Daytime soybean transcriptome fluctuations during water deficit stress

BMC Genomics

Alexandre Lima Nepomuceno⁴  Frank G Harmon¹  Francisco Pereira Lobo⁵  Hugo Bruno Correa Molinari⁴  Thiago Jonas Nakayama²  Juliana Marcolino-Gomes³  Renata Fuganti-Pagliarini³  Fabiana Aparecida Rodrigues³ 

[1] Department of Plant and Microbial Biology, UC Berkeley, Berkeley, CA, USA;Department of Crop Science, Federal University of Viçosa, Viçosa, MG, Brazil;Brazilian Agricultural Research Corporation- Embrapa Soybean, Embrapa Soybean- Rod. Carlos João Strass, s/n, Londrina 86001-970, PR, Brazil;Embrapa LABEX US Plant Biotechnology at ARS/USDA Plant Gene Expression Center, Albany, CA, USA;Brazilian Agricultural Research Corporation-Embrapa Agricultural Informatics, Campinas, SP, Brazil

关键词: Plant metabolism;    Glycine max;    Drought;    diel regulation;    Daily oscillation;    Abiotic stress;   

Others  :  1219119
DOI  :  10.1186/s12864-015-1731-x

 received in 2015-01-07, accepted in 2015-06-26,  发布年份 2015

【 摘 要 】

Background

Since drought can seriously affect plant growth and development and little is known about how the oscillations of gene expression during the drought stress-acclimation response in soybean is affected, we applied Illumina technology to sequence 36 cDNA libraries synthesized from control and drought-stressed soybean plants to verify the dynamic changes in gene expression during a 24-h time course. Cycling variables were measured from the expression data to determine the putative circadian rhythm regulation of gene expression.

Results

We identified 4866 genes differentially expressed in soybean plants in response to water deficit. Of these genes, 3715 were differentially expressed during the light period, from which approximately 9.55 % were observed in both light and darkness. We found 887 genes that were either up- or down-regulated in different periods of the day. Of 54,175 predicted soybean genes, 35.52 % exhibited expression oscillations in a 24 h period. This number increased to 39.23 % when plants were submitted to water deficit. Major differences in gene expression were observed in the control plants from late day (ZT16) until predawn (ZT20) periods, indicating that gene expression oscillates during the course of 24 h in normal development. Under water deficit, dissimilarity increased in all time-periods, indicating that the applied stress influenced gene expression. Such differences in plants under stress were primarily observed in ZT0 (early morning) to ZT8 (late day) and also from ZT4 to ZT12. Stress-related pathways were triggered in response to water deficit primarily during midday, when more genes were up-regulated compared to early morning. Additionally, genes known to be involved in secondary metabolism and hormone signaling were also expressed in the dark period.

Conclusions

Gene expression networks can be dynamically shaped to acclimate plant metabolism under environmental stressful conditions. We have identified putative cycling genes that are expressed in soybean leaves under normal developmental conditions and genes whose expression oscillates under conditions of water deficit. These results suggest that time of day, as well as light and temperature oscillations that occur considerably affect the regulation of water deficit stress response in soybean plants.

【 授权许可】

   
2015 Rodrigues et al.

【 图 表 】

Fig. 1.

Figure 11.

Fig. 3.

Fig. 4.

Fig. 5.

Fig. 6.

Fig. 7.

Fig. 8.

Fig. 9.

【 参考文献 】

• [1]Pedersen P, Kumudini S, Board J, Conley S. Soybean growth and development. In: Using foliar fungicides to manage soybean rust. Dorrance AE, Draper MA, Hershman DE, editors. Ohio State University, Ohio; 2005: p.41-7.
• [2]Xoconostle-Cázares B, Ramírez-Ortega FA, Flores-Elenes L, Ruiz-Medrano R. Drought tolerance in crop plants. Am J Plant Physiol. 2010; 5:241-56.
• [3]Osakabe Y, Yamaguchi-Shinozaki K, Shinozaki K, Tran L-SP. Sensing the environment: key roles of membrane-localized kinases in plant perception and response to abiotic stress. J Exp Bot. 2013; 64:445-58.
• [4]James AB, Syed NH, Bordage S, Marshall J, Nimmo GA, Jenkins GI, Herzyk P, Brown JWS, Nimmo HG. Alternative splicing mediates responses of the Arabidopsis circadian clock to temperature changes. Plant Cell. 2012; 24:961-81.
• [5]Khan S, Rowe SC, Harmon FG. Coordination of the maize transcriptome by a conserved circadian clock. BMC Plant Biol. 2010; 10:126. BioMed Central Full Text
• [6]Legnaioli T, Cuevas J, Mas P. TOC1 functions as a molecular switch connecting the circadian clock with plant responses to drought. EMBO J. 2009; 28:3745-57.
• [7]Covington MF, Maloof JN, Straume M, Kay SA, Harmer SL. Global transcriptome analysis reveals circadian regulation of key pathways in plant growth and development. Genome Biol. 2008; 9:R130. BioMed Central Full Text
• [8]Harmer SL, Hogenesch JB, Straume M, Chang H, Han B, Zhu T, Wang X. Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science. 2014; 290:2110-3.
• [9]Schaffer R, Landgraf J, Accerbi M, Simon V, Larson M, Wisman E. Microarray analysis of diurnal and circadian-regulated genes in Arabidopsis. Plant Cell. 2001; 13:113-23.
• [10]Pruneda-Paz JL, Kay SA. An expanding universe of circadian networks in higher plants. Trends Plant Sci. 2010; 15:259-65.
• [11]Covington MF, Harmer SL. The circadian clock regulates auxin signaling and responses in Arabidopsis. PLoS Biol. 2007; 5:1773-84.
• [12]Mizuno T, Yamashino T. Comparative transcriptome of diurnally oscillating genes and hormone-responsive genes in Arabidopsis thaliana: insight into circadian clock-controlled daily responses to common ambient stresses in plants. Plant Cell Physiol. 2008; 49:481-7.
• [13]Nováková M, Motyka V, Dobrev PI, Malbeck J, Gaudinová A, Vanková R. Diurnal variation of cytokinin, auxin, and abscisic acid levels in tobacco leaves. J Exp Bot. 2005; 56:2877-83.
• [14]Carbonell-Bejerano P, Rodríguez V, Royo C, Hernáiz S, Moro-González LC, Torres-Viñals M, Martínez-Zapater JM. Circadian oscillatory transcriptional programs in grapevine ripening fruits. BMC Plant Biol. 2014; 14:78. BioMed Central Full Text
• [15]Bieniawska Z, Espinoza C, Schlereth A, Sulpice R, Hincha DK, Hannah MA. Disruption of the Arabidopsis circadian clock is responsible for extensive variation in the cold-responsive transcriptome. Plant Physiol. 2008; 147:263-79.
• [16]Wilkins O, Waldron L, Nahal H, Provart NJ, Campbell MM. Genotype and time of day shape the Populus drought response. Plant J. 2009; 60:703-15.
• [17]Wilkins O, Bräutigam K, Campbell MM. Time of day shapes Arabidopsis drought transcriptomes. Plant J. 2010; 63:715-27.
• [18]Hudson KA. The circadian clock-controlled transcriptome of developing soybean seeds. Plant Genome J. 2010; 3:3.
• [19]He C-Y, Zhang J-S, Chen S-Y. A soybean gene encoding a proline-rich protein is regulated by salicylic acid, an endogenous circadian rhythm, and by various stresses. Theor Appl Genet. 2002; 104:1125-31.
• [20]Marcolino-Gomes J, Rodrigues FA, Fuganti-Pagliarini R, Bendix C, Nakayama TJ, Celaya B, Molinari HBC, Oliveira MCN, Harmon FG, Nepomuceno A. Diurnal oscillations of soybean circadian clock and drought responsive genes. PLoS One. 2014; 9:e86402.
• [21]Oya T, Nepomuceno AL, Neumaier N, Farias JRB, Tobita S, Ito O. Drought tolerance characteristics of Brazilian soybean cultivars: evaluation and characterization of drought tolerance of various Brazilian soybean cultivars in the field. Plant Prod Sci. 2004; 7:129-37.
• [22]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, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang X et al.. Genome sequence of the palaeopolyploid soybean. Nature. 2010; 463:178-83.
• [23]Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011; 6:e21800.
• [24]Hsu PY, Harmer SL. Circadian phase has profound effects on differential expression analysis. PLoS One. 2012; 7:e49853.
• [25]Raeini-Sarjaz M. Circadian rhythm leaf movement of Phaseolus vulgaris and the role of calcium ions. Plant Signal Behav. 2011; 6:962-7.
• [26]Frey A, Effroy D, Lefebvre V, Seo M, Perreau F, Berger A, Sechet J, To A, North HM, Marion-Poll A. Epoxycarotenoid cleavage by NCED5 fine-tunes ABA accumulation and affects seed dormancy and drought tolerance with other NCED family members. Plant J. 2012; 70:501-12.
• [27]Blanvillain R, Kim JH, Wu S, Lima A, Ow DW. OXIDATIVE STRESS 3 is a chromatin-associated factor involved in tolerance to heavy metals and oxidative stress. Plant J. 2009; 57:654-65.
• [28]de Carvalho MHC. Drought stress and reactive oxygen species. Plant Signal Behav. 2008; 3:156-65.
• [29]Borsani O, Díaz P, Agius MF, Valpuesta V, Monza J. Water stress generates an oxidative stress through the induction of a specific Cu/Zn superoxide dismutase in Lotus corniculatus leaves. Plant Sci. 2001; 161:757-63.
• [30]Fraire-Velázquez S, Balderas-Hernández V. Abiotic stress in plants and metabolic responses. In: Vahdati K, Leslie C. Rijeka. Abiotic stress: plant responses and applications in agriculture. InTech is registered [(MBS): 040244148; (OIB): 71692491655]; Rijeka, Croatia. 2013: 25-48.
• [31]Dhaubhadel S, Browning KS, Gallie DR, Krishna P. Brassinosteroid functions to protect the translational machinery and heat-shock protein synthesis following thermal stress. Plant J. 2002; 29:681-91.
• [32]Dhaubhadel S, Chaudhary S, Dobinson KF, Krishna P. Treatment with 24-epibrassinolide, a brassinosteroid, increases the basic thermotolerance of Brassica napus and tomato seedlings. Plant Mol Biol. 1999; 40:333-42.
• [33]Manosalva PM, Davidson RM, Liu B, Zhu X, Hulbert SH, Leung H, Leach JE. A germin-like protein gene family functions as a complex quantitative trait locus conferring broad-spectrum disease resistance in rice. Plant Physiol. 2009; 149:286-96.
• [34]Wang T, Chen X, Zhu F, Li H, Li L, Yang Q, Chi X, Yu S, Liang X. Characterization of peanut germin-like proteins, AhGLPs, in plant development and defense. PLoS One. 2013; 8:e61722.
• [35]Zimmermann G, Bäumlein H, Mock H-P, Himmelbach A, Schweizer P. The multigene family encoding germin-like proteins of barley: regulation and function in basal host resistance. Plant Physiol. 2006; 142:181-92.
• [36]Lu M, Han Y-P, Gao J-G, Wang X-J, Li W-B. Identification and analysis of the germin-like gene family in soybean. BMC Genomics. 2010; 11:620. BioMed Central Full Text
• [37]Membré N, Bernier F, Staiger D, Berna A. Arabidopsis thaliana germin-like proteins: common and specific features point to a variety of functions. Planta. 2000; 211:345-54.
• [38]Staiger D, Apel K, Trepp G. The Atger3 promoter confers circadian clock-regulated transcription with peak expression at the beginning of the night. Plant Mol Biol. 1999; 40:873-82.
• [39]Heintzen C, Fischer R, Melzer S, Kappeler S, Apel K, Staiger D. Circadian oscillations of a transcript encoding a germin-like protein that is associated with cell walls in young leaves of the long-day plant Sinapis alba. Plant Physiol. 1994; 106:905-15.
• [40]Carter C, Graham RA, Thornburg RW. Arabidopsis thaliana contains a large family of germin-like proteins: characterization of cDNA and genomic sequences encoding 12 unique family members. Plant Mol Biol. 1998; 38:929-43.
• [41]Hattori J, Boutilier KA, Campagne MML, Miki BL. A conserved BURP domain defines a novel group of plant proteins with unusual primary structures. Mol Gen Genet. 1998; 259:424-8.
• [42]Van Son L, Tiedemann J, Rutten T, Hillmer S, Hinz G, Zank T, Manteuffel R, Bäumlein H. The BURP domain protein AtUSPL1 of Arabidopsis thaliana is destined to the protein storage vacuoles, and overexpression of the cognate gene distorts seed development. Plant Mol Biol. 2009; 71:319-29.
• [43]Ding X, Hou X, Xie K, Xiong L. Genome-wide identification of BURP domain-containing genes in rice reveals a gene family with diverse structures and responses to abiotic stresses. Planta. 2009; 230:149-63.
• [44]Xu H, Li Y, Yan Y, Wang K, Gao Y, Hu Y. Genome-scale identification of soybean BURP-domain-containing genes and their expression under stress treatments. BMC Plant Biol. 2010; 10:197. BioMed Central Full Text
• [45]Wang H, Zhou L, Fu Y, Cheung M-Y, Wong F-L, Phang T-H, Sun Z, Lam H-M. Expression of an apoplast-localized BURP-domain protein from soybean (GmRD22) enhances tolerance toward abiotic stress. Plant Cell Environ. 2012; 35:1932-47.
• [46]Kant P, Gordon M, Kant S, Zolla G, Davydov O, Heimer YM, Chalifa-Caspi V, Shaked R, Barak S. Functional-genomics-based identification of genes that regulate Arabidopsis responses to multiple abiotic stresses. Plant Cell Environ. 2008; 31:697-714.
• [47]Urano K, Kurihara Y, Seki M, Shinozaki K. “Omics” analyses of regulatory networks in plant abiotic stress responses. Curr Opin Plant Biol. 2010; 13:132-8.
• [48]Spoel S, van Ooijen G. Circadian redox signaling in plant immunity and abiotic stress. Antioxid Redox Sign. 2013; 20:1-46.
• [49]Fehr WR, Caviness C, Burmood D, Pernnigton J. Stage of development description for soybeans (Glycine max L, Merrill). Crop Sci. 1971; 11:929-31.
• [50]Shukla A, Panchal H, Mishra M, Patel PR, Srivastava HS, Patel P, Shukla AK. Soil moisture estimation using gravimetric technique and FDR probe technique: a comparative analysis. Am Int J Res Formal, Appl Nat Sci. 2014; 8:89-92.
• [51]Tariq MA, Kim HJ, Jejelowo O, Pourmand N. Whole-transcriptome RNAseq analysis from minute amount of total RNA. Nucleic Acids Res. 2011; 39:e120.
• [52]Li H, Durbin R. Fast and accurate short-read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754-60.
• [53]Robinson MD, Mccarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26:139-40.
• [54]Benjamini Y, Rochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Methodol. 1995; 57:289-300.
• [55]Du Z, Zhou X, Ling Y, Zhang Z, Su Z. AgriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010; 38:64-70.
• [56]Hughes ME, Hogenesch JB, Kornacker K. JTK_CYCLE: an efficient nonparametric algorithm for detecting rhythmic components in genome-scale datasets. J Biol Rhythm. 2010; 25:372-80.
• [57]Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001; 25:402-8.
• [58]Jian B, Liu B, Bi Y, Hou W, Wu C, Han T. Validation of internal control for gene expression study in soybean by quantitative real-time PCR. BMC Mol Biol. 2008; 9:59. BioMed Central Full Text
• [59]Stolf-Moreira R, Lemos EGM, Abdelnoor RV, Beneventi MA, Rolla AAP, Pereira SS, Oliveira MCN, Nepomuceno AL, Marcelino-Guimarães FC. Identification of reference genes for expression analysis by real time. Pesqui Agropec Bras. 2011; 46:58-65.