De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification
© Yates et al.; licensee BioMed Central Ltd. 2014
Received: 29 November 2013
Accepted: 4 June 2014
Published: 9 June 2014
Red clover (Trifolium pratense L.) is a versatile forage crop legume, which can tolerate a variety of soils and is suitable for silage production for winter feed and for grazing. It is one of the most important forage legumes in temperate livestock agriculture. Its beneficial attributes include ability to fix nitrogen, improve soil and provide protein rich animal feed. It is however, a short-lived perennial providing good biomass yield for two or three years. Improved persistency is thus a major breeding target. Better water-stress tolerance is one of the key factors influencing persistency, but little is known about how red clover tolerates water stress.
Plants from a full sib mapping family were used in a drought experiment, in which the growth rate and relative water content (RWC) identified two pools of ten plants contrasting in their tolerance to drought. Key metabolites were measured and RNA-Seq analysis was carried out on four bulked samples: the two pools sampled before and after drought. Massively parallel sequencing was used to analyse the bulked RNA samples. A de novo transcriptome reconstruction based on the RNA-Seq data was made, resulting in 45181 contigs, representing ‘transcript tags’. These transcript tags were annotated with gene ontology (GO) terms. One of the most striking results from the expression analysis was that the drought sensitive plants were characterised by having approximately twice the number of differentially expressed transcript tags than the tolerant plants after drought. This difference was evident in most of the major GO terms. Before onset of drought the sensitive plants overexpressed a number of genes annotated as senescence-related. Furthermore, the concentration of three metabolites, particularly pinitol, but also proline and malate increased in leaves after drought stress.
This de novo assembly of a red clover transcriptome from leaf material of droughted and non-droughted plants provides a rich source for gene identification, single nucleotide polymorphisms (SNP) and short sequence repeats (SSR). Comparison of gene expression levels between pools and treatments identified candidate genes for further analysis of the genetic basis of drought tolerance in red clover.
KeywordsDrought stress Polymorphism Red clover RNA-Seq Transcriptome assembly Trifolium pratense
Red clover (Trifolium pratense L.) is a versatile forage crop legume, which can tolerate a variety of soils and is suitable for silage production for winter feed and for grazing. The benefits of red clover to farming include nitrogen fixation, soil improvement and high nutritive value in terms of protein-rich feed for livestock. Like other legumes, nitrogen fixation is facilitated by nodulation via symbiosis with the soil microbe Rhizobium leguminarosum[1, 2]. The nutritional benefits are attributable to easy digestibility, high voluntary intake by livestock and high protein content during ensiling . Red clover can be grown alone or in a sward mixed with grasses and other legumes, where it has been shown to be more productive than monocultures . Red clover is a short lived perennial that usually persists for two to three years, although more persistent varieties are available. Improving its longevity is a major breeding target [2, 5]. Persistency is a complex trait, and both biotic and abiotic stresses are known to have a major effect on longevity in the field . Red clover has a tap root system in the first year, but in subsequent years a more fibrous root system develops at the expense of the tap root, which senesce . This is believed to increase susceptibility to diseases of the crown, such as crown rot . Abiotic stresses are other factors considered to be of major importance for persistency [1, 3, 5]. Cold and drought stresses are important components of abiotic stress [6–10]. Even temperate forage crops are exposed to periods of drought that have adverse effects on yield and persistency . However, little is known about its genetic basis in red clover. Two studies have identified quantitative trait loci (QTLs) for traits related to persistence in red clover. Using full sib mapping families  found a major QTL for all persistency indices measured on linkage group 3, with heritabilites varying between 0.28 and 0.66. Klimenko et al.  found QTLs related to disease resistence and winter hardiness particularly on linkage groups 3 and 6. These studies are useful for future marker assisted selection. RNA-Seq approaches have potential to provide further valuable information about the molecular mechanisms underlying plant responses to abiotic stresses.
Effects of drought on gene expression have been studied using microarrays in rice , wheat  and Arabidopsis thaliana. These and other studies have provided information about the signalling pathways involved in the response to water stresses [17–19]. A typical response to water stress is to prevent water loss by closing the stomata, and producing abscisic acid (ABA). ABA initiates a signalling cascade which modifies the transcriptome and up-regulates genes encoding a number of proteins and enzymes involved in drought response [20, 21]. These can be divided into signalling and response proteins. Signalling proteins include transcription factors, protein kinases and ubiquitin ligases. Drought response proteins include dehydrins, heat shock proteins, aquaporins and scavengers of reactive oxygen species (ROS) [17, 22–25]. Additionally, metabolic pathways can be modified to compensate for reduced water and CO2 uptake . Osmolytes such as proline and pinitol are often produced to lower the osmotic potential in order to maintain water uptake and protect against build-up of toxic ion levels .
Next generation sequencing (NGS) can provide a depth of sequencing that is sufficient to cover the transcriptome of an organism many fold and allow quantification of the detected transcripts. The Illumina platform has been used for transcriptome analysis in several plant species [27–32]. An Illumina HiSeq2000 instrument can produce hundreds of millions of paired end reads per flow cell [33–35]. Gene expression is quantified by counting the number of reads per kilobase of transcript per million mapped reads (RPKM), enabling quantitative comparisons to be made . A complication arises with short reads if there is no reference sequence to map the reads onto. In the case of model organisms extensive genomic and transcript sequence data are available, but for non-model organisms the amount of expressed sequence tags (ESTs) and other transcript sequence data available to use as a reference can vary dramatically. Paucity of data in the reference transcriptome can potentially undermine the unbiased potential of NGS transcriptome analysis. Fortunately, software specifically designed for de novo reconstruction of the transcriptome from short reads, is available. A variety of methods using short read assemblers have been described in a number of plant species including Eucalyptus, Sonneratia alba, chickpea , sweet potato , alfalfa  and lupin .
The primary aim of this work was to use NGS technology to study changes in transcriptome patterns in pooled samples of red clover genotypes, contrasting in their phenotypic response to drought. Previous examples of transcriptome analysis have often described effects on seedlings [14, 16, 40] or other shorter term drought treatment. In this work we used mature plants in an experimental set-up designed to mimic field conditions as closely as possible. The plants were F1 progeny of a full sib mapping family. Based on growth rates and relative water content (RWC) we identified two bulks of genotypes contrasting in their response to drought stress. The Illumina platform was used to sequence four pools of RNA samples, two from the drought tolerant bulk before and after drought stress, and two from the sensitive bulk. It enabled us to assemble de novo a transcriptome library of red clover from paired end reads. As this study was made on genotypes from a mapping population, it also provided an opportunity to mine the pooled samples for putative SNP and SSR polymorphisms. Finally, information about differentially expressed genes informs our understanding of the fundamental aspects of drought response and identifies potential targets for improved abiotic stress tolerance. As part of this study we have also measured the concentration of a few key metabolites, known from other studies to undergo major changes in concentration in response to water stress.
The mapping family described in Methods was used to select a subset of genotypes that contrasted phenotypically in their response to drought stress. Leaf material from the two contrasting pools were used as a source for RNA-Seq data to generate a comprehensive library of transcript tags, and to explore how the two pools responded to drought stress in terms of differences in transcript expression patterns.
Selection of drought tolerant and sensitive genotypes
Data for plant growth rate, key metabolites, relative water content and osmotic potential
Sensitive (mg day-1)
Tolerant (mg day-1)
P = 0.119
Source of variation
-1.09 (-1.20-(-0.99)) MPa
-1.14 (-1.20-(-1.09)) MPa
-3.60 (-4.14-(-3.06)) MPa
-2.44 (-2.86-(-2.03) MPa
Red clover transcriptome assembly
The total RNA from both the drought and control treatments for individual plants of the two selected groups was extracted and quality checked. Out of the 40 RNA samples, two failed RNA quality check, one from the sensitive control pool and one from the sensitive drought treatment pool. The RNA samples which passed quality check were then combined in equal amounts in their respective pools and sent to Eurofins Ltd Genomic Service for sequencing. Over 100 M of 100 bp paired end reads were generated after quality checking by Eurofins. De novo transcriptome reconstruction of the data was made using RNA-Seq reads and publically available red clover EST sequences. Given the mixed sequencing data, i.e. the short reads produced in this work and the publically available red clover EST sequences (mostly Sanger sequences), this constitutes a hybrid assembly. This approach generated 45181 contigs (denoted ‘transcript tags’ hereafter) constituting 42 Mbp. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession number GAOU00000000. The version described in this paper is the first one, GAOU01000000. The full list of all transcript tags with their annotation and assignment to GO terms is available in Additional file 3. Other relevant statistics include an average transcript tag length of 933 bp, N50 of 622 bp, 196.4× coverage (actual number of reads mapped/transcriptome length) and longest sequence of 13855 bp. To further evaluate the quality of the assembly, all of the reads were mapped back to the assembled transcriptome using CLC Genomics Workbench v4.0.0. Using a paired end distance of up to 1000 bp 69% of the reads mapped back. A total of 34534 transcript tags (75%) were assigned a functional name by BLAST against Medicago, Arabidopsis or UNIPROT, and 29189 transcript tags (63%) were assigned at least one GO term.
Exploration of gene expression
Differentially expressed transcript tags in drought stressed red clover leaves compared with control conditions
While the effect of drought on gene expression was the main focus, we also looked at differences in gene expression between tolerant and sensitive pools under control conditions. The results are shown in Additional file 4 in the ‘control’ tab. The tolerant pool over-expressed 163 (>2 fold) transcript tags compared to the sensitive, which in turn had 436 up-regulated transcript tags. We inspected the data manually for gene functions or pathways and in the tolerant pool we found no consistent patterns in the data and interpret this as natural variation due to transcriptome plasticity. In the sensitive pool we identified eight senescence related proteins (RC.3083, 4561, 7070, 7319, 7433, 10857, 19071, 40938), which were up-regulated compared to the tolerant genotypes prior to drought stress.
SNP and SSR mining
Breakdown of SSR repeats present in the red clover transcriptome assembly
Number of seqs searched
Total size of sequences (bp)
41 943 532
Number of SSRs
Number of SSR containing seqs
Number of seqs containing more than one SSR
Number of compound SSRs
Distribution of SSR repeat types
Di-nucleotides (≥8 repeats)
Tri-nucleotides (≥6 repeats)
Tetra-nucleotides (≥5 repeats)
Penta-nucleotides (≥5 repeats)
Hexa-nucleotides (≥5 repeats)
The metabolite analysis data summarised in Table 1C show that pinitol, proline and malate all increased in concentration following exposure to drought, while glucose, fructose and myo-inositol remained unchanged. Particularly pinitol concentration rose dramatically from a little over 1.5% of dry matter to over 10% during the course of drought exposure. While proline concentration also rose significantly, it was from a very low base (undetectable), and the overall concentration of this metabolite was much lower than pinitol. Malate concentration was also measured. It can be formed by combining CO2 with phosphoenol pyruvate under conditions of limiting CO2 availability, such as stomatal closure [42, 43]. In C4 plants malate is shuttled between mesophyll and bundle sheath cells, where it is decarboxylated, and recycled. Given that water stress leads to stomatal closure and decreased CO2 availability, it seemed reasonable to look for changes in malate concentration in response to drought stress. While the concentration of malate increased (Table 1C), we did not observe convincing evidence for significant changes in expression of transcript tags encoding any of the key enzymes in CO2 metabolism which would be necessary for a malate shuttle. This includes phosphoenolpyruvate carboxylase, NADP-dependent malic enzyme, malate dehydrogenase and pyruvate-orthophosphate dikinase.
This work has generated a large amount of sequence data with NGS technology in a non-model crop species. Although red clover previously had a number of resources available, such as an EST database  the work described here, has greatly increased the number of transcript sequences available. The high number of sequenced reads resulted in significant coverage of the transcriptome in leaves, allowing gene discovery, quantification of transcripts in four pools and the identification of putative SNPs and SSRs.
In this experiment we used mature plants in a long term (60 days) experiment to simulate the effects of drought as realistically as possible in a greenhouse environment. Many drought experiments described in the literature involving transcriptomics have used seedlings and/or short term drought treatment which minimises variation due to environmental factors (e.g. [14, 16]). While long term experiments are likely to introduce more variation between genotypes and replicate plants, it may also provide potentially novel information on how the plant responds to prolonged periods of water deficiency. A large number of transcript tags were differentially expressed in response to the drought stress in both the sensitive and tolerant plants (Table 2), but variation in expression levels (the average CV of RPKM values of pairs of tags from the sensitive and tolerant pools) was slightly larger before than after drought. We used a pooling strategy to reduce within and between pool differences. This approach has been widely used in RNA-Seq experiments to study for example cell wall composition in Alfalfa , berry development in Vitis viniferia, tissue specific expression in Eucalyptus and floral sex determination in cucumber . Barcoding of the four pools enabled multiplexing libraries on a single lane of the flow-cell. This reduced possible confounding effects arising from between-flow-cell differences . It is also interesting to note that expression levels before drought were similar in both pools for 99% of the transcript tags. This suggests that the differences in gene expression which we observed in the two sets of genotypes at the end of the drought stress period are likely to represent a response to the drought stress rather than a difference between the two pools per se.
In the qRT-PCR experiment three out of four of the transcript tags tested showed similar changes in expression as was observed in the RNA-Seq experiment. The reasons for the discrepancy between the two methods for RC.31500 are not clear to us. While our data represent a small sample, the results are however of a similar nature to those reported in . They found that 12 out of 15 genes matched the expected RPKM values when quantified using qRT-PCR.
The results reported here represent a platform for further examination of candidate genes with a potential role in plant response to drought stress. The total size of the transcriptome assembled was 42 Mb. This constitutes around 10% of the predicted genome size (440 Mb) . Approximately 46 k transcript tags were generated in the transcriptome reconstruction. This is similar to those found in Medicago truncatula (44124 gene loci) , Soybean (46430 protein coding genes)  and Arabidopsis (37019 gene models) . The number and total size of transcripts naturally vary for de novo assemblies of other species, such as in chickpea (55 k) transcripts (28 Mb) , Sonneratia alba, 30 k transcripts (18 Mb)  and Eucalyptus, 18 k transcripts (22 Mb) . An 8 fold coverage filter on contigs was used in Eucalyptus, and a 2× minimum coverage in S. alba. This would explain the smaller numbers found in those transcriptomes. Using the same filters for RPKM a similar number of transcript tags would be realised in the present experiment. Due to the heterozygous and allogamous nature of red clover it is likely, that there is some redundancy in our transcriptome assembly. The transcriptome size may thus be over-estimated. On the other hand, erroneously merged assemblies would lead to underestimation of the transcript numbers. The SNP detection process showed that 462 transcript tags had greater than 10 SNPs (Additional file 5), which may suggest that the assembly process in those cases merged paralogues or that the reads could be miss-aligned. The allogamous nature of red clover and the pooling strategy used in the experiment makes it challenging to distinguish with certainty between allelic variants and paralogous genes. Therefore it is reasonable to assume there are a few inconsistencies in the assembled transcriptome. Nevertheless, the majority of the transcript tags do have high coverage and have high similarity to known transcripts.
In this experiment the level of water stress as measured by soil water content was identical across bins at the time of sampling during the drought treatment (Additional file 1). Nevertheless, the sensitive and tolerant plants showed contrasting responses to the same level of soil moisture. The sensitive genotypes had a RWC of approximately 43% (Figure 1). In contrast, the tolerant genotypes had a RWC of 78%, which is only slightly below what is found in non-droughted plants (~85%) (Table 1C). This difference manifested itself by the observation that more than twice the number of transcript tags was differentially expressed in the sensitive pool compared to its control, than was the case in the tolerant pool compared to its control. It would be difficult to infer which differentially expressed transcript tags qualify as genes directly involved in response to the drought stress, and as candidates for targeting in genetic improvement programmes. Nevertheless, the results do provide a comprehensive overview of gene expression changes after drought stress. The number of differentially expressed transcript tags which were detected as a result of the drought treatment (Additional file 4) are similar to those reported in other species such as rice (3097 up, 2391 down) , wheat (3056 differential)  and Arabidopsis (2059 up and 2075 down) , albeit for shorter term drought stress.
There was a large degree of overlap between GO terms of differentially expressed transcript tags between the sensitive and tolerant pools in their response to drought. The main difference is that a larger number of transcript tags from the sensitive pool with the same GO terms were either up-regulated or down-regulated in response to the drought, compared to the tolerant pool (Figure 4). The reason for this is not clear, but could simply be a consequence of the degree of severity of the stress experienced by the plants in the two pools. Down-regulation of the photosynthetic apparatus also involves more genes in the sensitive pool than the tolerant, possibly also caused by the difference in sensitivity to the stress (Figure 4). The results contrast with those of , who found that more genes were up-regulated in tolerant barley genotypes than in sensitive. Comparisons with short term drought experiments are however, difficult to make. The gradual onset of drought conditions in the present experiment may elicit the differential expression of more genes to counteract the severity of the stress. Nevertheless, even the relatively modest reduction of ~10% in RWC observed in the tolerant pool, elicited a significant transcriptome response albeit after prolonged water limitation.
Another finding is the high number of senescence-related transcript tags that were up-regulated in the sensitive pool in non-droughted plants. Perhaps the sensitive plants have a predisposition for reduced growth and less stress tolerance . Chaves et al. suggested that plants not predisposed to early senescence may be of value for breeding purposes. Identification of molecular markers or other biomarkers associated with these differences could potentially assist in the removal of such sensitive plants from the gene-pool in a population-based breeding programme. In two of those transcripts (RC.7319 and RC.19071) SNPs were identified (Additional file 5).
The transcript tags which were differentially expressed are similar to those from other drought studies where large transcriptome profiling was used [14–16]. They can be divided into three main groups: functional, regulatory and photosynthesis related genes. Functional drought response genes include protection enzymes (LEA, heat shock proteins, ROS scavenging), osmolyte biosynthesis (proline, trehalose, pinitol/ononitol), transporters (aquaporins, sugar transporters). Regulatory proteins include transcription factors (MYB, NAC, DREB and zinc finger), and those involved in post-transcriptional modification (splice factors, DEAD-box) and post-translational modification, (ubiquitin ligases, protein kinases) and epigenetic factors (histones) [15–17]. When plants are unable to use or dissipate absorbed light energy, the excess energy leads to production of ROS which can cause oxidative damage to the photosynthetic machinery . Under prolonged drought stress, the photosynthetic machinery is down-regulated . This is also the case in this study, where particularly transcript tags encoding Rubisco and polypeptides involved in PSII and PSI were down-regulated (Additional file 4). The rate of recovery of the plant after drought is tightly linked to photosynthetic recovery . Therefore plants which maintain larger portions of their photosynthesising apparatus would have an advantage over those with a lower portion once the stress is alleviated.
Three metabolites known to accumulate in plants in response to drought were measured in this work. Pinitol accumulation was most dramatic and statistically significant (Table 1C). We therefore tried to identify a candidate gene likely to represent the O-methyl transferase activity required for the key step in the biosynthesis of pinitol from myo-inositol . The greatest homology with the orthologue in the ice-plant (Mesembryanthemum crystallinum) was RC.4122, but when we did a BLAST search with that as the query, most of the best hits were to genes annotated as caffeic acid-O-methyl transferase. Furthermore, expression of RC.4122 did not increase in response to drought. One explanation could be that the accumulation of pinitol took place early in the drought stress period, and expression of the gene had declined by the end of the drought period. That seems unlikely, since transcript tags annotated as two of the key enzymes in proline biosynthesis were over-expressed in both the sensitive and the tolerant pool in response to drought (e.g. RC.34920: pyrroline-5-carboxylate reductase (P ≤ 7.99 × 10-5), and RC.35692: Δ1-pyrroline-5-carboxylate synthetase (P < 10-200) (Additional file 4). However, the similarity between most O-methyl transferases makes accurate annotation more challenging. It is surprising to us though, that this gene does not appear to have been identified in any legume despite the high concentration of pinitol and other polyols in this family [53–55].
This work makes available the first transcriptome library of red clover using next generation sequencing technology. The library was constructed using two pairs of pools of ten genotypes from an F1 mapping family. RNA from the two pairs of pools (drought sensitive and drought tolerant, as determined by phenotypic analysis) was extracted before and at the end of a 60 day period of drought. The RNA-Seq analysis provided a total of 45181 transcripts longer than 200 bp. A large number of transcript tags involved in the photosynthetic apparatus were down-regulated in both the sensitive and tolerant pools, but particularly in the former. A great many transcription factors, genes involved in protein modification and degradation were up-regulated in both pools after drought, but again, particularly in the sensitive pool. As expected, a number of genes involved in ABA biosynthesis and abiotic stress signal transduction were also up-regulated. The general trend was that a larger number of genes were up-regulated in the sensitive pool than in the tolerant. Despite the highly significant increase in the concentration of the compatible solute pinitol, we were unable to definitively identify the gene encoding the O-methyl transferase responsible for the key step in pinitol biosynthesis from myo-inositol. We recognise the limitations of the quantitative aspects of this work, and that further validation of the assembly is needed, but the results presented here provide a valuable resource for future work in terms of SNP identification, and annotation of genome sequence assembly currently in progress.
Plant material and experimental setup
The mapping family used throughout consisted of progeny from a F1 pseudo-testcross derived from a single genotype from each of two varieties, Milvus and Britta. The plant material was generated initially using a clonal cutting method from the original plant material involving excision of stem regions with nodes followed by auxin treatment and rooting in soil. Subsequent clones were generated by clonal splitting . We made use of 64 genotypes from the mapping family. Three blocks were used in the experimental design, where each block constituted one clone of each of the 64 genotypes. Each block consisted of two drought bins in a randomised block design, i.e. 32 plants per bin and 64 plants per block. The ‘drought bins’ measured 80 cm by 110 cm, with a depth of 80 cm. The bins were filled with gravel up to 30 cm: the remaining 50 cm was filled with John Innes Number 3 compost. Prior to the start of the experiment the bins were watered to field capacity.
Plant material was grown to mature size in 15 cm pots before transfer to drought bins on 06/03/2010. Plants were then allowed to adjust 2 days to the new environment before cutting back to approximately 6 cm height above soil level. Then the bins were watered daily to maintain field capacity (~55% soil moisture content) until the plants reached a mature size on 15/04/2010, at which time they were sampled for multiple phenotypic traits as described below. This was designated as DW0. From then onwards all watering ceased. Plants were cut again on 05/05/2010, and dry weight was recorded. This was designated as DW1. The subsequent sampling time for individual bins (designated DW2) was done when soil moisture content was down to 14% (see next section for details). This happened between 13-18/06/2010. All phenotypic traits were sampled at this point. Subsequently, the bins were fully hydrated and continually watered until the last cutting on 16/07/2010. This was designated DW3. Growth rate was quantified as dry weight divided by the number of days of growth since the previous cut.
Soil moisture, relative water content and osmotic potential
Two access tubes for Delta-T PR1 Profile probes were inserted into each drought bin. They were placed equidistant from the edge and 30 cm apart. The soil moisture content was measured using Delta-T PR1 Profile probe and Delta-T HH2 moisture meter at soil depths of 10, 20, 30 and 40 cm. The mean of the 4 measurements was used to indicate the water content of the soil in the bin. In addition, five plants per bin were sampled randomly at 8:30 am for measurement of relative water content (RWC). A single leaf sample was excised from each plant and placed in a sealable zip-lock bag on ice before quantifying RWC using the method described by Smart and Bingham . The soil moisture content and RWC were used as indicators of the drought status of the bins and plants, respectively. The final sampling was done when the soil moisture average fell below 14% and/or the RWC was below 60%.
Samples for measurement of osmotic potential (OP) were collected between 15:00 and 16:00, by excising mature leaf material and storing at -20°C. OP was then quantified using a Wescor vapour pressure osmometer 5520 (Wescor Inc., USA) as described by the manufacturer.
Samples for RNA extraction were collected between 13:00 and 15:00. The youngest fully open trifoliate leaf was used for this. The excised material was stored in 2 ml safe-lock Eppendorf tubes and then flash frozen in liquid N, before storing at -80°C. The leaf sample was ground to a fine powder using a Retsch MM300 mixer mill. RNA was then extracted using the Trizol (Invitrogen, USA) method as described by the manufacturer, with an additional chloroform extraction step added. Total RNA was solubilised in 50 μl DEPC-treated H2O. A 2 μl aliquot was taken and diluted 1:5 in DEPC treated water and used for RNA quantification. Quantification was done using BioRad Experion Automated Electrophoresis System with RNA StdSens chips (Bio-RAD Laboratories Inc., USA), as described by the manufacturer using 3 μl of the diluted RNA extract. The data were analysed using Experion™ Software System operation and data analysis tools, version 22.214.171.124 copyright © 2010 Bio-Rad laboratories. Four pooled samples were generated (see Results section), two consisting of 9-10 drought sensitive genotypes sampled before and at the end of the drought treatment, and two consisting of 9-10 drought tolerant genotypes exposed to the same treatments. For each pool, 2 μg of total RNA from each sample was included. Then 500 μl of 100% ethanol was added to each pool, and kept at -20°C during transport.
Sequencing and transcriptome construction
The four RNA pools were sent to Eurofins Genetic Services (Eurofins, Germany), for sequencing. This included preparation of 4 barcoded libraries, which were multiplexed on a single flow-cell of an Illumina HiSeq2000 (Illumina, USA). The sequencing was done as paired-end reads, 2×100 bp in length. The data was quality checked at Eurofins and provided as four libraries in FASTq format. Sequence quality controls were applied by Eurofins and the read libraries were checked using fastqc (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). After assembly the contigs were checked for adaptor sequences and none were found.
For the initial construction of the reference transcriptome Velvet  was used. Assemblies were created using k-mer sizes of 31, 41, 51 and 61. The 61-k-mer length assembly was chosen for further use because it gave the largest N50, the largest average scaffold size, the largest size of any single scaffold, and the smallest number of scaffolds. This generated 93.4 k contigs. To increase the length of contigs and reduce the number of redundant contigs the Velvet assemblies were combined with existing sequences from the databases from the publically available Trifolium pratense EST library, downloaded from http://www.plantgdb.org/download/download.php?dir=/Sequence/ESTcontig. The red clover ESTs were combined with the Velvet assemblies by using a second de novo assembly on CLCbio Genomics Workbench version 4.0.0 (31 k-mer). This created 47229 contigs of at least 200 bases in length. In order to retain polymorphism data, conflicting bases found during the assembly process were called ‘ambiguously’, rather than using the ‘vote’ option. Next, bacterial contigs were removed, by compiling over 4000 bacterial genomes from the NCBI RefSeq database. Then the transcriptome contigs were mapped against the bacterial genomes and un-mapped (non-bacterial) contigs were retained. This led to the removal of 634 contigs. However, to remove the ambiguity for submission to NCBI the reads were mapped back to the contigs to give a ‘vote’ scoring.
Transcript tag expression was determined using the method described by  and executed in CLC. For identification of differentially expressed transcript tags statistical analysis described by  was used. A false discovery rate (FDR) of 0.05 was used as the threshold for significance.
Annotation and gene ontology analysis
BLASTX was used with minimum E-value 1e-15 for annotation of transcript tags. Sequentially, contigs were searched using BLASTX against the following protein databases Medicago truncatula (v3.5, ), Arabidopsis thaliana (TAIR10, ) and the Universal Protein Resource (UniProt, [60, 61]). From each BLAST output the best match was used to annotate the transcript tags. After each BLAST search annotation tags with no matches and ones with ‘hypothetical’, ‘predicted’ or ‘unknown’ annotations were extracted for next sequential BLAST search (order: Medicago - Arabidopsis – UniProt). For GO assignment corresponding GO terms were downloaded from relevant databases. A small number of tags annotated as fungal were removed by manual curation.
TopGO from Bioconductor in R (http://www.r-project.org/) was used to identify enriched GO terms. This was done for all differentially expressed genes, split between over- and under-expressed transcript tags identified by CLC, for each pairwise comparison. The Fishers test implemented in TopGO was used to identify enriched GO terms relative to its expectation. GO terms with a P-value < 0.01 were selected. Both common and unique GO terms were identified by intersection of the data, including those unique to individual pools. Cytoscape  was used to visualise the selected GO terms. The GO map for A. thaliana was downloaded and installed in Cytoscape, the enriched GO terms were selected and a sub-network was created for visualisation of the GO network.
Aliquots of the RNA extracts were DNAse treated with the Ambion DNA-Free™ kit (Ambion, USA), according to the manufacturer instructions. The cDNA synthesis was done using Invitrogen™ Superscript™ II reverse transcriptase (RT) (Invitrogen, USA) as described by the manufacturer using an Oligo(dT)25 Primer with 50 ng total RNA. Samples were then treated with RNaseH (Promega, USA) and stored at -20°C. All qRT-PCR work was done using an Applied Biosystems (USA) 7500 Real time PCR system. The 7500 system sequence detection software version 1.2.3 was used for analysis including auto Ct/CP (cross threshold/crossing point), mean CP and outlier removal at 95% confidence). All samples were run in triplicate using the following PCR concentrations in a volume of 10 μl: 5 μl SYBR Green master mix (Applied Biosystems, USA), 900nM forward and reverse primers (Biolegio, Netherlands), and 1 μl 10× diluted cDNA. All primers were designed using the Roche Universal ProbeLibrary Assay Design Center at https://www.roche-applied-science.com/sis/rtpcr/upl/index.jsp?id=UP030000. The following amplification program was used: 95°C for 10 min, 40 cycles at 92°C for 15 s, annealing temperature for 30s, 72°C for 45 s followed by a final elongation at 72°C for 10 min. The specific annealing temperatures and nucleotide sequence of the primers are shown in Additional file 8. The relative expression software tool (REST-384) was used to calculate the relative fold change in gene expression (http://www.gene-quantification.info/)  using the mathematical model described by Pfaffl . Primer efficiencies were measured using a serial dilution of stock cDNA at 1:1, 1:10 and 1:100. Every plate included three negative controls for detection of contamination. RC.44146 was used as the reference gene. This was the Actin gene, which was found not to be differentially expressed in the RNA-seq data analysis. For comparison between RNA-seq and qRT-PCR data, the fold change in expression of drought stressed samples was calculated relative to control samples for both RNA-seq and qRT-PCR data.
SNP and SSR mining
The SNP detection tool in CLC was used to identify SNPs. All reads were included in this process. The minimum count for the variant allele was 200 and window length was 17 bp. For identification of SNPs inherited from one heterozygous parent and one homozygous parent, the minimum variant frequency was set to 20% to allow variation from the expected haplotype frequency of 25/75%.
The MISA (Microsatellite) Perl script (http://pgrc.ipk-gatersleben.de/misa) was used for identification of SSRs. The settings for minimum number of repeats were as follows: Di-nucleotide 8; Tri-nucleotide 6; and Tetra, Penta and Hexa-nucleotide repeats 5. For compound SSRs the maximum distance between two SSR runs was 50 bp. The BatchPrimer3 programme  was used to design primers pairs for amplification of the SSR motifs. The default settings were used except for the annealing temperature, which was set for an optimum of 60°C.
Metabolite sampling and extraction
Whole young trifoliate leaf samples were collected for metabolite quantification between 11:00 and 13:30. They were stored on ice in re-sealable plastic bags, until an amount of 20-60 mg was weighed. This was stored in 2 ml safe-lock microfuge tubes (Eppendorf, Germany) before they were flash frozen in liquid N2 and then stored at -20°C. Before extraction a metal ball bearing was added to each tube. During this process the samples were kept frozen by brief exposures to liquid N2. This was also the case when the samples were homogenized to a fine powder in a Retsch-mill MM300, mixer mill (Retch, Germany) for 2 min (20 impacts/s). Afterwards 300 μl of 100% ethanol was added and mixed for a further 15 minutes (20 impacts/s). Next, 200 μl of chloroform (CHCl3) was added and mixing continued for 5 minute (20 impacts/s). Then, 400 μl H2O was added and the samples were mixed by vortexing, before centrifugation for 5 min at 14 k rpm to separate the aqueous ethanol phase from the chloroform phase. The total volume of the aquous ethanol phase was approximately 525 μl. Two 75 μl-aliquots were transferred to new microfuge tubes (1.5 ml). The aliquots were dried using a speed vac concentrator for four hours (without heating), with the lids open under vacuum. All samples were then stored at 4°C.
Determination of metabolite concentration using gas chromatography mass spectroscopy
Prior to derivatization of the samples two internal standards were added: castanospermine (Enzo life sciences, USA) (20 mg) and cyclo-leucine (Sigma-Aldrich, UK) (20 mg). The samples were dried down completely before derivatizing. Derivatization of the samples was achieved using a 2-step process. First, 30 μl of freshly prepared MOX solution (O-methylhydroxylamine hydrochloride in pyridine, Sigma-Aldrich, UK) was added to the dried samples, and kept at room temperature for 2 hours to allow solubilisation. The samples were then transferred to crimp sealed 0.3 ml GC vials (Chromacol/Fisher, UK), and the vials were heated for 15 minutes at 90°C to convert ketone groups to the oxime derivative. The vials were then de-capped and 20 μl BSTFA (N,O-Bis(trimethylsilyl) trifluoroacetamide) (Sigma-Aldrich, UK) was added. The vials were recapped and heated for a further 15 minutes at 90°C to convert remaining polar groups to their trimethylsilyl derivatives. This method is similar to those described by Parveen et al..
Analysis was carried out using an Agilent GC-MS system (Agilent Technologies UK Ltd.) comprising a 5973 network Mass Selective detector, an 6890 Series GC and a 7683 series autosampler. The GC was fitted with a Varian FactorFour™ VF-5 ms capillary column (30 m × 0.25 mm ID × 0.25 μm film thickness). The injection volume was set to 1 μl with a 1:50 split ratio. The oven was set at an initial temperature of 80°C and increased to 280°C at a rate of 10°C per minute. The inlet temperature was set to 280°C and the transfer line to 320°C. The mass spectrometer was set to scan between 50 and 600 m/z. For quantification of metabolites standard curves were prepared using standards (Sigma-Aldrich, UK with the exception of pinitol which was purchased from ACROS Organics, Belgium). Quantification was achieved using Enhanced Data Analysis software (ChemStation 1701 CA version 00.00, Agilent Technologies UK Ltd.). Additional file 9 shows the metabolites that were quantified and relevant information for quantification. The quantification of glucose, fructose, myo-inositol and pinitol in non-droughted samples required standard curves ranging from 1-50 μg. In droughted samples the higher concentrations of pinitol required the range to be increased to 500 μg. An additional standard curve was set up for pinitol concentrations above 200 μg since the curve was non-linear in the high range.
Statistical analysis of plant growth and metabolites
The plant growth data were analysed using a repeated measures ANOVA approach with the control data (DW0) as a co-variate. The data were log10 transformed. Because there were only two timepoints for the RWC, OP and metabolite data, we used a split-plot in time ANOVA design for the analysis of those data. The data for glucose, fructose, myo-inositol, malate and pinitol were log10 transformed. The analysis was performed using Genstat, 15th edition (VSN International Ltd (http://www.vsni.co.uk)).
Availability of supporting data
All the sequencing data have been deposited in the NCBI short read archive under the Bioproject PRJNA219226 (http://www.ncbi.nlm.nih.gov/bioproject/PRJNA219226). Other supporting data are included as additional files.
Basic local alignment search tool
Dehydration responsive element binding
Late embryogenesis abundant
Nicotinamide adenine dinucleoside phosphate
- PSI & PSII:
Photosystem I & II
Quantitative real time polymerase chain reaction
Reactive oxygen species
Ribulose bisphosphate carboxylase
The Arabidopsis Information Resource.
This work was funded by the Biotechnology and Biological Sciences Research Council through the Crop Science 2007 Initiative, and an Institute Strategic Programme Grant (IBERS/08/2). We wish to thank Rhys Kelly for technical assistance, and Ruth Sanderson for advice on the statistical analysis of the data presented in Table 1 and Additional file 2.
- Taylor NL, Quesenberry KH: Red Clover Science. 1996, Dordrecht: Kluwer Academic PublishersView ArticleGoogle Scholar
- Abberton MT, Marshall AH: Progress in breeding perennial clovers for temperate agriculture. J Agric Sci. 2005, 143: 117-135. 10.1017/S0021859605005101.View ArticleGoogle Scholar
- Frame J, Charlton JFL, Laidlaw AS: Temperate Forage Legumes. 1998, Wallingford: CAB InternationalGoogle Scholar
- Nyfeler D, Huguenin-Elie O, Suter M, Frossard E, Connolly J, Lüscher A: Strong mixture effects among four species in fertilized agricultural grassland led to persistent and consistent transgressive overyielding. J Appl Ecol. 2009, 46: 683-691. 10.1111/j.1365-2664.2009.01653.x.View ArticleGoogle Scholar
- Taylor NL: A century of clover breeding developments in the United States. Crop Sci. 2008, 48: 1-13. 10.2135/cropsci2007.08.0446.View ArticleGoogle Scholar
- Thomashow MF: molecular basis of plant cold acclimation: insights gained from studying the CBF cold response pathway. Plant Physiol. 2010, 154: 571-577. 10.1104/pp.110.161794.PubMed CentralPubMedView ArticleGoogle Scholar
- Thomashow MF: Plant Cold Acclimation: Freezing Tolerance Genes and Regulatory Mechanisms. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 571-599. 10.1146/annurev.arplant.50.1.571.PubMedView ArticleGoogle Scholar
- Thomashow MF: Role of Cold-Responsive Genes in Plant Freezing Tolerance. Plant Physiol. 1998, 118: 1-8. 10.1104/pp.118.1.1.PubMed CentralPubMedView ArticleGoogle Scholar
- Rizhsky L, Liang H, Mittler R: The combined effect of drought stress and heat shock on gene expression in tobacco. Plant Physiol. 2002, 130: 1143-1151. 10.1104/pp.006858.PubMed CentralPubMedView ArticleGoogle Scholar
- Rasmussen S, Barah P, Suarez-Rodriguez MC, Bressendorff S, Friis P, Costantino P, Bones AM, Nielsen HB, Mundy J: Transcriptome responses to combinations of Stresses in Arabidopsis. Plant Physiol. 2013, 161: 1783-1794. 10.1104/pp.112.210773.PubMed CentralPubMedView ArticleGoogle Scholar
- Peterson PR, Sheaffer CC, Hall MH: Drought effects on perennial forage legume yield and quality. Agron J. 1992, 84: 6-Google Scholar
- Herrmann D, Boller B, Studer B, Widmer F, Kölliker R: Improving persistence in red clover: insights from QTL analysis and comparative phenotypic evaluation. Crop Sci. 2008, 48: 269-277. 10.2135/cropsci2007.03.0143.View ArticleGoogle Scholar
- Klimenko I, Razgulayeva N, Gau M, Okumura K, Nakaya A, Tabata S, Kozlov N, Isobe S: Mapping candidate QTLs related to plant persistency in red clover. Theor Appl Genet. 2010, 120: 1253-1263. 10.1007/s00122-009-1253-5.PubMed CentralPubMedView ArticleGoogle Scholar
- Lenka SK, Katiyar A, Chinnusamy V, Bansal KC: Comparative analysis of drought-responsive transcriptome in Indica rice genotypes with contrasting drought tolerance. Plant Biotechnol J. 2011, 9: 315-327. 10.1111/j.1467-7652.2010.00560.x.PubMedView ArticleGoogle Scholar
- Aprile A, Mastrangelo A, De Leonardis A, Galiba G, Roncaglia E, Ferrari F, De Bellis L, Turchi L, Giuliano G, Cattivelli L: Transcriptional profiling in response to terminal drought stress reveals differential responses along the wheat genome. BMC Genomics. 2009, 10: 279-10.1186/1471-2164-10-279.PubMed CentralPubMedView ArticleGoogle Scholar
- Matsui A, Ishida J, Morosawa T, Mochizuki Y, Kaminuma E, Endo TA, Okamoto M, Nambara E, Nakajima M, Kawashima M, et al: Arabidopsis Transcriptome Analysis under Drought, Cold, High-Salinity and ABA Treatment Conditions using a Tiling Array. Plant Cell Physiol. 2008, 49: 1135-1149. 10.1093/pcp/pcn101.PubMedView ArticleGoogle Scholar
- Shinozaki K, Yamaguchi-Shinozaki K: Gene networks involved in drought stress response and tolerance. J Exp Bot. 2007, 58: 221-227.PubMedView ArticleGoogle Scholar
- Chaves MM, Flexas J, Pinheiro C: Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot. 2009, 103: 551-560.PubMed CentralPubMedView ArticleGoogle Scholar
- Chaves MM, Maroco JP, Pereira JS: Understanding plant responses to drought - from genes to the whole plant. Funct Plant Biol. 2003, 30: 239-264. 10.1071/FP02076.View ArticleGoogle Scholar
- Schroeder JI, Allen GJ, Hugouvieux V, Kwak JM, Waner D: Guard cell signal transduction. Annu Rev Plant Physiol Plant Mol Biol. 2001, 52: 627-658. 10.1146/annurev.arplant.52.1.627.PubMedView ArticleGoogle Scholar
- Agarwal PK, Jha B: Transcription factors in plants and ABA dependent and independent abiotic stress signalling. Biol Plant. 2010, 54: 201-212. 10.1007/s10535-010-0038-7.View ArticleGoogle Scholar
- Mittler R, Vanderauwera S, Gollery M, Van Breusegem F: Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9: 490-498. 10.1016/j.tplants.2004.08.009.PubMedView ArticleGoogle Scholar
- Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K: Monitoring the Expression Pattern of 1300 Arabidopsis Genes under Drought and Cold Stresses by Using a Full-Length cDNA Microarray. Plant Cell. 2001, 13: 61-72. 10.1105/tpc.13.1.61.PubMed CentralPubMedView ArticleGoogle Scholar
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, Satou M, Akiyama K, Taji T, Yamaguchi-Shinozaki K, Carninci P, Kawai J, Hayashizaki Y, Shinozaki K: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31: 279-292. 10.1046/j.1365-313X.2002.01359.x.PubMedView ArticleGoogle Scholar
- Seki M, Umezawa T, Urano K, Shinozaki K: Regulatory metabolic networks in drought stress responses. Curr Opin Plant Biol. 2007, 10: 296-302. 10.1016/j.pbi.2007.04.014.PubMedView ArticleGoogle Scholar
- Ashraf M, Foolad MR: Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot. 2007, 59: 206-216. 10.1016/j.envexpbot.2005.12.006.View ArticleGoogle Scholar
- Garg R, Patel RK, Tyagi AK, Jain M: De Novo Assembly of Chickpea Transcriptome Using Short Reads for Gene Discovery and Marker Identification. DNA Res. 2011, 18: 53-63. 10.1093/dnares/dsq028.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang SS, Tu Z, Cheung F, Xu W, Lamb J, Jung H-J, Vance C, Gronwald J: Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems. BMC Genomics. 2011, 12: 199-10.1186/1471-2164-12-199.PubMed CentralPubMedView ArticleGoogle Scholar
- Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of Transcriptional Complexity during Berry Development in Vitis vinifera Using RNA-Seq. Plant Physiol. 2010, 152: 1787-1795. 10.1104/pp.109.149716.PubMed CentralPubMedView ArticleGoogle Scholar
- Mizrachi E, Hefer CA, Ranik , Joubert F, Myburg AA: De novo assembled expressed gene catalogue of a fast growing Eucalyptus tree produced by Illumina mRNA-Seq. BMC Genomics. 2010, 11: 681-10.1186/1471-2164-11-681.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu T, Qin Z, Zhou X, Feng Z, Du Y: Transcriptome profile analysis of floral sex determination in cucumber. J Plant Physiol. 2010, 167: 905-913. 10.1016/j.jplph.2010.02.004.PubMedView ArticleGoogle Scholar
- Wang L, Li P, Brutnell TP: Exploring plant transcriptomes using ultra high-throughput sequencing. Briefings Funct Genom. 2010, 9: 118-128. 10.1093/bfgp/elp057.View ArticleGoogle Scholar
- Kircher M, Kelso : High-throughput DNA sequencing – concepts and limitations. Bioessays. 2010, 32: 524-536. 10.1002/bies.200900181.PubMedView ArticleGoogle Scholar
- Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92 (5): 255-264. 10.1016/j.ygeno.2008.07.001.PubMedView ArticleGoogle Scholar
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26: 1135-1145. 10.1038/nbt1486.PubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008, 5: 621-628. 10.1038/nmeth.1226.View ArticleGoogle Scholar
- Chen S, Zhou R, Huang Y, Zhang M, Yang G, Zhong C, Shi S: Transcriptome sequencing of a highly salt tolerant mangrove species Sonneratia alba using Illumina platform. Mar Genome. 2011, 4: 129-136. 10.1016/j.margen.2011.03.005.View ArticleGoogle Scholar
- Wang Z, Fang B, Chen J, Zhang X, Luo Z, Huang L, Chen X, Li Y: De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato (Ipomoea batatas). BMC Genomics. 2010, 11: 726-10.1186/1471-2164-11-726.PubMed CentralPubMedView ArticleGoogle Scholar
- O'Rourke JA, Yang SA, Miller SS, Bucciarelli B, Liu J, Rydeen A, Bozsoki Z, Uhde-Stone C, Tu ZJ, Allan D, Gronwald JW, Vance CP: An RNA-Seq transcriptome analysis of orthophosphate-deficient white lupin reveals novel insights into phosphorus acclimation in plants. Plant Physiol. 2013, 161: 705-724. 10.1104/pp.112.209254.PubMed CentralPubMedView ArticleGoogle Scholar
- Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Monitoring Expression Profiles of Rice Genes under Cold, Drought, and High-Salinity Stresses and Abscisic Acid Application Using cDNA Microarray and RNA Gel-Blot Analyses. Plant Physiol. 2003, 133: 1755-1767. 10.1104/pp.103.025742.PubMed CentralPubMedView ArticleGoogle Scholar
- Torii KU: Leucine-Rich Repeat Receptor Kinases in Plants: Structure, Function, and Signal Transduction Pathways. 2004, 234-Google Scholar
- Kanai R, Edwards GE: The biochemistry of C4 photosynthesis. C4 Plant Biology. Edited by: Sage RF, Monson RK. 1999, San Diego: Academic, 49-87.View ArticleGoogle Scholar
- Hatch MD: C4 photosynthesis: A historical perspective. C4 Plant Biology. Edited by: Sage RF, Monson RK. 1999, San Diego: Academic, 17-46.View ArticleGoogle Scholar
- Sato S, Isobe S, Asamizu E, Ohmido N, Kataoka R, Nakamura Y, Kaneko T, Sakurai N, Okumura K, Klimenko I, Sasamoto S, Wada T, Watanabe A, Kohara M, Fujishiro T, Tabata S: Comprehensive Structural Analysis of the Genome of Red Clover (Trifolium pratense L.). DNA Res. 2005, 12: 301-364.PubMedView ArticleGoogle Scholar
- Auer PL, Doerge RW: Statistical Design and Analysis of RNA Sequencing Data. Genetics. 2010, 185: 405-416. 10.1534/genetics.110.114983.PubMed CentralPubMedView ArticleGoogle Scholar
- Young ND, Debelle F, Oldroyd GED, Geurts R, Cannon SB, Udvardi MK, Benedito VA, Mayer KFX, Gouzy J, Schoof H, Van de Peer Y, Proost S, Cook DR, Meyers BC, Spannagl M, Cheung F, De Mita S, Krishnakumar V, Gundlach H, Zhou S, Mudge J, Bharti AK, Murray JD, Naoumkina MA, Rosen B, Silverstein KAT, Tang H, Rombauts S, Zhao PX, Zhou P, et al: The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011, 480: 520-524.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barrie K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, et al: Genome sequence of the palaeopolyploid soybean. Nature. 2010, 463: 178-183. 10.1038/nature08670.PubMedView ArticleGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, Radenbaugh A, Singh S, Swing V, Tissier C, Zhang P, Huala E: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucl Acids Res. 2008, 36 (suppl 1): D1009-D1014.PubMed CentralPubMedGoogle Scholar
- Guo P, Baum M, Grando S, Ceccarelli S, Bai G, Li R, von Korff M, Varshney RK, Graner A, Valkoun J: Differentially expressed genes between drought-tolerant and drought-sensitive barley genotypes in response to drought stress during the reproductive stage. J Exp Bot. 2009, 60: 3531-3544. 10.1093/jxb/erp194.PubMed CentralPubMedView ArticleGoogle Scholar
- Rivero RM, Kojima M, Gepstein A, Sakakibara H, Mittler R, Gepstein S, Blumwald E: Delayed leaf senescence induces extreme drought tolerance in a flowering plant. Proc Natl Acad Sci U S A. 2007, 104: 19631-19636. 10.1073/pnas.0709453104.PubMed CentralPubMedView ArticleGoogle Scholar
- Flexas J, Medrano H: Drought-inhibition of photosynthesis in C-3 plants: Stomatal and non-stomatal limitations revisited. Ann Bot. 2002, 89: 183-189. 10.1093/aob/mcf027.PubMed CentralPubMedView ArticleGoogle Scholar
- Vernon DM, Bohnert HJ: A novel methyl transferase induced by osmotic stress in the facultative halophyte Mesembryanthemum crystallinum. EMBO J. 1992, 11: 2077-2085.PubMed CentralPubMedGoogle Scholar
- Minchin FR, Minguez MI, Sheehy JE, Witty JF, Skøt L: Relationships between nitrate and oxygen supply in symbiotic nitrogen fixation by white clover. J Exp Bot. 1986, 37: 1103-1113. 10.1093/jxb/37.8.1103.View ArticleGoogle Scholar
- Skøt L, Egsgaard H: Identification of ononitol and O-methyl-scyllo-inositol in pea root nodules. Planta. 1984, 161: 32-36. 10.1007/BF00951457.PubMedView ArticleGoogle Scholar
- Ford CW: Accumulation of low molecular weight solutes in water-stressed tropical legumes. Phytochemistry. 1984, 23: 1007-1015. 10.1016/S0031-9422(00)82601-1.View ArticleGoogle Scholar
- Herrmann D, Boller B, Studer B, Widmer F, Kölliker R: QTL analysis of seed yield components in red clover (Trifolium pratense L.). Theor Appl Genet. 2006, 112: 536-545. 10.1007/s00122-005-0158-1.PubMedView ArticleGoogle Scholar
- Smart RE, Bingham GE: Rapid Estimates of Relative Water Content. Plant Physiol. 1974, 53: 258-260. 10.1104/pp.53.2.258.PubMed CentralPubMedView ArticleGoogle Scholar
- Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Baggerly KA, Deng L, Morris JS, Aldaz CM: Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics. 2003, 19: 1477-1483. 10.1093/bioinformatics/btg173.PubMedView ArticleGoogle Scholar
- Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LL: UniProt: the Universal Protein knowledgebase. Nucl Acids Res. 2004, 32 (suppl 1): D115-D119.PubMed CentralPubMedView ArticleGoogle Scholar
- Consortium TU: Update on activities at the Universal Protein Resource (UniProt) in 2013. Nucl Acids Res. 2013, 41 (D1): D43-D47.View ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucl Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT–PCR. Nucl Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMed CentralPubMedView ArticleGoogle Scholar
- You F, Huo N, Gu Y, Luo M-c, Ma Y, Hane D, Lazo G, Dvorak J, Anderson O: BatchPrimer3: A high throughput web application for PCR and sequencing primer design. BMC Bioinform. 2008, 9: 253-10.1186/1471-2105-9-253.View ArticleGoogle Scholar
- Parveen I, Moorby JM, Fraser MD, Allison GG, Kopka J: Application of gas chromatography − mass spectrometry metabolite profiling techniques to the analysis of heathland plant diets of sheep. J Agric Food Chem. 2007, 55: 1129-1138. 10.1021/jf062995w.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.