ORIGINAL RESEARCH article
Sec. Plant Physiology
Seasonal Variation of Carbon Metabolism in the Cambial Zone of Eucalyptus grandis
- 1Laboratório Max Feffer de Genética de Plantas, Departamento de Genética, Escola Superior de Agricultura Luiz de Queiroz, Universidade de São Paulo, Piracicaba, Brazil
- 2Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden
Eucalyptus species are the most widely hardwood planted in the world. It is one of the successful examples of commercial forestry plantation in Brazil and other tropical and subtropical countries. The tree is valued for its rapid growth, adaptability and wood quality. Wood formation is the result of cumulative annual activity of the vascular cambium. This cambial activity is generally related to the alternation of cold and warm, and/or dry and rainy seasons. Efforts have focused on analysis of cambial zone in response to seasonal variations in trees from temperate zones. However, little is known about the molecular changes triggered by seasonal variations in trees from tropical countries. In this work we attempted to establish a global view of seasonal alterations in the cambial zone of Eucalyptus grandis Hill ex Maiden, emphasizing changes occurring in the carbon metabolism. Using transcripts, proteomics and metabolomics we analyzed the tissues harvested in summer-wet and winter-dry seasons. Based on proteomics analysis, 70 proteins that changed in abundance were successfully identified. Transcripts for some of these proteins were analyzed and similar expression patterns were observed. We identified 19 metabolites differentially abundant. Our results suggest a differential reconfiguration of carbon partioning in E. grandis cambial zone. During summer, pyruvate is primarily metabolized via ethanolic fermentation, possibly to regenerate NAD+ for glycolytic ATP production and cellular maintenance. However, in winter there seems to be a metabolic change and we found that some sugars were highly abundant. Our results revealed a dynamic change in E. grandis cambial zone due to seasonality and highlight the importance of glycolysis and ethanolic fermentation for energy generation and maintenance in Eucalyptus, a fast growing tree.
Eucalyptus is the most widely planted hardwood genus in the world. The trees are valued for their fast growth, high adaptability to different climatic conditions and multiple uses of their wood (e.g., pulp and paper industries, charcoal-based steel, wood panels and potential feedstock for lignocellulosic biofuels; Albaugh et al., 2013; Nogueira et al., 2015). Wood formation (xylogenesis) is a complex and highly dynamic process. It is the result of cumulative annual activity of the vascular cambium (Li et al., 2010), a secondary meristem which repeatedly induces its cell division, originating xylem and phloem and is responsible for self-maintenance and signal transfer via translocation of growth regulators (Larson, 1994). The annual course of cambial activity is generally related to the alternation of cold and warm, and/or dry and rainy seasons (Lachaud et al., 1999). The knowledge about cambial activity is fundamental since it is the time in which trees receive environmental signals directly responsible for their growth (Callado et al., 2014). Perennial woody plants from temperate zones have developed mechanisms, that undergo seasonal cycles of activity and dormancy, which are collectively known as annual periodicity (Ko et al., 2011; Begum et al., 2013). This periodicity plays an important role in the formation of wood and reflects the environmental adaptation of woody species. Therefore, the quantity and quality of wood depend on the division of cambial cells and the differentiation of cambial derivatives (Begum et al., 2013). In tropical regions (e.g., Brazil) seasonal changes are less pronounced than in temperate regions, cambial activity is relatively longer and may continue throughout the year (Prislan et al., 2013). It is suggested that in tropical regions water availability is the main factor that induces cambial seasonality. An annual dry season with a length of 2–3 months and less than 60 mm monthly precipitation induces the reduction in cambial activity, which is reestablished in the seasons where monthly precipitations are higher (Worbes, 1995). Oliveira et al. (2012) investigated the relationship between precipitation and wood production in a 23-years old Eucalyptus grandis in Brazil, using x-ray densitometry. The results showed a positive correlation between precipitation data and annual increment of wood. The molecular and physiological mechanisms that enable trees to survive and maintain themselves under limiting conditions, such as winter and limited water availability, are crucial to woody plants (Ko et al., 2011). However, the molecular mechanisms that occur in the cambial zone during seasonal changes are largely unknown. Furthermore, despite the influence that seasonality has on cambial activity, molecular studies regarding the changes that occur in these tissues have only been carried out in species from temperate zones, emphasizing cold acclimation (Schrader et al., 2004; Yang et al., 2004; Gricar and Cufar, 2008; Li et al., 2010; Ko et al., 2011; Galindo-González et al., 2012).
Glycolysis and Tricarboxylic Acid (TCA) cycle are known as central backbones of plant primary metabolism. Under aerobic conditions, pyruvate is transported into mitochondria and oxidized through TCA cycle into organic acids, CO2 and water, via aerobic respiration. During these steps the reducing equivalent, NADH is formed and used by the mitochondrial electron transport chain to power the synthesis of ATP (Fernie et al., 2004). However, under oxygen-limiting conditions (hypoxia), fermentative metabolism is activated to recycle NAD+ from NADH in order to avoid depletion of the cytosolic NAD pool (Zabalza et al., 2009) and to keep glycolysis running in the absence of oxidative phosphorylation by the mitochondrial electron transport chain (van Dongen et al., 2011). Oxygen availability is a key factor in respiration metabolism because it is the final electron acceptor in the oxidative phosphorylation (Moller, 2001). However, under oxygen-limiting conditions (hypoxia), plants have developed alternative pathways to aerobic respiration (e.g., fermentative reactions) in order to maintain glycolysis through the regeneration of NAD+ and thus obtaining the energy necessary to maintain key process of metabolism (membrane integrity, protein synthesis and turnover, regulation on cytosolic pH, among others (Tadege et al., 1999; Rocha et al., 2010). Thus, pyruvate might be converted to ethanol by pyruvate decarboxylase (PDC) and alcohol dehydrogenase (ADH; Zabalza et al., 2009). In trees from tropical areas this hypoxic condition is probably more intense during summer, when trees (e.g., Eucalyptus) are growing fast. Moreover, it is known that the bark has a low permeability to gaseous diffusion (Waisel, 1995; Sorz and Hietz, 2006). The reduced bark permeability may limit the atmospheric oxygen diffusion to the cambium zone, causing a reduction in respiration, by both, limiting the regenerating capacity of NAD+ and increasing the accumulation of CO2 (Pfanz et al., 2002; Aschan and Pfanz, 2003). Such changes in primary metabolism could lead to a reconfiguration of carbon allocation during summer and winter. Celedon et al. (2007) and Carvalho et al. (2008) reported the identification of proteins and transcripts related to carbon metabolism in the cambial zone of E. grandis trees at different ages, including alcohol dehydrogenase and pyruvate decarboxylase. Leonardi et al. (2015) identified stem proteins (e.g., energy metabolism) potentially related to cold stress responses in Eucalyptus urograndis. On the other hand, Durand et al. (2011) found two sucrose synthase more abundant in the cambial zone of Populus tremula L. × P. alba L. submitted to drought stress.
Given the importance of Eucalyptus wood to tropical and subtropical regions, the objective of the present study was to identify changes in transcripts, proteome and metabolome that occur in E. grandis cambial zone during summer (wet season and cambial active growth) and winter (dry season). We hypothesize that during the summer season, the rapid growth and intense metabolic activity in E. grandis trees is compensated by an increase in glycolysis and ethanolic fermentation, in order to maintain the carbon metabolism in the cambial zone. To our knowledge this is the first study of cambial zone metabolism analyzed at three different molecular levels for a tropical and subtropical cultivated tree. This approach allows for a more comprehensive characterization of the metabolic processes that are occurring during summer and winter. We found significant changes in the metabolism of the cambial zone, suggesting that during summer, pyruvate produced by glycolysis is primarily diverted to ethanolic fermentation, to regenerate NAD+ and maintain glycolysis producing ATP. However, during winter we found some sugars and transcripts related to sugar biosynthesis highly abundant. We also hypothesized that these metabolic changes between summer and winter are mainly triggered by higher water availability during the fast growing period contrasting with diminished water in the dry season (winter). Taken together, these results show that different pathways are actively expressed in summer/winter tree and provide new insights into molecular basis of cambial zone differential responses to seasonality in tropical trees.
Materials and Methods
Plant Material and Experimental Conditions
Tissue samples were harvested from clonal trees of 6 years-old E. grandis W. Hill ex Maiden, kindly provided by Suzano Papel e Celulose. The field-trial was situated in Itapetininga, State of São Paulo, Brazil (23°35′20″ S, 48°03′11″ W) at an altitude of 656 m. The cambial zone of the trees was harvested during the summer (26/January/2009), the actively growing season and during the winter (19/August/2009), less-actively growing season. The term cambial zone adopted here describes the entire group of cells which includes cambial initials, xylem and phloem mother cells and developing tissues (differentiating xylem and phloem). To collect the cambial zone, bark of each tree was removed at chest height (1.30 m, exposing an area of ~20 × 15 cm) and the differentiating phloem and xylem were scraped with a razor blade from the inner side of the bark and the outermost side of the stem, respectively, as shown in Supplementary Figure 1. Tissues harvested were immediately frozen in liquid nitrogen. The stem tissue was scraped until the fibrous material below the differentiating cells was reached as described by Celedon et al. (2007) and Carvalho et al. (2008). Samples were harvested in the morning, between 9 a.m. and 10 a.m. The average temperature and precipitation data based on the previous 15 days before sampling were 22.5°C and 99.0 mm, respectively, in summer/January 2009. During the month of winter/August 2009 these parameters were 16.8°C and 15.3 mm, respectively. The field-trial was a completely randomized design. Six bulked samples (10 trees each) were made by grinding and mixing the sampled material. Three bulks represented summer and three represented winter trees.
RNA Extraction and mRNA Isolation
Total RNA was extracted from the cambial zone tissue using the protocol described by Zeng and Yang (2002). Total RNA concentration was measured spectrophotometrically at 260/280 nm, using a U-3300 spectrophotometer (Hitachi, Tokyo, Japan). The absence of RNA degradation was verified by electrophoresis on a formamide-formaldehyde denaturing agarose gel (1%). mRNA was isolated from total RNA (50 μg) using Dynabeads® mRNA purification kit (Invitrogen Dynal, Oslo, Norway), according to the manufacturer's instructions.
Gene-specific primers pairs were designed with the Primer 3 software (Supplementary Table 1). Primer pairs were designed as follows: primer length between 18 and 25 bp, product length of 100–250 bp, melting temperatures 55–60°C, GC% between 40 and 60%. First and second strand cDNA synthesis were performed using the SuperScript™One-Step RT-PCR Platinum® Taq kit (Invitrogen, Carlsbad, CA, USA) with RT/Platinum® Taq (Invitrogen, Carlsbad, CA, USA) and using primers specific for the genes of interest (Supplementary Table 1). The cDNAs were produced in a Gene Amp® PCR System 9700 thermocycler (Applied Biosystems, Foster City, CA, USA) using as annealing temperatures 57°C. The cDNAs were used as a template for RT-qPCR assays, carried out in an iQ5 instrument (BioRad) to obtain de threshold quantification cycle (Cq) and the amplification efficiencies (E). At the end of the PCR cycles, the thermocycler was programmed to perform a denaturation curve. The final volume of each reaction was 20 μL, including cDNA, 10 mM of each primer and 1x Supermix SYBR Green real-time RT-PCR (Invitrogen). The reaction condition was as follows: 95°C for 30 s, 35 cycles of 95°C for 10 s, 57°C for 10 s, and 72°C for 10 s. A negative control (no cDNA template) was included for every gene. All RT-qPCR reactions were carried out in triplicate. The software LinReg (Ramakers et al., 2003) was used to calculate the PCR efficiencies and the Cq values of each gene analyzed. Reference genes were identified using NormFinder (Andersen et al., 2004) which gave the best stability value (0.4) for the combination of alpha-tubulin and malate dehydrogenase (MDH) genes. The calculation of relative expression ratios was carried out with the Relative Expression Software Tool (REST) using the pairwise fixed reallocation randomization test for the statistical significance (P ≤ 0.05; Pfaffl et al., 2002).
Protein Extraction and IEF-SDS PAGE
Total protein of the cambial zone was extracted from 6 grams of frozen ground tissue using the phenolic method according to Hurkman and Tanaka (1986), with minor modifications described in Celedon et al. (2007). The protein pellet was dried under vacuum at 4°C and suspended in 1 mL of solubilization buffer [7 M urea, 2 M thiourea, 0.4% v/v Triton X-100, 50 mM dithiothreitol (DTT)]. Proteins were quantified using the Bradford method (Bradford, 1976). Protein samples (500 μg) were mixed with buffer (340 μL) containing 10 mM DTT, 4% (w/v) CHAPS, 1% IPG buffer (GE Healthcare, Chalfont St. Giles, UK) and 1% (w/v) bromophenol blue and used to rehydrate Immobiline IPG strips (18 cm, 4–7 linear immobilized pH, GE Healthcare, Chalfont St. Giles, UK) for 12 h at 20°C using 50 V. Rehydrated strips were isoelectrofocused (IEF) in an Ettan™IPGphor II™(GE Healthcare) at 100 V for 1 h and then 500 V for 1 h, 1000 V for 1 h, 5000 V for 1 h, and 8000 until reaching a total of 80000 V-h. After IEF, the strips were kept at −80°C until needed. Before the second dimension, strips were kept at room temperature for 15 min in equilibration buffer (6 M urea, 2% w/v SDS, 50 mM Tris-HCl, pH 6.8, 30% v/v glycerol) firstly, with 1% w/v DTT and then with 2.5% w/v iodoacetamide (IAA) and 0.001% bromophenol blue. The second dimension was performed in 12% (w/v) polyacrylamide gels, using Protean II XI 2-D cell electrophoresis system (GE Healthcare), at 30 mA per gel. Three biological replicates were performed for each treatment. Proteins were detected using Coomassie Brilliant Blue G-250 (Candiano et al., 2004), with modifications. Gels were incubated for 1 h in a solution containing 40% (v/v) ethanol and 10% (v/v) acetic acid, in water. For protein detection the gels were left overnight in staining solution (20% (v/v) methanol, 10% (w/v) ammonium sulfate, 10% v/v phosphoric acid, and 0.1% (w/v) Coomassie G-250). After three washes in water (2 h each), the gels were stored in 5% (w/v) ammonium sulfate until image analysis, spot selection and picking.
Gel Image Analysis
Gels were imaged using Image scanner III and Labscan v 7.0 software (GE Healthcare). Image analysis was performed automatically using the Image Master 2D Platinum software v 7.0 (GE Healthcare). Gels from three independent biological replicates were used. Spots were detected using a smoothness of 8, minimum area of 15 and a saliency of 40, and spots across gels were matched using 5 landmarks per gel. Matching was performed automatically and systematically confirmed after one-by-one visual checking: artifacts, or spots that could not be confidently validated as true matches were disregarded and misalignments were corrected manually when appropriate. Spots were considered reproducible when well-resolved in at least two of the three biological replicates. The abundance of each protein spot was estimated by the percentage of volume (% vol). For each treatment analyzed the average spot volume of the three replicate gels was determined, followed by normalization (individual spot volume/total spot volume × 100). The normalized volumes (% vol.) of the corresponding spots from summer and winter samples were compared to estimate differential expression of proteins during different seasons. To ensure the reproducibility between biological replicates, spots with coefficient of variation higher than 30% were excluded. To identify spots that changed in abundance between the two seasons, the data collected from protein spots volumes were subjected to Students t-test (P ≤ 0.05) in Image Master v 7.0 software. The fold-change (summer/winter) value of each identified spot were also given (Table 1). However, no fold-change cut-off was established and their values were not taken into account in data discussion.
Protein Identification by Mass Spectrometry
In-gel digestion of proteins was done as described in Celedon et al. (2007). After, peptide mixtures were sequenced by online chromatography using a nano-Acquity UPLC (Waters®) system coupled to a Q-TOF Ultima API mass spectrometer (Waters, UK). Mass spectrometer parameters were set according to Celedon et al. (2007). Ten microliters of sample were loaded onto a pre-column Symmetry C18 5 μm, 5 × 30 mm (Waters) for sample preconcentration and desalination, followed by peptide separation on an LC column Symmetry C18, 5 μm, 32 × 150 mm (Waters). Peptides were eluted using a10–45% linear gradient of solvent B [95% (v/v) ACN, 0.1% (v/v) FA in water], starting at a flow rate of 5 μL/min for the first 15 min, then changing to 2 μL/min for the next 25 min, and back to 5 μL/min in the last 5 min. Solvent A consisted of 5% v/v ACN, and 0.1% v/v FA in water. All analyses were performed using a positive ion mode at 3 kV needle voltage. The mass range was set from m/z 300 to 2000, and the MS/MS spectra were acquired for the most intense peaks having at least 15 counts. The LC-MS/MS were processed using ProteinLynx v 2.0 (Waters) and Mascot Daemon software (Matrix Science, Boston, MA), and the sequences searched against an in-house Eucalyptus database from Phytozome v1.1 (www.phytozome.net/eucalyptus.php) and NCBI. Combined MS/MS search criteria used were as follows: trypsin digestion; fixed modification set as carbamidomethylation of cysteine; variable modification set as methionine oxidation; mass accuracy of 50 ppm for the parent ion and MS/MS mass tolerance of 0.1 Da. According to MASCOT probability analysis, only significant hits (P ≤ 0.05) were accepted. A match was considered significant if the peptide had a score higher than 70, based on Perkins et al. (1999). The mass spectrometry proteomics data have been deposited to the ProteomeXchangeConsortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository (Vizcaino et al., 2013) with the dataset identifier (PXD003504).
Metabolite Extraction, Data Processing, and Statistical Analysis
The cambial zone tissue of each bulk was ground into powder in liquid N2 and dried for 24 h using the Modulyod-320 freeze dry system (Thermo Scientific). Metabolites were extracted from powder tissue and analyzed according to the methods described by Hoffman et al. (2010), with minor changes. Approximately 5 mg of dried tissue was mixed with 1 mL of a chloroform-methanol-water mix (6:2:2) containing stable isotope reference compounds [15 ng mL−1 each of (13C3)-myristic acid, (13C4)-hexadecanoic acid, (2H4)-succinic acid, (13C5, 15N)-glutamic acid, (2H7)-cholesterol, (13C5)-proline, (13C4)-disodium-ketoglutarate (13C12)-sucrose, (2H4)-putrescine, (2H6)-salicylic acid and (13C6)-glucose)]. The metabolite extraction proceeded using a vibration mill set to a frequency of 30 Hz s−1 for 3 min, with 3 mm tungsten carbide beads added to each extraction tube to increase the extraction efficiency. The extracts were then centrifuged for 10 min at 14.000 xg in an Eppendorf centrifuge (model 54178). After, 100 μL of each supernatant was transferred to a GC-vial and evaporated to dryness. The samples were then derivatized with 30 μL of methoxyamine hydrochloride (15 mg mL−1) in pyridine for 16 h at room temperature. Trimethylsilation was performed by adding 30 μL of n-methyl-n-(trimethylsilyl) trifluoroacetamide (MSTFA) with 1% TMCS to the samples and incubating them for 1 h at room temperature. After silylation, 30 μL of heptane was added. Samples were analyzed, according to Gullberg et al. (2004), using gas chromatography with time-of-flight mass spectrometry (GC-TOF/MS) together with blank control samples and a series of n-alkanes (C12–C40), which allowed retention indices to be calculated (Schauer et al., 2005). Three biological replicates with three technical replicates were used for each treatment. One microliter of each derivatized sample was injected splitless by a CTC Combi Pal Xt Duo autosampler (CTC Analytics AG, Switzerland) into an Agilent 7890A gas chromatograph equipped with a 30 m × 0.25 mm i.d. fused-silica capillary column with a chemically bonded 0.25 μm DB 5-MS UI stationary phase (J&W Scientific, Folsom, CA). The injector temperature was 260°C, the septum purge flowrate was 20 mL min−1 and the purge was turned on after 75 s. The gas flow rate through the column was 1 mL min−1, the column temperature was held at 70°C for 2 min, then increased by 15°C min−1 to 320°C, and held there for 4 min. The column effluent was introduced into the ion source of a Pegasus HT time-of-flight mass spectrometer (Leco Corporation, St. Joseph, MI, USA). The transfer line and the ion source temperatures were 250 and 200°C, respectively. Ions were generated by a 70 eV electron beam at an ionization current of 2.0 mA, and 20-30 spectra s−1 (30 spectras−1 run 1, 20 spectra s−1 run 2) were recorded in the mass range m/z 50−800. The acceleration voltage was turned on after a solvent delay of 290 s. The detector voltage was 1450-1490 V (1450 V run 1, 1490 V run 2). All non-processed MS-files from the metabolic analysis were exported into Chroma TOF 2.12 software (Leco Corporation), in which all manual integrations and metabolite identification were done. The identities of the compounds were determined by performing database searches based on the mass spectra and the compounds retention indices (RIs). The databases used were the NIST mass spectra library, an in-house database established by Umeå Plant Science Centre (UPSC), and the mass spectra library maintained by the Max Planck Institute in Golm, Germany. All data treatment procedures (smoothing, baseline correction and chromatogram alignment) were performed using custom scripts (Jonsson et al., 2005) in MATLAB. To compare the metabolite changes between summer and winter seasons the normalized data set (to tissue dry weight and internal standards) was Pareto scaled, log transformed and applied to multivariate and univariate analytical methods using the MetaboAnalyst software (Xia et al., 2012). The multivariate supervised classification method PLS-DA (Partial least squares-discriminant analysis) was carried out to discriminate between different groups (summer and winter). The variable importance in the projection (VIP) was used to rank the predictor variables, the metabolites, according to their contribution to the response in the respective PLS-DA model (Korn et al., 2010). VIP is a weighted sum of squares of the PLS loadings that takes into account the amount of explained Y-variance of each component (Xia and Wishart, 2011). PLS-DA model fit was evaluated using the R2 and Q2 cross-validation performance measures (Xia et al., 2012), both of which vary between 0 and 1. R2, the squared correlation coefficient between the dependent variable and the PLS-DA prediction, provides an indication of the “goodness of fit” (a value between zero and one, where one is a perfect correlation) from the model. Q2 provides an indication of “goodness-of-prediction” and is the averaged correlation coefficient between the dependent variable and the PLS-DA predictions. Differential metabolites were selected from the PLS-DA model using a combination of VIP value >1 and P ≤ 0.05, by the univariate unpaired, two-tailed Student's t-test.
Carbohydrate Extraction and HPLC Analysis
Glucose, sucrose and fructose were extracted from cambial zone samples, three biological replicates with three technical replicates were used for each treatment. Tissues were grounded and freeze dried for 48 h, after 1 mL of water was added in 0.2 g of dry powder and samples were kept in bath (80°C) for 1 h. Then, samples were centrifuged for 10 min, 16.000 xg, the supernatant was recovered and stored at −4°C. Summer and winter samples were analyzed using a high-performance liquid chromatography (ICS 2500, HPLC Dionex) with amperometric detection (ED50) equipped with an autosampler AS50 (Dionex). Sugars were assigned according to the retention times of standards (sucrose, glucose and fructose). A Carbopac PA-1 column (4 × 250 mm, Dionex) and a guard Carbopac PA-10 column (4 × 50 mm, Dionex) were used. To identify statistical differences between summer and winter samples two-tailed Student's t-test (P ≤ 0.05) was performed.
Changes in mRNA Levels during Summer and Winter Seasons
As a first step to understand the molecular mechanisms underlying the seasonal reorganization in the cambial zone, we examined the transcriptional pattern of 26 genes (Figure 1 and Supplementary Table 1) related to specific pathways of plant primary metabolism: glycolysis, ethanolic fermentation, TCA cycle and carbon fixation. These pathways were chosen because they are crucial to plant respiration, maintenance and ATP generation. Of the 14 genes analyzed from glycolysis, three were up-regulated in summer (PGK, phosphoglycerate kinase; PK, pyruvate kinase; PGAM, phosphoglycerate mutase) and five were up-regulated during winter (GPI, glucose-6-phosphate isomerase; GAPDH, glyceraldehyde-3-phosphate; PGM, phosphoglucomutase; Susy1, sucrose synthase 1; PFP, PPi-dependent phosphofructokinase; Figure 1A). PGK and PGAM genes catalyze subsequent reversible reactions in glycolysis while PK catalyzes an irreversible reaction. Among the five highly expressed genes in winter, four (SUSY1, GPI, PFP, and GAPDH) catalyze reversible reactions in glycolysis, whereas PGM catalyzes an irreversible reaction. PFP, GPI, and PGM participate in subsequent steps in glycolysis and were all up-regulated in winter. Considering the genes participating in ethanolic fermentation we analyzed three alcohol dehydrogenases (ADH) alternative transcripts, all of them from ADH class I, and one piruvate descarboxylase (PDC). Only ADH3 and PDC were statistically significant and both were up-regulated in summer (Figure 1B). In addition, the transcriptional levels of two TCA cycle genes (SCL, succinyl CoA ligase; MDHm, malate dehydrogenase) were analyzed (Figure 1C). SCL and MDHm showed opposing expression patterns, SCL was up-regulated in summer while MDHm was up-regulated in winter. Six genes involved in carbon fixation were analyzed (Figure 1D), with three being up-regulated during winter (CA, carbonic anhydrase; RbcL, rubisco large subunit; FBAcl, fructose bisphosphate aldolase).
Figure 1. Comparison of transcripts abundances of genes involved in primary metabolism: (A) glycolysis, (B) Alcoholic fermentation, (C) TCA cycle, and (D) Photosynthesis, in E. grandis cambial zone by RT-qPCR. Data are expressed as log and winter values were used as a control. Expression was determined relative to alpha-tubulin and MDH, as described in Materials and Methods. Asterisks indicates genes that are significantly expressed (P ≤ 0.05). Abbreviations: ALDO, fructose bisphosphate aldolase; GPI, glucose 6 phosphate isomerase; PGK, phosphoglycerate kinase; GAPDH, glyceraldehyde 3 phosphate dehydrogenase; PK, pyruvate kinase; PEPC, phosphoenolpyruvate carboxylase; PFK, ATP-dependent phosphofructokinase; ENO, enolase; PGM, phosphoglucomutase; PGAM, phosphoglycerate mutase; PDH, pyruvate dehydrogenase; SuSy1, sucrose synthase 1; SuSy3, sucrose synthase 3; PFP, PPi-dependent phosphofructokinase; ADH1, alcohol dehydrogenase 1; ADH2, alcohol dehydrogenase 2; ADH3, alcohol dehydrogenase 3; PDC, pyruvate decarboxylase; SCL, succinyl-coa ligase; MDHm, malate dehydrogenase mitochondrial; NADPme, NADP-malic enzyme; CA, carbonic anhydrase; RbcS, rubisco small subunit; RbcL, rubisco large subunit; FBAcl, fructose bisphosphate aldolase chloroplastidial; RPI, ribose-5-phosphateisomerase. Three biological replicates, each with three technical replicates were analyzed per sample and error bars are standard errors of mean.
Protein Profile in the Cambial Zone during Summer and Winter Seasons
Protein profile of cambial zone (summer vs. winter) were characterized using IEF-SDS-PAGE in conjunction with LC-MS/MS (Supplementary Figure 2). More than 400 proteins were detected by Image Master Platinum (GE Healthcare). Of these, 129 spots had significant (P ≤ 0.05) and reproducible changes in abundance between samples harvested in summer and winter. Among these proteins, ~60% (70 spots) were successfully identified (Table 1). From these identified proteins spots, 74.3% (52 spots) were up-regulated in summer and 25.7% (18 spots) in winter. The others proteins spots were disregarded because they did not show similarities when queried against Eucalyptus and NCBI databases or had scores lower than 70 (Perkins et al., 1999). Many analyzed peptides yielded hits from the same protein accession. Similar observations were made in other studies which provided possible explanations such as protein degradation and post-translational modification (Jiang et al., 2007; Shu et al., 2011; Galindo-González et al., 2012). The proteins identified were classified into six main biological categories (Figure 2 and Table 1) as described by Rison et al. (2000).
Figure 2. Categorization of proteins identified in E. grandis cambial zone during summer and winter.
Proteins in the category of Metabolism and Energy were the most abundant (42%) and were distributed into five classes: Energetic Metabolism, Energy Transfer, Amino acid Metabolism, Nucleotide Metabolism and Secondary Metabolism. This category included mainly proteins related to glycolysis (aldolase, PGI, PGK, enolase, and GAPDH), ethanolic fermentation (ADH and PDC) and TCA cycle (citrate synthase, NADP-isocitrate dehydrogenase (NADP-IDH) and Pyruvate dehydrogenase). Most of these proteins were analyzed by RT-qPCR, showing similar expression patterns in many cases. The identification of proteins related to glycolysis, ethanolic fermentation and TCA cycle is interesting, because they are part of an interconnected pathway fulfilling the cambial region needs for ATP, NAD, NADH, and pyruvate for mitochondrial respiration (Celedon et al., 2007). In addition, three ADH isoforms were found, with ADH1 and 2 (spots 5 and 17) being highly expressed in summer, while ADH3 (spot 48) was highly expressed in winter. Another important protein identified was UGPase, for which two isoforms (spots 12 and 25) were mainly expressed in summer.
Contrary to results obtained with gene expression, we did not find proteins from carbon fixation pathways. This can be explained by limitations in 2-DE technique. The second most abundant category was the Information Pathway (39%), mainly represented by nine heat-shock proteins isoforms (six highly expressed in summer), besides two TCP-1/cpn60 chaperones (with opposite pattern of expression), two HSP20-like chaperones (highly expressed in winter) and three ATPases (highly expressed in summer). These are stress response proteins, acting in protein folding, assembly, translocation and degradation, maintenance of proteins functional conformation and preventing the aggregation of non-native proteins (Wang et al., 2004). Structure and Organization category was represented by proteins involved in the lignification process [caffeic acid O-methyltransferase (COMT) and caffeoyl CoA O-methyltransferase (CCoAOMT)] and proteins related to cytoskeleton formation (two alpha tubulins and one tubulin folding factor), all of them up-regulated in summer.
Soluble Sugars Availability and Metabolic Changes in Cambial Zone during Summer and Winter
To investigate the soluble sugar accumulation in response to seasonal changes we measured glucose, fructose and sucrose contents from summer and winter samples. The concentrations of all sugars were similar in both seasons. However, sucrose showed a higher concentration, compared to glucose and fructose, probably because it is the main sugar translocated in plants (Geiger, 2011). GC/TOF-MS was performed to identify changes in metabolites levels in the cambial zone due to seasonal variation. A total of 53 metabolites were identified and any differences between summer and winter samples were explored by multivariate and univariate methods (Supplementary Table 2). A PLS-DA model with 5 components was conducted to reveal the differences between the two groups. Based on this we found a clear separation between summer and winter groups (Figure 3). The two first principal components from PLS-DA accounted for 60% of the total variance and the model showed a significant R2 and Q2 values of 0.92 and 0.95, respectively. Based on VIPS >1 and P < 0.05 we found 19 differentially abundant metabolites between the summer and the winter groups (Table 2). The 19 metabolites were classified into five categories: amino acid (21%), fatty acid (5, 26%), organic acid (31, 57%), sugar (36, 84%), and unknown (5, 26%; Table 2). Most of the metabolites identified (68, 42%) were more abundant in summer samples, including all the amino acids.
Figure 3. PLS-DA score plots of eucalyptus cambial zone showing a significant separation (R2 = 0.99 and Q2 = 0.98) between summer (square) and winter (circle) samples. The x-axis represented the modeled co-variation and the y-axis is the predictive component. Component 1 and 2 contributes with 43.4 and 16.89%, respectively, of the total variance. Three biological replicates, each with three technical replicates were analyzed per sample.
Growth is controlled by both, endogenous factors, such as plant genotype or physiological processes (Schrader et al., 2004), and by exogenous factors, such as temperature, light and water availability (Deslauriers and Morin, 2005). Regular fluctuations of these limiting factors promote a periodic rhythm in plant development, reproduced at the cambium level by an active and a dormant stage, even in tropical trees (Callado et al., 2013). Cambial activity is genetically controlled but the rhythm of new cells production is determined by environmental limitations (Kozlowski and Pallardy, 1997). However, cambial responses and variations in wood properties are often difficult to interpret and generalize, due to the complexity and co-variance of factors involved in controlling xylem development, including tree hormonal balance, carbon fixation and allocation and tree water status (Chaffey et al., 2002). Water availability is the major climate driver of wood production. As a consequence, months with lowest water availability are less favorable for cell expansion and can affect cell wall thickness thus, properties such as wood density may be markedly altered (Balducci et al., 2014; Wagner et al., 2016). In the current work we aimed to identify molecular changes in the cambial zone of E. grandis due to seasonal variations (summer and winter) observed trees growing in the southeast of Brazil. Summer is characterized as a raining season, in which cambium is actively growing. On the other hand, winter is a dry season and this low water availability is less favorable for cell expansion (Wagner et al., 2016). We hypothesized that the intense metabolic activity in cambium, during summer, is probably maintained by glycolysis and ethanolic fermentation because of the intense aerobic respiration that might deplete cambial O2 supplies. Thus, alcoholic fermentation is required to generate NAD+ so that glycolysis can produce limited amount of ATP for cellular maintenance until the re-establishment of oxidative phosphorylation (Drew, 1997). We also hypothesized that cambium seasonal variations observed in our work were mainly trigged by water availability in soil instead of low temperatures. We found that metabolic pathways related to primary metabolism are differentially regulated in summer and winter, resulting in different carbon partitioning in the two seasons. As we expected we found transcripts and proteins from alcoholic fermentation pathway highly abundant in summer, the season which coincides with the rainfall period. In this season, with plenty of available water, the tree is growing rapidly, the cambium zone is metabolically active and cell division must be occurring in high rates. Thus, NAD+ must be continuously regenerated by fermentative reactions, maintaining the primary metabolism in the cambial zone, corroborating our hypothesis. On the other hand, during winter, when it is dry, tree's growth rates diminished and cell walls are ticker, the flux of carbon seems to favoring carbohydrate accumulation, which will probably be required for the next growing season.
Effects of Seasonal Variation in the Gene Regulation of Metabolism, Influencing Carbon Partioning in the Cambial Zone of E. grandis
Transcriptomic analyzes of different Eucalyptus tissues have been reported in an attempt to understand wood quality (Salazar et al., 2011, 2013), cold acclimation (Liu et al., 2014), frost tolerance (Fernández et al., 2015), and water stress (Thumma et al., 2012). In this study we analyzed a set of transcripts related to primary metabolism by RT-qPCR. We found differential transcription patterns for glycolytic genes in response to summer and winter seasons (Figure 1A). During winter the up-regulation of PFP, GPI, and PGM transcripts, which catalyze sequential step in glycolysis, might be indicating that the C flux is driving toward glucose-6-phosphate formation (G-6P), while the up-regulation of SUSY1 and GAPDH indicates that both enzymes are working to promote the formation of UDP-glucose (UDP-G) and glyceraldehyde-3-P (G-3P), respectively. In the summer condition, the up-regulation of PGK, PGAM, and PK suggests that C flux goes in the direction of pyruvate formation. As mentioned before, the final destination of the pyruvate formed in glycolysis, depends on the oxygen availability. In the presence of O2, pyruvate is completely oxidized to CO2 and water via aerobic respiration (Kimmerer and Stringer, 1988). However, when oxygen availability decreases, below the level at which oxygen becomes limiting for oxidative phosphorylation in TCA cycle, ethanolic fermentation is required. In our work ADH3 and PDC were up-regulated in the cambial zone during summer (Figure 1B) indicating that fermentation might be high in order to supply glycolysis with NAD+ to maintain the supply of ATP, as we hypothesized. In addition, this energy produced might be utilized for cell division and biosynthesis of new cell wall in the cambium. Despite the importance of ethanolic fermentation for plant metabolism under O2 limitation, the TCA cycle is essential to maintain respiration even under anaerobic condition. The first study about alcohol dehydrogenase in tree stems were done by Kimmerer and Stringer (1988). They hypothesized that tree with thick bark would have greater ADH activity and ethanol content than thinner-barked trees of the same species, however, they did not observe any relation between bark thickness and ADH activity.
We found two transcripts of the TCA cycle (SCL and MDHm) differentially expressed and they showed opposite expression patterns (Figure 1C). The TCA cycle is a universal feature of the aerobic organism metabolism, thus eucalyptus probably has alternative transcripts of TCA cycle enzymes that are seasonally regulated. In RNA-seq experiments, Xu et al. (2014) found that 20.7% of highly expressed transcripts in the developing xylem of Eucalyptus were affected by alternative splicing events. More detailed studies are necessary to understand how seasonal changes regulate this cycle in woody plants. Glycolysis and TCA cycle pathways have several bypass reactions, allowing carbon flux to take different routes. The intermediates of central C metabolism also serve as precursor for synthesis of a plethora of primary and secondary metabolites (Plaxton and Podesta, 2006). Our results suggest that at the transcriptional level, cambial glycolysis and TCA cycle were differentially regulated, resulting in different carbon partitioning between the two seasons.
Changes in Proteins Mainly Related to Primary Metabolism in the Cambial Zone during Summer and Winter Seasons
Through changes in protein profile, based on Student's t-test analysis (P ≤ 0.05), between summer and winter, we found five proteins from the glycolytic pathway (GPI, FBA, GAPDH, PGK and ENO), two proteins involved in ethanolic fermentation (ADH and PDC) and three proteins from TCA cycle (NADP - isocitrate dehydrogenase, citrate synthase and pyruvate dehydrogenase). Not all expression patterns of proteins were correlated to gene expression. These data suggest that in the cambial zone, a complex regulatory network occurred due to seasonal variation also at proteomic level.
In accordance with gene expression, during summer the cambial zone favors pyruvate production, which is distributed mainly to ethanolic fermentation. We found three ADH isoforms; ADH 1 and ADH 2 were highly abundant in summer while ADH 3 was highly abundant in winter. These ADH expression pattern suggest a temporal expression of the isoforms in the cambial zone and need to be further characterized. Another interesting finding were enzymes related to carbon metabolism, two UGPase isoforms, both highly expressed in summer. UGPase activity is related to cell wall synthesis and it promotes carbohydrate accumulation in storage tissues under low temperatures (Kleczkowski et al., 2004). Proteins of carbohydrate and energy metabolism have been reported in the stems of white spruce (Galindo-González et al., 2012), Populus (Plomion et al., 2006; Kalluri et al., 2009; Durand et al., 2011), Eucalyptus (Celedon et al., 2007; Leonardi et al., 2015), and Pinus (Gion et al., 2005). In summer samples, we also identified an upregulation of proteins involved in cell wall synthesis and cytoskeleton (Table 1). These data indicate that cell wall synthesis and deposition occurs during the growing season. Thus, remodeling of the cell wall could facilitate resilience against water deficit and low temperature conditions (Galindo-González et al., 2012).
Seasonal Variation of Photosynthetic Transcripts and Proteins in E. grandis Cambial Zone
Although green leaves are regarded as the main source of photosynthates, the stem can be photosynthetically active (Wiebe, 1975; Pfanz et al., 2002; Aschan and Pfanz, 2003; Cernusak and Cheesman, 2015; Vandegehuchte et al., 2015). Eucalyptus and other trees experience the stem recycling photosynthesis (Ávila et al., 2014), in which the main function is to recycle some portion of the CO2 produced inside the stem as a result of respiratory activity in the woody tissue (Saveyn et al., 2010; Cernusak and Cheesman, 2015). Some authors have been investigating this issue from an ecological and physiological perspective (Berveiller et al., 2007; Eyles et al., 2009; Saveyn et al., 2010; Cernusak and Hutley, 2011), nevertheless photosynthesis in stem is receiving little attention at the molecular level. We checked the expression of six genes involved in carbon fixation and found three genes up-regulated in winter: CA, RbcL, and FBAcl (Figure 1D). Transcripts related to chloroplasts or photosynthetic activity, including RbcS, were found up-regulated in ray cambial cells of hybrid poplar (Populus trichocarpa T.&G. × Populus deltoïdes Marsh var. Boelare) compared to fusiform cambial cells (Goué et al., 2008). CA, RbcL and RbcS transcripts were previously identified in the cambial region of poplar trees (Hertzberg et al., 2001). It is interesting to note that we found transcripts of RbcL and RbcS, however, only RbcL expression was up-regulated in winter, suggesting that this subunit of Rubisco may be temporally regulated being highly required during winter. Rubisco is also considered to be an unconventional storage protein (Cooke and Weih, 2005). The up-regulation of RbcL may also indicate that perhaps, this subunit could play a storage function during winter, once trees accumulate carbon and nitrogen reserves that can be mobilized upon the next growing season. Similar result was observed by Galindo-González et al. (2012).
Contrary to our gene expression results, we did not find proteins from carbon fixation pathways. This can be explained by limitations in 2-DE technique such as, low detection sensitivity and linearity, limited loading capacity of gradient pH strips, poor solubility of membrane proteins among others (Monteoliva and Albar, 2004). CO2 fixation proteins such as sedoheptulose-bisphosphatase, ribose-5-phosphate isomerase and rubisco (large and small subunit) were found in Picea glauca stem (Galindo-González et al., 2012). Rubisco large subunit was found in eucalyptus cambium (Celedon et al., 2007) and both subunits were found in Populus cambial zone (Juan et al., 2006). Our findings are in accordance with these observations, giving molecular support for the concept that photosynthesis occurs in the stem. Nevertheless, more detailed studies are necessary to understand the role of photosynthetic transcripts and proteins in the cambial zone.
Metabolic Profiling Changes in Cambial Zone during Summer and Winter
Most of the differentially abundant metabolites were more abundant in the summer season. This is the case for the amino acids glycine, serine, valine and threonine. In plants it has been proposed that the cytosolic conversion of serine into glycine, could be an important one-carbon source in plant C1 metabolism. It is possible that this metabolism is required during xylem maturation for the production of glycine, which is abundant in the cell wall proteins (Moreau et al., 2005). Oleic acid, was one of the three free fatty acids (FA) found. This saturated metabolite and others free FA were identified in the sapwood of Pinus sylvelstris, their levels were greatest in the beginning and end of the growing season (Saranpaa and Nyberg, 1987).
Xylose, xylitol and erythrose-4-phosphate (E4P) were abundantly found in summer. Xylitol, also called wood sugar, is made from xylose, which is found in the cell walls of most land plants (Nigam and Singh, 1995). Both sugars are important for the food and beverage industries (Nigam and Singh, 1995; Prakasham et al., 2009). We believed that xylose and xylitol were highly abundant in summer due to cell wall division and expansion during the tree active growth. Ko et al. (2011) found that xylose content increased in winter/dormancy stems of Populus trees. The E4P is an intermediate product of glucose oxidation through the pentose phosphate pathway (Bochkov et al., 2012). E4P and phosphoenolpyruvate (PEP) are precursors of the shikimate pathway (Zulak et al., 2008). This pathway provides an alternative route to aromatic compounds (Dewick, 2002) and is responsible to link the lignin biosynthesis to primary metabolism (Boerjan et al., 2003). Besides, the shikimate pathway is also related to the tree active growth (Larisch et al., 2012). Quinic acid, an organic acid precursor of aromatic amino acids (Minamikawa, 1976), was also highly abundant in summer. We also found the shikimic acid metabolite, but it did not show statistical differences between seasons. The organic acid metabolism is fundamental for several biochemical pathways, including energy production, formation of precursors for amino acid biosynthesis and at the whole plant level in modulating adaptation to the environment (López-Bucio et al., 2000). The organic compounds fumarate, maleic acid and threonic acid-1,4-lactone were more abundant during summer. Fumarate is a component of the TCA cycle and can be metabolized to yield energy and carbon skeletons of production of other compounds (Chia et al., 2000), its abundance in summer indicates that TCA cycle is also being required in addition to ethanolic fermentation. Pipecolic acid was the unique organic acid highly abundant during winter, it is a non-protein amino acid, homologous to proline that results from the lysine degradation, in plants. Pipecolic acid accumulation in higher plants seems to be closely related to metabolic responses induced by environmental stress (Moulin et al., 2002).
The three soluble sugars (glucose, sucrose and fructose) measured by HPLC were not statistically different between summer and winter. Since the cambial tissue is not recognized as a storage organ, carbohydrates present in this tissue are being supplied from the bark, which is the storage compartment. Similar results were reported in Populus trichocarpa (Ko et al., 2011). The authors measured soluble sugars from winter/dormancy and summer stem samples and observed a significant increase in total sugars content in the winter/dormancy stem. But, the concentrations of sucrose, fructose and glucose were similar in the two stem types from different seasons.
However, sugars identified by GC-MS were the metabolite category which was the most abundant in winter, having interesting metabolites such as glyceric acid, raffinose, and galactinol. Carbohydrate pools are mobilized in trees to supply respiration during low photosynthesis (Ögren, 2000), leaf development, bud-break, early spring growth (Barbaroux and Breda, 2002), when trees experiences stress such as water stress and pest damage (Canham et al., 1999). Different authors associate the carbohydrate accumulation with seasonal changes and physiological mechanisms in response to cold acclimation (Travert et al., 1997; Welling and Palva, 2006; Bonhomme et al., 2009; Turhan and Ergin, 2012). Raffinose and galactinol are known to act as cryoprotectors and were previously identified in the cambial region, accumulating in response to low temperatures (Guy et al., 2008; Ko et al., 2011; Plavcová et al., 2013). Raffinose levels increased in the cambium of mature conifer trees (Picea abies Karst and Larix decidua Mil) during winter, protecting cell membranes from damage during frost-induced dehydration by detoxifying radical oxygen species (ROS), that accumulates at low temperatures (Simard et al., 2013). On the other hand, both sugars have been hypothesized to be osmoprotectants in drought-stress conditions, and have frequently been implicated in drought response in plants (Taji et al., 2002; Nishizawa et al., 2008). Thus, we believed that in our work raffinose and galactinol might be accumulated during winter in response to diminished water availability instead of low temperatures. Raffinose and galactinol were some of the most highly accumulated metabolites in response to water-deficit conditions in Populus balsamifera leaves (Hamanishi et al., 2015). Deslauriers et al. (2014) found raffinose in low concentrations in the cambium and xylem of black spruce during the growing season. However, the authors reported an increase in raffinose levels in trees under water stress. All these results indicate a complex and dynamic metabolic change in cambial zone.
Overview of the E. grandis Primary Metabolism in the Cambial Zone
Carbon partitioning is a fundamental process in plant physiology. At the most basic level, plants assimilate inorganic carbon and nitrogen into reduced compounds required for growth. At the molecular level the variety of produced metabolites is large and the pathways that they feed into are complex and interconnected (Tegeder and Weber, 2007). A more comprehensive knowledge about a biological process can be achieved by the analysis of transcripts, proteins and metabolites profiles. This could improve and facilitate the characterization of changes in levels of important compounds during cellular regulation. Using data from multiple platforms Srivastava et al. (2013) presented a system response of Populus to oxidative stress integrating data from analysis of the cambial region of wild-type controls and trees expressing high-isoeletric-point superoxide dismutase. However, most of the studies have focused on acquiring and integrating data at two omics levels (Alvarez et al., 2008; Hoffman et al., 2010; Galindo-González et al., 2012; Hamanishi et al., 2015; Tuttle et al., 2015; Yu et al., 2015). To have a better understanding about the dynamic changes occurring in the E. grandis cambial zone, due to summer and winter seasons, the RT-qPCR, proteomics and metabolomics data related to carbon metabolism were grouped (Figure 4). We observed differential abundance of molecules related to different pathways/metabolisms such as sugar, glycolysis, TCA cycle, amino acids, phenylpropanoids. We found significant changes in all levels (transcripts, proteins and metabolites), indicating a clear occurrence of a metabolic reconfiguration triggered by seasonal variation. This tight correlation between different levels of response in primary metabolism implies that transcriptional regulation plays an important role in carbon metabolism. The transcriptional and proteomic levels suggest different carbon partitioning between the two seasons. In summer, metabolism seems to favor pyruvate formation, which is metabolized by the fermentative enzymes. Thus, energy and reducing power are produced to supply the energetic demands required during the high metabolic activity of the season. Transcripts accumulation of genes involved in energy metabolism, cell division, differentiation and expansion were observed in the xylem tissues of Pinus radiata during its greater growth rate (Li et al., 2010). In this study, we showed a reconfiguration in the cambial zone metabolism of E. grandis due to seasonal changes in the tropic region. Although in the tropics plants do not experience extremely low temperatures similarly to those of temperate regions, our data indicate a reconfiguration in carbon partitioning from summer to winter. Different metabolic pathways were favored in each season. With this comprehensive dataset in hand, we can now explore the role of new genes/protein isoforms, mainly related to glycolysis and ethanolic fermentation, to study carbon partitioning during seasonal variation in trees.
Figure 4. Differentially abundant transcripts, proteins and metabolites involved in the cambial zone primary metabolism. Transcripts are indicated by italic names and proteins are indicated by name. Subscribed letters S (summer) or W (winter) just below transcripts or proteins names indicates the season in which it was highly expressed. Asterisks in transcripts names indicates that the correspondent protein was also found as differentially expressed. Metabolites are indicated by boxes, metabolites in dashed boxes were detected in our work but were not significantly affected by seasonal changes (summer/winter). Metabolites highly abundant in summer or winter are indicates by white and gray boxes, respectively. RbcL, rubisco large subunit; FBAcl, fructose bisphosphate aldolase chloroplastidial; SuSy, sucrose synthase; GPI, glucose 6 phosphate isomerase; PGM, phosphoglucomutase; PFP, PPi-dependent phosphofructokinase; FBA, fructose bisphosphate aldolase citoplasmatic; GAPDH, glyceraldehyde 3 phosphate dehydrogenase; PGK, phosphoglycerate kinase; PGAM, phosphoglycerate mutase; ENO, enolase; PK, pyruvate kinase; PDC, pyruvate decarboxylase; ADH, alcohol dehydrogenase; PDH, pyruvate dehydrogenase; SCL, succinyl-coa ligase; MDHm, malate dehydrogenase mitochondrial; ATP-Citrate synthase.
IB conceived of the study, carried out the experiments, analyzed the data and wrote the manuscript. DM assisted in the RT-qPCR analysis and manuscript revision. PL and TM assisted with metabolomics studies. CL conceived of the study, participated in its design and coordination, and reviewed the manuscript. All author's approved the final manuscript.
Conflict of Interest Statement
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.
This work was supported by the FAPESP (São Paulo Research Foundation) (Grant 2009/00808-2). We would like to thank Dr. Esteban Gonzalez (FuturaGene Corporation) for providing the biological material. We also would like to thank José Mateus Wisniewski Gonsalves (FuturaGene) and Dr. Luis Felipe Boaretto for the help with the sample harvesting.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00932
Supplementary Figure 1. Bark was removed and the tissues representing the cambial zone (differentiating xylem) were scrapped with a razor blade and immediately frozen with liquid nitrogen. (A) Bark removal; (B) and (C) scrapping the tissue of cambial zone; (D) and (E) transversal section of eucalyptus wood before removal of the bark and after removal of the bark, respectively. Bars = 100 μm.
Supplementary Figure 2. Representative 2-DE gels of eucalyptus cambial region during summer (A) and winter (B). Arrows indicate spots collected for protein identification.
Supplementary Table 1. Primers used for RT-qPCR and the respective amplification size, efficiency (%) and P(H1) values.
Supplementary Table 2. Matrix containing the original metabolite data.
Andersen, C. L., Jensen, J. L., and Orntoft, T. F. (2004). Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 64, 5245–5250. doi: 10.1158/0008-5472.CAN-04-0496
Balducci, L., Deslauriers, A., Giovannelli, A., Beaulieu, M., Delzon, S., Rossi, S., et al. (2014). How do drought and warming influence survival and wood traits of Picea mariana saplings? J. Exp. Bot. 4, 2–13. doi: 10.1093/jxb/eru431
Barbaroux, C., and Breda, N. (2002). Contrasting distribution and seasonal dynamics of carbohydrate reserves in stem wood of adult-ring porous sessile oak and diffuse-porous beech trees. Tree Physiol. 22, 1201–1210. doi: 10.1093/treephys/22.17.1201
Begum, S., Nakaba, S., Yamagishi, Y., Oribe, Y., and Funada, R. (2013). Regulation of cambial activity in relation to environmental conditions: understanding the role of temperature in wood formation of trees. Physiol. Plant. 147, 46–54. doi: 10.1111/j.1399-3054.2012.01663.x
Bochkov, D. V., Sysolyatin, S. V., Kalashnikov, A. I., and Surmacheva, I. A. (2012). Shikimic acid: review of its analytical, isolation, and purification techniques from plant and microbial sources. J. Chem. Biol. 5, 5–17. doi: 10.1007/s12154-011-0064-8
Bonhomme, M., Peuch, M., Ameglio, T., Rageau, R., Guilliot, A., Decourteix, M., et al. (2009). Carbohydrate uptake from xylem vessels and its distribution among stem tissues and buds in walnut (Juglans regia L.). Tree Physiol. 30, 89–102. doi: 10.1093/treephys/tpp103
Bradford, M. M. (1976). A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye biding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3
Callado, C. H., Roig, F. A., Tomazello-Filho, M., and Barros, C. B. (2013). Cambial growth periodicity studies of south american woody species – a review. IAWA J. 34, 213–230. doi: 10.1163/22941932-00000019
Callado, C. H., Vasconsellos, T. J., Costa, M. S., Barros, C. F., Roig, F. A., and Tomazello-Filho, M. (2014). Studies on cambial activity: advances and challenges in the knowledge of growth dynamics of Brazilian woody species. An. Acad. Bras. Cien. 86, 278–223. doi: 10.1590/0001-3765201320130033
Candiano, G., Bruschi, M., Musante, L., Santucci, L., Ghiggeri, G. M., Carnemolla, B., et al. (2004). Blue silver: a very sensitive colloidal Coomassie G-250 staining for proteome analysis. Electrophoresis 25, 1327–1333. doi: 10.1002/elps.200305844
Canham, C. D., Kobe, R. K., Latty, E. F., and Chazdon, R. L. (1999). Interspecific and intraspecific variation in tree seedling survival: effects of allocation to roots versus carbohydrate reserves. Oecologia 121, 1–11. doi: 10.1007/s004420050900
Carvalho, M. C. C. G., Gomes, D. G. C., Carneiro, R. T., Moon, D. H., Salvatierra, G. R., Franceschin, L. M., et al. (2008). SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis. Tree Physiol. 28, 905 -918. doi: 10.1093/treephys/28.6.905
Celedon, P. A., De Andrade, A., Meireles, K. G., Carvalho, M. C. C. G., Caldas, D. G., Moon, D. H., et al. (2007). Proteomic analysis of the cambial region in juvenile Eucalyptus grandis at three ages. Proteomics 13, 2258–2274. doi: 10.1002/pmic.200600989
Cernusak, L., and Hutley, L. B. (2011). Stable isotopes reveal the contribution of corticular photosynthesis to growth in branches of Eucalyptus miniata. Plant Physiol. 155, 515–523. doi: 10.1104/pp.110.163337
Chaffey, N., Cholewa, E., Regan, S., and Sundberg, B. (2002). Secondary xylem development in Arabidopsis: a model for wood formation. Physiol. Plant. 114, 594–600. doi: 10.1034/j.1399-3054.2002.1140413.x
Chia, D. W., Yoder, T. J., Reiter, W. D., and Gibson, S. I. (2000). Fumaric acid: an overlooked form of fixed carbon in Arabidopsis and other plant species. Planta 211, 743–751. doi: 10.1007/s004250000345
Cooke, J. E., and Weih, M. (2005). Nitrogen storage and seasonal nitrogen cycling in Populus: bridging molecular physiology and ecophysiology. New Phytol. 167, 19–30. doi: 10.1111/j.1469-8137.2005.01451.x
Deslauriers, A., Beaulieu, M., Balducci, L., Giovannelli, A., Gagnon, M. J., and Rossi, S. (2014). Impact of warming and drought on carbon balance related to wood formation in black spruce. Ann. Bot. 114, 335–345. doi: 10.1093/aob/mcu111
Drew, M. C. (1997). Oxygen deficiency and root metabolism: injury and acclimation under hypoxia and anoxia. Annu. Rev. Plant Physiol. Plant Mol. Biol. 48, 223–250. doi: 10.1146/annurev.arplant.48.1.223
Durand, T. C., Sergeant, K., Renaut, J., Planchon, S., Hoffmann, L., Carpin, S., et al. (2011). Poplar under drought: comparison of leaf and cambial proteomic responses. J. Proteomics 74, 1396–1410. doi: 10.1016/j.jprot.2011.03.013
Eyles, A., Pinkard, E. A., O'Grady, A. P., Worledge, D., and Warren, C. R. (2009). Role of corticular photosynthesis following defoliation in Eucalyptus globulus. Plant Cell Environ. 32, 1004–1014. doi: 10.1111/j.1365-3040.2009.01984.x
Fernández, M., Troncoso, V., and Valenzuela, S. (2015). Transcriptome profile in response to frost tolerance in Eucalyptus globulus. Plant Mol. Biol. Rep. 33, 1472–1485. doi: 10.1007/s11105-014-0845-7
Fernie, A. L., Carrari, F., and Sweetlove, L. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr. Opin. Plant Biol. 7, 254–261. doi: 10.1016/j.pbi.2004.03.007
Galindo-González, L. M., El-Kayal, W., Ju, C. J. T., Allen, C. C. G., King-Jones, S., and Cooke, J. E. K. (2012). Integrated transcriptomic and proteomic profiling of white spruce stems during the transition from active growth to dormancy. Plant Cell Environ. 35, 682–701. doi: 10.1111/j.1365-3040.2011.02444.x
Gion, J. M., Lalanne, C., Le Provost, G., Ferry-Dumazet, H., Paiva, J., Chaumeil, P., et al. (2005). The proteome of maritime pine wood forming tissue. Proteomics 5, 3731–3751. doi: 10.1002/pmic.200401197
Goué, N., Lesage-Descauses, M. C., Mellerowicz, E. J., Magel, E., Label, P., and Sundberg, B. (2008). Microgenomic analysis reveals cell type-specific gene expression patterns between ray and fusiform initials within the cambial meristem of Populus. New Phytol. 180, 46–56. doi: 10.1111/j.1469-8137.2008.02556.x
Gricar, J., and Cufar, K. (2008). Seasonal dynamics of phloem and xylem formation in silver fir and Norway spruce as affected by drought. Russ. J. Plant Physiol. 55, 538–543. doi: 10.1134/S102144370804016X
Gullberg, J., Jonsson, P., Nordström, A., Sjöström, M., and Moritz, T. (2004). Design of experiments: an eficient strategy to identify factors influencing extraction and derivatization of Arabidopsis thaliana samples in metabolomic studies with gas chromatography/mass spectrometry. Anal. Biochem. 331, 283–295. doi: 10.1016/j.ab.2004.04.037
Hamanishi, E. T., Barchet, G. L. H., Dauwe, R., Mansfield, S. D., and Campbell, M. M. (2015). Poplar trees reconfigure the transcriptome and metabolome in response to drought in a genotype- and time-of-day-dependent manner. BMC Genomics 16:329. doi: 10.1186/s12864-015-1535-z
Hertzberg, M., Aspeborg, H., Schrader, J., Andersson, A., Erlandsson, R., Blomqvist, K., et al. (2001). A transcriptional roadmap to wood formation. Proc. Natl. Acad. Sci. U.S.A. 98, 14732–14737. doi: 10.1073/pnas.261293398
Hoffman, D. E., Jonsson, P., Bylesjö, M., Trygg, J., Antti, H., Eriksson, M. E., et al. (2010). Changes in diurnal patterns within the Populus transcriptome and metabolome in response to photoperiod variation. Plant Cell Environ. 23, 1298–1313. doi: 10.1111/j.1365-3040.2010.02148.x
Jiang, X., Han, G., Ye, M., and Zou, H. (2007). Optimization of filtering criterion for SEQUEST database searching to improve proteome coverage in shotgun proteomics. BMC Bioinformatics 8:323. doi: 10.1186/1471-2105-8-323
Jonsson, P., Johansson, A. I., Gullberg, J., Trygg, J., Grung, B., Marklund, S., et al. (2005). High-throughput data analysis for detecting and identifying differences between samples in GC/MS-based metabolomic analyses. Anal. Chem. 77, 5635–5642. doi: 10.1021/ac050601e
Juan, D., Hong-Li, X., De-Qiang, Z., Xin-Qiang, H., Min-Jie, W., Ying-Zhang, L., et al. (2006). Regeneration of the secondary vascular system in poplar as a novel system to investigate gene expression by a proteomic approach. Proteomics 6, 881–889. doi: 10.1002/pmic.200401348
Ko, J. H., Prassinos, C., Keathley, D., and Han, K. H. (2011). Novel aspects of transcriptional regulation in the winter survival and maintenance mechanism of poplar. Tree Physiol. 31, 208–225. doi: 10.1093/treephys/tpq109
Korn, M., Gartner, T., Erban, A., Kopka, J., Selbig, J., and Hincha, D. K. (2010). Predicting Arabidopsis freezing tolerance and heterosis in freezing tolerance from metabolite composition. Mol. Plant. 3, 224–235. doi: 10.1093/mp/ssp105
Larisch, C., Dittrich, M., Wildhagen, H., Lautner, S., Fromm, J., Polle, A., et al. (2012). Poplar wood rays are involved in seasonal remodeling of tree physiology. Plant Physiol. 160, 1515–1529. doi: 10.1104/pp.112.202291
Leonardi, G. A., Carlos, N. A., Mazzafera, P., and Balbuena, T. S. (2015). Eucalyptus urograndis stem proteome is responsive to short-term cold stress. Genet. Mol. Biol. 39, 191–198. doi: 10.1590/S1415-475738220140235
Li, X., Wu, H. X., and Southerton, S. G. (2010). Seasonal reorganization of the xylem transcriptome at different tree ages reveals novel insights into wood formation in Pinus radiata. New Phytol. 187, 764–776. doi: 10.1111/j.1469-8137.2010.03333.x
Liu, Y., Jiang, Y., Lan, J., Zou, Y., and Gao, J. (2014). Comparative transcriptomic analysis of the response to cold acclimation in Eucalyptus dunnii. PLoS One 9:e113090. doi: 10.1371/journal.pone.0113091
López-Bucio, J., Nieto-Jacobo, M. F., Ramírez-Rodríguez, V., and Herrera-Estrella, L. (2000). Organic acid metabolism in plants: from adaptive physiology to transgenic varieties for cultivation in extreme soils. Plant Sci. 160, 1–13. doi: 10.1016/S0168-9452(00)00347-2
Moller, I. M. (2001). Plant mitochondria and oxidative stress: Electron transport, NADPH turnover, and metabolism of reactive oxygen species. Annu. Rev. Plant Physiol. Plant Mol. Biol. 52, 561–591. doi: 10.1146/annurev.arplant.52.1.561
Moreau, C., Aksenov, N., Lorenzo, M. G., Segerman, B., Funk, C., Nilsson, P., et al. (2005). A genomic approach to investigate developmental cell death in woody tissues of Populus trees. Genome Biol. 6:R64. doi: 10.1186/gb-2005-6-4-r34
Moulin, M., Deleu, C., Larher, F. R., and Bouchereau, A. (2002). High-performance liquid chromatography determination of pipecolic acid after precolumn ninhydrin derivatization using domestic microwave. Anal. Biochem. 308, 320–327. doi: 10.1016/S0003-2697(02)00202-6
Nishizawa, A., Yabuta, Y., and Shigeoka, S. (2008). Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 147, 1251–1263. doi: 10.1104/pp.108.122465
Nogueira, G. S., Marshall, P. L., Leite, H. G., and Campos, J. C. C. (2015). Thinning intensity and pruning impacts on Eucalyptus plantations in Brazil. Int. J. For. Resear. 2015, 1–10. doi: 10.1155/2015/168390
Oliveira, L. A., Breton, M. C., Bastolla, F. M., Camargo, S. S., Margis, R., Frazzon, J., et al. (2012). Reference genes for the normalization of gene expression in Eucalyptus species. Plant Cell Physiol. 53, 405–422. doi: 10.1093/pcp/pcr187
Pfaffl, M. W., Horgan, G. W., and Dempfle, L. (2002). Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 30:e36. doi: 10.1093/nar/30.9.e36
Pfanz, H., Aschan, G., Langenfeld-Heyser, R., Wittmann, C., and Loose, M. (2002). Ecology and ecophysiology of tree stems: corticular and wood photosynthesis. Naturwissenschaften 89, 147–162. doi: 10.1007/s00114-002-0309-z
Plavcová, L., Hacke, U. G., Almeida-Rodriguez, A. M., Li, E., and Douglas, J. C. (2013). Gene expression patterns underlying changes in xylem structure and function in response to increased nitrogen availability in hybrid poplar. Plant Cell Environ. 36, 186–199. doi: 10.1111/j.1365-3040.2012.02566.x
Plomion, C., Lalanne, C., Claverol, S., Meddour, H., Kohler, A., and Bogeat-Triboulot, M. (2006). Mapping the proteome of poplar and application to the discovery of drought-stress responsive proteins. Proteomics 6, 6509–6527. doi: 10.1002/pmic.200600362
Ramakers, C., Ruijtera, J. M., Depreza, R. H. L., and Moorman, A. F. M. (2003). Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 339, 62–66. doi: 10.1016/S0304-3940(02)01423-4
Rocha, M., Licausi, F., Araújo, W. L., Nunes-Nesi, A., Sodek, L., Fernie, A. L., et al. (2010). Glycolysis and the tricarboxylic acid cycle are linked by alanine aminotransferase during hypoxia induced by waterlogging of Lotus japonicas. Plant Physiol. 152, 1501-1513. doi: 10.1104/pp.109.150045
Salazar, M. M., Nascimento, L. C., Camargo, E. L. O., Goncalves, D. C., Neto, J. L., Marques, W. L., et al. (2013). Xylem transcription profiles indicate potential metabolic responses for economically relevant characteristics of Eucalyptus species. BMC Genomics 14:201. doi: 10.1186/1471-2164-14-201
Salazar, M. M., Nascimento, L. C., Camargo, E. L. O., Vidal, R. O., Lepikson-Neto, J., Goncalves, D. C., et al. (2011). Comparative transcriptome analysis of tree Eucalyptus species using RNAseq technology: analysis of genes interfering in wood quality aspects. BMC Proc. 5(Suppl. 7):P175. doi: 10.1186/1753-6561-5-S7-P175
Saveyn, A., Steppe, K., Ubierna, N., and Dawson, T. E. (2010). Woody tissue photosynthesis and its contribution to trunk growth and bud development in young plants. Plant Cell Environ. 33, 1949–1958. doi: 10.1111/j.1365-3040.2010.02197.x
Schauer, N., Steinhauser, D., Strelkov, S., Schomburg, D., Allison, G., Moritz, T., et al. (2005). GC-MS libraries for the rapid identification of metabolites in complex biological samples. FEBS Lett. 579, 1332–1337. doi: 10.1016/j.febslet.2005.01.029
Schrader, J., Nilsson, J., Mellerowicz, E., Berglund, A., Nilsson, P., Hertzberg, M., et al. (2004). A high-resolution transcript profile across the wood-forming meristem of poplar identifies potential regulators of cambial stem cell identity. Plant Cell 16, 2278–2292. doi: 10.1105/tpc.104.024190
Shu, L., Qiaojun, L., Ma, C., Ding, W., Zhou, J., and Wu, J. (2011). Genetic, proteomic and metabolic analysis of the regulation of energy storage in rice seedlings in response to drought. Proteomics 11, 4122–4138. doi: 10.1002/pmic.201000485
Simard, S., Giovannelli, A., Treydte, K., Traversi, M. L., King, G. M., Frank, D., et al. (2013). Intra-annual dynamics of non-structural carbohydrates in the cambium of mature conifer trees reflects radial growth demands. Tree Physiol. 33, 913–923. doi: 10.1093/treephys/tpt075
Srivastava, V., Obudulu, O., Bygdell, J., Löfstedt, T., Rydén, P., Nilsson, R., et al. (2013). OnPLS integration of transcriptomic, proteomic and metabolomic data shows multi-level oxidative stress responses in the cambium of transgenic hipI- superoxide dismutase populus plants. BMC Genomics 14:893. doi: 10.1186/1471-2164-14-893
Taji, T., Ohsumi, C., Iuchi, S., Seki, M., Kasuga, M., Kobayashi, M., et al. (2002). Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 29, 417–426. doi: 10.1046/j.0960-7412.2001.01227.x
Tegeder, M., and Weber, A. P. M. (2007). “Metabolite transporters in the control of plant primary metabolism,” in Control of Primary Metabolism in Plants, eds W. C. Plaxton and M. T. McManus (Oxford, UK: Blackwel Publushing), 85–120.
Thumma, B. R., Sharma, N., and Southerton, S. G. (2012). Transcriptome sequencing of Eucalyptus camaldulensis seedlings subjected to water stress reveals functional single nucleotide polymorphisms and genes under selection. BMC Genomics 13:364. doi: 10.1186/1471-2164-13-364
Travert, S., Valerio, L., Fourasté, I., Boudet, A. M., and Teulieres, C. (1997). Enrichment in specific soluble sugars of two Eucalyptus cell-suspension cultures by various treatments enhances their frost tolerance via a non colligative mechanism. Plant Physiol. 114, 1333–1442.
Turhan, E., and Ergin, S. (2012). Soluble sugars and sucrose-metabolizing enzymes related to cold acclimation of sweet cherry cultivars grafted on different rootstocks. Sci. World J. 2012, 1–7. doi: 10.1100/2012/979682
Tuttle, J. R., Nah, G., Duke, M. V., Alexander, D. C., Guan, X., Song, Q., et al. (2015). Metabolomic and transcriptomic insights into how cotton fiber transitions to secondary wall synthesis, represses lignification, and prolongs elongation. BMC Genomics 16:477. doi: 10.1186/s12864-015-1708-9
van Dongen, J. T., Gupta, K. J., Ramírez-Aguilar, S. J., Araújo, W. L., Nunes-Nesi, A., and Fernie, A. R. (2011). Regulation of respiration in plants: a role for alternative metabolic pathways. J. Plant Physiol. 168, 1434–1443. doi: 10.1016/j.jplph.2010.11.004
Vizcaino, J. A., Cote, R. G., Csordas, A., Dianes, J. A., Fabregat, A., Foster, J. M., et al. (2013). The Proteomics Identifications (PRIDE) database and associated tools: status in 2013. Nucleic Acids Res. 41, D1063–D1069. doi: 10.1093/nar/gks1262
Wagner, F. H., Hérault, B., Bonal, D., Stahl, C., and Andreson, L. O. (2016). Climate seasonality limits carbon assimilation and storage in tropical forests. Biogeosci. Discuss. 2016, 1–50. doi: 10.5194/bg-2015-619
Wang, W., Vinocur, B., Shoseyov, O., and Altman, A. (2004). Role of plant heat- hock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 9, 244–252. doi: 10.1016/j.tplants.2004.03.006
Xia, J., Mandal, R., Sinelnikov, I., Broadhurst, D., and Wishart, D. S. (2012). MetaboAnalyst 2.0 - a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 40, W127–W133. doi: 10.1093/nar/gks374
Xu, P., Kong, Y., Song, D., Huang, C., and Li, X. (2014). Conservation and functional influence of alternative splicing in wood formation of Populus and Eucalyptus. BMC Genomics 15:780. doi: 10.1186/1755-8794-7-S2-S1
Yang, J., Kamdem, D. P., Keathley, D. E., and Han, K. H. (2004). Seasonal changes in gene expression at the sapwood–heartwood transition zone of black locust (Robinia pseudoacacia) revealed by cDNA microarray analysis. Tree Physiol. 24, 461–474. doi: 10.1093/treephys/24.4.461
Yu, Q., Plotto, A., Baldwin, E. A., Bai, J., Huang, M., Yu, Y., et al. (2015). Proteomic and metabolomic analyses provide insight into production of volatile and non-volatile flavor components in mandarin hybrid fruit. BMC Plant Biol. 15:76. doi: 10.1186/s12870-015-0466-9
Zabalza, A., Van Dongen, J. T., Froehlich, A., Oliver, S. N., Faix, B., Gupta, K. J., et al. (2009). Regulation of respiration and fermentation to control the plant internal oxygen concentration. Plant Physiol. 149, 1087–1098. doi: 10.1104/pp.108.129288
Zulak, K. G., Weljie, A. M., Vogel, H. J., and Facchini, P. (2008). Quantitative 1H-NMR metabolomics reveals extensive metabolic reprogramming of primary and secondary metabolism in elicitor-treated opium poppy cell cultures. BMC Plant Biol. 8:5. doi: 10.1186/1471-2229-8-5
Keywords: Eucalyptus grandis, primary metabolism, RT-qPCR, 2-DE Proteome, metabolome, ethanolic fermentation
Citation: Budzinski IGF, Moon DH, Lindén P, Moritz T and Labate CA (2016) Seasonal Variation of Carbon Metabolism in the Cambial Zone of Eucalyptus grandis. Front. Plant Sci. 7:932. doi: 10.3389/fpls.2016.00932
Received: 02 February 2016; Accepted: 11 June 2016;
Published: 28 June 2016.
Edited by:Jean Rivoal, Université de Montréal, Canada
Reviewed by:Atsushi Fukushima, RIKEN, Japan
Abir U. Igamberdiev, Memorial University of Newfoundland, Canada
Jorge Almiro P. Paiva, Institute of Plant Genetics Polish Academy of Sciences, Poland
Copyright © 2016 Budzinski, Moon, Lindén, Moritz and Labate. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Carlos A. Labate, email@example.com