Open Access Highly Accessed Research article

A transcriptional analysis of carotenoid, chlorophyll and plastidial isoprenoid biosynthesis genes during development and osmotic stress responses in Arabidopsis thaliana

Stuart Meier1, Oren Tzfadia23, Ratnakar Vallabhaneni23, Chris Gehring14 and Eleanore T Wurtzel23*

Author Affiliations

1 Division of Chemistry, Life Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia

2 Department of Biological Sciences, Lehman College, The City University of New York, 250 Bedford Park Blvd. West, Bronx, NY 10468, USA

3 The Graduate School and University Center-CUNY, 365 Fifth Ave., New York, NY 10016-4309, USA

4 Department of Biotechnology, University of the Western Cape, Private Bag X17, Cape Town - Bellville 7535, South Africa

For all author emails, please log on.

BMC Systems Biology 2011, 5:77  doi:10.1186/1752-0509-5-77

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/5/77


Received:14 December 2010
Accepted:19 May 2011
Published:19 May 2011

© 2011 Meier et al; licensee BioMed Central Ltd.

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

Abstract

Background

The carotenoids are pure isoprenoids that are essential components of the photosynthetic apparatus and are coordinately synthesized with chlorophylls in chloroplasts. However, little is known about the mechanisms that regulate carotenoid biosynthesis or the mechanisms that coordinate this synthesis with that of chlorophylls and other plastidial synthesized isoprenoid-derived compounds, including quinones, gibberellic acid and abscisic acid. Here, a comprehensive transcriptional analysis of individual carotenoid and isoprenoid-related biosynthesis pathway genes was performed in order to elucidate the role of transcriptional regulation in the coordinated synthesis of these compounds and to identify regulatory components that may mediate this process in Arabidopsis thaliana.

Results

A global microarray expression correlation analysis revealed that the phytoene synthase gene, which encodes the first dedicated and rate-limiting enzyme of carotenogenesis, is highly co-expressed with many photosynthesis-related genes including many isoprenoid-related biosynthesis pathway genes. Chemical and mutant analysis revealed that induction of the co-expressed genes following germination was dependent on gibberellic acid and brassinosteroids (BR) but was inhibited by abscisic acid (ABA). Mutant analyses further revealed that expression of many of the genes is suppressed in dark grown plants by Phytochrome Interacting transcription Factors (PIFs) and activated by photoactivated phytochromes, which in turn degrade PIFs and mediate a coordinated induction of the genes. The promoters of PSY and the co-expressed genes were found to contain an enrichment in putative BR-auxin response elements and G-boxes, which bind PIFs, further supporting a role for BRs and PIFs in regulating expression of the genes. In osmotically stressed root tissue, transcription of Calvin cycle, methylerythritol 4-phosphate pathway and carotenoid biosynthesis genes is induced and uncoupled from that of chlorophyll biosynthesis genes in a manner that is consistent with the increased synthesis of carotenoid precursors for ABA biosynthesis. In all tissues examined, induction of β-carotene hydroxylase transcript levels are linked to an increased demand for ABA.

Conclusions

This analysis provides compelling evidence to suggest that coordinated transcriptional regulation of isoprenoid-related biosynthesis pathway genes plays a major role in coordinating the synthesis of functionally related chloroplast localized isoprenoid-derived compounds.

Background

The carotenoids are pure isoprenoids that are synthesized in chloroplasts from geranylgeranyl diphosphate (GGPP) which additionally serves as an immediate precursor for other chloroplastic localized isoprenoid biosynthesis pathways including plastoquinone (PQ), the phytol tail of chlorophylls, phylloquinones (PhQ) and tocopherols as well as the phytohormone gibberellic acid (GA). While the biochemistry of carotenoid biosynthesis (CrtBS) has been extensively studied and most genes encoding enzymes that function in the CrtBS pathway have been identified, little is known about how the synthesis of these enzymes is coordinated and additionally how this synthesis is coordinated with that of other interdependent and interrelated isoprenoid-derived compounds. We have performed a global in-silico expression correlation analysis using microarray experimental data to identify genes that share a high level of co-expression and thus may share closely associated functional relationships with phytoene synthase (PSY). Comprehensive expression profiling of chloroplastic isoprenoid-related biosynthesis pathway genes was performed over a range of developmental and stress-related conditions in order to identify important regulatory components such as phytohormones and transcriptional regulatory factors that are important in coordinating their collective expression.

A number of chloroplast localized isoprenoid-derived compounds constitute important components of the photosynthetic apparatus. The carotenoids perform a range of functions including the acquisition of light energy and photoprotection [1] and additionally serve as precursors for abcisic acid (ABA) biosynthesis [2]. The chlorophylls are the main light absorbing pigments of the photosynthetic apparatus while PhQ and PQ function in photosynthetic electron transfer reactions. Plastoquinone additionally functions as an essential electron carrier in CrtBS desaturation reactions mediated by phytoene desaturase (PDS) and ζ-carotene desaturase (ZDS) [3]. As the GGPP molecule is an immediate precursor for the biosynthesis of these functionally related molecules, it serves as an important metabolic hub in the biosynthesis of essential components of the photosynthetic apparatus (see Figure 1, dark text).

thumbnailFigure 1. Diagram of the plastidial isoprenoid biosynthesis pathway detailing the level of co-expression that PSY shares with individual interrelated isoprenoid biosynthesis pathway genes. The pathways represented include the Calvin Cycle, MEP, Carotenoid, Chlorophyll, Phylloqinone, Plastoquinone, ABA and Gibberellins and are collectively referred to as the PSY-correlated isoprenoid interrelated genes (PSY-CIIG). Reaction substrates and products are represented in bold black letters while genes that encode pathway enzymes are in black italic letters. Numbers in red parentheses represent expression correlation r-values and numbers in square brackets indicate the number of paralog genes that are annotated to encode the respective enzymes. All r-values > 0.5 had p-values and e-values <1.0-15. Non-significant r-values are indicated as n.s. Only the highest correlated member of the paralog gene family and those that have a co-expression value > 0.6 are listed. See Additional File 1 for list of corresponding gene IDs, details of statistics for individual genes and an extended list including additional paralog gene family members.

The synthesis of GGPP in plastids starts from pyruvate and glyceraldehyde 3-phosphate (GAP), that can be generated directly from the Calvin cycle (photosynthesis) or glycolysis [4], and serve as precursor molecules for the methylerythritol 4-phosphate (MEP) pathway [5]. The MEP pathway consists of a series of seven enzymes that function sequentially to catalyze the synthesis of the prenyl diphosphate precursors, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). GGPP synthase (GGPPS) then catalyzes the sequential addition of three molecules of IPP to one molecule of DMAPP resulting in the formation of the poly-isoprenoid, GGPP [4].

The PSY gene encodes the first dedicated enzyme of the CrtBS pathway which catalyzes the condensation of two molecules GGPP to form 15-cis-phytoene (Figure 1) [6-8]. Phytoene is metabolized to lycopene in a linear series of desaturation and isomerisation reactions that involves four enzymes [9,10]. CrtBS then branches into two distinct pathways (Figure 1), the β-cyclase (LCY-β, also known as LCYB or LYC-β) enzyme converts lycopene into β-carotene while the dual action of epsilon cyclase (LCY-ε, also known as LCYE or LYC-ε) and LCY-β results in the formation of α-carotene. The α- and β- carotenes can then be hydroxylated to form α- and β- branch xanthophylls which are essential components of the photosynthetic apparatus in higher plants where they function in photosystem assembly, light harvesting and photoprotection [11,12]. In addition, violaxanthin and neoxanthin also serve as direct precursors for ABA biosynthesis and are alternative substrates for plastid localized nine cis-epoxycarotenoid dioxygenases (NCEDs) [2,13] (Figure 1).

The carotenoids have been shown to have important functional roles in early stages of post-germination development. In etioplasts of dark grown seedlings, lutein and violaxanthin biosynthesis is required for assembly of the prolamellar body (PLB) [14,15]; a lattice of tubular membranes composed primarily of lipids, carotenoids and a ternary complex of NADPH, protochlorophyllide oxidoreductase (POR) and the chlorophyll precursor, protochlorophyllide (Pchlide) [16]. The synthesis of carotenoids in PLBs is thought to optimize the transition of etiolated plants to photomorphogenic development since it has been shown to enhance chlorophyll accumulation and greening upon light-induced de-etiolation [14,17].

Light activates the differentiation of etioplasts into chloroplasts in a process that is accompanied by a large and coordinated increase in the biosynthesis and accumulation of carotenoids, chlorophylls and pigment-binding proteins; this accumulation supports the development of a functional photosynthetic apparatus [14,15,18]. The light-induced synthesis of carotenoids is characterized by an increase in expression of PSY and select MEP pathway genes [18,19] as well as an increase in PSY enzymatic activity [15]. The expression of PSY has been shown to be elevated in response to a broad spectrum of continuous (c) light wavelengths including far-red (cFR), red (cR), blue (cB) and white (cW) [18,20,21]. While PSY transcript levels have been reported to increase in response to cFR [18], only light wavelengths that activate POR - which catalyzes the light-dependent conversion of Pchlide to chlorophyllide - cause the decay of PLBs, the synthesis of chlorophylls and the transition of etioplasts into chloroplasts [15,22]. These studies demonstrate that the coordinated and co-localized synthesis of carotenoids with chlorophyll precursors and chlorophylls in etiolated and de-etiolated plants respectively is required for normal photomorphogenic development.

The light-induced increase in PSY expression has been shown to be mediated by the phytochrome (PHYs) photoreceptors. Mutant studies have shown that the induction of PSY expression in response to cFR is dependent on the light-labile PHY-A while the cR-induction is thought to be mediated by light-stable PHYs other than PHY-B [18]. Upon light-induced activation, the cytoplasmic localized PHYs are translocated to the nucleus where they interact with and mediate the degradation of the Phytochrome Interacting transcription Factors (PIFs); these factors bind to G-boxes in the promoters of light-induced genes and negatively regulate their expression [23]. Recently, the PIFs have been shown to have an important role in regulating the transcription of PSY and other carotenoid and chlorophyll biosynthesis genes during light-induced de-etiolation [24]. The PHYs and PIFs are therefore interesting candidate regulatory factors that may function to coordinate the transcription of genes that encode enzymes that function in the interrelated and interdependent chloroplastic isoprenoid biosynthesis pathways during early development.

Transcriptional co-regulation has been shown to play a major role in coordinating cellular responses that involve multiple genes and their products. A number of studies have shown that genes that have been confirmed to be co-expressed in response to a range of conditions have correlated functional relationships, including physical interactions between their encoded proteins [25-28]. These finding also extend to metabolic pathways where it has been shown that many genes encoding metabolic enzymes that function within the same or functionally related pathways form co-expression modules [29,30]. Thus, it is conceivable, that the synthesis of functionally related chloroplast localized isoprenoid molecules is mediated by their transcriptional co-regulation.

The model plant species Arabidopsis thaliana is ideal for studying global transcriptional responses since there are thousands of publicly available full-genome microarray experiments that encompass a broad range of experimental conditions including different developmental stages, stress, chemical and hormone treatments and mutants. In addition, analysis tools are available to identify modules of co-expressing genes and genome sequence data allows analysis of promoter regulatory regions and the identification of putative regulatory elements.

It is pertinent to acknowledge that changes in gene transcription do not necessarily translate to changes in protein abundance and functional activity due to post-transcriptional regulatory mechanisms. However, as these mechanisms rely on a gene being transcribed in the first instance, gene transcription can be considered the primary level of regulation of protein synthesis. While cells can alter the activity of specific proteins/enzymes to fine tune cellular responses, the protein must be synthesized and present at appropriate quantities for this to occur. Changes in gene transcription in response to specific stimuli can be considered a primary regulatory response that reflects a change in requirement for a specific protein(s) at a specific point in time. In addition, in comparison to a single gene, when the expression of a large group of functionally related genes is altered in a uniform manner in response to a specific stimulus, it is a stronger indicator that the transcriptional response is representative of a cell's intent to change the associated functional activity in response to the stimuli.

Here we aim to elucidate the role transcriptional regulation plays in coordinating CrtBS and the synthesis of other functionally related isoprenoid-derived compounds during early development and in response to osmotic stress. A global co-expression analysis revealed that PSY is highly co-expressed with many photosynthesis-related genes including, those involved in chlorophyll, PQ and PhQ biosynthesis as well as genes that function in the upstream Calvin cycle and MEP pathway that synthesize the commonly required GGPP precursor. Stimulus-specific transcription profiling revealed that expression of the isoprenoid biosynthesis genes is almost universally activated following germination, during both etiolated and de-etiolated growth and the induction during early development is positively regulated by BRs and GA and inhibited by ABA. During etiolated growth, the PIFs appear to suppress the expression of the genes while PHYs mediate their photoactivation. An enrichment in putative BR-auxin response elements and G-boxes (which bind PIFs) in the promoter of PSY and the co-expressed genes further supports a role for BRs and PIFs in regulating expression of the genes. In osmotically stressed root tissue, transcription of CrtBS-related genes is induced in a manner that is consistent with the increased synthesis of carotenoid precursors for ABA biosynthesis. In all tissues examined, induction of β-carotene hydroxylase transcript levels is linked to increased demand for ABA. We therefore conclude that transcriptional regulation plays a major role in coordinating the synthesis of functionally related isoprenoid-derived compounds in chloroplasts.

Results and Discussion

PSY co-expression analysis

In order to elucidate the role transcriptional regulation plays in coordinating CrtBS and the synthesis of other functionally related isoprenoid-derived compounds; an expression correlation analysis was undertaken using PSY as the driver gene in order to determine the level of co-expression that PSY shares with all of the other genes represented on the ATH1 microarray (22 K) chip. Key to the accuracy of this analysis is that co-expression is measured over a large number of diverse experimental conditions (see methods and Ref [31]). PSY was selected as the driver gene for this analysis as it is the first dedicated enzyme of carotenogenesis and its transcription is known to be positively correlated with and a major driving force for carotenoid production [9,14,20,32]. It is thought that genes that are highly co-expressed with PSY will have closely associated functional roles.

The expression of PSY was shown to be highly correlated with many genes in the genome with the top 50 expression correlated genes (PSY-ECG50) having a Pearson correlation coefficient (r-value) ranging from 0.91 to 0.84 (Table 1). In total, approximately 1000 genes (4.3%) had an r-value > 0.6 while around 600 (2.6%) had an r-value > 0.7 supporting the specificity of the analysis since it indicates that PSY is co-expressed with only small percentage of select genes in the Arabidopsis thaliana genome. All genes in the PSY-ECG50 had highly significant p-values (<1-35) and e-values (<1-35) supporting the biological significance of the results.

Table 1. List of the 50 genes that are most highly co-expressed with PSY

Functional enrichment analysis of the PSY-ECG50

The high expression correlation of the PSY-ECG50 is a strong indicator that these genes may function in common biological processes. The PSY-ECG50 was therefore subjected to a functional enrichment analysis using "Fatigoplus" [33] which identified a number of significant enrichments in functional terms associated with the group (Table 1). In the biological process category, significant enrichments are found with genes associated with the terms photosynthesis, plastid organization and biogenesis, PQ biosynthetic process, and carotenoid and tetraterpenoid metabolic processes. In the cellular component category at level nine, genes associated with the terms plastid parts, thylakoid parts and chloroplasts are enriched.

Specifically, a number of genes in the PSY-ECG50 encode enzymes that directly function in the synthesis of chloroplastic localized isoprenoids. This includes ZDS (At3G04870) and LCY-β (At3G10230; r = 0.86 for both; Table 1, Figure 1) which function in the CrtBS pathway and the PHYTOENE DESATURATION 2 (PDS2, At3g11950, r = 0.86) [3,34,35] and ALBINO OR PALE GREEN MUTANT 1 (APG1; AT3G63410, r = 0.86) [36] genes that both function in the PQ biosynthesis pathway. In addition to its function as an electron carrier in PSII light-dependent photosynthesis reactions, PQ is also an essential compound in the synthesis of carotenoids where it has a role as a hydrogen acceptor in the desaturation reactions mediated by PDS and ZDS [3]. The GLCYERALDEHYDE-3-PHOSPHATE DEHYDROGENASE B SUBUNIT (GAPDβ) gene (At1g42970, r = 0.89) functions in the Calvin cycle to synthesize GAP which is a direct substrate for the MEP pathway [5,37].

In general, PSY is co-expressed with genes that encode proteins that have critical functional roles in the photosynthetic machinery; these proteins include enzymes that function in the biosynthesis of carotenoids, chlorophylls and components of the photosynthetic electron transport chain including PQ, PhQ, plastidial NAD(P)H dehydrogenase complex, thioredoxin, ferredoxin, plastocyanin and the cytochrome b6/f complex as well as proteins that form structural components of photosystem I and II. The high co-expression of PSY with genes that encode proteins that have important functional roles in the photosynthetic machinery, including a number of isoprenoid biosynthesis genes, illustrates that PSY is indeed highly co-expressed with functionally related genes and this gives confidence in the accuracy of the analysis.

PSY co-expression with functionally related isoprenoid biosynthesis genes

The expression correlation values were next extracted for all known Arabidopsis genes that encode enzymes that function in plastidial isoprenoid biosynthesis; this included Calvin cycle and MEP pathway genes as well as carotenoid, chlorophyll, PQ, PhQ, ABA, and GA biosynthesis genes (Figure 1 - see Additional File 1 for full list of genes). These genes will collectively be referred to as the PSY-correlated interrelated isoprenoid biosynthesis genes (PSY-CIIG).

Additional file 1. Additional Table 1. Extended list of genes in the PSY-correlated interrelated isoprenoid biosynthesis genes and their expression correlation relative to PSY.

Format: XLS Size: 84KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

This analysis revealed that the expression of all nuclear genes that encode enzymes that are known or predicted to function at each of the individual steps in the CrtBS pathway are highly correlated with PSY (Figure 1). In addition, PSY is also highly co-expressed with many isoprenoid-related biosynthesis pathway genes including, Calvin cycle and MEP pathway genes as well as chlorophyll, PQs and PhQs biosynthesis genes (Figure 1). It is noteworthy that the expression of PSY was found to be highly correlated with genes that function in different branches of chlorophyll biosynthesis; this includes chlorophyll synthetase (ChlSyn, At3g51820, r = 0.82) that functions in phytol side chain biosynthesis, as well as glutamyl tRNA reductase (GluTR/HEMA1, r = 0.77, e-value <1-35) and glutamate 1-semialdehyde aminotransferase (GSA2, r = 0.72, e-value <1-35) that function in the upstream tetrapyrrole branch of chlorophyll biosynthesis. Significantly, the GluTR and GSA2 genes encode enzymes that catalyze the biosynthesis of 5-aminolevulinic acid (ALA) which is the rate-limiting step for this pathway [38,39]. The high degree of co-expression of these genes strongly suggests that their transcription is regulated by a common mechanism. In contrast, none of the ABA biosynthesis genes that operate downstream of ABA1/ZEP, or any GA biosynthesis genes are positively expression correlated with PSY (Figure 1).

This analysis also revealed that for the carotenoid and chlorophyll biosynthesis-related pathway enzymes that are encoded by multiple genes, only specific family members displayed high co-expression levels; this may imply their functional importance in their respective biosynthesis pathways (Figure 1 and Additional File 1). In the MEP pathway, 1-deoxy-D-xylulose 5-phosphate synthase (DXPS) is the only enzyme that is encoded for by multiple (three) nuclear genes in Arabidopsis [40] and of these, only the functionally determined DXPS2 (At4g15560, r = 0.69, e-value <1-35) displays a high level of co-expression with PSY [41]. The two Arabidopsis IPP isomerase (IPPI) genes, show very little correlation with PSY (Additional File 1) and this is consistent with, and lends support to a recent study that reported that these enzymes have minor functional roles in plastidial isoprenoid biosynthesis, since IPPI and DMAPP are directly synthesized by the MEP pathway in plastids [42].

Of the family of 12 annotated GGPPS genes in Arabidopsis [41], only GGPPS1 (At4g36810; r = 0.74, e-value <1-35), that encodes a functionally active and plastid localized enzyme [43], displays a high level of expression correlation. The GGPPS-like protein, geranylgeranyl reductase (GGR, At4g38460, r = 0.64, e-value = 8.6-35) also shows some expression correlation, however, GGR does not have GGPPS activity in vitro and its function remains unknown although it has been suggested to encode a GPPS subunit [43].

Since GGPP is a common substrate, and thus an important metabolic link in the synthesis of multiple isoprenoid-derived compounds, it is significant that in addition to PSY, GGPPS1 is the most highly co-expressed GGPPS in Arabidopsis with a number of other genes that encode chloroplast localized enzymes that directly use GGPP as a substrate. These genes include GGR2 (At1g74470, GGR2 to GGPPS1, r = 0.53 (e-value = 2.2-20), data not shown) and two genes that encode solanesyl diphosphate synthase enzymes, SPS-1 (At1g78510, SPS1 to GGPPS1, r = 0.51 (e-value = 2.7-18), data not shown) and SPS-2 (At1g17050, SPS2 to GGPPS1, r = 0.58, (e-value = 6.8-26), data not shown). The GGR2 enzyme reduces GGPP to phytyl pyrophosphate [44] which forms essential phytol side chains for both chlorophyll and PhQ (vitamin K1) biosynthesis while the chloroplast localized SPS-2 enzyme catalyzes the synthesis of solanesyl diphosphate (SPP) which is thought to be a precursor of the PQ side-chain in Arabidopsis [45,46]. In addition, PSY is also highly co-expressed with all enzymes that function downstream in these pathways including, as mentioned ChlSyn for chlorophyll biosynthesis, C-methyltransferase for PhQ synthesis (At1g23360, r = 0.74) and, as mentioned above, APG1 and PDS2 for PQ synthesis (Figure 1).

A co-correlation scatterplot between PSY and GGR2 (Figure 2) illustrates that both genes have a high level of co-expression with GGPPS1 and many genes that function in the chlorophyll, PQ and PhQ biosynthesis pathways thus providing strong evidence that the GGPPS1 enzyme plays a major role in generating a common pool of GGPP substrate that is used in the biosynthesis of these compounds. This interpretation is supported by a recent study that shows a reduction in carotenoid and chloroplast levels in a ggpps1 knock out mutant [47] and suggests that transcriptional regulation of the GGPPS1 gene serves as an important regulatory node in coordinating carotenoid, chlorophyll, PhQ and PQ biosynthesis.

thumbnailFigure 2. Co-correlation scatter plot illustrating the level of co-expression of all Arabidopsis genes relative to PSY and GGR2. Genes that function in defined biosynthesis pathways (PSY-CIIG ) are color highlighted as indicated in the legend. Select individual genes of interest are highlighted. All genes listed in Figure 1 are represented.

The scatterplot also shows that expression of carotenoid and chlorophyll biosynthesis genes is not correlated with any of the plastid localized GA biosynthesis genes, and in particular, GGPPS1 showed no co-expression with GA1(r = -0.26, e-value = 0.04) or GA2 (r = - 0.03, e-value > 1), which directly use plastidial GGPP pools as substrates for GA synthesis. This implies that GGPPS1 does not function in the synthesis of GGPP for GA biosynthesis and this in turn is supported by reports that state that ggpps1 mutants show no signs of GA deficiency [47] and that the generation of GGPP pools for GA biosynthesis is dependent on the action of GPPS (At2g34630) [48].

Of the carotenoid cleavage dioxygenase gene (CCD) family, only CCD1 (At3g63520, r = 0.76, e-value <1-35), and to a lesser degree CCD4 (At4g19170, r = 0.51, e-value = 1.4-18), showed any degree of expression correlation with PSY (Additional File 1). In Arabidopsis, the CCD1 protein is localized to the cytoplasm and thus inaccessible to plastid localized carotenoid substrates [49]. It has been proposed to function in the metabolism of carotenoids that are localized in the chloroplast envelope [50] or present in dry seeds which lack well defined organelles such as chloroplasts [49]. Significantly, PSY is not co-expressed with any of the five plastid localized ABA biosynthesis NCED genes (-2, -3, -5, -6) that catalyze the synthesis of xanthoxin from β,β-xanthophylls [50,51], or other downstream enzymes that function in ABA biosynthesis (Figure 2 and 3), indicating that transcription of the CrtBS genes is not directly coupled to ABA biosynthesis. Indeed, previous reports have shown that ABA biosynthesis is correlated with the expression of the NCED genes suggesting that their expression is important in regulating ABA-biosynthesis [2].

Stimulus specific expression analysis

The expression correlation analysis provides a generalized measure of how PSY is co-expressed with other genes in the genome since it is performed across multiple tissues and in response to a broad range of experimental conditions. In order to determine if expression of the genes corresponds to the known timing of carotenoid and chlorophyll biosynthesis and development of the photosynthetic apparatus, the expression of individual isoprenoid biosynthesis genes and genes in the PSY-ECG50 were examined throughout key developmental stages that are known to involve coordinated changes in the synthesis of carotenoids and chlorophylls.

Many previous studies on CrtBS gene expression have focused primarily on the transcriptional responses of a small subset of genes during de-etiolation since this is when large increases in CrtBS occur concomitant with the development of a functional photosynthetic apparatus [15,18,19,52]. However, carotenoids also have important functional roles during seed development, maturation and germination since they serve as precursors for ABA biosynthesis in developing seeds and CrtBS in dark-grown seedlings has been shown to be essential for PLB formation in etioplasts [14,17]. Thus, expression of the PSY-CIIG was examined throughout developmental stages encompassing seed development and maturation, imbibition and germination as well as etiolated and de-etiolated growth. Inhibitor and mutant experiments were also examined in order to determine the role that the early developmental-related phytohormones, ABA, GA and BR have in regulating the expression of these genes. Details of the experimental conditions for the microarray data examined are provided in Additional File 2.

Additional file 2. Additional Text 1. Description of experimental conditions in the microarray data-sets examined.

Format: DOC Size: 34KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Isoprenoid gene expression during seed and seedling development

The results in Figure 3 are presented as signal values since this provides information regarding the relative expression levels of individual genes at specific developmental stages and can provide insights into genes that may be rate-limiting due to low expression levels. While fold change ratios identify changes in expression in response to different conditions, they do not provide information of the relative abundance of transcripts.

thumbnailFigure 3. Expression heat map illustrating the relative expression levels of the PSY-CIIG during early developmental stages. The experimental conditions examined are listed across the top of the heat map and include, developing seed (GSE5634); stratified seed, etiolated and de-etiolated growth (GSE17159); imbibition in wild-type (wt), aba2 (ABA-deficient), and the cyp707a1,a2,a3 triple mutant (elevated ABA levels, GSE15700); post-germinative growth (24h) in the presence and absence of exogenously applied ABA or PAC (GSE5751). Individual genes included in the analysis are listed on the right and are arranged in sequential pathway order. Arrows indicate branch points where reaction products are used in multiple biosynthesis pathways. Results are presented as normalized signal values to reveal the relative expression levels of individual genes at conditions examined. Genes that function in biosynthesis pathway represented in Figure 1 are analysed. Details of the microarray experimental conditions are presented in Additional File 2 (Supporting Information).

The heat maps generated from the microarray expression analysis revealed that transcription of the PSY-ECG50 is modulated in a largely uniform manner in response to a range of different experimental conditions (Additional File 3) which is consistent with the high expression correlation of these genes. In general, the expression of PSY, other select CrtBS genes, chlorophyll, PQ and PlQ biosynthesis genes decline progressively throughout seed development, remained very low in dry and stratified seeds before being induced during imbibition, germination, skotomorphogenic and photomorphogenic growth (Figure 3). In addition, most genes that function in the upstream biosynthesis pathways to synthesize the commonly required GGPP precursor also follow a very similar expression profile, including the Calvin cycle GAPD subunit encoding genes (GAPD -β, A-1 and A-2); the MEP pathway genes DXPS2, MCT and GGPPS1. The very low expression level of specific genes in mature (stage 10), stratified and dry seeds may translate to very low protein/enzyme levels which would be rate-limiting for the biosynthesis of their respective molecules. While the activity of enzymes may be modulated to fine-tune biosynthesis rates, if there is no enzyme present then transcription of the gene will become the rate-limiting factor.

Additional file 3. Additional Figure 1. Heatmaps illustrating the expression of the PSY-ECG50 in response to the range of experimental conditions examined.

Format: PPT Size: 289KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

The observed induction in PSY expression in two day dark grown plants and in response to cR (Figure 3) in the absence of increased expression of PDS and ZDS is consistent with previous reports which additionally show that this response is sufficient to activate carotenoid biosynthesis [14,18,20]. The absence of a requirement for PDS and ZDS induction at these stages may be explained by the relatively high expression levels of these genes in mature (stage 10) and stratified seeds, when compared to PSY, which will presumably maintain relatively high enzyme levels. The activity of the PDS and ZDS enzymes however may be regulated by the abundance of the PDS2 enzyme; PDS2 functions in the biosynthesis of PQ [3,34] which is essential for the desaturation reactions mediated by PDS and ZDS [53]. Since the transcription of PDS2 closely mirrors that of PSY throughout these stages, it may be important in regulating the activities of the PDS and ZDS enzymes.

The collective moderate induction of these genes in dark grown plants and the strong induction in response to light coincide with the timing of carotenoid and chlorophyll biosynthesis during these developmental stages. As previously stated, carotenoids, the chlorophyll precursor Pchlide and the light-activated POR enzyme have been shown to accumulate in PLB of dark grown plants and carotenoid and chlorophyll biosynthesis is strongly activated in response to light [16]. This illustrates a correlation between biosynthesis pathway gene expression and carotenoid and chlorophyll biosynthesis at these developmental stages which implies that the transcriptional coordination of these pathway genes plays a major role in coordinating the synthesis of carotenoids and chlorophylls. Further, since PSY expression has been shown to be rate-determining and a major driving force for carotenoid biosynthesis [9,14], the closely coupled expression of other carotenoid and chlorophyll biosynthesis genes during these developmental stages suggests that expression of these genes may also be important in regulating and coordinating the biosynthesis of carotenoids and chlorophylls.

CrtBS gene expression and ABA biosynthesis in developing seeds

The expression profile of a number of CrtBS genes, including, βCHY1 and -2 and ZDS/ABA1, is in stark contrast with that of PSY during seed development in that their transcript levels remain elevated or increase during seed maturation, remain high in dry seeds and sharply decline during imbibition and in dark grown seedlings. These expression profiles are strikingly similar to a number of ABA biosynthesis genes, including NCED5, -6 and -9 that function directly downstream of ZEP/ABA1. Indeed, the induction of the ABA-responsive genes EM6 and RD29B [54] at latter stages of seed development (stage 7) strongly supports that an increase in endogenous ABA biosynthesis and accumulation occurs at this stage [55]. The coupled induction of βCHY1 and -2, and ABA1 with ABA biosynthesis genes during later stages of seed development that coincide with the accumulation of ABA indicates that βCHY1 and -2, and ABA1 may function to drive carotenoid intermediates towards β-xanthophyll and ultimately ABA biosynthesis during these stages. This is consistent with the observed reduced expression of LCY-ε in dry seeds and reports that βchy1 βchy2 double mutants have a reduced ability to synthesize ABA during drought stress [56]. Interestingly, while the expression of βCHY1 and -2 and ABA1 is reduced during imbibition and in dark grown plants, their expression is rapidly induced by light with βCHY1 being expressed at levels greater than two fold above that of βCHY2. The coupling of expression of these genes with other CrtBS genes in response to light may be indicative of their essential role in synthesizing β-xanthophylls which in turn are required for photoprotection.

Phytohormone regulation of isoprenoid gene expression in early development

Abscisic acid, GA and BRs have been shown to have important roles in regulating germination and post-germinative development and gene expression. Abscisic acid is known to inhibit germination and the expression of photosynthesis-related genes in imbibed seeds [54] while GA acts as essential hormone in promoting germination and etiolated development while negatively regulating ABA levels [57,58] in a process that is dependent on BRs [59]. Given the cross-talk between these important developmental-related phytohormones, we examined their role in regulating transcription of the genes using mutant and chemical treatment studies.

The induction of PSY and other photosynthesis-related genes following imbibition in continuous light is negatively regulated by ABA since in the ABA deficient mutant (aba2), induction is enhanced, while in the cyp707a1, -a2 and -a3 triple mutant, which has elevated ABA levels, the induction is almost completely abolished (Figure 3). ABA-mediated suppression is consistent with the known inhibitory role of ABA in germination and reports that ABA inhibits the expression of photosynthesis-related genes at this stage [54]. The observed ABA-mediated alteration of gene expression only occurs post-germination since gene expression levels in dry mutant seeds are not different from wild type (data not shown).

The presence of exogenous ABA or the GA biosynthesis inhibitor PAC in the growth media of light germinating seeds had very similar effects in that they strongly suppress the induction of PSY and other photosynthesis-related genes while maintaining expression of βCHY2, ABA1 and other ABA biosynthesis genes (Figure 3). The high expression level of EM6 and RD29B in PAC treated seeds indicates that these plants maintain high levels of ABA in the absence of GA. This is in line with reports that GA acts as essential hormone in promoting germination and negatively regulating ABA levels [57,58]. These results demonstrate that GA is required to activate the expression of PSY and the photosynthesis-related genes during germination in a process that is likely to involve a reduction in endogenous ABA levels. Indeed, ABA and GA are known to antagonistically regulate each others levels. GA levels in developing seeds follow an opposite trend to ABA in that they decrease progressively during seed maturation and increase sharply during germination [55].

The PAC-mediated repression of gene expression in light-germinated seeds observed here is in disagreement with a report that PAC increases expression of PSY and CrtBS genes in dark grown seedlings [14]. The discrepancy may be explained by differences in the growth conditions since in the above study [14], seeds were germinated in light for two to six hours before being grown in the dark for three days in the presence of PAC, whereas in the study analyzed here, seeds were stratified and germinated in the presence of light and PAC or ABA [60]. Thus, it appears that GA is essential for the early induction of these genes immediately following germination but may inhibit their expression at later developmental stages. This confirms reports that GA is required for germination and involved in the establishment of etiolated seedling development in darkness while repressing photomorphogenesis in a process that is dependent on BR [59].

The BRs also appear to have a positive role in regulating the expression of the genes. The expression of the PSY-CIIG was strongly reduced to non-detectable levels in both root (six day old) and whole shoot tissue (four day old) in the BREVIS RADIX (brx) loss-of-function mutant; brx has an impaired root development phenotype due to a root-specific BR deficiency [61] (Figure 4). The transcription of the genes in brx was rapidly restored to control levels following three hour brassinolide (BL) treatment strongly indicating that the miss-regulated expression was due to BL deficiency [61]. The addition of BL to wild type plants, however, failed to alter expression of the genes, indicating that while optimal levels of BL are required for correct expression of the genes, excess BL does not induce further expression. In contrast to the reduced expression of photosynthesis-related genes in the 'minimal' shoot tissue of young brx seedlings [61], the shoot system morphology, including leaves of mature brx plants have been reported to resemble that of wild-type plants [62]. This suggests that the repression of photosynthesis-related genes in the shoot system of brx may only be temporary in young developing seedlings. The effect may result from ABA-mediated inhibition since BRs have been shown to positively regulate germination by reducing ABA sensitivity [63]. In addition, along with GA, BRs have been found to function in the establishment of etiolated development in Arabidopsis seedlings while repressing photomorphogenesis [59]. Thus, like GA, BL appears to be required for the normal expression of carotenoid and chlorophyll biosynthesis genes in young post-germination tissue in a process that may involve inhibition of ABA sensitivity and biosynthesis.

thumbnailFigure 4. Heat map illustration the fold change in expression of the PSY-CIIG in developmental-related mutants. The experimental conditions included, etiolated and de-etiolated growth in pifq mutants (GSE17159), effect of constitutively active PHYBY276H allele in dark-grown phyAphyB mutant seedlings (GSE8951), exposure of dark grown phyA and phyB mutants to cFR and cR light respectively (Tepperman 2001, 2004) and effect of brassinolide (BL) in roots (6 day) and shoots (4 day) of young brx loss-of-function plants (E-MEXP-635) Details of the microarray experimental conditions are presented in Additional File 2 (Supporting Information).

GA biosynthesis genes

The expression of GPPS that functions in GA biosynthesis, increased progressively during seed development, remained high in dry and imbibed seeds before being reduced in dark- and light-grown plants (Figure 3). This is consistent with the established role of GPPS in the synthesis of GGPP pools for GA biosynthesis, but not carotenoid and chlorophyll biosynthesis [48]. However, the GA1 and GA2 genes that encode plastid localized enzymes that directly use GGPP as a substrate showed minimal differential expression throughout. The expression profile of GA3 was similar to GPPS which is in line with their functional roles in the early steps in the GA biosynthesis pathway. The high expression level of GPPS in maturing, dry and imbibed seeds is likely to contribute to the increased synthesis of GGPP pools that is required for the increase in GA that occurs during germination [55,58]. The expression profile of these GA biosynthesis genes is in marked contrast to those of GGPPS1, PSY and GGR2 that decrease during seed maturation and increase in dark- and light-grown seedlings. Thus, these distinct expression profiles are entirely consistent with GPPS functioning in the synthesis of GGPP pools for GA biosynthesis and GGPPS1 catalyzing the synthesis of GGPP precursors for carotenoid and chlorophyll biosynthesis. A number of the late GA biosynthesis genes including the GA20 and GA3 oxidases that encode enzymes catalyzing the final steps in the synthesis of bioactive GAs are also strongly induced during imbibition (data not shown) in line with elevated GA levels that occur at this stage [58].

Regulation of isoprenoid gene expression by PHYs and PIFs during early development

The expression of the gene groups was next examined in a time course experiment where dark grown seedlings were exposed to cFR and cR. This experiment was performed on the 8 K Arabidopsis microarray chip that does not include all genes used in our analysis; it does however include PSY and a number of other genes under investigation and thus provides a reasonable representation of the biosynthesis pathways being examined [64,65]. It is not unexpected that exposure to both cFR and cR induced a largely universal increase in expression of the isoprenoid biosynthesis genes including, PSY, GAPDβ, GGR2, GluTR, GSA2 and ChlSyn (Figure 4, for 6 h time point). These light responses were additionally examined in phy-A and -B mutants since the PHYs are considered to be the predominant photoreceptors that mediate light-induced germination [66]. The cFR-induction was largely abolished for most genes in the phyA mutant while the cR-induction remained largely unaltered in the phyB mutant; these studies illustrate that phyA is required for early cFR-induced gene expression while the early cR induction can be mediated by PHYs other than PHYB. The importance of PHYA as an essential signaling component of cFR-regulated gene expression is well documented with a number of studies reporting that phyA mutants are disrupted in cFR-induced expression [64,65,67]. In addition, PHYA has been shown to be the dominant PHY in mediating the induction of early-response genes to cR [64,68,69] and to exert an early functional role in inhibiting hypocotyl growth [70]. However, while PHYB does not appear necessary to activate early cR-induced gene expression, phyB mutants have been reported to display a distinct morphological phenotype in cR, including long hypocotyls and small cotyledons, pointing to an important functional role for PHYB in plant photomorphogenesis [67].

While it is evident that cFR activates expression of many carotenoid and chlorophyll biosynthesis genes, and can induce de-etiolation (repress hypocotyl elongation) via a PHYA-dependent mechanism [67], it does not activate chlorophyll biosynthesis or chloroplast development which is dependent on light-induced activation of POR that catalyzes the conversion of Pchlide to chlorophyllide [22]. Indeed, it has been demonstrated that dark-grown seedlings exposed to cFR have carotenoid and chlorophyll contents that are around 80% and 20% respectively of the levels present in seedlings exposed to cR [18]. Further, cFR grown plants have a phenotype that is an intermediate between dark and cR grown plants [71] and cFR induces a PHYA-dependent growth pattern essential for soil emerging seeds or seedling survival in conditions of deep canopy shade which are characterized by reduced ratios of R:FR [72]. Hence, while expression of all three Arabidopsis POR genes (-A, -B and -C) is high in dark grown seedlings, the activation of their enzymatic activity and the induction of chlorophyll biosynthesis is ultimately light-dependent.

The PHYs are known to activate gene expression following their light-induced translocation from the cytoplasm to the nucleus where they specifically interact with PIFs and mediate their degradation [73,74]. The PIFs are a subset of basic helix-loop-helix (bHLH) transcription factors (TFs) that bind to the promoters of light-induced genes and function somewhat redundantly to repress their expression and photomorphogenesis in dark-grown seedlings [23,73]. The PHY-mediated degradation of PIFs allows activation of light-induced genes and pif mutants have been shown to have a constitutive photomorphogenic (cop)-like phenotype in true dark-grown seedlings [75]. A recent study showed that the PIF1 TF binds specifically to G-box cis motifs present in the PSY promoter and mutant studies revealed that PIF1 and other members of the PIF family function to inhibit PSY expression and carotenoid and chlorophyll biosynthesis in dark-grown seedlings [24]. The expression of the gene sets was thus examined in response to a number of PIF loss-of-function mutants in order to provide a broader systems perspective of the role that PIFs have in regulating the synthesis of chloroplast localized isoprenoid derived compounds.

The expression of the gene sets was not substantially altered in dark grown pif single and double mutants including pif1, pif3 and pif4,5 (data not shown) and this is probably a reflection of their redundant functions. In the pif -1,-3,-4,-5 quadruple mutant (pifq), however, expression of the CrtBS genes and other genes involved in the synthesis of the photosynthetic apparatus reached quantitatively similar levels to that observed in 2 day cR-exposed wild-type plants (Figure 4). As reported previously for most dark grown pifq differentially expressed genes [23], the induction of the CrtBS genes occurs post-germination since their transcript levels are similar in pifq and wild type seeds. In this study [23], PSY and a number of other genes investigated, including DXPS2, βCHY2, ABA1, GluTR, GUN5, chlorophyllide a oxygenase (CH1) and CHL1 were identified as direct target candidates of PIF-mediated repression in the dark based on their expression being, firstly, elevated in dark grown pifq mutants compared to wild types, secondly, rapidly elevated after one hour Rc exposure (stimulates rapid PHY-induced PIF degradation ) and thirdly, sustained after germination in two days cR. In addition, 84% of genes in the PSY-ECG50 are induced >1.5 fold in dark-grown pifq mutants (Additional File 3).

While the PIFs clearly appear to negatively regulate expression of the genes in dark-grown seedlings, it is noted that expression of many genes including PSY, Calvin cycle genes, MEP pathway and chlorophyll biosynthesis genes (Figure 3) as well as most genes in the PSY-ECG50 (Additional File 3) were previously shown to be strongly induced in dark-grown wild-type seedlings when compared to stratified seeds. Although not as great as in response to light, for some genes, including, the three GAPD subunit encoding genes, PSY, LCY-ε, DXPS2, MCT, GGR2, ChlSyn and many other chlorophyll biosynthesis genes, the increase was greater than two-fold illustrating that expression of these genes is positively regulated in dark-grown wild-type seedlings. This documents that while the PIFs limit gene expression in dark grown seedlings, the inhibition is not absolute and that increases in expression do occur in the dark in the presence of PIFs. This is consistent with studies which show that increases in the biosynthesis of carotenoids and chlorophyll precursors in the dark is required for optimal greening upon light exposure [14].

In another related mutant experiment, dark-grown phyA phyB double mutant seedlings expressing the constitutively active Y276H missense allele of Arabidopsis PHYB (PHYBY276H) [76], were similarly able to mimic the cR-induced transcriptional activation of the gene sets as observed in dark-grown pifq mutants. In this mutant, PHYBY276H undergoes light-independent nuclear localization which may mediate degradation of PIFs and allow the induction of light-inducible genes [77,78]. Thus, while the phyB mutant experiment indicates that cR induction of the gene sets can occur independently of PHYB, this experiment clearly shows that active and nuclear localized PHYB can induce expression of the genes in the dark.

In summary, these results strongly support that PHYA and PHYB have important functional roles in coordinating the transcription of the interrelated isoprenoid carotenoid and chlorophyll biosynthesis genes during de-etiolation. This process most likely involves the PHY-mediated degradation of PIFs, thus enabling light-induced gene expression.

Carotenoid gene expression and ABA biosynthesis in response to osmotic stress

The carotenoids are precursors for ABA biosynthesis and we reported in this study that expression of some late CrtBS genes is induced at a time that coincides with increased ABA biosynthesis in maturing seeds (Figure 3). We therefore next examined expression of the PSY-CIIG in shoot and root tissue in a time course response to osmotic stress (mannitol) which induces the synthesis of ABA and can thus help resolve how expression of CrtBS genes is coordinated with that of ABA biosynthesis in these tissues. The experimental results reveal some interesting tissue specific expression response patterns (Figure 5). Not surprisingly, responses were more immediate in root tissue where the stress was applied, resulting in an early and sustained increase in expression of a number of the genes including, Calvin cycle genes, MEP pathway genes and dedicated CrtBS genes including, PSY, ZDS, βCHY1 and -2, ABA1 and VDE. In a similar manner to maturing seeds, this increase was paralleled with a strong increase in expression of a number of ABA biosynthesis genes including NCED3, and the ABA-responsive genes, EM6 and RD29B, suggesting an increase in endogenous ABA levels [79]. It is noted that there was little change in the expression of chlorophyll biosynthesis genes here.

thumbnailFigure 5. Expression heat map time-course illustrating the effect of osmotic stress on expression of the PSY-CIIG in root and shoot tissue. Fold-change (log2) in gene expression was measured in root and shoot tissue at the indicated time points following continuous osmotic stress application (mannitol) to root tissue (ME00327). Details of the microarray experimental conditions are presented in Additional File 2 (Supporting Information).

The increase in expression of CrtBS pathway genes in root tissue is in contrast to that in the shoot where there is a general reduction in expression of carotenoid and chlorophyll biosynthesis genes from around 3-6 h which progressively decreases up to 24 h. However, there is a strong and transient induction of βCHY2 and ABA1 between 3-12 h in shoot tissue while NCED3 expression is induced early and sustained for the duration. As observed previously, βCHY2 and ABA1 expression is also strongly induced independent of other CrtBS genes during seed maturation, a process that also requires increased ABA biosynthesis and illustrates that expression of these genes can be uncoupled from other CrtBS genes, in both non-photosynthetic (seeds) and photosynthetic tissues, under conditions that require increased ABA biosynthesis.

The NCEDs have been proposed to be key regulators of ABA synthesis since their increased expression is correlated with increased endogenous ABA concentrations [80]. Notably here, NCED2 and -3 were the predominant NCEDs induced in root and shoot tissue, which is in contrast to developing seeds where induction of the NCED-5, -6, and -9 genes parallels the increase in ABA production.

The more universal induction of carotenoid-related biosynthesis genes in osmotically stressed roots may reflect the lower concentration of carotenoid precursors that are present in this tissue. Since photosynthetic tissue contains high concentrations of epoxycarotenoids, it appears that transcription, and presumably translation of only the late CrtBS genes, βCHY2 and ABA1 are required to increase violaxanthin precursor levels for ABA biosynthesis. In contrast, in root tissue which is a major site of ABA biosynthesis, low concentrations of carotenoids may be rate-limiting for ABA biosynthesis. We therefore propose that an increase in Calvin cycle and MEP pathway genes and a more universal induction of the CrtBS genes is required in root tissue to generate violaxanthin precursors for ABA biosynthesis, a hypothesis that is supported by studies in maize [81] and rice [82]. The absence of any change in the expression of chlorophyll biosynthesis genes documents that the expression of the carotenoid and chlorophyll genes can be regulated independently, at least in root tissue.

Promoter enrichment analysis

The high expression correlation values of the genes within the PSY-ECG50 points to the possibility that their expression is coordinately regulated. A promoter content analysis was therefore performed in an attempt to identify the presence of enriched putative regulatory elements that may be causative for their co-expression. The analysis was performed examining regions 2000 base pairs (bp) upstream of the coding region/translation start sites (TlSS) of genes since regulatory elements have previously been identified in 5 prime untranslated regions (5' UTR) of light-induced genes [83,84]. This is particularly relevant to PSY which is annotated in TAIR to have a 779 bp sequence upstream of the coding region that includes two 5'UTRs and an intron. A number of elements were found to be significantly enriched in the promoters of the co-expressed PSY-ECG50 and are thus considered candidate regulatory elements that may coordinate their transcription. In addition, a number of these elements correspond to known plant cis regulatory elements including, a slightly degenerate G-box (CACGNG (p-value = 9.8-03) compared to CACGTG) and the auxin-responsive element (AuxRE, TGTCTC (p-value = 0.02), Additional File 4).

Additional file 4. Additional Table 2. Enriched motifs identified in the promoters of genes in the PSY-ECG50.

Format: DOC Size: 70KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The G-box is known to be present in the promoters of many light-regulated genes [85-87] and is enriched in the promoters of genes that are rapidly-induced by PHYA [88]. As previously mentioned, G-boxes present in the promoter of PSY have been shown to specifically bind the PIF1 TF resulting in inhibition of PSY expression [24]. Thus, the identification of an enrichment of G-boxes in the promoters of the PSY-ECG50 is consistent with the observed PHYA dependency for induction of these genes and the inhibitory effect of PIFs on their expression in dark grown plants. We noted that one of the two G-boxes identified in the PSY promoter in this analysis is positioned in the 5' UTR in close proximity to the TlSS (-21 to -16 and -919 to -914) and both differ from those identified previously [67] which examined promoters more than 2000 bp upstream of the TlSS.

The enrichment of AuxREs in the promoters of the genes is consistent with the observed BR-dependency for CrtBS gene expression in young tissues. While the AuxRE was initially believed to confer auxin responsiveness to promoters [89,90], more recent studies have indicated that this element is also a target of BR signaling. It has been suggested that the AuxRE should in fact be considered a BR-AuxRE [91] since it has been found to be enriched in auxin- and BR-responsive genes rather than genes specifically regulated by auxin [92]. Thus, the enrichment of BR-AuxRE in the promoters of the genes is in line with the BR-dependent gene expression observed in young tissues in this study and additionally adds strength to studies that have shown BRs have a role in the establishing the etiolated development program in dark-grown Arabidopsis seedlings [59].

In summary, the GA- and BR-dependent induction of the carotenoid and chlorophyll biosynthesis genes following germination is entirely consistent with the role of these hormones in establishing etioplast development [59] and the requirement for carotenoid and chlorophyll precursor accumulation in developing etioplasts [14]. The expression level of the genes appears to be restricted by the PIF TFs in dark grown plants which are subsequently degraded by light activated PHY molecules allowing a strong and coordinated induction of the genes and a subsequent increase in carotenoid and chloroplasts biosynthesis. The identification of an enrichment in putative BR-AuxRE and G-boxes in promoters of the PSY-ECG50 complements the observed transcriptional regulatory roles of BRs and PIFs respectively.

Conclusions

The tightly coupled expression and induction of PSY and many other isoprenoid biosynthesis genes throughout key developmental stages that correspond to the timing of increased carotenoid and chlorophyll synthesis and development of the photosynthetic apparatus strongly suggests that the coordinated transcription of these biosynthesis genes is critical in regulating and coordinating the biosynthesis of the functionally related carotenoid, chlorophyll, PQ and PhQ molecules. The phytohormones GA, BR and ABA as well as the transcriptional-related PHYs and PIFs appear to have important roles in regulating and coordinating the transcription of these isoprenoid-derived compounds.

Methods

PSY expression correlation analysis

An expression correlation analysis was performed for PSY using the freely available Arabidopsis co-expression tool (ACT) (http://www.arabidopsis.leeds.ac.uk/ webcite)[31]. This particular tool uses hybridization signal intensities from microarray experiments to calculate a Pearson correlation coefficient (r-value), which is a scale-invariant measure of expression similarity. The analysis was performed across all of the 322 available Ath1 22 K microarrays from the NASC/GARNet dataset which contain probe sets that recognize 21,891 Arabidopsis genes. The arrays included in this analysis are derived from a broad range of experimental samples including specific tissue types, developmental stages, abiotic and biotic treatments, and a range of mutants. Importantly, the ACT tool uses NASC/GARNet data sets that were labeled, hybridized and analyzed using a standardized procedure thus providing a homogeneous and readily comparable data set.

The analysis was performed leaving the gene list limit blank resulting in the return of a global correlation analysis of all probe sets relative to PSY ranging from the most positive to the most negatively expression correlated genes (total over 22,500 probe IDs). The top 50 genes that had the highest expression correlation with PSY were extracted from the list as were genes that were included in the PSY-CIIG. Both lists were filtered to include only genes that were represented by a unique probe on the microarray chip.

The co-correlation analysis was performed using the 2D scatter plot tool present on the ACT website. Probe IDs for PSY and GGR2 was inserted into the X and Y axis and all other specifically highlighted genes were inserted in the highlight option.

Functional enrichment analysis

A gene ontology (GO) analysis was performed using the "Fatigo plus" (version 3.1) compare tool in the Babelomics suite (http://babelomics.bioinfo.cipf.es/functional.html webcite) [33,93] to determine if there was any statistically enriched terms associated with the PSY-ECG50 expression correlated genes compared to the expected frequency in the complete genome. The top 50 genes were selected for this analysis since their expression was highly correlated with PSY (r-value range 0.91-0.84). All the available functional annotation options for Arabidopsis were selected which include the three GO categories of biological process (BP), cellular component (CC) and molecular function (MF) as well as KEGG pathways. The tool uses a Fisher's exact test and returns adjusted p-values (Family Wise Error Rate) to accounting for multiple testing to determine statistical significance.

Microarray stimuli specific transcription analysis

An in silico global expression analysis was subsequently performed for both gene sets in response to specific stimuli and in selected mutants to identify conditions that induce differential expression of the genes (Figures 3 and 5 and Additional File 3). The expression profiles of PSY and its positively correlated gene sets were initially screened over all of the available ATH1: 22 K array Affymetrix public microarray data in the gene response viewer tool (GRV) in Genevestigator [94]. Normalized microarray data were downloaded for experiments that were found to induce differential expression of the genes from the following sites:

NASCArrays (http://affymetrix.arabidopsis.info/narrays/experimentbrowse.pl webcite)[95], TAIR-ATGenExpress (http://www.ebi.ac.uk/microarray-as/ae/ webcite), GEO (NCBI) (http://www.ncbi.nlm.nih.gov/geo/ webcite) [96].

(see attached file for experiment descriptions).

Promoter enrichment analysis

A number of tools in the POXO (http://ekhidna.biocenter.helsinki.fi/poxo webcite) [97] promoter analysis suite were used to analyze promoter regions 2000 bp upstream of the coding regions of the genes in the PSY-ECG50. The POCO tool was used to identify enriched elements and the POBO tool was used to verify the presence of identified elements in the PSY promoter. The identified significantly enriched motifs were filtered to ensure that they were present and enriched in the PSY promoter and were present in greater than 70% of the genes in the PSY-ECG50.

Authors' contributions

ETW and CG conceived the initial project. OT and RV performed an initial analysis on the carotenoid biosynthesis pathway genes and promoters. SM expanded the project to include additional isoprenoid biosynthesis pathway genes, generated the results presented, interpreted data and wrote the manuscript with contributions by ETW, OT and CG. All authors read and approved the final manuscript.

Acknowledgements

This research was supported by grants (to ETW) from the US National Institutes of Health (GM081160) and New York State, and grants (to CG) from the South African National Research Fund and the Oppenheimer Memorial Trust (South Africa).

References

  1. Demmig-Adams B, Gilmore AM, Adams WW III: Carotenoids 3: In vivo function of carotenoids in higher plants.

    FASEB J 1996, 10:403-412. PubMed Abstract | Publisher Full Text OpenURL

  2. Nambara E, Marion-Poll A: Abscisic acid biosynthesis and catabolism.

    Annu Rev Plant Biol 2005, 56:165-185. PubMed Abstract | Publisher Full Text OpenURL

  3. Norris SR, Barrette TR, Dellapenna D: Genetic dissection of carotenoid synthesis in arabidopsis defines plastoquinone as an essential component of phytoene desaturation.

    Plant Cell 1995, 7:2139-2149. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Cordoba E, Salmi M, Leon P: Unravelling the regulatory mechanisms that modulate the MEP pathway in higher plants.

    J Exp Bot 2009, 60:2933-2943. PubMed Abstract | Publisher Full Text OpenURL

  5. Rodriguez-Concepcion M: Supply of precursors for carotenoid biosynthesis in plants.

    Arch Biochem Biophys 2010, 504:118-122. PubMed Abstract | Publisher Full Text OpenURL

  6. Beyer P, Weiss G, Kleinig H: Solubilization and reconstitution of the membrane-bound carotenogenic enzymes from daffodil chromoplasts.

    Eur J Biochem 1985, 153:341-346. PubMed Abstract | Publisher Full Text OpenURL

  7. Dogbo O, Laferriere A, d'Harlingue A, Camara B: Carotenoid biosynthesis: Isolation and characterization of a bifunctional enzyme catalyzing the synthesis of phytoene.

    Proc Natl Acad Sci USA 1988, 85:7054-7058. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Matthews PD, Wurtzel ET: Biotechnology of food colorant production. In Food Colorants: Chemical and Functional Properties. CRC Press Boca Raton; 2007:347-398. OpenURL

  9. Sandmann G, Romer S, Fraser PD: Understanding carotenoid metabolism as a necessity for genetic engineering of crop plants.

    Metab Eng 2006, 8:291-302. PubMed Abstract | Publisher Full Text OpenURL

  10. Chen Y, Li F, Wurtzel ET: Isolation and characterization of the Z-ISO gene encoding a missing component of carotenoid biosynthesis in plants.

    Plant Physiol 2010, 153:66-79. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Dall'Osto L, Fiore A, Cazzaniga S, Giuliano G, Bassi R: Different roles of alpha- and beta-branch xanthophylls in photosystem assembly and photoprotection.

    J Biol Chem 2007, 282:35056-35068. PubMed Abstract | Publisher Full Text OpenURL

  12. Fiore A, Dall'Osto L, Fraser PD, Bassi R, Giuliano G: Elucidation of the beta-carotene hydroxylation pathway in Arabidopsis thaliana.

    FEBS Lett 2006, 580:4718-4722. PubMed Abstract | Publisher Full Text OpenURL

  13. Wasilewska A, Vlad F, Sirichandra C, Redko Y, Jammes F, Valon C, Frei dit FN, Leung J: An update on abscisic acid signaling in plants and more.

    Mol Plant 2008, 1:198-217. PubMed Abstract | Publisher Full Text OpenURL

  14. Rodriguez-Villalon A, Gas E, Rodriguez-Concepcion M: Phytoene synthase activity controls the biosynthesis of carotenoids and the supply of their metabolic precursors in dark-grown Arabidopsis seedlings.

    Plant J 2009. OpenURL

  15. Welsch R, Beyer P, Hugueney P, Kleinig H, Von Lintig J: Regulation and activation of phytoene synthase, a key enzyme in carotenoid biosynthesis, during photomorphogenesis.

    Planta 2000, 211:846-854. PubMed Abstract | Publisher Full Text OpenURL

  16. Paddock TN, Mason ME, Lima DF, Armstrong GA: Arabidopsis protochlorophyllide oxidoreductase A (PORA) restores bulk chlorophyll synthesis and normal development to a porB porC double mutant.

    Plant Mol Biol 2010, 72:445-457. PubMed Abstract | Publisher Full Text OpenURL

  17. Park H, Kreunen SS, Cuttriss AJ, Dellapenna D, Pogson BJ: Identification of the carotenoid isomerase provides insight into carotenoid biosynthesis, prolamellar body formation, and photomorphogenesis.

    Plant Cell 2002, 14:321-332. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Von Lintig J, Welsch R, Bonk M, Giuliano G, Batschauer A, Kleinig H: Light-dependent regulation of carotenoid biosynthesis occurs at the level of phytoene synthase expression and is mediated by phytochrome in Sinapis alba and Arabidopsis thaliana seedlings.

    Plant J 1997, 12:625-634. PubMed Abstract | Publisher Full Text OpenURL

  19. Botella-Pavia P, Besumbes O, Phillips MA, Carretero-Paulet L, Boronat A, Rodriguez-Concepcion M: Regulation of carotenoid biosynthesis in plants: evidence for a key role of hydroxymethylbutenyl diphosphate reductase in controlling the supply of plastidial isoprenoid precursors.

    Plant J 2004, 40:188-199. PubMed Abstract | Publisher Full Text OpenURL

  20. Li F, Vallabhaneni R, Yu J, Rocheford T, Wurtzel ET: The maize phytoene synthase gene family: overlapping roles for carotenogenesis in endosperm, photomorphogenesis, and thermal stress tolerance.

    Plant Physiol 2008, 147:1334-1346. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Li F, Tzfadia O, Wurtzel ET: The phytoene synthase gene family in the Grasses: subfunctionalization provides tissue-specific control of carotenogenesis.

    Plant Signal Behav 2009, 4:208-211. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Philippar K, Geis T, Ilkavets I, Oster U, Schwenkert S, Meurer J, Soll J: Chloroplast biogenesis: the use of mutants to study the etioplast-chloroplast transition.

    Proc Natl Acad Sci USA 2007, 104:678-683. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Leivar P, Tepperman JM, Monte E, Calderon RH, Liu TL, Quail PH: Definition of early transcriptional circuitry involved in light-induced reversal of PIF-imposed repression of photomorphogenesis in young Arabidopsis seedlings.

    Plant Cell 2009, 21:3535-3553. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Toledo-Ortiz G, Huq E, Rodriguez-Concepcion M: Direct regulation of phytoene synthase gene expression and carotenoid biosynthesis by phytochrome-interacting factors.

    Proc Natl Acad Sci USA 2010, 107:11626-11631. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. 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:24. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Meier S, Ruzvidzo O, Morse M, Donaldson L, Kwezi L, Gehring C: The Arabidopsis wall associated kinase-like 10 gene encodes a functional guanylyl cyclase and is co-expressed with pathogen defense related genes.

    PLoS ONE 2010, 5:e8904. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Allocco DJ, Kohane IS, Butte AJ: Quantifying the relationship between co-expression, co-regulation and gene function.

    BMC Bioinformatics 2004, 5:18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Stuart JM, Segal E, Koller D, Kim SK: A gene-coexpression network for global discovery of conserved genetic modules.

    Science 2003, 302:249-255. PubMed Abstract | Publisher Full Text OpenURL

  29. Wei H, Persson S, Mehta T, Srinivasasainagendra V, Chen L, Page GP, Somerville C, Loraine A: Transcriptional coordination of the metabolic network in Arabidopsis.

    Plant Physiol 2006, 142:762-774. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Ihmels J, Bergmann S, Barkai N: Defining transcription modules using large-scale gene expression data.

    Bioinformatics 2004, 20:1993-2003. PubMed Abstract | Publisher Full Text OpenURL

  31. Manfield IW, Jen CH, Pinney JW, Michalopoulos I, Bradford JR, Gilmartin PM, Westhead DR: Arabidopsis Co-expression Tool (ACT): web server tools for microarray-based gene expression analysis.

    Nucleic Acids Res 2006, 34:W504-W509. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Gallagher CE, Matthews PD, Li F, Wurtzel ET: Gene duplication in the carotenoid biosynthetic pathway preceded evolution of the grasses.

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

  33. Al Shahrour F, Minguez P, Tarraga J, Medina I, Alloza E, Montaner D, Dopazo J: FatiGO +: a functional profiling tool for genomic data. Integration of functional annotation, regulatory motifs and interaction data with microarray experiments.

    Nucleic Acids Res 2007, 35:W91-W96. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Tian L, Dellapenna D, Dixon RA: The pds2 mutation is a lesion in the Arabidopsis homogentisate solanesyltransferase gene involved in plastoquinone biosynthesis.

    Planta 2007, 226:1067-1073. PubMed Abstract | Publisher Full Text OpenURL

  35. Sadre R, Gruber J, Frentzen M: Characterization of homogentisate prenyltransferases involved in plastoquinone-9 and tocochromanol biosynthesis.

    FEBS Lett 2006, 580:5357-5362. PubMed Abstract | Publisher Full Text OpenURL

  36. Motohashi R, Ito T, Kobayashi M, Taji T, Nagata N, Asami T, Yoshida S, Yamaguchi-Shinozaki K, Shinozaki K: Functional analysis of the 37 kDa inner envelope membrane polypeptide in chloroplast biogenesis using a Ds-tagged Arabidopsis pale-green mutant.

    Plant J 2003, 34:719-731. PubMed Abstract | Publisher Full Text OpenURL

  37. Munoz-Bertomeu J, Cascales-Minana B, Mulet JM, Baroja-Fernandez E, Pozueta-Romero J, Kuhn JM, Segura J, Ros R: Plastidial glyceraldehyde-3-phosphate dehydrogenase deficiency leads to altered root development and affects the sugar and amino acid balance in Arabidopsis.

    Plant Physiol 2009, 151:541-558. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Tanaka R, Tanaka A: Tetrapyrrole biosynthesis in higher plants.

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

  39. Masuda T, Fujita Y: Regulation and evolution of chlorophyll metabolism.

    Photochem Photobiol Sci 2008, 7:1131-1149. PubMed Abstract | Publisher Full Text OpenURL

  40. Floss DS, Hause B, Lange PR, Kuster H, Strack D, Walter MH: Knock-down of the MEP pathway isogene 1-deoxy-D-xylulose 5-phosphate synthase 2 inhibits formation of arbuscular mycorrhiza-induced apocarotenoids, and abolishes normal expression of mycorrhiza-specific plant marker genes.

    Plant J 2008, 56:86-100. PubMed Abstract | Publisher Full Text OpenURL

  41. Lange BM, Ghassemian M: Genome organization in Arabidopsis thaliana: a survey for genes involved in isoprenoid and chlorophyll metabolism.

    Plant Mol Biol 2003, 51:925-948. PubMed Abstract | Publisher Full Text OpenURL

  42. Okada K, Kasahara H, Yamaguchi S, Kawaide H, Kamiya Y, Nojiri H, Yamane H: Genetic evidence for the role of isopentenyl diphosphate isomerases in the mevalonate pathway and plant development in Arabidopsis.

    Plant Cell Physiol 2008, 49:604-616. PubMed Abstract | Publisher Full Text OpenURL

  43. Okada K, Saito T, Nakagawa T, Kawamukai M, Kamiya Y: Five geranylgeranyl diphosphate synthases expressed in different organs are localized into three subcellular compartments in Arabidopsis.

    Plant Physiol 2000, 122:1045-1056. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Keller Y, Bouvier F, d'Harlingue A, Camara B: Metabolic compartmentation of plastid prenyllipid biosynthesis--evidence for the involvement of a multifunctional geranylgeranyl reductase.

    Eur J Biochem 1998, 251:413-417. PubMed Abstract | Publisher Full Text OpenURL

  45. Hirooka K, Bamba T, Fukusaki E, Kobayashi A: Cloning and kinetic characterization of Arabidopsis thaliana solanesyl diphosphate synthase.

    Biochem J 2003, 370:679-686. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Hirooka K, Izumi Y, An CI, Nakazawa Y, Fukusaki E, Kobayashi A: Functional analysis of two solanesyl diphosphate synthases from Arabidopsis thaliana.

    Biosci Biotechnol Biochem 2005, 69:592-601. PubMed Abstract | Publisher Full Text OpenURL

  47. Ruppel N, Hangarter R: A mutant in geranylgeranyl diphosphate synthase 1 (GGPS1) of Arabidopsis thaliana that affects chloroplast development in adult leaves. In 19th International conference on Arabidopsis Research: 23 - 27July 2008. Montreal, Canada; 90-91. OpenURL

  48. van Schie CC, Ament K, Schmidt A, Lange T, Haring MA, Schuurink RC: Geranyl diphosphate synthase is required for biosynthesis of gibberellins.

    Plant J 2007, 52:752-762. PubMed Abstract | Publisher Full Text OpenURL

  49. Auldridge ME, Block A, Vogel JT, Dabney-Smith C, Mila I, Bouzayen M, Magallanes-Lundback M, Dellapenna D, McCarty DR, Klee HJ: Characterization of three members of the Arabidopsis carotenoid cleavage dioxygenase family demonstrates the divergent roles of this multifunctional enzyme family.

    Plant J 2006, 45:982-993. PubMed Abstract | Publisher Full Text OpenURL

  50. Tan BC, Joseph LM, Deng WT, Liu L, Li QB, Cline K, McCarty DR: Molecular characterization of the Arabidopsis 9-cis epoxycarotenoid dioxygenase gene family.

    Plant J 2003, 35:44-56. PubMed Abstract | Publisher Full Text OpenURL

  51. Iuchi S, Kobayashi M, Taji T, Naramoto M, Seki M, Kato T, Tabata S, Kakubari Y, Yamaguchi-Shinozaki K, Shinozaki K: Regulation of drought tolerance by gene manipulation of 9-cis-epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis.

    Plant J 2001, 27:325-333. PubMed Abstract | Publisher Full Text OpenURL

  52. Ghassemian M, Lutes J, Tepperman JM, Chang HS, Zhu T, Wang X, Quail PH, Markus Lange B: Integrative analysis of transcript and metabolite profiling data sets to evaluate the regulation of biochemical pathways during photomorphogenesis.

    Arch Biochem Biophys 2006, 448:45-59. PubMed Abstract | Publisher Full Text OpenURL

  53. Mayer MP, Beyer P, Kleinig H: Quinone compounds are able to replace molecular oxygen as terminal electron acceptor in phytoene desaturation in chromoplasts of Narcissus pseudonarcissus L.

    Eur J Biochem 1990, 191:359-363. PubMed Abstract | Publisher Full Text OpenURL

  54. Okamoto M, Tatematsu K, Matsui A, Morosawa T, Ishida J, Tanaka M, Endo TA, Mochizuki Y, Toyoda T, Kamiya Y, Shinozaki K, Nambara E, Seki M: Genome-wide analysis of endogenous abscisic acid-mediated transcription in dry and imbibed seeds of Arabidopsis using tiling arrays.

    Plant J 2010, 62:39. PubMed Abstract | Publisher Full Text OpenURL

  55. Braybrook SA, Harada JJ: LECs go crazy in embryo development.

    Trends Plant Sci 2008, 13:624-630. PubMed Abstract | Publisher Full Text OpenURL

  56. Tian L, Dellapenna D, Zeevaart JAD: Effect of hydroxylated carotenoid defi ciency on ABA accumulation in Arabidopsis.

    Physiologia Plantarum 2004, 122:314-320. Publisher Full Text OpenURL

  57. Seo M, Hanada A, Kuwahara A, Endo A, Okamoto M, Yamauchi Y, North H, Marion-Poll A, Sun TP, Koshiba T, Kamiya Y, Yamaguchi S, Nambara E: Regulation of hormone metabolism in Arabidopsis seeds: phytochrome regulation of abscisic acid metabolism and abscisic acid regulation of gibberellin metabolism.

    Plant J 2006, 48:354-366. PubMed Abstract | Publisher Full Text OpenURL

  58. Seo M, Nambara E, Choi G, Yamaguchi S: Interaction of light and hormone signals in germinating seeds.

    Plant Mol Biol 2009, 69:463-472. PubMed Abstract | Publisher Full Text OpenURL

  59. Alabadi D, Gil J, Blazquez MA, Garcia-Martinez JL: Gibberellins repress photomorphogenesis in darkness.

    Plant Physiol 2004, 134:1050-1057. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. 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:1887-1899. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Mouchel CF, Osmont KS, Hardtke CS: BRX mediates feedback between brassinosteroid levels and auxin signalling in root growth.

    Nature 2006, 443:458-461. PubMed Abstract | Publisher Full Text OpenURL

  62. Mouchel CF, Briggs GC, Hardtke CS: Natural genetic variation in Arabidopsis identifies BREVIS RADIX, a novel regulator of cell proliferation and elongation in the root.

    Genes Dev 2004, 18:700-714. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Steber CM, McCourt P: A role for brassinosteroids in germination in Arabidopsis.

    Plant Physiol 2001, 125:763-769. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Tepperman JM, Zhu T, Chang HS, Wang X, Quail PH: Multiple transcription-factor genes are early targets of phytochrome A signaling.

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

  65. Tepperman JM, Hudson ME, Khanna R, Zhu T, Chang SH, Wang X, Quail PH: Expression profiling of phyB mutant demonstrates substantial contribution of other phytochromes to red-light-regulated gene expression during seedling de-etiolation.

    Plant J 2004, 38:725-739. PubMed Abstract | Publisher Full Text OpenURL

  66. Mathews S: Phytochrome-mediated development in land plants: red light sensing evolves to meet the challenges of changing light environments.

    Mol Ecol 2006, 15:3483-3503. PubMed Abstract | Publisher Full Text OpenURL

  67. Ma L, Li J, Qu L, Hager J, Chen Z, Zhao H, Deng XW: Light control of Arabidopsis development entails coordinated regulation of genome expression and cellular pathways.

    Plant Cell 2001, 13:2589-2607. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Tepperman JM, Hwang YS, Quail PH: phyA dominates in transduction of red-light signals to rapidly responding genes at the initiation of Arabidopsis seedling de-etiolation.

    Plant J 2006, 48:728-742. PubMed Abstract | Publisher Full Text OpenURL

  69. Quail PH: Phytochrome-regulated Gene Expression.

    Journal of Integrative Plant Biology 2007, 49:11-20. Publisher Full Text OpenURL

  70. Parks BM, Spalding EP: Sequential and coordinated action of phytochromes A and B during Arabidopsis stem growth revealed by kinetic analysis.

    Proc Natl Acad Sci USA 1999, 96:14142-14146. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  71. Jiao Y, Ma L, Strickland E, Deng XW: Conservation and divergence of light-regulated genome expression patterns during seedling development in rice and Arabidopsis.

    Plant Cell 2005, 17:3239-3256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Yanovsky MJ, Casal JJ, Whitelam GC: Phytochrome A, phytochrome B and HY4 are involved in hypocotyl growth responses to natural radiation in Arabidopsis: weak de-etiolation of the phyA mutant under dense canopies.

    Plant, Celt and Environment 1995, 18:794. OpenURL

  73. Castillon A, Shen H, Huq E: Phytochrome Interacting Factors: central players in phytochrome-mediated light signaling networks.

    Trends Plant Sci 2007, 12:514-521. PubMed Abstract | Publisher Full Text OpenURL

  74. Nagatani A: Light-regulated nuclear localization of phytochromes.

    Curr Opin Plant Biol 2004, 7:708-711. PubMed Abstract | Publisher Full Text OpenURL

  75. Leivar P, Monte E, Oka Y, Liu T, Carle C, Castillon A, Huq E, Quail PH: Multiple phytochrome-interacting bHLH transcription factors repress premature seedling photomorphogenesis in darkness.

    Curr Biol 2008, 18:1815-1823. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Hu W, Su YS, Lagarias JC: A Light-Independent Allele of Phytochrome B Faithfully Recapitulates Photomorphogenic Transcriptional Networks.

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

  77. Oh E, Kim J, Park E, Kim JI, Kang C, Choi G: PIL5, a phytochrome-interacting basic helix-loop-helix protein, is a key negative regulator of seed germination in Arabidopsis thaliana.

    Plant Cell 2004, 16:3045-3058. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Al Sady B, Ni W, Kircher S, Schafer E, Quail PH: Photoactivated phytochrome induces rapid PIF3 phosphorylation prior to proteasome-mediated degradation.

    Mol Cell 2006, 23:439-446. PubMed Abstract | Publisher Full Text OpenURL

  79. Chaves MM, Flexas J, Pinheiro C: Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell.

    Ann Bot 2009, 103:551-560. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Schwartz SH, Qin X, Zeevaart JA: Elucidation of the indirect pathway of abscisic acid biosynthesis by mutants, genes, and enzymes.

    Plant Physiol 2003, 131:1591-1601. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. Li F, Vallabhaneni R, Wurtzel ET: PSY3, a new member of the phytoene synthase gene family conserved in the Poaceae and regulator of abiotic stress-induced root carotenogenesis.

    Plant Physiol 2008, 146:1333-1345. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  82. Welsch R, Wust F, Bar C, Al Babili S, Beyer P: A third phytoene synthase is devoted to abiotic stress-induced abscisic acid formation in rice and defines functional diversification of phytoene synthase genes.

    Plant Physiol 2008, 147:367-380. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  83. Ovadia A, Tabibian-Keissar H, Cohen Y, Kenigsbuch D: The 5'UTR of CCA1 includes an autoregulatory cis element that segregates between light and circadian regulation of CCA1 and LHY.

    Plant Mol Biol 2010, 72:659-671. PubMed Abstract | Publisher Full Text OpenURL

  84. Tompa M: Identifying functional elements by comparative DNA sequence analysis.

    Genome Res 2001, 11:1143-1144. PubMed Abstract | Publisher Full Text OpenURL

  85. Martinez-Garcia JF, Huq E, Quail PH: Direct targeting of light signals to a promoter element-bound transcription factor.

    Science 2000, 288:859-863. PubMed Abstract | Publisher Full Text OpenURL

  86. Huq E, Quail PH: PIF4, a phytochrome-interacting bHLH factor, functions as a negative regulator of phytochrome B signaling in Arabidopsis.

    EMBO J 2002, 21:2441-2450. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  87. Moon J, Zhu L, Shen H, Huq E: PIF1 directly and indirectly regulates chlorophyll biosynthesis to optimize the greening process in Arabidopsis.

    Proc Natl Acad Sci USA 2008, 105:9433-9438. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Hudson ME, Quail PH: Identification of promoter motifs involved in the network of phytochrome A-regulated gene expression by combined analysis of genomic sequence and microarray data.

    Plant Physiol 2003, 133:1605-1616. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  89. Ulmasov T, Ohmiya A, Hagen G, Guilfoyle T: The Soybean GH2/4 Gene That Encodes a Glutathione S-Transferase Has a Promoter That Is Activated by a Wide Range of Chemical Agents.

    Plant Physiol 1995, 108:919-927. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  90. Ulmasov T, Hagen G, Guilfoyle TJ: ARF1, a transcription factor that binds to auxin response elements.

    Science 1997, 276:1865-1868. PubMed Abstract | Publisher Full Text OpenURL

  91. Nemhauser JL, Mockler TC, Chory J: Interdependency of brassinosteroid and auxin signaling in Arabidopsis.

    PLoS Biol 2004, 2:E258. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  92. Goda H, Shimada Y, Asami T, Fujioka S, Yoshida S: Microarray analysis of brassinosteroid-regulated genes in Arabidopsis.

    Plant Physiol 2002, 130:1319-1334. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  93. Al Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes.

    Bioinformatics 2004, 20:578-580. PubMed Abstract | Publisher Full Text OpenURL

  94. Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox.

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

  95. Craigon DJ, James N, Okyere J, Higgins J, Jotham J, May S: NASCArrays: a repository for microarray data generated by NASC's transcriptomics service.

    Nucleic Acids Res 2004, 32:D575-D577. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  96. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data.

    Nucleic Acids Res 2009, 37:D885-D890. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  97. Kankainen M, Pehkonen P, Rosenstom P, Toronen P, Wong G, Holm L: POXO: a web-enabled tool series to discover transcription factor binding sites.

    Nucleic Acids Res 2006, 34:W534-W540. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL