Correlation Network Analysis reveals a sequential reorganization of metabolic and transcriptional states during germination and gene-metabolite relationships in developing seedlings of Arabidopsis
© Allen et al; licensee BioMed Central Ltd. 2010
Received: 12 October 2009
Accepted: 13 May 2010
Published: 13 May 2010
Holistic profiling and systems biology studies of nutrient availability are providing more and more insight into the mechanisms by which gene expression responds to diverse nutrients and metabolites. Less is known about the mechanisms by which gene expression is affected by endogenous metabolites, which can change dramatically during development. Multivariate statistics and correlation network analysis approaches were applied to non-targeted profiling data to investigate transcriptional and metabolic states and to identify metabolites potentially influencing gene expression during the heterotrophic to autotrophic transition of seedling establishment.
Microarray-based transcript profiles were obtained from extracts of Arabidopsis seeds or seedlings harvested from imbibition to eight days-old. 1H-NMR metabolite profiles were obtained for corresponding samples. Analysis of transcript data revealed high differential gene expression through seedling emergence followed by a period of less change. Differential gene expression increased gradually to day 8, and showed two days, 5 and 7, with a very high proportion of up-regulated genes, including transcription factor/signaling genes. Network cartography using spring embedding revealed two primary clusters of highly correlated metabolites, which appear to reflect temporally distinct metabolic states. Principle Component Analyses of both sets of profiling data produced a chronological spread of time points, which would be expected of a developmental series. The network cartography of the transcript data produced two distinct clusters comprising days 0 to 2 and days 3 to 8, whereas the corresponding analysis of metabolite data revealed a shift of day 2 into the day 3 to 8 group. A metabolite and transcript pair-wise correlation analysis encompassing all time points gave a set of 237 highly significant correlations. Of 129 genes correlated to sucrose, 44 of them were known to be sucrose responsive including a number of transcription factors.
Microarray analysis during germination and establishment revealed major transitions in transcriptional activity at time points potentially associated with developmental transitions. Network cartography using spring-embedding indicate that a shift in the state of nutritionally important metabolites precedes a major shift in the transcriptional state going from germination to seedling emergence. Pair-wise linear correlations of transcript and metabolite levels identified many genes known to be influenced by metabolites, and provided other targets to investigate metabolite regulation of gene expression during seedling establishment.
Germination is a phenomenon with complex regulation that is a balance between the release of dormancy and the promotion of germination. This reflects the relationship between the hormones gibberellic acid (GA) and abscisic acid (ABA), environmental cues [1, 2], and the spatial distribution of hormone action and gene expression [3–5]. Considerable effort has been put into elucidating the molecular mechanisms controlling seed germination with greater application of gene expression profiling. These studies have highlighted the roles of gene expression changes in mediating GA and ABA interactions in the control of dormancy and germination [6–10]. To complement the growing number of gene expression studies, Fait et al.  conducted an integrated metabolomic and gene expression study of various seed developmental stages from maturation through germination. They identified distinct metabolite profiles associated with the various developmental stages and suggested that seeds are metabolically primed for germination during desiccation and subsequent metabolic programming during imbibition and germination is essential for seedling establishment. An integrated metabolomic and transcriptomic study of photomorphogenesis in red light and far-red light treated seedlings showed that even though transcript profiles were relatively similar, phenotypic differences could be explained by significant differences at the level of the metabolome .
Prior to seed germination, the mobilization of stored triacylglycerol (TAG) begins in earnest in order to feed the developing seedling. The processes by which germination and lipid mobilization are regulated have been found to be distinct , and it is likely that reserve mobilization is governed by abscisic acid-related processes within the embryo . In Arabidopsis, stored sugars are consumed by the time the radicle has emerged, and within 48 h after germination lipid and protein stores have been consumed . At this point, the seedling must become photosynthetically competent. It has been suggested that metabolic signals may regulate the transition from heterotrophy to autotrophy in seedlings in order to maximize the use of storage compounds . Exploiting the altered behavior of seed germination and of seedling vigor for forward genetic screens of Arabidopsis mutants has been instrumental in revealing the potential signaling properties of metabolites, primarily sugars , and nutrients . Mutant studies have revealed the interaction of sugars and hormones [18, 19] and the concept of a carbon:nitrogen 'matrix effect' in metabolic regulation . Through a forward genetic screen using the toxic analogue monofluoroacetic acid, we identified mutants disrupted in their ability to metabolize exogenously supplied acetate through the glyoxylate cycle [21, 22]. A physiological analysis of the mutants provided evidence that carbohydrate responses of seedlings may be impaired within the mutants. This suggests a cross-talk between organic acid and carbohydrate signaling in developing seedlings  with the possibility of either acetate or down-stream metabolites influencing gene expression in developing seedlings.
Many forward genetic screens have relied on observing differential sensitivity of mutants to added compounds. This approach does not work for many metabolites, since artificially high concentrations must be used and undesired traits are selected. For example, organic acids pose this problem because they are weak acids, and mutant selection for specific responses may be confounded by responses to altered intracellular pH. Integrated analysis of metabolite and transcript data offers a way to identify co-regulatory networks of metabolites and genes [23, 24]. This has been applied successfully to identify potential genetic regulation of metabolite levels concerning sulfur stress [25–29], glucosinolate metabolism , and nitrogen responses [31, 32] in Arabidopsis and fruit development in tomato [28, 33–35]. The suggestion that strong correlations between metabolites and transcripts may reflect metabolite effects on gene expression [27, 28, 36], therefore, enables integrated analysis to be used to identify potential signaling metabolites for subsequent detailed studies. We obtained metabolite and transcript profiling data from a series of samples spanning germination and establishment, and analyzed the data to identify pair-wise combinations of genes and metabolites strongly correlated over this developmental transition. We discuss how analysis of metabolite-gene correlations provided evidence for differential regulation of a common ontological class of genes. Furthermore, the network correlation analysis approach can provide supplemental information on the progression of metabolic and transcriptional states during developmental transitions [27, 28]. Both types of profiling data were mined for interesting gene expression and metabolite patterns and relationships. Principle Component Analysis (PCA) and network correlation analysis based on spring-embedding  were used to integrate and visualize the data to obtain information about the metabolic and transitional states present during germination and seedling development.
Gene expression during seedling development
Behavior of metabolites during development
A total of 27 metabolites corresponding to a variety of known and unknown metabolites including four soluble carbohydrates, nine amino acids and five organic acids were quantified from the 1H-NMR spectra of seedling extracts. Although these metabolites comprise a small proportion of the total metabolic complement of a cell, these metabolites are the most abundant ones. They reflect the nutritional state of the tissue as an immediate source of carbon and/or nitrogen and serve as respiratory substrates for energy production. A direct comparison of data from our NMR profiling platform with GC-MS acquired data demonstrated a similar capacity to distinguish metabolic states . Additionally, a number of the metabolites are well known effectors of gene expression and some, such as sucrose, isoleucine and glutamine, have high regulatory potential as determined by correlation analysis . The values obtained for each metabolite are given in Additional File 1.
The relationships between individual metabolites are clearer when correlations are included in the visualization as shown in the spring embedding plots (Fig. 2B-C). Our threshold p-value produced a correlation coefficient cut-off value of 0.68. We observed significant correlations between 19 metabolites that appear to be separated into three clusters. The cluster containing glutamine, fructose, fumarate and unkD8.0 are linked to the lactate cluster (Fig. 2B) via sucrose. The lactate cluster contains those metabolites that are decreasing over time, such as trigonelline, threonine, citrate, and alanine (Fig. 2C). Malate has also been included within this cluster and has a relatively high correlation to alanine. This could reflect a partitioning of malate into alanine either via oxaloacetate or pyruvate. The third cluster (Fig. 2D) consists of the aliphatic amino acids and the compounds choline, and formate. It was expected that the three amino acids leucine, isoleucine and valine would be highly correlated as they share common synthetic and catabolic pathways.
Transcriptional states of developing seedlings
Metabolic states of developing seedlings
Metabolite and transcript co-analysis
Correlations between metabolites and transcript levels in developing seedlings.
Number of Relationships*
Direction of Correlations†
N (1), P (2)
N (6), P (1)
N (2), P (1)
N (1), P (3)
N (15), P (1)
N (2), P (4)
N (1), P (1)
N (18), P (1)
Identified regulatory genes correlated with metabolites.
Suppressor of Frigida 4
Myc-related bHLH transcription factor
C2H2 zinc finger protein FZF
Early Flowering 3
Auxin Response Factor 10
bZIP transcription factor POSF21
Nuclear localized bHLH protein that interacts with active PhyB
Arabidopsis thaliana TCP
Family transcription factor
Homologous to the maize transcription factor Viviparous-1.
IAA14 is a member of the Aux/IAA protein family. Solitary Root locus.
High Mobility Group B 5
Member of the R2R3 factor gene family
Callus Expression of RBCS
Vernalization Independence 5 putative heterogeneous nuclear ribonucleoprotein
DNA binding, Chloroplast
Arabidopsis Pumilio 24
Homologous to the DNA-directed RNA polymerase II subunit (At5g59180)
arginine/serine-rich splicing factor
Cold shock domain protein
ARABIDOPSIS MEI2-LIKE 4
Arabidopsis thaliana Bromodomain and Extra terminal Domain protein 9
Polypyrimidine tract-binding (PTB) Protein
An example of the relationships between a metabolite and a connected gene is given for lactate (Additional File 1). Lactate was negatively correlated with 21 genes, five of which had no assigned identity in the TAIR database. Of the 17 genes with an assigned identity, seven showed an involvement in photosynthetic-related functions. Eight of the genes within the group are also affected by abiotic or biotic stress. These include 1-aminocyclopropane-1-carboxylate oxidate (ACO4, At1g05010), calmodulin-like 9 (At3g51920) and the PIP2A aquaporin (At3g53420), which is induced during dehydration stress. Due to our stringent statistical cut-off, none of the 21 genes of the lactate cluster include any of the 19 hypoxia inducible genes reported by Loreti . However, the genes encoding alanine aminotransferase (At1g17290), alcohol dehydrogenase (At1g77120), Class I non-symbiotic hemoglobin (At2g16060), and pyruvate decarboxylase 1 (At4g33070) reported as anoxia inducible by Sachs et al.  - and which are included in the set of 19 inducible genes - were positively correlated with lactate at a p-value less than 0.03 (r > 0.72). Although a number of gene expression profiles were correlated with more than one metabolite concentration, it was observed that of the seven photosynthetic genes correlated to lactate, only three showed correlations with other metabolites. Besides the two TFs shown in Table 2, the unknown metabolite unkM1.85 was negatively correlated with 17 other genes, seven of which were assigned photosynthetic functions.
More and more, integrative approaches are being employed to describe the function of molecular systems in development (for reviews see [24, 46–49]. Whereas most metabolite-gene interaction studies have been from the point of view of understanding the genetic bases for changes in metabolism, such studies can be integral to understanding the control of gene expression by metabolic factors [28, 49]. In fact, strong correlations between metabolite and transcript levels more likely reflect metabolite regulation of transcription than vice versa. A recent study reported that a series of distinct metabolic switches were characteristic of the transition from dormant, dry seed to germinating embryo . The results presented in this work extend the analysis to provide an overview of metabolite and transcriptome profiles from imbibed, non-dormant seeds through to established seedlings with the aim of identifying potential targets of metabolic control during seedling development covering the heterotrophic to autotrophic transition.
Differential gene expression activity coincides with developmental stage
In the present work, we observed a relatively high level of transcriptional change occurring over the first three days of seedling development, which encompasses germination and emergence. Fifty percent of all DE events occurred during the first three days. Of these, more than 7% were genes categorized as either TF or signaling genes. This implies large changes in transcriptional activity during emergence. Alba et al.  made a similar conclusion from their analysis of gene expression in developing pericarp of ethylene treated tomatoes. Of 628 known DE genes during the 10 developmental stages they analyzed, 11% were either TF or signaling genes. We observed a substantial increase in DE at day 3 compared to day 2, with most genes comprising both chloroplast/plastid and TF/signaling genes. A high degree of transcriptional alteration may not be required for seedling development (i.e. cell division or differentiation) at this time, since they are geared for a constant rate of lipid degradation [50, 51], and cell expansion is the principal means of seedling growth . Subsequently, relative gene UR and DR appears to follow a cyclic pattern during the subsequent days, and it is interesting that this corresponds to likely transitional stages of development. Transcriptional transition is lowest between days 2 and day 4, since TF gene expression is DR. By day 4, all lipid reserves have been depleted and so by day 5 any potential catabolite repression would be eliminated to permit full development of autotrophy, which would be revealed by a relative increase in UR genes, such as those encoding chloroplast/photosynthetic proteins. Between days 6 and 7 rapid leaf growth begins  and a corresponding UR of gene expression may ensue. Accordingly, spikes of chloroplast/plastid and TF/signaling UR take place at days 5 and 7.
Metabolic state establishes prior to germination and the switch in transcriptional programming
The changes to levels of various metabolites going from imbibition to early germination follow similar patterns as reported previously [11, 13]. Fait et al.  observed a change of metabolic activity during post-imbibition germination 24 h after transfer of seedlings from cold to a germination inductive temperature. The grouping of the metabolite profiles produced during this experiment supports these findings, demonstrating that in the cold, a relatively stable metabolic state for the major metabolites is present and then changes relatively little for 24 h. A larger metabolic switch occurs from day 1 to day 2 but the metabolic state stabilizes during seedling establishment even though a number of metabolites show transient increases (Fig. 4).
Although gene expression profiles are changing from day 0 to day 2, there is a more dramatic change from day 2 to day 3 (Fig. 3). By this point, the seedling has emerged, but attainment of full photosynthetic competence does not appear to happen quickly. From the large DR of expression from day 3 to day 4, it is interesting to speculate that prior to the emergence of the radicle, i.e. by day 2, the embryo has attained a metabolic state that primes the seedling to reduce aspects of gene expression in preference for emergence and reserve mobilization.
Revealing potential metabolic signals by correlation analysis
Deciphering metabolic contributions to switches in transcriptional states, such as observed during seedling development, will entail identifying individual signaling metabolites, the genes they affect and the concerted degree of affect . Although any metabolic regulation during the heterotrophic to autotrophic transition would be complex, it should be possible to identify metabolites involved in signaling gene expression events by examining their behavior in relation to the expression of specific genes . We determined linear correlations between each metabolite-gene pair with the assumption that the strength of correlation would indicate the potential for a regulatory relationship to exist (Fig. 5). Sucrose levels showed positive correlations with 129 gene transcripts. Comparison of the gene transcripts correlated with sucrose levels with previous microarray experiments and online databases showed that 44 of the 129 genes had previously been identified as sucrose-responsive [36, 53–56]. The correlation of sucrose levels with a large proportion of previously identified sucrose-responsive gene transcripts reinforces the validity of the use of correlation coefficients to identify interesting relationships. Two well-studied TFs that were highly correlated with sucrose were PIF1 and ABI3. PIF1 has been identified as a negative regulator of photomorphogenesis in seedlings [57, 58] and ABI3 may play a role in sugar-induced seedling developmental arrest . The TF genes IAA14 and ARF10 showed high negative correlation with the unknown metabolites unkM1.85 and unkM7.9, respectively. In order for correlations to be identified a gene had to be expressed at each time point. Therefore, it may act in some capacity outside the developmental stage in which it is commonly associated. For example, Penfield et al.  concluded that factors controlling cotyledon expansion in imbibed seeds -- a gibberellin mediated processes -- continue well into seedling establishment. The CHO1 AP2 domain TF that functions in the glucose signaling pathway downstream of ABI4 also appears to function well into seedling establishment . In support of these views, it is interesting that we observe expression of genes known to be involved in the imposition of dormancy well into seedling establishment.
ABI3 expression decreases within two days of imbibition to remove dormancy and permit seed germination. ARF10 is believed to increase seed sensitivity to ABA  and its delayed suppression would likely contribute to a graded control of ABA responses while ABI3 transcript levels decline. In contrast, the levels of SLR/IAA14 transcripts are increasing until day 4-5. SLR/IAA14 is known to repress ARF7 and ARF19 during initiation of lateral roots  and it is interesting to speculate that it may play an additional role to restrict the expression of ARF10 to vascular tissue in cotyledons and roots in older seedlings . Even if metabolic regulation of highly correlated genes is shown not to occur or is minimal, observing the behavior of expression within a broader physiological and biochemical context through a network correlation analysis may reveal as of yet unknown interactions.
Identifying mechanisms of metabolic regulation
Imbibition results in seeds undergoing a period of anoxia during which lactate production occurs [1, 41]. In animals, elevated lactate has been shown to alter gene expression in certain tumor types [64–66] and may involve carbohydrate response-like elements . We looked for elements within the promoters of photosynthetic genes correlated to lactate and unkM1.85 as a start toward identifying regulatory mechanisms (Supplemental Table 2) in a manner analogous to co-regulated genes identified by microarray analysis . Genes correlated with lactate contained Ocs-like elements responsive to oxidative stress, auxin and salicylic acid [69, 70] and motifs showing similarity to those involved in light-responses. Since it is difficult to distinguish between the effects of light and metabolic stimulus [71, 72], it is possible that elements identified as light responsive might be metabolite responsive instead. Interestingly, the potential promoter motifs identified in the photosynthetic genes correlated with metabolite unkM1.85 predominantly included elements associated with response to various stresses and did not contain any light-responsive elements. The differences between the potential promoter motifs identified in the two sets of photosynthetic genes indicate that distinct regulatory mechanisms may be operate in groups of genes that may be considered initially as functionally similar through ontological classification. Identification of the TFs that bind to these motifs, and the characterization of identified, but unknown promoter elements will help elucidate the signaling pathways involved in the expression of these genes and potential involvement of metabolic regulation.
The pair-wise analysis of metabolite and transcript levels appears to be a useful investigatory tool to identify potential links between genes and metabolites, thereby providing a number of targets for further examination. However, the identification of a correlated gene and metabolite does not provide information relating to the causality within the relationship, i.e. metabolite affecting gene expression or vice versa. It is also difficult to determine whether the observed relationship results from a direct interaction between a gene and a metabolite, or whether a downstream signaling event is involved. Such questions can be addressed, in part, by repeating the metabolite measurements in the appropriate mutant and/or by direct measurement of transcript levels in rigorously controlled metabolite feeding experiments.
A systems biology approach was adopted to investigate the interactions of metabolites and gene expression during seedling development. Both transcript and metabolite data were analysed at various levels and the results visualized using PCA and correlation-based network cartography. The analysis of transcript data alone showed that germination and seedling development is marked by stages of differing gene expression activity. These stages fall at important developmental points, such as at the beginning of seedling emergence, the end of reserve mobilization and the onset of leaf formation. Metabolite levels were revealed to fall into two clusters that reflect the pattern and timing of change, principally those that decrease post-imbibition and those that show a transient increase after seedling emergence. Network cartography, whereby the degree of correlation between variables was used as the basis of sample comparison, provided a clearer picture of the relationship among samples than PCA. This network analysis indicates a shift in the state of nutritionally important metabolites precedes the major shift in transcriptional state going from germination to seedling emergence. Therefore, a suitable metabolic state achieved prior to germination may be necessary for the initiation of gene expression programs for efficient seedling development. Some aspects of gene expression may be regulated by specific metabolites. The key is to identify signaling metabolites and the genes they affect, which may be accomplished by holistic profiling and correlation analysis. In addition, network correlation analysis may reveal component interactions when visualised within the context of a dependent or regulatory process, such as we noted with potential TF relationships uncovered by metabolite correlations.
Plant material and growth conditions
All Arabidopsis thaliana L. (ecotype Col-0) seeds were surface sterilized and sown onto 0.8% agar media plates containing 1/2-strength MS salts, pH 5.7 . The sowing density was approximately 500 seeds evenly spread on a 9 cm Petri dish. Agar and MS salts were purchased from Fisher Scientific and Sigma, respectively. The edges of each plate were wrapped in 3 MM surgical tape and the plates were incubated in the dark at 4°C for 4 days before transferring to the growth room. Transfer of plates occurred at 09:00 h. Seeds were germinated at 20°C at 70 μmol of photons m-2 s-1 constant illumination using standard white fluorescent bulbs (General Electric). A drop of less 2 μmol of photons m-2 s-1 was observed going from the centre to the edges of the shelf. We used seeds from completely brown siliques that had been after-ripened for at least 3-weeks after harvesting. The seeds were imbibed for 4 days prior to transfer to the growth room with dormancy being maintained by incubation at low temperature. Images of the stages at which we selected seedlings for analysis is given in Rylott et al. .
Design of tissue sampling
Sample harvesting and preparation was conducted in three sets of nine samples with each set encompassing the time points Day 0 to Day 8 after transfer to the growth room . We alternated tissue harvesting regimes to obtain sets of tissue for total RNA and metabolite extraction, thus corresponding tissue samples were used to compare metabolite and transcript profiles. A shelf in the growth room was divided into three sections and plates for each set were arranged horizontally around the shelf, such that each section contained an equal number of plates. Tissue was harvested each day at 09:00 h and only one sample per day was harvested resulting in three biological replicates for each time point. Each tissue sample consisted of pooled seedlings from an equal number of plates from each section. The total number of plates selected for each sample varied depending on the developmental stage. Each plate was examined under a microscope to ensure that seedlings were harvested at the required stage of development. Only those plates with greater than 95% of seedlings at the appropriate developmental stage were used. For sample days 0-4, approximately 0.4 g of seeds or seedlings was washed from the surface of the agar petri dishes with distilled sterile water into a filtration unit. Once the water had passed through, the seedlings were washed in 10 ml more water, weighed and immediately frozen in liquid nitrogen. For sample days 5-8, approximately 0.4 to 0.5 g of seedlings were removed from plates by forceps, rinsed briefly, and immediately frozen in liquid nitrogen. The time from opening the petri dish to freezing the sample was at most 3 min with rinsing times for all samples being within one and a half to two minutes.
Transcript profiling and DE Analysis
Total RNA was isolated using a borate-based extraction protocol . Production of labelled cDNA, quality checking, and slide hybridization were conducted as described in by Armengaud et al. . For each labelling reaction 1 μg of Oligo dT20 primer was added to 100 μg of total RNA in a total reaction volume of 20.5 μl. Printed 70-mer oligonucleotide microarrays were obtained from the laboratory of Prof. David Galbraith at the University of Arizona. Versions 1 and 3 arrays containing 29 K elements were used in these experiments. The identity of each spot in the meta-grid was obtained from the Galbraith laboratory http://www.ag.arizona.edu/microarray/. Only one cDNA sample was hybridized per slide. Since transcriptome profiles were produced from more than one microarray print version, only those genes common to all microarrays were used in this analysis.
Hybridized slides were scanned using an Affymetrix 428 scanner set on a gain setting to yield no more than 10 saturated spots per slide and gain settings were varied to account for the quality of the hybridization. Spot checking and intensity determination were done using ImaGene™ (BioDiscovery Inc., CA, USA). The quantified gene expression data produced by ImaGene was normalized using GeneSight™ version 4.1 (BioDiscovery Ltd.). Background signals were subtracted and spots designated as poor hybridization events were discounted from future analysis. In order to address the problem of negative spots, signal intensities below a set value of 20 were raised to that value. The raw expression data after spot removal has been deposited into the ArrayExpress database under the accession number M-MEXP-2493 http://www.ebi.ac.uk/microarray-as/ae/. A standard normalization procedure was applied to the quantified gene expression values obtained for each printed microarray to facilitate comparisons among individual microarrays (Affymetrix GeneChip Expression Analysis technical Manual, 2004). In brief, the top and bottom 2% of the signal values were removed and the mean calculated for the 96% of the values remaining. A value, the scaling factor, was calculated to adjust the mean of the remaining values to 100. Each of the signal intensities on the array was then multiplied by the appropriate scaling factor to normalize signal intensities on an array-by-array basis. A file of combined normalized data is also available from ArrayExpress under the accession number E-MEXP-2493. A group of differentially up-regulated and down-regulated genes was identified between each day at a confidence interval of 99% based on the bootstrapping procedure reported by Kerr & Churchill , which is resident within Genesight™. The groups of differentially expressed genes were filtered further with a threshold cut-off value of 1.5-fold. Functional analysis by ontology was done using the information obtained from the TAIR database http://www.arabidopsis.org.
In order to verify the robustness of our hybridizations and data normalization, we compared the relative expression levels of the genes encoding ACT2 (ACT2, At3g18780) and ribosomal protein genes S9 (At1g74970) and L32 (At5g46430) from the arrays to expression values obtained by qRT-PCR. The array and RT-PCR values were normalized by the expression levels of ubiquitin 10 (At4g05320) to allow direct comparison of data from the two different quantitative techniques. The levels of transcripts in each extract were determined by real-time RT-PCR according to Love et al. . Transcript values for each reference gene were obtained using a standard curve produced from purified plasmid DNA containing the appropriate cDNA. All genes (except UBQ10) were chosen, because each was spotted between 13 and 27 times throughout the microarray depending on the version used, and thus the average intensity would provide a general indication of the quality of the hybridization. These genes were selected for qRT-PCR analysis prior to hybridization and processing of the arrays, therefore they represent an unbiased indication of the quality of the array hybridization. A plot of array versus qRT-PCR values produced lines with high correlation (near unity for ACT2 and RPS9) demonstrating data of suitable quality for subsequent statistical analysis (Additional File 1).
Metabolite extraction and quantification
Metabolites were extracted according to Weckwerth et al.  and quantified by 1H-NMR as described in Moing et al. . Briefly, the dried extracts were resuspended in 400 mM phosphate buffer pH 6.0 in D2O and analyzed at 500.162 MHz on a Bruker spectrometer (Bruker Biospin Avance). Special care was taken to allow absolute quantification of the individual metabolites through addition of ethylene diamine tetraacetic acid sodium salt solution (5 mM final concentration in the NMR tube) to improve the resolution and quantification of organic acids such as malate and citrate, adequate choice of the NMR acquisition parameters (pulse angle 90°, relaxation delay 10 s) and use of an electronic reference (ERETIC mode ) calibrated with glucose, fructose, glutamine and glutamic acid sodium salt solutions as described previously . Individual metabolites were identified using published data [78, 80, 81], acquisition of NMR spectra of reference compounds under exact solvent conditions, and spiking extracts with reference compounds. They were quantified using the metabolite mode of AMIX software (Bruker Biospin v. 3.5.6) based on the number of protons comprising the corresponding resonance. Concentrations in the NMR tube were converted to amounts per g fresh weight using the mass of sample extracted. Citrate, formate, fumarate, glutamate, lactate and malate are expressed as μg of the acid form. The concentration of NMR unknown compounds (named according to the form of the resonance, S for singlet, D for doublet, M for multiplet, and its frequency in ppm) was calculated on the assumption that the measured resonance corresponded to one proton and using an arbitrary molecular weight of 100 Da. We verified the robustness of the quantification procedure by observing a near 1 to 1 relationship between levels of metabolites when we compared those measured in a mixture of samples, one from each day, with the corresponding calculated theoretical mix (Additional File 1).
Data analysis for network cartography
For the network cartography and correlation analysis, only those genes which had at least one expression value for each of the sample days were included in the analysis. This pre-processing step produced a set of 10,005 genes. Similar data pre-processing was not required for metabolite levels as they had been quantified absolutely for each extract. Since material for both the transcript and metabolite profiling was collected in three independent groups, the values were averaged for each day. Hierarchical, K-means and 2D-SOM clustering were done using the metabolite data imported into Genesight ™. PCA was performed using MATLAB™ release 2007b ((The MathWorks, Natick, MA). Visualization of significant correlations, i.e. network cartography, was conducted using springScape . Pearson correlation coefficients (signed r value) were used to generate the similarity matrices for the spring embedding of metabolite, transcript and combined data. MATLAB™ release 2007a was the mathematics platform for the spring embedding and metabolite-gene correlation analysis. Depending on the data set being analysed, the initial similarity matrix was cut at a threshold value to facilitate the spring embedding and to enhance the significance of the output correlations. For the similarity matrix of correlations between the 27 individual metabolites, a Bonferroni adjustment was applied based on a significance value of 0.1 and 351 tests. The day-by-metabolite and the day-by-gene similarity matrices were adjusted strictly by correlation coefficient in order to compare directly clustering between the two data sets. A FDR  was used to threshold the level of significance for the metabolite by gene correlations based on a significance value of 0.1 and 276,345 (27 × 10,235) tests. Determination of metabolites differentially expressed between the two groups of days was conducted by Student's t-tests also using a FDR of 0.1 (6 × 20 tests) to determine significance.
differential expression or differentially expressed
false discovery rate
principal component analysis
quantitative real-time PCR
transcription factor: UR: up-regulated (regulation)
We would like to thank the laboratory of Professor David Galbraith at Arizona State University for producing and distributing the Arabidopsis microarrays, Dr. Joel Milner and Janet Laird at Glasgow University for their help with the quantitative RT-PCR, Dr Catherine Deborde at INRA-Bordeaux for her help with metabolite identification, and Dr. David Broadhurst at the Manchester Interdisciplinary Biocentre for his advice on data analysis. This work was supported by the Biotechnology and Biological Sciences Research Council [grant number P19408]. We would also like to thank the Sir William Roberts Trust and the School of Biological Sciences for funding the studentship.
- Bewley JD: Seed germination and dormancy. Plant Cell. 1997, 9 (7): 1055-1066. 10.1105/tpc.9.7.1055PubMed CentralView ArticlePubMedGoogle Scholar
- Finch-Savage WE, Leubner-Metzger G: Seed dormancy and the control of germination. New Phytologist. 2006, 171 (3): 501-523.View ArticlePubMedGoogle Scholar
- Ogawa M, Hanada A, Yamauchi Y, Kuwahara A, Kamiya Y, Yamaguchi S: Gibberellin biosynthesis and response during Arabidopsis seed germination. Plant Cell. 2003, 15 (7): 1591-1604. 10.1105/tpc.011650PubMed CentralView ArticlePubMedGoogle Scholar
- Penfield S, Rylott EL, Gilday AD, Graham S, Larson TR, Graham IA: Reserve mobilization in the Arabidopsis endosperm fuels hypocotyl elongation in the dark, is independent of abscisic acid, and requires PHOSPHOENOLPYRUVATE CARBOXYKINASE1. Plant Cell. 2004, 16 (10): 2705-2718. 10.1105/tpc.104.024711PubMed CentralView ArticlePubMedGoogle Scholar
- Penfield S, Li Y, Gilday AD, Graham S, Graham IA: Arabidopsis ABA INSENSITIVE4 regulates lipid mobilization in the embryo and reveals repression of seed germination by the endosperm. Plant Cell. 2006, 18 (8): 1887-1899. 10.1105/tpc.106.041277PubMed CentralView ArticlePubMedGoogle Scholar
- Nakabayashi K, Okamoto M, Koshiba T, Kamiya Y, Nambara E: Genome-wide profiling of stored mRNA in Arabidopsis thaliana seed germination: epigenetic and genetic regulation of transcription in seed. Plant Journal. 2005, 41 (5): 697-709. 10.1111/j.1365-313X.2005.02337.xView ArticlePubMedGoogle Scholar
- Cadman CSC, Toorop PE, Hilhorst HWM, Finch-Savage WE: Gene expression profiles of Arabidopsis Cvi seeds during dormancy cycling indicate a common underlying dormancy control mechanism. Plant Journal. 2006, 46 (5): 805-822. 10.1111/j.1365-313X.2006.02738.xView ArticlePubMedGoogle Scholar
- Cao DN, Cheng H, Wu W, Soo HM, Peng JR: Gibberellin mobilizes distinct DELLA-dependent transcriptomes to regulate seed germination and floral development in arabidopsis. Plant Physiology. 2006, 142 (2): 509-525. 10.1104/pp.106.082289PubMed CentralView ArticlePubMedGoogle Scholar
- Carrera E, Holman T, Medhurst A, Peer W, Schmuths H, Footitt S, Theodoulou FL, Holdsworth MJ: Gene expression profiling reveals defined functions of the ATP-binding cassette transporter COMATOSE late in phase II of germination. Plant Physiology. 2007, 143 (4): 1669-1679. 10.1104/pp.107.096057PubMed CentralView ArticlePubMedGoogle Scholar
- Bassel GW, Fung P, Chow TFF, Foong JA, Provart NJ, Cutler SR: Elucidating the germination transcriptional program using small molecules. Plant Physiology. 2008, 147 (1): 143-155. 10.1104/pp.107.110841PubMed CentralView ArticlePubMedGoogle Scholar
- Fait A, Angelovici R, Less H, Ohad I, Urbanczyk-Wochniak E, Fernie AR, Galili G: Arabidopsis seed development and germination is associated with temporally distinct metabolic switches. Plant Physiology. 2006, 142 (3): 839-854. 10.1104/pp.106.086694PubMed CentralView ArticlePubMedGoogle Scholar
- Ghassemian M, Lutes J, Tepperman JM, Chang HS, Zhu T, Wang X, Quail PH, Lange BM: Integrative analysis of transcript and metabolite profiling data sets to evaluate the regulation of biochemical pathways during photomorphogenesis. Archives of Biochemistry and Biophysics. 2006, 448 (1-2): 45-59. 10.1016/j.abb.2005.11.020View ArticlePubMedGoogle Scholar
- Pritchard SL, Charlton WL, Baker A, Graham IA: Germination and storage reserve mobilization are regulated independently in Arabidopsis. Plant Journal. 2002, 31 (5): 639-647. 10.1046/j.1365-313X.2002.01376.xView ArticlePubMedGoogle Scholar
- Mansfield SG, Briarty LG: The dynamics of seedling and cotyledon cell development in Arabidopsis thaliana during reserve mobilization. International Journal of Plant Sciences. 1996, 157 (3): 280-295. 10.1086/297347.View ArticleGoogle Scholar
- Sheen J: Metabolic Repression of Transcription in Higher-Plants. Plant Cell. 1990, 2 (10): 1027-1038. 10.1105/tpc.2.10.1027PubMed CentralView ArticlePubMedGoogle Scholar
- Rolland F, Baena-Gonzalez E, Sheen J: SUGAR SENSING AND SIGNALING IN PLANTS: Conserved and Novel Mechanisms. 2006, 57: 675-709.Google Scholar
- Rubio V, Bustos R, Irigoyen ML, Cardona-Lopez X, Rojas-Triana M, Paz-Ares J: Plant hormones and nutrient signaling. Plant Molecular Biology. 2009, 69 (4): 361-373. 10.1007/s11103-008-9380-yView ArticlePubMedGoogle Scholar
- Leon P, Sheen J: Sugar and hormone connections. Trends in Plant Science. 2003, 8 (3): 110-116. 10.1016/S1360-1385(03)00011-6View ArticlePubMedGoogle Scholar
- Carrari F, Fernie AR, Iusem ND: Heard it through the grapevine? ABA and sugar cross-talk: the ASR story. Trends Plant Sci. 2004, 9 (2): 57-59. 10.1016/j.tplants.2003.12.004View ArticlePubMedGoogle Scholar
- Coruzzi GM, Zhou L: Carbon and nitrogen sensing and signaling in plants: emerging 'matrix effects'. Curr Opin Plant Biol. 2001, 4 (3): 247-253. 10.1016/S1369-5266(00)00168-0View ArticlePubMedGoogle Scholar
- Turner JE, Greville K, Murphy EC, Hooks MA: Characterization of Arabidopsis fluoroacetate-resistant mutants reveals the principal mechanism of acetate activation for entry into the glyoxylate cycle. The Journal of biological chemistry. 2005, 280 (4): 2780-2787. 10.1074/jbc.M407291200View ArticlePubMedGoogle Scholar
- Hooks MA, Turner JE, Murphy EC, Graham IA: Acetate non-utilizing mutants of Arabidopsis: evidence that organic acids influence carbohydrate perception in germinating seedlings. Mol Genet Genomics. 2004, 271 (3): 249-256. 10.1007/s00438-004-0985-9View ArticlePubMedGoogle Scholar
- Urbanczyk-Wochniak E, Luedemann A, Kopka J, Selbig J, Roessner-Tunali U, Willmitzer L, Fernie AR: Parallel analysis of transcript and metabolic profiles: a new approach in systems biology. EMBO reports. 2003, 4 (10): 989-993. 10.1038/sj.embor.embor944PubMed CentralView ArticlePubMedGoogle Scholar
- Saito K, Hirai MY, Yonekura-Sakakibara K: Decoding genes with coexpression networks and metabolomics - 'majority report by precogs'. Trends Plant Sci. 2008, 13 (1): 36-43. 10.1016/j.tplants.2007.10.006View ArticlePubMedGoogle Scholar
- 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 (27): 10205-10210. 10.1073/pnas.0403218101PubMed CentralView ArticlePubMedGoogle Scholar
- Hirai MY, Klein M, Fujikawa Y, Yano M, Goodenowe DB, Yamazaki Y, Kanaya S, Nakamura Y, Kitayama M, Suzuki H, et al.: Elucidation of gene-to-gene and metabolite-to-gene networks in arabidopsis by integration of metabolomics and transcriptomics. The Journal of biological chemistry. 2005, 280 (27): 25590-25595. 10.1074/jbc.M502332200View ArticlePubMedGoogle Scholar
- Nikiforova VJ, Daub CO, Hesse H, Willmitzer L, Hoefgen R: Integrative gene-metabolite network with implemented causality deciphers informational fluxes of sulphur stress response. J Exp Bot. 2005, 56 (417): 1887-1896. 10.1093/jxb/eri179View ArticlePubMedGoogle Scholar
- Szymanski J, Bielecka M, Carrari F, Fernie AR, Hoefgen R, Nikiforova VJ: On the processing of metabolic information through metabolite-gene communication networks: an approach for modelling causality. Phytochemistry. 2007, 68 (16-18): 2163-2175. 10.1016/j.phytochem.2007.04.017View ArticlePubMedGoogle Scholar
- Hoefgen R, Nikiforova VJ: Metabolomics integrated with transcriptomics: assessing systems response to sulfur-deficiency stress. Physiol Plant. 2008, 132 (2): 190-198.View ArticlePubMedGoogle Scholar
- Hirai MY, Sugiyama K, Sawada Y, Tohge T, Obayashi T, Suzuki A, Araki R, Sakurai N, Suzuki H, Aoki K, et al.: Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proc Natl Acad Sci USA. 2007, 104 (15): 6478-6483. 10.1073/pnas.0611629104PubMed CentralView ArticlePubMedGoogle Scholar
- Gutierrez RA, Stokes TL, Thum K, Xu X, Obertello M, Katari MS, Tanurdzic M, Dean A, Nero DC, McClung CR, et al.: Systems approach identifies an organic nitrogen-responsive gene network that is regulated by the master clock control gene CCA1. Proc Natl Acad Sci USA. 2008, 105 (12): 4939-4944. 10.1073/pnas.0800211105PubMed CentralView ArticlePubMedGoogle Scholar
- Thum KE, Shin MJ, Gutierrez RA, Mukherjee I, Katari MS, Nero D, Shasha D, Coruzzi GM: An integrated genetic, genomic and systems approach defines gene networks regulated by the interaction of light and carbon signaling pathways in Arabidopsis. BMC Syst Biol. 2008, 2: 31- 10.1186/1752-0509-2-31PubMed CentralView ArticlePubMedGoogle Scholar
- Carrari F, Fernie AR: Metabolic regulation underlying tomato fruit development. J Exp Bot. 2006, 57 (9): 1883-1897. 10.1093/jxb/erj020View ArticlePubMedGoogle Scholar
- Schauer N, Semel Y, Balbo I, Steinfath M, Repsilber D, Selbig J, Pleban T, Zamir D, Fernie AR: Mode of inheritance of primary metabolic traits in tomato. Plant Cell. 2008, 20 (3): 509-523. 10.1105/tpc.107.056523PubMed CentralView ArticlePubMedGoogle Scholar
- Mounet F, Moing A, Garcia V, Petit J, Maucourt M, Deborde C, Bernillon S, Le Gall G, Colquhoun I, Defernez M, et al.: Gene and Metabolite Regulatory Network Analysis of Early Developing Fruit Tissues Highlights New Candidate Genes for the Control of Tomato Fruit Composition and Development. Plant Physiology. 2009, 149 (3): 1505-1528. 10.1104/pp.108.133967PubMed CentralView ArticlePubMedGoogle Scholar
- Gibon Y, Usadel B, Blaesing OE, Kamlage B, Hoehne M, Trethewey R, Stitt M: Integration of metabolite with transcript and enzyme activity profiling during diurnal cycles in Arabidopsis rosettes. Genome Biol. 2006, 7 (8):
- Ebbels TMD, Buxton BF, Jones DT: springScape: visualisation of microarray and contextual bioinformatic data using spring embedding and an 'information landscape'. Bioinformatics. 2006, 22 (14): E99-E107. 10.1093/bioinformatics/btl205View ArticlePubMedGoogle Scholar
- Alba R, Payton P, Fei ZJ, McQuinn R, Debbie P, Martin GB, Tanksley SD, Giovannoni JJ: Transcriptome and selected metabolite analyses reveal multiple points of ethylene control during tomato fruit development. Plant Cell. 2005, 17 (11): 2954-2965. 10.1105/tpc.105.036053PubMed CentralView ArticlePubMedGoogle Scholar
- Rylott EL, Hooks MA, Graham IA: Co-ordinate regulation of genes involved in storage lipid mobilization in Arabidopsis thaliana. Biochemical Society Transactions. 2001, 29: 283-287. 10.1042/BST0290283View ArticlePubMedGoogle Scholar
- Biais B, Allwood JW, Deborde C, Xu Y, Maucourt M, Beauvoit B, Dunn WB, Jacob D, Goodacre R, Rolin D, et al.: 1H NMR, GC-EI-TOFMS, and data set correlation for fruit metabolomics: application to spatial metabolite analysis in melon. Anal Chem. 2009, 81 (8): 2884-2894. 10.1021/ac9001996View ArticlePubMedGoogle Scholar
- Crawford RMM: Tolerance of Anoxia and Ethanol-Metabolism in Germinating Seeds. New Phytologist. 1977, 79 (3): 511-517. 10.1111/j.1469-8137.1977.tb02235.x.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57 (1): 289-300.Google Scholar
- Liu PP, Montgomery TA, Fahlgren N, Kasschau KD, Nonogaki H, Carrington JC: Repression of AUXIN RESPONSE FACTOR10 by microRNA160 is critical for seed germination and post-germination stages. Plant J. 2007, 52 (1): 133-146. 10.1111/j.1365-313X.2007.03218.xView ArticlePubMedGoogle Scholar
- Loreti E: A genome-wide analysis of the effects of sucrose on gene expression in Arabidopsis seedlings under anoxia. Plant Physiology. 2005, 137 (3): 1130-1138. 10.1104/pp.104.057299PubMed CentralView ArticlePubMedGoogle Scholar
- Sachs MM, Subbaiah CC, Saab IN: Anaerobic gene expression and flooding tolerance in maize. Journal of Experimental Botany. 1996, 47 (294): 1-15. 10.1093/jxb/47.1.1.View ArticleGoogle Scholar
- Fiehn O, Kloska S, Altmann T: Integrated studies on plant biology using multiparallel techniques. Current Opinion in Biotechnology. 2001, 12 (1): 82-86. 10.1016/S0958-1669(00)00165-8View ArticlePubMedGoogle Scholar
- Sweetlove LJ, Fernie AR: Regulation of metabolic networks: understanding metabolic complexity in the systems biology era. New Phytologist. 2005, 168 (1): 9-23. 10.1111/j.1469-8137.2005.01513.xView ArticlePubMedGoogle Scholar
- Schauer N, Fernie AR: Plant metabolomics: towards biological function and mechanism. Trends in Plant Science. 2006, 11 (10): 508-516. 10.1016/j.tplants.2006.08.007View ArticlePubMedGoogle Scholar
- Yuan JS, Galbraith DW, Dai SY, Griffin P, Stewart CN: Plant systems biology comes of age. Trends in Plant Science. 2008, 13 (4): 165-171. 10.1016/j.tplants.2008.02.003View ArticlePubMedGoogle Scholar
- Eastmond PJ, Germain V, Lange PR, Bryce JH, Smith SM, Graham IA: Postgerminative growth and lipid catabolism in oilseeds lacking the glyoxylate cycle. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (10): 5669-5674. 10.1073/pnas.97.10.5669PubMed CentralView ArticlePubMedGoogle Scholar
- Lawand S, Dorne AJ, Long D, Coupland G, Mache R, Carol P: Arabidopsis a bout de souffle, which is homologous with mammalian carnitine acyl carrier, is required for postembryonic growth in the light. Plant Cell. 2002, 14 (9): 2161-2173. 10.1105/tpc.002485PubMed CentralView ArticlePubMedGoogle Scholar
- Krouk G, Tranchina D, Lejay L, Cruikshank AA, Shasha D, Coruzzi GM, Gutierrez RA: A Systems Approach Uncovers Restrictions for Signal Interactions Regulating Genome-wide Responses to Nutritional Cues in Arabidopsis. Plos Computational Biology. 2009, 5 (3):Google Scholar
- Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiology. 2004, 136 (1): 2621-2632. 10.1104/pp.104.046367PubMed CentralView ArticlePubMedGoogle Scholar
- Blasing OE, Gibon Y, Gunther M, Hohne M, Morcuende R, Osuna D, Thimm O, Usadel B, Scheible WR, Stitt M: Sugars and circadian regulation make major contributions to the global regulation of diurnal gene expression in Arabidopsis. Plant Cell. 2005, 17 (12): 3257-3281. 10.1105/tpc.105.035261PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzali S, Loreti E, Solfanelli C, Novi G, Alpi A, Perata P: Identification of sugar-modulated genes and evidence for in vivo sugar sensing in Arabidopsis. Journal of Plant Research. 2006, 119: 115-123. 10.1007/s10265-005-0251-1View ArticlePubMedGoogle Scholar
- Muller R, Morant M, Jarmer H, Nilsson L, Nielsen TH: Genome-wide analysis of the Arabidopsis leaf transcriptome reveals interaction of phosphate and sugar metabolism. Plant Physiology. 2007, 143 (1): 156-171. 10.1104/pp.106.090167PubMed CentralView ArticlePubMedGoogle Scholar
- Huq E, Al-Sady B, Hudson M, Kim CH, Apel M, Quail PH: PHYTOCHROME-INTERACTING FACTOR 1 is a critical bHLH regulator of chlorophyll biosynthesis. Science. 2004, 305 (5692): 1937-1941. 10.1126/science.1099728View ArticlePubMedGoogle Scholar
- Oh E, Yamaguchi S, Hu JH, Yusuke J, Jung B, Paik I, Lee HS, Sun TP, Kamiya Y, Choi G: PIL5, a phytochrome-interacting bHLH protein, regulates gibberellin responsiveness by binding directly to the GAI and RGA promoters in Arabidopsis seeds. Plant Cell. 2007, 19 (4): 1192-1208. 10.1105/tpc.107.050153PubMed CentralView ArticlePubMedGoogle Scholar
- Lopez-Molina L, Mongrand B, McLachlin DT, Chait BT, Chua NH: ABI5 acts downstream of ABI3 to execute an ABA-dependent growth arrest during germination. Plant Journal. 2002, 32 (3): 317-328. 10.1046/j.1365-313X.2002.01430.xView ArticlePubMedGoogle Scholar
- Penfield S, Josse EM, Kannangara R, Gilday AD, Halliday KJ, Graham IA: Cold and light control seed germination through the bHLH transcription factor SPATULA. Current Biology. 2005, 15 (22): 1998-2006. 10.1016/j.cub.2005.11.010View ArticlePubMedGoogle Scholar
- Yamagishi K, Tatematsu K, Yano R, Preston J, Kitamura S, Takahashi H, McCourt P, Kamiya Y, Nambara E: CHOTTO1, a Double AP2 Domain Protein of Arabidopsis thaliana, Regulates Germination and Seedling Growth Under Excess Supply of Glucose and Nitrate. Plant and Cell Physiology. 2009, 50 (2): 330-340. 10.1093/pcp/pcn201View ArticlePubMedGoogle Scholar
- Brady SM, Sarkar SF, Bonetta D, McCourt P: The ABSCISIC ACID INSENSITIVE 3 (ABI3) gene is modulated by farnesylation and is involved in auxin signaling and lateral root development in Arabidopsis. Plant J. 2003, 34 (1): 67-75. 10.1046/j.1365-313X.2003.01707.xView ArticlePubMedGoogle Scholar
- Fukaki H, Taniguchi N, Tasaka M: PICKLE is required for SOLITARY-ROOT/IAA14-mediated repression of ARF7 and ARF19 activity during Arabidopsis lateral root initiation. Plant Journal. 2006, 48 (3): 380-389. 10.1111/j.1365-313X.2006.02882.xView ArticlePubMedGoogle Scholar
- Roth S, Gmunder H, Droge W: Regulation of intracellular glutathione levels and lymphocyte functions by lactate. Cellular immunology. 1991, 136 (1): 95-104. 10.1016/0008-8749(91)90384-NView ArticlePubMedGoogle Scholar
- Walenta S, Mueller-Klieser WF: Lactate: Mirror and motor of tumor malignancy. Seminars in Radiation Oncology. 2004, 14 (3): 267-274. 10.1016/j.semradonc.2004.04.004View ArticlePubMedGoogle Scholar
- Schmid SA, Gaumann A, Wondrak M, Eckermann C, Schulte S, Mueller-Klieser W, Wheatley DN, Kunz-Schughart LA: Lactate adversely affects the in vitro formation of endothelial cell tubular structures through the action of TGF-beta 1. Experimental Cell Research. 2007, 313 (12): 2531-2549. 10.1016/j.yexcr.2007.05.016View ArticlePubMedGoogle Scholar
- Walenta S, Schroeder T, Mueller-Klieser W: Lactate in solid malignant tumors: Potential basis of a metabolic classification in clinical oncology. Current Medicinal Chemistry. 2004, 11 (16): 2195-2204.View ArticlePubMedGoogle Scholar
- Haberer G, Mader MT, Kosarev P, Spannagl M, Yang L, Mayer KF: Large-scale cis-element detection by analysis of correlated expression and sequence conservation between Arabidopsis and Brassica oleracea. Plant Physiol. 2006, 142 (4): 1589-1602. 10.1104/pp.106.085639PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang B, Singh KB: Ocs Element Promoter Sequences Are Activated by Auxin and Salicylic-Acid in Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91 (7): 2507-2511. 10.1073/pnas.91.7.2507PubMed CentralView ArticlePubMedGoogle Scholar
- Chen W, Singh KB: The auxin, hydrogen peroxide and salicylic acid induced expression of the Arabidopsis GST6 promoter is mediated in part by an ocs element. Plant J. 1999, 19 (6): 667-677. 10.1046/j.1365-313x.1999.00560.xView ArticlePubMedGoogle Scholar
- Smeekens S: Sugar-induced signal transduction in plants. Annual Review of Plant Physiology and Plant Molecular Biology. 2000, 51: 49-81. 10.1146/annurev.arplant.51.1.49View ArticlePubMedGoogle Scholar
- Rolland F, Sheen J: Sugar sensing and signalling networks in plants. Biochemical Society Transactions. 2005, 33: 269-271. 10.1042/BST0330269View ArticlePubMedGoogle Scholar
- Murashige T, Skoog F: A revised medium for rapid growth and bioassay with tobacco tissue cultures. Physiologia Plantarum. 1962, 15: 473-496. 10.1111/j.1399-3054.1962.tb08052.x.View ArticleGoogle Scholar
- Armengaud P, Breitling R, Amtmann A: The potassium-dependent transcriptome of Arabidopsis reveals a prominent role of jasmonic acid in nutrient signaling. Plant Physiology. 2004, 136 (1): 2556-2576. 10.1104/pp.104.046482PubMed CentralView ArticlePubMedGoogle Scholar
- Kerr MK, Churchill GA: Bootstrapping cluster analysis: Assessing the reliability of conclusions from microarray experiments. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98 (16): 8961-8965. 10.1073/pnas.161273698PubMed CentralView ArticlePubMedGoogle Scholar
- Love AJ, Laval V, Geri C, Laird J, Tomos AD, Hooks MA, Milner JJ: Components of Arabidopsis defense- and ethylene-signaling pathways regulate susceptibility to Cauliflower mosaic virus by restricting long-distance movement. Mol Plant-Microbe Interact. 2007, 20 (6): 659-670.View ArticlePubMedGoogle Scholar
- Weckwerth W, Wenzel K, Fiehn O: Process for the integrated extraction identification, and quantification of metabolites, proteins and RNA to reveal their co-regulation in biochemical networks. Proteomics. 2004, 4 (1): 78-83. 10.1002/pmic.200200500View ArticlePubMedGoogle Scholar
- Moing A, Maucourt M, Renaud C, Gaudillere M, Brouquisse R, Lebouteiller B, Gousset-Dupont A, Vidal J, Granot D, Denoyes-Rothan B, et al.: Quantitative metabolic pro. ling by 1-dimensional H-1-NMR analyses: application to plant genetics and functional genomics. Functional Plant Biology. 2004, 31 (9): 889-902. 10.1071/FP04066.View ArticleGoogle Scholar
- Akoka S, Barantin L, Trierweiler M: Concentration measurement by proton NMR using the ERETIC method. Analytical Chemistry. 1999, 71 (13): 2554-2557. 10.1021/ac981422i.View ArticlePubMedGoogle Scholar
- Ward JL, Harris C, Lewis J, Beale MH: Assessment of H-1 NMR spectroscopy and multivariate analysis as a technique for metabolite fingerprinting of Arabidopsis thaliana. Phytochemistry. 2003, 62 (6): 949-957. 10.1016/S0031-9422(02)00705-7View ArticlePubMedGoogle Scholar
- Le Gall G, Metzdorff SB, Pedersen J, Bennett RN, Colquhoun IJ: Metabolite profiling of Arabidopsis thaliana (L.) plants transformed with an antisense chalcone synthase gene. Metabolomics. 2005, 1 (2): 181-198. 10.1007/s11306-005-4434-5.View ArticleGoogle Scholar