De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

BMC Genomics

Leif Skøt⁵  Glyn Jenkins⁴  Michael T Abberton³  Tom Ruttink²  Gordon G Allison⁵  Matthew Lowe⁵  Igor Chernukin¹  Matthew J Hegarty⁵  Martin T Swain⁴  Steven A Yates¹ 

[1] Present address: School of Biological Sciences, University of Essex, Wivenhoe Park, Colchester, Essex CO4 3SQ, UK;Institute for Agricultural and Fisheries Research (ILVO), Plant Sciences Unit – Growth and Development, Caritasstraat 21, 9090 Melle, Belgium;Present address: International Institute of Tropical Agriculture (IITA), PMB 5320, Oyo Road, Ibadan, Nigeria;Institute of Biological, Environmental and Rural Sciences, Aberystwyth University, Penglais, Aberystwyth, Ceredigion SY23 3FL, UK;Institute of Biological, Environmental and Rural Sciences, Aberystwyth University, Gogerddan, Aberystwyth, Ceredigion SY23 3 EB, UK

关键词: Trifolium pratense;    Transcriptome assembly;    RNA-Seq;    Red clover;    Polymorphism;    Drought stress;   

Others  :  1216628
DOI  :  10.1186/1471-2164-15-453

 received in 2013-11-29, accepted in 2014-06-04,  发布年份 2014

【 摘 要 】

Background

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.

Results

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.

Conclusions

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.

【 授权许可】

   
2014 Yates et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

【 参考文献 】

• [1]Taylor NL, Quesenberry KH: Red Clover Science. Dordrecht: Kluwer Academic Publishers; 1996.
• [2]Abberton MT, Marshall AH: Progress in breeding perennial clovers for temperate agriculture. J Agric Sci 2005, 143:117-135.
• [3]Frame J, Charlton JFL, Laidlaw AS: Temperate Forage Legumes. Wallingford: CAB International; 1998.
• [4]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.
• [5]Taylor NL: A century of clover breeding developments in the United States. Crop Sci 2008, 48:1-13.
• [6]Thomashow MF: molecular basis of plant cold acclimation: insights gained from studying the CBF cold response pathway. Plant Physiol 2010, 154:571-577.
• [7]Thomashow MF: Plant Cold Acclimation: Freezing Tolerance Genes and Regulatory Mechanisms. Annu Rev Plant Physiol Plant Mol Biol 1999, 50:571-599.
• [8]Thomashow MF: Role of Cold-Responsive Genes in Plant Freezing Tolerance. Plant Physiol 1998, 118:1-8.
• [9]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]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.
• [11]Peterson PR, Sheaffer CC, Hall MH: Drought effects on perennial forage legume yield and quality. Agron J 1992, 84:6.
• [12]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.
• [13]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.
• [14]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.
• [15]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. BioMed Central Full Text
• [16]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.
• [17]Shinozaki K, Yamaguchi-Shinozaki K: Gene networks involved in drought stress response and tolerance. J Exp Bot 2007, 58:221-227.
• [18]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.
• [19]Chaves MM, Maroco JP, Pereira JS: Understanding plant responses to drought - from genes to the whole plant. Funct Plant Biol 2003, 30:239-264.
• [20]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.
• [21]Agarwal PK, Jha B: Transcription factors in plants and ABA dependent and independent abiotic stress signalling. Biol Plant 2010, 54:201-212.
• [22]Mittler R, Vanderauwera S, Gollery M, Van Breusegem F: Reactive oxygen gene network of plants. Trends Plant Sci 2004, 9:490-498.
• [23]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.
• [24]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.
• [25]Seki M, Umezawa T, Urano K, Shinozaki K: Regulatory metabolic networks in drought stress responses. Curr Opin Plant Biol 2007, 10:296-302.
• [26]Ashraf M, Foolad MR: Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot 2007, 59:206-216.
• [27]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.
• [28]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. BioMed Central Full Text
• [29]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.
• [30]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. BioMed Central Full Text
• [31]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.
• [32]Wang L, Li P, Brutnell TP: Exploring plant transcriptomes using ultra high-throughput sequencing. Briefings Funct Genom 2010, 9:118-128.
• [33]Kircher M, Kelso : High-throughput DNA sequencing – concepts and limitations. Bioessays 2010, 32:524-536.
• [34]Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics 2008, 92(5):255-264.
• [35]Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol 2008, 26:1135-1145.
• [36]Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth 2008, 5:621-628.
• [37]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.
• [38]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. BioMed Central Full Text
• [39]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.
• [40]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.
• [41]Torii KU: Leucine-Rich Repeat Receptor Kinases in Plants: Structure, Function, and Signal Transduction Pathways. 2004, 234.
• [42]Kanai R, Edwards GE: The biochemistry of C4 photosynthesis. In C4 Plant Biology. Edited by Sage RF, Monson RK. San Diego: Academic; 1999:49-87.
• [43]Hatch MD: C4 photosynthesis: A historical perspective. In C4 Plant Biology. Edited by Sage RF, Monson RK. San Diego: Academic; 1999:17-46.
• [44]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.
• [45]Auer PL, Doerge RW: Statistical Design and Analysis of RNA Sequencing Data. Genetics 2010, 185:405-416.
• [46]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.
• [47]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.
• [48]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.
• [49]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.
• [50]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.
• [51]Flexas J, Medrano H: Drought-inhibition of photosynthesis in C-3 plants: Stomatal and non-stomatal limitations revisited. Ann Bot 2002, 89:183-189.
• [52]Vernon DM, Bohnert HJ: A novel methyl transferase induced by osmotic stress in the facultative halophyte Mesembryanthemum crystallinum. EMBO J 1992, 11:2077-2085.
• [53]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.
• [54]Skøt L, Egsgaard H: Identification of ononitol and O-methyl-scyllo-inositol in pea root nodules. Planta 1984, 161:32-36.
• [55]Ford CW: Accumulation of low molecular weight solutes in water-stressed tropical legumes. Phytochemistry 1984, 23:1007-1015.
• [56]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.
• [57]Smart RE, Bingham GE: Rapid Estimates of Relative Water Content. Plant Physiol 1974, 53:258-260.
• [58]Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 2008, 18:821-829.
• [59]Baggerly KA, Deng L, Morris JS, Aldaz CM: Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics 2003, 19:1477-1483.
• [60]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.
• [61]Consortium TU: Update on activities at the Universal Protein Resource (UniProt) in 2013. Nucl Acids Res 2013, 41(D1):D43-D47.
• [62]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.
• [63]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.
• [64]Pfaffl MW: A new mathematical model for relative quantification in real-time RT–PCR. Nucl Acids Res 2001, 29:e45.
• [65]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. BioMed Central Full Text
• [66]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.