|
|
||||||||
|
First published online July 28, 2006; 10.1104/pp.106.086256 Plant Physiology 142:135-147 (2006) © 2006 American Society of Plant Biologists OPEN ACCESS ARTICLE
The Effects of Elevated CO2 Concentration on Soybean Gene Expression. An Analysis of Growing and Mature Leaves1,[W],[OA]United States Department of Agriculture/Agricultural Research Service Photosynthesis Research Unit (E.A.A.), Department of Plant Biology (E.A.A.), and Department of Crop Sciences (A.R., L.O.V.), University of Illinois at Urbana-Champaign, Urbana, Illinois 61801; Department of Environmental Sciences, Brookhaven National Laboratory, Upton, New York 119735000 (A.R.); and ICG-III, Juelich Research Center, D52425 Juelich, Germany (E.A.A., A.W., U.S.)
Improvements in carbon assimilation and water-use efficiency lead to increases in maximum leaf area index at elevated carbon dioxide concentration ([CO2]); however, the molecular drivers for this increase are unknown. We investigated the molecular basis for changes in leaf development at elevated [CO2] using soybeans (Glycine max) grown under fully open air conditions at the Soybean Free Air CO2 Enrichment (SoyFACE) facility. The transcriptome responses of rapidly growing and fully expanded leaves to elevated [CO2] were investigated using cDNA microarrays. We identified 1,146 transcripts that showed a significant change in expression in growing versus fully expanded leaves. Transcripts for ribosomal proteins, cell cycle, and cell wall loosening, necessary for cytoplasmic growth and cell proliferation, were highly expressed in growing leaves. We further identified 139 transcripts with a significant [CO2] by development interaction. Clustering of these transcripts showed that transcripts involved in cell growth and cell proliferation were more highly expressed in growing leaves that developed at elevated [CO2] compared to growing leaves that developed at ambient [CO2]. The 327 [CO2]-responsive genes largely suggest that elevated [CO2] stimulates the respiratory breakdown of carbohydrates, which provides increased energy and biochemical precursors for leaf expansion and growth at elevated [CO2]. While increased photosynthesis and carbohydrate production at elevated [CO2] are well documented, this research demonstrates that at the transcript and metabolite level, respiratory breakdown of starch is also increased at elevated [CO2].
By 2050, soybean (Glycine max) will grow in an atmosphere with a 50% higher carbon dioxide concentration ([CO2]) (Prentice et al., 2001
Increased leaf growth, leading to larger individual leaf size, is one component of increased LAI measured at elevated [CO2] in field experiments (Taylor et al., 2003
Leaf growth is a spatially and temporally dynamic process (for review, see Schurr et al., 2006
The first objective of this research was to investigate molecular changes in growing and fully expanded soybean leaves developed at elevated [CO2] under fully open-air conditions. Research has shown that the response of soybean to elevated [CO2] in the field is less than predicted from chamber studies (Long et al., 2006
The transcriptome response of fully expanded trifoliate 4 (T4) and growing trifoliate 6 (T6) soybean leaves to elevated [CO2] was analyzed using cDNA microarrays (Fig. 1 ). On July 7, T4 leaflets were fully expanded and longer in elevated [CO2] compared to ambient [CO2] (Fig. 2 ). T6 leaflets were growing with a relative increase in length of 42% ± 6% per day in both ambient and elevated [CO2] (Fig. 2). We were specifically interested in how elevated [CO2] alters genes related to leaf development, so samples were taken between 1 and 2 AM, which corresponded to the time of maximum leaf expansion rate (Ainsworth et al., 2005
Genes Associated with Leaf Development
Expression of 1,146 genes was significantly different in growing versus fully expanded leaves, irrespective of growth [CO2] (for a complete list of transcripts, see Supplemental Table I). A total of 178 transcripts, encoding genes for a wide variety of functions, showed 1.5 times lower gene expression in growing leaves compared to fully expanded leaves. This group included genes involved in secondary metabolism, transport, stress and metal handling, and major and minor carbohydrate metabolism. Notably, starch phosphorylase (Gm-r1088-8633), a glucan-metabolizing enzyme (Smith et al., 2005
A total of 132 transcripts showed at least 1.5 times higher expression in growing leaves and, therefore, represent our best estimates of control points of leaf expansion. The major family of genes that were highly expressed in growing tissues included ribosomal proteins (Table I). Some cell-cycle genes (histones) and cell wall-loosening genes (expansins) were also included in this group (Table I). Tubulin genes, necessary for regulating the direction of diffuse growth in plants (Abe et al., 2004
We identified 139 transcripts with a CO2 x development interaction (P < 0.05). These were of particular interest because they represent potential genes involved in growth that may be altered by [CO2] treatments. These transcripts were clustered into four groups using k-means clustering (Saeed et al., 2003
The second cluster contained transcripts with lower expression in T4E compared to T6E (Fig. 4B) and rather subtle changes in other comparisons. This group contained a number of ribosomal proteins (Gm-r1070-3758, Gm-r1070-6640, Gm-r1070-8751, Gm-r1070-3694) involved in protein synthesis, as well as a binding protein (BiP; Gm-r1070-7989), a highly conserved endoplasmic-reticulum luminal protein that functions as a molecular chaperone (Kalinski et al., 1995
The third cluster included 43 transcripts that showed higher expression in T4E compared to T6E (Fig. 4C). This cluster included genes with a wide range of functions, including amino acid synthesis and transport, carbohydrate and cell wall metabolism, protein degradation, redox, and stress response (Fig. 4C). The fourth cluster included 23 genes that showed lower expression in T4A compared to both T4E and T6A (Fig. 4D), including two genes involved in cell wall metabolism, a putative NAD-dependent epimerase and a glycosyl hydrolase family 17 protein (Gm-r1070-5796 and Gm-r1070-767). Glycosyl hydrolase family 17 proteins hydrolyze 1,3-
The 327 CO2-responsive genes were assigned to different functional categories (Fig. 5 ). Many genes with roles in cellular functions (i.e. cell cycle, RNA regulation of transcription, DNA synthesis, and cell organization) showed higher expression in elevated [CO2]. Within this category, most of the genes with higher expression in elevated [CO2] were transcription factors (Table II). While increased expression of transcription factors suggests increased protein synthesis, most transcripts in the protein category (Fig. 5) were involved in protein degradation. These included ubiquitin-specific proteases, Cys proteinases, and different proteosome subunits (Table II). Therefore, we might hypothesize that growth at elevated [CO2] accelerates protein turnover. Other categories where genes were differentially expressed in elevated [CO2] included nitrogen (N) metabolism, hormone metabolism, secondary metabolism (in particular lignin biosynthesis), and transport (Fig. 5; Table II).
Growth of soybeans at elevated [CO2] stimulates photosynthesis during the day and results in marked and significant accumulations of soluble carbohydrates and starch at the end of the photoperiod (Fig. 3, A and C; Rogers et al., 2004
The transcripts of genes encoding enzymes of central C metabolism support this hypothesis. Figure 6
depicts a representation of central C metabolism, annotated with the steps where transcript levels indicated that they were up-regulated at elevated [CO2] (Heldt, 1997 While developing leaves at both ambient and elevated [CO2] had similarly high relative leaf expansion rates, mature leaves reach a larger final area at elevated [CO2] (Fig. 2). Data at the transcript level are consistent with the carbohydrate data (Fig. 3) and provide further evidence that biochemical precursors and energy from soluble carbohydrate and starch degradation may stimulate increased leaf growth and area at elevated [CO2]. Carbon from carbohydrate and starch degradation may be used along with other substrates to produce cell walls and phospholipid membranes. This was supported by the increased transcript levels of genes associated with fatty acid biosynthesis and desaturases (Fig. 6; Table II).
Soybeans get most of their N through their association with N-fixing bacteria (Ritchie et al., 1997
In this field study, we investigated the transcriptome response of soybean to elevated [CO2] in growing and fully expanded leaves. We tested the hypothesis that increased C assimilation in plants grown at elevated [CO2] altered pools of carbohydrates and transcripts that control growth and expansion of young leaves. It is well established that elevated [CO2] increases photosynthetic C fixation and carbohydrate synthesis (Long et al., 2004
Experimental Site
Soybeans (Glycine max cv 93B15; Pioneer Hi-Bred) were grown at the SoyFACE facility, located in Champaign, IL (40°02'N, 88°14'W, 228 m above sea level). SoyFACE was established on a tile-drained field that has been in continuous cultivation for more than 100 years. The 32-ha site has organically rich Flanagan/Drummer series soil. Following standard agronomic practice in the region, no fertilizer was applied. The crop was planted on May 28, 2004, and measurements were made on July 8, 2004, when the crop was in the vegetative growth phase (Ritchie et al., 1997
The length of T4 and T6 lateral leaflets was tracked with a ruler (±0.1 cm) approximately every other day from initiation of T4 until sampling of both developmental stages on July 7, 2004. Growth of 12 leaflets on six randomly selected plants per plot was followed. Leaf development in field-grown plants was similar to leaf development of plants raised in growth chambers, where a homogeneous distribution of growth along the leaf blade and a distinct diurnal rhythm of expansion were described for leaflets of a similar developmental stage (Ainsworth et al., 2005
Leaf discs from T4 and T6 middle leaflets of three plants within each plot were sampled for analysis of carbohydrates between 1 and 2 AM on July 8, 2004. Therefore, 12 leaflets per developmental stage and [CO2] treatment were sampled. Each disc (approximately 1.8 cm2) was removed from a vein-free area of a middle leaflet, wrapped in foil, and plunged immediately into liquid N. Samples were lyophilized prior to analysis.
Individual leaf discs were powdered in liquid N. Foliar contents of carbohydrates were extracted from ground leaf tissue in 80% (v/v), buffered (2 mM HEPES, pH 7.8) ethanol at 80°C. Four 20-min incubations were needed to recover the soluble carbohydrates. Glc, Fru, and Suc were determined using a continuous enzymatic substrate assay (Jones et al., 1977
T4 and T6 lateral leaflets from 12 individual soybeans within each plot were harvested between 1 and 2 AM. Entire leaflets were cut, wrapped in foil, plunged immediately into liquid N, and then lyophilized (Multi-Dry Lyophilizer; FTS Systems) and stored at 20°C. Total RNA was extracted from six pooled freeze-dried leaflets from each plot and developmental stage using a SDS/phenol chloroform method and lithium chloride precipitation (Wang and Vodkin, 1994
Spot intensities were quantified using Imagene 6.1 (Biodiscovery). The local background was subtracted for each spot, and spots were normalized to the median intensity of each dye on each slide. The natural log of the background-corrected median signal was used for all statistical analyses. Spots flagged by the Imagene image analysis software were removed from subsequent analyses (Prakash and Petrov, 2004
Biological and technical replications were averaged for each plot for statistical analysis. A mixed-model ANOVA was performed, with trifoliate and CO2 treatment as fixed effects and block as a random effect. The model was tested for conformation to the assumption of normality of the residuals using the Shapiro-Wilkes Test. A Bonferroni significance level was used as an initial criterion for rejecting the null hypothesis of a significant treatment effect (0.05/5,314 for Gm-1070 and 0.05/5,831 for Gm-1088). No genes were significant at the Bonferroni level, so we used a second nominal threshold of
We thank Steve Long and Tim Mies for management and maintenance of the SoyFACE facility. We thank L. McIntyre for assistance with statistical analysis of the experiment and insightful comments on an early draft of the manuscript. We thank R. Knepp, K. Gillespie, A.M. Boone, and S.I. Jones for technical help with RNA extractions, microarray protocols, and bioinformatics. Received July 3, 2006; accepted July 25, 2006.
1 This work was supported by the Illinois Council for Food and Agricultural Research, by the Archer Daniels Midland Company, and by the U.S. Department of Agriculture/Agricultural Research Service. E.A.A. was supported by an Alexander von Humboldt postdoctoral research fellowship. A.R. was supported by the U.S. Department of Energy Office of Science contract no. DEAC0298CH10886 to Brookhaven National Laboratory. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Elizabeth A. Ainsworth (ainswort{at}uiuc.edu).
[W] The online version of this article contains Web-only data.
[OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.106.086256 * Corresponding author; e-mail ainswort{at}uiuc.edu; fax 2172444419.
Abe T, Thitamadee S, Hashimoto T (2004) Microtubule defects and cell morphogenesis in the lefty1lefty2 tubulin mutant of Arabidopsis thaliana. Plant Cell Physiol 45: 211220 Ainsworth EA, Davey PA, Bernacchi CJ, Dermody OC, Heaton EA, Moore DJ, Morgan PB, Naidu SL, Ra HSY, Zhu XG, et al (2002) A meta-analysis of elevated [CO2] effects on soybean (Glycine max) physiology, growth and yield. Glob Change Biol 8: 695709[CrossRef] Ainsworth EA, Walter A, Schurr U (2005) Glycine max leaflets lack a base-tip gradient in growth rate. J Plant Res 118: 343346[CrossRef][ISI][Medline] Avery GS (1933) Structure and development of the tobacco leaf. Am J Bot 20: 565592[CrossRef][ISI] Bernacchi CJ, Leakey ADB, Heady LE, Morgan PB, Rogers A, Long SP, Ort DR (2006) Hourly and seasonal variation in photosynthesis and stomatal conductance of soybean grown at future CO2 and ozone concentrations for three years under fully open air conditions. Plant Cell Environ doi/10.1111/j.1365-3040.2006.01581.x Bunce JA (1977) Leaf elongation in relation to leaf water potential in soybean. J Exp Bot 28: 156161 Chatterjee M, Sparvoli S, Edmunds C, Garosi P, Finlay K, Martin C (1996) DAG, a gene required for chloroplast differentiation and palisade development in Antirrhinum majus. EMBO J 15: 41944207[ISI][Medline] Citterio S, Piatti S, Albertini E, Aina R, Varotta S, Barcaccia G (2006) Alfalfa Mob1-like proteins are involved in cell proliferation and are localized in the cell division plane during cytokinesis. Exp Cell Res 312: 10501064[CrossRef][ISI][Medline] Davey PA, Hunt S, Hymus GJ, DeLucia EH, Drake BG, Karnosky DF, Long SP (2004) Respiratory oxygen uptake is not decreased by an instantaneous elevation of [CO2], but is increased with long-term growth in the field at elevated CO2. Plant Physiol 134: 520527 Dennis DT, Blakeley SD (2000) Carbohydrate metabolism. In BB Buchanan, W Gruissem, RL Jones, eds, Biochemistry and Molecular Biology of Plants. American Society of Plant Physiologists, Rockville, MD, pp 631675 Dermody O, Long SP, Delucia EH (2006) How does elevated CO2 or ozone affect the leaf-area index of soybean when applied independently? New Phytol 169: 145155[CrossRef][ISI][Medline] Druart N, Rodríguez-Buey M, Barron-Gafford G, Sjödin A, Bhalerao R, Hurry V (2006) Molecular targets of elevated [CO2] in leaves and stems of Populus deltoides: implications for future tree growth and carbon sequestration. Funct Plant Biol 33: 121131 Ferris R, Sabatti M, Miglietta F, Mills RF, Taylor G (2001) Leaf area is stimulated in Populus by free air CO2 enrichment (POPFACE), through increased cell expansion and production. Plant Cell Environ 24: 305315[CrossRef] Gupta P, Duplessis S, White H, Karnosky DF, Martin F, Podila GK (2005) Gene expression patterns of trembling aspen trees following long-term exposure to interacting elevated CO2 and tropospheric O3. New Phytol 167: 129142[CrossRef][ISI][Medline] Heldt HW (1997) Plant Biochemistry and Molecular Biology. Oxford University Press, New York Hendriks JHM, Kolbe A, Gibon Y, Stitt M, Geigenberger P (2003) ADP-glucosepyro-phosphorylase is activated by posttranslational redox-modification in response to light and to sugars in leaves of Arabidopsis and other plant species. Plant Physiol 133: 838849 Kalinski A, Rowley DL, Loer DS, Foley C, Buta G, Herman EM (1995) Binding-protein expression is similar to temporal, developmental and stress-induced regulation in terminally differentiated soybean organs. Planta 195: 611621[ISI][Medline] Kerr MK, Churchill GA (2001) Statistical design and the analysis of gene expression microarray data. Genome Res 77: 123128 Jones MGK, Outlaw WH, Lowry OH (1977) Enzymic assay of 107 to 1014 moles of sucrose in plant tissues. Plant Physiol 60: 379383 Li H, Singh AK, McIntyre LM, Sherman LA (2004) Differential gene expression in response to hydrogen peroxide and the putative PerR regulon of Synechocystis sp. strain PCC 6803. J Bacteriol 186: 33313345 Long SP, Ainsworth EA, Leakey ADB, Nosberger J, Ort DR (2006) Food for thought: lower-than-expected crop yield stimulation with rising CO2 concentrations. Science 312: 19181921 Long SP, Ainsworth EA, Rogers A, Ort DR (2004) Rising atmospheric carbon dioxide: plants FACE the future. Annu Rev Plant Biol 55: 591628[CrossRef][Medline] Matsubara S, Hurry V, Druart N, Benedict C, Janzik I, Chavarría-Krauser A, Walter A, Schurr U (2005) Nocturnal changes in leaf growth of Populus deltoides are controlled by cytoplasmic growth. Planta 223: 13151328[Medline] Miglietta F, Peressotti A, Vaccari FP, Zaldei A, de Angelis P, Scarascia-Mugnozza G (2001) Free-air CO2 enrichment (FACE) of a poplar plantation: the POPFACE fumigation system. New Phytol 150: 465476[CrossRef] Morgan PB, Bollero GA, Nelson RL, Dohleman FG, Long SP (2005) Smaller than predicted increase in aboveground net primary production and yield of field-grown soybean under fully open-air [CO2] elevation. Glob Change Biol 11: 18561865[CrossRef] Ort DR, Ainsworth EA, Aldea M, Allen DJ, Bernacchi CJ, Berenbaum MR, Bollero GA, Cornic G, Davey PA, Dermody O, et al (2006) SoyFACE: the effects and interactions of elevated [CO2] and [O3] on soybean. In J Nösberger, SP Long, RJ Norby, M Stitt, GR Hendrey, H Blum, eds, Managed Ecosystems and CO2: Case Studies, Processes and Perspectives. Springer, Berlin, pp 7185 Prakash PJ, Petrov A (2004) Gene Flagging in ImaGene. Biodiscover Inc. Technical Bulletin. Biodiscover Inc., El Segundo, CA, pp 17 Prentice IC, Farquhar GD, Fasham MJR, Goulden ML, Heimann M, Jaramillo VJ, Kheshgi HS, LeQuere C, Scholes RJ, Wallace DWR, et al (2001) The carbon cycle atmospheric carbon dioxide. In JT Houghton, Y Ding, DJ Griggs, M Noguer, PJ Van der Linder, X Dai, K Maskell, CA Johnson, eds, Climate Change 2001: The Scientific Basis. Contributions of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, New York, pp 183239 Raines CA, Paul MJ (2006) Products of leaf primary carbon metabolism modulate the developmental programme determining plant morphology. J Exp Bot 57: 18571862 Ritchie SW, Hanaway JJ, Thompson HE, Benson GO (1997) How a Soybean Plant Develops. Special Report Number 53. Iowa State University, Ames, IA Rogers A, Allen DJ, Davey PA, Morgan PB, Ainsworth EA, Bernacchi CJ, Cornic G, Dermody O, Dohleman FG, Heaton EA, et al (2004) Leaf photosynthesis and carbohydrate dynamics of soybeans grown throughout their life-cycle under Free-Air Carbon dioxide Enrichment. Plant Cell Environ 27: 449458[Medline] Rogers A, Gibon Y, Stitt M, Morgan PB, Bernacchi CJ, Ort DR, Long SP (2006) Increased carbon availability at elevated carbon dioxide concentration improves N assimilation in a legume. Plant Cell Environ 29: 16511658[CrossRef][Medline] Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier TM, Sturn A, et al (2003) TM4: a free, open-source system for microarray data management analysis. Biotechniques 34: 374378[ISI][Medline] Sambrook J, Fritsch EF, Maniatis T (1989) Molecular Cloning: A Laboratory Manual, Ed 2. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY Schurr U, Walter A, Rascher U (2006) Functional dynamics of plant growth and photosynthesisfrom steady-state to dynamicsfrom homogeneity to heterogeneity. Plant Cell Environ 29: 340352[CrossRef][Medline] Smith AM, Zeeman SC, Smith SM (2005) Starch degradation. Annu Rev Plant Biol 56: 7398[CrossRef][Medline] Sugimoto-Shirasu K, Roberts K (2003) "Big it up": endoreduplication and cell-size control in plants. Curr Opin Plant Biol 6: 544553[CrossRef][ISI][Medline] Taylor G, Ranasinghe S, Bosac C, Gardner SDL, Ferris R (1994) Elevated CO2 and plant growth: cellular mechanisms and responses of whole plants. J Exp Bot 45: 17611774 Taylor G, Street NR, Tricker PJ, Sjödin A, Graham L, Skogström O, Calfapietra C, Scarascia-Mugnozza G, Jansson S (2005) The transcriptome of Populus in elevated CO2. New Phytol 167: 143154[Medline] Taylor G, Tricker PJ, Zhang FZ, Alston VJ, Miglietta F, Kuzminsky E (2003) Spatial and temporal effects of free-air CO2 enrichment (POPFACE) on leaf growth, cell expansion, and cell production in a closed canopy of poplar. Plant Physiol 131: 177185 Thomas BR, Romero GO, Nevins DJ, Rodriguez RL (2000) New perspectives on the endo-beta-glucanases of glycosyl hydrolase Family 17. Int J Biol Macromol 27: 139144[Medline] Trainotti L, Pavanello A, Casadoro G (2004) Differential expression of genes in apical and basal tissues of expanding tobacco leaves. Plant Sci 167: 679686 Tricker PJ, Calfapietra C, Kuzminsky E, Puleggi R, Ferris R, Nathoo M, Pleasants LJ, Alston V, de Angelis P, Taylor G (2004) Long-term acclimation of leaf production, development, longevity and quality following 3 yr exposure to free-air CO2 enrichment during canopy closure in Populus. New Phytol 162: 413426[CrossRef] Turgeon R (1989) The sink-source transition in leaves. Annu Rev Plant Physiol Plant Mol Biol 40: 119138[CrossRef][ISI] Vodkin LO, Khanna A, Shealy R, Clough SJ, Gonzalez DO, Philip R, Zabala G, Thibaud-Nissen F, Sidarous M, Stromvik M, et al (2004) Microarrays for global expression constructed with a low redundancy set of 27,500 sequenced cDNAs representing an array of developmental stages and physiological conditions of the soybean plant. BMC Genomics 5: 73[CrossRef][Medline] Walter A, Christ MM, Barron-Gifford GA, Grieve KA, Murthy R, Rascher U (2005) The effect of elevated CO2 on diel leaf growth cycle, leaf carbohydrate content and canopy growth performance of Populus deltoides. Glob Change Biol 11: 12581271 Wang C-S, Vodkin LO (1994) Extraction of RNA from tissues containing high levels of procyanidins that bind RNA. Plant Mol Biol Rep 12: 132145 Wayne ML, McIntyre LM (2002) Combining mapping and arraying: an approach to candidate gene identification. Proc Natl Acad Sci USA 99: 1490314906 Related articles in Plant Physiol.:
This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ASPB Publications | PLANT PHYSIOLOGY | THE PLANT CELL | |
|---|---|---|---|