Regulatory modules controlling early shade avoidance response in maize seedlings
BMC Genomics volume 17, Article number: 269 (2016)
Optimization of shade avoidance response (SAR) is crucial for enhancing crop yield in high-density planting conditions in modern agriculture, but a comprehensive study of the regulatory network of SAR is still lacking in monocot crops.
In this study, the genome-wide early responses in maize seedlings to the simulated shade (low red/far-red ratio) and also to far-red light treatment were transcriptionally profiled. The two processes were predominantly mediated by phytochrome B and phytochrome A, respectively. Clustering of differentially transcribed genes (DTGs) along with functional enrichment analysis identified important biological processes regulated in response to both treatments. Co-expression network analysis identified two transcription factor modules as potentially pivotal regulators of SAR and de-etiolation, respectively. A comprehensive cross-species comparison of orthologous DTG pairs between maize and Arabidopsis in SAR was also conducted, with emphasis on regulatory circuits controlling accelerated flowering and elongated growth, two physiological hallmarks of SAR. Moreover, it was found that the genome-wide distribution of DTGs in SAR and de-etiolation both biased toward the maize1 subgenome, and this was associated with differential retention of various cis-elements between the two subgenomes.
The results provide the first transcriptional picture for the early dynamics of maize phytochrome signaling. Candidate genes with regulatory functions involved in maize shade avoidance response have been identified, offering a starting point for further functional genomics investigation of maize adaptation to heavily shaded field conditions.
A major challenge in modern agriculture is to produce enough food with limited arable land to meet the demand from an ever-growing world population . One effective practice is to increase the crop planting density in order to enhance the yield on a per area unit basis, which has been successfully employed in maize production . A key determinant of the optimum planting density of field-grown maize is its response to shade. At dense stands, maize perceives the shade mainly as a decrease in the red (R):far-red (FR) light ratio caused by the depletion of photosynthetically active radiation (PAR) by neighboring plants. Such a signal serves as a warning of forthcoming competition from neighbors and triggers a number of physiological alterations of maize plants for them to outgrow their competitors, including reallocation of resources into stem elongation from other organs; longer leaves; lowered mechanical strength of stems; elevated ear height; higher lodging rate; early flowering and senescence; reduced photosynthesis efficiency and grain filling; and also active suppression of defense responses against pathogens and herbivores [1, 3–9]. The above-mentioned phenomena, collectively known as the shade avoidance response (SAR), is detrimental to the grain yield of maize as plant density increases.
The molecular components and genetic network controlling SAR have been extensively studied in the dicot model species Arabidopsis . The Arabidopsis genome encodes five phytochromes (phyA to E) serving as photoreceptors for R and FR light. Among them phyB plays a major role in shade detection . Phytochromes exist in two photo-convertible isoforms: R light triggers the photoconversion of phytochromes from their inactive Pr form to the active Pfr form, while FR light photoconverts the Pfr form back to the Pr form. The Pfr form of phytochormes translocates into the nucleus to trigger genome-wide transcriptional changes and subsequent photo-responses. Shade reduces the Pfr:Pr ratio, and also the nuclear accumulation, of phyB, leading to the accumulation of the E3 ligase CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1) in the nucleus  and enhanced 26S proteasome-mediated degradation of several transcription factors, including ELONGATED HYPOCOTYL5 (HY5), HY5-HOMOLOG (HYH), LONG HYPOCOTYL IN FAR-RED1 (HFR1), and LONG AFTER FAR-RED LIGHT1 (LAF1) [13–16]. Recent progress showed that phytochromes could function as transcriptional regulators by interacting directly with numerous transcription factors on the promoters of target genes, conferring rapid responses to light signals [17, 18]. More importantly, phytochromes suppress the shade response by antagonizing a group of bHLH transcription factors termed PHYTOCHROME-INTERACTING FACTORS (PIFs). Shade induces the dephosphorylation and activation of PIF7 and also promotes the protein stability of other PIFs (PIF1/3/4/5) as the result of deactivation of phyB [19–24]. PIFs regulate nearly all physiological aspects of SAR, such as flowering and auxin-dependent elongated growth [23, 25–27]. Shade also swiftly induces the transcript abundance of a number of transcription factors in Arabidopsis, including a group of homeobox genes (HAT1-4 and AtHB4) and bHLH family members (PAR1, PAR2, HFR1, PIL1, and PIF6/PIL2), which act as either positive (such as AtHB2, PIL1) or negative (such as HFR1, PAR1, and PAR2) regulators to fine-tune the progress of SAR [28–32].
Although the genetic and regulatory networks controlling SAR has been extensively studied in Arabidopsis, these are far less well defined in monocots, including maize . Maize as well as all other monocots only has three types of phytochromes (phyA to C) . Similar to Arabidopsis, the maize SAR is largely mediated by phyB, and the two copies of phyB (phyB1 and phyB2) derived from an ancient tetraploidization event contributes differently to distinct physiological aspects of SAR, indicating both redundancy and sub-functionalization of phyB activities . Dramatic physiological and morphological variations of shade response have been observed in natural populations of maize, and modern maize varieties are generally associated with enhanced shade tolerance under high-density planting conditions [2, 35, 36], thus attenuated SAR might have been historically selected by maize breeders for adaptation to high density planting . A number of domestication genes controlling important agronomic traits have been cloned in maize, such as teosinte branched1 (tb1, controls tillering) and teosinte glume architecture1 (tga1, controls glume development) [37, 38], but none has been identified to regulate SAR hitherto. It is highly likely that the improved shade response in modern maize is a relatively complex trait collectively controlled by multiple loci, considering the complexity of the regulatory network of SAR in Arabidopsis. Despite a lack of comprehensive understanding of the genetic architecture of SAR in maize and many other crop species, there has been significant interest in modifying photoperception for crop improvement. For example, overexpression or targeted modification of phyA or phyB has been used to improve plant architecture and increase yield in densely planted crops [39–45]. Other master regulators of shade response downstream of phytochromes such as PIFs, however, have not been utilized as targets in crop SAR engineering.
To improve maize SAR more efficiently, either by molecular engineering or traditional breeding methods, it is necessary to achieve a mechanistic understanding of the molecular components and regulatory networks controlling maize SAR . Here, the maize SAR controlled by phyB signaling, and also the FRc-mediated de-etiolation controlled by phyA signaling, was examined by time-series genome-wide expression profilling. The shade avoidance transcript profile significantly overlapped with that in de-etiolation, with contrasting directions of expression regulation for most overlapping differentially regulated genes in the two experiments. Enrichment analysis of cis-elements on light-responsive gene promoters identified binding sites for TFs involved in various signaling pathways that cross-talk with phytochrome signaling. Co-expression network analysis identified two transcription factor modules each functioning in SAR and de-etiolation, and may serve as regulatory hubs in each process. A cross-species comparison of the shade transcriptome between maize and Arabidopsis identified both conserved and divergent patterns of expression regulation on various biological processes such as flowering and elongated growth. Moreover, it was found that the two subgenomes of maize differentially retained shade/de-etiolation-responsive genes and cis-elements, indicating fractionation of the two subgenomes in light-signaling. This study will be helpful for future functional genomics investigation of maize adaptation to heavily shaded field conditions.
Results and discussion
The simulated shade environment
Genome-wide transcriptional dynamics of the shade avoidance response (SAR) has been extensively studied in the model dicot Arabidopsis [46–48], but less well defined in economically important monocots. This study aimed to dissect the regulatory network controlling the maize SAR, a process believed to be attenuated by intense artificial selection during domestication and breeding to avoid the adverse effects of SAR under densely planted field conditions . A simulated shade condition identical to that used before in Arabidopsis  was adopted to facilitate direct comparison between the two species (Fig. 1a). Maize seedlings grown under the simulated shade condition displayed phenotypes characteristic of SAR, including elongated leaves and less anthocyanin accumulation, compared to those grown under high R/FR (Additional file 1: Figure S1), indicating the validity of such an experimental condition.
Dynamic transcriptional reprogramming in the maize shade avoidance response
To analyze the dynamics of the SAR at a genome-wide scale, global gene expression profiles upon 1 h, 3 h and 6 h exposure to low R/FR were inspected. A transcriptional profiling on maize de-etiolation under FRc was also conducted (Fig. 1a), a process sharing common regulatory components with SAR . Principal component analysis (PCA) showed high degree of reproducibility among biological replicates within each group (Additional file 2: Figure S2). A total of 580 genes were identified as differentially regulated in response to the simulated shade within 6 h (absolute fold change > =2 and FDR adjusted p-value < =0.05), with 333 up-regulated genes and 247 down-regulated genes. Comparable numbers of DTGs were found in FR-mediated de-etiolation, with 336 genes up-regulated and 194 down-regulated (Fig. 1b). Due to the partial overlap between the two DTG sets, in total 988 early-responsive DTGs in response to shade and/or FR were observed, representing 2.5 % of the maize transcriptome (Summarized in Additional file 3: Table S1). Validation of the expression levels of 20 genes by qRT-PCR showed a high correlation (R2 = 0.938) between RNA-seq and qRT-PCR (Additional file 4: Table S2).
The number of DTGs increased with different dynamics in the two experiments. Most DTGs occurred between 3 h and 6 h in SAR, but between 1 h and 3 h in de-etiolation, suggesting faster responses of maize transcriptome in de-etiolation (Fig. 1c). There was only a limited overlap between shade- and FR-responsive gene sets (122 out of the 988 DTGs were present in both DTG sets), and most common DTGs (110 out of the 122 common DTGs) were regulated in opposite directions in SAR and de-etiolation (Fig. 1b). In agreement with this, there was a negative correlation in the expression pattern of these 110 common DTGs between SAR and de-etiolation, with Pearson correlation coefficients of −0.57, −0.53 and −0.34 for 1 h vs. 0 h, 3 h vs.0 h, and 6 h vs. 0 h, respectively. Genes with such expression patterns were also reported in Arabidopsis . Such phenomenon is expected since according to the model established in Arabidopsis, SAR and FR-mediated de-etiolation are initiated mainly by inactivation of phyB and activation of phyA, respectively, and these two classes of phytochromes share overlapping downstream signaling cascades . Therefore, DTGs regulated in opposite directions are likely primary targets downstream of phytochromes. In support of this notion, the aforementioned negative correlation weakened with time, indicating a gradual increase of secondary effects downstream of phytochrome signaling cascades.
To gain an insight into the biological processes involved in maize SAR and de-etiolation, the 988 DTGs were grouped using the K-Means clustering algorithm embedded in MapMan , and 12 gene clusters (K1-K12) were identified (Fig. 1d and Additional file 3: Table S1). All clusters showed enrichment for one or more MapMan bins (functional categories) (Fig. 2d). For example, genes involved in photosynthesis, major and minor CHO metabolism, glycolysis, transcriptional regulation and protein folding were enriched in K4, which contained shade-repressed genes. Notably, the K8 cluster, which represents genes that were highly induced by shade and therefore of particular interest, were enriched for genes unassigned with any functions. It was caused by the fact that all the 47 genes in K8 do not have functionally characterized orthologs in other organisms except two orthologs of Arabidopsis FLAVIN-BINDING, KELCH REPEAT, F-BOX1 (FKF1) (Additional file 3: Table S1). During FR-mediated de-etiolation, early-induced genes (K12, induced within 1 h FR treatment) and late-induced genes (K10) showed differential enrichment of MapMan bins. Although they were both enriched for genes involved in light signaling and hormone metabolism, K12 was particularly enriched with transcription factors while K10 was enriched for genes of phenylpropanoid and flavonoid metabolism pathways, suggesting the importance of transcription factors as early primary targets of phyA signaling in maize.
Differentially enriched cis-elements in distinctly light-regulated gene clusters
Genes with similar expression patterns are likely to share common regulatory cis-elements in their promoters, thus the enrichment of cis-elements for each DTG cluster was analyzed. This led to the identification of a total of 119 cis-elements enriched in at least one DTG cluster (Additional file 5: Table S3). Cis-elements enriched in more than 6 clusters or with functions associated with light signaling were shown in Fig. 2. K7, K8 and K1 represent shade-induced DTGs (Fig. 1d). One striking common feature of them is significant enrichment of I-box, the evening element, CCA1 binding site, and the light-responsive GATA motif. In addition, the K8 genes, which were most highly induced by shade, were enriched with the binding sites for a group of HB transcription factors, including HB1/5/6 and ARR10, consistent with the finding in Arabidopsis that HB family members were among the early shade-induced genes . Notably, several enriched motifs for K1, K7 and K8 genes, such as the light-responsive GATA motif, ARR10 binding site and GCbox, were not enriched on the promoters of shade-responsive Arabidopsis genes , suggestive of the unique functions of corresponding transcription factors in maize SAR. K12 and K10 genes, which represent early and late FRc-induced DTGs, respectively, were enriched with G-box motifs, especially binding sites for MYC2/3/4, implying a potential role of JA pathway components in FR-mediated de-etiolation. In addition to the binding sites of MYC2/3/4, cis-elements related to various other hormone pathways were found enriched in the DTG set, such as those related to auxin (ARF binding site motif), GA (ARR10 binding site motif) and BR (binding sites of BES1 and BZR1) (Additional file 5: Table S3 and Fig. 2), suggesting combinatorial complexity of these cis-elements together with light-regulated cis-elements on the transcriptional regulation of target genes. To investigate whether functionally related cis-elements might have similar enrichment patterns, cis-elements were hierarchically clustered based on their pattern of enrichment across all DTG clusters (Fig. 2). Transcription factor binding motifs with similar sequences tend to group together, such as PIL5/BZR1 and PIFs/MYCs. It was reported that over half of PIL5 target genes are also targets of BZR1 , but the integration of other transcription factors on light controlled promoters still awaits further investigation. Despite the abundance of G-box like elements in this analysis, enrichment of HY5 binding sites was not observed in any DTG cluster, consistent with the fact that HY5 transcript abundance was not significantly altered during SAR in this study. In Arabidopsis, HY5 was initially considered to play no roles in SAR under simulated shade conditions . But it was later found that HY5 was late induced by shade, and HY5 binding sites were weakly enriched in shade-responsive genes after prolonged (>1 d) shade treatment . And two recent studies identified HY5 as a repressor of SAR when a sunfleck interrupts shade light on a daily basis  and controls hypocotyl and leaf morphology in response to shade . A lack of HY5 binding site enrichment in this dataset suggests a limited role of HY5 in the early signaling events of SAR in maize, or distinct recognition sites of maize HY5 compared to Arabidopsis HY5.
Changing inventories of transcription factors in SAR and de-etiolation
Most functionally characterized domesticated genes are transcription factors [54, 55], and it was suggested that maize SAR has been targets of domestication in order for maize to adapt to the high-density stands in modern field conditions [1, 40]. Thus special attention was paid to TFs that are transcriptionally regulated, and are therefore potentially involved, in SAR. To identify TFs in our DTG sets, DTGs with homology to known TFs were combined with information from the GRASSIUS database, followed by manual removal of mis-annotated genes in GRASSIUS (Additional file 6: Table S4). Among the 580 shade-responsive DTGs, 27 up-regulated and 33 down-regulated differentially transcribed TFs (DTTFs) were identified. A larger proportion of DTGs in de-etiolation was identified as transcription factors (53 up-regulated and 32 down-regulated DTTFs among 530 DTGs) (Fig. 3a), but this set of DTTFs only shared limited overlap with that of SAR. More DTTFs were regulated in opposite directions than those regulated in the same direction in SAR and de-etiolation, a scenario also observed for DTGs (Fig. 3b). Interestingly, in contrast to the previous report that over 30 Arabidopsis TFs displayed over two fold induction or repression within 1 h of shade treatment, only 4 up-regulated TFs were found in maize at 1 h of the same treatment, including ZmHB53 (orthologous to Arabidopsis HB2/HAT4, a marker gene for SAR [56, 57]), ZmLBD6 (orthologous to Arabidopsis LBD37, which controls leaf morphogenesis and nitrogen metabolism ), ZmMYBR97 (a putative ortholog of Arabidopsis ENHANCER OF TRY AND CPC 3 (ETC3), a regulator of epidermal cell differentiation and nitrate signaling [59, 60]), and ZmGLK54 (a putative ortholog of Arabidopsis LUX). Nevertheless, maize TFs responded to FR rather swiftly: 25 up-regulated and one down-regulated transcription factors were identified from the 52 early-responsive DTGs at 1 h FR treatment (Fig. 3a). This is comparable to the scenario in Arabidopsis, in which over 30 % of the early-responsive genes in FR-mediated photomorphogenesis were transcription factors, representing a master set of genes that orchestrate the expression of the downstream targets in the phyA-directed transcriptional network .
To detect potential family-specific expression trends, DTTFs were clustered using the Self Organization Tree Algorithm (SOTA) based on their expression pattern in SAR and de-etiolation, and further categorized into 3 groups (G1-G3) (Fig. 3c). G1 and G2 represent TFs that were induced or repressed by shade, respectively. Some TF family members preferably or exclusively belong to G1 and/or G2, such as AT hook containing TFs, AP2/EREBP, MADS, GATA, PRR, MYB, SIG and DOF. TFs that were non-responsive to shade, but nevertheless FR-inducible, comprised G3. TIFY and WRKY TFs preferentially belong to G3, and are therefore more functionally related to phyA signaling during de-etiolation (Fig. 4d).
Homology-based function prediction of DTTFs in maize SAR showed many biological processes potentially regulated by these DTTFs, including light signaling, circadian clock, flowering, stress response, hormone signaling, as well as growth and development (Fig. 3e). Eleven shade-regulated BBX family members were identified, including 3 shade-induced (orthologous to Arabidopsis BBX12, BBX22, and BBX24) and 8 shade-repressed genes (orthologous to Arabidopsis BBX4, BBX5, BBX18, BBX25, and BBX32). In Arabidopsis, BBX22 negatively regulates SAR [62, 63], while BBX24 promotes SAR by impairing DELLA activity . Moreover, BBX18 and BBX25 promote hypocotyl elongation under a low R/FR ratio [65, 66]. The roles of other BBX genes in SAR have not been defined, although BBX4 and BBX32 were identified as positive and negative regulators of photomorphogenesis, respectively [67, 68].
Plastid was considered as an important player in light signaling , and two putative chloroplast localized sigma factors (SIGA/SIG2 and SIGE/SIG5) were found shade-repressed. SIGA, but not SIGE, was also reported as a shade-repressed gene in Arabidopsis . Moreover, Arabidopsis SIGA and SIGE expression displayed phytochrome-dependent induction by R or FR illumination during de-etiolation, and function as positive and negative regulators of photomorphogenesis, respectively [70–72].
TFs with functions in circadian clock and flowering were also among shade regulated DTTFs. Orthologs of ELF3, LUX, TOC1, and PRR5 were up-regulated while those of LHY, PRR7, RVE2, RVE8 were down-regulated. Two Dof-type zinc finger genes orthologous to Arabidopsis CDF2 and CDF3 were down-regulated by shade. Arabidopsis CDF2 and CDF3 repress CO expression to delay flowering and such repression is antagonized by the GI/FKF1 complex [73, 74]. Interestingly, in addition to the shade-repressed CDF2/3 expression, strong induction of the two GI and two FKF1 orthologs in maize in response to shade was observed. Therefore, the mechanism of shade-induced early flowering might be conserved between maize and Arabidopsis.
Shade also influenced many aspects of growth and development in maize. Three DTTFs potentially involved in nitrogen signaling were identified, including orthologs of Arabidopsis TCP20 , ETC3 , and LBD37 . Notably, ETC3 and LBD37 are two of the four extremely early shade-responsive DTTFs (induced within 1 h), suggesting the importance of nitrogen signaling in early maize SAR. Other shade-induced up-regulated genes include an ortholog of Arabidopsis WRINKLED1, which controls biosynthesis of triacylglycerol, an important storage compound ; an ortholog of AtNAC1, which controls SAM formation ; and an HB2 ortholog potentially involved in cell expansion and proliferation. Moreover, shade suppressed the expression of ZmMYB88, an ortholog of ATMYB59, which functions to regulate cell cycle progression ; ZmMYBR35, an ortholog of KUA1, which is essential for cell expansion and ROS homeostasis ; and four single-MYB transcription factor genes orthologous to Arabidopsis RADIALIS-LIKE SANT/MYB 2 (RSM2) and RSM3. Although RSM2 and RSM3 have not been functionally studied, RSM1 overexpression plants are hypersensitive to red light during de-etiolation, characterized by a shorter hypocotyl than wild type plants .
A number of TFs related to stress response were also identified, including the upregulation of two Oxidation-related Zinc Finger (OZF) genes orthologous to Arabidopsis OZF1 and OZF2, positive regulators of the oxidative stress response and salt stress response, respectively [81, 82]; and also ZmbZIP112, an ortholog of Arabidopsis bZIP60, which regulates unfolded protein response (UPR) in ER stress . Moreover, two JA signaling components were differentially regulated, including up-regulation of a JAZ1 ortholog and down-regulation of a MYC4 ortholog.
In Arabidopsis, a large number of shade-induced TFs have been identified and functionally characterized, including a group of homeobox genes (HAT1-4 and ATHB4) and bHLH family members (PAR1, PAR2, HFR1, PIL1, and PIF6/PIL2), which function as either positive (such as ATHB2, PIL1) or negative (such as HFR1, PAR1, and PAR2) regulators of SAR [28–31]. Phylogenetic analysis of homeobox TFs in Arabidopsis and maize identified ZmHB53 and ZmHB78 as two orthologs of the five shade-induced Arabidopsis HB genes, indicating an expansion of this small HB family clade in Arabidopsis after its divergence from maize (Additional file 7: Figure S3). ZmHB53 was an extremely early shade-induced gene up-regulated by shade within 1 h, and it was found down-regulated in de-etiolation. ZmHB78 was unresponsive to any light treatment in this study. Arabidopsis PAR1 and PAR2 are two atypical bHLH TFs, and only one maize ortholog was found (GRMZM2G364528) (Additional file 8: Figure S4). However, this gene was not up-regulated by shade, although its down-regulation in de-etiolation was observed. PIFs and PILs, together with an atypical bHLH gene HFR1, comprise a small bHLH clade that are crucial components of SAR. Phylogenetic analysis showed that seven maize bHLH proteins were orthologous to Arabidopsis PIF1/3/4/5, all characterized by an N-terminal APB motif responsible for their interaction with phyB (Additional file 9: Figure S5). Although regulation of PIF activity mainly occurs at the protein level in Arabidopsis, two light-regulated maize PIFs, ZmbHLH47 (shade-suppressed) and ZmbHLH60 (FR-induced), were identified. For HFR1, PIL1 and PIL2, however, no obvious orthologs were detected in maize, suggesting their fast evolution and expansion in Arabidopsis or a loss of these TFs in maize.
The dynamics of DTTFs during FRc-mediated de-etiolation also suggests complex cascades of transcriptional control. Interestingly, up-regulation of a number of JA pathway signaling components, including 15 JAZs, 1 MYC2 ortholog, and 1 JAM3 ortholog, was observed. The integration of phyA signaling and JA pathway has been a field of intense investigation . MYC2 functions synergistically with SPA1 to suppress FR-mediated photomorphogenesis . The role of JAZ genes in photomorphogenesis is unclear, but JAZ1 protein was reported to be destablized in FRc-mediated de-etiolation and stablized during SAR, and ectopic overexpression of JAZ1 resulted in exaggerated SAR [7, 86]. In Arabidopsis, JA was reported to positively regulate the phenylpropanoid pathway and promote biosynthesis of flavonoids (especially anthocyanin) in a phyA-dependent manner , and genes involved in phenylpropanoid and flavonoid metabolism were enriched in K2, K3 and K10, clusters containing FR-induced genes in de-etiolation (Fig. 1e). Consistently, several TFs involved in phenylpropanoid biosynthesis were induced during de-etiolation, including the maize r1/colored 1 gene orthologous to Arabidopsis GLABROUS 3 (GL3), a maize ortholog of Arabidopsis IAA26/phytochrome-associated protein1 (PAP1), three orthologs of HY5, and ZmMYB40/ZmMYB-IF35, an ortholog of Arabidopsis MYB12. In Arabidopsis, PAP1 is a phytochrome interacting protein and regulates phyA-induced chalcone synthase (CHS, EC: 184.108.40.206) expression through HY5. ZmMYB-IF35 has not been functionally characterized, but its ortholog MYB12 is a transcriptional activator of CHS, flavanone 3-hydroxylase (F3H, EC: 220.127.116.11), flavonol synthase (FLS, EC: 18.104.22.168) and chalcone flavanone isomerase (CHI, EC: 22.214.171.124) . In Arabidopsis, the promotion of MYB12 expression is dependent on HY5 , consistent with the role of HY5 as a positive regulator of anthocyanin biosynthesis . It’s still unknown whether the HY5-MYB12 axis is conserved in maize, despite the up-regulation of both HY5 and MYB12 orthologs during de-etiolation in our dataset.
Co-expression analysis shows possible modularized regulatory hubs in SAR and de-etiolation
To identify co-regulated transcription factor modules that may function as core regulators of maize light signaling, a co-expression correlation network was constructed according to Pearson correlation values between TFs, and densely connected TFs (modules) were identified using the Markov cluster algorithm (MCL) (Fig. 4a). Two modules were identified at a cut-off of 10 for node numbers, and they were found each involved in SAR and de-etiolation (Fig. 4b,c). Module I consists of 13 TFs whose expression levels gradually decreased throughout the shade treatment, including four BBX genes, three orthologs of Arabidopsis LHY, two sigma factors, as well as orthologs of Arabidopsis PIF5, MYC4, RL6, and OFP6; Module II contains 13 TFs that are early-induced (within 1 h) during FRc-mediated de-etiolation, but the inductive effect of FRc on these TFs was later alleviated. Interestingly, most of them have been implicated in stress responses, such as several JA-pathway components (five JAZ genes and one JAM1 ortholog), and also orthologs of Arabidopsis WRKY40, ZAT6, ATAF2, bHLH92, all involved in defense response to biotic and/or abiotic stresses [91–94].
Co-regulated genes may function in concert in specific biological processes, making it possible to divulge the biological functions of the above mentioned TF modules by placing them in the context of a genome-wide co-expression network. The co-expression analysis was extended to all 988 DTGs identified in this study (Fig. 4d). To identify biological processes most closely linked to TF modules, the modules were considered as single nodes in the graph and their first neighbors were subjected to GO enrichment analysis. Ninety-three first neighbors of Module I were enriched with diverse GO terms including “plastid”, “photosynthesis”, terms related to primary and secondary metabolism (such as flavonoid, anthocyanin, cellular amino acid derivatives and myo-inositol hexakisphosphate biosynthesis), as well as terms involved in responses to environment changes (“response to light stimulus” and “response to sucrose stimulus”), suggesting that Module I TFs play regulatory roles in these processes (Fig. 4e). The Module II, however, appeared rather isolated in the co-expression network, and only harbors two first neighbors: GRMZM2G031432 (encoding GA2OX1) and GRMZM2G325683 (encoding an unknown protein). The scarcity of first neighbors was probably due to the fact that Module II TFs responded to FRc much faster than most other genes, leading to relatively low correlation coefficients between these TFs and other genes. The role of GA2OX genes in de-etiolation has been reported in Pisum sativum, in which HY5 positively regulates the GA catabolism gene GA2OX2 to reduce the levels of active GAs and suppress the elongated growth in skotomorphogenesis . The maize GA2OX1 gene is another GA catabolism gene, and the maize HY5 gene in Module II possibly up-regulated the GA2OX1 gene expression. Moreover, although HY5 was largely regulated at the protein level, up-regulation of maize HY5 expression during de-etiolation was observed.
Cross-species comparison of transcriptional profiles in SAR identifies auxin biosynthesis as a potential target of maize domestication
Although monocots and dicots have split around 140 ~ 200 million years ago, they share a number of common morphological and physiological responses in their shade avoidance syndrome, such as elongated growth and accelerated flowering. Thus we compared shared patterns of expression across all orthologous maize-Arabidopsis gene pairs in response to shade. Putative orthologs of shade-responsive maize genes in Arabidopsis were obtained from the Gramene database, and their expression patterns under shade were retrieved from two previous studies [46, 48]. 174 of the 333 shade-induced maize genes have orthologs in Arabidopsis. Among these 174 maize genes, 20 have orthologous Arabidopsis genes that are also induced by shade. 182 out of 247 shade-repressed maize genes have their orthologous counterparts in Arabidopsis. However, only 11 of these maize genes have Arabidopsis orthologs that are shade-repressed (Additional file 10: Table S5). These results indicate that the regulation of shade-induced genes is more conserved than that of shade-repressed genes. Interestingly, 19 orthologous gene pairs displayed contrasting expression patterns in response to shade, suggesting that the biological processes controlled by these genes may be distinctly regulated in maize and Arabidopsis during SAR (Additional file 11: Table S6).
Special attention was paid to genes controlling elongated growth and accelerated flowering, two hallmarks of SAR. Two flowering promoting genes, GI and FT, were shade-induced in both species (Additional file 10: Table S5). In Arabidopsis, GI and FKF1 promote the CO/FT module activity by multiple mechanisms (Fig. 5) [74, 96–100]. Although FKF1 was not a shade-induced gene in Arabidopsis, strong shade induction of the two FKF1 orthologs in maize was observed (Additional file 3: Table S1, Fig. 5a). Shade also repressed the expression of two maize CDFs, which are suppressed by FKF1 and negatively regulate CO expression in Arabidopsis. These results suggested a conserved GI/FKF1-CO/FT axis promoting flowering in SAR, although the regulation of FKF1 and CDFs may have distinct mechanisms in the two species.
Elongated growth is another hallmark of SAR. In Arabidopsis, this is achieved by shade-induced rapid biosynthesis of auxin dependent on TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS 1 (TAA1)  and YUCCA genes . A recently revised model of auxin biosynthesis pathway placed YUCCA genes downstream of TAA1 . In Arabidopsis, shade strongly induces YUCCAs but also exerts a weak suppression effect on TAA1 expression . As YUCCAs catalyze the rate-limiting step, the biosynthesis of auxin is strongly enhanced in response to shade. In maize, although up-regulation of a TAA1 ortholog was observed, the expression of YUCCA genes was not obviously altered (Fig. 5b). Furthermore, in contrast to the previous report that shade-responsive genes in Arabidopsis were particularly enriched with the GO term “response to auxin stimulus” , this and other auxin-related GO terms were not enriched in shade-responsive maize genes (Additional file 12: Table S7), consistent with the distinct regulatory patterns of auxin biosynthetic genes in the two species. It is possible that the unresponsiveness of YUCCAs to shade is due to artificial selection for attenuated SAR in maize, but confirmation of such a hypothesis still awaits further investigation of YUCCA expression patterns in the undomesticated ancestors of maize.
Biased distribution of DTGs toward maize1 is associated with fractionation of cis-elements between the two maize subgenomes
An ancient tetraploidy event 5–12 million years ago gave rise to the two subgenomes of maize, with the maize1 subgenome experiencing less gene loss and displaying generally higher levels of gene expression compared to maize2 , but it’s still unknown whether the two subgenomes respond to environmental stimuli differently. To study the response of the two subgenomes in SAR and de-etiolation, the distribution of DTGs in the two subgenomes was investigated. 212 and 105 DTGs in SAR could be unambiguously assigned to maize1 and maize2, respectively, indicating more dynamic transcriptomic changes in maize1 (p < 0.001 in chi-squared test, with the null hypothesis that DTGs fall into maize1 and maize2 with an equal probability) (Fig. 6a, Additional file 13: Table S8). In each subgenome, induced or repressed DTGs were categorized according to their counterparts (homologs) in the other subgenome: 1) counterpart not differentially expressed; 2) counterpart with a similar induction/repression pattern; or 3) counterpart lost after tetraploidy (fractionation). The number of genes falling into the second group should be equal between the two subgenomes. Genes in the third group in maize1 outnumbered those in maize2, which is expected since maize2 suffered from more severe gene loss than maize1. Interestingly, distribution of the first group genes significantly biased toward maize1, for both up- and down-regulated DTGs. Similar gene distribution was also observed for de-etiolation under FRc (Fig. 6b, Additional file 13: Table S8). In sum, the higher responsiveness to shade/de-etiolation for maize1 is due to both higher gene retention rate in maize1, and also higher retention rate of light-responsiveness in existing genes.
The reduced responsiveness of maize2 to light stimuli might be indicative of fractionation of light-regulated cis-elements between the two subgenomes. To test this hypothesis, it was investigated whether cis-elements enriched in our DTG sets were differentially retained in maize1-maize2 gene pairs. A “shade-induced gene pair” was defined as a pair of homologs each belonging to a subgenome, with at least one of them shade-induced. Enrichment of cis-elements was evaluated separately for maize1 and maize2 members of all “shade-induced gene pairs”. As shown in Table 1, distribution of 3 cis-elements were found strongly biased toward maize1, including the Evening Element (EE), FUS3 binding site, and SORLREP3. The strongest bias was observed for EE, indicating a significant loss of EE in maize2. In de-etiolation, enrichment of Z-box, AGL3 binding site, ABRE-like binding site, PIF4 binding site, and BZR1 binding site was significantly biased toward maize1. Interestingly, enhanced enrichment of several cis-elements in maize2 was also observed, such as the CDC5 binding site in shade-induced gene pairs and the CCA1 binding site in FR-induced gene pairs, suggesting that these cis-elements were not involved in gene-induction in this study, or that they function to suppress gene-induction in SAR and de-etiolation, respectively. Differential distribution of cis-elements was also found for shade-suppressed and FR-suppressed genes (Table 1). Taken together, these results indicated that fractionation of cis-elements contributed to the differential responses of the two subgenomes of maize to light signals.
In this study, maize SAR and FRc-mediated de-etiolation was transcriptionally profiled, leading to the identification of important biological processes and cis-elements involved in the response to the two treatments. Co-expression network analysis identified transcription factor modules that may serve as pivotal regulators in SAR and de-etiolation. Cross-species comparison of transcriptional regulation in SAR between Arabidopsis and maize identified shade-induced auxin biosynthesis as a possible target of domestication. Moreover, it was found that the distribution of shade- and FR-responsive genes was biased toward the maize1 subgenome, and this was associated with differential retention of light-responsive cis-elements between the two subgenomes, indicating that genome fractionation can occur both on genes and non-coding sequences with regulatory functions.
Plant material for RNA-seq experiments
To dissect the SAR in maize, the plant growth condition and treatment is the same as that used before to study Arabidopsis SAR . Briefly, B73 maize plants were grown in a light-emitting diode growth chamber (Percival Scientific, Perry, IA) at 12 h light/12 h dark cycles at 24 °C under the following light outputs: Blue, 15 umol/m2s; Red, 96 umol/m2s; Far-red, 21 umol/m2s. On the fourth day after germination (at V1 stage), shoots were harvested above the upper node of the mesocotyl 3 h after lights were on. Then light outputs were immediately adjusted to: Blue, 15 umol/m2s; Red, 12 umol/m2s; Far-red, 105 umol/m2s. Samples were then harvested 1 h, 3 h, and 6 h after the onset of the simulated shade treatment. For each time point, six plants were pooled for each of the three biological replicates. To dissect the de-etiolation under FR, the plant growth condition and treatment is similar to that used in two previous studies [47, 61]. Briefly, B73 plants were grown in darkness at 24 °C. Four days after germination (etiolated V1 stage seedlings), shoots were harvested above the upper node of the mesocotyl. Remaining seedlings were immediately irradiated with 6 umol/m2s FR. Samples were then harvested after 1 h, 3 h, and 6 h treatment. For each time point, six plants were pooled for each of the three biological replicates.
RNA-seq library construction, sequencing and analysis
RNA-seq libraries were generated from 2–5 μg total RNA and size-selected for a 250–300 bp insert for paired-end (PE) sequencing (100 bp for each end). Libraries were quantified on an Agilent bioanalyzer (Agilent, Santa Clara, CA) and sequenced using the Illumina HiSeq2000 system using standard Illumina protocols (Illumina Inc., San Diego, CA). The raw reads were trimmed and filtered using Fastx-toolkit (version 0.0.13). Tophat (version 2.0.12) was used to align reads to the maize reference genome (AGP v3.23). Gene-level expression values are represented by Reads Per Kilobase per Million (RPKM). Differentially expressed genes were identified using the exact test of Robinson and Smyth for two-group comparisons implemented in the EdgeR package . A False Discovery Rate (FDR) corrected p-value < =0.05 and a threshold fold change > =2 were used to call differentially transcribed genes.
Total RNA was isolated using Trizol (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions. Isolated RNA was treated with RNase-free DNase and purified using RNeasy mini columns (Qiagen, Hilden, Germany). First-strand cDNA was synthesized using SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA). The tublin transcript was used as an internal control to normalize the RNA quantity. Three biological replicates were included in quantitative RT–PCR analysis.
Maize gene annotations using GO terms were acquired using the BLAST2GO pipeline, combined with existent GO annotation retrieved from Gramene BioMart. Maize genes were also assigned with MapMan bins based on their homology with rice and Arabidopsis genes using the Mercator pipeline .
Cluster and functional enrichment analysis
Cluster analysis was applied to genes showing differential expression in at least one comparison in SAR or de-etiolation. Log2-transformed ratios were subjected to cluster analysis. For K-means clustering, the K-Means/K-Medians Support (KMS) module in the MEV program (http://www.tm4.org/mev.html) was used for clustering of DTGs. The number of clusters was determined by the Figures of Merit (FOM) module. Genes in each cluster were then classified into Mapman functional categories, and the Fisher’s exact test was used to test for enrichment of each functional category. FDR is controlled by Benjamini and Hochberg’s procedure. Differentially expressed transcription factors were clustered using the Self Organizing Tree Algorithm (SOTA) embedded in the MEV program.
Cis-regulatory motif enrichment analysis
To determine enriched cis-regulatory elements on promoters of co-expressed genes, proximal promoter was defined as 2 kbp upstream and 500 bp downstream of transcription start site (TSS) since this region is adequate to capture the 5′ UTR and the first intron of most maize gene models. A comprehensive collection of plant position weight matrices (PWMs)  and also PWMs deposited in the JASPAR database were used. Promoters were scanned for significantly enriched cis-elements using the PWMEnrich package.
The expression correlation network was calculated with Pearson correlation values between genes as edge weights. The resulting weighted network was clustered with the Markov cluster algorithm (MCL) using clusterMaker2 with granularity parameter 2 and default advanced settings . The graph was visualized in Cytoscape version 3.2.1.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the SRA repository (accession SRP068070).
chalcone flavanone isomerase
- COP1 :
differentially transcribed gene
differentially transcribed TF
- ETC3 :
enhancer of try and CPC3
false discovery rate
- FKF1 :
flavin-binding, kelch repeat, F-BOX1
- GL3 :
- HFR1 :
long hypocotyl in far-red1
- HY5 :
- HYH :
- LAF1 :
long after far-red light1
Markov cluster algorithm
- OZF :
oxidation-related zinc finger
- PAP1 :
photosynthetically active radiation
- PIF :
position weight matrices
reads per kilobase per million
- RSM2 :
radialis-like SANT/MYB 2
shade avoidance response
self organization tree algorithm
- TAA1 :
tryptophan aminotransferase of arabidopsis1
- tb1 :
- tga1 :
teosinte glume architecture1
transcription start site
unfolded protein response
Kebrom TH, Brutnell TP. The molecular analysis of the shade avoidance syndrome in the grasses has begun. J Exp Bot. 2007;58(12):3079–89.
Duvick D. What is yield? In: Edmeades GO, Banziger M, Mickelson HR, Pena-Valdivia CB, editors. Developing drought- and low N-tolerant maize. El Batan: CIMMYT: Proceedings of CIMMYT Symposium, 25–29 March 1996; 1997. p. 332–5.
Smith H, Whitelam GC. The shade avoidance syndrome: multiple responses mediated by multiple phytochromes. Plant Cell Environ. 1997;20(6):840–4.
Ballare CL. Light regulation of plant defense. Annu Rev Plant Biol. 2014;65:335–63.
Shyu C, Brutnell TP. Growth-defence balance in grass biomass production: the role of jasmonates. J Exp Bot. 2015.
Ugarte CC, Trupkin SA, Ghiglione H, Slafer G, Casal JJ. Low red/far-red ratios delay spike and stem growth in wheat. J Exp Bot. 2010;61(11):3151–62.
Robson F, Okamoto H, Patrick E, Harris SR, Wasternack C, Brearley C, Turner JG. Jasmonate and phytochrome A signaling in Arabidopsis wound and shade responses are integrated through JAZ1 stability. Plant Cell. 2010;22(4):1143–60.
Weijschede J, Martinkova J, de Kroon H, Huber H. Shade avoidance in Trifolium repens: costs and benefits of plasticity in petiole length and leaf size. New Phytol. 2006;172(4):655–66.
Schmitt J, Stinchcombe JR, Heschel MS, Huber H. The adaptive evolution of plasticity: phytochrome-mediated shade avoidance responses. Integr Comp Biol. 2003;43(3):459–69.
Wang H, Wang H. Phytochrome signaling: time to tighten up the loose ends. Mol Plant. 2015;8(4):540–51.
Franklin KA, Praekelt U, Stoddart WM, Billingham OE, Halliday KJ, Whitelam GC. Phytochromes B, D, and E act redundantly to control multiple physiological responses in Arabidopsis. Plant Physiol. 2003;131(3):1340–6.
Pacin M, Legris M, Casal JJ. COP1 re-accumulates in the nucleus under shade. Plant J. 2013;75(4):631–41.
Seo HS, Yang JY, Ishikawa M, Bolle C, Ballesteros ML, Chua NH. LAF1 ubiquitination by COP1 controls photomorphogenesis and is stimulated by SPA1. Nature. 2003;423(6943):995–9.
Jang IC, Yang JY, Seo HS, Chua NH. HFR1 is targeted by COP1 E3 ligase for post-translational proteolysis during phytochrome A signaling. Genes Dev. 2005;19(5):593–602.
Hardtke CS, Gohda K, Osterlund MT, Oyama T, Okada K, Deng XW. HY5 stability and activity in arabidopsis is regulated by phosphorylation in its COP1 binding domain. EMBO J. 2000;19(18):4997–5006.
Holm M, Ma LG, Qu LJ, Deng XW. Two interacting bZIP proteins are direct targets of COP1-mediated control of light-dependent gene expression in Arabidopsis. Genes Dev. 2002;16(10):1247–59.
Chen F, Li B, Li G, Charron JB, Dai M, Shi X, et al. Arabidopsis phytochrome A directly targets numerous promoters for individualized modulation of genes in a wide range of pathways. Plant Cell. 2014;26(5):1949–66.
Martinez-Garcia JF, Huq E, Quail PH. Direct targeting of light signals to a promoter element-bound transcription factor. Science. 2000;288(5467):859–63.
Lorrain S, Allen T, Duek PD, Whitelam GC, Fankhauser C. Phytochrome-mediated inhibition of shade avoidance involves degradation of growth-promoting bHLH transcription factors. Plant J. 2008;53(2):312–23.
Shen H, Zhu L, Castillon A, Majee M, Downie B, Huq E. Light-induced phosphorylation and degradation of the negative regulator PHYTOCHROME-INTERACTING FACTOR1 from Arabidopsis depend upon its direct physical interactions with photoactivated phytochromes. Plant Cell. 2008;20(6):1586–602.
Shen Y, Khanna R, Carle CM, Quail PH. Phytochrome induces rapid PIF5 phosphorylation and degradation in response to red-light activation. Plant Physiol. 2007;145(3):1043–51.
Leivar P, Monte E, Al-Sady B, Carle C, Storer A, Alonso JM, Ecker JR, Quail PH. The Arabidopsis phytochrome-interacting factor PIF7, together with PIF3 and PIF4, regulates responses to prolonged red light by modulating phyB levels. Plant Cell. 2008;20(2):337–52.
Li L, Ljung K, Breton G, Schmitz RJ, Pruneda-Paz J, Cowing-Zitron C, Cole BJ, Ivans LJ, Pedmale UV, Jung HS,et al. Linking photoreceptor excitation to changes in plant architecture. Genes Dev. 2012;26(8):785–90.
Al-Sady B, Ni WM, Kircher S, Schafer E, Quail PH. Photoactivated phytochrome induces rapid PIF3 phosphorylation prior to proteasorne-mediated degradation. Mol Cell. 2006;23(3):439–46.
Kumar SV, Lucyshyn D, Jaeger KE, Alos E, Alvey E, Harberd NP,Wigge PA. Transcription factor PIF4 controls the thermosensory activation of flowering. Nature. 2012;484(7393):242–5.
Leivar P, Quail PH. PIFs: pivotal components in a cellular signaling hub. Trends Plant Sci. 2011;16(1):19–28.
Leivar P, Monte E. PIFs: systems integrators in plant development. Plant Cell. 2014;26(1):56–78.
Zhou P, Song M, Yang Q, Su L, Hou P, Guo L, et al. Both PHYTOCHROME RAPIDLY REGULATED1 (PAR1) and PAR2 promote seedling photomorphogenesis in multiple light signaling pathways. Plant Physiol. 2014;164(2):841–52.
Bou-Torrent J, Roig-Villanova I, Galstyan A, Martinez-Garcia JF. PAR1 and PAR2 integrate shade and hormone transcriptional networks. Plant Signal Behav. 2008;3(7):453–4.
Shi H, Zhong S, Mo X, Liu N, Nezames CD, Deng XW. HFR1 sequesters PIF1 to govern the transcriptional network underlying light-initiated seed germination in Arabidopsis. Plant Cell. 2013;25(10):3770–84.
Hornitschek P, Lorrain S, Zoete V, Michielin O, Fankhauser C. Inhibition of the shade avoidance response by formation of non-DNA binding bHLH heterodimers. EMBO J. 2009;28(24):3893–902.
Sessa G, Carabelli M, Sassi M, Ciolfi A, Possenti M, Mittempergher F, Becker J, Morelli G, Ruberti I. A dynamic balance between gene activation and repression regulates the shade avoidance response in Arabidopsis. Genes Dev. 2005;19(23):2811–5.
Sheehan MJ, Farmer PR, Brutnell TP. Structure and expression of maize phytochrome family homeologs. Genetics. 2004;167(3):1395–405.
Sheehan MJ, Kennedy LM, Costich DE, Brutnell TP. Subfunctionalization of PhyB1 and PhyB2 in the control of seedling and mature plant traits in maize. Plant J. 2007;49(2):338–53.
Dubois PG, Olsefski GT, Flint-Garcia S, Setter TL, Hoekenga OA, Brutnell TP. Physiological and genetic characterization of end-of-day far-red light response in maize seedlings. Plant Physiol. 2010;154(1):173–86.
Maddonni GA, Otegui ME, Andrieu B, Chelle M, Casal JJ. Maize leaves turn away from neighbors. Plant Physiol. 2002;130(3):1181–9.
Studer A, Zhao Q, Ross-Ibarra J, Doebley J. Identification of a functional transposon insertion in the maize domestication gene tb1. Nat Genet. 2011;43(11):1160–3.
Wang H, Nussbaum-Wagler T, Li B, Zhao Q, Vigouroux Y, Faller M, Bomblies K, Lukens L, Doebley JF. The origin of the naked grains of maize. Nature. 2005;436(7051):714–9.
Garg AK, Sawers RJ, Wang H, Kim JK, Walker JM, Brutnell TP, Parthasarathy MV, Vierstra RD, Wu RJ. Light-regulated overexpression of an Arabidopsis phytochrome A gene in rice alters plant architecture and increases grain yield. Planta. 2006;223(4):627–36.
Warnasooriya SN, Brutnell TP. Enhancing the productivity of grasses under high-density planting by engineering light responses: from model systems to feedstocks. J Exp Bot. 2014;65(11):2825–34.
Boccalandro HE, Ploschuk EL, Yanovsky MJ, Sanchez RA, Gatz C, Casal JJ. Increased phytochrome B alleviates density effects on tuber yield of field potato crops. Plant Physiol. 2003;133(4):1539–46.
Thiele A, Herold M, Lenk I, Quail PH, Gatz C. Heterologous expression of Arabidopsis phytochrome B in transgenic potato influences photosynthetic performance and tuber development. Plant Physiol. 1999;120(1):73–82.
Robson PR, McCormac AC, Irvine AS, Smith H. Genetic engineering of harvest index in tobacco through overexpression of a phytochrome gene. Nat Biotechnol. 1996;14(8):995–8.
Rao AQ, Irfan M, Saleem Z, Nasir IA, Riazuddin S, Husnain T. Overexpression of the phytochrome B gene from Arabidopsis thaliana increases plant growth and yield of cotton (Gossypium hirsutum). J Zhejiang Univ Sci B. 2011;12(4):326–34.
Schittenhelm S, Menge-Hartmann U, Oldenburg E. Photosynthesis, carbohydrate metabolism, and yield of phytochrome-B-overexpressing potatoes under different light regimes. Crop Sci. 2004;44(1):131–43.
Ciolfi A, Sessa G, Sassi M, Possenti M, Salvucci S, Carabelli M, Morelli G, Ruberti I. Dynamics of the shade-avoidance response in Arabidopsis. Plant Physiol. 2013;163(1):331–53.
Roig-Villanova I, Bou J, Sorin C, Devlin PF, Martinez-Garcia JF. Identification of primary target genes of phytochrome signaling. Early transcriptional control during shade avoidance responses in Arabidopsis. Plant Physiol. 2006;141(1):85–96.
Devlin PF, Yanovsky MJ, Kay SA. A genomic analysis of the shade avoidance response in Arabidopsis. Plant Physiol. 2003;133(4):1617–29.
Li J, Li G, Wang H, Wang Deng X. Phytochrome signaling mechanisms. Arabidopsis Book. 2011;9, e0148.
Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37(6):914–39.
Sun Y, Fan XY, Cao DM, Tang W, He K, Zhu JY, He JX, Bai MY, Zhu S, Oh E, et al. Integration of brassinosteroid signal transduction with the transcription network for plant growth regulation in Arabidopsis. Dev Cell. 2010;19(5):765–77.
Sellaro R, Yanovsky MJ, Casal JJ. Repression of shade-avoidance reactions by sunfleck induction of HY5 expression in Arabidopsis. Plant J. 2011;68(5):919–28.
Nozue K, Tat AV, Kumar Devisetty U, Robinson M, Mumbach MR, Ichihashi Y, Lekkala S, Maloof JN. Shade avoidance components and pathways in adult plants revealed by phenotypic profiling. PLoS Genet. 2015;11(4), e1004953.
Gross BL, Olsen KM. Genetic perspectives on crop domestication. Trends Plant Sci. 2010;15(9):529–37.
Ferrier T, Matus JT, Jin J, Riechmann JL. Arabidopsis paves the way: genomic and network analyses in crops. Curr Opin Biotechnol. 2011;22(2):260–70.
Turchi L, Carabelli M, Ruzza V, Possenti M, Sassi M, Penalosa A, Sessa G, Salvi S, Forte V, Morelli G, et al. Arabidopsis HD-Zip II transcription factors control apical embryo development and meristem function. Development. 2013;140(10):2118–29.
Kunihiro A, Yamashino T, Nakamichi N, Niwa Y, Nakanishi H, Mizuno T. Phytochrome-interacting factor 4 and 5 (PIF4 and PIF5) activate the homeobox ATHB2 and auxin-inducible IAA29 genes in the coincidence mechanism underlying photoperiodic control of plant growth of Arabidopsis thaliana. Plant Cell Physiol. 2011;52(8):1315–29.
Rubin G, Tohge T, Matsuda F, Saito K, Scheible WR. Members of the LBD family of transcription factors repress anthocyanin synthesis and affect additional nitrogen responses in Arabidopsis. Plant Cell. 2009;21(11):3567–84.
Nemie-Feyissa D, Olafsdottir SM, Heidari B, Lillo C. Nitrogen depletion and small R3-MYB transcription factors affecting anthocyanin accumulation in Arabidopsis leaves. Phytochemistry. 2014;98:34–40.
Wang S, Hubbard L, Chang Y, Guo J, Schiefelbein J, Chen JG. Comprehensive analysis of single-repeat R3 MYB proteins in epidermal cell patterning and their transcriptional regulation in Arabidopsis. BMC Plant Biol. 2008;8:81.
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 U S A. 2001;98(16):9437–42.
Chang CS, Maloof JN, Wu SH. COP1-mediated degradation of BBX22/LZF1 optimizes seedling development in Arabidopsis. Plant Physiol. 2011;156(1):228–39.
Crocco CD, Holm M, Yanovsky MJ, Botto JF. AtBBX21 and COP1 genetically interact in the regulation of shade avoidance. Plant J. 2010;64(4):551–62.
Crocco CD, Locascio A, Escudero CM, Alabadi D, Blazquez MA, Botto JF. The transcriptional regulator BBX24 impairs DELLA activity to promote shade avoidance in Arabidopsis thaliana. Nat Commun. 2015;6:6202.
Gangappa SN, Crocco CD, Johansson H, Datta S, Hettiarachchi C, Holm M, Botto JF. The Arabidopsis B-BOX protein BBX25 interacts with HY5, negatively regulating BBX22 expression to suppress seedling photomorphogenesis. Plant Cell. 2013;25(4):1243–57.
Crocco CD, Holm M, Yanovsky MJ, Botto JF. Function of B-BOX under shade. Plant Signal Behav. 2011;6(1):101–4.
Datta S, Hettiarachchi GH, Deng XW, Holm M. Arabidopsis CONSTANS-LIKE3 is a positive regulator of red light signaling and root growth. Plant Cell. 2006;18(1):70–84.
Holtan HE, Bandong S, Marion CM, Adam L, Tiwari S, Shen Y, Maloof JN, Maszle DR, Ohto MA, Preuss S, et al. BBX32, an Arabidopsis B-Box protein, functions in light signaling by suppressing HY5-regulated gene expression and interacting with STH2/BBX21. Plant Physiol. 2011;156(4):2109–23.
Ruckle ME, Burgoon LD, Lawrence LA, Sinkler CA, Larkin RM. Plastids are major regulators of light signaling in Arabidopsis. Plant Physiol. 2012;159(1):366–90.
Khanna R, Shen Y, Toledo-Ortiz G, Kikis EA, Johannesson H, Hwang YS, Quail PH. Functional profiling reveals that only a small number of phytochrome-regulated early-response genes in Arabidopsis are necessary for optimal deetiolation. Plant Cell. 2006;18(9):2157–71.
Oh S, Montgomery BL. Phytochrome-induced SIG2 expression contributes to photoregulation of phytochrome signalling and photomorphogenesis in Arabidopsis thaliana. J Exp Bot. 2013;64(18):5457–72.
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(5):728–42.
Fornara F, Panigrahi KC, Gissot L, Sauerbrunn N, Ruhl M, Jarillo JA, Coupland G. Arabidopsis DOF transcription factors act redundantly to reduce CONSTANS expression and are essential for a photoperiodic flowering response. Dev Cell. 2009;17(1):75–86.
Sawa M, Nusinow DA, Kay SA, Imaizumi T. FKF1 and GIGANTEA complex formation is required for day-length measurement in Arabidopsis. Science. 2007;318(5848):261–5.
Guan P, Wang R, Nacry P, Breton G, Kay SA, Pruneda-Paz JL, Davani A, Crawford NM. Nitrate foraging by Arabidopsis roots is mediated by the transcription factor TCP20 through the systemic signaling pathway. Proc Natl Acad Sci U S A. 2014;111(42):15267–72.
Baud S, Wuilleme S, To A, Rochat C, Lepiniec L. Role of WRINKLED1 in the transcriptional regulation of glycolytic and fatty acid biosynthetic genes in Arabidopsis. Plant J. 2009;60(6):933–47.
Guo HS, Xie Q, Fei JF, Chua NH. MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for arabidopsis lateral root development. Plant Cell. 2005;17(5):1376–86.
Mu RL, Cao YR, Liu YF, Lei G, Zou HF, Liao Y, et al. An R2R3-type transcription factor gene AtMYB59 regulates root growth and cell cycle progression in Arabidopsis. Cell Res. 2009;19(11):1291–304.
Lu D, Wang T, Persson S, Mueller-Roeber B, Schippers JH. Transcriptional control of ROS homeostasis by KUODA1 regulates cell expansion during leaf development. Nat Commun. 2014;5:3767.
Hamaguchi A, Yamashino T, Koizumi N, Kiba T, Kojima M, Sakakibara H, Mizuno T. A small subfamily of Arabidopsis RADIALIS-LIKE SANT/MYB genes: a link to HOOKLESS1-mediated signal transduction during early morphogenesis. Biosci Biotechnol Biochem. 2008;72(10):2687–96.
Huang P, Chung MS, Ju HW, Na HS, Lee DJ, Cheong HS, Kim CS. Physiological characterization of the Arabidopsis thaliana oxidation-related zinc finger 1, a plasma membrane protein involved in oxidative stress. J Plant Res. 2011;124(6):699–705.
Huang P, Ju HW, Min JH, Zhang X, Chung JS, Cheong HS, Kim CS. Molecular and physiological characterization of the Arabidopsis thaliana oxidation-related zinc finger 2, a plasma membrane protein involved in ABA and salt stress response through the ABI2-mediated signaling pathway. Plant Cell Physiol. 2012;53(1):193–203.
Song ZT, Sun L, Lu SJ, Tian Y, Ding Y, Liu JX. Transcription factor interaction with COMPASS-like complex regulates histone H3K4 trimethylation for specific gene expression in plants. Proc Natl Acad Sci U S A. 2015;112(9):2900–5.
Hsieh HL, Okamoto H. Molecular interaction of jasmonate and phytochrome A signalling. J Exp Bot. 2014;65(11):2847–57.
Gangappa SN, Prasad VB, Chattopadhyay S. Functional interconnection of MYC2 and SPA1 in the photomorphogenic seedling development of Arabidopsis. Plant Physiol. 2010;154(3):1210–9.
Chico JM, Fernandez-Barbero G, Chini A, Fernandez-Calvo P, Diez-Diaz M, Solano R. Repression of jasmonate-dependent defenses by shade involves differential regulation of protein stability of MYC transcription factors and their JAZ repressors in Arabidopsis. Plant Cell. 2014;26(5):1967–80.
Li T, Jia KP, Lian HL, Yang X, Li L, Yang HQ. Jasmonic acid enhancement of anthocyanin accumulation is dependent on phytochrome A signaling pathway under far-red light in Arabidopsis. Biochem Biophys Res Commun. 2014;454(1):78–83.
Mehrtens F, Kranz H, Bednarek P, Weisshaar B. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005;138(2):1083–96.
Stracke R, Favory JJ, Gruber H, Bartelniewoehner L, Bartels S, Binkert M, Funk M, Weisshaar B, Ulm R. The Arabidopsis bZIP transcription factor HY5 regulates expression of the PFG1/MYB12 gene in response to light and ultraviolet-B radiation. Plant Cell Environ. 2010;33(1):88–103.
Zhang Y, Zheng S, Liu Z, Wang L, Bi Y. Both HY5 and HYH are necessary regulators for low temperature-induced anthocyanin accumulation in Arabidopsis seedlings. J Plant Physiol. 2011;168(4):367–74.
Shi H, Wang X, Ye T, Cheng F, Deng J, Yang P, Zhang Y, Chan Z. The Cys2/His2-type zinc finger transcription factor ZAT6 modulates biotic and abiotic stress responses by activating salicylic acid-related genes and CBFs in Arabidopsis. Plant Physiol. 2014.
Jiang Y, Yang B, Deyholos MK. Functional characterization of the Arabidopsis bHLH92 transcription factor in abiotic stress. Mol Genet Genomics. 2009;282(5):503–16.
Wang X, Culver JN. DNA binding specificity of ATAF2, a NAC domain transcription factor targeted for degradation by Tobacco mosaic virus. BMC Plant Biol. 2012;12:157.
Geilen K, Bohmer M. Dynamic subnuclear relocalisation of WRKY40 in response to Abscisic acid in Arabidopsis thaliana. Sci Rep. 2015;5:13369.
Weller JL, Hecht V, Vander Schoor JK, Davidson SE, Ross JJ. Light regulation of gibberellin biosynthesis in pea is mediated through the COP1/HY5 pathway. Plant Cell. 2009;21(3):800–13.
Mizoguchi T, Wright L, Fujiwara S, Cremer F, Lee K, Onouchi H, Mouradov A, Fowler S, Kamada H, Putterill J et al: Distinct roles of GIGANTEA in promoting flowering and regulating circadian rhythms in Arabidopsis. Plant Cell. 2005;17(8):2255–70.
Song YH, Estrada DA, Johnson RS, Kim SK, Lee SY, MacCoss MJ, Imaizumi T. Distinct roles of FKF1, Gigantea, and Zeitlupe proteins in the regulation of Constans stability in Arabidopsis photoperiodic flowering. Proc Natl Acad Sci U S A. 2014;111(49):17672–7.
Song YH, Smith RW, To BJ, Millar AJ, Imaizumi T. FKF1 conveys timing information for CONSTANS stabilization in photoperiodic flowering. Science. 2012;336(6084):1045–9.
Imaizumi T, Schultz TF, Harmon FG, Ho LA, Kay SA. FKF1 F-box protein mediates cyclic degradation of a repressor of CONSTANS in Arabidopsis. Science. 2005;309(5732):293–7.
Sawa M, Kay SA. GIGANTEA directly activates flowering locus T in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2011;108(28):11698–703.
Tao Y, Ferrer JL, Ljung K, Pojer F, Hong F, Long JA, et al. Rapid synthesis of auxin via a new tryptophan-dependent pathway is required for shade avoidance in plants. Cell. 2008;133(1):164–76.
Sun J, Qi L, Li Y, Chu J, Li C. PIF4-mediated activation of YUCCA8 expression integrates temperature into the auxin pathway in regulating arabidopsis hypocotyl growth. PLoS Genet. 2012;8(3), e1002594.
Mashiguchi K, Tanaka K, Sakai T, Sugawara S, Kawaide H, Natsume M, Hanada A, Yaeno T, 963 Shirasu K, Yao H, et al. The main auxin biosynthesis pathway in Arabidopsis. Proc Natl Acad Sci U S A. 2011;108(45):18512–7.
Schnable JC, Springer NM, Freeling M. Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc Natl Acad Sci U S A. 2011;108(10):4069–74.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Lohse M, Nagel A, Herter T, May P, Schroda M, Zrenner R, Tohge T, Fernie AR, Stitt M, Usadel B. Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ. 2014;37(5):1250–8.
Eveland AL, Goldshmidt A, Pautler M, Morohashi K, Liseron-Monfils C, Lewis MW, Kumari S, Hiraga S, Yang F, Unger-Wallace E, et al. Regulatory modules controlling maize inflorescence architecture. Genome Res. 2014;24(3):431–43.
Morris JH, Apeltsin L, Newman AM, Baumbach J, Wittkop T, Su G, Bader GD, Ferrin TE. clusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinformatics. 2011;12:436.
We thank Dr. Sunita Kumari and Dr. David Jackson at Cold Spring Harbor Laboratory for kindly providing a comprehensive collection of PWMs of cis-elements in plants. This work was supported by National Natural Science Foundation of China grant 91435112 (to Chunyi Zhang) and 31430008 (to Haiyang Wang). No conflict of interest was declared.
The authors declare that they have no competing interests.
Author contributions are as follows: HW, HYW, CYZ and ZHL designed the experiments. HW collected tissues and extracted RNA. HW built the analysis pipeline and analyzed the sequence data. GXW, BBZ and BBW verified gene expression data by qRT-PCR. HW, HYW and CYZ interpreted the results of all analysis and wrote the manuscript. All authors have read and approved the manuscript.
Maize seedlings (inbred line B73) were grown under high R/FR and low R/FR (simulated shade) for 10 days. Note that seedlings grown under the simulated shade were with longer leaves and accumulated less anthocyanin. (PDF 131 kb)
Principle Component Analysis (PCA) of samples using genome-wide gene expression values. (PDF 192 kb)
K-means clustering of DTGs. (XLSX 298 kb)
Confirmation of RNA-Seq data by qRT-PCR. (XLSX 18 kb)
Enriched cis-elements in each DTG cluster. (XLSX 29 kb)
A summary of differentially transcribed transcription factors. (XLSX 28 kb)
Phylogenetic analysis of a group of early shade-induced Arabidopsis homeodomain-leucine zipper genes and their two putative maize orthologs. The evolutionary tree was inferred with amino acid sequences using the Neighbor-Joining method implemented in MEGA6. Bootstrap values (with 100 replicates) were shown next to the branches. The tree was rooted on a close homolog in P. patens. (PDF 100 kb)
Phylogenetic analysis showed a single maize ortholog of the early shade-inducible Arabidopsis atypical basic helix-loop-helix (bHLH) genes PAR1 and PAR2. The evolutionary tree was inferred with amino acid sequences using the Neighbor-Joining method implemented in MEGA6. Bootstrap values (with 100 replicates) were shown next to the branches. The tree was rooted on a close homolog in P. patens. (PDF 97 kb)
Phylogeny and structure of putative maize PIFs. Seven maize bHLH family transcription factors were identified as PIFs, due to their close relationship with Arabidopsis PIFs in the phylogenetic tree (a), and also the conserved ABP motif shared among all PIFs responsible for their physical interaction with phytochrome B (b). Phylogenetic analysis was conducted by MEGA6 using the neighbor-joining method. Bootstrap values were obtained using 100 bootstrap replicates and are shown next to the branches. The tree was rooted on the closest homolog of Arabidopsis PIFs in P. patens. Note that HFR1 and PIL1 are on a branch (shaded in blue) without any maize orthologs. (PDF 129 kb)
Conserved inter-species expression patterns in SAR. (XLSX 11 kb)
Contrasting inter-species expression patterns in SAR. (XLSX 10 kb)
Enriched GO terms for shade-responsive DTGs. (XLSX 13 kb)
A comparison of light-responsive genes between the two subgenomes of maize. (XLSX 24 kb)
About this article
Cite this article
Wang, H., Wu, G., Zhao, B. et al. Regulatory modules controlling early shade avoidance response in maize seedlings. BMC Genomics 17, 269 (2016). https://doi.org/10.1186/s12864-016-2593-6