ORIGINAL RESEARCH article

Front. Plant Sci., 19 September 2017

Sec. Plant Breeding

Volume 8 - 2017 | https://doi.org/10.3389/fpls.2017.01506

Metabolic Pathways Involved in Carbon Dioxide Enhanced Heat Tolerance in Bermudagrass

  • 1. College of Agro-grassland Science, Nanjing Agricultural University Nanjing, China

  • 2. Department of Plant Biology and Pathology, Rutgers, The State University of New Jersey, New Brunswick NJ, United States

Abstract

Global climate changes involve elevated temperature and CO2 concentration, imposing significant impact on plant growth of various plant species. Elevated temperature exacerbates heat damages, but elevated CO2 has positive effects on promoting plant growth and heat tolerance. The objective of this study was to identify metabolic pathways affected by elevated CO2 conferring the improvement of heat tolerance in a C4 perennial grass species, bermudagrass (Cynodon dactylon Pers.). Plants were planted under either ambient CO2 concentration (400 μmol⋅mol-1) or elevated CO2 concentration (800 μmol⋅mol-1) and subjected to ambient temperature (30/25°C, day/night) or heat stress (45/40°C, day/night). Elevated CO2 concentration suppressed heat-induced damages and improved heat tolerance in bermudagrass. The enhanced heat tolerance under elevated CO2 was attributed to some important metabolic pathways during which proteins and metabolites were up-regulated, including light reaction (ATP synthase subunit and photosystem I reaction center subunit) and carbon fixation [(glyceraldehyde-3-phosphate dehydrogenase, GAPDH), fructose-bisphosphate aldolase, phosphoglycerate kinase, sedoheptulose-1,7-bisphosphatase and sugars) of photosynthesis, glycolysis (GAPDH, glucose, fructose, and galactose) and TCA cycle (pyruvic acid, malic acid and malate dehydrogenase) of respiration, amino acid metabolism (aspartic acid, methionine, threonine, isoleucine, lysine, valine, alanine, and isoleucine) as well as the GABA shunt (GABA, glutamic acid, alanine, proline and 5-oxoproline). The up-regulation of those metabolic processes by elevated CO2 could at least partially contribute to the improvement of heat tolerance in perennial grass species.

Introduction

Global climate changes involve elevated temperature and CO2 concentration, imposing significant impact on plant growth (). During this century, global temperatures are predicted to rise by 2–5°C; atmospheric CO2 concentration has increased by 100 μmol mol-1 since the beginning of the industrialized era and the concentration is predicted to continue rising at a rate of approximately 2 μmol mol-1 per year (). Previous research has shown that elevated CO2 promotes plant growth under optimal growing temperatures in various plant species (; ; ). Recent research also found that elevated CO2 has positive effects on promoting heat tolerance in terms of vegetative growth in C3 species, such as rice (Oryza sativa) (; ; ), wheat (Triticum aestivum) (; ), and cool-season perennial grass species (, ) and C4 plant species, such as Bouteloua gracilis (), peanut (Arachis hypogaea) (), grain sorghum (Sorghum bicolor) () and maize (Zea mays) (). The mechanisms regulating elevated CO2 effects on C3 plant species have been reported, which have been associated with enhanced cellular expansion and cell division resulted from increased carbohydrate availability and changes in proteins and gene transcript levels (; ; ; ). However, metabolic factors underlying elevated CO2 improvement of heat tolerance in C4 perennial grass species are not well understood.

Metabolic and proteomic analysis mostly in C3 plant species demonstrated that elevated CO2 causes changes in various metabolic processes or pathways such as photosynthetic carbon fixation, respiratory metabolism, cellular growth, and stress defense (; , ; , ). The improved heat tolerance by doubling ambient CO2 concentration in C3 grass species, such as tall fescue (Festuca arundinacea), has been attributed to increases in the accumulation of metabolites, such as organic acids (shikimic acid, malonic acid, glyceric acid, threonic acid, galactaric acid, and citric acid), sugars (sucrose and maltose) and amino acids (valine, serine, and 5-oxoproline) involved in photosynthesis, respiration and amino acid metabolism (). In addition, doubling ambient CO2 concentration significantly increased the accumulation of soluble leaf carbohydrates and activity of adenosine-5′-diphosphoglucose pyrophosphorylase under high temperature in kidney bean (Phaseolus vulgaris) (). Proteomic profiling of tall fescue exposed to elevated CO2 concentration under heat stress found increased abundance of proteins associated with functions of photosynthetic light reaction, electron transport carrier molecule, ATP generation enzyme and antioxidant system (). It has been reported that C4 plant species are generally less responsive to elevated CO2 than C3 species when they are exposed to their respective optimal temperature conditions (; ). Mechanisms of elevated CO2-induced stimulation of photosynthesis in C3 plants were mainly associated with changes in electron transport during in light reaction as well as capacity for carbon fixation and assimilation during dark respiration (, ; ). However, the key changes in metabolites and proteins induced by elevated CO2 in C4 plants under heat stress have not yet to be determined.

The objective of the current study was to identify metabolic pathways affected by elevated CO2 conferring the improvement of heat tolerance in a C4 perennial grass species, bermudagrass (Cynodon dactylon) widely used as forage and turfgrass species. Understanding changes of metabolites and proteins in C4 species in response to elevated CO2 concentration will provide new insights to mechanisms about elevated CO2-mitigated effects on heat stress.

Materials and Methods

Plant Materials and Growth Conditions

Stolons of bermudagrass (cv. ‘Tifway’) plants were collected from the research farm at Nanjing Agricultural University in Nanjing, China, and transplanted into pots (20 cm in diameter and 20 cm in depth) filled with a mixture of soil and sand (soil: sand = 1:1, v/v). Plants were grown in a greenhouse with average temperature of 30/22°C (day/night), natural sunlight and irrigated once a week with half-strength Hoagland’s nutrient solution () to establish canopy and roots for 2 months. During this period, plants were trimmed once a week to keep a canopy height of 4–5 cm. After establishment, plants were transferred to growth chambers (Xubang, Jinan, Shandong province, China) with the temperature of 30/25°C (day/night), 70% relative humidity, photosynthetically active radiation of 650 μmol⋅m-2⋅s-1 and a 12-h photoperiod.

Experimental Design and Treatments

The CO2 concentrations set-up and control in growth chambers followed the same designed as described in ,). Each CO2 treatment was imposed in four growth chambers on September 1, 2015. In order to evaluate the long-term effects of elevated CO2, plants were grown under the two CO2 concentrations for 70 days prior to the exposure to heat stress. Plants grown under either CO2 treatment was then exposed to (45/40°C) (heat stress) or 30/25°C (non-stress control) in two growth chambers on November 9, 2015 until December 8, 2015. Plants were randomly relocated in each chamber twice per week to avoid confounding effects of environmental variation between different chambers. The CO2 concentration inside each growth chamber was controlled by an automated, open-chamber CO2 control system connected to a gas tank containing 100% CO2 ().

The experiment was arranged as factorial design with two CO2 concentrations (ambient CO2 concentration at 400 ± 10 μmol⋅mol-1 and elevated CO2 concentration at 800 ± 10 μmol⋅mol-1) and two temperature treatments [30/25°C (day/night, optimal temperature control) and 45/40°C (day/night, heat stress)]. Each treatment was repeated in four pots of plants (four replicates).

Measurements of Physiological Indexes

Leaf net photosynthetic rate (Pn) was determined by inserting 4–5 individual leaves (second full-expanded from the top) collected from each pot to a 6 cm2 cuvette with a portable infrared gas analyzer (Li-6400, LI-COR, Inc., Lincoln, NB, United States). Leaves were placed in a leaf chamber with a built-in red and blue light source of the Li-6400 with the light intensity of 800 μmol photon⋅m-2⋅s-1.

For leaf chlorophyll content (Chl), 0.2 g of fresh leaves were detached from plants and then immersed in dimethyl sulfoxide (DMSO) in dark for at least 72 h for a complete extraction of total chlorophyll. The absorbance of the Chl extract was measured at wavelengths of 663 and 645 nm, respectively, by a spectrophotometer (Ultrospec 2100 pro, Biochrom Ltd., Cambridge, England) to calculate Chla and Chlb content. Chl was determined as described by . For photochemical efficiency (Fv/Fm), chlorophyll fluorescence (the ratio of variable to maximum fluorescence as Fv/Fm) was measured by a fluorescence induction monitor (Bioscientific Ltd., Herts, United Kingdom) following 30 min dark acclimation through leaf tips.

Metabolites Extraction and Quantification

The extraction procedure was conducted following the method of and . Leaf samples collected at 28 days of treatment were collected and immediately frozen in liquid nitrogen, then stored at -80°C for metabolic profiling analysis. For each sample, frozen dry leaves were ground to a fine powder with liquid nitrogen, and then 25 mg of powder was transferred into a 10 mL microcentrifuge tubes, and extracted in 1.4 mL of 80% (v/v) aqueous methanol at 23°C for 2 h. Ribitol solution of 10 μL (2 mg⋅mL-1 water) as an internal standard was added prior to incubation. Then, extraction was performed in a water bath at 70°C for 15 min. Tubes were centrifuged for 30 min at 9660 gn and the supernatant was decanted into new tubes, 1.4 mL of water and 0.75 mL of chloroform were added. The mixture was vortexed thoroughly and centrifuged for 15 min at 5025 gn and then 1 mL of the polar phase (methanol/water) was pipetted into HPLC vials and dried in a centrifugal concentrator (Centrivap, Labconco Corporation, Kansas City, MO, United States). The dried polar phase was methoximated with 80 μL of 20 mg⋅mL-1 methoxyamine hydrochloride at 30°C for 90 min and then was trimethylsilylated with 80 μL N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) (with 1% TMCS) for 60 min at 70°C.

The Gas Chromatography-Mass Spectrometer (GC-MS) analysis was modified from . The derivatized extracts were analyzed with a GC coupled with a TurboMass-Autosystem XL MS (Perkin Elmer Inc., Waltham, MA, United States). A 1 μL extracts was injected into a DB-5MS capillary column (30 m × 0.25 mm × 0.25 μm, Agilent J & W Scientific, Folsom, CA, United States). The inlet temperature was held at 260°C. After a 6.5 min solvent delay, initial GC oven temperature was maintained at 60°C; 1 min after injection, the GC oven temperature was raised to 280°C at a rate of 5°C⋅in-1, and finally maintained at 280°C for 15 min. The injection temperature was set at 280°C and the ion source temperature was adjusted to 200°C. Helium was used as the carrier gas with a constant flow rate of 1 mL⋅min-1. The measurements were performed through electron impact ionization (70 eV) in the full scan mode (m/z 30–550). The detected metabolites were identified with Turbomass 4.1.1 software (PerkinElmer Inc., Waltham, MA, United States). For GC/MS results, compounds were identified based on retention time (RT) and comparison with reference spectra in mass spectral libraries.

Protein Extraction and Quantification

Leaf samples were collected from each tube at 28 days, immediately frozen in liquid nitrogen, then ground into fine powder and stored at -80°C until analysis. Proteins were extracted using the trichloroacetic acid (TCA)/Acetone method described from . Leaf powder samples (0.5 g) were homogenized on ice in precipitation solution (10% TCA and 0.07% 2-mercaptoethanol in acetone) for 10 min and then incubated at -20°C for 2 h. The protein pellet was collected and washed with cold acetone containing 0.07% 2-mercaptoethanol until the supernatant was colorless. The pellet was then vacuum-dried and suspended in resolubilization solution [8 M urea, 2 M thiourea, 2% CHAPS, 1% dithiothreitol (DTT), and 1% pharmalyte]. The suspension was centrifuged at 21000 g for 20 min and the supernatant was collected for further protein quantification. Protein content was determined using the method of . A 10 μL aliquot of protein extract was mixed with 0.5 mL of a commercial color reagent (Bio-Rad Laboratories, Hercules, CA, United States) by a bovine serum albumin (BSA) standard. The absorbance was measured spectrophotometrically at 595 nm between 5 and 30 min after reaction.

Two-Dimensional PAGE and Protein Analysis

An IPGPhor apparatus (GE Healthcare, Waukesha, WI, United States) was used for the first isoelectric focusing (IEF) described by . The extracts containing 300 μg of sample protein were used for IEF in immobilized pH gradient (IPG) strips (pH 3.0–10.0, linear gradient, 13 cm), formed by rehydrating strips for 12 h at room temperature in 250 μL of rehydration buffer (8 M urea, 2 M thiourea, 2% CHAPS, 1% DTT, 1% v/v IPG buffer, and 0.002% bromophenol blue). Following IEF, the IPG strips were equilibrated for 15 min twice at room temperature in equilibration buffer (50 mM Tris–HCl pH 8.8, 6 M urea, 30% glycerol, 2% SDS, and 1% DTT), then transferred to the same equilibration buffer containing 2.5% iodoacetamide instead of 1% DTT. The second dimension electrophoresis was run on a 12.5% SDS–polyacrylamide gel with a Hoefer SE 600 Ruby electrophoresis apparatus (GE Healthcare, Waukesha, WI, United States). The running conditions were 5 mA per strip for 30 min followed by 20 mA per strip for about 5 h. Gels were stained with Coomassie brilliant blue G-250 and scanned using a Personal Densitometer SI (63-0016-46, GE Healthcare, Waukesha, WI, United States).

Gel images analysis was performed by Progenesis software (Nonlinear Dynamics, Durham, NC, United States). Automatic default spot analysis settings were coupled with manual correction and editing of spot features. The spot volumes were normalized as a percentage of the total volume of all spots on the gel to correct the variability due to staining. Variance analysis of data was used to test the treatment effects on each transgenic line.

Selected protein spots were manually excised from gels and subjected to a trypsin digestion. The peptides were identified by MALDI-TOF-MS as described by . Data were searched against the National Center for Biotechnology Information (NCBI) database. Proteins containing at least two peptides with a confidence interval value >95% were considered to be successfully identified ().

Protein functional classification was performed by Mapman software () in combination with the criteria proposed by . The identified proteins were distributed to different subcellular location by SUBA (). Gene ontology (GO) for biological process, molecular function and cellular component was conducted by the agrigo database1; threshold was -log10 > 4 ().

Statistical Analysis

Data were analyzed using statistics software (SPSS 13.0; SPSS Inc., Chicago, IL, United States). Analysis of variance (ANOVA) was used to determine differences among treatment effects at a given treatment time. The means ± SE were calculated for each parameter. When a particular F-test was significant, means were tested with least significant difference (LSD) at a confidence level of 0.05.

Results

Physiological Effects of Elevated CO2

Under normal temperature, elevated CO2 significantly increased Pn and Chl (Figures 1A,B) while it had no significant effects on Fv/Fm (Figure 1C). Under heat stress, plants grown at elevated CO2 had significantly higher Pn (Figure 1A), Chl (Figure 1B), and Fv/Fm (Figure 1C) than that at ambient CO2 concentration.

FIGURE 1

Effects of Elevated CO2 on Metabolic Profiles

A total of 53 metabolites, including 18 organic acids and phosphoric acid, 12 amino acids, 18 sugars and 4 sugar alcohols, in responsive to elevated CO2 and heat stress were identified and quantified by GC-MS. The name, RT, derivative and mass to charge (m/z) as well as the relative expression of each metabolite was presented in Figure 2 and Table 1.

FIGURE 2

Table 1

CompoundRTDerivativem/zCompoundRTDerivativem/z
Pyruvic acid8.8O-TMSa,MEOX1b174Erythrose30.306O-3TMS,MEOX1205
Lactic acid9.091O-2TMS147Tagatose31.15O-5TMS,MEOX1103
Propenoic acid9.819O-2TMS147Pimelic acid32.978O-3TMS300
Alanine10.24N,O-TMS116Myo-Inositol33.289O-6TMS305
Oxalic acid11.267O-2TMS147Allose33.641O-5TMS,MEOX1319
Valine13.331N,O-TMS144Glucopyranose33.802O-6TMS389
Glycerol14.977O-2TMS205Cellobiose35.218O-8TMS204
Isoleucine15.47O-2TMS117Gulose35.932O-5TMS204
Proline15.547N,O-TMS142Maltose40.354O-8TMS204
Glycine15.785N,N,O-TMS174Galacturonic acid40.611O-5TMS204
Succinic acid16.059O-2TMS147Mannobiose40.972O-8TMS204
Glyceric Acid16.4167O-3TMS147Sucrose42.414O-8TMS361
Serine17.274N,O,O-TMS204Galactinol47.19O-9TMS204
Threonine17.929N,O,O-TMS218Gentiobiose49.095O-8TMS204
Malic acid20.565O-3TMS147Psicose26.779O-5TMS,MEOX1103
5-Oxoproline21.276O-2TMS156Mannopyranose36.968O-4TMS204
Aspartic acid21.328O-3TMS232Fucose37.747O-4TMS204
GABA21.528N,N,O-TMS174Glucuronic acid31.323O-3TMS317
Lysine21.611N,N,O-TMS174Turanose31.491O-8TMS361
Threonic acid22.286O-4TMS292Gluconic acid31.59O-6TMS333
α-Ketoglutaric acid22.666O-2TMS,MEOX1198Palmitic acid32.566O-TMS313
Glutamic acid23.699N,O,O-TMS246Oxaloacetic acid32.701O-3TMS147
Lyxose24.649O-4TMS103Phosphoric acid35.269O-5TMS357
Mucic acid27.687O-6TMS333
Shikimic acid27.912O-4TMS204
Citric acid28.12O-4TMS273
Fructose29.109O-5TMS,MEOX1307
Galactose29.47O-5TMS,MEOX1319
Glucose29.611O-5TMS,MEOX1319
Mannitol30.236O-6TMS319

Metabolites identified by GC-MS in response to different CO2 concentrations and temperatures in leaves of bermudagrass at 28 days of treatments.

RT, retention time; m/z, mass to charge ratio; aTrimethylsilyl derivative(s), bMethoxime derivative(s).

Total content of organic acids, amino acids, sugars, and sugar alcohols were presented in Figure 3. Under normal temperature, no effects of elevated CO2 were detected on total content of organic acids, amino acids, sugars and sugar alcohols compared with ambient CO2. Under heat stress, elevated CO2 resulted in significant increases in the content of organic acids, amino acids, sugars, and sugar alcohols by 52%, 2.79-fold, 29% and 30%, respectively (Figure 3).

FIGURE 3

For organic acids, under heat stress, plants grown at elevated CO2 exhibited significantly lower content of mucic acid, galacturonic acid, lactic acid, but higher content of pyruvic acid, α-ketoglutaric acid, citric acid, glyceric acid, pimelic acid, malic acid and threonic acid compared to those with ambient CO2 (Figures 4AC). Plants exposed to elevated CO2 had significantly lower content of pyruvic acid and α-ketoglutaric acid than those with ambient CO2 under normal temperature (Figure 4B).

FIGURE 4

For amino acids, under heat stress, the content of all amino acids were increased by elevated CO2 except glycine compared with ambient CO2 (Figure 5). The content of alanine, GABA, and serine (Figure 5C) was significantly lower and the content of aspartic acid, isoleucine, lysine and glutamic acid (Figure 5D) was significantly higher in plants exposed to elevated CO2 compared with ambient CO2 treatments under normal temperature.

FIGURE 5

For sugars and sugar alcohols, the content of gentiobiose, gulose, lyxose (Figure 6A), and myo-inositol (Figure 7) was decreased while that of erythrose and glucopyranose (Figure 6B) was increased by elevated CO2 compared to plants grown at ambient CO2 concentration under normal temperature. Under heat stress, 8 out of 18 sugars, and two sugar alcohols (mannitol and galactinol) exhibited increases in the content in plants exposed to elevated CO2 compared with ambient CO2 (Figures 6, 7).

FIGURE 6

FIGURE 7

Out of 53 identified metabolites, 43 were placed into the metabolic pathways associated with GABA shunt, TCA cycle, sugar and amino acid metabolism (Figure 8). These 43 metabolites included 16 organic acids, 12 amino acids, 11 sugars and 4 sugar alcohols. Under heat stressed conditions, elevated CO2 enhanced the accumulation of metabolites associated with GABA shunt, sugar and amino metabolisms.

FIGURE 8

Proteomic Responses to Elevated CO2

A total of 70 and 53 protein spots were differentially expressed in leaves of bermudagrass due to elevated CO2 compared to those at ambient CO2 under normal and high temperature, respectively (Figure 9 and Table 2).

FIGURE 9

Table 2

Spot no.Unipro. IDPro. name [species]pIMWPM
n38A0A0A9D510Uncharacterized protein [Arundo donax]11.47278620592.0182
n196C0PFV4Cytokinin inducible protease1 [Zea mays]6.235711020414
n211A0A0A9R3Q6Uncharacterized protein [Arundo donax]9.639793410972.5892
n238C1K9J1Heat shock protein 90 [Zea mays]5.004387800905
n248Q8W0Q7Methionine synthase protein [Sorghum bicolor]5.93044383788.726
n250Q8W0Q7Methionine synthase protein [Sorghum bicolor]5.93044383788.726
n300X4Z319Heat shock protein 70 [Saccharum hybrid cultivar]5.12775471034.477
n303A4ZYQ0Chloroplast heat shock protein 70 [Pennisetum americanum]5.23328473010.55
n311K4AEH8Glutathione S-transferase [Arabidopsis thaliana]5.505119325757.4872
n324Q7SIC9Transketolase, chloroplastic [Zea mays]5.46634772993.412
n382K3XFX0Phosphoglycerate mutase [Arabidopsis thaliana]5.738609363645.4724
n411A0A096PMM2Chaperonin-60 alpha [Arabidopsis thaliana]5.053413463186.3972
n417C0PHP3Putative TCP-1/cpn60 chaperonin family protein [Zea mays]4.75060344074.174
n447A0A059PYZ3Catalase [Saccharum hybrid cultivar]6.579429656439.8512
n476A0A059Q9W7ATP synthase subunit alpha, chloroplastic [Neyraudia reynaudiana]5.723014855674.8047
n486A0A024GW45ATP synthase subunit alpha, chloroplastic [Lecomtella madagascariensis]5.86518155704.876
n522K3Z2G6ATP synthase subunit alpha [Setaria italica]5.702720655314.3913
n525A0A024GW49ATP synthase subunit beta, chloroplastic [Lecomtella madagascariensis]5.30698453954.826
n530A0A059Q9X1ATP synthase subunit beta, chloroplastic [Neyraudia reynaudiana]5.301643453997.8438
n537A0A024BLC0ATP synthase subunit beta [Pennisetum americanum]5.306983953910.7657
n551A0A0G2UKF5Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit [Orinus thoroldii]6.233680751506.5686
n561A0A059Q9V4Ribulose-1,5-bisphosphate carboxylase/oxygenas [Neyraudia reynaudiana]6.036079452724.9276
n574A0A0U5GUY4Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit [Neostapfiella perrieri]6.339851449871.6557
n598A0A059Q008Elongation factor 1-alpha [Saccharum hybrid cultivar]9.139488249276.9933
n616C0P699Elongation factor Tu [Zea mays]6.195442250776.3463
n619A0A077JG84S-adenosylmethionine synthase [Andropogon virginicus]5.564079343045.7854
n660K4AA01Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.100700446993.463
n666K4AA01Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.100700446993.462
n673K3Z5U9Phosphoglycerate kinase [Setaria italica]6.072074949686.2174
n675K3Z5U9Phosphoglycerate kinase [Setaria italica]6.07207549686.223
n678K3XH82Phosphoglycerate kinase [Setaria italica]8.488258450239.9924
n705B6T2L2Sedoheptulose-1,7-bisphosphatase [Zea mays]6.07453241816.73
n706K3YTN2Glutamine synthetase [Setaria italica]5.512168939158.0482
n736K3XHJ0Aspartate aminotransferase [Setaria italica]8.802711550210.4554
n760K3XHJ0Aspartate aminotransferase [Setaria italica]8.802711550210.4554
n762A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
n769P0C1M0ATP synthase subunit gamma, chloroplastic [Zea mays]8.437202539789.8075
n770A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
n776K3YS38Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]9.386756951746.2142
n781K3ZIS7Fructose-bisphosphate aldolase [Setaria italica]6.265724242104.9795
n784C4JBS8Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]6.45926736494.673
n793C0PD30Fructose-bisphosphate aldolase [Zea mays]6.373924338146.5524
n794K3ZIS7Fructose-bisphosphate aldolase [Setaria italica]6.265724242104.9795
n795A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
n808K3YIG5Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.972602836597.8313
n814A0A140GYJ8Cysteine synthase C1 [Arabidopsis thaliana]7.728721640549.9372
n816K3YIG5Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.972602836597.8313
n827K3Z7Q4Malate dehydrogenase [Setaria italica]8.22945435523.893
n865B6TEW2Ferredoxin–NADP reductase, leaf isozyme [Zea mays]8.37290237506.142
n941B4F9R9Oxygen-evolving enhancer protein 1 [Zea mays]5.53951335079.655
n942B4F9R9Oxygen-evolving enhancer protein 1 [Zea mays]5.53951335079.655
n968B6SQQ0Inorganic pyrophosphatase [Zea mays]5.78838331736.972
n1019J9QDZ6Ascorbate peroxidase [Saccharum hybrid cultivar]5.17603327159.673
n1040B4FT85Isochorismate synthase 1 [Zea mays]7.85582729470.342
n1066B4FNR1Chlorophyll a-b binding protein 2 [Zea mays]5.14025127815.743
n1073B6UG30Triosephosphate isomerase [Zea mays]6.13968732392.872
n1076B4FNR1Chlorophyll a-b binding protein 2 [Zea mays]5.14025127815.742
n1117B6SS26Adenylate kinase [Zea mays]6.79027631139.682
n1118K3YIW7Adenosine monophosphate kinase [Arabidopsis thaliana]7.687599231556.0692
n1127B6SUC4Chlorophyll a-b binding protein 8 [Zea mays]8.94081928984.342
n1173C4J9M72-cys peroxiredoxin BAS1 [Zea mays]5.80782328272.363
n1195B6SSN3Chlorophyll a-b binding protein 6A [Zea mays]6.21477526309.183
n1364B4F9N4Cytochrome b6-f complex iron-sulfur subunit [Zea mays]8.579048224054.5083
n1403A0A0A6Z9F5Photosystem I reaction center subunit II [Saccharum hybrid cultivar]9.91397921844.14
n1410A0A0A6Z9F5Photosystem I reaction center subunit II [Saccharum hybrid cultivar]9.91397921844.15
n1473B6SPC1Photosystem I reaction center subunit IV A [Zea mays]9.78665914893.94
n1507A0A024GWT9ATP synthase epsilon chain, chloroplastic [Lecomtella madagascariensis]5.02799215245.62
n1511B4G259Photosystem II Subunit Q [Arabidopsis thaliana]9.77191923132.724
n1650A0A0A9IAK2Ribulose bisphosphate carboxylase small chain [Arundo donax]6.30663314832.182
n1802O65101Photosystem I reaction center subunit VI, chloroplastic [Zea mays]10.0983414929.32
h405B5AMJ8Alpha-1,4 glucan phosphorylase [Zea mays]6.856071594452.8242
h481Q8W0Q7Methionine synthase protein [Sorghum bicolor]5.930442883788.7256
h634A0A096QX48Succinate dehydrogenase [ubiquinone] flavoprotein subunit [Zea mays]6.042381363934.2212
h712A0A096PMM2CPN60A [Arabidopsis thaliana]5.053413463186.3972
h716A0A096RAX3Malic enzyme [Zea mays]8.000450167809.0362
h724A0A096RAX3Malic enzyme [Zea mays]8.000450167809.0364
h730C0PHP3Putative TCP-1/cpn60 chaperonin family protein [Zea mays]4.750602744074.1734
h736C0PHP3Putative TCP-1/cpn60 chaperonin family protein [Zea mays]5.469337564030.345
h743A0A059PYZ3Catalase [Saccharum hybrid cultivar R570]6.579429656439.8512
h826A0A024GW49ATP synthase subunit beta, chloroplastic [Lecomtella madagascariensis]5.306983953954.826
h839C5Z2J6Catalase [Sorghum bicolor]6.615745556841.3323
h843A0A0G2UKF5Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit [Orinus thoroldii]6.233680751506.5686
h844C5Z2J6Catalase [Sorghum bicolor]6.615745556841.3323
h845A0A059Q0R4NADP-dependent glyceraldehyde-3-phosphate dehydrogenase [Saccharum hybrid cultivar R570]6.800315953254.562
h903A0A059Q008Elongation factor 1-alpha [Saccharum hybrid cultivar R570]9.139488249276.9933
h951A0A077JG84S-adenosylmethionine synthase [Andropogon virginicus]5.564079343045.7854
h964K3ZIK0Rubisco activase [Arabidopsis thaliana]6.198219347531.2554
h1015K4AA01Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.100700446993.463
h1022K4AA01Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.100700446993.462
h1067K3XH82Phosphoglycerate kinase [Setaria italica]8.488258450239.992
h1150A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
h1154B6T2L2Sedoheptulose-1,7-bisphosphatase [Zea mays]6.074531641816.74
h1157B6T2L2Sedoheptulose-1,7-bisphosphatase [Zea mays]6.074531641816.73
h1162A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
h1164K3YS38Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]9.386756951746.2142
h1167K3YIG5Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.972602836597.8313
h1176P0C1M0ATP synthase subunit gamma, chloroplastic [Zea mays]8.437202539789.8075
h1182K3YIG5Glyceraldehyde-3-phosphate dehydrogenase [Setaria italica]6.972602836597.8313
h1192A0A096TAE3Glyceraldehyde-3-phosphate dehydrogenase [Zea mays]7.001228342856.792
h1211K3ZIS7Fructose-bisphosphate aldolase [Setaria italica]6.265724242104.9792
h1230A0A140GYJ8Cysteine synthase C1 [Arabidopsis thaliana]7.728721640549.9372
h1233C0PD30Fructose-bisphosphate aldolase [Zea mays]6.373924338146.5524
h1237K3Z7Q4Malate dehydrogenase [Setaria italica]8.22945435523.893
h1242K3Z7Q4Malate dehydrogenase [Setaria italica]8.22945435523.893
h1258K3XJN7Malate dehydrogenase [Setaria italica]7.671791135479.8543
h1303B6TEW2Ferredoxin–NADP reductase, leaf isozyme [Zea mays]8.372901937506.142
h1384C0PK05Lactoylglutathione lyase [Zea mays]5.825767532344.83
h1389C0PK05Lactoylglutathione lyase [Zea mays]5.825767532344.8043
h1398B4F9R9Oxygen-evolving enhancer protein 1 [Zea mays]5.539512635079.655
h1438K3YUP7Pyrophosphorylase 6 [Arabidopsis thaliana]5.782615731748.9442
h1504B4FT85Isochorismate synthase 1 [Zea mays]7.855827329470.342
h1528K3Y8C6Ascorbate peroxidase 4 [Arabidopsis thaliana]8.165473938125.0193
h1533B6TVL8APx2-Cytosolic Ascorbate Peroxidase [Zea mays]5.282951427201.753
h1584B4FQW0Stem-specific protein TSJT1 [Zea mays]5.23328424666.2512
h1586B6SS26Adenylate kinase [Zea mays]6.790275631139.684
h1655B6SUC4Chlorophyll a-b binding protein 8 [Zea mays]8.940818828984.343
h1854A0A0B4J349Peptidyl-prolyl cis-trans isomerase[Zea mays]9.401176526443.9813
h1900B4F9N4Cytochrome b6-f complex iron-sulfur subunit [Zea mays]8.579048224054.5083
h1905A0A0A6Z9F5Photosystem I reaction center subunit II [Saccharum hybrid cultivar]9.913978621844.14
h1930B4F9N4Cytochrome b6-f complex iron-sulfur subunit [Zea mays]8.579048224054.5083
h2067B6SPC1Photosystem I reaction center subunit IV A [Zea mays]9.786659214893.93
h2087B6ST36Chloroplast oxygen-evolving complex/thylakoid lumenal 25.6kDa protein [Zea mays]9.344459526165.122
h2563B4FAC2Photosystem I reaction center subunit N [Zea mays]9.211692815485.723

Differentially expressed proteins in response to different CO2 concentrations and temperatures by comparison between elevated CO2 and ambient CO2 in leaves of bermudagrass at 28 days of treatments.

n, normal temperature (30°C); h, heat stress (45°C); pI, isoelectric point; MW (kDa), molecular weight; PW, the number of unique peptides matched.

Those 70 proteins up- or down-regulated by elevated CO2 under normal temperature were found in plastid (67.1%), cytosol (18.6%), mitochondrion (5.7%), peroxisome (4.3%), cytosol and plasma membrane (1.4%) and unknown locations (2.9%) (Figure 10). GO category enrichment showed that 70 proteins participated in various biological processes (metabolic process, response to stress, generation of precursor metabolites and energy, photosynthesis, carbohydrate metabolic process, glycolysis, protein folding, and carbohydrate biosynthetic process), molecular functions (catalytic activity and oxidoreductase activity) and cellular components (intracellular, cell, cytoplasm, organelle, plastid, chloroplast, membrane, protein complex, thylakoid, mitochondrion, envelope, plastoglobule, cytosol, stromule, and photosystem) (Figure 11A).

FIGURE 10

FIGURE 11

The majority of these 53 proteins up- or down-regulated by elevated CO2 under heat stress mainly distributed in plastid (62.3%) followed by cytosol (17%) (Figure 10). GO category enrichment indicated that the biological processes regulated by CO2 included cellular metabolic process, responses to stress, response to abiotic stimulus, generation of precursor metabolites and energy, photosynthesis, response to inorganic substance, electron transport chain, carbohydrate metabolism, molecular functions (catalytic activity and oxidoreductase activity) and cellular components (cytoplasm, intracellular, cell, organelle, plastid, chloroplast, membrane, mitochondrion, thylakoid, protein complex, envelope, photosystem, and stromule) (Figure 11B).

Based on the Venn analysis, 18 and 19 differential proteins were up-regulated by elevated CO2 only under either normal temperature or heat stress, respectively (Figure 12A). Elevated CO2 caused 12 proteins to be up-regulated regardless of temperature (Figure 12A). A total of 40 proteins were down-regulated by elevated CO2 under normal and high temperature (Figure 12B). There were 27 proteins down-regulated under normal temperature and 7 proteins down-regulated under heat stress alone due to elevated CO2 treatment. The differentially expressed proteins in responses to elevated CO2 and heat stress were classified into different functional categories (Figure 13).

FIGURE 12

FIGURE 13

Under normal temperature condition, the differential proteins caused by elevated CO2 compared with ambient CO2 were involved in photosynthesis (55.7%), followed by amino acid metabolism (8.6%), glycolysis (7.1%), protein synthesis and degradation (4.3%), stress defense (8.6%), nucleotide metabolism (4.3%), and the remaining (11.4%) including those unknown functions (Figure 13A). For proteins related to photosynthesis, significant increases in the relative fold change were found in ATP synthase subunit (ATPA, n476) by 1.5-fold, rubisco large subunit (RBCL, n551, n561, n574) by 1.1- to 2.6-fold, glyceraldehyde-3-phosphate dehydrogenase (GAPDH, n660, n762, n770, n795) by 1.1- to 1.4-fold, ATP synthase subunit gamma (ATPC, n769) by 1.6-fold, fructose-bisphosphate aldolase (FBA, n781) by 1.4-fold, ferredoxin-NADP reductase (FNR, n865) by 1.3-fold, oxygen-evolving enhancer protein (OEE, n941) by 1.6-fold, chlorophyll a-b binding protein (LHC, n1066, n1076, n1127, n1195) by 1.3- to 2.4-fold, cytochrome b6-f complex iron-sulfur subunit (PGR, n1364) by 2.0-fold, rubisco small chain (RBCS, n1650) by 1.5-fold (Figure 14A). Other 21 proteins involved in photosynthesis [(chaperonin-60 alpha, CPN60A), n411 by 1.14-fold; (cpn60 chaperonin family protein, CPN60B), n417 by 1.52-fold; n486, ATPA by 1.53-fold; n525, n530 and n537, ATPB by 1.17- to 1.23-fold; n666, GAPDH by 1.23-fold; (Phosphoglycerate kinase, PGK), n673, n675 and n678 by 1.17- to 1.21-fold; (Sedoheptulose-1,7-bisphosphatase, SBPase), n705 by 1.18-fold; n793 and n794, FBA by 1.22 – 1.35; n942, OEE by 1.37-fold; n1073, TIM by 1.36-fold; (Photosystem I reaction center subunit, Psa), PsaD, n1403 and n1410 by 1.47- to 2.11-fold; n1473, PsaE by 1.23-fold; n 1507, ATPE by 1.3-fold; n1511, PsbQ by 2.02-fold; n1802, PsaH by 1.89-fold] were significantly down-regulated by elevated CO2 under normal temperature compared with ambient CO2 (Figure 14A).

FIGURE 14

Among proteins associated with the function of protein synthesis and degradation, cytokinin inducible protease (CLPC, n196) had a 1.9-fold up-regulation and the other two [(elongation factor, EF), n598 and n616] with 1.2- to 3.0-fold down-regulation compared with ambient CO2 under normal temperature. Transketolase (n324) involved in oxidative pentose phosphate showed a 1.8-fold increase in response to elevated CO2 (Figure 14A). Proteins involved in amino acid metabolism exhibited increases in methionine synthase protein (MS, n248, n250) by 1.4- to 1.7-fold and cysteine synthase C1 (CSase, n814) by 1.6-fold as well as decreases in aspartate aminotransferase (ASP, n736, n760) by 1.3- to 1.5-fold and S-adenosylmethionine synthase (SAMS, n619) by 1.3-fold in plants grown at elevated CO2 compared with ambient CO2. GAPDH associated with glycolysis were all down-regulated by 1.4- to 1.9-fold (n382, n776, n784) under elevated CO2 except n808, which was increased by 2.4-fold under elevated CO2 concentration. There were six proteins associating with stress defense of which four proteins [Heat shock protein (Hsp) 90, n238; Hsp 70, n300; Catalase, n447; Ascorbate peroxidase, n1019] were up-regulated by 1.2- to 1.4-fold and tow proteins (Hsp 70, n303; 2-cys peroxiredoxin BAS1, n1173) were down-regulated by 1.1- to 1.2-fold under elevated CO2 concentration (Figure 14B).

Under heat stress, elevated CO2-regulated proteins were classified the into following functional categories: photosynthesis (52.8%), TCA cycle (11.3%), stress defense (9.4%), glycolysis (7.5%), amino acid metabolism (5.7%), nucleotide metabolism (3.8%), metal handling (1.9%), major CHO metabolism (1.9%), protein synthesis (1.9%), cell cycle (1.9%) and transport (1.9%) (Figure 13B). Proteins associated with photosynthesis were mainly up-regulated by elevated CO2 under heat stress, including CPN60A (h712) by 1.2-fold, CPN60B (h730, h736) by 1.2- to 1.3-fold, ATP synthase subunit beta (ATPB, n826) by 1.4-fold, RBCL (h843) by 1.5-fold, PGK, (h1067) by 1.3-fold, GAPDH (h1150, h1162, h1162, h1192) by 1.3-fold, SBPase (h1154, h1157) by 1.3-fold, ATPC (h1176) by 1.3-fold, FBA (h1211, h1233) by 1.3- to 1.5-fold, OEE (h1398) by 2.1-fold, LHC (h1655) by 1.6-fold, PGR (h1900, h1930) by 1.2- to 1.7-fold, PsaD (h1905) by 2.5-fold, photosystem I reaction center subunit N (PsaN, h2563) by 1.5-fold (Figure 15A). One protein associated with protein synthesis was upregulated by elevated CO2 under heat stress (Figure 15A). All proteins involved in amino acid metabolism and nucleotide metabolism were up-regulated by elevated CO2 under heat stress, including MS (1.6-fold), SAMS (1.5-fold), CSase (1.4-fold), pyrophosphorylase 6 (PPa6, 1.7-fold), adenylate kinase (ADK, 1.5-fold). Most TCA cycle related proteins [(Malic enzyme, ME), h716 and h724] by 1.4- to 1.5-fold; [(Malate dehydrogenase, MDH), h1237, h1242 and h1258] by 1.2- to 1.5-fold)] were up-regulated by elevated CO2 compared with ambient CO2 during heat stress (Figure 15B).

FIGURE 15

Discussion

Previous studies have shown positive effects of elevated CO2 on plant growth of C4 species under optimal temperature conditions (). In this study, elevated CO2 significantly improved physiological activities of C4 bermudagrass under heat stress or mitigated heat stress damages, as manifested by physiological indexes, including higher leaf Pn, Fv/Fm and Chl. The positive physiological effects were associated with changes in various metabolic pathways regulated by elevated CO2. Metabolic and proteomic analysis in this study indicated that the underlying mechanisms of elevated CO2-mitigation of heat stress were mainly related to photosynthesis, respiration (glycolysis and TCA cycle), amino acid metabolism, and GABA shunt, and some of the metabolic factors regulated by elevated CO2 in the C4 grass species, bermudagrass, in this study are in common and some are different from those previously found in C3 grass species (; ). Due to the large number of metabolites and the complexity of metabolic pathways involved in CO2 effects, the following sections focused on the discussion of unique or different metabolic pathways found in bermudagrass in our study from those findings previously reported in other C3 plant species.

Proteins and Metabolites in Photosynthesis Regulated by Elevated CO2 under Heat Stress

In our present study, 67 out of 123 proteins (54.5%) associated with photosynthetic pathways, including proteins involved in electron transport chain and Calvin cycle, were responsive to elevated CO2 concentration under normal and high temperatures, as shown by the decrease or increase in their abundance (Figure 16). Under heat stress, the majority of proteins involved in photosynthesis exhibited accumulation in response to elevated CO2 concentration, such as ATP synthase subunit (h826, h1176), photosystem I reaction center subunit [(PsaD, h1905) and (PsaN, h2563)] in light reactions of photosynthesis, and fructose-bisphosphate aldolase (FBA, h781, h1211, h1233), phosphoglycerate kinase (PGK, h1067) and sedoheptulose-1,7-bisphosphatase (SBPase, h1154, h1157) in Calvin cycle (Figure 15). FBA is a primary enzyme involved in the sixth reaction of Calvin cycle to convert fructose 1,6-bisphosphate into glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate as well as ATP (). FBA content at the level of protein significantly increased under elevated CO2 in C3 tall fescue under heat stress () and creeping bentgrass (Agrostis stolonifera) under both well-water and drought stress (). SBPase functions as a bisphosphatase enzyme catalyzing sedoheptulose 1,7-bisphosphate dephosphorylation to sedoheptulose-7-phosphate during the regeneration phase of Calvin cycle (). Overexpression of SBPase in C3 tobacco (Nicotiana tabacum) had higher photosynthesis at elevated CO2 compared with that at ambient CO2 under field conditions (). The benefits of SBPase on the stimulation of photosynthesis depended on light intensity (; ). Therefore, our study and another case in creeping bentgrass, conducted in light saturated growth chambers, the abundance of SBPase was enhanced by elevated CO2 under abiotic stresses (). FBA is a primary enzyme involved in the sixth reaction of Calvin cycle to convert fructose 1,6-bisphosphate into glyceraldehyde-3-phosphate (G3P), dihydroxyacetone phosphate and ATP (). In addition, FBA could directly affect ribulose-1,5-bisphosphate (RuBP) regeneration which actions as substrate of carbon fixation (). FBA content at the level of protein showed the grater accumulation in elevated CO2 in C3 tall fescue under heat stress () and creeping bentgrass (Agrostis stolonifera) under both well-water and drought stress (). PGK is a major enzyme catalyzing the phosphorylation of 3-phosphoglycerate to produce 1, 3-bisphosphoglycerate and ADP which is one of vital steps regenerating RuBP during Calvin cycle (). The regulation of FBA and PGK induced by elevated CO2 indicated that elevated CO2 availability in atmosphere could be helpful for sustaining ATP supply and RuBP regeneration for plant growth under heat stress. ATP synthase is a critical enzyme for creating energy storage molecule ATP. Under high CO2 availability, ATP synthase was found to decline in wheat grain (). To our knowledge, our case is the first report on the abundance of ATP synthase and PGK in response to elevated CO2 were found in C4 plant species grown under heat stress. Our previous study in C3 plant species found differential responses of photosynthesis-related proteins to elevated CO2 different from found in bermudagrass in our study. In tall fescue, the abundance of ATP synthase subunit and PGK did not change in response to elevated CO2 under heat stress (). The increase in abundance and activity of a single or some enzyme(s) during photosynthesis could enhance carbon assimilation (). Taken together, the enhanced accumulation of proteins involved in photosynthesis by elevated CO2 under heat stress in bermudagrass suggested that elevated CO2 could help to maintain photosynthesis to withstand the adverse environments as various proteins are involved in the light harvesting, electron transport, and carbohydrate assimilation processes of photosynthesis.

FIGURE 16

Other proteins related to photosynthesis such as GAPDH, OEE, PGR exhibited the enhanced expression in plants grown at elevated CO2 concentration under both temperatures in our study. GAPDH could convert G3P to D-glycerate 1,3-bisphosphate as well as mediating the formation of NADH and ATP (). It has multiple functions, such as two chloroplastic forms playing photosynthetic function locating in chloroplast and one cytosolic form participating in glycolysis in higher plants (; ). In chloroplasts, GAPDH catalyzes a reaction of NADPH-consuming which is regulated by light utilizing thioredoxins and metabolites during Calvin cycle (). Various stresses caused the decline in chloroplastic GAPDH whereas stress-tolerant species exhibited higher GAPDH abundance than stress-sensitive plants, such as creeping bentgrass under heat stress (; ), salinity stress () and drought stress (). Plants with lower GAPDH abundance were generally associated with decreased photosynthetic capacity resulted from reduced RuBP regeneration rate, followed with the decline in accumulation of photosynthetic products (; ). However, elevated CO2 had no effects on chloroplastic GAPDH abundance under heat stressed condition but caused significant decrease under non-stressed control plants (). Overall, our study suggested that enhanced abundance of photosynthesis-related proteins could contribute to the improved photosynthetic activities by elevated CO2, particularly under heat stress, which could be reflected with improved Pn and increased content of sugars, such as fructose, glucose, sucrose, erythrose, and glucopyranose.

Proteins and Metabolites in Respiration Regulated by Elevated CO2 under Heat Stress

It has been widely known that glycolysis and TCA cycle are vital pathways for energy supply, amino acid synthesis and various other biological processes in plants (; ). As substrate of photosynthesis for carboxylation, plants grown at elevated CO2 tended to accumulate the larger amount of non-structural carbohydrates (; ). Most monosaccharides (glucose, fructose, galactose, etc.) as substrate or intermediates play vital roles during glycolysis. Glycolysis pathway could convert glucose into pyruvate via a series of intermediate metabolites and cytosolic GAPDH is one of essential enzymes catalyzing the sixth step of respiratory glycolysis to convert G3P to 1, 3-bisphosphateglycerate (1, 3-BPG) which is one of the most important reactions during the glycolytic pathway (). The increase of pyruvic acid (pyruvate) as the product of glycolysis, followed by the enhanced content of valine, isoleucine and alanine, was partly due to elevated CO2-caused accumulation of glucose under heat stress in our study, since those metabolites are all derived from glucose. In C3 tall fescue, we also observed the significant increases in valine and alanine but not for isoleucine resulted from elevated CO2 under heat stress (). During the pathway of glycolysis, the abundance of GAPDH in cytosol (n776, n784, n816 except n808) and phosphoglycerate mutase (PGAM, n382) exhibited the down-regulation in response to elevated CO2 rather than ambient CO2 under normal temperature, while under heat stressed conditions elevated CO2 caused up-regulation in GAPDH (h845, h1164, h1167 except h1182) in bermudagrass (Figure 16). GAPDH might serve as a provider of additional energy for plant growth and development under stressed conditions and stress tolerance could be enhanced by improved abundance of GAPDH to cope with environmental stresses (; ). In C3 plants, no consistent changes were found due to variations in plant species. For example, in tall fescue and creeping bentgrass, the abundance of cytosolic GAPDH exhibited either no changes or decrease under elevated CO2 and heat stressed condition (; ). found that GAPDH gene-silenced lines showed more sensitive traits to drought stress than non-silenced lines in potato (Solanum tuberosum). By contrast, the higher shoot length and weight were detected in GAPDH overexpression transgenic plants compared with wild-type plants (). In the level of transcription in potato, cytosolic GAPDH RNA accumulation was also increased under biological stress (). Therefore, in our study, the higher abundance of GAPDH caused by elevated CO2 was beneficial for energy supply to support plant growth under heat stress.

Malate dehydrogenase (MDH) acts as an enzyme to catalyze the oxidation of malate to oxaloacetate via the reduction of NAD+ to NADH in mitochondrial matrix during TCA cycle (). Environmental stresses including drought (), heat (), salinity () and Al-stress () have been shown to decrease the level of MDH in various plant species. However, limited studies about MDH were found in plants grown at elevated CO2 concentrations, especially under stressed conditions (). In this study, elevated CO2-responsive MDH (n827, h1237, h1242, h1258) involved in TCA cycle exhibited up-regulated expression regardless of temperatures, suggesting that CO2 inhibited the heat-induced reduction in MDH to catalyze the enhanced malate (malic acid) to oxaloacetate (oxaloacetic acid) during malate metabolism.

Amino Acid Metabolism and GABA Shunt Regulated by Elevated CO2 under Heat Stress

In addition to function in TCA cycle, MDH also participates in the process of amino acid synthesis due to the relations among malate, oxaloacetate and aspartate (; ). Several amino acids including aspartate (aspartic acid), methionine, threonine, isoleucine, lysine derived from oxaloacetate and aspartate is the precursor of methionine, threonine, isoleucine and lysine (). Along with the significant increase in malic acid and aspartic acid, the content of threonine, isoleucine and lysine were stimulated by elevated CO2 during heat stress. Furthermore, the content of alanine, valine and serine were also enhanced by elevated CO2 compared with ambient CO2 under heat stress. Alanine, valine and serine are used for synthesis of several proteins and associated with many metabolic processes (). The stimulation of elevated CO2 concentration on the content of alanine, valine and serine was found in other species under abiotic stresses, as previously reported in C3 grass species under heat stress () and tree seedlings under drought stress (). The increase in synthesis of both alanine and valine in present study is directly associated with the higher content of pyruvate (pyruvate acid) which is the final product of glycolysis (). Superior stress tolerance has been reported with the higher content of alanine, valine and serine as well as other amino acids such as GABA, glutamic acid, proline and 5-oxoproline involved in the GABA shunt pathway in plant species, including perennial grasses (; ; ; ,). The GABA shunt was considered to be a part of the TCA cycle during respiration besides its central role in primary carbon and nitrogen metabolism (). In tall fescue, the content of GABA was significantly decreased by elevated CO2 under high temperature (). While, in bermudagrass of this case, GABA and glutamic acid exhibited the opposite response to elevated CO2 under heat stress. Increased GABA caused the enhanced content of alanine and pyruvate which was turned into TCA cycle and proline metabolism. The content of all amino acids except arginine during GABA shunt was increased by elevated CO2 under heat stress suggesting a predominant role of elevated CO2 in carbon and nitrogen metabolism in C4 bermudagrass.

Proteins including methionine synthase (MS), cysteine synthase (CSase) and S-adenosylmethionine synthase (SAMS) associated with amino acid metabolism were up-regulated by 1.4- to 1.6-fold by elevated CO2 under heat stress. MS and SAMS serve as regulators in the synthesis and degradative pathways of various amino acids (). It was detected by the same proteomic analysis that many proteins involved in amino acid metabolism accumulated more or degraded less in stress-tolerant plants, such as MS and SAMS (). CSase functions in the final strep in cysteine synthesis in plants. Plants with overexpressing CSase gene displayed high tolerance to toxic environmental pollutants, such as sulfur dioxide and sulfite (), cadmium toxicity () in tobacco and aluminum toxicity in rice (). The accumulation of many amino acids as well as proteins involved in amino acid metabolism in this study could contribute to elevated CO2-improved heat tolerance.

In summary, elevated CO2 concentration suppressed heat-induced damages in bermudagrass, as shown by the increased Pn, Chl and Fv/Fm. The improvement of heat tolerance under elevated CO2 could be associated with some important metabolic pathways during which proteins and metabolites were up-regulated, including proteins, sugars and/or amino acids involved in light reaction (ATP synthase subunit and photosystem I reaction center subunit) and carbon fixation of photosynthesis (GAPDH, FBA, PGK, SBPase and sugars), glycolysis (GAPDH, glucose, fructose and galactose) and TCA cycle (pyruvic acid, malic acid and MDH) of respiration, amino acid metabolism (aspartic acid, methionine, threonine, isoleucine, lysine, valine, alanine and isoleucine) as well as the GABA shunt (GABA, glutamic acid, alanine, proline and 5-oxoproline). The molecular factors and mechanisms underlying the metabolic changes caused by elevated CO2 during plant responses to heat stress require further investigation.

Statements

Author contributions

JY and BH designed the experiments and wrote the manuscript. JY and RL conducted the experiments. NF helped with the sample analysis. ZY arranged the experiments and did the data analysis.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (31301799) and the Fundamental Research Funds for the Central Universities (KYZ201673).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AbbasiF. M.KomatsuS. (2004). A proteomic approach to analyze salt-responsive proteins in rice leaf sheath.Proteomics420722081. 10.1002/pmic.200300741

  • 2

    AbebeA.PathakH.SinghS. D.BhatiaA.HaritR. C.KumarV. (2016). Growth, yield and quality of maize with elevated atmospheric carbon dioxide and temperature in north-west India.Agric. Ecosyst. Environ.2186672. 10.1016/j.agee.2015.11.014

  • 3

    AlonsoA.PérezP. A.Martinez-CarrascoR. (2009). Growth in elevated CO2 enhances temperature response of photosynthesis in wheat.Physiol. Plant.135109120. 10.1111/j.1399-3054.2008.01177.x

  • 4

    ArnonD. I. (1949). Copper enzymes in isolated chloroplasts. Polyphenoloxidase in Beta vulgaris.Plant Physiol.24115. 10.1104/pp.24.1.1

  • 5

    BenczeS.VeiszO.BedoZ. (2005). Effect of elevated CO2 and high temperature on the photosynthesis and yield of wheat.Cereal Res. Commun.33385388. 10.1556/CRC.33.2005.1.95

  • 6

    BernsteinB. E.MichelsP. A.HolW. G. (1997). Synergistic effects of substrate-induced conformational changes in phosphoglycerate kinase activation.Nature385275278. 10.1038/385275a0

  • 7

    BertrandA.PrévostD.BigrasF. J.CastonguayY. (2007). Elevated atmospheric CO2 and strain of rhizobium alter freezing tolerance and cold-induced molecular changes in alfalfa (Medicago sativa).Ann. Bot.99275284. 10.1093/aob/mcl254

  • 8

    BevanM.BancroftI.BentE.LoveK.GoodmanH.DeanC.et al (1998). Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana.Nature391485488. 10.1038/35140

  • 9

    BohnertH. J.JensenR. G. (1996). Strategies for engineering water-stress tolerance in plants.Trends Biotechnol.148997. 10.1016/0167-7799(96)80929-2

  • 10

    BourguignonJ.RebeilleF.DouceR. (1999). “Serine and glycine metabolism in higher plants,” in Plant Amino Acids, ed.SinghB. K. (New York, NY: Marcel Dekker), 111146.

  • 11

    BradfordM. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding.Anal. Biochem.72248254. 10.1016/0003-2697(76)90527-3

  • 12

    BurgessP.HuangB. (2014). Root protein metabolism in association with improved root growth and drought tolerance by elevated carbon dioxide in creeping bentgrass.Field Crops Res.1658091. 10.1016/j.fcr.2014.05.003

  • 13

    BurgessP.HuangB. (2016). Leaf protein abundance associated with improved drought tolerance by elevated carbon dioxide in creeping bentgrass.J. Am. Soc. Hortic. Sci.1418596.

  • 14

    FaitA.FrommH.WalterD.GaliliG.FernieA. R. (2008). Highway or byway: the metabolic role of the GABA shunt in plants.Trends Plant Sci.131419. 10.1016/j.tplants.2007.10.005

  • 15

    FernieA. R.CarrariF.SweetloveL. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport.Curr. Opin. Plant Biol.7254261. 10.1016/j.pbi.2004.03.007

  • 16

    FigueiredoN.CarrancaC.TrindadeH.PereiraJ.GoufoP.CoutinhoJ.et al (2015). Elevated carbon dioxide and temperature effects on rice yield, leaf greenness, and phenological stages duration.Paddy Water Environ.13313324. 10.1007/s10333-014-0447-x

  • 17

    FukayamaH.FukudaT.MasumotoC.TaniguchiY.SakaiH.ChengW.et al (2009). Rice plant response to long term CO2 enrichment: gene expression profiling.Plant Sci.177203210. 10.1016/j.plantsci.2009.05.014

  • 18

    HamerlynckE. P.HuxmanT. E.LoikM. E.SmithS. D. (2000). Effects of extreme high temperature, drought and elevated CO2 on photosynthesis of the Mojave Desert evergreen shrub, Larrea tridentata.Plant Ecol.148183193. 10.1023/A:1009896111405

  • 19

    HaradaE.ChoiY. E.TsuchisakaA.ObataH.SanoH. (2001). Transgenic tobacco plants expressing a rice cysteine synthase gene are tolerant to toxic levels of cadmium.J. Plant Physiol.158655661. 10.1078/0176-1617-00314

  • 20

    HoaglandD. R.ArnonD. I. (1950). The water-culture method for growing plans without soil.Calif. Agric. Exp. Stn. Circ.347132.

  • 21

    HögyP.ZörbC.LangenkämperG.BetscheT.FangmeierA. (2009). Atmospheric CO2 enrichment changes the wheat grain proteome.J. Cereal Sci.50248254. 10.1016/j.jcs.2009.06.002

  • 22

    HuangB.XuY. (2015). Cellular and molecular mechanisms for elevated CO2-regulation of plant growth and stress adaptation.Crop Sci.551405. 10.2135/cropsci2014.07.0508

  • 23

    Intergovernmental Panel on Climate Change [IPCC] (2007). Climate Change: Fourth Assessment Report.London: Cambridge University Press.

  • 24

    KappacheryS.Baniekal-HiremathG.YuJ. W.ParkS. W. (2015). Effect of over-and under-expression of glyceraldehyde 3-phosphate dehydrogenase on tolerance of plants to water-deficit stress.Plant Cell Tissue Organ Cult.12197107. 10.1007/s11240-014-0684-0

  • 25

    KirkhamM. B. (2011). Elevated Carbon Dioxide: Impacts on Soil and Plant Water Relations.Boca Raton, FL: CRC Press. 10.1201/b10812

  • 26

    LaiS. K.ZhuangS. T.WuY. Z.WangY. X.ZhuJ. G.YangL. X.et al (2015). Impact of elevated atmospheric CO2 concentration and temperature on growth and development of super rice.Am. J. Roentgenol.3412531262.

  • 27

    LaxaltA. M.CassiaR. O.SanllorentiP. M.MadridE. A.AndreuA. B.DaleoG. R.et al (1996). Accumulation of cytosolic glyceraldehyde-3-phosphate dehydrogenase RNA under biological stress conditions and elicitor treatments in potato.Plant Mol. Biol.30961972. 10.1007/BF00020807

  • 28

    LefebvreS.LawsonT.ZakhleniukO. V.LloydJ. C.RainesC. A.FryerM. (2005). Increased sedoheptulose-1,7-bisphosphatase activity in transgenic tobacco plants stimulates photosynthesis and growth from an early stage in development.Plant Physiol.138451460. 10.1104/pp.104.055046

  • 29

    LiZ.YuJ.PengY.HuangB. (2016a). Metabolic pathways regulated by γ-aminobutyric acid (GABA) contributing to heat tolerance in creeping bentgrass (Agrostis stolonifera).Sci. Rep.6:30338. 10.1038/srep30338

  • 30

    LiZ.YuJ.PengY.HuangB. (2016b). Metabolic pathways regulated by abscisic acid, salicylic acid, and γ-aminobutyric acid in association with improved drought tolerance in creeping bentgrass (Agrostis stolonifera).Physiol. Plant.1594258. 10.1111/ppl.12483

  • 31

    MaX.XuQ.MeyerW. A.HuangB. (2016). Hormone regulation of rhizome development in tall fescue (Festuca arundinacea) associated with proteomic changes controlling respiratory and amino acid metabolism.Ann. Bot.118481494. 10.1093/aob/mcw120

  • 32

    MerewitzE. B.DuH.YuW.LiuY.GianfagnaT.HuangB. (2012). Elevated cytokinin content in ipt transgenic creeping bentgrass promotes drought tolerance through regulating metabolite accumulation.J. Exp. Bot.6313151328. 10.1093/jxb/err372

  • 33

    MerewitzE. B.GianfagnaT.HuangB. (2011). Protein accumulation in leaves and roots associated with improved drought tolerance in creeping bentgrass expressing an ipt gene for cytokinin synthesis.J. Exp. Bot.6253115333. 10.1093/jxb/err166

  • 34

    MijeongJ.SoochulP.HawkbinK.MyungokB. (2000). Isolation and characterization of the gene encoding glyceraldehyde-3-phosphate dehydrogenase.Biochem. Biophys. Res. Commun.278192196. 10.1006/bbrc.2000.3732

  • 35

    MorganJ. A.LecainD. R.PendallE.BlumenthalD. M.KimballB. A.CarrilloY.et al (2011). C4 grasses prosper as carbon dioxide eliminates desiccation in warmed semi-arid grassland.Nature476202205. 10.1038/nature10274

  • 36

    MuehlbauerG. J.GengenbachB. G.SomersD. A. (1994). Genetic and amino-acid analysis of two maize threonine-overproducing, lysine-insensitive aspartate kinase mutants.Theor. Appl. Gen.89767774. 10.1007/BF00223717

  • 37

    MusratiR. A.KollárováM.MernikN.MikulásováD. (1998). Malate dehydrogenase: distribution, function and properties.Gen. Physiol. Biophys.17193210.

  • 38

    NojiM.AonoM.SaitoK. (2001). Cysteine synthase overexpression in tobacco confers tolerance to sulfur-containing environmental pollutants.Plant Physiol.126973. 10.1104/pp.126.3.973

  • 39

    PrasadP. V. V.BooteK. J.AllenL. H.Jr. (2006). Adverse high temperature effects on pollen viability, seed-set, seed yield and harvest index of grain-sorghum [Sorghum bicolor (L.) Moench] are more severe at elevated carbon dioxide due to higher tissue temperatures.Agric. For. Meteorol.139237251. 10.1016/j.agrformet.2006.07.003

  • 40

    PrasadP. V. V.BooteK. J.AllenL. H.Jr.ThomasJ. M. G. (2002). Effects of elevated temperature and carbon dioxide on seed-set and yield of kidney bean (Phaseolus vulgaris L.).Glob. Change Biol.8710721. 10.1046/j.1365-2486.2002.00508.x

  • 41

    PrasadP. V. V.BooteK. J.AllenL. H.ThomasJ. M. G. (2010). Super-optimal temperatures are detrimental to peanut (Arachis hypogaea L.) reproductive processes and yield at both ambient and elevated carbon dioxide.Glob. Change Biol.917751787. 10.1046/j.1365-2486.2003.00708.x

  • 42

    PrasadP. V. V.BooteK. J.JcvV.LhjrA. (2004). The carbohydrate metabolism enzymes sucrose-P synthase and ADG-pyrophosphorylase in phaseolus bean leaves are up-regulated at elevated growth carbon dioxide and temperature.Plant Sci.16615651573. 10.1016/j.plantsci.2004.02.009

  • 43

    PriceG. D.EvansJ. R.VonC. S.YuJ. W.BadgerM. R. (1995). Specific reduction of chloroplast glyceraldehyde-3-phosphate dehydrogenase activity by antisense RNA reduces CO2 assimilation via a reduction in ribulose bisphosphate regeneration in transgenic tobacco plants.Planta195369378. 10.1007/BF00202594

  • 44

    PritchardS. H. G.RogersH. O. H.PriorS. A.PetersonC. T. M. (1999). Elevated CO2 and plant structure: a review.Glob. Change Biol.5807837. 10.1046/j.1365-2486.1999.00268.x

  • 45

    QaderiM. M.KurepinL. V.ReidD. M. (2006). Growth and physiological responses of canola (Brassica napus) to three components of global climate change: temperature, carbon dioxide and drought.Physiol. Plant.128710721. 10.1111/j.1399-3054.2006.00804.x

  • 46

    QiuY.SuM.LiuY.ChenM.GuJ.ZhangJ.et al (2007). Application of ethyl chloroformate derivatization for gas chromatography-mass spectrometry based metabonomic profiling.Anal. Chim. Acta583277283. 10.1016/j.aca.2006.10.025

  • 47

    RainesC. A.LloydJ. C.DyerT. A. (1999). New insights into the structure and function of sedoheptulose-1,7-bisphosphatase; an important but neglected Calvin cycle enzyme.J. Exp. Bot.5018.

  • 48

    Ramírez-BenítezJ. E.Chee-GonzálezL.Hernandez-SotomayorS. M. T. (2008). Aluminium induces changes in organic acids metabolism in Coffea arabica suspension cells with differential Al-tolerance.J. Inorg. Biochem.10216311637. 10.1016/j.jinorgbio.2008.03.002

  • 49

    ReadJ. J.MorganJ. A. (1996). Growth and partitioning in Pascopyrum smithii (C3) and Bouteloua gracilis (C4) as influenced by carbon dioxide and temperature.Ann. Bot.77487496. 10.1006/anbo.1996.0059

  • 50

    RizhskyL.LiangH.ShumanJ.ShulaevV.DavletovaS.MittlerR. (2004). When defense pathways collide. The response of Arabidopsis to a combination of drought and heat stress.Plant Physiol.13416831696. 10.1104/pp.103.033431

  • 51

    RoessnerU.WagnerC.KopkaJ.TretheweyR. N.WillmitzerL. (2000). Simultaneous analysis of metabolites in potato tuber by gas chromatography-mass spectrometry.Plant J.23131142. 10.1046/j.1365-313x.2000.00774.x

  • 52

    RosenthalD. M.LockeA. M.KhozaeiM.RainesC. A.LongS. P.OrtD. R. (2011). Over-expressing the C3 photosynthesis cycle enzyme Sedoheptulose-1-7 Bisphosphatase improves photosynthetic carbon gain and yield under fully open air CO2 fumigation (FACE).BMC Plant Biol.11:123. 10.1186/1471-2229-11-123

  • 53

    SchulzesiebertD.HeinekeD.ScharfH.SchultzG. (1984). Pyruvate-derived amino acids in spinach chloroplasts 1.Plant Physiol.76465471. 10.1104/pp.76.2.465

  • 54

    ShiH.JiangC.YeT.TanD. X.ReiterR. J.ZhangH.et al (2014). Comparative physiological, metabolomic, and transcriptomic analyses reveal mechanisms of improved abiotic stress resistance in bermudagrass [Cynodon dactylon (L). Pers.] by exogenous melatonin.J. Integr. Plant Biol.66681694. 10.1093/jxb/eru373

  • 55

    SongY. L.YuJ. J.HuangB. (2014). Elevated CO2-mitigation of high temperature stress associated with maintenance of positive carbon balance and carbohydrate accumulation in Kentucky bluegrass.PLOS ONE9:e89725. 10.1371/journal.pone.0089725

  • 56

    SparlaF.ZaffagniniM.WedelN.ScheibeR.PupilloP.TrostP. (2005). Regulation of photosynthetic GAPDH dissected by mutants.Plant Physiol.13822102219. 10.1104/pp.105.062117

  • 57

    SujathaK. B.UpretyD. C.RaoD. N.RaoP. R.DwivediN. (2008). Up-regulation of photosynthesis and sucrose-P synthase in rice under elevated carbon dioxide and temperature conditions.Plant Soil Environ.54155162.

  • 58

    TaizL.ZeigerE. (2010). Plant Physiology, 5th Edn. Sunderland, MA: Sinauer Associates.

  • 59

    TanzS. K.CastledenI.HooperC. M.VacherM.SmallI.MillarH. A. (2013). SUBA3: a database for integrating experimentation and prediction to define the SUBcellular location of proteins in Arabidopsis.Nucleic Acids Res.4111851191. 10.1093/nar/gks1151

  • 60

    TarzeA.DeniaudA.BrasM. L.MaillierE.MolleD.LarochetteN.et al (2007). GAPDH, a novel regulator of the pro-apoptotic mitochondrial membrane permeabilization.Oncogene2626062620. 10.1038/sj.onc.1210074

  • 61

    ThimmO.BläsingO.GibonY.NagelA.MeyerS.KrügerP.et al (2004). Mapman: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes.Plant J.37914939. 10.1111/j.1365-313X.2004.02016.x

  • 62

    TristanC.ShahaniN.SedlakT. W.SawaA. (2011). The diverse functions of GAPDH: views from different subcellular compartments.Cell. Signal.23317323. 10.1016/j.cellsig.2010.08.003

  • 63

    TschaplinskiT. J.StewartD. B.NorbyR. J. (1995). Interactions between drought and elevated CO2 on osmotic adjustment and solute concentrations of tree seedlings.New Phytol.131169177. 10.1111/j.1469-8137.1995.tb05718.x

  • 64

    WenW.LiK.AlseekhS.OmranianN.ZhaoL.ZhouY.et al (2015). Genetic determinants of the network of primary metabolism and their relationships to plant performance in a maize recombinant inbred line population.Plant Cell2718391856. 10.1105/tpc.15.00208

  • 65

    XuC.HuangB. (2008). Root proteomic responses to heat stress in two Agrostis grass species contrasting in heat tolerance.J. Exp. Bot.5941834194. 10.1093/jxb/ern258

  • 66

    XuC.HuangB. (2010a). Differential proteomic response to heat stress in thermal Agrostis scabra and heat-sensitive Agrostis stolonifera.Physiol. Plant.139192204. 10.1111/j.1399-3054.2010.01357.x

  • 67

    XuC.HuangB. (2010b). Comparative analysis of drought responsive proteins in Kentucky bluegrass cultivars contrasting in drought tolerance.Crop Sci.5025432552. 10.2135/cropsci2010.03.0152

  • 68

    XuC.SibickyT.HuangB. (2010). Protein profile analysis of salt-responsive proteins in leaves and roots in two cultivars of creeping bentgrass differing in salinity tolerance.Plant Cell Rep.29595615. 10.1007/s00299-010-0847-3

  • 69

    XuY.DuH.HuangB. (2013). Identification of metabolites associated with superior heat tolerance in thermal bentgrass through metabolic profiling.Crop Sci.5316261635. 10.2135/cropsci2013.01.0045

  • 70

    YangQ.WangY.ZhangJ.ShiW.QianC.PengX. (2007). Identification of aluminum-responsive proteins in rice roots by a proteomic approach: cysteine synthase as a key player in Al response.Proteomics7737749. 10.1002/pmic.200600703

  • 71

    YuJ.SunL.FanN.YangZ.HuangB. (2015). Physiological factors involved in positive effects of elevated carbon dioxide concentration on bermudagrass tolerance to salinity stress.Environ. Exp. Bot.1152027. 10.1016/j.envexpbot.2015.02.003

  • 72

    YuJ. J.DuH. M.XuM.HuangB. R. (2012a). Metabolic responses to heat stress under elevated atmospheric CO2 concentration in a cool-season grass species.J. Am. Soc. Hortic. Sci.137221228.

  • 73

    YuJ. J.ChenL. H.XuM.HuangB. R. (2012b). Effects of elevated CO2 on physiological responses of tall fescue to elevated temperature, drought stress, and the combined stresses.Crop Sci.5218481858. 10.2135/cropsci2012.01.0030

  • 74

    YuJ. J.YangZ. M.JespersenD.HuangB. R. (2014). Photosynthesis and protein metabolism associated with elevated CO2-mitigation of heat stress damages in tall fescue.Environ. Exp. Bot.997585. 10.1016/j.envexpbot.2013.09.007

Summary

Keywords

bermudagrass, elevated CO2, heat stress, metabolites, protein

Citation

Yu J, Li R, Fan N, Yang Z and Huang B (2017) Metabolic Pathways Involved in Carbon Dioxide Enhanced Heat Tolerance in Bermudagrass. Front. Plant Sci. 8:1506. doi: 10.3389/fpls.2017.01506

Received

31 March 2017

Accepted

15 August 2017

Published

19 September 2017

Volume

8 - 2017

Edited by

Luis A. J. Mur, Aberystwyth University, United Kingdom

Reviewed by

Lucia Guidi, University of Pisa, Italy; M. B. Kirkham, Kansas State University, United States

Updates

Copyright

*Correspondence: Zhimin Yang, Bingru Huang,

These authors have contributed equally to this work.

This article was submitted to Crop Science and Horticulture, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics