Unraveling multifaceted contributions of small regulatory RNAs to photomorphogenic development in Arabidopsis

Background 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. Results 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. Conclusion 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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3937-6) contains supplementary material, which is available to authorized users.


Background
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 [9]. 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 lightresponsive genes during de-etiolation [10]. 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].
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 [26]. 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 [25]. 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 [27]. 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.

Results
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 [28], 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) [14] were expressed in de-etiolating seedlings (Table 1 and Additional file 1: Table S3). In addition to miRNAs and pha-siRNAs, 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 [32] 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 Fig. 1 Pipeline for investigating sRNA-mediated regulation in de-etiolating Arabidopsis seedlings. Four-d-old dark-grown seedlings (W0 h) were exposed to continuous white light (W; 100 μE) for 1, 3, 6, 12 and 24 h; 18-22 million (M) reads per library were acquired. Degradome sequencing for dark-(Dark) and light-(equally pooled from the five light-treated time points) grown seedlings were used for identifying sRNA-mediated cleavage of target mRNA. Colored numbers in parentheses indicate the number of sRNAs/targets passing the corresponding filtering criteria  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 [33]. 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 [34]. 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 [19]. The availability of transcriptome data for mRNAs [12], 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. 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 [25]. 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 [35]. 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 a c b 0 0 inhibition of hypocotyl elongation [36]. Cluster II miRNAs are also upregulated by light, and the up regulation is more prominent after prolonged light exposure. The miR-NAs in cluster II include miR157d, reported to target HY5 during photomorphogenic development [25]. The miR-NAs 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 [25], 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 [39], whereas GRF7 is a transcriptional repressor of abscisic acid and osmotic stress-responsive genes [40]. 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 [38], 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) [38] 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   [17,18]. Our results suggest that in addition to RdDM, siRNAmediated 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 [25], 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 [42]; 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 [25] 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-25, 30, 43-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 [47] and a tRNA-derived small RNA via the association with AGO1 in pollens [48]. 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. a Expressed miRNA families with predicted (black) or identified (red) target cleavages in de-etiolating Arabidopsis. b Expression of miR396 is transiently upregulated and gradually decreased on W exposure. Data are mean ± SD from three biological replicates of sRNA sequencing. c Northern blot analyses confirmed the expression of miR396 in de-etiolating Arabidopsis. One representative gel image is shown. Data are mean ± SD for the relative expression from three biological replicates. SYBR-Gold-stained rRNA/tRNA was a loading control. d Degradome T-plot marked the miR396-mediated GRF1/ GRF2/GRF3/GRF7 mRNA cleavage in de-etiolating seedlings. Red asterisks indicate the degradome signatures detected at expected cleavage sites for GRF1, GRF2, GRF3 and GRF7. e Light regulation of GRF1, GRF2, GRF3 and GRF7. Data are mean ± SD from three technical replicates of one representative qRT-PCR experiment. * p < 0.01 in Student's t-test. Three biological replicates were performed with similar results

Conclusions
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 posttranscriptional 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  [51]. 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 [52]. 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 [52]. 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 darkgrown 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 [53] with the parameters -f -n 0 -e 80 -l 18 -a -m 5 -beststrata. 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 [54]. 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 [14]. 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 [12]. Potential targets included those predicted by use of psRNATarget [55] (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 γ-32 P-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 [14]. Images were quantified by using ImageJ v1.47 [51]. For qRT-PCR, cDNA was synthesized from 2 μg total RNA from 4-d-old deetiolating 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 [60]. 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.

Additional files
Additional file 1: 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) Additional file 2: 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.