Integrated analysis of gene expression from carbon metabolism, proteome and metabolome, reveals altered primary metabolism in Eucalyptus grandis bark, in response to seasonal variation
© The Author(s). 2016
Received: 1 December 2015
Accepted: 22 June 2016
Published: 1 July 2016
Seasonal variation is presumed to play an important role in the regulation of tree growth, especially for Eucalyptus grandis, a fast-growing tree. This variation may induce changes in the whole tree at transcriptional, protein and metabolite levels. Bark represents an important group of tissues that protect trees from desiccation and pathogen attack, and it has been identified as potential feedstock for lignocellulosic derived biofuels. Despite the growing interest, little is known about the molecular mechanisms that regulates bark metabolism, particularly in tropical countries.
In this study we report the changes observed in the primary metabolism of E. grandis bark during two contrasting seasons in Brazil, summer (wet) and winter (dry), through the combination of transcripts (RT-qPCR), proteome (2-DE gels) and metabolome (GC-MS) analysis, in an integrated manner. Twenty-four genes, involved in carbon metabolism, were analyzed in the two seasons. Eleven were up-regulated in summer, three were up-regulated in winter and ten did not show statistical differences in the expression pattern. The proteomic analysis using 2-DE gels showed 77 proteins expressing differences in abundance, with 38 spots up-regulated in summer and 37 in winter. Different metabolites significantly accumulated during winter.
This study revealed a metabolic reconfiguration in the primary metabolism of E. grandis bark, triggered by seasonal variation. Transcripts and protein data suggests that during winter carbohydrate formation seems to be favored by tree metabolism. Glucose, fructose and sucrose accumulated at significant levels during the winter.
Eucalyptus species are the most widely planted hardwood due to the quality of its wood. These fast-growing trees are cultivated under a range of different climates and can be destined to different industrial processes (e.g. pulp and paper, charcoal, fuel wood, and solid wood products). Most of the current Eucalyptus production in Brazil is cultivated in an area over 5.4 million hectares . Given its fast growing rates and coppicing ability, eucalyptus has also been identified as a potential feedstock for biofuels . Besides, bark is also a source of nutrients, carbon as well as being used to form a protective covering of soils in commercial plantations. Bark comprises all the tissues outside the vascular cambium and it includes primary and secondary phloem, cortex, first periderm, rhytidome and tissues formed by dilatation growth . These tissues also protect woody plant organs and healing tissues from dehydration, solar irradiation and pathogens . Bark formation is initiated by the process of cell division at the cambium, which produces xylem on the inner woody side and phloem, the primary bark tissue, on the exterior bark side. The phloem tissue contains phloem parenchyma, phloem fibers, companion cells and sieve cells . Although tree development is regulated by seasonal periods, little is known about the underlying molecular processes related to growth, especially in the bark. Soler et al.  analyzed the seasonal variation in mRNA abundance in cork tissue from Quercus suber. They found transcripts for structural genes involved in suberin production accumulating in late spring; this accumulation was significantly correlated with temperature and relative humidity. The increased expression of genes involved in stress was also strongly correlated to temperature. Using the proteomic approach, Pagter et al.  observed distinct seasonal protein patterns in bark of Hydrangea macrophylla and Hydrangea paniculata during cold acclimation and de-acclimation.
In the present study, we investigated the metabolic response of E. grandis bark in two contrasting seasons: summer/wet and winter/dry, using RT-qPCR, proteomic and metabolomic analyses, with emphasis on carbon metabolism. Despite the importance of bark metabolism to the whole tree and its potential for biofuel production , to our knowledge this is the first molecular study showing changes in eucalyptus bark metabolism, in response to seasonal variation.
Results and discussion
Seasonality influences mRNA expression in bark
Most of the genes involved in carbohydrate metabolism and glycolysis were up-regulated in summer (Fig. 1a), suggesting that during this season, primary metabolism is being required for the production of reducing power and ATP, necessary for the development and growth of the trees. During winter, the metabolic activity decreases and, as a consequence, trees reduce their growth rates, thus reducing their energy consumption. It is interesting to note that each SuSy analyzed showed a different expression pattern; SuSy1 was highly expressed in winter while SuSy3 were highly expressed in summer. Susy plays a crucial role in the sucrose metabolism. The enzyme catalyzes the reversible conversion of sucrose and UDP to UDP-glucose and fructose, but it also degrades sucrose . Thus, each alternative transcript could be required for a different molecular function. Moreover, it has been proposed that there are two forms of Susy in higher plants , one is a soluble enzyme found in the cytoplasm (S-Susy) and the other is a membrane associate enzyme (P-Susy). Both forms are probably regulated by the phosphorylation status of the enzyme . Further investigations are necessary to completely understand the mechanism that regulates Susy isoforms in eucalyptus bark.
The genes ADH3 and PDC involved in the fermentative metabolism showed statistically significant differences in expression (Fig. 1b). Curiously, each transcript was up-regulated in a different season. As we expected ADH3 was up-regulated in summer, during the growing season. We also expected that PDC was up-regulated in summer; however, it showed an opposite expression pattern. As a fast growing tree, eucalyptus requires high levels of ATP in summer to maintain this process. At the same time, the cambial and bark tissues are under hypoxic conditions due to anatomical barriers to gas exchange. Thus, respiration might shift from the aerobic to the ethanolic fermentation mode, as a means to maintain substrate-level for ATP production. This requires transcriptional activation of the essential genes of ethanolic fermentation, PDC and ADH . Our data suggest that ethanolic fermentation is probably required in summer and winter, with the participation of different PDC alternative transcripts. The regulatory role of PDC in ethanolic fermentation has not yet been fully explained . In Arabidopsis there are four genes encoding PDC and a microarray dataset related to low oxygen conditions, revealed that PDC1 and PDC2 were strongly up-regulated under low oxygen, whereas PDC3 and PDC4 mRNA levels were not induced by anoxia, suggesting that these two genes were unlikely to play a role during anoxic stress .
Yang et al.  found one ADH in the bark of Robinia pseudoacacia, however the authors did not discuss the transcript function. The study of the ethanolic fermentation process in trees started around the late 1980s [15–17] although until now little information is available and more research in this field is necessary.
IDH and SCL were the genes from TCA analyzed in our work, and both were up-regulated in summer (Fig. 1c). IDH is responsible to the oxidative decarboxylation of isocitrate to 2-oxoglutarate, the function of this enzyme has been associated with the maintenance of the 2-oxoglutarate level and the regulation on nitrogen assimilation . SCL catalyzes the reversible interconversion of succinyl-CoA to succinate, characterization of the regulatory properties of this enzyme suggests that it may represent an adaptive mechanism in the attempt to maintain the rate of respiration under suboptimal condition [19, 20]. Our results may suggest that in summer, the pyruvate produced via glycolysis is mainly sent to mitochondria, where it will be used by the TCA cycle, for energy production, instead of being used in the fermentative metabolism. It’s important to mention that as described above, most of the transcripts analyzed in carbohydrate metabolism and glycolysis were also up-regulated in summer, indicating that this pathway is working in the direction of pyruvate formation. Trees show higher metabolic activity in summer and during this period, active growth, glycolysis and the TCA cycle are fundamental to maintain high metabolic rate. Curiously, we observed that PEPC was up-regulated in summer. The export of TCA cycle intermediates requires the importation of substrates that can generate both acetyl-CoA and oxaloacetate (OAA) . If pyruvate is used as a unique substrate, export of TCA cycle intermediates, would reduce OAA regeneration and bring the TCA cycle to a halt. In such condition PEPC plays an important role in the anaplerotic fixation of CO2 and in the route to replenish TCA intermediates that are withdrawn from the pool [22, 23]. PEPC activity was reported as being more than ten times higher in stem of Fagus sylvatica than in leaves . This high PEPC activity in stem could be explained by the anaplerotic roles of PEPC, common in C3 plants or, as in C4 metabolism, PEPC could be supplying malate for NADP-ME, which supplies CO2 for Rubisco by malate descarboxylating .
In trees several physiological process such as growth, respiration, the lack of stomata and the low permeability of stem peridermal layers to gaseous diffusion, result in a high internal CO2 concentration ([CO2]) (range <1 to 26 %) and thus 500–800 times higher CO2 levels than in ordinary plant organs or ambient air [23, 24]. High CO2 concentrations effectively eliminate photorespiration, and enhance photosynthetic potential within branches and bark . In accordance with this, we found in eucalyptus bark the accumulation of mRNAs related to CO2 fixation (Fig. 1). Among the analyzed transcripts only RbcL and RPI were differentially expressed, they were up-regulated in summer. The NADP-ME transcript was not differentially expressed between summer and winter. RbcL was up-regulated in summer, the season in which trees have the highest respiration levels, leading to an increase in CO2 release in the interior of the stem. The RPI was also up-regulated in summer, it catalyzes the conversion of ribose 5-phosphate into ribulose 5-phosphate in the Calvin cycle and pentose phosphate pathways . Considering the high CO2 availability inside the bark, our results suggests that RbcL and RPI are acting in the refixation of metabolically produced CO2. The stems of woody plants possess greenish tissues, that contain chlorophylls (the chlorenchymes) and are localized below the outer peridermal or rhytidomal layers. These tissues are able to use the stem internal CO2 and the light penetrating the rhytidome to fixate carbon , explaining the presence of rubisco in the bark. Thus, the high CO2 concentration and the low oxygen availability can explain the fact that RbcL was differentially expressed and RbcS was not. Under such rich CO2 environment, inside the bark, maybe it is unnecessary for the plant to over produce RbcS to promote CO2 fixation. Higher plant Rubisco is composed of eight large subunits coded for by a single gene, the RbcL, and eight small subunits coded for by the nuclear RbcS multigene family . The RbcL contain the catalytic site of the enzyme and is responsible for the carboxylase and oxygenase reactions but the RbcS, whose precise role in structure and function of Rubisco remains poorly understood, contributes to the differences in kinetic properties among Rubisco enzyme . Transcripts of enzymes that participates in the photosystems I and II were also found in Robinia pseudoacacia bark .
Proteomic Analysis of E. grandis bark during summer and winter
Identification of differentially expressed proteins spots from 2-DE gels
N° of Peptides
Fold change (summer/winter)
22.214.171.124 C1 Metabolism
RuBisCO large subunit-binding
Ribulose bisphosphate carboxylase large chain
RuBisCO large subunit-binding
1.1.2 Energy Metabolism (Carbon)
1.1.3 Energy Transfer/Atp-Proton Motive Force
HSP20-like chaperones superfamily protein
ATPase, V1 complex, subunit B
ATP synthase alpha/beta family protein
1.2.1 Amino AcidMetabolism
1.2.3 Nucleotide/Nucleoside and Nucleotide Sugar Metabolism
1.2.8 Secondary Metabolism
2.1 Cell Processes
Vacuolar H + -ATPasecatalyticsubunit
2.2.2 Protection Responses/Detoxification
Glutathione S-transferase, C-terminal-like
126.96.36.199 Expansins, Xetand Extensin
Major Latex protein MLP-like
Major Latex protein (MLP-like)
5.2 Rna Related
Glycine-rich RNA-binding protein
5.3.3 Translation Related
Elongationfactor 1 beta
5.3.5 Protein Folding/Chaperoning
Heat shock cognate 70 kDa protein
Heatshock 70 kDaprotein
Heatshockcognate 70 kDa
Heat shock 70 kDa protein, mitochondrial
5.3.6 Protein Cleavage and Turnover
Proteasomesubunit beta type
26S proteasome non-ATPase regulatory particle
6.2 Putative Protein
NAD(P)-binding Rossmann-fold superfamily protein
Functional classification of proteins differentially expressed found in bark
Proteins identified related to the primary metabolism
Surprisingly we found only few proteins differentially expressed acting in primary metabolism. This result can be due to limitations of 2-DE gels technique or by the fact that only a reduced number of proteins related to primary metabolism changed in abundance between the two seasons. We found different proteins acting in carbohydrate metabolism, cell wall biosynthesis, glycolysis and TCA: UDP-glucose pyrophosphorylase (UGPase), UDP-glucose dehydrogenase (UGDH), phosphoglycerate kinase (three probable isoforms), triosephosphate isomerase, enolase (four probable isoforms), citrate synthase (two probable isoform) and malate dehydrogenase, all of them belong to different subcategories from the category “1. Metabolism and Energy”.
Four probable isoforms of UGPase were found, showing different expression patterns. Three of them were up-regulated in winter and the last one was up-regulated in summer, indicating temporal expression regulation. UGPase is a key enzyme in sucrose metabolism that catalyzes the reversible production of glucose-1-phosphate and UTP to UDP-glucose and pyrophosphate, depending on the metabolic status of the tissue. In photosynthetic tissues UGPase converts glucose-1-phosphate to UDP-glucose, which can be utilized for sucrose synthesis, or cell wall polysaccharides . In non-photosynthetic sink tissues, UGPase is linked to sucrose degradation pathways by converting UDP-glucose produced by sucrose synthase to glucose-1-phosphate . UGDH was up-regulated in winter and converts UDP-glucose to UDP-glucuronate, which is a precursor of hemicellulose and pectin. In woody tissues, the role of UGPase and UGP is poorly understood. Two UGDH genes mainly expressed in roots, stem and bark of 6-month-old E. grandis were cloned .
We found tree probable isoforms of phosphoglycerate kinase and four possible isoforms of triosephosphate isomerase, all of them up-regulated in winter. Two enolases were also up-regulated in winter and two were up-regulated in summer. Changes in abundance of proteins related to carbohydrate and energy metabolism during cold acclimation were observed in the Hydrangea paniculata bark . In the subcategory “Carbon Metabolism” we found tree RbcL isoforms, all of them up-regulated during winter. The RbcS was not found in our work, it is possible that RbcS was not detected as a differentially expressed spot. To our knowledge there are no data in literature related to the rubisco subunits expression pattern (transcripts and proteins) in trees bark. In the E. grandis cambial zone of trees with different ages, it has been demonstrated the presence of RbcL, by immunoblotting, the RbcS subunit was only detected in the leaf control . Proteins related to carbon fixation were found in the bark of Prunus persica  and Picea sitchensis . The protein profile of RbcL observed in our work was the opposite of the transcriptional pattern that we found, as the transcript was up-regulated in summer. The discrepancy between transcriptomic and proteomic data is widely discussed in the literature and probably indicates the occurrence of post-transcriptional and/or post-translational modifications [38–40].
Proteins identified related to other biological pathways
Proteomics provides the identification of a set of proteins expressed at a specific time, tissue or condition. Thus, we identified a diverse range of proteins involved in other biological process, besides primary metabolism. The proteins implicated in Protein folding/Chaperoning were the most abundant, among them 17 heat shock (HSP) were found and two were up-regulated in winter. Some HSPs are molecular chaperones that regulate the folding, localization, accumulation, and degradation of protein . Thus HSPs play a crucial role in protecting plants against multiple environmental stresses by re-establishing normal protein conformation and homeostasis. Other proteins related to stress were found, such as: one late embryogenic protein (LEA), seven ascorbate peroxidases (APX) and two 14-3-3 proteins. LEA protein was up-regulated in winter; this protein is found in plant seeds and also in vegetative tissues under stress conditions such as cold, drought, or high salinity . Among the seven APX identified, five were up-regulated in summer and two in winter. APX expression is induced in response to different forms of stress that results in the accumulation of reactive oxygen species (ROS) . 14-3-3 proteins were also found. These are phosphoserine-binding proteins that regulate the activities of a wide array of targets via direct protein–protein interactions. In plants 14-3-3 regulates the plasma membrane H+-ATPase and enzymes of carbon and nitrogen metabolism . Three translation-related proteins were identified, two initiation factors (up regulated in summer) and one elongation factor (up regulated in winter). Two caffeic acid o-methyltransferase (COMT) were up regulated in winter. COMT is one of the most important enzymes controlling lignin monomers production in plant cell wall synthesis.
Soluble sugars in Bark
Metabolic profiling during seasonal variation
Metabolites differentially abundant in bark
Integrated analyses of E. grandis bark in response to seasonal variations
We did not find proteins related to ethanolic fermentation, probably due to limitations in 2D-gels, on the other hand we found differential expression for ADH and PDC transcripts. ADH3 was up-regulated in summer and PDC was up-regulated in winter. Although we did not find any metabolite directly related to carbon fixation, transcripts and proteins acting in Calvin-Benson Cycle were identified, especially RbcL, strongly suggesting that carbon fixation is occurring in bark.
Studies about the metabolism of bark trees are scarce and the effect of seasonal variation is also poorly understood. Given the importance of Eucalyptus, especially in Brazil, for different industrial applications, it is important to comprehend how seasonal variations affect the whole tree and tissues. Our results strongly suggest a metabolic reconfiguration triggered by the shift between summer and winter periods, as we found significant differences in all levels investigated (transcripts, proteins and metabolites). During summer, when trees are fast growing, energy compounds are necessary to support glycolysis and mitochondrial electron transport chain. However, the high respiration rates associated with anatomical barriers generate a hypoxic environment inside the bark. Thus, ethanolic fermentation is an important pathway regenerating NAD+ to the maintenance of glycolysis and plant metabolism. It is known that bark is a storage tissue and that during winter tree growth diminished. In winter, soluble sugars accumulate probably because of the diminished water available in tropical regions, and not because of lower temperatures (cold acclimation) as observed in temperate trees. An interesting data we observed was the identification of RbcL transcripts and proteins in bark. This study provides important data to understand seasonal variation in Eucalyptus bark. Future studies are necessary to identify which isoforms are involved in ethanolic fermentation and carbohydrate metabolism and also to clarify the precise function of RbcL in bark.
Plant material and experimental conditions
Tissue samples were harvested from commercial clonal trees of six years-old Eucalyptus grandis, kindly provided by Suzano Papel e Celulose. The field-trial was situated in Itapetininga city, State of São Paulo, Brazil (23°35′20″ S, 48°03′11″ W) at an altitude of 656 m. To analyze the changes in transcripts, proteins and metabolites during summer/active growth compared to winter/diminished growth, bark samples were harvested in January/2009 and July/2009 for summer and winter, respectively. Samples were harvested in the morning, between 9 and 10:00 am. The average temperature and precipitation during January/summer were 23.7 °C and 213.2 mm, respectively. During the month of July/winter these parameters were 16.8 °C and 47.7 mm, respectively. The bark of each tree was removed at chest height (1.30 m, exposing an area of approximately 20 × 15 cm2). The cambial zone tissues in the inner surface of the bark were scraped with a razor blade and discarded, the bark samples were immediately frozen in liquid nitrogen. The field trial was a completely randomized design. Six bulked samples (10 trees each) were used as biological replicates. Three bulks represented summer and three represented winter trees. The plant material was used with the permission of Suzano Papel e Celulose S/A.
RNA extraction and mRNA isolation
Total RNA was extracted from bark samples using the protocol described by Zeng and Yang . Total RNA concentration was measured spectrophotometrically at 260/ 280 nm, in 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 using Dynabeads® mRNA purification kit (Invitrogen Dynal, Oslo, Norway), according to the manufacturer’s instructions.
Gene-specific primers were designed with Primer 3 software (Additional file 3: Table S1). Primer pairs were designed as follows: primer length between 18 pb - 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 SuperScriptTM One-Step RT-PCR Platinum® Taq (Invitrogen, Carlsbad, CA, USA) kit with RT/Platinum®Taq (Invitrogen, Carlsbad, CA, USA) and using primers specific for the genes of interest (Additional file 3: Table S1). 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 (10-1) 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). A negative control (no cDNA template) was included for every gene. Three biological replicates, each with three technical replicates, were analyzed. 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) . Reference genes (α-tubulin and citoplasmatic malate dehydrogenase (MDHc)) were identified using NormFinder . The software LinReg  was used to calculate the PCR efficiencies and the Cq values of each gene analyzed.
Total protein from bark was extracted by grinding the frozen tissue (4 g) and using the phenolic method according to Hurkman and Tanaka , with minor modifications described in Celedon et al. . The tissues were homogenized in 15 mL of extraction buffer (0.7 M sucrose, 0.5 M Tris-HCl, pH 7.5, 50 mM EDTA, 0.1 M KCl, 1 % w/v polyvinylpolypirrolidone (PVPP), 2 % v/v 2-mercaptoethanol, and 2 mM PMSF), by shaking for 30 min at 4 °C. An equal volume of Tris-HCl saturated phenol pH 8.5 was added to the protein suspension. After an additional 30 min of shaking at 4 °C, the phases were separated by centrifugation (10,000 × g for 30 min at 4 °C). The phenol phase was recovered and re-extracted with an equal volume of extraction buffer. Proteins were precipitated from the phenol phase by adding 5 vol. of 0.1 M ammonium acetate in methanol and incubated overnight at -20 °C. The samples were then centrifuged (10,000 × g, 30 min at 4 °C) and the resulting pellets were washed three times with 0.1 M ammonium acetate in methanol, followed by a wash with acetone. 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 . Protein samples (700 μg) were mixed with buffer (340 μL) containing 10 mM DTT, 4 % (w/v) CHAPS, 1 % IPG buffer (GE Healthcare, Chalfont St. Giles, UK). Bromophenol blue (1 % w/v) was used to rehydrate for 12 h (20 °C and 50 V) the strips of Immobiline IPG (18 cm -pH 4–7, linear gradient, GE Healthcare, Chalfont St. Giles, UK). Rehydrated strips were isoelectrofocused in an Ettan™ IPGphor II™ (GE Healthcare) for 1 h, starting at 100 V and then 500 V for 1 h, 1000 V for 1 h, 5000 V for 1 h and 8000 until reaching a total of 80,000 V-h. 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 a Protean II XI 2-D cell electrophoresis system (GE Helathcare), at 30 mA per gel until the dye reached the bottom of the gel. Three biological replicates were performed for each treatment. Proteins were detected using Coomassie Brilliant Blue G-250 . 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). Gels were imaged using an 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 Amersham Bioscience). 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: artefacts, 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 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 technical replicates, spots with coefficient of variation higher than 30 % were excluded. To identify spots significantly expressed the data collected from protein spot volumes were subjected to Student’s t-test (P ≤ 0,05) in Image Master v 7.0 software. In-gel digestion of proteins was performed as described in Celedon et al.  After, peptide mixtures were sequenced by online chromatography using a nano-Acquity UPLC (Waters®) sistem coupled to a Q-TOF Ultima API mass spectrometer (Waters, UK). Mass spectrometer parameters were adjusted according Celedon et al. . Ten microliters of sample were loaded onto a pre-column Symetry MCA C18 5 mm, 5630 mm (Waters) for sample pre-concentration and desalination, followed by peptide separation on an LC column Symmetry C185 mm, 32 × 150 mm (Waters). Peptides were eluted using a linear gradient (10–45 %) of solvent B (95 % (v/v) acetonitrile, 0.1 % (v/v) formic acid in water). The flow rate started with 5 mL/min for the first 15 min, then changing to 2 mL/min for the next 25 min, and back to 5 mL/min in the last 5 min. Solvent A consisted of 5 % v/v acetonitrile, and 0.1 % v/v formic acid in water. All analyses were performed using a positive ion mode at 3 kV needle voltage. The mass range was set from 300 to 2000 m/z, and the MS/MS spectra were acquired for the most intense peaks having at least 15 counts.
Protein identification/MSMS-data analysis
The LC-MS/MS were processed using ProteinLynx v 2.0 (Waters) and Mascot Daemon (Matrix Science, Boston, MA) software, and the sequences searched against an in-house Eucalyptus database from Phytozome v1.1 (www.phytozome.net/eucalyptus.php) and NCBI. Combined MS-MS/MS search criteria used were as follows: trypsin digestion; fixed modification set as carbamidomethylation of cystein; 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. .
Bark samples were ground into powder in liquid N2 and freeze-dried. Metabolites were extracted according to the method described by Hoffman et al. , 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)-disodiumketoglutarate, (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 rpm in an Eppendorf centrifuge (model 54178).
After, 100 mL 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. Trimethylsilylation 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. , 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 index to be calculated . One microliter of each derivatized sample was injected in splitless mode by a CTC Combi Pal Xt Duo (CTC Analytics AG, Switzerland) auto-sampler in 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 flow rate 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 50 − 800 m/z. 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 ChromaTOF 2.12 software (Leco Corporation), in which all manual integrations and metabolite identification were done. All data treatment procedures (base-line correction and chromatogram alignment) were performed using customs scripts  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 . The supervised classification method PLS-DA was carried out to discriminate between different groups (summer and winter). PLS-DA model fit was evaluated using the R2 and Q2 cross-validation performance measures , 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. To identify the metabolites that contributed to the separation between the two groups we used the VIP from the PLS-DA model. VIP is a weighted sum of squares of the PLS loadings which indicates the importance of the variable to the whole model. Differential metabolites were selected based on PLS-DA model using a combination of VIP value > 1.0 and p-value (P ≤ 0.05), by the univariate unpaired, two-tailed Student’s t-test.
Carbohydrate extraction and HPLC analysis
Soluble sugars (glucose, sucrose and fructose) were extracted from bark samples. Tissues were grinded 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 g, 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. 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 a Tukey’s test (P ≤ 0.05) was performed in SAS (SAS Institute, Cary, NC, USA).
2-DE, two-dimensional chain electrophoresis; ADH2, alcohol dehydrogenase 2; ADH3, alcohol dehydrogenase 3; ATP, Adenosine triphosphate; CA, carbonic anhydrase; Cq, threshold quantification cycle; ENO, enolase; FBAcl, fructose bisphosphate aldolase chloroplast; FBAcyt, fructose bisphosphate aldolase cytoplasmatic; GC-MS, gas chromatography mass spectrometry; GPI, glucose 6-phosphate isomerase; IDH, isocitrate dehydrogenase; NADP-ME, NADP Malic enzyme; PDC, pyruvate decarboxylase; PDH, pyruvate dehydrogenase; PEPC, phosphoenolpyruvate carboxylase; PFK, ATP-dependent phosphofructokinase; PFP, PPi-dependent phosphofructokinase; PGAM, phosphoglycerate mutase; PGK, phosphoglycerate kinase; PGM, phosphoglucomutase; PK, pyruvate kinase; PLS-DA, partial lest squares discriminant analysis; RbcL, rubisco large subunit; RbcS, Rubisco small subunit; RPI, ribose 5-phosphate isomerase; RT-qPCR, reverse transcription quantitative real time polymerase chain reaction; SCL, succinyl-CoA ligase; SuSy1, sucrose synthase 1; SuSy3, sucrose synthase 3; TCA, tricarboxylic acid cycle; UDP, uridine diphosphate; UTP, uridine triphosphate; VIP, variable importance in the projection
We would like to thank Dr. Esteban Gonzalez (FuturaGene Corporation) for providing the biological material and José Mateus Wisniewski Gonsalves (FuturaGene) and Dr. Luis Felipe Boaretto for the help with the samples harvesting. We would like to thank FAPESP for the fellowship to J.S.M. (Proc. no 2011/01236-2).
This work was supported by the FAPESP (São Paulo Research Foundation) (Grant 2009/00808-2 and 2008/50361-1).
Availability of data and materials
The data sets supporting the results of this article are included within the article and its additional files.
IGF conceived of the study, carried out the experiments, analyzed the data and wrote the manuscript. DHM assisted in the RT-qPCR analysis and manuscript revision. JSM participate in the RT-qPCR and proteomic analysis. JB assisted in the carbohydrate extraction and HPLC analysis. PL and TM assisted with metabolomics studies. CAL conceived of the study, participated in its design and coordination, and reviewed the manuscript. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Brazilian Association of Planted Forest Producers (ABRAF): Annual report. Brazil. 2013. http://www.ipef.br/estatisticas/relatorios/anuario-ABRAF13-EN.pdf. Accessed 29 Jan 2014.
- Albaugh JM, Dye PJ, King JS. Eucalyptus and water use in souths Africa. Int J For Res. 2013;2013:1–11.Google Scholar
- Lev-Yadun, Simcha: Bark - Encyclopedia of Life Sciences (ELS). Chichester, UK: John Wiley & Sons; 2011.Google Scholar
- Soler M, Serra O, Molinas M, Berthou EG, Caritat A, Figueras M. Seasonal variation in transcript abundance in cork tissue analyzed by real time RT-PCR. Tree Physiol. 2008;28:743–51.View ArticlePubMedGoogle Scholar
- Ogunwusi AA. Potentials of industrial utilization of bark. J Nat Sci Res. 2013;3:106–15.Google Scholar
- Pagter M, Sergeant K, Moller SM, Bertram HC. Changes in the proteome and water state in bark and xylem ofHydrangea paniculata during loss of freezing tolerance. Environ Exp Bot. 2013;106:99–111.View ArticleGoogle Scholar
- Lima MA, Lavorente GB, Silva HK, Bragatto J, Rezende CA, Bernardinelli OD, Azevedo ER, Gomez LD, McQueen-Mason SJ, Labate CA, Polikarpov I. Effects of pretreatment on morphology, chemical composition and enzymatic digestibility of eucalyptus bark: a potentially valuable source of fermentable sugars for biofuel production - part 1. Biotechnol Biof. 2013;6:75–92.View ArticleGoogle Scholar
- Kleczkowski LA. Glucose activation and metabolism through UDPglucose pyrophosphorylase in plants. Phytochemistry. 1994;37:1507–15.View ArticleGoogle Scholar
- Amor Y, Haigler CH, Johnson S, Wainscott M, Delmer DP. A membrane-associated form of sucrose synthase and its potential role in synthesis of cellulose and callose in plants. Proc Natl Acad Sci. 1995;92:9353–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Haigler CH, Ivanova-Datcheva M, Hogan PS, Salnikov VV, Hwang S, Martin LK, Delmer DP. Carbon partitioning to cellulose synthesis. Plant Mol Biol. 2001;47:29–51.View ArticlePubMedGoogle Scholar
- Giuntoli B, Lee SC, Licausi F, Kosmacz M, Oosumi T, van Dongen JT, Bailey-Serres J, Perata P. A trihelix DNA binding protein counterbalances hypoxia-responsive transcriptional activation in Arabidopsis. PLoS Biol. 2014;12:e1001950.View ArticlePubMedPubMed CentralGoogle Scholar
- Banti V, Giuntoli B, Gonzali S, Loreti E, Magneschi L, Novi G, Paparelli E, Parlanti S, Pucciariello C, Santaniello A, Perata P. Low oxygen response mechanisms in green organisms. Int J Mol Sci. 2013;14:4734–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Mithran M, Paparelli E, Novi G, Perata P, Loreti E. Analysis of the role of the pyruvate decarboxylase gene family in Arabidopsis thaliana under low oxygen conditions. Plant Biol. 2014;16:28–34.View ArticlePubMedGoogle Scholar
- Yang J, Kamdem DP, Keathley DE, Han KH. Seasonal changes in gene expression at the sapwood–heartwood transition zone of black locust (Robinia pseudoacacia) revealed by cDNA microarray analysis. Tree Physiol. 2004;24:461–74.View ArticlePubMedGoogle Scholar
- Kimmerer TW, Stringer MA. Alcohol dehydrogenase and ethanol in the stems of trees. Plant Physiol. 1988;87:693–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Harry DE, Kimmerer TW. Molecular genetics and physiology of alcohol dehydrogenase in woody plants. For Ecol Manage. 1991;43:251–72.View ArticleGoogle Scholar
- MacDonald RC, Kimmerer TW. Ethanol in the stems of trees. Physiol Plant. 1991;4:582–8.View ArticleGoogle Scholar
- Nunes-Nesi A, Fernie AR, Stitt M. Metabolic and signaling aspects underpinning the regulation of plant carbon nitrogen interactions. Mol Plant. 2010;3:973–96.View ArticlePubMedGoogle Scholar
- Studart-Guimarães C, Fait A, Nunes-Nesi A, Carrari F, Usadel B, Fernie AR. Reduced expression of Succinyl CoA ligase can be compensated for by an upregulation of the γ-amino-butyrate (GABA) shunt in illuminated tomato leaves. Plant Physiol. 2007;3:626–339.View ArticleGoogle Scholar
- Cavalcanti JHF, Esteves-Ferreira AA, Quinhones CGS, Pereira-Lima IA, Nunes-Nesi A, Fernie AR, Araújo WL. Evolution and functional implications of the tricarboxylic acid cycle as revealed by Phylogenetic Analysis. Genome Biol Evol. 2014;10:2830–48.View ArticleGoogle Scholar
- Tegeder M, Weber A. Metabolite transporters in the control of plant primary metabolism. In: Control of primary metabolism in plants. Volume 1st ed. Oxford: Willey; 2006. p. 85–120.View ArticleGoogle Scholar
- Plaxton WC, Podestá FE. The functional organization and control of plant respiration. Crit Rev Plant Sci. 2006;25:159–98.View ArticleGoogle Scholar
- Berveiller D, Damesin C. Carbon assimilation by tree stems: potential involvement of phosphoenolpyruvate carboxylase. Trees. 2008;22:149–57.View ArticleGoogle Scholar
- Aschan G, Pfanz H. Non-foliar photosynthesis – a strategy of additional carbon acquisition. Flora. 2003;198:81–97.View ArticleGoogle Scholar
- Wittmann C, Pfanz H, Loreto F, Centritto M, Pietrini F, Alessio G. Stem CO2 release under illumination: corticular photosynthesis, photorespiration or inhibition of mitochondrial respiration? Plant Cell Environ. 2006;29:1149–58.View ArticlePubMedGoogle Scholar
- Martin W, Henze K, Kellerman J, Flechner A, Schnarrenberger C. Microsequencing and cDNA cloning of the Calvin cycle/OPPP enzyme ribose-S-phosphate isomerase (EC 188.8.131.52) from spinach chloroplasts. Plant Mol Biol. 1996;30:795–805.View ArticlePubMedGoogle Scholar
- Pfanz H, Aschan G, Langenfeld-Heyser R, Wittmann C, Loose M. Ecology and ecophysiology of tree stems: corticular and wood photosynthesis. Naturwissenschaften. 2002;89:147–62.View ArticlePubMedGoogle Scholar
- Spreitzer RJ. Role of the small subunit in ribulose-1,5-bisphosphate carboxylase/oxygenase. Arch Biochem Biophis. 2003;2003:141–9.View ArticleGoogle Scholar
- Andersson I, Backlund A. Structure and Function of rubisco. Plant Physiol Bioch. 2008;2008:275–91.View ArticleGoogle Scholar
- Rison SC, Hodgman TC, Thornton JM. Comparison of functional annotation schemes for genomes. Funct Integr Genomics. 2000;1:56–69.View ArticlePubMedGoogle Scholar
- Carvalho MCCG, Gomes DGC, Carneiro RT, Moon DH, Salvatierra GR, Franceschin LM, De Andrade A, Celedon PAF, Oda S, Labate CA. SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis. Tree Physiol. 2008;28:905–18.View ArticleGoogle Scholar
- Kleczkowski LA, Geisler M, Ciereszko I, Johansson H. UDP-glucose pyrophosphorylase. An old protein with new tricks. Plant Physiol. 2004;134:912–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Winter H, Huber SC. Regulation of sucrose metabolism in higher plants: Localization and regulation of activity of key enzymes. Crit Rev Biochem Mol Biol. 2000;35:253–89.View ArticlePubMedGoogle Scholar
- Labate MTV, Bertolo ALF, Nascimento DD, Gutmanis G, Andrade A, Rodrigues MJC, Camargo ELO, Boaretto LF, Moon DH, Bragatto J. Labate CA. Cloning and endogenous expression of a Eucalyptus grandis UDP-glucose dehydrogenase cDNA. Genet Mol Biol. 2010;33:686–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Celedon PA, De Andrade A, Meireles KG, Carvalho MCCG, Caldas DG, Moon DH, Carneiro RT, Franceschini LM, Oda S, Labate CA. Proteomic analysis of the cambial region in juvenile Eucalyptus grandis at three ages. Proteomic. 2007;13:2258–74.View ArticleGoogle Scholar
- Renaut J, Hausman JF, Bassett C, Artlip T, Cauchie HM, Witters E, Wisniewski M. Quantitative proteomic analysis of short photoperiod and low-temperature responses in bark tissues of peach (Prunus persica L. Batsch). Tree Genet Genomes. 2008;4:589–600.View ArticleGoogle Scholar
- Lippert D, Chowrira S, Ralph SR, Zhuang J, Aeschliman D, Ritland C, Ritland K, Bohlmann, J. Conifer defense against insects: Proteome analysis of Sitka spruce (Picea sitchensis) bark induced by mechanical wounding or feeding by white pine weevils (Pissodesstrobi). Proteomics. 2007;7:248–70.View ArticlePubMedGoogle Scholar
- Yin L, Tao Y, Zhao Z, Shao J, Li X, Liu G, Liu S, Zhu L. Proteomic and transcriptomic analysis of rice mature seed-derived callus differentiation. Proteomics. 2007;7:755–68.View ArticlePubMedGoogle Scholar
- Fernie AR, Stitt M. On the discordance of metabolomics with proteomics and transcriptomics: coping with increasing complexity in logic, chemistry and network interactions. Plant Physiol. 2012;158:1139–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13:227–32.PubMedPubMed CentralGoogle Scholar
- Feder ME, Hofmann GE. Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999;61:243–82.View ArticlePubMedGoogle Scholar
- Hundertmark M, Hincha DK. LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008;9:118–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Mittler R, Lam E, Shulaev V, Cohen M. Signals controlling the expression of cytosolic ascorbate peroxidase during pathogen-induced programmed cell death in tobacco. Plant Mol Biol. 1999;39:1025–35.View ArticlePubMedGoogle Scholar
- Roberts MR. 14-3-3 Proteins find new partners in plant cell signaling. Trends Plant Sci. 2003;8:218–23.View ArticlePubMedGoogle Scholar
- Chantuma P, Lacointe A, Kasemsap P, Thanisawanyangkura S, Gohet E, Clement A, Guilliot A, Ameglio T, Thaler P. Carbohydrate storage in wood and bark of rubber trees submitted to different level of C demand induced by latex tapping. Tree Physiol. 2009;43:1–11.Google Scholar
- Ashworth EN, Stirm VE, Volenec JJ. Seasonal variations in soluble sugars and starch within woody stems of Corms sericea L. Tree Physiol. 1993;13:379–88.View ArticlePubMedGoogle Scholar
- Caldeira MC, Fernandéz V, Tomé J, Pereira JS. Positive effect of drought on longicorn borer larval survival and growth on eucalyptus. Ann For Sci. 2002;59:99–106.View ArticleGoogle Scholar
- Epron D, Nouvellon Y, Ryan MG. Introduction to the invited issue on carbon allocation of trees and forests. Tree Physiol. 2012;32:639–43.View ArticlePubMedGoogle Scholar
- Bonhomme M, Peuch M, Ameglio T, Rageau R, Guilliot A, Decourteix M, Alves G, Sakr S, Lacointe A. Carbohydrate uptake from xylem vessels and its distribution among stem tissues and buds in walnut (Juglans regia L.). Tree Physiol. 2009;30:89–102.View ArticlePubMedGoogle Scholar
- Travert S, Valerio L, Fourasté I, Boudet AM, Teulieres C. Enrichment in specific soluble sugars of two eucalyptus cell-suspension cultures by Various treatments enhances their frost tolerance via a noncolligative mechanism. Plant Physiol. 1997;114:1333–442.Google Scholar
- Herrmann KM, WEAVER LM. The shikimate pathway. Annu Rev Plant Physiol Plant Mol Biol. 1999;50:473–503.View ArticlePubMedGoogle Scholar
- Parson HT, Fry SC. Oxidation of dehydroascorbic acid and 2,3-diketogulonate under plant apoplastic conditions. Phytochemistry. 2011;75:41–9.View ArticleGoogle Scholar
- Fahnenstich H, Saigo M, Niessen M, Zanor MI, Andreo CS, Fernie AR, Drincovich MF, Flugge UI, Maurino VG. Alteration of organic acid metabolism in A. thaliana overexpressing the maize C4-NADPmalic enzyme causes accelerated senescence during extended darkness. Plant Physiol. 2007;145:640–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Hibberd JM, Quick P. Characteristics of C4 photosynthesis in stems and petioles of C3 flowering plants. Nature. 2002;415:451–4.View ArticlePubMedGoogle Scholar
- Pennycooke JC, Vepachedu R, Stushnoff C, Jones M. Expression of an alpha-galactosidase gene in petunia is upregulated during low-temperature deacclimation. J Am Soc Hortic Sci. 2004;129:491–6.Google Scholar
- Luo ZB, Calfapietra C, Liberloo M, Scarascia-Mugnozza G, Polle A. Carbon partitioning to mobile and structural fractions in poplar wood under elevated CO2 (EUROFACE) and N fertilization. GCB Bioenergy. 2006;12:272–83.Google Scholar
- Zeng Y, Yang T. Rna isolation from highly viscous samples rich in polyphenols and polysaccharides. Plant Mol Biol Rep. 2002;20:417–7.View ArticleGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30:e36.View ArticlePubMedPubMed CentralGoogle Scholar
- Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001;29:2003–7.View ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF. 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. 2004;64:5245–50.View ArticlePubMedGoogle Scholar
- Ramakers C, Ruijtera JM, Depreza RHL, Moorman AFM. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339:62–6.View ArticlePubMedGoogle Scholar
- Hurkman WJ, Tanaka CK. Solubilization of plant membrane proteins for analysis by two-dimensional gel electrophoresis. Plant Physiol. 1986;81:802–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bradford MM. A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye biding. Anal Biochem. 1976;72:248–54.View ArticlePubMedGoogle Scholar
- Candiano G, Bruschi M, Musante L, Santucci L, Ghiggeri GM, Carnemolla B, Orecchia P, Zardi L, Righetti, PG. Blue silver: a very sensitive colloidal Coomassie G-250 staining for proteome analysis. Electrophoresis. 2004;25:1327–33.View ArticlePubMedGoogle Scholar
- Perkins DN, Pappin DJ, Creasy DM, Cottrell JS. Probability-based protein identifications by searching sequence using mass spectrometry data. Electrophoresis. 1999;20:3551–67.View ArticlePubMedGoogle Scholar
- Hoffman DE, Jonsson P, Bylesjö M, Trygg J, Antti H, Eriksson ME, Moritz T. Changes in diurnal patterns within the Populus transcriptome and metabolome in response to photoperiod variation. Plant Cell Environ. 2010;23:1298–313.Google Scholar
- Gullberg J, Jonsson P, Nordström A, Sjöström M, Moritz T. Design of experiments: an eficient strategy to identify factors inXuencing extraction and derivatization of Arabidopsis thaliana samples in metabolomic studies with gas chromatography/mass spectrometry. Anal Biochem. 2004;331:283–95.View ArticlePubMedGoogle Scholar
- Jonsson P, Johansson AI, Gullberg J, Trygg J, Grung B, Marklund S, Sjostrom M, Antti H, Moritz T. High-throughput data analysis for detecting and identifying differences between samples in GC/MS-based metabolomic analyses. Anal Chem. 2005;77:5635–42.View ArticlePubMedGoogle Scholar
- Xia J, Wishart DS. Web-based inference of biological patterns, functions and pathways from metabolomic data using MetaboAnalyst. Nat Protoc. 2011;6:743–60.View ArticlePubMedGoogle Scholar