Comparative transcriptomics of drought responses in Populus: a meta-analysis of genome-wide expression profiling in mature leaves and root apices across two genotypes
© Cohen et al; licensee BioMed Central Ltd. 2010
Received: 8 June 2010
Accepted: 12 November 2010
Published: 12 November 2010
Comparative genomics has emerged as a promising means of unravelling the molecular networks underlying complex traits such as drought tolerance. Here we assess the genotype-dependent component of the drought-induced transcriptome response in two poplar genotypes differing in drought tolerance. Drought-induced responses were analysed in leaves and root apices and were compared with available transcriptome data from other Populus species.
Using a multi-species designed microarray, a genomic DNA-based selection of probesets provided an unambiguous between-genotype comparison. Analyses of functional group enrichment enabled the extraction of processes physiologically relevant to drought response. The drought-driven changes in gene expression occurring in root apices were consistent across treatments and genotypes. For mature leaves, the transcriptome response varied weakly but in accordance with the duration of water deficit. A differential clustering algorithm revealed similar and divergent gene co-expression patterns among the two genotypes. Since moderate stress levels induced similar physiological responses in both genotypes, the genotype-dependent transcriptional responses could be considered as intrinsic divergences in genome functioning. Our meta-analysis detected several candidate genes and processes that are differentially regulated in root and leaf, potentially under developmental control, and preferentially involved in early and long-term responses to drought.
In poplar, the well-known drought-induced activation of sensing and signalling cascades was specific to the early response in leaves but was found to be general in root apices. Comparing our results to what is known in arabidopsis, we found that transcriptional remodelling included signalling and a response to energy deficit in roots in parallel with transcriptional indices of hampered assimilation in leaves, particularly in the drought-sensitive poplar genotype.
Water deficit is recognised as one of the main environmental constraints restricting natural and agro-ecosystem productivity [1, 2]. The influence of water availability on plant productivity suggests that water limitation has shaped the natural variation and evolution of many physiological traits . Biotechnology has investigated the genetic basis of drought tolerance by targeting relevant genes [4, 5]. However, manipulating a single gene at a time, even genes encoding transcription factors, has proved insufficient to maintain productivity under drought; the tendency of cell systems is to restore homeostasis and mutants showing improved drought-or salt-tolerance are often stunted [6, 7]. Transcriptome meta-analyses and regulon detection could contribute to identifying physiological processes relevant to drought response. The availability of genomic tools such as high density arrays or whole-genome microarrays has led to an increasing number of studies examining drought-induced transcriptional remodelling [8–15].
Drought tolerance, which is obviously a multigenic trait, clearly depends on genome x environment interactions . Comparative genomics has emerged as a promising means of unravelling the molecular networks underlying complex traits [17, 18]. Studies analysing gene expression in parallel with quantitative trait loci (QTL) have shown that drought-induced changes in gene expression in specific genotypes were consistent with the observed physiological responses [11, 19]. Wilkins et al. highlighted the genotype specificity of the transcriptome response occurring in poplar leaves under drought conditions . Two drought-tolerant genotypes of maize showed more rapid and drastic changes in gene expression under drought, especially during recovery, than a drought-sensitive genotype . Moreover, drought tolerance relies on physiological adjustments occurring in distinct organs, involves the interplay of signalling/sensing cascades, and supposes integrated responses from molecular to whole-plant level. The few studies that have compared leaf and root transcriptomes have highlighted the organ specificity of drought responses [9, 20, 21]. Roots sense the edaphic water deficit, send chemical signals to shoots, and maintenance of root growth despite reduced water availability can contribute to drought tolerance through water foraging . Species-dependent features also shape the transcriptome response; almost none of the 27 genes reliably responsive to water stress in arabidopsis were regulated under drought in poplar and pine [7, 23]. Besides, many genes "inducing drought tolerance" have been identified under an abrupt and/or severe stress, which is a far cry from a realistic, slow developing and long-lasting drought [8, 24]. By quantifying and controlling water deficit by physical variables such as soil relative water content, the two genotypes experienced a similar degree of water deficit, allowing a robust comparison of their physiological and molecular responses . Such an ecophysiological approach has proved to be an efficient means of comparing drought response across genotypes [26, 27].
Since the publication of the Populus trichocarpa genome in 2006, poplar has become the model species for trees and the most studied tree species overall . In addition, poplars are also ecologically and economically important. Poplar inter-specific hybrids (Populus spp.) are among the fastest-growing trees under temperate latitudes and are grown for pulp, paper wood and fuel production purposes . While poplars are known to be sensitive to water deprivation as compared to other trees, drought tolerance varies considerably between genotypes of Populus, both inter-and intra-specifically [30, 31]. The purpose of the present study was to analyse the transcriptome responses induced by mild-to-moderate water deficit in two poplar genotypes. On the one hand, we applied a short-term water deficit to access sensing and signalling events; on the other, we imposed prolonged water deficit to reveal the molecular controls of plant performance under steady state stress. Transcriptome responses were analysed in mature leaves and in growing root apices in order to gain a wide assessment of the response and to draw an integrative picture of molecular responses to drought. The comparison of two genotypes known to differ in their drought tolerance revealed not only reliable drought markers but also the divergences and similarities in transcriptional networks, highlighting candidate genes for future diversity screening.
We focused on two hybrid genotypes ('Carpaccio' and 'Soligo') exhibiting contrasting tolerance to drought in field experiments, Carpaccio productivity being less hampered by drought than that of Soligo . Young trees were submitted either to a short-term water deficit by withholding irrigation 36 hours before harvest [early response (EAR)] or to a 10 day-long treatment, inducing a mild drought [long-term response to mild stress (LMI)] or a moderate drought [long-term response to moderate stress (LMO)].
Genomic DNA-based selection of probesets
Number of probesets
Present only on genomic DNA arrays
Present only on Carpaccio arrays
Present only on genomic DNA arrays
Expressed in root and leaf
Present only on Soligo arrays
Present only on genomic DNA arrays
Expressed in root and leaf
2-Probesets matching on both genotype arrays
Both genotypes expression arrays
Expressed in root and leaf
Both genotypes genomic DNA arrays
Only on Carpaccio expression arrays
Only on Soligo expression arrays
Transcriptome response in root apices is consistent across treatments and genotypes but contrasted in leaves
Correlation between changes in gene expression in the twelve conditions
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
EAR vs CTL
LMI vs CTL
LMO vs CTL
Functional categories involved in drought-induced transcriptome responses
Drought-responsive probesets were assigned to poplar gene models and annotated using web-based queries. Screening allowed the annotation of almost all probesets. Only 89 out of the 6,725 probesets were neither assigned to a poplar gene model nor functionally annotated. For the sake of brevity, and given the assumption that best homology implies a true ortholog, gene names of the closest arabidopsis genes are used to describe poplar genes. An overview of functional groups involved in drought responses was obtained through singular enrichment analysis (SEA, Additional file 2). Concerning the Cellular Component ontology, "endomembrane system" (GO:0012505) and "cell wall" (GO:0005618) were sorted in response to all drought conditions and were highly enriched for Soligo roots (EAR down-regulation). Another term, "nucleus" (GO:0005634), was enriched in the up-regulated responses to EAR for both organs and to all drought conditions for Soligo roots, whereas it was found among the down-regulated responses for LMO-treated Carpaccio roots. The most enriched terms of Biological Processes ontology were "response to stimulus" (GO:0006950), "response to abiotic stimulus" (GO:0009628), "response to stress" (GO:0006950) and "response to endogenous stimulus" (GO:0009719). These enrichments were seen for both organs under all conditions and indicated not only that the constraint was perceived--although the treatments were only mild-to-moderate--but also that prolonged drought remained a stressful growing condition. Interestingly, "growth" (GO:0040007) was specifically enriched for Soligo roots (EAR down-regulation). For mature leaves, metabolisms responded differentially in the two genotypes. "Primary metabolic process" (GO:0044238) was enriched for Carpaccio (LMI up-regulation). In contrast, for Soligo, "photosynthesis" (GO:0015979) and "secondary metabolic process" (GO:0019748) were detected among the down-regulations occurring in response to short and prolonged drought, respectively.
In terms of Molecular Function ontology, the most significant enriched term was "binding" (GO:0005488). Consistent with this, "transcription factor activity" (GO:0003700) and "transcription regulation activity" (GO:0030528) were sorted under all treatments (up and down regulations). "Transporter activity" (GO:000521) was enriched for Soligo leaves (LMI up and down regulations) and for Carpaccio roots (LMI and LMO up-regulations). Other GO terms, such as "catalytic activity" (GO:0003824), "hydrolase activity" (GO:0016787), and "transferase activity" (GO:0016740) were found in most lists.
Drought-dependent enrichment of functional groups in leaf arrays
CL-EAR vs CL-CTL
Major intrinsic protein
Unknown protein/ABI5 binding
Bet v I allergen/PYL
Carotene dioxygenase activity/NCED
CL-LMI vs CL-CTL
Metal ion binding/SOD, metallothionein
ATP dependant helicase/DEAD-box
CL-LMO vs CL-CTL
Nutrient reservoir activity/Germin, Extensin-like
Protein amino-acid phosphorylation/Protein kinase
Flavonoid 3'-monooxygenase activity
Calcium ion binding/EF-hand
Cell redox homeostasis/Glutaredoxin
Metal ion binding/SOD, metallothionein
Zinc ion binding
SL-EAR vs SL-CTL
Calcium ion binding/EF-hand
Regulation of transcription, DNA-dependent
Protein amino acid phosphorylation
Ankyrin repeat family protein
ATP synthesis coupled proton transport
SL-LMI vs SL-CTL
Cysteine-type peptidase activity/Papain
Drug transporter activity/MatE
SL-LMO vs SL-CTL
Protein ser/thr phosphatase activity/PP2C
O-glucosyl hydrolase activity/β-glucosidase
ABA metabolic process/Xerico
Amino acid transport
No apical meristem (NAM) protein
Two-component signal transduction
In roots (Additional file 3), many functional groups were detected in accordance with the strong responsiveness to drought of the root transcriptome. In these large lists, the iGA procedure revealed a striking conservation of the concerted changes in most conditions, highlighting a common response to drought. Whatever the condition, we detected enrichment in groups of genes known to be responsive to abiotic stress and/or drought. The generic response involved genes that were related to i) ABA biosynthesis/signalling (up-regulated: "NCED", "PP2C"; down-regulated: "PYL"), ii) cell rescue and/or cell redox homeostasis (up-regulated: "dnaJ", "heat shock protein", "glutathion-S-transferase", "metallothionein"), and iii) the response to hypoxia (down-regulated: "alcohol dehydrogenase"; "pyruvate decarboxylase", "LOB domain-containing protein"). As expected for actively growing organs, stress impacted recurrent groups of genes that were involved either in expansion (up-regulated: "aquaporins"; down-regulated: "pectinesterase", "L-ascorbate oxidase") or in meristematic activity and cell cycle (down-regulated: "chromosome organization", "DNA replication"). EAR-treated Soligo roots underwent an extensive metabolism-related response (up-regulated: "raffinose synthase", "asparagine synthase", "trehalose phosphatase"; down-regulated: "nucleotide-sugar metabolism", "fatty acid desaturation"). In addition, an erosion of the "transcription factors" group was detected across treatments. In Carpaccio roots, a "transcription factors" group of nine up-regulated genes was detected in response to EAR (putative "Homeobox proteins" (ATHB12, ATHB40, ATHB6), "Responsive to desiccation" (RD26), "Nuclear transcription factor Y, alpha" (NF-YA), "Heat shock transcription factor C1" (AtHSFC1), "ABA-Repressor1" (ABR1)), whereas five up-regulated genes were detected under LMI (putative "ATHB12", "ATHB40", "WRKY"), and only one under LMO. Concerning Soligo, we detected five "transcription factors" groups in response to EAR (38 genes including putative "ATHB12", "RD26", "ABR1", "AtHSFC1"), two groups under LMI (3 genes including putative "ATHB12", "RD26") and no enrichment under LMO. As highlighted by this functional categorization, the applied treatments clearly drove transcriptome responses in both organs and genotypes. Both procedures provided consistent results, enabling the extraction of processes physiologically relevant to drought responses. Several genes of interest can be discriminated on the basis of their contribution to enriched functional categories.
Analysis of drought-responsive gene networks based on gene co-expression relationships provides robust drought markers and candidate genes
In a joint analysis of drought-driven regulation in leaves and roots, 16 clusters of co-expressed genes were first identified in Carpaccio (Figure 3a). The DCA procedure revealed 10 transcriptional modules as fully conserved between the two genotypes (clusters 1, 3, 4, 5, 7, 8, 9, 10, 12 and 15). Delineation of these transcriptional groups confirmed the contrasting gene regulation in mature leaves and growing root apices. Four clusters (13-16) appeared to be relatively distant from each other compared to two homogeneous sets of clusters (clusters 1-8, 1341 up-regulated genes; clusters 9-12, 565 down-regulated genes). In between, some co-regulation relationships appeared to be partially conserved (clusters 6 and cluster 11), meaning that a subset of the genes that are co-expressed in the reference genotype lost their co-expression relationships in the target genotype. This phenomenon applied to 39 genes (with transporters of different substrates included in cluster 6, subset b, see Additional file 4) and 28 genes (with transcription factors and hormone biosynthetic enzymes, such as putative Allene oxide synthase included in cluster 11, subset b, see Additional file 4). Finally, the large divergence of co-expression relationships between genes in Carpaccio and Soligo were highlighted (clusters 2, 13, 14 and 16, totalling about 600 genes). Up-regulated cluster 2 included several Aquaporins and genes related to ABA signalling or the response to oxidative stress (Additional file 4). Loss of co-expression relationships between genes was due mainly to almost invariant gene expression in Carpaccio (cluster 16) or in Soligo (clusters 13 and 14). These two latter transcriptional modules collected genes with a similar annotation (Unknown proteins, putative Cytochrome P450, Leucin-rich repeat proteins) or associated with "cell wall" or "transport activity" (putative Wall-associated kinases, Amino-acid permeases).
In mature leaves, drought-driven regulation of gene expression was weak (Figure 2) and led to the clear distinction of 12 transcriptional modules (Figure 3b). Drought-driven transcriptome responses appeared strongly dependent on stress duration, splitting clearly into early and long-term responses in accordance with global gene expression analyses (Table 2). Most transcriptional modules of co-expressed genes were conserved. On the one hand, some conserved modules were regulated exclusively in response to prolonged (clusters 1, 6 and 7) or short (cluster 9) drought. On the other hand, conserved clusters 8 and 11 encompassed 81 genes that were responsive to all drought conditions but that were inversely regulated in short and in prolonged drought. The contrast between short and prolonged drought responses in mature leaves relied on the regulation of 301 distinct genes (Additional file 4). The DCA procedure also revealed drought marker genes that were up-regulated in both genotypes and under all drought conditions (conserved cluster 2; including putative RCI2A, ATHB12, Galactinol synthase AtGolS2, ABC transporter). The sole transcriptional module identified as not-conserved between the two genotypes (cluster 12) resulted partly from the genotype-specific response to LMI (which repressed a lower number of genes in Soligo than in Carpaccio, Figure 2a). Far more informative were the partially conserved clusters (clusters 3, 5, and clusters 4, 10, Additional file 4). Genes were almost invariant in Soligo but strongly regulated in Carpaccio in sub-group 3b (24 up-regulated genes, e.g. putative Aquaporin, Flavonol synthase, five transcription factors) and in sub-group 5b (nine down-regulated genes including three putative Wall-associated kinases). In cluster 4b, genes were strongly repressed by LMI in Carpaccio leaves but were not drought-responsive in Soligo (15 genes including putative β-xylosidase1, Pectinesterase, 5 genes related to defence response). Conversely, genes that were almost invariant in Carpaccio and drought-responsive in Soligo were split into up-regulated sub-group 10a (including hormone signalling of jasmonic acid, auxin and ABA) and down-regulated sub-group 10b (including putative ATP sulfurylase, Squalene epoxidase, Gibberellin-2-β-dioxygenase2).
In root apices, the strong intensity of drought-driven expression patterns allowed 12 transcriptional modules to be defined unambiguously (Figure 3c, left panel). Two distant sets of clusters could be discriminated (down-regulated clusters 1-4, totalling 612 genes, and up-regulated clusters 5-12, totalling about 1300 genes). The down-regulated set included genes related to cell cycle and DNA processes, notably putative cyclin-dependent protein kinases and DEAD-box RNA helicases. The up-regulated groups included genes related to metabolic processes (putative Phosphoenolpyruvate carboxylase; Dihydrodipicolinate reductase1, Alternative oxidase), and to catabolism (putative Ubiquitin-protein ligases). As shown in Figure 3c, Carpaccio transcriptional modules were delineated according to differences in intensity of gene regulation across treatments. This clustering outline was valid for both genotypes, which is consistent with the high correlation observed between treatments and genotypes (Table 2). However, most of the transcriptional modules were labelled "not conserved". In these cases, the homogeneous response, in terms of up or down regulation, masked the fine genotype-specific tuning of transcriptome responses that was revealed by the DCA procedure. In cluster 8 only, the expression pattern differed widely between genotypes. Besides, some modules were labelled "fully" or "partially" conserved, indicating that genes conserved their co-expression properties (clusters 5 and 7, sub-groups 6a, 9a and 11a, totalling about 472 genes, Additional file 4). In root apices, drought-driven regulation was highly consistent in both genotypes, and it was the differential tuning across drought conditions that accounted specifically for the loss of conservation in gene co-expression relationships.
Common and specific components of drought-induced response in different organs
Our experimental design did not allow unequivocal testing of organ-preferred expression since tissue specificity (root versus leaf) and maturity level (growing versus mature) were confounded. To our knowledge, only one study has compared leaf and root transcriptomes under drought, thus precluding a robust meta-analysis of the differential impact of drought according to the organ. Given that gene expression patterns between leaf and root have been compared more extensively under optimal conditions, we delineated two subsets of drought-regulated genes on the basis of consistent expression pattern in control and drought treatment (Figure 4B). Among drought-regulated genes with a leaf-preferred expression pattern, one-fourth was confirmed to be expressed preferentially in leaf in other Populus species. A similar result was found for genes expressed preferentially in root. The expression patterns of 775 genes were not confirmed in other Populus species and the patterns of 475 genes could not be tested due to missing information. To our knowledge, the interaction of organ maturity and drought on the transcriptome responses has not been documented in poplar. Although maturity level could affect the sensitivity of the transcriptome response to drought, we searched for expression profiles in young and mature organs under optimal conditions in another Populus species (exPlot) . From the two subsets of drought-regulated genes in our study, we extracted six contrasting expression profiles under optimal conditions (Figure 4C). S1 gathers potential markers of growing tissues and could also be drought-regulated in growing leaves. In contrast, genes in S2 and S3 were confirmed as being expressed preferentially in root, either in old roots (S2) or in young roots (S3). Genes that are expressed preferentially in young leaves under optimal conditions were drought-regulated in growing roots or mature leaves (S4). Among genes expressed preferentially in mature leaves in our experiment, 20 were confirmed to be expressed preferentially in mature tissue, either in leaf (S5) or in both leaf and root (S6). Although indirect (i.e. extrapolated from other Populus species grown under optimal conditions), these arguments strengthen the hypothesis of an organ-preferred expression of some drought-regulated genes. However, missing information precludes a robust meta-analysis allowing the drought responses within the Populus genus to be unravelled. Given that most drought transcriptome studies have focused on mature leaves, the comparative information provided here will benefit future integrative approaches.
Maintenance of water status, a key process in plant functioning, is actively regulated on the whole plant scale (from root uptake to stomata) in response to variations in water availability. Given the known difference in drought sensitivity of the two genotypes, severe drought could yield two contrasting physiological states . In our study, the physiological and transcriptional responses clearly indicate that both genotypes perceived water deficit as stressful. However, leaf water status was maintained and growth similarly hampered in both genotypes. The moderate stress levels applied induced similar physiological responses in both genotypes, allowing genotype-dependent transcriptional responses to be considered as intrinsic divergences in genome functioning rather than the result of the interaction between genome and physiological status. Drought sensing and metabolic adjustments involve tight molecular control . We examined this control in poplar, analysing transcriptional remodelling in response to short and prolonged water deficits and in parallel in root apices and mature leaves. Given that poplar pathway information was inferred mainly from the arabidopsis literature, common gene names of the closest arabidopsis homolog were used to describe poplar genes (see Additional file 1 for correspondence) .
Controlling energy and drought signalling under drought: a candidate process related to productivity
ABA biosynthesis and signalling under drought: consistent transcriptome responses in leaf and root
The up-regulation of 9-cis-epoxycarotenoid dioxygenase (NCED), Phytoene synthase and β-carotene hydroxylase (BETA-OHASE1) suggests enhanced ABA biosynthesis in roots under all drought conditions (Additional file 7). Consistently, this regulation has been previously associated with a peak in ABA content during poplar bud development . BETA-OHASE could be involved in ABA biosynthesis and its over-expression enhanced stress tolerance in arabidopsis [43, 44]. In leaves, transcript profiling suggested that ABA biosynthesis was no longer activated under prolonged drought--only NCED was up-regulated under short drought stress--although ABA signalling was still detected, suggesting long-distance transport of the phytohormone.
Concerning ABA-responsive genes, Pyarabactin resistance-like (PYL) were down-regulated and those of Protein phosphatase type-2C (PP2C) were up-regulated. PYLs encode ABA receptors that interact with PP2C as regulators of the ABA-mediated signalling pathway . The opposite regulation of PYL and PP2C, either in roots or in both organs, was described in arabidopsis in response to ABA treatment . As shown in Additional file 7, we detected several orthologs of Dehydration-responsive-element binding protein (DREB2A, 1A, 1D up regulated in roots and DREB3 repressed in leaves) and four predicted ABA-RESPONSIVE-ELEMENT binding factors known to be involved in signalling cascades . Interestingly, a putative transcription factor, close to Homeobox protein ATHB7 and ATHB12, was up-regulated under all drought conditions. These ABA-induced growth mediators were up-regulated in response to water-deficit in arabidopsis . Among up-regulated transcription factors were two putative Responsive to desiccation 26-- a transcriptional activator in ABA signal transduction . Eight orthologs of Arabidopsis nuclear factor-Y (NF-YA or NF-YB) were drought-regulated in poplar . In arabidopsis, drought-driven induction of NF-YA5 and of NF-YB1 was ABA-dependent or ABA-independent, respectively [51, 52]. Here, putative NF-YA and NF-YB were both up-regulated in Carpaccio roots. However, only putative NF-YA transcripts accumulated in EAR-treated leaves, which is in line with ABA-dependent activation. The ZFP family, encoding C2H2-type zinc finger proteins, was found to be drought-responsive in poplar. Three genes similar to Salt-tolerance zinc fingers and one similar to Arabidopsis zinc finger protein 2 were regulated by drought in poplar, in accordance with their responsiveness to ABA and abiotic stress in arabidopsis . In addition, three drought-responsive genes in poplar were similar to ZAT12, which is involved in ROS and abiotic stress signalling in arabidopsis . Among them, ZFP2 was up-regulated in roots under all drought conditions, in line with its induction by abiotic stress in poplar . Three putative RARE-COLD-INDUCIBLE (RCI2) were up-regulated in poplar. In arabidopsis, most RCI2 genes are induced by ABA and abiotic stress, and are potentially involved in the regulation of plasma membrane potential . RCI2-A contributes to salt tolerance by preventing over-accumulation of K+ and Na+. The stress-responsive plasma membrane protein COR413-PM is potentially involved in signal transduction . Two putative COR413-PM were induced in roots in response to drought. As highlighted in Figure 4, drought repressed several genes in poplar that are known to be up-regulated in response to hypoxia in arabidopsis and in poplar [59–61]. Such down-regulation has been reported previously in ABA and/or drought responses and could reflect antagonism between ABA and ethylene signalling [62–65]. Accordingly, poplar response to drought implied cross-talk between hormonal pathways (hormone metabolism or/and signalling, Figure 4, Additional file 1) in accordance with the literature [66, 67].
Multi-stress responsive genes in poplar
WRKY transcription factors constitute a large family of plant-specific regulators controlling senescence and responses to stress and ABA [66, 68–70]. In poplar, we detected 31 drought-responsive WRKY. In arabidopsis, AtWRKY-53, -54 and -70 were found to be structurally related, AtWRKY-18 and -40 could interact each other, AtWRKY-53 and -70 exhibited partial overlapping functions, and AtWRKY-70 and -54 counteracted accumulation of salicylic acid . In silico analysis of poplar WRKY-40 and -53 (eugene3.00061944 and grail3.0007034202) revealed an EAR-motif--a potential signature of transcriptional repressors . Putative AtWRKY-53, -40, -18, -54 and -70 exhibited similar expression patterns, being up-regulated in roots under prolonged drought and in leaves in the early response, and/or being repressed in leaves under prolonged drought. Conserved clusters 5 and 6a (Figure 3a) gathered nine WRKY and some co-regulated genes such as NIM1-interacting1 and a putative Pathogenesis-related protein1. In poplar, WRKY and a putative Plant natriuretic peptide (AtPNP-A) were co-regulated by drought as already observed in arabidopsis . AtPNP-A, an extracellular signalling molecule, could affect water and solute transport in response to stress .
Reactive oxygen species (ROS) have been associated with stress sensing/signalling, and have emerged as important general signals . Although the drought applied here was mild-to-moderate in degree, we highlighted the activation of oxidative detoxification processes. In maize under mild water deficit, ABA accumulation triggered the generation of ROS, which up-regulated the antioxidant system . Glutathion-S-transferase, Superoxide dismutase, Lactoylglutathione lyase, Catalase, Glutathione peroxidase, and Ascorbate peroxidase were up-regulated in roots under all conditions, and in leaves more strongly under prolonged stress. Drought also induced the up-regulation of several genes involved in oxidative stress tolerance/response in arabidopsis (Additional file 7), such as Alternative oxidase (AOX prevents mitochondrial ROS formation), Temperature-induced lipocalin, AT1G68440, Senescence-associated protein1 or Senescence-associated gene21[77–81]. In poplar roots, drought positively regulated two putative heat shock-like protein (HSPRO2), which, in arabidopsis, were involved in tolerance to oxidative stress and were drought-and ABA-responsive [23, 82]. The ChaC-like family protein, which was repressed in response to oxidative stress, were strongly down-regulated in poplar leaves . In accordance with ROS production and detoxification processes, we detected up-regulation of four raffinose synthases in roots, of a stachyose synthase in leaves and of three galactinol synthases (two in roots and one in leaves). An increase in galactinol, raffinose and stachyose content could improve osmoprotection and ROS scavenging . The poplar drought response also implied the induction of genes related to cell rescue, including detoxification and chaperone-like activities (peptidyl-prolyl-cis-trans-isomerase, heat shock protein, dnaJ).
In poplar, predicted lipid transfer proteins (LTP) and remorin (REM) were induced strongly in roots, and late embryogenesis proteins (LEA) were up-regulated in both organs. Similar inductions in response to drought were detected in arabidopsis . Non-specific LTPs could be calmodulin-binding, implying a possible Ca2+/CaM signalling function . AtREM4.1 and AtREM4.2 were induced strongly in response to osmotic, salt, drought, ABA and brassinosteroid treatments . Stress responsive LEAs are suspected to act as chaperones and/or ROS scavengers, to bind metal ions or divalent cations . Among the four Heat-stable protein1 (HSP1) up-regulated in roots, eugene3.00101442 was specifically regulated in Carpaccio. In Populus tremula, SP1 was identified as a new stress responsive protein . In poplar, we detected several genes that are commonly down-regulated by drought in arabidopsis . Among them, two putative Germin-like (AtGER1 and AtGER3) were repressed under short water deficit but induced in response to LMO. AtGER1 and AtGER3 may be involved in the control of synthesis of cell wall polysaccharides and/or in scavenging of extracellular nucleotide-sugars [89, 90]. In poplar, the two ESKIMO orthologs, namely gw1.VIII.1375.1 and gw1.X.1696.1, were repressed under most drought conditions. In arabidopsis, ESKIMO1, a positive regulator of transcription, negatively regulates cold acclimation . The eskimo1 mutant was not drought-or salt-tolerant although ATHB7, ATHB12 or PP2C were constitutively up-regulated . However, Eskimo1 is suspected to play a role in whole-plant water economy in arabidopsis . While poplar gene functions are extrapolated mainly from sequence similarities, the transcriptome analysis performed here has given new insights into their involvement in cell physiology.
Comparative genomics is a powerful tool that can help decipher the molecular basis of drought responses and reveal physiologically relevant processes. Reliable stress markers were extracted as well as genes whose expression differed in tolerant and sensitive genotypes. Similarly, analysis of variance detected a strong genotype effect in the transcriptome responses of poplar leaves to drought . However, when using multi-species designed arrays, the risk of misinterpreting divergent signals has to be acknowledged and controlled for. We used a genomic DNA-based selection strategy to improve the detection of differentially expressed transcripts. Hybridizing genomic DNA was previously used in genotyping arabidopsis accessions and for analysing transcriptomes by cross-hybridization (banana/rice, chimpanzee/human) [93–95]. Our quantitative analysis of gene expression in poplar provides an unambiguous comparison of two hybrid transcriptomes. Although requiring to be tested further (wider range of organs, other poplar species, field conditions) our meta-analysis has revealed several candidate genes and processes that are differentially regulated in root and leaf, potentially under developmental control, and preferentially involved in rapid and long-term response to drought. Since most of these genes were not previously ascribed to poplar drought response, our work provides expression data that will enrich our knowledge of gene function in Populus.
Cuttings of two Populus deltoides W. Bartram ex Marshall x Populus nigra L., namely 'Carpaccio' and 'Soligo', were planted in 2L-pots filled with a peat-sand mix (50/50 V/V) amended with 1 g L-1 of CaMg(CO3)2 and 4 g L-1 of fertilizer (Nutricote 711, Fertil; http://www.fertil.fr/) and were grown for two months in a greenhouse. In order to favour the development of a dense root system, the initial stem was cut at a few centimetres above its base. This "detopped cutting" was then transplanted into a 10L-pot filled with the same substrate. A new stem was allowed to grow for 10 weeks without water limitation. Trees were assigned randomly to 4 modalities of water supply. For controls (CTL), evaporative demand was compensated by 4 to 6 waterings to field capacity per day. For short-term water deficit (EAR), irrigation was withheld for 36 hours prior to harvest bringing soil relative extractable water (REW) into the range of 20-35%. For long-term drought, soil REW was controlled by water supply 4 times a day as detailed in . Soil REW was maintained for 10 days either at 20-35% (LMI) or at 10-20% (LMO). Ambient conditions depended on outside weather but temperature was maintained in the range of 19-26°C, humidity varied between 55 and 85% (day/night) and PAR between 400 and 950 μmol m-2 s-1 (cloudy versus sunny days).
For each genotype x treatment combination, six trees were assigned to ecophysiological monitoring. Leaf predawn water potential, leaf relative water content, leaf full turgor osmotic pressure, height growth rate, and gas exchange were measured as previously described . Instantaneous water use efficiency (WUEi) was calculated as the ratio of net CO2 assimilation rate to stomatal conductance for water vapour. Six other trees were devoted to molecular analyses. Controls and treated plants were harvested simultaneously following random sampling, between 11:00 am and 3:00 pm. Mature leaves and root apices of each tree were harvested in parallel. Two mature leaves were cut and immediately frozen in liquid nitrogen. In less than 30 seconds, about ten 1-cm-long apices were sampled in the whole root system, frozen immediately in liquid nitrogen and stored at -80°C.
DNA and RNA extraction, RNA amplification and array hybridization
Total genomic DNA from leaves was extracted, fragmented, labelled and hybridized as described in  with the following modifications. Fifty μg of gDNA were partially digested with DNAse1 (Promega, http://www.promega.com/). DNAse1 was heat-inactivated with 2 μL inactivation buffer (10 min at 65°C). gDNA fragments were labelled by adding 200 U terminal deoxynucleotidyl transferase (90 min at 37°C) and hybridized for 20 h at 45°C.
Total RNA of each sample was extracted separately. Total RNA was extracted from 100 mg of leaves and 30 mg of roots with an Rneasy Plant Mini kit, using a DNAse1 treatment (Qiagen, http://www1.qiagen.com/). After checking integrity (2100 Bioanalyzer, Agilent, http://www.home.agilent.com/), RNA was quantified (RiboGreen RNA Quantification Reagent, http://www.promega.com/). Amplification and hybridization on Affymetrix GeneChip Poplar Genome Arrays were performed according to the manufacturer's protocol (Affymetrix, http://www.affymetrix.com/). Arrays were scanned with the GeneChip Scanner 3000-7G piloted by the GeneChip Operating Software (GCOS).
Transcriptome analysis was conducted using 36 Affymetrix GeneChip Poplar Genome Arrays. Four arrays were devoted to genomic DNA hybridization (two technical replicates per genotypes). For expression arrays and for all conditions (2 genotypes × 4 treatments × 2 organs), the six trees were assigned randomly to two biological replicates. In each biological replicate, the same individuals were pooled for both organs and the three individuals contributed equally to the pool of total RNA. All raw and normalized data are available through both the CATdb database [AFFY_POPSEC_Nancy_Roots_poplar, AFFY_POPSEC_Nancy_Leaves_poplar and AFFY_genomic_Poplar, http://urgv.evry.inra.fr/CATdb] and the Gene Expression Omnibus repository at the NCBI [GSE17223, GSE17226, GSE17230 and GSE21334; http://www.ncbi.nlm.nih.gov/geo/]. The 16 arrays devoted to root samples, 16 arrays devoted to leaf samples and 4 arrays hybridized with genomic DNA were normalized separately with the gcrma algorithm available in the Bioconductor package [96, 97]. To determine which genes were differentially expressed between two given conditions, we performed a two group t-test assuming equal variance between groups. To fit the assumption of equal variance of gene expression per group, genes displaying extreme variation (too small or too large) were excluded from the analysis. The raw P- values were adjusted by the Bonferroni method, which controls the Family Wise Error Rate . A gene was declared differentially expressed if the Bonferroni P- value was below 0.05. SR-EAR vs SR-CTL comparison had a lower variance than the others (0.037 versus a mean variance of 0.052 ± 0.001). Seven genes were selected for RT-qPCR validation (Additional file 8). For each condition, four out of the six samples were chosen randomly for the RT-qPCR procedure (as described in ; with 500 ng RNA). Similar gene expression patterns were obtained with both methods (Additional file 9).
In the probeset selection procedure, the background level was set to 3.1 on the basis of the mean signal intensities of 62 reporter probesets (i.e. several controls that are not in the investigated poplar genomes; mean value 2.5 and mean maximum value 3.5). For a given array, any probeset with a signal intensity below this cut-off value was labelled "absent". In the expression arrays, when a probeset was labelled "present" in the two biological replicates of a condition, the targeted transcript was considered expressed. For genomic DNA, the analysis was less stringent. We considered that hybridization was possible when the probeset was labelled "present" in at least one of two technical replicates.
GO enrichment and detection of differentially expressed gene groups
Probesets were assigned to gene model (poplar genome v1.1, http://genome.jgi-psf.org/poplar/poplar.home.html) using the batch query of NetAffx Analysis Center http://www.affymetrix.com/analysis/index.affx, poplar database query at JGI and similarity researches at NCBI http://www.ncbi.nlm.nih.gov. Sequence alignments were performed using BLASTn (default parameters and a maximal e- value of 10-5) and only the best homologies were considered further. Following the release of phytozome v5.0 in January 2010, the best homology was confirmed using the annotation of the poplar genome v2.0 http://www.phytozome.net/poplar. We employed singular enrichment analysis (SEA, http://bioinfo.cau.edu.cn/agriGO) on up-regulated or down-regulated gene lists sorted for each condition . SEA using plant GO slim was performed independently for each condition with the Populus Affymetrix Genome array as a background list, followed by correction for multiple testing. Functional annotation and Gene Ontology were retrieved by querying the annotation browser in agriGO with Affymetrix probesets, by querying the annotation batch function at PopGenie http://188.8.131.52/popgenie1 with the poplar gene model, and by querying the annotation tool at The Arabidopsis Information Resource (TAIR) http://www.arabidopsis.org [37, 101]. Annotations from all origins were compiled to determine functional class enrichments using Iterative Group Analysis (iGA) . The iGA procedure, which is based on hypergeometric statistic calculations, detects concerted changes in functional classes and assigns a probability of change (PC-values) to each functional class. For each condition, differentially expressed genes were sorted by their mean normalized expression ratio in ascending or descending order. The iGA procedure was applied separately for up-and down regulation to determine which functional groups are most enriched at the top of the sorted gene lists .
Differential clustering algorithm
The DCA, first described by Ihmels et al., was performed using an R script http://www.R-project.org developed by Lelandais et al. [35, 36]. Transcriptional modules in the reference genotype were detected using a hierarchical clustering algorithm (with 'hclust' function and with the 'ward' method for probeset agglomeration) and, for each module, the corresponding gene expression patterns in the target genotype were segmented into two different sub-clusters (labelled as "a" and "b") using the same hierarchical clustering algorithm. The DCA results are presented as a distance matrix between gene expression measurements (reference genotype in rows and target genotype in columns). Transcriptional modules of co-expressed genes were first defined in Carpaccio (reference genotype), and their corresponding probesets in Soligo (target genotype) were next clustered into two groups according to their expression measurements across Soligo arrays. Clusters were automatically assigned to four categories ("full", "partial", "split" or "no" conservation), calculating the mean correlation of probesets within and between sub-groups "a" and "b", and comparing the values obtained to a specific threshold T (T = 0.4 in this study).
A set of unique gene models was delineated from the list of drought-responsive probesets. When gene annotation was multiple (cross-hybridization), data were discarded. Each gene was described by Affymetrix probeset identifiers, the Populus genome v1.1 gene name and the AGI code. For genes matching multiple probesets, the expression pattern was assessed by a "present call" in at least one probeset. Using Venn diagrams, our list was compared with previous poplar studies [11, 14, 17, 18, 21, 37, 61, 103–106]. Impact of maturity level on gene expression under optimal conditions was assessed using exPlot and data referred as UMA-0030 in UPSC-BASE and detailed in [17, 37, 106, 107].
This work was funded and DC was supported by the project ANR GENOPLANTE "POPSEC" (GPLA06028G). Part of this work was also funded by a Region Lorraine grant (# 12000160). We would like to thank P. Courtois, B. Clerc, C. Antoine, JM. Gioria, M. Granier, P. Gross, D. Himbert and B. Pêtre for their help during the experiment and the genomic platform of UMR1136 INRA-Nancy Université for hosting the molecular experiments. We thank Pr I. Couée for critical reading of the manuscript, P. Montpied for his advice, H. Rothnie for English editing and two anonymous reviewers for useful comments.
- Aussenac G: Interactions between forest stands and microclimate: Ecophysiological aspects and consequences for silviculture. Ann For Sci. 2000, 57: 287-301. 10.1051/forest:2000119.Google Scholar
- Passioura JB: Drought and drought tolerance. Plant Growth Regul. 1996, 20: 79-83. 10.1007/BF00024003.Google Scholar
- Dudley SA: Differing selection on plant physiological traits in response to environmental water availability: A test of adaptive hypotheses. Evolution. 1996, 50: 92-102. 10.2307/2410783.Google Scholar
- Salmeron J, Herrera-Estrella LR: Plant biotechnology-Fast-forward genomics for improved crop production. Curr Opin Plant Biol. 2006, 9: 177-179. 10.1016/j.pbi.2006.01.018.Google Scholar
- Century K, Reuber TL, Ratcliffe OJ: Regulating the regulators: The future prospects for transcription-factor-based agricultural biotechnology products. Plant Physiol. 2008, 147: 20-29. 10.1104/pp.108.117887.PubMed CentralPubMedGoogle Scholar
- Vinocur B, Altman A: Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol. 2005, 16: 123-132. 10.1016/j.copbio.2005.02.001.PubMedGoogle Scholar
- Polle A, Altman A, Jiang XN: Towards genetic engineering for drought tolerance in trees. Tree Transgenesis: Recent Developments. 2006, Berlin: Springer-Verlag, 275-297.Google Scholar
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, et al: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31: 279-292. 10.1046/j.1365-313X.2002.01359.x.PubMedGoogle Scholar
- Kreps JA, Wu YJ, Chang HS, Zhu T, Wang X, Harper JF: Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress. Plant Physiol. 2002, 130: 2129-2141. 10.1104/pp.008532.PubMed CentralPubMedGoogle Scholar
- Brosché M, Vinocur B, Alatalo ER, Lamminmaki A, Teichmann T, Ottow EA, Djilianov D, Afif D, Bogeat-Triboulot MB, Altman A, et al: Gene expression and metabolite profiling of Populus euphratica growing in the Negev desert. Genome Biol. 2005, 6: R101-10.1186/gb-2005-6-12-r101.PubMed CentralPubMedGoogle Scholar
- Street NR, Skogstrom O, 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. 10.1111/j.1365-313X.2006.02864.x.PubMedGoogle Scholar
- Wang HG, Zhang HL, Gao FH, Li JX, Li ZC: Comparison of gene expression between upland and lowland rice cultivars under water stress using cDNA microarray. Theor Appl Genet. 2007, 115: 1109-1126. 10.1007/s00122-007-0637-7.PubMedGoogle Scholar
- Hayano-Kanashiro C, Calderon-Vazquez C, Ibarra-Laclette E, Herrera-Estrella L, Simpson J: Analysis of Gene Expression and Physiological Responses in Three Mexican Maize Landraces under Drought Stress and Recovery Irrigation. PLoS One. 2009, 4: 10.1371/journal.pone.0007531.Google Scholar
- Wilkins O, Waldron L, Nahal H, Provart NJ, Campbell MM: Genotype and time of day shape the Populus drought response. Plant J. 2009, 60: 703-715. 10.1111/j.1365-313X.2009.03993.x.PubMedGoogle Scholar
- Berta M, Giovannelli A, Sebastiani F, Camussi A, Racchi ML: Transcriptome changes in the cambial region of poplar (Populus alba L.) in response to water deficit. Plant Biol. 2010, 12: 341-354. 10.1111/j.1438-8677.2009.00320.x.PubMedGoogle Scholar
- Tardieu F, Tuberosa R: Dissection and modelling of abiotic stress tolerance in plants. Curr Opin Plant Biol. 2010, 13: 206-212. 10.1016/j.pbi.2009.12.012.PubMedGoogle Scholar
- Street NR, Sjodin A, Bylesjo M, Gustafsson P, Trygg J, Jansson S: A cross-species transcriptomics approach to identify genes involved in leaf development. BMC Genomics. 2008, 9: 10.1186/1471-2164-9-589.Google Scholar
- Hamanishi ET, Raj S, Wilkins O, Thomas BR, Mansfield SD, Plant AL, Campbell MM: Intraspecific variation in the Populus balsamifera drought transcriptome. Plant Cell Environ. 2010, 33: 1742-1755. 10.1111/j.1365-3040.2010.02179.x.PubMedGoogle Scholar
- Hazen SP, Pathan MS, Sanchez A, Baxter I, Dunn M, Estes B, Chang H-S, Zhu T, Kreps JA, Nguyen HT: Expression profiling of rice segregating for drought tolerance QTLs using a rice genome array. Funct Integr Genomics. 2005, 5: 104-116. 10.1007/s10142-004-0126-x.PubMedGoogle Scholar
- Zhou JL, Wang XF, Jiao YL, Qin YH, Liu XG, He K, Chen C, Ma LG, Wang J, Xiong LZ, et al: Global genome expression analysis of rice in response to drought and high-salinity stresses in shoot, flag leaf, and panicle. Plant Mol Biol. 2007, 63: 591-608. 10.1007/s11103-006-9111-1.PubMed CentralPubMedGoogle Scholar
- Bogeat-Triboulot MB, Brosche M, Renaut J, Jouve L, Le Thiec D, Fayyaz P, Vinocur B, Witters E, Laukens K, Teichmann T, et al: 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. 10.1104/pp.106.088708.PubMed CentralPubMedGoogle Scholar
- Tardieu F, Simonneau T: Variability among species of stomatal control under fluctuating soil water status and evaporative demand: modelling isohydric and anisohydric behaviours. J Exp Bot. 1998, 49: 419-432. 10.1093/jexbot/49.suppl_1.419.Google Scholar
- Bray EA: Genes commonly regulated by water-deficit stress in Arabidopsis thaliana. J Exp Bot. 2004, 55: 2331-2341. 10.1093/jxb/erh270.PubMedGoogle Scholar
- Nanjo T, Futamura N, Nishiguchi M, Igasaki T, Shinozaki K, Shinohara K: Characterization of full-length enriched expressed sequence tags of stress-treated poplar leaves. Plant Cell Physiol. 2004, 45: 1738-1748. 10.1093/pcp/pci009.PubMedGoogle Scholar
- Granier C, Aguirrezabal L, Chenu K, Cookson SJ, Dauzat M, Hamard P, Thioux JJ, Rolland G, Bouchier-Combaud S, Lebaudy A, et al: PHENOPSIS, an automated platform for reproducible phenotyping of plant responses to soil water deficit in Arabidopsis thaliana permitted the identification of an accession with low sensitivity to soil water deficit. New Phytol. 2006, 169: 623-635. 10.1111/j.1469-8137.2005.01609.x.PubMedGoogle Scholar
- Aguirrezabal L, Bouchier-Combaud S, Radziejwoski A, Dauzat M, Cookson SJ, Granier C: Plasticity to soil water deficit in Arabidopsis thaliana: dissection of leaf development into underlying growth dynamic and cellular variables reveals invisible phenotypes. Plant Cell Environ. 2006, 29: 2216-2227. 10.1111/j.1365-3040.2006.01595.x.PubMedGoogle Scholar
- Reymond M, Muller B, Leonardi A, Charcosset A, Tardieu F: Combining quantitative trait loci analysis and an ecophysiological model to analyze the genetic variability of the responses of maize leaf growth to temperature and water deficit. Plant Physiol. 2003, 131: 664-675. 10.1104/pp.013839.PubMed CentralPubMedGoogle Scholar
- Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313: 1596-1604. 10.1126/science.1128691.PubMedGoogle Scholar
- Heilman PE, Hinckley TM, Roberts DA, Ceulemans R: Production Physiology. Biology of Populus and its implications for management and conservation. Edited by: Stettler RF, Bradshaw HDJ, Heilman PE, Hinckley TM. 1996, Ottawa: NRC Research Press, 459-490.Google Scholar
- Tschaplinski TJ, Tuskan GA, Gunderson CA: Water-stress tolerance of black and eastern cottonwood clones and four hybrid progeny. I. Growth, water relations, and gas exchange. Can J For Res. 1994, 24: 364-371. 10.1139/x94-049.Google Scholar
- Wullschleger SD, Weston DJ, Davis JM: Populus Responses to Edaphic and Climatic Cues: Emerging Evidence from Systems Biology Research. Crit Rev Plant Sci. 2009, 28: 368-374. 10.1080/07352680903241246.Google Scholar
- Monclus R, Dreyer E, Villar M, Delmotte FM, Delay D, Petit JM, Barbaroux C, Le Thiec D, Brechet C, Brignolas F: Impact of drought on productivity and water use efficiency in 29 genotypes of Populus deltoides x Populus nigra. New Phytol. 2006, 169: 765-777. 10.1111/j.1469-8137.2005.01630.x.PubMedGoogle Scholar
- Fossati T, Zapelli I, Bisoffi S, Micheletti A, Vietto L, Sala F, Castiglione S: Genetic relationships and clonal identity in a collection of commercially relevant poplar cultivars assessed by AFLP and SSR. Tree Genet Genomes. 2005, 1: 11-19. 10.1007/s11295-004-0002-9.Google Scholar
- Breitling R, Amtmann A, Herzyk P: Iterative Group Analysis (iGA): A simple tool to enhance sensitivity and facilitate interpretation of microarray experiments. BMC Bioinformatics. 2004, 5:Google Scholar
- Ihmels J, Bergmann S, Berman J, Barkai N: Comparative gene expression analysis by a differential clustering approach: Application to the Candida albicans transcription program. PLoS Genet. 2005, 1: 380-393. 10.1371/journal.pgen.0010039.Google Scholar
- Lelandais G, Tanty V, Geneix C, Etchebest C, Jacq C, Devaux F: Genome adaptation to chemical stress: clues from comparative transcriptomics in Saccharomyces cerevisiae and Candida glabrata. Genome Biol. 2008, 9: 10.1186/gb-2008-9-11-r164.Google Scholar
- Sjodin A, Street NR, Sandberg G, Gustafsson P, Jansson S: The Populus Genome Integrative Explorer (PopGenIE): a new resource for exploring the Populus genome. New Phytol. 2009, 182: 1013-1025. 10.1111/j.1469-8137.2009.02807.x.PubMedGoogle Scholar
- Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, Boerjan W, Rohde A: A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007, 19: 2370-2390. 10.1105/tpc.107.052811.PubMed CentralPubMedGoogle Scholar
- Baena-Gonzalez E, Rolland F, Thevelein JM, Sheen J: A central integrator of transcription networks in plant stress and energy signalling. Nature. 2007, 448: 938-U910. 10.1038/nature06069.PubMedGoogle Scholar
- Zhang YH, Primavesi LF, Jhurreea D, Andralojc PJ, Mitchell RAC, Powers SJ, Schluepmann H, Delatte T, Wingler A, Paul MJ: Inhibition of SNF1-Related Protein Kinase1 Activity and Regulation of Metabolic Pathways by Trehalose-6-Phosphate. Plant Physiol. 2009, 149: 1860-1871. 10.1104/pp.108.133934.PubMed CentralPubMedGoogle Scholar
- Baena-Gonzalez E, Sheen J: Convergent energy and stress signaling. Trends Plant Sci. 2008, 13: 474-482. 10.1016/j.tplants.2008.06.006.PubMed CentralPubMedGoogle Scholar
- Kant P, Kant S, Gordon M, Shaked R, Barak S: STRESS RESPONSE SUPPRESSOR1 and STRESS RESPONSE SUPPRESSOR2, two DEAD-box RNA helicases that attenuate Arabidopsis responses to multiple abiotic stresses. Plant Physiol. 2007, 145: 814-830. 10.1104/pp.107.099895.PubMed CentralPubMedGoogle Scholar
- Thompson AJ, Mulholland BJ, Jackson AC, McKee JMT, Hilton HW, Symonds RC, Sonneveld T, Burbidge A, Stevenson P, Taylor IB: Regulation and manipulation of ABA biosynthesis in roots. Plant Cell Environ. 2007, 30: 67-78. 10.1111/j.1365-3040.2006.01606.x.PubMedGoogle Scholar
- Davison PA, Hunter CN, Horton P: Overexpression of beta-carotene hydroxylase enhances stress tolerance in Arabidopsis. Nature. 2002, 418: 203-206. 10.1038/nature00861.PubMedGoogle Scholar
- Park SY, Fung P, Nishimura N, Jensen DR, Fujii H, Zhao Y, Lumba S, Santiago J, Rodrigues A, Chow TFF, et al: Abscisic Acid Inhibits Type 2C Protein Phosphatases via the PYR/PYL Family of START Proteins. Science. 2009, 324: 1068-1071.PubMed CentralPubMedGoogle Scholar
- Santiago J, Dupeux F, Round A, Antoni R, Park SY, Jamin M, Cutler SR, Rodriguez PL, Marquez JA: The abscisic acid receptor PYR1 in complex with abscisic acid. Nature. 2009, 462: 665-U143. 10.1038/nature08591.PubMedGoogle Scholar
- Yamaguchi-Shinozaki K, Shinozaki K: Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol. 2006, 57: 781-803. 10.1146/annurev.arplant.57.032905.105444.PubMedGoogle Scholar
- Olsson ASB, Engstrom P, Soderman E: The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol Biol. 2004, 55: 663-677. 10.1007/s11103-004-1581-4.PubMedGoogle Scholar
- Fujita M, Fujita Y, Maruyama K, Seki M, Hiratsu K, Ohme-Takagi M, Tran LSP, Yamaguchi-Shinozaki K, Shinozaki K: A dehydration-induced NAC protein, RD26, is involved in a novel ABA-dependent stress-signaling pathway. Plant J. 2004, 39: 863-876. 10.1111/j.1365-313X.2004.02171.x.PubMedGoogle Scholar
- Siefers N, Dang KK, Kumimoto RW, Bynum WE, Tayrose G, Holt BF: Tissue-Specific Expression Patterns of Arabidopsis NF-Y Transcription Factors Suggest Potential for Extensive Combinatorial Complexity. Plant Physiol. 2009, 149: 625-641. 10.1104/pp.108.130591.PubMed CentralPubMedGoogle Scholar
- Li WX, Oono Y, Zhu JH, He XJ, Wu JM, Iida K, Lu XY, Cui XP, Jin HL, Zhu JK: The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell. 2008, 20: 2238-2251. 10.1105/tpc.108.059444.PubMed CentralPubMedGoogle Scholar
- Nelson DE, Repetti PP, Adams TR, Creelman RA, Wu J, Warner DC, Anstrom DC, Bensen RJ, Castiglioni PP, Donnarummo MG, et al: Plant nuclear factor Y (NF-Y) B subunits confer drought tolerance and lead to improved corn yields on water-limited acres. Proc Natl Acad Sci USA. 2007, 104: 16450-16455. 10.1073/pnas.0707193104.PubMed CentralPubMedGoogle Scholar
- Sakamoto H, Maruyama K, Sakuma Y, Meshi T, Iwabuchi M, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2-type zinc-finger proteins function as transcription repressors under drought, cold, and high-salinity stress conditions. Plant Physiol. 2004, 136: 2734-2746. 10.1104/pp.104.046599.PubMed CentralPubMedGoogle Scholar
- Davletova S, Schlauch K, Coutu J, Mittler R: The zinc-finger protein Zat12 plays a central role in reactive oxygen and abiotic stress signaling in Arabidopsis. Plant Physiol. 2005, 139: 847-856. 10.1104/pp.105.068254.PubMed CentralPubMedGoogle Scholar
- Martin L, Leblanc-Fournier N, Azri W, Lenne C, Henry C, Coutand C, Julien JL: Characterization and expression analysis under bending and other abiotic factors of PtaZFP2, a poplar gene encoding a Cys2/His2 zinc finger protein. Tree Physiol. 2009, 29: 125-136. 10.1093/treephys/tpn011.PubMedGoogle Scholar
- Medina J, Ballesteros ML, Salinas J: Phylogenetic and functional analysis of Arabidopsis RCI2 genes. J Exp Bot. 2007, 58: 4333-4346. 10.1093/jxb/erm285.PubMedGoogle Scholar
- Mitsuya S, Taniguchi M, Miyake H, Takabe T: Disruption of RCI2A leads to over-accumulation of Na+ and increased salt sensitivity in Arabidopsis thaliana plants. Planta. 2005, 222: 1001-1009. 10.1007/s00425-005-0043-9.PubMedGoogle Scholar
- Breton G, Danyluk J, Charron JBF, Sarhan F: Expression profiling and bioinformatic analyses of a novel stress-regulated multispanning transmembrane protein family from cereals and Arabidopsis. Plant Physiol. 2003, 132: 64-74. 10.1104/pp.102.015255.PubMed CentralPubMedGoogle Scholar
- Liu FL, VanToai T, Moy LP, Bock G, Linford LD, Quackenbush J: Global transcription profiling reveals comprehensive insights into hypoxic response in Arabidopsis. Plant Physiol. 2005, 137: 1115-1129. 10.1104/pp.104.055475.PubMed CentralPubMedGoogle Scholar
- van Dongen JT, Frohlich A, Ramirez-Aguilar SJ, Schauer N, Fernie AR, Erban A, Kopka J, Clark J, Langer A, Geigenberger P: Transcript and metabolite profiling of the adaptive response to mild decreases in oxygen concentration in the roots of arabidopsis plants. Ann Bot. 2009, 103: 269-280. 10.1093/aob/mcn126.PubMed CentralPubMedGoogle Scholar
- Kreuzwieser J, Hauberg J, Howell KA, Carroll A, Rennenberg H, Millar AH, Whelan J: Differential Response of Gray Poplar Leaves and Roots Underpins Stress Adaptation during Hypoxia. Plant Physiol. 2009, 149: 461-473. 10.1104/pp.108.125989.PubMed CentralPubMedGoogle Scholar
- Matsui A, Ishida J, Morosawa T, Mochizuki Y, Kaminuma E, Endo TA, Okamoto M, Nambara E, Nakajima M, Kawashima M, et al: Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array. Plant Cell Physiol. 2008, 49: 1135-1149. 10.1093/pcp/pcn101.PubMedGoogle Scholar
- Gimeno J, Gadea J, Forment J, Perez-Valle J, Santiago J, Martinez-Godoy MA, Yenush L, Belles JM, Brumos J, Colmenero-Flores JM, et al: Shared and novel molecular responses of mandarin to drought. Plant Mol Biol. 2009, 70: 403-420. 10.1007/s11103-009-9481-2.PubMedGoogle Scholar
- Rubio S, Rodrigues A, Saez A, Dizon MB, Galle A, Kim TH, Santiago J, Flexas J, Schroeder JI, Rodriguez PL: Triple Loss of Function of Protein Phosphatases Type 2C Leads to Partial Constitutive Response to Endogenous Abscisic Acid. Plant Physiol. 2009, 150: 1345-1355. 10.1104/pp.109.137174.PubMed CentralPubMedGoogle Scholar
- Cheng WH, Chiang MH, Hwang SG, Lin PC: Antagonism between abscisic acid and ethylene in Arabidopsis acts in parallel with the reciprocal regulation of their metabolism and signaling pathways. Plant Mol Biol. 2009, 71: 61-80. 10.1007/s11103-009-9509-7.PubMed CentralPubMedGoogle Scholar
- Huang DQ, Wu WR, Abrams SR, Cutler AJ: The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. J Exp Bot. 2008, 59: 2991-3007. 10.1093/jxb/ern155.PubMed CentralPubMedGoogle Scholar
- Popko J, Hansch R, Mendel RR, Polle A, Teichmann T: The role of abscisic acid and auxin in the response of poplar to abiotic stress. Plant Biol. 2010, 12: 242-258. 10.1111/j.1438-8677.2009.00305.x.PubMedGoogle Scholar
- Xie Z, Zhang ZL, Zou XL, Huang J, Ruas P, Thompson D, Shen QJ: Annotations and functional analyses of the rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells. Plant Physiol. 2005, 137: 176-189. 10.1104/pp.104.054312.PubMed CentralPubMedGoogle Scholar
- Kankainen M, Brader G, Toronen P, Palva ET, Holm L: Identifying functional gene sets from hierarchically clustered expression data: map of abiotic stress regulated genes in Arabidopsis thaliana. Nucleic Acids Res. 2006, 34:Google Scholar
- Jiang YQ, Deyholos MK: Functional characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol Biol. 2009, 69: 91-105. 10.1007/s11103-008-9408-3.PubMedGoogle Scholar
- Eulgem T: Dissecting the WRKY web of plant defense regulators. PLoS Pathog. 2006, 2: 1028-1030. 10.1371/journal.ppat.0020126.Google Scholar
- Kagale S, Links MG, Rozwadowski K: Genome-Wide Analysis of Ethylene-Responsive Element Binding Factor-Associated Amphiphilic Repression Motif-Containing Transcriptional Regulators in Arabidopsis. Plant Physiol. 2010, 152: 1109-1134. 10.1104/pp.109.151704.PubMed CentralPubMedGoogle Scholar
- Meier S, Bastian R, Donaldson L, Murray S, Bajic V, Gehring C: Co-expression and promoter content analyses assign a role in biotic and abiotic stress responses to plant natriuretic peptides. BMC Plant Biol. 2008, 8: 10.1186/1471-2229-8-24.Google Scholar
- Morse M, Pironcheva G, Gehring C: AtPNP-A is a systemically mobile natriuretic peptide immunoanalogue with a role in Arabidopsis thaliana cell volume regulation. FEBS Lett. 2004, 556: 99-103. 10.1016/S0014-5793(03)01384-X.PubMedGoogle Scholar
- Ma SS, Bohnert HJ: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol. 2007, 8: 10.1186/gb-2007-8-4-r49.Google Scholar
- Jiang MY, Zhang JH: Water stress-induced abscisic acid accumulation triggers the increased generation of reactive oxygen species and up-regulates the activities of antioxidant enzymes in maize leaves. J Exp Bot. 2002, 53: 2401-2410. 10.1093/jxb/erf090.PubMedGoogle Scholar
- Umbach AL, Ng VS, Siedow JN: Regulation of plant alternative oxidase activity: A tale of two cysteines. Biochim Biophys Acta-Bioenerg. 2006, 1757: 135-142. 10.1016/j.bbabio.2005.12.005.Google Scholar
- Charron JBF, Ouellet F, Houde M, Sarhan F: The plant apolipoprotein D ortholog protects Arabidopsis against oxidative stress. BMC Plant Biol. 2008, 8: 10.1186/1471-2229-8-86.Google Scholar
- Baxter CJ, Redestig H, Schauer N, Repsilber D, Patil KR, Nielsen J, Selbig J, Liu JL, Fernie AR, Sweetlove LJ: The metabolic response of heterotrophic Arabidopsis cells to oxidative stress. Plant Physiol. 2007, 143: 312-325. 10.1104/pp.106.090431.PubMed CentralPubMedGoogle Scholar
- Woo HR, Kim JH, Nam HG, Lim PO: The delayed leaf senescence mutants of Arabidopsis, ore1, ore3, and ore9 are tolerant to oxidative stress. Plant Cell Physiol. 2004, 45: 923-932. 10.1093/pcp/pch110.PubMedGoogle Scholar
- Mowla SB, Cuypers A, Driscoll SP, Kiddle G, Thomson J, Foyer CH, Theodoulou FL: Yeast complementation reveals a role for an Arabidopsis thaliana late embryogenesis abundant (LEA)-like protein in oxidative stress tolerance. Plant J. 2006, 48: 743-756. 10.1111/j.1365-313X.2006.02911.x.PubMedGoogle Scholar
- Luhua S, Ciftci-Yilmaz S, Harper J, Cushman J, Mittler R: Enhanced tolerance to oxidative stress in transgenic Arabidopsis plants expressing proteins of unknown function. Plant Physiol. 2008, 148: 280-292. 10.1104/pp.108.124875.PubMed CentralPubMedGoogle Scholar
- Nishizawa A, Yabuta Y, Shigeoka S: Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008, 147: 1251-1263. 10.1104/pp.108.122465.PubMed CentralPubMedGoogle Scholar
- Bray EA: Classification of genes differentially expressed during water-deficit stress in Arabidopsis thaliana: an analysis using Microarray and differential expression data. Ann Bot. 2002, 89: 803-811. 10.1093/aob/mcf104.PubMed CentralPubMedGoogle Scholar
- Wang Z, Xie WQ, Chi F, Li CF: Identification of non-specific lipid transfer protein-1 as a calmodulin-binding protein in Arabidopsis. FEBS Lett. 2005, 579: 1683-1687. 10.1016/j.febslet.2005.02.024.PubMedGoogle Scholar
- Raffaele S, Mongrand S, Gamas P, Niebel A, Ott T: Genome-wide annotation of remorins, a plant-specific protein family: Evolutionary and functional perspectives. Plant Physiol. 2007, 145: 593-600. 10.1104/pp.107.108639.PubMed CentralPubMedGoogle Scholar
- Hundertmark M, Hincha DK: LEA (Late Embryogenesis Abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008, 9: 10.1186/1471-2164-9-118.Google Scholar
- Wang WX, Pelah D, Alergand T, Shoseyov O, Altman A: Characterization of SP1, a stress-responsive, boiling-soluble, homo-oligomeric protein from aspen. Plant Physiol. 2002, 130: 865-875. 10.1104/pp.002436.PubMed CentralPubMedGoogle Scholar
- Rodriguez-Lopez M, Baroja-Fernandez E, Zandueta-Criado A, Moreno-Bruna B, Munoz FJ, Akazawa T, Pozueta-Romero J: Two isoforms of a nucleotide-sugar pyrophosphatase/phosphodiesterase from barley leaves (Hordeum vulgare L.) are distinct oligomers of HvGLP1, a germin-like protein. FEBS Lett. 2001, 490: 44-48. 10.1016/S0014-5793(01)02135-4.PubMedGoogle Scholar
- Davidson RM, Reeves PA, Manosalva PM, Leach JE: Germins: A diverse protein family important for crop improvement. Plant Sci. 2009, 177: 499-510. 10.1016/j.plantsci.2009.08.012.Google Scholar
- Xin ZG, Mandaokar A, Chen JP, Last RL, Browse J: Arabidopsis ESK1 encodes a novel regulator of freezing tolerance. Plant J. 2007, 49: 786-799. 10.1111/j.1365-313X.2006.02994.x.PubMedGoogle Scholar
- Bouchabke-Coussa O, Quashie ML, Seoane-Redondo J, Fortabat MN, Gery C, Yu A, Linderme D, Trouverie J, Granier F, Teoule E, Durand-Tardif M: ESKIMO1 is a key gene involved in water economy as well as cold acclimation and salt tolerance. BMC Plant Biol. 2008, 8: 10.1186/1471-2229-8-125.Google Scholar
- Borevitz JO, Liang D, Plouffe D, Chang HS, Zhu T, Weigel D, Berry CC, Winzeler E, Chory J: Large-scale identification of single-feature polymorphisms in complex genomes. Genome Res. 2003, 13: 513-523. 10.1101/gr.541303.PubMed CentralPubMedGoogle Scholar
- Davey MW, Graham NS, Vanholme B, Swennen R, May ST, Keulemans J: Heterologous oligonucleotide microarrays for transcriptomics in a non-model species; a proof-of-concept study of drought stress in Musa. BMC Genomics. 2009, 10: 10.1186/1471-2164-10-436.Google Scholar
- Caceres M, Lachuer J, Zapala MA, Redmond JC, Kudo L, Geschwind DH, Lockhart DJ, Preuss TM, Barlow C: Elevated gene expression levels distinguish human from non-human primate brains. Proc Natl Acad Sci USA. 2003, 100: 13030-13035. 10.1073/pnas.2135499100.PubMed CentralPubMedGoogle Scholar
- 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. 10.1093/biostatistics/4.2.249.PubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge YC, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: 10.1186/gb-2004-5-10-r80.Google Scholar
- Ge X, Tsutsumi S, Aburatani H, Iwata S: Reducing false positives in molecular pattern recognition. Genome Inform. 2003, 14: 34-43.PubMedGoogle Scholar
- Merret R, Moulia B, Hummel I, Cohen D, Dreyer E, Bogeat-Triboulot MB: Monitoring the regulation of gene expression in a growing organ using a fluid mechanics formalism. BMC Biol. 2010, 8: 10.1186/1741-7007-8-18.Google Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: W64-W70. 10.1093/nar/gkq310.PubMed CentralPubMedGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, et al: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008, 36: D1009-D1014. 10.1093/nar/gkm965.PubMed CentralPubMedGoogle Scholar
- Armengaud P, Breitling R, Amtmann A: The potassium-dependent transcriptome of Arabidopsis reveals a prominent role of jasmonic acid in nutrient signaling. Plant Physiol. 2004, 136: 2556-2576. 10.1104/pp.104.046482.PubMed CentralPubMedGoogle Scholar
- Caruso A, Chefdor F, Carpin S, Depierreux C, Delmotte FM, Kahlem G, Morabito D: Physiological characterization and identification of genes differentially expressed in response to drought induced by PEG 6000 in Populus canadensis leaves. J Plant Physiol. 2008, 165: 932-941. 10.1016/j.jplph.2007.04.006.PubMedGoogle Scholar
- Quesada T, Li Z, Dervinis C, Li Y, Bocock PN, Tuskan GA, Casella G, Davis JM, Kirst M: Comparative analysis of the transcriptomes of Populus trichocarpa and Arabidopsis thaliana suggests extensive evolution of gene expression regulation in angiosperms. New Phytol. 2008, 180: 408-420. 10.1111/j.1469-8137.2008.02586.x.PubMedGoogle Scholar
- Segerman B, Jansson S, Karlsson J: Characterization of genes with tissue-specific differential expression patterns in Populus. Tree Genet Genomes. 2007, 3: 351-362. 10.1007/s11295-006-0077-6.Google Scholar
- Sjodin A, Wissel K, Bylesjo M, Trygg J, Jansson S: Global expression profiling in leaves of free-growing aspen. BMC Plant Biol. 2008, 8: 10.1186/1471-2229-8-61.Google Scholar
- Sjodin A, Bylesjo M, Skogstrom O, Eriksson D, Nilsson P, Ryden P, Jansson S, Karlsson J: UPSC-BASE-Populus transcriptomics online. Plant J. 2006, 48: 806-817. 10.1111/j.1365-313X.2006.02920.x.PubMedGoogle Scholar
- Baena-Gonzalez E: Energy Signaling in the Regulation of Gene Expression during Stress. Mol Plant. 2010, 3: 300-313. 10.1093/mp/ssp113.PubMedGoogle Scholar
- Seo M, Koshiba T: Complex regulation of ABA biosynthesis in plants. Trends Plant Sci. 2002, 7: 41-48. 10.1016/S1360-1385(01)02187-2.PubMedGoogle Scholar
- Ko JH, Yang SH, Han KH: Upregulation of an Arabidopsis RING-H2 gene, XERICO, confers drought tolerance through increased abscisic acid biosynthesis. Plant J. 2006, 47: 343-355. 10.1111/j.1365-313X.2006.02782.x.PubMedGoogle Scholar
- Himmelbach A, Hoffmann T, Leube M, Hohener B, Grill E: Homeodomain protein ATHB6 is a target of the protein phosphatase ABI1 and regulates hormone responses in Arabidopsis. Embo J. 2002, 21: 3029-3038. 10.1093/emboj/cdf316.PubMed CentralPubMedGoogle Scholar
- Kariola T, Brader G, Helenius E, Li J, Heino P, Palva ET: EARLY RESPONSIVE TO DEHYDRATION 15, a negative regulator of abscisic acid responses in Arabidopsis. Plant Physiol. 2006, 142: 1559-1573. 10.1104/pp.106.086223.PubMed CentralPubMedGoogle Scholar
- Pandey GK, Grant JJ, Cheong YH, Kim BG, Li LG, Luan S: ABR1, an APETALA2-domain transcription factor that functions as a repressor of ABA response in Arabidopsis. Plant Physiol. 2005, 139: 1185-1193. 10.1104/pp.105.066324.PubMed CentralPubMedGoogle Scholar
This article is published under license to 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.