Email updates

Keep up to date with the latest news and content from BMC Plant Biology and BioMed Central.

Open Access Highly Accessed Research article

Pathway analysis of the transcriptome and metabolome of salt sensitive and tolerant poplar species reveals evolutionary adaption of stress tolerance mechanisms

Dennis Janz1, Katja Behnke2, Jörg-Peter Schnitzler2, Basem Kanawati3, Philippe Schmitt-Kopplin3 and Andrea Polle1*

Author Affiliations

1 Forstbotanik und Baumphysiologie, Büsgen-Institut, Georg-August-Universität Göttingen, Büsgenweg 2, 37077 Göttingen, Germany

2 Institut für Meteorologie und Klimaforschung, Bereich Atmosphärische Umweltforschung, Karlsruher Institut für Technologie (KIT), Kreuzeckbahnstraße 19 82467 Garmisch-Partenkirchen, Germany

3 Institut für Ökologische Chemie, Deutsches Forschungszentrum für Gesundheit und Umwelt, Helmholtz Zentrum München, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany

For all author emails, please log on.

BMC Plant Biology 2010, 10:150  doi:10.1186/1471-2229-10-150

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/10/150


Received:9 December 2009
Accepted:17 July 2010
Published:17 July 2010

© 2010 Janz et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Populus euphratica is a salt tolerant and Populus × canescens a salt sensitive poplar species. Because of low transcriptional responsiveness of P. euphratica to salinity we hypothesized that this species exhibits an innate activation of stress protective genes compared with salt sensitive poplars. To test this hypothesis, the transcriptome and metabolome of mature unstressed leaves of P. euphratica and P. × canescens were compared by whole genome microarray analyses and FT-ICR-MS metabolite profiling.

Results

Direct cross-species comparison of the transcriptomes of the two poplar species from phylogenetically different sections required filtering of the data set. Genes assigned to the GO slim categories 'mitochondria', 'cell wall', 'transport', 'energy metabolism' and 'secondary metabolism' were significantly enriched, whereas genes in the categories 'nucleus', 'RNA or DNA binding', 'kinase activity' and 'transcription factor activity' were significantly depleted in P. euphratica compared with P. × canescens. Evidence for a general activation of stress relevant genes in P. euphratica was not detected. Pathway analyses of metabolome and transcriptome data indicated stronger accumulation of primary sugars, activation of pathways for sugar alcohol production, and faster consumption of secondary metabolites in P. euphratica compared to P. × canescens. Physiological measurements showing higher respiration, higher tannin and soluble phenolic contents as well as enrichment of glucose and fructose in P. euphratica compared to P. × canescens corroborated the results of pathway analyses.

Conclusion

P. euphratica does not rely on general over-expression of stress pathways to tolerate salt stress. Instead, it exhibits permanent activation of control mechanisms for osmotic adjustment (sugar and sugar alcohols), ion compartmentalization (sodium, potassium and other metabolite transporters) and detoxification of reactive oxygen species (phenolic compounds). The evolutionary adaptation of P. euphratica to saline environments is apparently linked with higher energy requirement of cellular metabolism and a loss of transcriptional regulation.

Background

Salinization due to sea-level rise, decreased ground water recharge, inappropriate irrigation regimes, and increases in evapotranspiration is a major threat to sustainable land-use [1]. To prevent soil degradation and to maintain productivity in areas affected by salinity, cultivation of salt tolerant plants is required. This applies not only to agricultural crops, but also to woody plants that are gaining importance as renewable resources.

Over the past years, poplar (Populus sp.), a fast growing species with high biomass production, has been established as a model organism for tree research [2-4]. Among the species of this genus Populus euphratica Olivier occurs naturally in salt afflicted areas [5] and exhibits high salt tolerance [6]. Under saline conditions, P. euphratica is able to maintain higher growth rates and higher photosynthetic rates than other poplar species [7]. To gain insight into the molecular basis of its ability to resist salt stress, transcriptional changes were investigated during short-term salt exposure under controlled conditions as well as in long-term salt acclimated mature trees grown in natural environments [8,9]. Previous microarray analyses based on a stress-responsive EST collection of approximately 6340 unique P. euphratica genes revealed surprisingly few, i.e., only 22 genes whose transcript abundances were affected by salinity under field conditions [9]. We, therefore, hypothesize that P. euphratica developed innate mechanisms to tolerate salt stress that may require no immediate gene regulation.

Salt cress (Thellungiella halophila), a close relative of the herbaceous model plant Arabidopsis thaliana, occurs in harsh environments. When salt cress was exposed to excess salinity only relatively few stress-responsive genes were detected compared with Arabidopsis that exhibited a global defense strategy [10]. Transcriptome and metabolome analyses suggested stress-anticipatory preparedness in salt cress [10,11]. High levels of stress tolerance are the result of evolutionary adaptation. In addition to other possibilities, plants thriving in extreme environments may achieve protection by constitutive activation of stress-related gene networks. As a consequence these regulons can not be detected by classical microarray approaches analyzing transcriptional changes in response to stress.

To elucidate differences in constitutive gene regulation of salt tolerant and sensitive species, direct cross-species comparisons are required. A difficulty, usually limiting straightforward whole genome comparisons between different species, is that for non-model species microarrays are not available and hybridization efficiencies of different species are affected by sequence dissimilarities of bait set (= probe set of the model species on the microarray) and prey set (cDNAs of transcripts of either model or non-model species). In spite of the above problems, cross-species hybridizations have recently been successfully applied to a number of species when their evolutionary distances were close [12]. The problem has also been partially circumvented by disregarding genes with lower transcript abundances because 'apparent decreases' might have been caused by lower hybridization efficiencies, whereas stronger responses between non-model prey and model bait set are indicative for higher transcript abundance [13].

The major goal of our study was to investigate whether P. euphratica possesses a preventive stress defense strategy involving constitutive activation of protective pathways. To this end we compared available sequence data of P. euphratica genes with those of Populus × canescens (syn. P. tremula × P. alba [Aiton] Sm.), a salt sensitive species [14-16]. Since this analysis revealed high nucleotide identity we used Affymetrix GeneChip Poplar Genome Array to examine transcriptional differences of P. euphratica and P. × canescens. The Poplar Genome Array contains 56,055 probe sets (= baits) that can be interrogated, each composed of 11 individual probes. The array covers the whole genome of P. trichocarpa (31,999 probe sets based on gene predictions and EST evidence) and contains additionally 23,657 probes sets of 13 further poplar species (based on mRNA or EST sequence data) including P. euphratica as well as P. × canescens [17]. Here, we provide evidence that evolutionary distance of prey and bait sets affect the results but that a condensed data set can be used after application of appropriate quality filters. Key pathways identified by this approach to differ between the two poplar species were validated by targeted biochemical and physiological analyses as well as by non-targeted metabolite profiling applying Fourier Transform-Ion Cyclotron Mass Spectrometry (FT-ICR MS) [18-20]. Metabolome and transcriptome data were combined by pathway analysis tools [21,22]. We provide evidence that energy metabolism, ion compartmentalization, and phenylpropanoid biosynthetic pathways constitute major differences between P. euphratica and P. × canescens.

Results

Sequence identity of structural genes of P. trichocarpa, P. euphratica and P. × canescens is high

A prerequisite for a transcriptional comparison of P. euphratica and P. × canescens is a high sequence identity of their mRNAs with the probe set oligonucleotides on the microarray, which are mainly based on P. trichocarpa genome data [3]. Therefore, we compiled available sequence information for cDNAs of P. × canescens (Pc), P. euphratica (Pe), and P. trichocarpa (Pt). We found full-length and partial sequences of 20 genes for all three poplar species, which were used to assess their relationship by analysis of nucleotide identity (Table 1). Identity was usually high with mean values between 95.6 and 97.2%, and no significant differences between the pair wise comparisons of Pc/Pt and Pe/Pt, between Pe/Pc and Pc/Pt or between Pe/Pc and Pe/Pt were found at P ≤ 0.05 (Bonferoni test). However, there are individual examples for low sequence identity, e.g., NdID (NAD+ dependent isocitrate dehydrogenase in Table 1).

Table 1. Comparison of sequence data of 20 genes from P. × canescens (Pc), P. euphratica (Pe) and P. trichocarpa (Pt)

Detection and evaluation of transcripts differentially expressed between P. euphratica and P. × canescens

The transcriptomes of mature, fully light-exposed leaves of soil grown P. euphratica and P. × canescens plants were compared by microarray analysis on GeneChip Poplar Genome arrays. In the first step, differentially expressed transcripts were identified by established methods employing significance analysis of microarrays (SAM) and removing transcripts with low signal intensities (Table 2, see Materials and Methods). As cross-species comparisons may be affected by sequence divergence [23,24], we also removed genes whose standard deviation (SD) of mean signal intensities was in the upper 5% quantile. This procedure was based on the following reasoning: since each mean signal is the result of signal intensities of 11 oligonucleotides from different regions of a gene (= 1 probe set), it is likely that sequence divergence will lead to irregular binding of cDNAs within the probe set and thus cause higher SD. High sequence divergence may affect binding efficiencies and thus lead to erroneous results regarding transcript abundance. The validity of this filtering step was corroborated by qRT-PCRs conducted for genes with high and low SDs. Indeed, only those with SDs within the 90 to 95% range showed a significant correlation with the probe sets on the array (R = 0.909, P = 0.032), whereas those within the upper 5% quantile showed no correlation of qRT-PCR and array signal intensities (R = -0.376, P = 0.536, Figure 1).

Table 2. Overview of the different filtering steps applied to the list of significant transcripts and the gene universe

thumbnailFigure 1. Correlation of qRT-PCR expression data and microarray signal ratios. Data show means (SE) for qRT-PCR expression ratios and for signal ratios of microarray data for P. euphratica/P. × canescens. Transcripts with probe set SDs in the upper 5% quantile are indicated by open symbols, those with lower probe set SDs by black symbols. Star: outlier. The outlier has an SD in the 9% range, which is close to the chosen threshold of 5%. The following genes were included (putative function, Affymetrix probe set ID, JGI gene model for analyzed genes): 1, Gibberellin regulated protein, Ptp.6252.1.S1_a_at, estExt_Genewise1_v1.C_LG_V1745; 2, MADS-Box protein, Ptp.5993.1.S1_a_at, eugene3.00150771; 3, Mitochondrial carrier protein, Ptp.5103.1.S1_at, grail3.0008039502; 4, Lil3 protein, Ptp.4571.1.S1_at, eugene3.01180096; 5, Aquaporin, Ptp.5700.1.S1_s_at, eugene3.00280238; 6, GTP-binding protein, PtpAffx.25286.1.S1_at, estExt_fgenesh4_pg.C_LG_I1368; 7, Nitrogen fixation protein, PtpAffx.1459.1.A1_s_at, estExt_fgenesh4_pm.C_LG_XII0286; 8, Ubiquitin-like protein, PtpAffx.157059.1.S1_s_at, estExt_fgenesh4_pg.C_LG_XIV1291; 9, 1-Aminocyclopropane-1-carboxylate oxidase, Ptp.5158.1.S1_at, estExt_Genewise1_v1.C_1660131; 10, Glycine dehydrogenase, PtpAffx.19705.1.A1_at, estExt_fgenesh4_pm.C_LG_VI0678.

Despite these efforts, three different groups of transcriptional responses of the Pe/Pc ratio were detected depending on the species used to produce the array probe set (Figure 2). The first and largest group was based on sequences from poplars of the Tacamahaca and Aigeiros sections (P. trichocarpa, P. trichocarpa × deltoides, P. trichocarpa × nigra, P. nigra, P. deltoides and P. euramericana). The transcript ratios of Pe/Pc in this group showed a homogenous distribution of increased and decreased transcript ratios (Figure 2). This result is expected if the binding efficiencies of transcripts of the test species P. euphratica and P. × canescens with genes used to construct the probe sets on the array are similar. Only genes of this group were used for GO-term enrichment analyses.

thumbnailFigure 2. Transcript ratios of P. euphratica and P. × canescens sorted according to the original poplar species used to construct the probe set on the microarray. All transcripts with significantly different signal intensities that were detected as present by the MAS 5.0 algorithm and had probe set SDs in the lower 95% quantile were used (f2, Table 2). Signal log ratios were calculated as log-signal(P. euphratica) - log-signal(P. × canescens), i.e. a positive value denotes a higher apparent expression in P. euphratica, a negative value correspondingly a higher apparent expression in P. × canescens. Sample number of significant transcripts (n) is displayed for each species. Boxes represent the interquartile length (IQL) and median, whiskers extend to the most extreme data point with a maximum length of 1.5 times IQL; outliers are shown as circles;

The second group consisted of probe sets based on species of the Leuce section of which P. × canescens is a member (P. tremula × tremuloides, P. tremula, P. × canescens [= P. tremula × P. alba], P. alba × tremula var. glandulosa, P. tremuloides and P. tomentiglandulosa). In this group, the distribution of Pe/Pc transcript ratios showed a clear shift towards lower values. Since it is unlikely that the Leuce section contains an overabundance of down-regulated genes this result is probably influenced by sequence divergence resulting in relatively stronger binding of P. × canescens-derived cDNAs to the probe sets than those of P. euphratica, consequently leading to an apparent decrease of the Pe/Pc ratio. We have therefore excluded the genes with decreased Pe/Pc ratio in the Leuce group from further analyses. Vice versa, in the last group, the Turanga section which is formed by P. euphratica, the Pe/Pc ratio was increased. Along the same line of argument, the genes with increased Pe/Pc ratios were excluded (Table 2, filter f4*).

Since the goals of this study were functional categorization and pathway analyses, all transcripts without annotation were removed from the list of significant genes as proposed by Falcon and Gentleman [25]. The final list of differentially expressed transcripts contained 2503 genes, of which 1645 had higher transcript abundances in P. euphratica and 858 in P. × canescens than the respective other species (Table 2, Additional data Table S1).

The transcriptomes of P. euphratica and P. × canescens exhibit differences in transport, energy metabolism and biosynthesis of secondary metabolites

To characterize differences between the transcriptomes of the two poplar species, we used the GO slim terms of the complex Gene Ontology annotations [26]. A 'gene universe' of 14,254 eligible genes was compiled after appropriate filtering from all 61,413 probe sets on the microarray (Table 2). Significant over- and underrepresentation of GO slim terms for P. euphratica and P. × canescens were determined by a hypergeometric test. In P. euphratica GO slim terms associated with transcription and regulation like 'nucleus', 'DNA or RNA binding', 'transcription factor activity', 'kinase activity' and 'transcription' were significantly underrepresented (Figure 3). This suggests that mechanisms to react to alterations of external conditions and to adapt the transcriptome accordingly were less active in P. euphratica than in P. × canescens. Overrepresented GO slim terms in the group of genes with higher expression in P. euphratica were 'other membranes' (i.e. other than plasma membrane), 'mitochondria', 'cell wall', 'hydrolase activity' and 'transporter activity' (Figure 3). In P. × canescens, the only significantly overrepresented GO slim term was 'plasma membrane'. Neither of the species showed an over- or underrepresentation of the GO slim term 'response to stress' contrary to our expectation (Figure 3).

thumbnailFigure 3. Over- or underrepresented GO slim categories in P. euphratica (red) and P. × canescens (blue). Bars show frequencies of GO slim annotations of genes with significantly higher relative transcript abundance in P euphratica (red) or P. × canescens (blue) compared with the 'gene universe' as defined in table 2 (white). All GO slim categories significantly over- or underrepresented as calculated by a hypergeometric test after Benjamini-Hochberg correction are shown except 'unknown biological processes' (underrepresented in P. × canescens), 'other enzyme activity' (overrepresented in P. euphratica) and 'unknown cellular components' (underrepresented in P. euphratica). Data set after filter 5 (table 2) was used. Significant over- or underrepresentation of categories are indicated by * for P ≤ 0.05, ** for P ≤ 0.01 and for P ≤ 0.001.

For more detailed analysis, the best match of the poplar gene models with Arabidopsis thaliana was searched, AGI numbers were assigned to the poplar genes, and Mapman was used for a functional classification [27]. In the category 'transport' a number of genes for cation transport with higher transcript abundance in P. euphratica than in P. × canescens were especially conspicuous with respect to the increased salt tolerance of P. euphratica (Additional file 1, Table S2). In the category 'cell wall' numerous glycosylases and hydrolase such as putative xyloglucan endotransglycosylases, glycosyl hydrolase family proteins in addition to pectin esterases and fasciclin-like proteins had higher transcript abundance in P. euphratica than in P. × canescens (Additional file 1, Table S2). Analysis of the category 'mitochondria' revealed that most components of the mitochondrial electron transport chain were up-regulated in P. euphratica compared with P. × canescens (Figure 4). Furthermore, out of twelve significantly differentially expressed genes encoding mitochondrial metabolite transporters, ten showed higher transcript abundance in P. euphratica than in P. × canescens (Figure 4).

thumbnailFigure 4. Mitochondrial electron transport chain and metabolite transporters. Mapman analysis of the category mitochondria with relative increased transcript abundance in P. euphratica (red) and in P. × canescens (blue).

Additional file 1. Table S1: List of genes with significantly different signal intensities in P. euphratica compared with P. × canescens. The table contains fold-change of P. euphratica/P. × canescens, filter settings according to Table 2 in the manuscript and genes present after application of filters. Table S2: Mapman analysis of genes with significant differences in transcript levels P. euphratica compared with P. × canescens. The gene list with AGI annotations obtained after filer 5 was used and genes were to the 35 categories present under "Mapman overview". Table S3: List of metabolites with significant differences between P. euphratica compared with P. × canescens. Positive peak ratios indicate significant enrichment in P. euphratica and negative ratios indicate significant enrichment in P. × canescens, respectively. Table S4: Primers used to determine gene expression by qRT-PCR. Gene names, accession numbers and primer sequences are shown.

Format: XLS Size: 3.3MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Among all genes with significant changes that were uploaded into Mapman (n = 2216), the category 'secondary metabolism' was significantly enriched (P = 0.02 after Benjamini Hochberg correction). This result was corroborated using the original GO annotation of the poplar gene models. Both GO terms 'phenylpropanoid' and 'flavonoid biosynthetic process' were significantly enriched in P. euphratica compared with P. × canescens with P = 0.042 and P = 0.046, respectively, after Benjamini Hochberg correction.

Integration of metabolomic and transcriptomic data on pathway maps and validation by physiological and biochemical analyses

The metabolomes of leaves of the two poplar species were analyzed by Fourier transform-ion cyclotron resonance mass spectrometry (FT-ICR-MS). Identification of ions was performed using 'Mass Translator into Pathways' [MassTrix; 22]. 789 ions were found with had significantly higher peak intensities in P. euphratica than in P. × canescens. To 265 of these ions, a bulk chemical formula could be assigned, of which 75 were annotated in poplar. Where metabolites share an identical bulk formula, multiple alternative annotations may be possible. Therefore, the 75 ions represent 39 different bulk chemical formulas for metabolites with higher abundance in P. euphratica than in P. × canescens (Additional file 1, Table S3). In P. × canescens, 983 ions with significantly higher peak intensities than in P. euphratica were detected. Assignment of a bulk chemical formula was possible for 313 ions, of which 96 with 51 different bulk formulas were annotated in poplar (Additional file 1, Table S3).

Identified and annotated metabolites were automatically mapped to KEGG pathways using MassTrix, simultaneously implementing the transcriptomic data. For this, Enzyme Commission (EC) numbers were obtained for significant genes from the JGI P. trichocarpa project [28]. Of the differentially expressed genes, 289 genes with higher expression in P. euphratica than in P. × canescens and 132 genes with higher expression in P. × canescens than in P. euphratica could be annotated by EC numbers and were added to the MassTrix pathway analysis query.

Based on the results of the GO term enrichment analysis that indicated significant effects in secondary metabolism and mitochondrial energy metabolism, five pathways were chosen from the 124 KEGG pathways available for poplar: 'phenylpropanoid biosynthesis', 'flavonoid biosynthesis', 'citrate cycle' with the connected pathways 'glycolysis/gluconeogenesis' and 'starch and sucrose metabolism'.

The pathways of phenylpropanoid and flavonoid biosynthesis showed higher transcript abundances of enzymes from P. euphratica than for those from P. × canescens, whereas the intermediate metabolites, with the exception of sinapyl alcohol, were usually higher in P. × canescens than in P. euphratica (Figure 5). Since this result may indicate higher flux through the phenolic pathway in P. euphratica than in P. × canescens, which would be expected to result in increased concentrations of end products, we determined soluble and non-soluble phenolics and tannins. In fact, the concentrations of the final products, especially condensed non-soluble tannins, were significantly higher in P. euphratica than in P. × canescens (Figure 6a-d).

thumbnailFigure 5. Analysis of pathways related to phenolics metabolism. Map displays selected steps from KEGG pathways dpop00940 'Phenylpropanoid biosynthesis' and dpop00941 'Flavonoid biosynthesis'. Colors indicate significant expression, respective metabolite content ratios between P. euphratica and P. × canescens, red indicates higher relative levels in P. euphratica, blue in P. ×canescens. Enzymes are given as EC numbers: 1.1.1.195, cinnamyl-alcohol dehydrogenase; 1.1.1.219, dihydroflavonol 4-reductase; 1.11.1.7, peroxidase; 1.14.11.9, flavanone 3-hydroxylase; 1.2.1.44, cinnamoyl-CoA reductase; 3.2.1.21, beta-glucosidase; 5.5.1.6, chalcone isomerase; 6.2.1.12, 4-coumarate:CoA ligase.

thumbnailFigure 6. Phenolic compounds in leaves of P. euphratica (red) and P. canescens (blue). Soluble (a) and non-soluble phenolic compounds (b) were determined as catechin equivalents; soluble (c) and non-soluble condensed tannins (d) were determined as extinction units per gram g of fresh weight (n = 5, mean ± SD). Data in (b) are expressed mg per g of cell wall material (CW). Significant differences between the two poplar species are indicated by * for P ≤ 0.05, ** for P ≤ 0.01 and ** for P ≤ 0.001.

Since our transcriptome analysis pointed to higher activity of the mitochondrial electron transport chain in P. euphratica than in P. × canescens, we mapped enzymes and metabolites of the glycolysis and the tricarboxylic acid cycle, necessary to supply the respiratory chain with reductant and localized in the cytosol and the mitochondrial matrix, respectively (Figure 7). While many of the enzymes involved in these pathways were up-regulated at the transcriptional level in P. euphratica, only few differences were found regarding the participating metabolites. Among these metabolites most were increased in P. × canescens compared with P. euphratica (Figure 7). Interestingly, the increases of metabolites were observed at the entry point of the TCA cycle (phosphoenol pyruvate, citrate, cis-aconitate and isocitrate), which may indicate lower demand and thus, slower turn over of this pathway in P. × canescens. This assumption was corroborated by lower dark respiration in P. × canescens than in P. euphratica (Table 3).

Table 3. Photosynthesis, respiration and carbohydrates in P. euphratica and P. × canescens

thumbnailFigure 7. Analysis of pathways related to energy metabolism. Map displays selected steps from KEGG pathways dpop00010 'Glycolysis/Gluconeogenesis' and dpop00020 'Citrate cycle (TCA cycle)'. Colors indicate significant expression, respective metabolite content ratios between P. euphratica and P. × canescens, red indicates higher relative levels in P. euphratica, blue in P. × canescens. Enzymes are given as EC numbers: 1.1.1.37, malate dehydrogenase; 1.1.1.42, isocitrate dehydrogenase; 1.2.1.12, glyceraldehyde-3-phosphate dehydrogenase; 1.2.4.1, pyruvate dehydrogenase; 1.3.5.1, succinate dehydrogenase; 2.3.1.12, dihydrolipoamide S-acetyltransferase; 2.7.1.40, pyruvate kinase; 3.1.3.11, fructose-bisphosphatase; 4.1.2.13, fructose bisphosphate aldolase; 4.2.1.11, 2-phospho-D-glycerate hydrolyase; 5.4.2.1, phosphoglycerate mutase; 5.4.2.2, phosphoglucomutase.

Further analyses of carbohydrate metabolism indicated up-regulation of transcripts for genes involved in starch and carbohydrate metabolism in P. euphratica compared with P. × canescens (Figure 8). In addition, pathway analysis of metabolome data pointed to increased concentrations of glucose and fructose in P. euphratica relative to P. × canescens (Figure 8). Independent biochemical carbohydrate analyses revealed that P. euphratica contained double the glucose concentrations of P. × canescens (Table 3). Higher activity of sugar metabolism can also be inferred from six carbohydrate transporters with increased transcript abundance in P. euphratica compared with P. × canescens, whereas no transporters of this class were found to be increased in P. × canescens (Additional file 1, Table S2). The net CO2 assimilation rates were, however, lower in P. euphratica than in P. × canescens (Table 3). It is likely that this difference was mainly caused by Rubisco activase, a chaperone that maintains Rubisco activity [29], which showed significantly lower transcript abundance in P. euphratica than in P. × canescens (Additional file 1, Table S2)

thumbnailFigure 8. Pathway analysis of starch and sucrose metabolism. Map displays selected steps from KEGG pathway dpop00500 'Starch and Sucrose metabolism'. Colors indicate significant expression, respective metabolite content ratios between P. euphratica and P. × canescens, red indicates higher relative levels in P. euphratica, blue in P. × canescens. Two-colored enzymes can occur where different genes models are annotated by the same enzymatic function. Enzymes are given as EC numbers: 2.4.1.13, sucrose synthase; 2.4.1.18, 1,4-α-glucan branching enzyme; 2.4.1.34, 1,3-β-glucan synthase; 3.2.1.2, β-amylase; 3.2.1.3, glucan 1,4-α-glucosidase; 3.2.1.4, cellulase; 3.2.1.20, α-glucosidase; 3.2.1.21, β-glucosidase; 3.2.1.39, glucan endo-1,3-β-D-glucosidase.

Since carbohydrates also play roles in salt tolerance by counteracting osmotic stress, transcriptional and metabolomic differences for 'minor carbohydrates' according to Mapman categories were also studied. The raffinose and myo-inositol pathways with enzymes leading to formation of galactinol, raffinose and stachyose as well as to myo-inositol were transcriptionally increased in P. euphratica compared with P. × canescens (Additional file 1, Table S2). Increased concentrations of myo-inositol and of raffinose were detected in P. euphratica compared with P. × canescens (Additional file 1, Table S3).

Discussion

Direct transcriptome comparisons of different poplar species are affected by their phylogenetic distances

In spite of its regular application, correct use of cross-species hybridization on microarrays is still an intensively discussed topic. Mismatch as a result of single nucleotide polymorphisms can decrease the signals of single probes and influence the mean probe set signal leading to a higher standard deviation as shown here and in cases of strong sequence divergence to systematically lower signal intensities [23,30]. To avoid species-related bias when comparing the transcriptomes of Arabidopsis thaliana and its close relative the salt tolerant cress Thellungiella halophila on an Arabidopsis cDNA microarray, only genes with higher signal ratios of Thellungiella/Arabidopsis were considered [13]. However, this may lead to loss of important information with respect to the biological differences.

Recently, Renn et al. [31] showed for different cichlid fish species that the consistency of cross-species hybridization microarray results increased with decreasing phylogenetic distance of the analyzed species. For species with a time of genetic divergence of less than 10 million years (Mya) ago, transcription profiling gave most robust results, but divergence times of greater than 65 Mya were still acceptable. The genus of Populus diverged from Salix around 60 to 65 Mya ago and evolution of Populus proceeds relatively slowly, approximately at one-sixth of the rate of Arabidopsis [3]. In line with this notion sequence identity of structural genes between P. × canescens, a species from the oldest section Leuce and P. trichocarpa, which is member of the relatively young Tacamahaca section [32], was not significantly different from that between P. euphratica in the section Turanga, which assumes a phylogenetically intermediate position, and P. trichocarpa. Still, a clear species-related bias for signal ratios was observed that reflecting the genetic phylogeny of different poplar species. We, therefore, conclude that interspecific transcriptome comparisons using Poplar Genome Array between species from different sections require the exclusion of false positives by appropriate filtering procedures, whereas comparisons within a phylogenetic section are unproblematic.

Pathway mapping reveals differences in energy and secondary metabolism between P. euphratica and P. × canescens

Applying tools of systems biology, we demonstrated that pathways of energy metabolism and for the production of secondary metabolites constituted major differences between P. euphratica and P. × canescens. It has been demonstrated for species in different kingdoms including poplars that transcript-level variation is under genetic control and heritable [33-37]. The identified differences between the transcriptomes of a salt tolerant and a susceptible poplar species reflect their inherited molecular phenotypes since developmental and environmental factors, that are also influencing gene regulation, were excluded. However, in contrast to our expectation no evidence for general increases in typical pathways mediating stress tolerance were found in P. euphratica compared with P. × canescens. Similarly, transcriptional comparisons of the salt tolerant T. halleri and the susceptible A. thaliana did not show activation of general stress responses but rather pointed to a precise defence strategy [10,38,39].

It is debatable if all observed differences between P. euphratica and P. × canescens were specifically related to their differences in salt susceptibility. For example, increases in phenolics induced by chemical effectors that resulted in ameliorative effects on salt stressed plants also correlated with increases in other antioxidative compounds [40]. Therefore, the contribution of phenolic compounds to improve salt persistence is still elusive and requires further investigation. However, it is well known that phenolic compounds have antioxidant properties [41]. Since excess salinity causes oxidative stress [6], it is likely that phenol-based compounds prevent negative effects of salt-related oxidative stress by scavenging free radicals, thereby, contributing to maintain cell and tissue integrity. Notably, the phenylpropanoid metabolism was also among the significantly increased pathways identified in a comparison of drought susceptible and tolerant poplar species [37]. Since P. euphratica is not drought tolerant [42] we suggest that phenolic compounds do not constitute a primary line of defence to prevent osmotic imbalance but that their enrichment may enhance stress tolerance by downstream effects such as affording increased radical scavenging.

Further major transcriptional differences between the two species were related to carbohydrate metabolism and confirmed by analysis of carbohydrate concentrations. Enhanced glucose, fructose, raffinose, and myo-inositol levels in P. euphratica leaves suggest that this species is pre-adapted to osmotic stress. An interesting feature is that the elevated carbohydrate concentrations do not necessitate further increases in osmolytes upon exposure to excess salinity. On the contrary, we have previously shown that osmotic adjustment, which is required when sodium accumulates, is achieved by decreases in carbohydrates, especially that of glucose [43]. Furthermore, some classes of osmolytes such as sugar alcohols ameliorate the negative consequences of high ionic strength by preventing unfolding and precipitation of proteins [44-46]. Further support for the importance of chaperones for salt preparedness comes from transcriptome analyses showing that transcript levels of cyclophilins (category 'cell organisation'), heat shock proteins and chaperones (category stress) were increased in comparison with P. × canescens (Additional file 1, Table S3). It is probably an advantage of P. euphratica that protective proteins and metabolites are already in place before the stress starts. Up-regulation of these pathways is found in P. × canescens only upon salt exposure [47] and, therefore, the protective effect in salt sensitive species may lag too much behind to prevent effectively salt imposed injury.

Salt preparedness involves constitutively high expression of ion carriers

A striking and important finding is the overrepresentation of transport systems in P. euphratica compared to P. × canescens (Figure 3). A closer inspection of this category shows the presence of ABC transporters, metabolite transporters, sugar transporters and transporters with known functions in salt tolerance such as aquaporins, Ca2+, Na+, and K+ transporters [48]. For instance, gating of aquaporins is important to maintain the water balance under stress [49]. In several species, including field grown P. euphratica, the expression of aquaporins decreased under salt stress [9,50,51], whereas overexpression of aquaporins increased plant vigour only under non-stressed conditions [52]. The activation of water channels in non-stressed P. euphratica suggests higher water supply under normal conditions. This conclusion is corroborated by a higher ratio of transpiration to photosynthesis in P. euphratica than in P. × canescens supporting that P. euphratica is a phreatophytic species [53].

Maintenance of calcium and potassium supply are crucial features of salt tolerance [48]. In P. euphratica, K+ efflux was lower than in a salt sensitive poplar species involving Na+-induced down regulation of potassium outward rectifying channels and nonselective cation channels and regulated by Ca2+ [54,55]. In our study, putative Ca2+-activated outward rectifying K+ channels showed higher expression in P. × canescens than in P. euphratica, whereas several systems for K+-influx were increased in P. euphratica. The expression and response to salt stress of a large range of K+ channels has been screened in P. × canescens [16]. In leaves, none of these transport systems increased in response to salinity and some were not detectable such as the stelar K+-outward rectifying channel (SKORC) that mediates the delivery of K+ from stelar cells to the xylem in the roots towards the shoot. SKORC was neither expressed in leaves Arabidopsis thaliana [56], whereas considerable activation was found in P. euphratica leaves. The most important K+-transporter with respect to salt tolerance is probably HKT1, which was among the most highly upregulated genes (92 fold) in P. euphratica compared with P. × canescens. HKT1 is a high affinity K+ uptake transporter that facilitates Na+ homeostasis [57,58]. In Arabidopsis thaliana AtHKT1;1 controls the rate of Na+ transport from root to shoot by the retrieval of Na+ from the xylem in the roots [59] and therefore is a key factor in determining salt tolerance in various species [60]. In P. euphratica sodium transport from root to shoot is also more strongly restricted than in salt sensitive species but the molecular basis of this effect is not yet known [61]. The strong transcriptional activation of PeHKT1 suggests that this transporter is crucial for the control of Na+ in P. euphratica leaves.

NHX-type antiporters, which are involved in Na-extrusion from the cytosol back into the apoplast, or into endosomal compartments or the vacuole [62,63], were not identified in our analysis. However, NaHD1 was increased in P. euphratica compared with P. × canescens. NaHD1 is a sodium proton antiporter, whose heterolog expression in E. coli that was lacking its two sodium/proton antiporters NahA and NahB restored salt tolerance [64]. The activity of PeNaHD1 is pH dependent and functions only at low pH [64]. Maintenance of its activity therefore requires a sufficiently large proton gradient. Our comparative analysis revealed increased transcript levels of H+-ATPases suggesting maintenance of a higher proton gradient in P. euphratica. Since permanent activation of proton pumps is ATP-consuming, a higher energy demand is expected for P. euphratica. It is, therefore, reasonable to assume that the observed increases in mitochondrial respiration in P. euphratica compared with P. × canescens are required as driving force for a tighter control of ion compartmentalization.

Conclusions

We showed that direct comparisons of the transcriptomes of different poplar species on Poplar Gene Array are influenced by their phylogenetic distance, if the species are from different sections as applies for the salt sensitive P. × canescens (Leuce) and the salt tolerant P. euphratica (Turanga). This influence can be reduced by application of appropriate filters. Using filtered transcriptome data, pathways for the production of secondary metabolites, and carbohydrates as well as a significant enrichment of genes in the categories 'transporters' and 'mitochondria' were identified that were activated in P. euphratica compared with P. × canescens. This interpretation was corroborated by a combination of metabolomic analyses and targeted biochemical and physiological measurements suggesting that higher metabolic flux through the pathways resulted in the observed increases of phenolic end products and carbohydrates in P. euphratica compared with P. × canescens. Tolerance mechanisms for excess salinity in P. euphratica are obviously relying on stronger pre-formed osmolyte and chaperone abundance and on tighter control of ion compartmentalization requiring increased energy allocation.

Underrepresentation of GO terms related to signaling and transcription suggest that mechanisms responsible for gene regulation and transcription were less active in P. euphratica. Since previous analyses revealed only moderate transcriptomic changes in response to stress [9], we propose that the underside to a high level of stress preparedness as a result of adaptation to a specific ecological niche is a loss of flexibility and adjustability of the transcriptome to environmental cues.

The results of our study may also have implications for wood production. Poplars are fast fast-growing tree species that constitute an important resource for woody biomass. However, expansion of tree cultivation will only be possible on marginal or polluted soils that are not suitable for food production [65]. For this purpose stress tolerant tree species are required. The results of this study laid a foundation for understanding tree salt tolerance at the level of system biology. This knowledge can be used for tree improvement by biotechnological methods or conventional breeding.

Methods

Plant culture

Plantlets from P. × canescens (P. alba × tremula) clone INRA717 1-B4 and P. euphratica clone B2 from the Ein Avdat valley in Israel were multiplied by in vitro micropropagation [66], transferred to aerated hydroponics using Long Ashton nutrient solution [67] for acclimation to greenhouse conditions, subsequently potted into soil (Frühstorfer Erde, Lauterbach-Wallenrod, Germany) and grown for 3 months in a greenhouse at 20 to 25°C, with a relative air humidity of 40 to 60% and additional 16 h of light (neon lamps: L58W/25 and 58W/840, Osram, Munich, Germany, and TLD 58W/840 Philips, Amsterdam, Netherlands) between 6:00 a.m. and 10:00 p.m. to reach a continuous photosynthetic active radiation (PAR) of 150 μmol m-2s-1. Fully expanded leaves were harvested at a stem height of approximately 70 cm (upper third of the stem) after four h in the light period. The leaves were shock-frozen in liquid nitrogen and kept at -80°C until further analysis.

Gas exchange

Gas exchange (net photosynthesis, respiration, and transpiration) was measured (1 min means) with a portable gas exchange system (HCM-1000 or GSF-3000, Walz, Effeltrich, Germany) under ambient conditions. For analysis of dark respiration plants were connected to the system and kept for 30 min in darkness prior to analysis. For each plant (n = 5) the gas exchange values of 5 min were averaged.

FT-ICR-MS measurements

Of each species, 5 biological replicates were measured. Frozen leaf tissue was ground in a ball mill (Retsch, Haan, Germany). Fine powder was diluted in methanol to a methanolic concentration of 70% to give highest ion density inside the electrospray, without eliminating those neutrals which are highly water soluble. High-resolution mass spectra were acquired on a Fourier Transform Ion Cyclotron Resonance Mass Spectrometer (APEX Qe, Bruker, Bremen, Germany) equipped with a 12-T superconducting magnet and an Apollo II Electrospray (ESI) source. The ionization source was run in the negative operation mode to generate mono charged negative ions. Thus, all important multi-functional organic compounds which bear at least one acidic site could be deprotonated to give rise to a detectable anion in the mass spectrometer [18]. Each sample was introduced into the ionization source at a flow rate of 2 μL min-1 by a microliter pump with a nebulizer gas pressure of 20 psi and a drying gas pressure of 15 psi (heated to 200°C). Each sample was measured three times and the signal intensities of each detected ion were averaged before running the multi-dimensional statistical data analysis.

Spectra were externally calibrated on clusters of arginine (10 mg l-1 in methanol); calibration errors in the relevant mass range were always below 0.05 ppmv. The spectra were acquired with a time domain of 1 Megaword with a mass-to-charge ratio (m/z) range of 146-2000. The spectra were zero filled to a processing size of 2 Megawords. Before Fourier transformation of the time-domain transient, a sine apodization was performed. No fragmentation experiments were performed in this study. Thus, the whole mass range could be scanned, and 300 scans were summed on in each acquisition. The ion accumulation time in the ion source was set to 0.1 s.

FT-ICR mass spectra were exported to peak lists at a signal to noise ratio (S/N) of two. From those lists, High to Low signal intensity and Low to High signal intensity profiles were obtained by use of the Hierarchical Clustering Explorer HCE Version 3.0 [68]. Ions with sharp differences between the acquisitions of the two species were highlighted as High to Low and Low to High signal intensity profiles. The m/z values of all ions were clustered by the use of the average linkage method with similarity/difference measure type: Euclidean distance. The search method is model-based and the used distance measure was Pearson's correlation coefficient of 0.8. For all obtained hits, t-tests were calculated to validate a statistical significance at α = 0.1.

Metabolite identification and annotation of the m/z values was performed via the MassTRIX web site with the following parameter settings: 'Scan mode negative ionisation (correct for H+ loss)', 'Max. error 1.0 ppm', 'Database KEGG with isotopes', 'Organism Populus trichocarpa' [22]. Identified compounds and differentially expressed genes added to the query as EC numbers were automatically mapped to KEGG pathways using KEGG/API [69]via MassTrix. Metabolites with the same total formula could not be distinguished.

Analysis of soluble and cell wall bound phenolics

Frozen leaf tissue was ground in a ball mill (Retsch, Haan, Germany). Fine powder (60 mg per sample) was extracted with 2 ml of 50% methanol in an ultrasonic bath (60 min, 40°C; Sonorex Super RK 510 H, Bandelin electronics, Berlin, Germany). The extract was centrifuged, the pellet was extracted once again and the supernatants were combined for analysis of soluble phenolics photometrically with the Folin Ciocalteus method [70]. Catechin (Sigma-Aldrich, Deisenhofen, Germany) was used for calibration and the phenolic concentrations were expressed as catechin equivalents.

To determine cell wall bound phenolics, the pellet was washed with n-hexane, dried, and weighed. This pellet representing cell walls was extracted with 1 M NaOH as described earlier [70]. The extracted phenolic compounds were measured with the Folin Ciocalteous reagent as above.

Analysis of soluble and non-soluble condensed tannins

Condensed tannins (proanthocyanidin) were determined after [71]. For extraction of soluble tannins, plant materials ground to a fine powder were extracted with 50% methanol, centrifuged, the pellet was re-extracted and the extracted solutions were combined. Aliquots were mixed with a reagent (1.6% methanol, acid butanol [5% concentrated HCl in butanol (v/v)] and 0.2 ml of 2% FeNH4(SO4)2·12H2O in 2 M HCl [FeNH4(SO4)2·12H2O/HCl = w/v], incubated at 95°C and used to measure the extinction at 550 nm. The amount of tannins was expressed as absorbance per dry mass.

To determine non-soluble tannins, the pellet was lyophilized, weighed and resuspended in 1 ml of 100% methanol. The slurry was combined with acid butanol reagent and processed as above.

Soluble carbohydrates

Frozen tissues were ground to a fine powder with a ball mill (Retsch, Haan, Germany) pre-cooled in liquid nitrogen and extracted in DMSO/HCl (dimethylsulfoxid:25% HCl = 80:20 (v:v)) at 60°C for 30 min. The mixture was centrifuged and aliquots of the supernatant were mixed with 0.2 M citrate buffer (pH 10.6). This extract was mixed with a reaction mixture (0.75 M triethanolamine-buffer, 4 mM NADP, 10 mM ATP, 9 mM MgSO4, pH 7.6) for enzymatic determination of soluble sugars as described previously [72]. Glucose was measured after addition of hexokinase (3 mg ml-1, Roche Diognostics GmbH, Mannheim, Germany), fructose after addition of phosphoglucose isomerase (10 mg ml-1, Roche Diognostics GmbH, Mannheim, Germany) and sucrose after incubation with β-fructosidase (5 mg ml-1, Roche Diognostics GmbH, Mannheim, Germany). Calibration curves were produced with glucose, fructose and sucrose. The sucrose content of the extracts was below the detection limit of the method.

In silico sequence analysis

DNA sequences of twenty genes were compared between P. × canescens, P. euphratica and P. trichocarpa. For five genes, full length sequences were available for P. euphratica and P. × canescens cDNA ('IAA-amino acid hydrolase', Ill3 [[73] GenBank:AJ744952, GenBank:AJ744953]; 'Na/H antiporter type D', NhaD1 [[74] GenBank:AJ561195]; 'salt induced serine-rich' and 'temperature induced lipocalin', SIS and TIL [[74] GenBank:FJ238515, GenBank:FJ238514 and GenBank:FJ238513, GenBank:FJ238512]; 'gibberellin regulated protein', GAST, [GenBank:FJ238511, GenBank:CF231013]). Criterion for the selection of further genes was their representation by probe sets on the microarray and the availability of sequence data for P. euphratica and P. × canescens in public databases. Probe sets both with and without significantly different signal intensities between P. euphratica and P. × canescens were randomly chosen, and sequence data for the corresponding P. trichocarpa gene model retrieved from the Populus genome project (Table 1) [28]. Because most sequence data is available from EST datasets and thus as cDNA sequences, open reading frame (ORF) sequences of the P. trichocarpa gene models were used to screen the GenBank database for homologs in P. euphratica and P. × canescens by using the Megablast algorithm [75]. Pair wise alignments and base identity calculations of the cDNA sequences were generated with the GeneDoc software [76]. An ANOVA was conducted to test for significant differences between the sequence identity values at α = 0.05.

RNA extraction

Leaf material of 3 plants was pooled. Three biological replicates per poplar species were analyzed (i.e. a total of 9 nine plants per species was included). Frozen leaf tissue was milled in a ball mill (Retsch, Haan, Germany). Total RNA was extracted from 500 mg of frozen plant powder according to Chang et al. [77] with minor modifications: No spermidine was applied in the extraction buffer, and 2% β-mercaptoethanol was used. RNA was additionally purified using an RNeasy Mini Kit (Qiagen, Valencia, CA). Total RNA yield and purity were determined spectrophotometrically (BioPhotometer, Eppendorf, Hamburg, Germany) at A260 and A280. RNA integrity was assessed on an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) at the Microarray Facility (Tübingen, Germany).

Microarray analysis

Of each species, 3 biological replicates were analyzed on the GeneChip® Poplar Genome Array (Affymetrix, Santa Clara, CA). Synthesis of one-cycle cDNA and biotin-labelled cRNA, fragmenting of cRNA, hybridization to the Poplar Genome Array, washing, staining and scanning was performed as stated by Affymetrix (GeneChip® Expression Analysis Technical Manual) at the Microarray Facility (Tübingen, Germany). Data were deposited at ArrayExpress [EMBL: E-MEXP-1928].

Statistical analysis of the raw signal intensity data was conducted using the following functions from packages released by the bioconductor project, implemented in R [78,79]. Background correction, quantile normalization and summarization of the Affymetrix CEL output files resulting in a raw list of normalized genes was computed by robust multi-array average (RMA) using the 'rma' function from the 'affy' package [23]. Statistical testing for differentially expressed genes was performed on this raw list with the 'sam' function from the 'siggenes' package [80]. In the SAM analysis, Δ was chosen to obtain FDR = 0.05. Annotation of the genes was carried out via the PopArray-Database [81]. Gene Ontology (GO) terms were matched via The Arabidopsis Information Resource (TAIR) [82].

Gene lists were subjected to several filtering processes after [25]: genes that were not expressed were removed by calculating Affymetrix' MAS 5.0 change calls using the 'mas5calls' function from the 'affy' package; only genes with two or three 'present' calls in any one of the two species were considered. Probe sets presumably affected by differences in binding efficiencies were removed. For this purpose, SDs of prey signals were calculated for each bait set on a microarray and then averaged over all six microarrays for each gene. The resulting data were compared with the range of the SDs of all prey set signals. A gene was removed from the list when its prey signal SD was in the upper 5% quantile of all 'present' bait set SDs on the microarray. The setting of the threshold for this filter affects the stringency of the results. For the purpose of this study the upper 5%-quantile was used.

Probe sets for which no annotation was available were removed. In cases where multiple probe sets were matched with the same JGI P. trichocarpa gene model, duplicate probe sets were removed so that each gene model was represented only once. The probe set with the lower p-value in the SAM analysis was used for further analysis if both genes showed the same direction of change in gene expression. Probe sets annotated by the same gene model with contrasting expression ratios were removed.

Bioinformatic analyses

For evaluation of over- and underrepresented GO slim terms, the frequencies of annotations within list of significant genes were compared to the frequencies within the defined 'gene universe'. Statistical testing was conducted by calculating the cumulative hypergeometric distribution function, using the 'phyper' function in R. Benjamini-Hochberg correction was applied to the resulting p-values using the 'p.adjust' function.

For statistical analysis of overrepresentation with 'The Ontologizer' [83], a gene ontology file for Populus was adapted. In this file, all JGI Populus trichocarpa gene models represented on the Affymetrix GeneChip were linked to the GO identifier of their closest Arabidopsis matches as listed by TAIR [82]. The 'population' of the Ontologizer conforms with our 'gene universe', whereas the two lists of genes higher expressed in one of the poplar species each correspond to a 'study set'. As parameter settings, term-for-term analysis with Benjamini-Hochberg correction was used.

qRT-PCR

For quantitative Real-Time PCR (qRT-PCR), primer pairs were designed for actin 9 [GenBank:AJ778775 (P. euphratica) and GenBank:CX656348 (P. × canescens)] as a reference gene and 10 transcripts that were differentially expressed in P. euphratica and P. × canescens (Additional file 1, Table S4). All primers met the following conditions: 1) 100% sequence identity with both P. euphratica and P. × canescens cDNA; 2) fragment lengths between 90 and 180 base pairs (bp); 3) calculated salt-adjusted melting temperatures Tm(salt) between 58 and 62°C. Primer design was performed with the Oligo Explorer, and suitable primers were tested for similar Tm(salt), primer dimers and primer loops by Oligo Analyzer (both Gene Link, Hawthorne, NY, USA).

Total RNA was DNAse treated with a Turbo DNA-free kit (Ambion, Austin, TX) and transcribed to cDNA with a RevertAid™ First Strand cDNA Synthesis Kit (MBI Fermentas, St. Leon-Rot, Germany). For each gene, three technical repeats were analyzed for three biological samples. The qRT-PCR was performed on an iCycler (Bio-Rad, Hercules, CA). 250 ng cDNA were used in a 25 μl reaction with 1 × ABsolute qPCR SYBR Green Fluorescein Mix (ABgene, Surrey, UK; including Thermo-Start DNA Polymerase) and 10 μM primer. First denaturation and activation of the Taq-polymerase occurred at 95°C for 15 min, followed by 45 cycles of denaturation at 95°C for 15 s, annealing at 56°C for 30 s and elongation at 72°C for 30 s.

qRT-PCR outputs of three biological replicates per species and gene were analyzed using the MyiQ software (Bio-Rad, Hercules, CA). Primer specificity was assessed by melting curve analysis. Statistical analysis of expression ratios and standard error was conducted using a pair wise fixed reallocation randomization test implemented in Excel with the Relative Expression Software Tool (REST)-384 [84,85]. All primer pairs showed less than 5% differences in PCR efficiency between P. euphratica and P. × canescens.

Statistical analyses

Data are shown as means (± SE). Five to six plants were analyzed per species. Statistical analysis (ANOVA) was performed using Statgraphics (Centurion XV, St Louis, MO). Significant differences at P ≤ 0.05 were calculated by ANOVA followed by a Tukey-Kramer test [86].

Abbreviations

ANOVA - analysis of variance, bp - base pairs, cDNA - complementary DNA, DNA - deoxyribonucleic acid, EST - expressed sequence tag, FT-ICR/MS - Fourier transform-ion cyclotron resonance mass spectrometry, GO - gene ontology, IQL - interquartile length, JGI - Joint Genome Institute, KEGG - Kyoto Encyclopedia of Genes and Genomes, m/z - mass-to-charge ratio, Mya - million years, ORF - open reading frame, P. - Populus, ppm - parts per million, qRT-PCR - quantitative real-time polymerase chain reaction, REST - relative expression software tool, RMA - robust multi-array average, SAM - significance analysis of microarrays, TAIR - The Arabidopsis Information Resource

Authors' contributions

DJ conducted this work as part of his PhD thesis. He performed the experiments, did the bioinformatic analyses and drafted the manuscript. PSK and BK provided FT-ICR-MS data and conducted the hierarchical clustering analysis. KB and JPS conducted gas exchange analyses. AP designed the experiments and finalized the study. All authors read and commented on the manuscript.

Acknowledgements

We are grateful to the German Science foundation for funding Poplar research Group Germany (FOR 496, PO362-12, PO361-13, SCHN653-3, SCHN653-4). We acknowledge excellent technical support by S. Elend, C. Kettner and T. Klein (Laboratory for Radioisotopes).

References

  1. Kundzewicz ZW, Mata LJ, Arnell NW, Döll P, Kabat P, Jiménez B, Miller KA, Oki T, Sen Z, Shiklomanov IA: Freshwater resources and their management. In Climate Change 2007: Impacts, Adaptation and Vulnerability.Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Edited by Parry ML, Canziani OF, Palutikof JP, van der Linden PJ, Hanson CE. Cambridge, UK: Cambridge University Press; 2007:173-210. OpenURL

  2. Brunner AM, Busov VB, Strauss SH: Poplar genome sequence: functional genomics in an ecologically dominant plant species.

    Trend Plant Sci 2004, 9:9-56. Publisher Full Text OpenURL

  3. Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, et al.: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray).

    Science 2006, 313:1596-1604. PubMed Abstract | Publisher Full Text OpenURL

  4. Jansson S, Douglas CJ: Populus: A model system for plant biology.

    Annu Rev Plant Biol 2007, 58:435-458. PubMed Abstract | Publisher Full Text OpenURL

  5. Feng Q, Endo KN, Cheng GD: Towards sustainable development of the environmentally degraded arid rivers of China - a case study from Tarim River.

    Environm Geo 2001, 41:229-238. OpenURL

  6. Chen SL, Polle A: Salt tolerance in Populus.

    Plant Biol 2010, 12:317-333. PubMed Abstract | Publisher Full Text OpenURL

  7. Wang RG, Chen SL, Deng L, Fritz E, Huttermann A, Polle A: Leaf photosynthesis, fluorescence response to salinity and the relevance to chloroplast salt compartmentation and anti-oxidative stress in two poplars.

    Trees 2007, 21:581-591. Publisher Full Text OpenURL

  8. Gu RS, Fonseca S, Puskas LG, Hackler L, Zvara A, Dudits D, Pais MS: Transcript identification and profiling during salt stress and recovery of Populus euphratica.

    Tree Physiol 2004, 24:265-276. PubMed Abstract | Publisher Full Text OpenURL

  9. Brosché M, Vinocur B, Alatalo ER, Lamminmaki A, Teichmann T, Ottow EA, Djilianov D, Afif D, Bogeat-Triboulot MB, Altman A, Polle A, Dreyer E, Rudd S, Lars P, Auvinen P, Kangasjarvi J: Gene expression and metabolite profiling of Populus euphratica growing in the Negev desert.

    Genome Biol 2005, 6:R101. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  10. Gong Q, Li P, Ma S, Rupassara SU, Bohnert HJ: Salinity stress adaptation competence in the extremophile Thellungiella halophila in comparison with its relative Arabidopsis thaliana.

    Plant J 2005, 44:826-839. PubMed Abstract | Publisher Full Text OpenURL

  11. Amtmann A: Learning from evolution: Thellungiella generates new knowledge on essential and critical components of abiotic stress tolerance in plants.

    Mol Plant 2009, 2:3-12. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Bar-Or C, Czosnek H, Koltai H: Cross-species microarray hybridizations: a developing tool for studying species diversity.

    Trend Genet 2007, 23:200-207. Publisher Full Text OpenURL

  13. Taji T, Seki M, Satou M, Sakurai T, Kobayashi M, Ishiyama K, Narusaka Y, Narusaka M, Zhu JK, Shinozaki K: Comparative genomics in salt tolerance between Arabidopsis and Arabidopsis-related halophyte salt cress using Arabidopsis microarray.

    Plant Physiol 2004, 135:1697-1709. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Ehlting B, Dluzniewska P, Dietrich H, Selle A, Teuber M, Hänsch R, Nehls U, Polle A, Schnitzler JP, Rennenberg H, Gessler A: Interaction of nitrogen nutrition and salinity in grey poplar (Populus tremula × alba).

    Plant Cell Environm 2007, 30:796-811. Publisher Full Text OpenURL

  15. Teuber M, Zimmer I, Kreuzwieser J, Ache P, Polle A, Rennenberg H, Schnitzler JP: VOC emission of Grey poplar leaves as affected by salt stress and different N sources.

    Plant Biol 2008, 10:86-96. PubMed Abstract OpenURL

  16. Escalante-Pérez M, Lautner S, Nehls U, Selle A, Teuber M, Schnitzler JP, Teichmann T, Fayyaz P, Hartung W, Polle A, Fromm J, Hedrich R, Ache P: Salt stress affects xylem differentiation of grey poplar (Populus × canescens).

    Planta 2009, 229:299-309. PubMed Abstract | Publisher Full Text OpenURL

  17. Affymetrix GeneChip Poplar Genome Array Data Sheet [http://www.affymetrix.com/support/technical/datasheets/poplar_datasheet.pdf] webcite

  18. Schmitt-Kopplin P, Hertkorn N: Ultrahigh resolution mass spectrometry.

    Anal Bioanal Chem 2007, 389:1309-1310. Publisher Full Text OpenURL

  19. Gougeon RD, Lucio M, Frommberger M, Peyron D, Chassagne D, Alexandre H, Feuillat F, Cayot Ph, Voilley A, Gebefügi I, Hertkorn N, Schmitt-Kopplin Ph: Fom the cover: The chemodiversity of wines can reveal a metabologeography of cooperage oak wood.

    Proc Natl Acad Sci USA 2009, 106:9174-9179. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Leon C, Rodriguez-Meizoso I, Lucio M, Garcia-Cañas V, Ibañez E, Schmitt-Kopplin Ph, Cifuentes A: Metabolomics of transgenic foods combining Fourier Transform-Ion Cyclotron Resonance-Mass Spectrometry, Capillary Electrophoresis-Mass Spectrometry and Pressurized Liquid Extraction.

    J Chrom A 2009, 1216:7314-7323. Publisher Full Text OpenURL

  21. Hirai MY, Yano M, Goodenowe DB, Kanaya S, Kimura T, Awazuhara M, Arita M, Fujiwara T, Saito K: Integration of transcriptomics and metabolomics for understanding of global responses to nutritional stresses in Arabidopsis thaliana.

    Proc Natl Acad Sci USA 2004, 101:10205-10210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Suhre K, Schmitt-Kopplin P: MassTRIX: mass translator into pathways.

    Nuc Acid Res 2008, 36:W481-W484. Publisher Full Text OpenURL

  23. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data.

    Biostatistics 2003, 4:249-264. PubMed Abstract | Publisher Full Text OpenURL

  24. Gilad Y, Rifkin SA, Bertone P, Gerstein M, White KP: Multi-species microarrays reveal the effect of sequence divergence on gene expression profiles.

    Genome Res 2005, 15:674-680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association.

    Bioinformatics 2007, 23:257-258. PubMed Abstract | Publisher Full Text OpenURL

  26. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology.

    Nature Gen 2000, 25:25-29. Publisher Full Text OpenURL

  27. Thimm O, Blaesing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes.

    Plant J 2004, 37:914-939. PubMed Abstract | Publisher Full Text OpenURL

  28. JGI Populus trichocarpa genome release 1.1 [http://genome.jgi-psf.org/Poptr1_1/Poptr1_1.home.html] webcite

  29. Salvucci ME, Osteryoung KW, Crafts-Brandner SJ, Vierling E: Exceptional sensitivity of Rubisco activase to thermal denaturation in vitro and in vivo.

    Plant Physiol 2001, 127:1053-1064. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Rouchka EC, Phatak AW, Singh AV: Effect of single nucleotide polymorphisms on Affymetrix(R) match-mismatch probe pairs.

    Bioinformation 2008, 2:405-411. PubMed Abstract | PubMed Central Full Text OpenURL

  31. Renn SCP, Aubin-Horth N, Hofmann HA: Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray.

    BMC Genomics 2004, 5:42. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  32. Cervera MT, Storme V, Soto A, Ivens B, Van Montagu M, Rajora OP, Boerjan W: Intraspecific and interspecific genetic and phylogenetic relationships in the genus Populus based on AFLP markers.

    Theor Appl Gen 2005, 111:1440-1456. Publisher Full Text OpenURL

  33. Brem RB, Yvert G, Clinton R, Kruglyak L: Genetic dissection of transcriptional regulation in budding yeast.

    Science 2002, 296:752-755. PubMed Abstract | Publisher Full Text OpenURL

  34. Wayne ML, McIntyre LM: Combining mapping and arraying: an approach to candidate gene identification.

    Proc Natl Acad Sci USA 2002, 99:14903-14906. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, Colinayo V, Ruff TG, Milligan SB, Lamb JR, Cavet G, Linsley PS, Mao M, Stoughton RB, Friend SH: Genetics of gene expression surveyed in maize, mouse and man.

    Nature 422:297-302. PubMed Abstract | Publisher Full Text OpenURL

  36. Kirst M, Basten CJ, Myburg AA, Zeng ZB, Sederoff RR: Genetic architecture of transcript-level variation in differentiating xylem of a Eucalyptus hybrid.

    Genetics 2005, 169:2295-2303. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Street NR, Skogstrom O, Sjodin A, Tucker J, Rodriguez-Acosta M, Nilsson P, Jansson S, Taylor G: The genetics and genomics of the drought response in Populus.

    Plant J 2006, 48:321-341. PubMed Abstract | Publisher Full Text OpenURL

  38. Ma S, Bohnert H: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell specific expression.

    Genome Biol 2007, 8:R49. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  39. Ma S, Gong Q, Bohnert H: Dissecting salt stress pathways.

    J Exp Bot 2006, 57:1097-1107. PubMed Abstract | Publisher Full Text OpenURL

  40. Muthukumarasamy M, Dutta Gupta S, Panneerselvam R: Enhancement of peroxidase, polyphenol oxidase and superoxide dismutase activities by Triadimefon in NaCl stressed Raphanus sativus L.

    Biol Plant 2000, 43:317-320. Publisher Full Text OpenURL

  41. Rice-Evans C, Miller N, Paganga G: Antioxidant properties of phenolic compounds.

    Trend Plant Sci 1997, 2:152-159. Publisher Full Text OpenURL

  42. Bogeat-Triboulot M, Brosché M, Renaut J, Jouve L, Le Thiec D, Fayyaz P, Vinocur B, Witters E, Laukens K, Teichmann T, Hausman JF, Polle A, Kangasjärvi J, Dreyer E: Gradual soil water depletion results in reversible changes of gene expression, protein profiles, ecophysiology and growth performance in Populus euphratica, a poplar growing in arid regions.

    Plant Physiol 2007, 143:876-892. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Ottow EA, Brinker M, Teichmann T, Fritz E, Kaiser W, Brosché M, Kangasjärvi J, Jiang X, Polle A: Populus euphratica displays apoplastic sodium accumulation, osmotic adjustment by decreases in calcium and soluble carbohydrates, and develops leaf succulence under salt stress.

    Plant Physiol 2005, 139:1762-1772. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Welch WJ, Brown CR: Influence of molecular and chemical chaperones on protein folding.

    Cell Stress Chaperones 1996, 1:109-115. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Vinocur B, Altman A: Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations.

    Curr Op Biotechn 2005, 16:123-132. Publisher Full Text OpenURL

  46. Bressan RA, Bohnert HJ, Hasegawa PM: Genetic engineering for salinity stress tolerance.

    Adv Plant Biochem Mol Biol 2008, 1:347-384. Publisher Full Text OpenURL

  47. Luo ZB, Janz D, Jiang XN, Göbel C, Wildhagen H, Tan YP, Rennenberg H, Feussner I, Polle A: Upgrading root physiology for stress tolerance by ectomycorrhizas: insights from metabolite and transcriptional profiling into reprogramming for stress anticipation.

    Plant Physiol 2009, 151(4):1902-1817. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Munns R, Tester M: Mechanisms of salinity tolerance.

    Annu Rev Plant Biol 2008, 59:651-681. PubMed Abstract | Publisher Full Text OpenURL

  49. Chaumont F, Moshelion M, Daniels MJ: Regulation of plant aquaporin activity.

    Biol Cell 2005, 97:749-764. PubMed Abstract | Publisher Full Text OpenURL

  50. Ouziad F, Wilde P, Schmelzer E, Hildebrandt U, Bothe H: Analysis of expression of aquaporins and Na+/H+ transporters in tomato colonized by arbuscular mycorrhizal fungi and affected by salt stress.

    Environm Exp Bot 2006, 57:177-186. Publisher Full Text OpenURL

  51. López-Berenguer C, Martínez-Ballesta MA, García-Viguera C, Carvajal M: Leaf water balance mediated by aquaporins under salt stress and associated glucosinolate synthesis in broccoli.

    Plant Sci 2008, 174:321-328. Publisher Full Text OpenURL

  52. Aharon R, Shahak Y, Wininger S, Bendov R, Kapulnik Y, Galili G: Overexpression of a plasma membrane aquaporin in transgenic tobacco improves plant vigor under favorable growth conditions but not under drought or salt stress.

    Plant Cell 2003, 15:439-447. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Gries D, Zeng F, Foetzki A, Arndt SK, Bruelheide H, Thomas FM, Zhang X, Runge M: Growth and water relations of Tamarix ramosissima and Populus euphratica on Taklamakan desert dunes in relation to depth to a permanent water table.

    Plant Cell Environm 2003, 26:725-736. Publisher Full Text OpenURL

  54. Sun J, Chen S, Dai S, Wang R, Li N, Shen X, Zhou X, Lu C, Zheng X, Hu Z, Zhang Z, Song J, Xu Y: NaCl-Induced alternations of cellular and tissue ion fluxes in roots of salt-resistant and salt-sensitive poplar species.

    Plant Physiol 2009, 149:1141-1153. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Sun J, Dai S, Wang R, Chen S, Li N, Zhou X, Lu C, Shen X, Zheng X, Hu Z, Zhang Z, Song J, Xu Y: Calcium mediates root K+/Na+ homeostasis in poplar species differing in salt tolerance.

    Tree Physiol 2009, 29:1175-1186. PubMed Abstract | Publisher Full Text OpenURL

  56. Gaymard F, Pilot G, Lacombe B, Bouchez D, Bruneau D, Boucherez J, Michaux-Ferrière N, Thibaud JB, Sentenac H: Identication and disruption of a plant shaker-like outward channel involved in K+ release into the xylem sap.

    Cell 1998, 94:647-655. PubMed Abstract | Publisher Full Text OpenURL

  57. Mäser P, Eckelman B, Vaidyanathan R, Horie T, Fairbairn D, Kubo M, Yamagami M, Yamaguchi K, Nishimura M, Uozumi N: Altered shoot/root Na+ distribution and bifurcating salt sensitivity in Arabidopsis by genetic disruption of the Na+ transporter AtHKT1.

    FEBS Lett 2002, 531:157-161. PubMed Abstract | Publisher Full Text OpenURL

  58. Rus A, Lee BH, Munoz-Mayor A, Miura K, Bressan RA, Hasegawa PM: AtHKT1 facilitates Na+ homeostasis and K+ nutrition in planta.

    Plant Physiol 2004, 136:2500-2511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Davenport RJ, Munoz-Mayor A, Jha D, Essah PA, Rus A, Tester M: The Na+ transporter AtHKT1;1 controls retrieval of Na+ from the xylem in Arabidopsis.

    Plant Cell Environm 2007, 30:497-507. Publisher Full Text OpenURL

  60. Huang S, Spielmeyer W, Lagudah ES, Munns R: Comparative mapping of HKT genes in wheat, barley, and rice, key determinants of Na+ transport, and salt tolerance.

    J Exp Bot 2008, 59:927-937. PubMed Abstract | Publisher Full Text OpenURL

  61. Chen S, Li J, Fritz E, Wang S, Hüttermann A: Sodium and chloride distribution in roots and transport in three poplar genotypes under increasing NaCl stress.

    For Ecol Manag 2002, 168:217-230. Publisher Full Text OpenURL

  62. Quintero FJ, Blatt MR, Pardo JM: Functional conservation between yeast and plant endosomal Na+/H+ antiporters.

    FEBS Lett 2000, 471:224-228. PubMed Abstract | Publisher Full Text OpenURL

  63. Zhu JK: Salt and drought stress signal transduction in plants.

    Annu Rev Plant Biol 2002, 53:247-273. PubMed Abstract | Publisher Full Text OpenURL

  64. Ottow E, Polle A, Brosche M, Kangasjarvi J, Dibrov P, Zörb C, Teichmann T: Molecular characterisation of PeNHaD1: the first member of the NaHD Na+/H+ antiporter family of plant origin.

    Plant Mol Biol 2005, 58:75-88. PubMed Abstract | Publisher Full Text OpenURL

  65. Polle A, Douglas C: The molecular physiology of poplars: paving the way for knowledge-based biomass production.

    Plant Biol 2010, 12:239-241. PubMed Abstract | Publisher Full Text OpenURL

  66. Leplé JCC, Brasileiro ACM, Michel MF, Delmotte F, Jouanin L: Transgenic poplars - expression of chimeric genes using 4 different constructs.

    Plant Cell Rep 1992, 11:137-141. Publisher Full Text OpenURL

  67. Hewitt EJ, Smith TA: Plant mineral nutrition. London, UK: English U.P.; 1975.

  68. Seo J, Shneiderman B: A rank-by-feature framework for interactive exploration of multidimensional data.

    Information Visualization 2005, 4:96-113. Publisher Full Text OpenURL

  69. KEGG API [http://www.genome.jp/kegg/soap/] webcite

  70. Luo ZB, Calfapietra C, Scarascia-Mugnozza G, Liberloo M, Polle A: Carbon-based secondary metabolites and internal nitrogen pools in Populus nigra under Free Air CO2 Enrichment (FACE) and nitrogen fertilisation.

    Plant Soil 2008, 304:45-57. Publisher Full Text OpenURL

  71. Booker FL, Anttonen S, Heagle AS: Catechin, proanthocyanidin and lignin concentrations of loblolly pine (Pinus taeda) needles after chronic exposure to ozone.

    New Phytol 1996, 132:483-492. Publisher Full Text OpenURL

  72. 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.

    Global Change Biol 2006, 12:272-283. Publisher Full Text OpenURL

  73. Junghans U, Polle A, Duchting P, Weiler E, Kuhlman B, Gruber F, Teichmann T: Adaptation to high salinity in poplar involves changes in xylem anatomy and auxin physiology.

    Plant Cell Environm 2006, 29:1519-1531. Publisher Full Text OpenURL

  74. Fayyaz P: Effects of salt stress on ecophysiological and molecular characteristics of Populus euphratica Oliv., Populus × canescens (Aiton) Sm. and Arabidopsis thaliana L. In Ph.D. Thesis. Georg-August-Universität Göttingen; 2007. OpenURL

  75. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW: GenBank.

    Nucl Acid Res 2009, 37:D26-D31. Publisher Full Text OpenURL

  76. GeneDoc Homepage [http://www.nrbsc.org/gfx/genedoc/index.html] webcite

  77. Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees.

    Plant Mol Biol Rep 1993, 11:113-116. Publisher Full Text OpenURL

  78. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge YC, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang JH: Bioconductor: open software development for computational biology and bioinformatics.

    Genome Biol 2004, 5:R80. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  79. The R project for statistical computing [http://www.r-project.org/] webcite

  80. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response.

    Proc Natl Acad Sci USA 2001, 98:5116-5121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. PopARRAY: A database for cross-reference of multiple microarray platforms in Populus [http://aspendb.uga.edu/poparray] webcite

  82. The Arabidopsis Information Resource [http://www.arabidopsis.org/] webcite

  83. Bauer S, Grossmann S, Vingron M, Robinson PN: Ontologizer 2.0 - a multifunctional tool for GO term enrichment analysis and data exploration.

    Bioinformatics 2008, 24:1650-1651. PubMed Abstract | Publisher Full Text OpenURL

  84. Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR.

    Nucleic Acid Res 2001, 29:e45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  85. Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST (c)) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.

    Nucleic Acid Res 2002, 30:e36. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  86. Kramer CY: Extension of multiple range tests to group means with unequal numbers of replications.

    Biometrics 1956, 12:307-310. Publisher Full Text OpenURL

  87. Ottow EA: Molecular and ecophysiological responses of Populus euphratica (Oliv.) and Arabidopsis thaliana (L.) to salt stress. PhD thesis. Georg-August-Universität Göttingen; 2005. OpenURL