- Open Access
Unraveling multifaceted contributions of small regulatory RNAs to photomorphogenic development in Arabidopsis
BMC Genomicsvolume 18, Article number: 559 (2017)
Post-transcriptional control of gene expression mediated by small regulatory RNAs (sRNAs) is vital for growth and development of diverse organisms. The biogenesis of sRNAs is regulated by both positive and negative regulators known to regulate photomorphogenic development. Two microRNAs (miRNAs), miR157 and miR319, also regulate photomorphogenesis. However, genome-wide profiling of sRNAs and their regulation of target genes during photomorphogenesis has been missing. We provide a comprehensive view of sRNA-controlled gene expression in this developmental process.
By profiling sRNAs and the 5′ ends of degraded mRNAs during the first 24 h of photomorphogenic development in Arabidopsis, we identified 335 sRNA-mediated mRNA cleavage events in de-etiolating seedlings. These cleavage events are primarily resulted from actions of highly expressed miRNAs and irrelevant to the abundance of target mRNAs. In the light, the expression of the slicer protein gene ARGONAUTE1 in the miRNA functioning pathway could be fine-tuned by miRNA168a/b. We also found that miR396a/b positively regulates de-etiolation by suppressing GROWTH REGULATING FACTORs. Our results suggest that the miRNAs are required to tune down the target mRNAs and regulate photomorphogenesis.
sRNAs may have a broad impact on gene expression regulation for optimized photomorphogenic development. With both positive and negative regulators under the control of sRNAs, young Arabidopsis seedlings can have a timely but not exaggerated developmental adaptation to light.
Plants have evolved a plethora of morphological alterations to adapt to their surroundings. Photomorphogenesis, or de-etiolation, is one such process when seedlings first experience light irradiation. The rate of hypocotyl elongation decreases in seedlings under light exposure, which allows for the formation of firm structural support for seedlings emerging from the soil surface. Also, the cotyledons open and expand to maximize the area of light perception and to expose the shoot apical meristem for the development of true leaves. Light also triggers the development of chloroplasts for photosynthesis so that plants can utilize light energy for autotrophic growth and development [1,2,3].
In Arabidopsis, photomorphogenesis is under the control of at least three types of photoreceptors, including the red (R)–far-red (FR) light photoreceptor phytochromes (phys), blue light (B) photoreceptor cryptochromes (crys) and the UV-B photoreceptor, UVR8 [3,4,5,6,7,8]. The perceived light signals trigger signaling cascades that reprogram gene expression for photomorphogenic development. Transcriptional profiling for Arabidopsis seedlings exposed to B, FR, R light and the light–dark transition have revealed differential expression of approximately one-third of the genome . The light-regulated genome-wide transcriptomic adjustment requires the actions of transcription factors. One of the most well-characterized transcription factors conveying light signals to changes of gene expression is ELONGATED-HYPOCOTYL 5 (HY5). HY5 is a light-regulated bZIP transcription factor that upregulates the expression of many light-responsive genes during de-etiolation . In addition to activating transcription, light also enhances the translational efficiency of thousands of genes, especially those committed to the translation apparatus and chloroplast functions [11, 12].
Plant small regulatory RNAs (sRNAs) are 20 to 24 nt long and can be classified into microRNAs (miRNAs) and small interfering RNAs (siRNAs) primarily according to different modes of biogenesis. MiRNAs originate from stem-loop structures of primary transcripts, and siRNAs are mostly derived from double-stranded RNAs . Phased siRNAs (phasiRNAs) are a special group of siRNAs generated from mRNAs cleaved by 22-nt miRNAs or siRNAs [14,15,16]. Plant miRNAs can mediate the cleavage or translation inhibition of target mRNAs, whereas siRNAs function via RNA-dependent DNA methylation (RdDM) for transcriptional gene silencing or post-transcriptional target mRNA cleavage [17,18,19,20,21,22].
Previous studies have implied that sRNAs are involved in gene expression regulation during de-etiolation. Mutants defective in genes for miRNA biogenesis and functions have altered light responses. For example, in Arabidopsis, light hypersensitive phenotypes have been observed to carry mutations in the miRNA processor HYPONASTIC LEAVES 1 (HYL1), the sRNA methyltransferase HUA ENHANCER1 (HEN1), the sRNA transporter HASTY (HST), and the slicer protein ARGONAUTE1 (AGO1) [23,24,25]. A negative regulator of photomorphogenesis, CONSTITUTIVE PHOTOMORPHOGENESIS 1 (COP1), can protect HYL1 against degradation, thereby leading to a stabilized miRNA pool . Transcripts of the positive regulator HY5 and negative regulator TEOSINTE BRANCHED 1, CYCLOIDEA AND PCF TRANSCRIPTION FACTORs (TCPs) of photomorphogenesis were shown to be under regulation by miR157d and miR319, respectively . In addition, HY5 was found to bind to promoter regions of at least 8 miRNAs (MIR) loci and required for the accumulation of miR156d, miR402, miR408, miR775 and miR858 . These studies provide a glimpse into photomorphogenic development shaped by the actions of a few sRNAs. A global investigation of sRNAs and their targets would greatly help in assessing the impact of post-transcriptional regulation in photomorphogenic development. However, such information is currently missing in de-etiolating seedlings.
In this study, we profiled sRNAs at 6 times during the first 24 h of Arabidopsis photomorphogenic development. We also sequenced 5′ ends of degraded mRNAs (degradome) in both dark- and light-grown seedlings to reveal sRNA-mediated cleavage of mRNAs during de-etiolation. Pairwise studies of sRNAs and their target mRNAs indicated that a high sRNA-to-target ratio is a key determinant for successful mRNA target repression by sRNAs. The high ratio is mainly contributed by the abundance of sRNAs. A total of 335 sRNA-target mRNA regulatory pairs were identified in de-etiolating seedlings, with several sRNAs demonstrated to regulate photomorphogenesis. The action of miR168 leads to reduced expression of AGO1 under light, thereby offering a feedback regulation of miRNA functioning during de-etiolation. miR396 were identified to act as positive regulators of photomorphogenesis. In addition, we revealed that some 24-nt siRNAs had potential to cause target cleavage in de-etiolating seedlings. Our data indicate that sRNAs function in multiple regulatory circuits for optimized seedling growth under light illumination.
The expression and actions of small RNAs in de-etiolating Arabidopsis seedlings
We first used deep sequencing to survey the sRNAs in de-etiolating Arabidopsis seedlings. The sRNAs were isolated from 4-d-old dark-grown (W0 h) seedlings and seedlings that were further treated with continuous white light irradiation for 1 to 24 h (W1 h to W24 h) and subjected to deep sequencing (Fig. 1). Approximately 18–22 million reads were obtained for each sample in 3 biological replicates. For each dataset, 94–98% of the filtered reads (see methods) could be mapped to the TAIR10 genome (Additional file 1: Table S1). We first analyzed miRNAs and phasiRNAs, as they are frequently studied groups of plant sRNAs. Among the 427 Arabidopsis miRNAs annotated in miRBase 21 , 207 (48.5%) are considered expressed (see methods for criteria) (Additional file 1: Table S2). Overall 58 phasiRNAs derived from 12 phasiRNA-generating loci (PHAS, or trans-acting siRNA-generating loci, TAS)  were expressed in de-etiolating seedlings (Table 1 and Additional file 1: Table S3). In addition to miRNAs and phasiRNAs, 4255 20–24 nt siRNAs were defined expressed in this developmental stage (Additional file 1: Table S4).
Since both miRNAs and siRNAs can target mRNAs for cleavage [19, 29,30,31], we aimed to identify sRNA–target pairs that may be involved in photomorphogenic development. We performed degradome sequencing followed by CleaveLand analyses  to identify target mRNAs cleaved by expressed sRNAs in de-etiolating seedlings. The degraded mRNA samples were obtained from seedlings grown under dark (W0 h) or light (equally mixed from samples treated with W for 1, 3, 6, 12 and 24 h). Approximately 50 million reads for each of the libraries were obtained; 81–85% could be mapped to TAIR10 cDNAs (Additional file 1: Table S1). Our analyses suggested that 262 non-redundant sRNAs could mediate the cleavage of 306 Arabidopsis mRNAs (a total of 335 target cleavage sites). Among them, 142 cleavage events were mediated by miRNAs, 13 by phasiRNAs and 180 by other siRNAs (Fig. 1; Additional file 1: Table S5 and S6). These newly identified sRNA-target pairs are potential players in post-transcriptional gene expression regulation in Arabidopsis photomorphogenesis.
sRNA abundance determines the likelihood of target mRNA cleavage
When analyzing the mRNA degradomes, we noticed that although 90 miRNA families were expressed in de-etiolating seedlings (Additional file 1: Table S2), target cleavage was detected for only members of 49 miRNA families (Additional file 1: Table S5). The results prompted us to investigate factors affecting target cleavage or the identification of degradome signatures. Previous reports have shown that high target abundance will compromise the repression activity of miRNAs and siRNAs when introduced via transfection in animal cell lines . A computational model based on fixed concentration of miRNAs and mRNAs implicated that the concentration of miRNAs has a greater effect on miRNA-mRNA interaction in Drosophila melanogaster and in human . In contrast to the seed pairing seen in most animal miRNA-mRNA interactions, most plant miRNAs interact with their target mRNAs at high complementarity that leads to cleavage of target mRNAs . The availability of transcriptome data for mRNAs , sRNAs and degradome signatures in this study allowed us to investigate whether miRNAs/siRNAs or target abundance is important for effective miRNA/siRNA-mediated target cleavage in de-etiolating Arabidopsis seedlings.
Our analysis indicated that miRNAs causing target cleavage tended to have higher abundance, as compared with miRNAs that failed to generate detectable target degradome signatures (Fig. 2a). The results were similar for miRNAs expressed in both the dark (p = 5.2 × 10−8, D = 0.4327 in Kolmogorov-Smirnov test, K-S test) and light (p = 3.0 × 10−8, D = 0.4394, K-S test) (Fig. 2a). In contrast, for mRNAs with predicted miRNA target sites, the transcript abundance was comparable for mRNAs with or without degradome signatures identified (Fig. 2b). Pair-wise examination of miRNA-to-target ratios revealed relatively higher ratios for miRNA-mRNA pairs with observed cleavage under both dark (p = 1.3 × 10−12, D = 0.4108, K-S test) and light (p = 2.5 × 10−12, D = 0.4067, K-S test) (Fig. 2c). Therefore, miRNAs with high abundance may give rise to high miRNA-to-target ratios, thereby leading to successful target mRNA cleavage.
The above notion remains true when applied to siRNA-mediated mRNA cleavage (Additional file 2: Figure S1), although not as significant as for miRNAs (Fig. 2a). Among the 4255 expressed siRNAs, only 180 have degradome signatures identified for their target mRNAs (Additional file 1: Table S6), as compared with 155 targets resulting from 265 miRNA/phasiRNA-mediated cleavages (Additional file 1: Table S5), possibly because of the significantly lower expression of most siRNAs than miRNAs (Fig. 2d).
Although the expression of most of the miRNAs remained unchanged before and after light treatment, light appears to down-regulate the expression of target mRNAs with degradome signatures but not that of mRNAs without evidence of cleavage (Fig. 2e). Therefore, instead of regulating miRNA expression, light signals may potentiate the target-cleavage activities of miRNAs to tune down the expression of their target genes during de-etiolation. Whether this reduction is achieved by regulating the expression or enzymatic activities of slicer complexes remains to be investigated.
Light optimizes miRNA functioning via the action of miR168
Previously, we reported a feedback regulation between HY5 and HEN1, which indicates that an sRNA equilibrium is required during photomorphogenic development . Here, we sought to identify whether light regulates steps in addition to HEN1 in sRNA biogenesis and functioning. In Arabidopsis, miR168 targets the sRNA slicer gene AGO1 . Deep sequencing and northern blot results indicated that the expression pattern of miR168 only slightly fluctuated under light (Fig. 3a, b). However, the mean abundance of miR168a/b ranged from 880 to 1270 read per million reads (RPM) (Additional file 1: Table S2), which is more than 10 times greater than the median level of miRNAs overall (Fig. 2a). Thus, miR168a/b has high potential in mediating the cleavage of AGO1 transcript. Indeed, under light, AGO1 cleavage signatures could be detected (Additional file 1: Table S5), which led to the down-regulation of AGO1 (Fig. 3c). The detection of AGO1 cleavage signature only under light is also consistent with preferential light-mediated downregulation of miRNA target mRNAs (Fig. 2e). Thus, miR168a/b have potential to desensitize the sRNA actions by targeting AGO1 for degradation during photomorphogenesis.
Light regulates the expression of some miRNAs and phasiRNAs
Although most miRNA levels were unchanged before and after light treatment in young Arabidopsis seedlings (Fig. 2d; Additional file 1: Table S2), we still observed that 32% (67 of 207) of expressed miRNAs were regulated by light (Fig. 4; Additional file 1: Table S2). Only 8 of 58 expressed phasiRNAs were differentially regulated by light (Additional file 1: Table S3). Because sRNA abundance is a major determinant for target cleavage in seedlings (Fig. 2), any changes in sRNA levels under light may alter their target suppression capacity. Thus, the light regulation of miRNAs and phasiRNAs may provide a timely control of target mRNAs to shape photomorphogenic development. The light responsiveness of the 67 light-regulated miRNAs could be classified into 3 major clusters by k-mean clustering (Fig. 4; Table 2). miR163 belongs to cluster I, whose expression is barely detectable in the dark but is rapidly induced by light. miR163 has been shown to promote seed germination and primary root growth during early seedling development, but not involved in light-induced inhibition of hypocotyl elongation . Cluster II miRNAs are also upregulated by light, and the up regulation is more prominent after prolonged light exposure. The miRNAs in cluster II include miR157d, reported to target HY5 during photomorphogenic development . The miRNAs in cluster III were down regulated by light, especially after 6-h light exposure.
miR396 promote photomorphogenesis by tuning GRF levels
In de-etiolating seedlings, the degradome signatures were most frequently found for mRNAs targeted by members in the miR156/157 and miR396 families (Fig. 5a). miR157d directly targets HY5 to desensitize the light signals during photomorphogenesis , but the functions of miR396 in photomorphogenic development remain obscure.
The expression of miR396 was slightly decreased upon light treatment (Fig. 5b and c). miR396 can target 7 GROWTH REGULATING FACTORs (GRFs) [37, 38], and the cleaved signatures of all 7 GRFs were detected in our degradome analysis (Additional file 1: Table S5). GRF1, GRF2, GRF3 and GRF7 showed relatively higher degradome signature reads amongst the GRFs (Fig. 5d), and all showed clear down regulation under light (Fig. 5e). GRF1, GRF2 and GRF3 cooperatively regulate leaf and cotyledon development , whereas GRF7 is a transcriptional repressor of abscisic acid and osmotic stress-responsive genes . However, their functions in photomorphogenesis remain unknown.
To assay the regulatory roles of miR396–GRF pairs in photomorphogenesis, we first isolated and analyzed the mir396a single mutant (SALK_064047; Additional file 2: Figure S2a). Although overall miR396 levels were reduced, the phenotypes of the mir396a mutant were indistinguishable from that of the wild type under dark or light (Additional file 2: Fig. S2a and b). Possibly, the residual amount of miR396b in mir396a is sufficient for normal seedling development (Additional file 2: Figure S2a).
Because both miR396a and miR396b can suppress GRFs , we sought to simultaneously sequester these two miRNAs by generating target mimicry lines (MIM396) (Fig. 6a). The levels of GRF1, GRF2, GRF3 and GRF7 were indeed increased in two independent MIM396 lines (Fig. 6a). The MIM396 lines showed elongated hypocotyl length under 50 μE white light (Fig. 6b), which suggests that functional miR396 can positively regulate photomorphogenesis. We also examined the hypocotyl lengths of the MIR396A overexpression line (MIR396Aox)  and found that light sensitivity was not further exaggerated (Additional file 2: Figure S2c and d), so the endogenous miR396 pool may be at a saturated level for its functions in light responses.
To further understand the mechanistic roles of miR396 in photomorphogenesis, we examined the hypocotyl lengths of grf1 grf2 grf3 and grf7 mutants. Compared to wild-type Ws, the grf1 grf2 grf3 triple mutant showed short hypocotyl under 50 μE white light, which indicates that the three GRFs act as negative regulators of photomorphogenesis. In contrast, grf7 showed a short hypocotyl only under dark (Fig. 6c). Together with the quick repression of GRF7 expression by light (Fig. 5e), the major function of GRF7 is likely to promote hypocotyl elongation under dark. In sum, miR396 can act as a positive regulator of hypocotyl elongation by suppressing the negative regulator GRFs.
Varied length and target properties of siRNAs in Arabidopsis seedlings
Among the expressed siRNAs in de-etiolating seedlings, 164 are capable of mediating target cleavage (Additional file 1: Table S6). Among them, 70 (> 40%) are 24 nt long (Fig. 7a), which is quite different from the typical 21- or 22-nt miRNAs/phasiRNAs that mediate target mRNA cleavage [14, 19, 41]. Most target genes with degradome signatures for both miRNAs/phaiRNAs and siRNAs are protein coding genes (Fig. 7b). Intriguingly, 30 cleavage events from the actions of siRNAs were identified from 18 transposable elements (TEs), which is significantly higher than the number targeted by miRNAs/phasiRNAs (p = 3.2 × 10−3 by Fisher’s exact test). Most siRNA-targeted transposons are in the gypsy-like retrotransposon and CACTA-like transposase family. Eleven sRNAs that mediate TE mRNA cleavage are also derived from annotated TEs, with 4 sRNAs derived from their target loci (Additional file 1: Table S6), so these TEs may be capable of self-suppressing through TE-derived siRNAs and self-targeted cleavage. Among the 30 cleavage events derived from TE mRNAs, 23 potentially resulted from cleavage mediated by 24-nt siRNAs. The 24-nt siRNAs derived from transposable elements can mediate silencing of their original transposable elements via RNA-dependent DNA methylation (RdDM) [17, 18]. Our results suggest that in addition to RdDM, siRNA-mediated cleavage may function as an additional mechanism to prevent TE mRNA accumulation, which may escape from incomplete RdDM.
Light regulates miRNA-target pairs in photomorphogenic development
In addition to the previously discovered miR157d–HY5 and miR319–TCP pairs , we found additional miRNA–target pairs that positively regulate photomorphogenesis. In this study, miR168 could tune down AGO1 level under light (Fig. 3), which can counteract the actions of sRNAs stabilized by the light-induced HEN1 expression. miR396 can act as positive regulators of photomorphogenesis by suppressing GRF1, GRF2, GRF3 and GRF7. The grf7 mutant shows shorter hypocotyl length under dark, so it may positively regulate hypocotyl elongation under dark. Under light conditions, the grf7 mutant phenotype is comparable to that of the wild type possibly because light also markedly represses GRF7 expression (Fig. 5e, Additional file 1: Table S8).
We cannot rule out that light down-regulates the expression of GRFs at the transcriptional level. However, the detection of the miR396-mediated cleavage events on GRFs (Additional file 2: Table S5) suggested that miR396 indeed functions to optimize the GRF mRNA levels during de-etiolation. GRFs are known as transcription activators ; hence, future investigation of GRF downstream genes will help demystify genes regulated by the miR396-GRF lineage and provide a future research direction for their contribution in photomorphogenic development.
sRNAs regulate photomorphogenesis from multiple angles
Our results in Fig. 2 and Additional file 2: Figure S1 showed that abundant sRNAs have a better likelihood of mediating target mRNA cleavage during photomorphogenic development. Also, despite no negative correlation between the expression of miR168/miR396 and their target mRNAs (Figs. 3 and 5), degradome signatures from their target mRNAs were observed (Additional file 1: Table S5, Fig. 5). This finding indicated that although light does not affect the accumulation of miR168 and miR396, these miRNAs can contribute to the expression repression of their target mRNAs in de-etiolating seedlings. The steady-state mRNA transcriptome during photomorphogenic growth likely is a finely orchestrated balance between the well-studied transcriptional regulation by light signals and target mRNA cleavage mediated by small regulatory RNAs as examples shown in this study.
Combined with our previous  and current discovery, sRNAs could fine-tune the expression of both positive (HY5) and negative (TCPs, AGO1, GRFs) regulators of photomorphogenesis (Fig. 8). Clearly, as key regulators of these complex and interlocked regulatory circuits, the whole plethora of sRNAs is crucial for an optimal transcriptome during photomorphogenesis. This observation also explains why mutations of single MIR or target gene usually show less prominent phenotypic changes (Fig. 6), as compared with mutants with a defective miRNA pathway [23,24,25, 30, 43,44,45,46]. We have observed a considerable amount of degradome signatures that were predicted to be results of siRNA-mediated mRNA cleavage. This suggests that, in addition to miRNAs, siRNAs also contribute considerably to down regulate their target mRNAs in de-etiolating seedlings. Further investigation of the miRNA- and siRNA-target pairs will continue to shed light on post-transcriptional regulation of photomorphogenic growth.
Finally, our observation suggests that siRNA-mediated TE mRNA cleavage may serve as an additional mode of TE silencing (Fig. 7). In Arabidopsis, TE mRNAs could also be cleaved by miR859  and a tRNA-derived small RNA via the association with AGO1 in pollens . In Drosophila germ cells, Piwi-interacting RNAs (piRNAs) can interact with Aubergine (Aub) or AGO3 for the cleavage of TE mRNAs [49, 50]. It remains to be clarified with which AGO protein(s) siRNAs interact for silencing plant TE mRNAs.
Photomorphogenesis is a coordinated result of gene expression regulation at multiple levels. Our analyses revealed multiple sRNA–mRNA pairs contributing to this important development process. We also confirmed a comprehensive impact of sRNAs on regulating post-transcriptional gene expression during de-etiolation in Arabidopsis. sRNAs target multiple positive and negative regulators of photomorphogenesis, offering sophisticated fine-tuning power for regulating gene expression during de-etiolation. The potency of an sRNA in target cleavage is primarily determined by its abundance, adding an extra regulation dimension in addition to target recognition.
Plant materials and growth conditions
Seeds of wild-type Arabidopsis, Col-0, Ler, Ws, T-DNA insertion lines SALK_064047 (mir396a) and SAIL_1256_F08 (grf7) were acquired from stock centers, ABRC or NASC. Homozygous lines of T-DNA insertion lines were screened and confirmed for phenotype observation. The grf1 grf2 grf3 triple mutant and 35S:MIR396aox lines were kindly provided by Drs. Jeong Hoe Kim and DiQiu Yu, respectively. For phenotype observation, Arabidopsis seeds were surface-sterilized with 30% bleach and sown on half-strength Murashige and Skoog medium (Duchefa) without supplementing vitamin or sucrose, with 0.8% phyto agar at tissue culture grade (Duchefa, CAS number 9002–18-0). Seeds were stratified (4 °C for 4 days in the dark) to synchronize germination, then exposed to white light for 1 h to stimulate germination, and transferred to different light (W) conditions at 22 °C (Dark, W 12.5 and 50 μE for 4 days). The white light source was a PHILIPS LIFEMAX T-LD 18 W/840 T25 cool white tube. Hypocotyl lengths of seedlings were measured by using ImageJ v1.47 . The means and SEM were calculated from more than 30 seedlings. At least 3 biological replicates for each line were used for each experiment.
Construction of MIM396 lines
Primer sequences used in this study are listed in Additional file 1: Table S10. 35S:MIM396 target mimicry lines were generated as described . Briefly, the genomic fragment of IPS1 was amplified by using the iProof High-fidelity PCR kit (Bio-Rad) and cloned into the pGEMT-easy vector (Promega). The miR399 target site on original IPS1 sequence was modified to sequester miR396a/b, as shown in Fig. 6a, by overlapping PCR during construction . All constructs were then subcloned into the pCambia-1390 binary vector (CSIRO) digested with SalI and SacI. The constructs were transformed into Agrobacterium tumefaciens GV3101 strain, and introduced into Arabidopsis Col-0 by floral dipping. Two independent homozygous transgenic lines per construct were used for further analyses.
RNA sequencing and data analyses
For sRNA sequencing and data analyses, 4-d-old dark-grown Arabidopsis Col-0 seedlings were exposed to white light (100 μE) for 1, 3, 6, 12 and 24 h. The aerial tissues of approximately 5000 seedlings were collected for RNA isolation. Ten to 15 μg total RNA was size fractionated on 15% Tris-Borate-EDTA-Urea gel. sRNAs ranging from 17 to 30 nucleotides were gel-purified and used for cDNA library construction (Illumina Truseq for replicates 1 and 2, Small RNA v1.5 for replicates 3) and sequencing with the use of an Illumina HiSeq 2500. Twelve barcoded samples were sequenced in one single flowcell (a total of 240 M reads output per flowcell) at a read length of 50 nt. The adaptor-trimmed reads with size >18 nt were mapped to the Arabidopsis TAIR10 genome by using Bowtie  with the parameters -f –n 0 –e 80 –l 18 –a –m 5 –best –strata. For miRNA profiling, reads that perfectly matched to mature miRNA sequences were counted, normalized to total mapped reads of 20–24 nt and were shown as reads per million reads (RPM). Reads that mapped to miRNA families (e.g., miR156) were weighted by dividing the read count and equally assigning to each miRNA family member. For siRNA quantification, the Bowtie parameters were -f –n 0 –e 80 –l 18 –a –v 2 –best –strata. For phasiRNAs, the prediction of PHAS loci involved use of the UEA sRNA Workbench . In brief, adaptor-trimmed reads longer than 16 nt were mapped to the TAIR10 genome, and the 21-nt phasing register was set to detect phasing within a 251-nt window, based on hypergeometric distributions described previously . Among the uniquely mapped miRNA/siRNAs, only those with read counts ≥5 in ≥1 time point for all 3 biological replicates were considered expressed. Light-regulated sRNAs were defined as sRNAs with p < 0.05 on Student’s t-test compared to dark treatment (W0 h) for all 3 biological replicates. For mRNA transcriptome analysis, the RPKM for sRNA target genes were analyzed by using datasets published previously . Potential targets included those predicted by use of psRNATarget  (with UPE = 25, expectation = 3), miRNA targets identified in previous studies [25, 37] and miRNA targets detected in our degradome analysis (see Additional file 1: Table S5). Expressed genes had RPKM > 0.01 in at least one time point in both biological replicates. The transcript levels of target genes with degradome signatures are in Additional file 1: Table S5 and S6.
For degradome sequencing, 100 μg total RNA was isolated from 4-d-old dark-grown seedlings and mixtures of 4-d-old dark-grown seedlings exposed to 1, 3, 6, 12 and 24 h of light. Degradome sequencing was performed as described [56,57,58]. Putative cleavage sites were identified by using Cleaveland v4.4.3 [32, 59]. Those with CleaveLand category ≤ 2, p ≤ 0.05 and at least 5 reads at the predicted cleavage site were reported as valid targeting events.
Northern blot analysis and qRT-PCR
In total, 20 to 50 μg total RNA was separated on 15% TBE-Urea gel (Invitrogen). SYBR-Gold (Life Technologies) was used for visualizing RNAs on gels. RNAs were then transferred to Hybond-N+ Nylon membrane (GE Healthcare), by using Transblot SD Semi-Dry Transfer Cell (Bio-Rad) and hybridized with γ-32P-labeled miRNA probes as indicated at 37° Celsius overnight in UltraHyb Oligo buffer (Ambion). Hybridized blots were washed and exposed to Phosphoimager (GE Healthcare), then analyzed by using Typhoon FLA 7000 (GE Healthcare Life Sciences), as described . Images were quantified by using ImageJ v1.47 . For qRT-PCR, cDNA was synthesized from 2 μg total RNA from 4-d-old de-etiolating Arabidopsis seedlings. The SuperScript II RT kit (Invitrogen) was used for reverse transcription of mRNA. For qRT-PCR, cDNA with 0.25 ng equivalence of mRNA was used as a template for each sample. PCR amplification and detection was as described . Primers are in Additional file 1: Table S8. Data for one representative biological replicate were shown in Figs. 3c, 5e and 6a. Results for 2 additional biological replicates were shown in Additional file 2: Figure S3.
CONSTITUTIVE PHOTOMORPHOGENESIS 1
GROWTH REGULATING FACTOR
HYPONASTIC LEAVES 1
- K-S test:
- miRNA or miR:
- PHAS :
RNA-dependent DNA methylation
read per million reads
small interfering RNA
small regulatory RNA
- TAS :
trans-acting siRNA-generating loci
TEOSINTE BRANCHED 1
CYCLOIDEA AND PCF TRANSCRIPTION FACTOR
Sullivan JA, Deng XW. From seed to seed: the role of photoreceptors in Arabidopsis development. Dev Biol. 2003;260:289–97.
Kami C, Lorrain S, Hornitschek P, Fankhauser C. Light-regulated plant growth and development. Curr Top Dev Biol. 2010;91:29–66.
Wu SH. Gene expression regulation in photomorphogenesis from the perspective of the central dogma. Annu Rev Plant Biol. 2014;65:311–33.
Yeh KC, Wu SH, Murphy JT, Largarias JC. A cyanobacterial Phytochrome two-component light sensory system. Science. 1997;277:1505–8.
Ahmad M, Jarillo JA, Cashmore AR. Chimeric proteins between cry1 and cry2 Arabidopsis blue light photoreceptors indicate overlapping functions and varying protein stability. Plant Cell. 1998;10:197–207.
Rizzini L, Favory JJ, Cloix C, Faggionato D, O'Hara A, Kaiserli E, Baumeister R, Schäfer E, Nagy F, Jenkins GI, Ulm R. Perception of UV-B by the Arabidopsis UVR8 protein. Science. 2011;332:103–6.
Wang H, Wang H. Phytochrome signaling: time to tighten up the loose ends. Mol Plant. 2015;8:540–51.
Yin R, Arongaus AB, Binkert M, Ulm R. Two distinct domains of the UVR8 photoreceptor interact with COP1 to initiate UV-B signaling in Arabidopsis. Plant Cell. 2015;27:202–13.
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–607.
Chattopadhyay S, Ang LH, Puente P, Deng XW, Wei N. Arabidopsis bZIP protein HY5 directly interacts with light-responsive promoters in mediating light control of gene expression. Plant Cell. 1998;10:673–83.
Liu MJ, Wu SH, Chen HM, Wu SH. Widespread translational control contributes to the regulation of Arabidopsis photomorphogenesis. Mol Syst Biol. 2012;8:566.
Liu MJ, Wu SH, Wu JF, Lin WD, Wu YC, Tasi TY, Tsai HL, Wu SH. Translational landscape of Photomorphogenic Arabidopsis. Plant Cell. 2013;25:3699–710.
Meister G, Tuschl T. Mechanisms of gene silencing by double-stranded RNA. Nature. 2004;431:343–9.
Chen HM, Li YH, Wu SH. Bioinformatic prediction and experimental validation of a microRNA-directed tandem trans-acting siRNA cascade in Arabidopsis. Proc Natl Acad Sci U S A. 2007;104:3318–23.
Fei Q, Xia R, Meyers BC. Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks. Plant Cell. 2013;25:2400–15.
Allen E, Xie Z, Gustafson AM, Carrington JC. microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005;121:207–21.
Law JA, Jacobsen SE. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11:204–20.
Matzke MA, Mosher RA. RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat Rev Genet. 2014;15:394–408.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Chen X. Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009;35:21–44.
Guleria P, Mahajan M, Bhardwaj J, Yadav SK. Plant small RNAs: biogenesis, mode of action and their roles in abiotic stresses. Genomics Proteomics Bioinform. 2011;9:183–99.
Ruiz-Ferrer V, Voinnet O. Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol. 2009;60:485–510.
Bollman KM. HASTY, the Arabidopsis ortholog of exportin 5/MSN5, regulates phase change and morphogenesis. Development. 2003;130:1493–504.
Sorin C, Bussell JD, Camus I, Ljung K, Kowalczyk M, Geiss G, McKhann H, Garcion C, Vaucheret H, Sandberg G, Bellini C. Auxin and light control of adventitious rooting in Arabidopsis require ARGONAUTE1. Plant Cell. 2005;17:1343–59.
Tsai HL, Li YH, Hsieh WP, Lin MC, Ahn JH, Wu SH. HUA ENHANCER1 is involved in posttranscriptional regulation of positive and negative regulators in Arabidopsis photomorphogenesis. Plant Cell. 2014;26:2858–72.
Cho SK, Ben Chaabane S, Shah P, Poulsen CP, Yang SW. COP1 E3 ligase protects HYL1 to retain microRNA biogenesis. Nat Commun. 2014;5:5867.
Zhang H, He H, Wang X, Wang X, Yang X, Li L, Deng XW. Genome-wide mapping of the HY5-mediated gene networks in Arabidopsis that involve both transcriptional and post-transcriptional regulation. Plant J. 2011;65:346–58.
Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152–7.
Cheng Y, Chen X. Posttranscriptional control of plant development. Curr Opin Plant Biol. 2004;7:20–5.
Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert JL, Bartel DP, Crete P. Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell. 2004;16:69–79.
Baumberger N, Baulcombe D. Arabidopsis ARGONAUTE1 is an RNA slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci U S A. 2005;102:11928–33.
Brousse C, Liu Q, Beauclair L, Deremetz A, Axtell MJ, Bouche N. A non-canonical plant microRNA target site. Nucleic Acids Res. 2014;42:5270–9.
Arvey A, Larsson E, Sander C, Leslie CS, Marks DS. Target mRNA abundance dilutes microRNA and siRNA activity. Mol Syst Biol. 2010;6:363.
Ragan C, Zuker M, Ragan MA. Quantitative prediction of miRNA-mRNA interaction based on equilibrium concentrations. PLoS Comput Biol. 2011;7:e1001090.
Vaucheret H, Vazquez F, Crete P, Bartel DP. The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 2004;18:1187–97.
Chung PJ, Park BS, Wang HI, Jang IC, Chua NH. Light-inducible MiR163 targets PXMT1 transcripts to promote seed germination and primary root elongation in Arabidopsis. Plant Physiol. 2016;170:1772–82.
Jones-Rhoades MW, Bartel DP. Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004;14:787–99.
Liu D, Song Y, Chen Z, Yu D. Ectopic expression of miR396 suppresses GRF target gene expression and alters leaf growth in Arabidopsis. Physiol Plant. 2009;136:223–36.
Kim JH, Choi D, Kende H. The AtGRF family of putative transcription factors is involved in leaf and cotyledon growth in Arabidopsis. Plant J. 2003;36:94–104.
Kim JS, Mizoi J, Kidokoro S, Maruyama K, Nakajima J, Nakashima K, Mitsuda N, Takiguchi Y, Ohme-Takagi M, Kondou Y, et al. Arabidopsis growth-regulating factor7 functions as a transcriptional repressor of abscisic acid- and osmotic stress-responsive genes, including DREB2A. Plant Cell. 2012;24:3393–405.
Chen HM, Chen LT, Patel K, Li YH, Baulcombe DC, Wu SH. 22-nucleotide RNAs trigger secondary siRNA biogenesis in plants. Proc Natl Acad Sci U S A. 2010;107:15269–74.
Lee BH, Ko JH, Lee S, Lee Y, Pak JH, Kim JH. The Arabidopsis GRF-INTERACTING FACTOR gene family performs an overlapping function in determining organ size as well as multiple developmental properties. Plant Physiol. 2009;151:655–68.
Li J, Yang Z, Yu B, Liu J, Chen X. Methylation protects miRNAs and siRNAs from a 3'-end uridylation activity in Arabidopsis. Curr Biol. 2005;15:1501–7.
Park MY, Wu G, Gonzalez-Sulser A, Vaucheret H, Poethig RS. Nuclear processing and export of microRNAs in Arabidopsis. Proc Natl Acad Sci U S A. 2005;102:3691–6.
Vaucheret H, Mallory AC, Bartel DP. AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell. 2006;22:129–36.
Yang L, Liu Z, Lu F, Dong A, Huang H. SERRATE is a novel nuclear regulator in primary microRNA processing in Arabidopsis. Plant J. 2006;47:841–50.
Creasey KM, Zhai J, Borges F, Van EF, Regulski M, Meyers BC, Martienssen RA. miRNAs trigger widespread epigenetically activated siRNAs from transposons in Arabidopsis. Nature. 2014;508:411–5.
Martinez G, Choudury SG, Slotkin RK. tRNA-derived small RNAs target transposable element transcripts. Nucleic Acids Res. 2017;45:5142–52.
Aravin AA, Sachidanandam R, Girard A, Fejes-Toth K, Hannon GJ. Developmentally regulated piRNA clusters implicate MILI in transposon control. Science. 2007;316
Sato K. The epigenetic regulation of transposable elements by PIWI-interacting RNAs in drosophila. Genes Genet Syst. 2013;88:9–17.
Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.
Franco-Zorrilla JM, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, Leyva A, Weigel D, Garcia JA, Paz-Ares J. Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet. 2007;39:1033–7.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Stocks MB, Moxon S, Mapleson D, Woolfenden HC, Mohorianu I, Folkes L, Schwach F, Dalmay T, Moulton V. The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics 2012;28:2059–2061.
Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39:W155–9.
Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ. Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008;18:758–62.
German MA, Pillay M, Jeong DH, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, et al. Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008;26:941–6.
Gregory BD, O'Malley RC, Lister R, Urich MA, Tonti-Filippini J, Chen H, Millar AH, Ecker JR. A link between RNA metabolism and silencing affecting Arabidopsis development. Dev Cell. 2008;14:854–66.
Addo-Quaye C, Miller W, Axtell MJ. CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics. 2009;25:130–1.
Wang Y, Wu JF, Nakamichi N, Sakakibara H, Nam HG, Wu SH. LIGHT-REGULATED WD1 and PSEUDO-RESPONSE REGULATOR9 form a positive feedback regulatory loop in the Arabidopsis circadian clock. Plant Cell. 2011;23:486–98.
We specially thank Ming-Tsung Wu and Shu-Jen Chou for establishing the degradome sequencing system, Bo-Jun Liu and Chi-Yu Chen for assisting in experiments related to MIM396 under the supervision of MCL and SHW. We also thank Jeong Hoe Kim at Kyungpook National University and Diqiu Yu at Xishuangbanna Tropical Botanical Garden, Chinese Academiy of Sciences, for kindly providing the grf mutants and 35S:MIR396aox line, respectively. Finally, we thank Chih-Ching Lin, Chiung-swey Joanne Chang, Ho-Ming Chen, Ho-Wei Wu, Wen-Ping Hsieh, Yi-Chen Wu and Ying Wang for inspiring discussion and suggestions.
This research, including design and collection of experimental results, data analyses and preparation of the manuscript, was supported by the Investigator Award to S.-H. Wu from Academia Sinica, Taiwan.
Availability of data and materials
The datasets generated and analyzed in this manuscript can be accessed from Gene Expression Omnibus (GEO) under accession number GSE83646. The other supporting data are included as Additional files.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Table S1. Sequencing and mapping statistics. Description: Contains the number and percentage of mapped sRNA reads in sRNA and degradome sequencing. Table S2. Expressed miRNAs in de-etiolating Arabidopsis seedlings. Description: Contains the expression level of expressed miRNAs, with RPM and fold changes. Table S3. Expressed phasiRNAs in de-etiolating seedlings. Description: Contains the expression level of expressed phasiRNAs, with the corresponding PHAS locus, RPM and fold changes. Table S4. Expressed siRNAs in de-etiolating seedlings. Description: Contains the expression level of expressed siRNAs, with the sequence, length, RPM and fold changes. Table S5. miRNA-mRNA and phasiRNA-mRNA pairs from cross comparisons of sRNA transcriptome and mRNA degradome. Contains the information of miRNA/phasiRNA-mediated target cleavage, including locus number, CleaveLand categories, p-values and target mRNA levels. Table S6. siRNA-mRNA pairs from cross comparisons of sRNA transcriptome and mRNA degradome. Contains the information of siRNA-mediated target cleavage, including locus number, CleaveLand categories, p-values and target mRNA levels. Table S7. Expression levels of GRFs. Contains the expression level of GRFs under dark and light. Table S8. Primers used in this study. Contains sequences of the primers that are used in this study. (XLSX 2667 kb)
Figure S1. miRNAs/siRNAs with target cleavage tend to be more abundant than miRNAs/siRNAs that did not show target cleavage signatures. Contains supplemental figure and legend showing K-S test results of siRNA abundance. Figure S2. Molecular and phenotypic analyses of mir396a mutant and MIR396aox lines. Contains supplemental figure and legend showing the examination of mir396a mutant and MIR396aox lines. Figure S3. qRT-PCR results of two additional biological replicates for a AGO1 and b, c GRF shown in Figs. 3c, 5e and 6a, respectively. Contains supplemental figure and legend showing the qRT-PCR results of two additional biological replicates in this study. (PDF 1074 kb)