Skip to main content

Plastid retrograde regulation of miRNA expression in response to light stress



MicroRNAs (miRNAs) are a class of endogenous noncoding RNAs that play a pivotal role in the regulation of plant development and responses to the surrounding environment. Despite the efforts made to elucidate their function in the adaptation of plants to many abiotic and biotic stresses, their role in high light (HL) stress is still vague. HL stress often arises upon plant exposure to full sunlight. Subsequent changes in nuclear gene expression are triggered by chloroplast-derived retrograde signals.


In this study, we show that HL is involved in miRNA-dependent regulation in Arabidopsis thaliana rosettes. Microtranscriptomic screening revealed a limited number of miRNAs reacting to HL. To explain the miRNA regulation mechanisms at the different biogenesis stages, chemical and genetic approaches were applied. First, we tested the possible role of plastoquinone (PQ) redox changes using photosynthetic electron transport chain inhibitors. The results suggest that increased primary transcript abundance (pri-miRNAs) of HL-regulated miRNAs is dependent on signals upstream of PQ. This indicates that such signals may originate from photosystem II, which is the main singlet oxygen (1O2) source. Nevertheless, no changes in pri-miRNA expression upon a dark–light shift in the conditional fluorescent (flu) mutant producing 1O2 were observed when compared to wild-type plants. Thus, we explored the 1O2 signaling pathway, which is initiated independently in HL and is related to β-carotene oxidation and production of volatile derivatives, such as β-cyclocitral (β-CC). Pri-miRNA induction by β-CC, which is a component of this 1O2 pathway, as well as an altered response in the methylene blue sensitivity 1 (mbs1) mutant support the role of 1O2 signaling in miRNA regulation.


We show that light stress triggers changes in miRNA expression. This stress response may be regulated by reactive oxygen species (ROS)-related signaling. In conclusion, our results link ROS action to miRNA biogenesis, suggesting its contribution to inconsistent pri- and mature miRNA dynamics.

Peer Review reports


Plants take fundamental advantage of light absorption, but constant light fluctuations often result in episodes of excess light energy [1]. Many physiological and molecular processes are engaged in adjusting the plant response to light stress because of its frequent occurrence and possible consequences in photoinhibition [2]. Light-originated stimuli are perceived by chloroplasts and trigger retrograde signaling, resulting in nuclear gene expression changes [3,4,5]. MicroRNAs (miRNAs) are nucleus-encoded molecules that play an effective fine-tuning role in the plant response to environmental stresses [6]. Thus, the study of their contribution to the light stress response involving chloroplast-derived signals is extremely important. Moreover, miRNAs largely target transcription factors (TFs), which have paramount importance for plant growth, reproduction, and defense [7,8,9]. Surprisingly, knowledge concerning light-regulated miRNAs is limited. For example, far-red light-responsive miRNAs were described in soybean [10]. Nevertheless, there is almost no information about miRNA expression changes under high light (HL) conditions except studies on Dendrocalamus latiflorus, an important Asian bamboo species [11], a recent report describing small RNAs during high light acclimation [12] and the systemic miRNA response in Arabidopsis thaliana roots [13]. Additionally, there are some reports on UV-A-, UV-B-, and gamma radiation-regulated miRNAs [14,15,16,17,18].

Almost all MIR genes (encoding miRNAs) are transcribed as independent transcriptional units by RNA polymerase II (PolII) and may be regulated through cis regulatory promoter elements [9, 19,20,21]. Mature miRNAs are processed from longer primary transcripts (pri-miRNAs) and cleaved by a core microprocessor complex consisting of type III RNAse, DICER-LIKE1 (DCL1), zinc finger protein SERRATE (SE), and dsRNA binding protein HYPONASTIC LEAVES1 (HYL1) [22,23,24]. Stepwise cleavage generates precursor miRNAs (pre-miRNAs) at the first step and mature miRNAs in the second cleavage reaction step. miRNAs are then transported to the cytoplasm and loaded onto ARGONAUTE proteins to form the sequence-specific RNA-induced silencing complex (RISC) [25]. Once programmed with an miRNA, RISC can silence target genes by translational inhibition, mRNA cleavage, or heterochromatin formation [26]. As trans-acting regulators, foremost efforts were put into elucidating the temporal and spatial expression changes of individual miRNAs and targeted genes, while the exact regulatory mechanism is often obscure. Therefore, transcriptional and posttranscriptional control of miRNA expression is the subject of intensive research, yet there is a knowledge gap with respect to its crosstalk with chloroplast-derived signals [27,28,29].

Light propels photosynthesis, which is a well-established source of reactive oxygen species (ROS) in plants [30]. Perturbations in photosynthesis during light stress conditions lead to intensified production of ROS—1O2 (singlet oxygen) at PSII and H2O2 and O2•− at PSI. To date, only a few 1O2-dependent retrograde signaling pathways have been reported. In the grana margins (GM), 1O2 is mostly produced from tetrapyrrole biosynthesis intermediates and activates three signaling pathways: 1) EX1/EX2-dependent programmed cell death, 2) E3 ubiquitin ligase plant U-box 4 (PUB4)-dependent selective chloroplast degradation, and 3) a newly described SAFEGUARD1 (SAFE1)-dependent pathway [31,32,33]. Nevertheless, 1O2 is also produced in the grana core (GC), and the signal is transduced through the 1) β-cyclocitral (β-CC)-dependent or 2) OXI1 kinase-mediated signaling pathways [34, 35].

In addition to ROS production, electron flux between PSII and PSI also results in the redox regulation of photosynthetic electron transport chain (PET) components. Accordingly, HL causes plastoquinone (PQ) reduction, located downstream of PSII, to plastoquinol (PQH2). These PQ/PQH2 redox state changes are responsible for regulating at least 750 nuclear genes [3, 36]. In summary, the PQ pool redox state and 1O2-related signals may contribute to retrograde communication because they originate in chloroplasts and induce nuclear gene expression changes during HL conditions. There is some evidence for a PQ pool oxidation/reduction role in miRNA processing [37], as well as an influence on alternative splicing [38]. Nevertheless, the coupling of 1O2 signaling with miRNA abundance is currently an unexplored part of miRNA and retrograde signaling cross-talk.

Here, we demonstrate that HL causes changes in the miRNA levels in A. thaliana rosettes. The possible role of retrograde signals, especially 1O2, in miRNA expression was studied to verify the putative link to specific 1O2 signaling pathways. To achieve that, analyses of an A. thaliana flu mutant, conditionally producing 1O2, and mutants impaired in 1O2-mediated retrograde communications exposed to HL were conducted. β-CC application, a volatile retrograde signaling mediator, followed by monitoring of pri-miRNA expression, was also implemented. Taken together, our study is the first report linking ROS production and miRNA expression in the context of light stress. Nevertheless, future work should be undertaken to fully understand this mechanism and its biological function.


High light induces miRNA expression changes

To assess the impact of HL on miRNA expression changes, low light (LL)-acclimated Arabidopsis thaliana plants grown in hydroponic conditions were subjected to HL treatment for 2 h (HL), which was 10 times greater than that under the LL growth conditions (see section Methods). Such HL treatment was shown to cause a photoinhibitory effect and promote chloroplast to nucleus signaling [39,40,41]. The photoinhibitory effect was confirmed by measuring PSII activity, which was determined as the maximum quantum efficiency of PSII expressed as the Fv/Fm parameter, the ratio of variable to maximum chlorophyll fluorescence (Fig. 1a). Additionally, the induction of APX2, ELIP1, and RRTF1, known as HL response marker genes [1, 42, 43], confirmed the stress effect of the HL treatment (Additional file 1: Fig. S1). Since we hypothesized that miRNAs are involved in light-triggered gene regulation, we first searched the list of Arabidopsis nuclear genes coding for chloroplast proteins using available small RNA target prediction software and tested several miRNA candidates [44, 45] (Fig. 1b; Additional file 1: Table. S1). The expression of miR163 and miR5021 was upregulated in plants treated with HL (HL) versus that in control plants (LLc) with a simultaneous decrease in miR395 expression (Fig. 1b). To identify other miRNAs connected with the HL response, we performed microtranscriptomic sequencing of LLc-, HL- and LLr (LLr—HL-treated plants followed by 4 h recovery in LL)-treated plants, which allowed us to identify 21 miRNA candidates regulated by HL. Of these, 7 were up- and 14 were downregulated, with confirmed upregulation of miR163 under HL (Fig. 1c; Additional file 1: Fig. S2; Fig. S3). The observed HL-triggered miRNA expression changes were limited, dynamic, and rather subtle (fold changes ranging from approximately 0.4 to 2.8). Selected candidates were validated by the two-tailed RT-qPCR method (TT-RT qPCR), [46] and two of them, miR163 and miR840, which were upregulated just after HL stress (in the microtranscriptomic screen), were chosen for further analysis (Additional file 1: Fig.S4; Fig S5; Fig. 1c, d). HL induction was also visible 4 h after HL stimuli (LLr), although only in the case of miR840 it was confirmed by both RNAseq and qPCR analysis (Additional file 1: Fig. S2; Fig. S4; Fig. 1d). The involvement of miR163 and miR840 in the HL response is not surprising because miR163 was previously found to be induced by light during seedling de-etiolation [47, 48] after 6 h of HL [12] or red light treatment [49], and miR840 was described as a gamma-ray-responsive miRNA [16].

Fig. 1
figure 1

HL induces miRNA expression changes in Arabidopsis thaliana shoots. a Fv/Fm measured in 4-week-old Arabidopsis plants grown in LL under hydroponic conditions. LLc – control plants; HL—plants exposed to HL for 2 h (n = 18). b TT-qRT PCR for selected miRNA-targeted transcripts of nuclear encoded proteins localized in chloroplasts under LLc and HL. Transcript levels were normalized with respect to small nucleolar85 RNA (sno85) and small nucleolar101 RNA (sno101). c Results of a microtranscriptomic screening for miRNA expression changes in LLc, HL and LLr – plants exposed to HL for 2 h and a subsequent recovery for 4 h in LL. d TT-qRT PCR for miR163 and miR840 and e qRT-PCR for pri-miR163 and pri-miR840 in Col-0 LLc, HL and LLr plants. Transcript levels were normalized with respect to sno85 and sno101 (for miRNAs) or the PP2A and UPL7 genes (for pri-miRNAs), respectively. f Occupancy of PolII on MIR genes. Line charts present ChIP profile of total PolII on examined genes. Grey lines represent results for LLc plants and blue lines represent results for HL plants. Above each chart, gene structure is shown with black boxes representing miR, and grey boxes representing primary transcripts (pri-miRs).Orange lines show amplified regions (primer localization used for qRT-PCR analysis). Above each gene structure 100 bp scale is shown. g RNA stability assay performed on Arabidopsis seedlings in control plants (LLc) and plant exposed to HL (HL). Degradation curves after cordycepin treatment (plot) were used to calculate half-life of pri-miR840, pri-miR163, UBC (control) transcripts and short-lived mRNA transcribed from gene At3G45970 (table h). Presented values are averages from three biological replicates. For better clarity of chart only pri-miRNAs data were presented. Asterisks indicate significant differences according to Tukey’s HSD test (panel: a,b,d,e) or t-test (panel f) at the level of ** ≤ 0.01 and *** ≤ 0.001. Mean values ± SDs (n = 3) were provided

The maturation of functional miRNAs is a multistep process; consequently, a better understanding of their regulation cannot be limited to the analysis of mature forms of miRNAs. Therefore, we also monitored the expression of pri-miR163 and pri-miR840 (Fig. 1e). Pri-miR163 was highly accumulated under HL, while the expression of pri-miR840 was slightly increased under both HL and LLr. Although the expression of pri-miR163 is elevated 5 times under HL, changes in miR163 level do not exceed twofold (Fig. 1d) On the other hand, miR840 exhibits comparable level of expression fold changes at the analyzed stages of miRNA biogenesis. This divergence between miR163 and miR840 may be caused by different maturation process of analyzed miRNAs. It was proven that physical interactions between the DCL1 and HYL1 proteins are necessary for precise miR163 precursor processing [23]. Generally, most miRNAs require HYL1 for their processing; therefore, they are defined as HYL1-dependent miRNAs, and their precursors overaccumulate in hyl1 mutant plants. Nevertheless, HYL1 activity is not crucial for the maturation of all miRNAs [9, 23, 50, 51]. Moreover, several miRNAs may become HYL1/SE conditionally independent at decreased temperatures [52]. Therefore we investigated pri-miRNA levels in the Col-0 and hyl1 mutant (Additional file 1: Fig. S6). Pri-miR163 overaccumulated in the hyl1 mutant, while the pri-miR840 level was similar in Col-0 plants, which clearly indicated that pri-miR163 maturation is HYL1-dependent (Additional file 1: Fig. S6). Remarkably, the mature miR840 level in hyl1 was several-fold greater than that in Col-0 (Additional file 1: Fig. S6). Since the level of pri-miRNAs is an outcome of the transcription rate and stability of transcript we performed PolII-ChIP assay (Fig. 1f). This experiment show that no statistical differences in PolII occupancy was detected at tested miRNA genes between LL and HL conditions. On the other hand, GUS staining revealed that HL enhanced pri-miR163 expression (Additional file 1: Fig. S7). Probably, observed changes are connected with regulation of pri-miRNA stability rather than the transcription rate or PolII-ChIP experimental setup missed shorter gene activation periods. To test first hypothesis cordycepin assay was carried out. Stability of pri-miR163 increases after HL treatment in opposite to pri-miR840 whose half-life is substantially reduced (Fig. 1g-h).

In summary, miRNA sequencing followed by RT-qPCR analysis confirmed the miRNA response to HL. Moreover, this response was regulated at the different stages of miRNA biogenesis.

The role of PQ redox status in miRNA expression is not conclusive

To test the engagement of chloroplast-derived signals, we used inhibitors of PET because its components act as excess light messengers and nuclear gene expression regulators [1, 53,54,55]. The PQ redox status can be easily modulated using PET inhibitors. DCMU (3-3,4-dichlorophenyl-1,1-dimethylurea) blocks the PSII PQ binding site, oxidizing the PQ pool [56]. DBMIB (2,5-dibromo-6-isopropyl-3-methyl-1,4-benzoquinone) has been used as a specific inhibitor of plastoquinol oxidation at the Qo binding site of cytochrome b6f, causing reduction of the PQ pool [57] (Fig. 2a). Plants were kept in the dark (dark) or LL (LLc), and plants were treated with DCMU or DBMIB for 4 h in LL (LLtrt) (Fig. 2b,c; Additional file 1: Fig. S8). The effective inhibitor concentrations were determined by monitoring the Fv/Fm parameter, which decreased 4 h after chemical application (Fig. 2b). Since we showed that miRNAs can be regulated at different biogenesis stages, pri-miRNAs and mature miRNAs were monitored. Pri-miR163 and pri-miR840 were upregulated in light-treated plants compared to dark-treated plants (Fig. 2c). DCMU further enhanced pri-miR163 induction, whereas light-dependent pri-miR840 induction was almost completely abolished. In the DBMIB experiment, pri-miR163 was strongly upregulated, while the pri-miR840 level was downregulated compared to that in the LL control plants but remained significantly elevated compared to that in dark-treated plants. Since the dominant effect was related to the dark–light switch, we deduce that the increase in pri-miRNAs may not be dependent on PQ because it occurred in the DCMU treatment (PQ oxidized) as well as DBMIB treatment (PQ reduced). However, in the case of pri-miR840, DCMU seemed to block the effect of light to some extent (similar results were also observed in pri-miR319b, which was upregulated in the LLr treatment in the microtranscriptomic screening (Additional file 1: Fig. S2; S9). Since the analysis with inhibitors was not conclusive, we decided to check the expression level of pri-miR163 and pri-miR840 in protein kinase STATE TRANSITION 7 (stn7) and SALICYLIC ACID INDUCTION DEFICIENT 2 protein (sid2) mutant plants, which were previously shown to have a reduced PQ level compared to that of Col-0 plants (Additional file 1: Fig. S10) [58]. We observed no differences in pri-miRNA levels (and miRNA) in the tested mutants, which suggests that PQ redox status is not responsible for changes in pri-miRNAs.

Fig. 2
figure 2

The role of PQ redox status in miRNA expression. a Scheme representing the action of DCMU and DBMIB in the photosynthetic electron transport chain. b Treatment with DCMU and DBMIB influences the maximum efficiency of PSII expressed as the Fv/Fm parameter. Plant material: dark—plants kept in darkness for 4 h; LLc—control plants in LL; LLtrt—plants treated with DCMU or DBMIB for 4 h and kept in LL. Asterisks indicate significant differences according to Tukey’s HSD test at the level of *** ≤ 0.001. Mean values ± SDs (n = 9) were provided. c qRT-PCR for pri-miR163 and pri-miR840 after using DCMU or DBMIB (upper panel). Transcript levels were normalized with respect to the PP2A and UPL7 genes. TT-qRT PCR for miR163 and miR840 after treatment with DCMU or DBMIB (bottom panel). Transcript levels were normalized with respect to sno85 and sno101. Asterisks indicate significant differences according to Tukey’s HSD test at the level of * ≤ 0.05, ** ≤ 0.01 and *** ≤ 0.001. Mean values ± SDs (n = 3) were provided

Next, we tested the abundance of mature miRNAs after treatment with DCMU and DBMIB. Surprisingly, changes exhibited in pri-miRNAs were largely abolished in the mature forms (Fig. 2c). Although we still observed statistically significant differences in the miR163 expression level, the magnitude of the presented changes was scaled down (Fig. 2c, DBMIB panel) or abolished completely (Fig. 2c, DCMU panel). Because pri-miR163 induction in the DBMIB treatment was greater than that in DCMU, it can be assumed that the miRNA level response was almost equally reduced in both treatments. Simultaneously, the ratio and pattern of miR840 changes were maintained from the pri- to mature miRNAs, indicating different miR163 and miR840 maturation processes.

In summary, the presented results are not conclusive in defining the role of the PQ pool redox state in the light-triggered regulation of miRNAs. Thus, the impact of other possible retrograde signals should be considered.

Regulation of pri-miRNA expression is not dependent on EX1-dependent 1O2 signaling

Pri-miRNA expression changes induced by light were similar for both the DCMU and DBMIB treatments, which suggests that transcription of miRNAs may be induced upstream from PQ. The PSII light-harvesting antenna complex is the main place where the production of highly reactive 1O2, a photosynthesis byproduct, occurs [59,60,61,62,63]. Under mild stress conditions, 1O2 may promote programmed cell death by activating two nucleus-encoded proteins, EXECUTER1 (EX1) and EX2, which are located in the chloroplast GM, where chlorophyll is synthesized and the PSII repair cycle takes place [31, 64,65,66,67]. Under light stress, various ROS can be generated simultaneously; therefore, it is impossible to analyze the specific biological activity of 1O2. Fortunately, a flu conditional mutant that selectively overproduces 1O2 from the photosensitizer protochlorophyllide (Pchlide) initiated numerous studies on the signaling functions of this ROS [31, 68, 69]. The flu mutant in continuous light displays the wild-type phenotype because in these conditions, Pchlide is immediately photoreduced to chlide and consequently does not reach the level necessary for elevated 1O2 production [70]. By transferring light-grown flu plants to darkness for a period and then re-exposing them to light, we can easily modulate the Pchlide level because it accumulates proportionally over time and consequently produces 1O2. Interestingly, these symptoms are abrogated in flu/ex1 double mutants, demonstrating that Pchlide accumulation is not enough to trigger 1O2 signaling and requires EX1 protein function [67, 70, 71].

The potential role of 1O2 in the induction of MIR expression was verified by analysis of pri-miR163 and pri-miR840 expression changes in wild-type Col-0 and flu, ex1, and flu/ex1 double mutant seedlings grown in LL for 2 weeks and transferred to darkness for 12 h, then subsequently returned to LL for 2 h for 1O2 generation. After that time, Fv/Fm was measured in treated (trt) and control (ctrl, kept continuously in LL) plants (Fig. 3a). Treated flu seedlings exhibited a stress response (decrease in the Fv/Fm value) because of 1O2 overproduction in PSII. This decline in photosynthetic parameters was not exhibited in Col-0, ex1, and flu/ex1 plants, which implies that 1O2-mediated and EX1-dependent signaling takes place under noninhibitory light and that 1O2 produced in flu background seedlings does not directly damage PSII [67]. Moreover, 1O2 release was shown by elevated DRP expression, a known 1O2 marker gene [72, 73], in flu background plants (Fig. 3b). Pri-miRNA163 and pri-miR840 were upregulated in flu (pri-miR840 changes were statistically significant); however, this elevated level was also observed in Col-0, suggesting that it is not connected with EX1-dependent 1O2 signaling. Although Pchlide is first synthesized in the GM, when the darkness time exceeds 8 h (in the case of the flu and flu/ex1 plants), after reillumination, Pchlide also accumulates in GC and slightly accumulates in stroma lamellae [32]. Thus, the 12 h of darkness used in our experiments may also activate other EX1-independent 1O2 signaling pathways, which would explain the greater pri-miR840 level. No changes in pri-miRNAs in flu/ex1 plants indicate that regulation of these miRNAs is not dependent on EX1 (Fig. 3c). To confirm these observations, we exposed the ex1 mutant to HL stress for 2 h, and subsequent analysis of the pri-miRNA expression levels was carried out (Fig. 3d). Pri-miR163 expression was induced by HL in both Col-0 and ex1 plants. Although the pri-miR163 expression level in HL-treated mutant plants was reduced, it was already reduced in ex1 control plants, and the fold changes for ex1 were even greater than those in Col-0 (8.6 in ex1 versus 6.8 in wild-type plants) (Fig. 3d; Additional file 1: Fig. S11). We observed a slight increase in pri-miR840 expression in HL-treated Col-0 plants, while expression changes in LLr plants were nearly identical for Col-0 and ex1. Thus, our results suggest that the EX1-1O2 signaling pathway is not engaged in HL-induced miRNA regulation.

Fig. 3
figure 3

Regulation of pri-miRNA expression is not dependent on EX1-dependent 1O2 signaling. a The Fv/Fm parameter measured in control (ctrl) plants (grown in constant light (CL)) and treated (trt) plants (plants grown for 2 weeks in CL, placed for 12 h in darkness, and re-exposed for 2 h to LL) of the Col-0, ex1, flu, and flu/ex1 genotypes (n = 9–10). b qRT-PCR of the DRP gene in the Col-0, ex1, flu, and flu/ex1 genotypes in ctrl and trt plants. c qRT-PCR for pri-miR163 and pri-miR840 in the Col-0, ex1, flu, and flu/ex1 genotypes in ctrl and trt plants. d qRT-PCR for pri-miR163 and pri-miR840 in Col-0 and ex1 plants LLc-control plants; HL-plants exposed to HL for 2 h; LLr-plants exposed to HL for 2 h and subsequent recovery in LL for 4 h. Transcript levels were normalized with respect to the PP2A and UPL7 genes. Asterisks indicate significant differences according to Tukey’s HSD test at the level of * ≤ 0.05, ** ≤ 0.01 and *** ≤ 0.001. In (b) and (c), * indicates significance within the same genotype, while ^ indicates comparison to Col-0 within the same conditions. Mean values ± SDs (n = 3) were provided

pri-miRNA expression is regulated by β-CC-dependent 1O2 signaling

1O2 formation may also activate EX1-independent signaling, which occurs in photoinhibitory conditions. In such circumstances, light stress leading to β-carotene oxidative breakdown results in the release of small volatile compounds such as β-CC, which is known to induce 1O2-responsive genes (Singlet Oxygen Responsive Genes—SORGs, Fig. 4a) [34, 64]. Therefore, we tested the possible role of the β-CC-1O2-dependent pathway in miRNA expression. β-CC treatment (for details see Methods section) was utilized, and then we examined the DRP expression level to confirm that 1O2 signaling was activated (Fig. 4a-b). Significantly increased DRP expression was detected after β-CC (1 ml) application (Fig. 4b). Subsequently, we analyzed pri-miR163 and pri-miR840 expression levels in the same experimental setup. Pri-miR163 abundance was elevated at both β-CC doses (i.e., 50 µl and 1 ml), whereas pri-miR840 abundance was elevated at the greater dose only (Fig. 4c). Surprisingly, the level of both miRNAs decreased drastically after application of 1 ml of this chemical (Fig. 4d). Next, we considered the possible role of this retrograde signaling pathway in miRNA regulation. The β-CC signaling pathway requires the MBS1 protein, which is positioned downstream of β-CC, to mediate 1O2 signal transfer to the nucleus [68]. MBS1 was identified in a Chlamydomonas reinhardtii (green alga) mutant genetic screen displaying a defect in response to 1O2 generated by the photosensitizer methylene blue [74]. Moreover, an A. thaliana mbs1 mutant and mbs1/RNAi-mbs2 double mutants were shown to be more susceptible to HL conditions. Thus, while in wild-type plants β-CC causes induction of SORGs, resulting in greater photooxidative stress tolerance [34], the insertional knockdown mbs1 mutant exhibits deregulation in SORG expression after β-CC treatment and consequently does not achieve β-CC-induced HL tolerance [75]. In line with this idea, we exposed wild-type and mbs1 mutant plants to an episode of HL, and then we analyzed pri-miRNA expression levels (Fig. 4e; Additional file 1: Fig. S12). HL induces the expression of both pri-miR163 and pri-miR840 in Col-0, and the scale of this induction is much greater in pri-miR163. This response was noticeably reduced in mbs1 plants for both pri-miRNAs (Fig. 4e). The pri-miR163 induction was almost 40% weaker (13.7- versus 8.9-fold changes for Col-0 and mbs1, respectively), whereas pri-miR840 induction lost statistical significance in mutant plants (Fig. 4e). However, when we compared these data within the treatment, not within the genotype, we still observed a significant reduction in the expression level for pri-miR163 but not for pri-miR840 in HL (Additional file 1: Fig. S12). Only a partial effect of the mbs1 mutation on pri-miRNA induction suggests the existence of other parallel regulatory mechanisms or, more likely, the redundant role of a close MBS1 homolog – MBS2 [75]. Thus, chemical (β-CC) and partial genetic premises (mbs1 study) imply the potential role of 1O2 in the regulation of miRNA expression.

Fig. 4
figure 4

Regulation of pri-miRNA by β-CC-dependent 1O2 signaling. a Scheme represents the β-CC-dependent 1O2 signaling pathway induced in HL conditions. β-CC is formed in PSII as a result of β-carotene oxidative breakdown under HL conditions. MBS1 is a downstream component that transduces stress information to the nucleus, where it influences the expression of singlet oxygen responsive genes (SORGs). b qRT-PCR for the DRP gene and c pri-miR163 and pri-miR840 in control plants (0.00) and after β-CC treatment (0.05 ml and 1.00 ml) for 4 h. d TT-qRT PCR for miR163 and miR840 level in plants treated by β-cc. ctrl- control plants, trt—plants treated with 1.0 ml of β-CC. Transcript levels were normalized with respect to sno85 and sno101. e qRT-PCR for pri-miR163 and pri-miR840 in Col-0 and mbs1 plants. LLc-control plants; HL-plants exposed to HL for 2 h; LLr-plants exposed to HL for 2 h and subsequent recovery in LL for 4 h. Transcript levels were normalized with respect to the PP2A and UPL7 genes. Asterisks indicate significant differences according to the Tukey HSD test at the level of * ≤ 0.05, ** ≤ 0.01, *** ≤ 0.001. Mean values ± SDs (n = 3), were provided


Light affects posttranscriptional gene regulation

The ability to transfer stress signals from chloroplasts to the nucleus and consequently influence nuclear gene expression is indispensable for plant adaptive strategies. Thus, effective retrograde communication and its linkages with RNA metabolism have been extensively studied [76]. For example, it was demonstrated that changes in the PQ redox pool may regulate the alternative splicing process [38, 77, 78]. Recently, the basis of these mechanisms was revealed, showing that light increases the PolII elongation rate, while in darkness, elongation is reduced, which is consistent with the kinetic coupling model also found in mammalian cells [79, 80]. Light also influences other posttranscriptional gene expression regulation mechanisms, including the regulation of miRNAs [81, 82].

HL reveals incosistency in pri-miRNA and miRNA abundance

Light may affect miRNA expression at different stages of biogenesis. For example, light has a negative impact on miRNA processing (described as the phenomenon of inconsistency in miRNA levels) and is manifested by a reduction in miRNA response compared to significantly greater changes in pri-miRNA abundance, which was observed in the context of de-etiolation [83]. This light-dependent mechanism decreases the capability of miRNA processing, and activity of SMALL RNA DEGRADING NUCLEASE 1 (SDN1), which reduces the half-life of several miRNAs in de-etiolated seedlings, conferring relatively constant miRNA levels [83]. Recently, it was demonstrated that the forkhead-associated domain 2 (FHA2) protein is the light-stabilized suppressor of miRNAs through interactions with DCL1 and HYL1 [84, 85]. Briefly, FHA2 inhibits DCL1 enzyme activity and facilitates the accessibility of HYL1 to pri-miRNAs. Despite all of these reports, there is a knowledge gap concerning the role of light stress in miRNA processing, and the majority of current observations are based on nonstressed conditions. Certainly, there is some kind of negative regulation of miRNA processing during HL because we also observed inconsistency in miRNA abundance after HL treatment. For instance, upon HL, the fold changes for pri-miR163 ranged from 5 to 15 in Col-0 (Fig. 1e, Fig. 3d, Fig. 4e). Simultaneously, the changes of the mature miRNA did not exceed twofold, regardless of the experimental setup or methodology used (Fig. 1b,d, Additional file 1: Fig. S4). These results are compatible with those of experiments that used PET inhibitors for pri-miR163 (Fig. 2). Conversely, in the case of pri-miR840, mature miRNA changes were almost identical to those at the pri-miRNA level.

Functioning of the core components of miRNA biogenesis is disturbed by HL

We concluded that the discrepancy between the level of pri-miRNAs and their mature forms may depend on HYL1 activity because we observed accumulation of primary transcripts in the case of pri-miR163 in the hyl1 mutant but not for pri-miR840 (Additional file 1: Fig. S6). Interestingly, the mature miR840 level in hyl1 was several-fold greater than that in Col-0, indicating that miR840 biogenesis was more efficient in the absence of HYL1 (in the middle of the photoperiod light phase). This highlights the putative role of HYL1 as the microprocessor component underlying the regulation of miRNA expression by light. Since light may affect HYL1 localization and activity [37], such a conclusion seems justifiable. Nevertheless, there are some contradictory reports referring to light-dependent HYL1 regulation. For example, during the dark to red light transition, phytochrome-interacting factor 4 (PIF4) destabilized HYL1 and DCL1 [82, 86]. However, COP1 (constitutive photomorphogenic 1), a RING-finger E3 ligase and negative photomorphogenesis regulator, contributed to light-dependent proteolytic stabilization of the HYL1 protein during the daytime [87]. Additionally, HYL1 activity was maintained by its nucleo-cytoplasmic relocation and phosphorylation/dephosphorylation changes [82, 88, 89]. The nuclear pool of HYL1 in prolonged darkness is phosphorylated, which protects it from degradation. After plant reillumination, HYL1 is dephosphorylated to rapidly restore miRNA production [37]. In addition, the possible influence of chloroplast signals with respect to HYL1 posttranscriptional regulation was suggested because DCMU almost completely blocks HYL1 dephosphorylation after prolonged darkness [37]. Although there is some evidence of light-triggered HYL1 regulation, there is a lack of data about HL-governed changes in microprocessor functioning. Moreover, in addition to the importance of HYL1, DCL1 plays a dominant role in miRNA processing. Like HYL1, DCL1 is subject to multifaceted regulation [29]. For example, it was demonstrated that DCL1 highly accumulates in de-etiolation, but the processing activity decreases; consequently, the level of mature miRNAs is nearly unaltered [83]. Although we observed no changes in DCL1 accumulation after HL exposure (Additional file 1: Fig. S13), its activity is likely disturbed by this stress. This suspicion was supported by a dramatic decrease in miRNA levels after β-CC treatment (Fig. 4d). One of the plausible explanations is the negative regulation of microprocessors by overproduced ROS (major engagement of 1O2). Probably the applied dose of β-CC causes not only stronger pri-miRNAs induction but at the same time has much more severe impact on further steps of miRNAs biogenesis than HL. The exact mechanism by which HL stress adjusts the functionality of microprocessor components to suppress miRNA biogenesis still remains unknown. Further study on the role of retrograde signals may contribute to a deeper understanding of HL-regulated miRNA expression changes.

Retrograde 1O2, signaling is important in miRNA regulation

To date, knowledge about the role of plastid-derived signals in miRNA regulation in the context of light signaling, in addition to the possible role of DCMU in changing the HYL1 phosphorylation status, is still vague and limited to a few but significant reports. For example, the involvement of the plastid-nucleus located DNA/RNA binding protein WHIRLY1 in miRNA expression regulation during light stress in barley was postulated [90]. sRNA sequencing combined with mRNA/lncRNA sequencing on Arabidopsis wild-type plants and two retrograde signaling mutants, gun1 and gun5, treated with the herbicide norflurazon shed light on cross-talk between sRNAs (including miRNAs) and retrograde signaling [91]. Additionally, the functional role of chloroplast-derived signals in miRNA regulation was postulated recently by Fang et al., 2019 [92], in the context of heat stress. They showed that tocopherols enhance the accumulation of 3’-phosphoadenosine 5’ phosphate (PAP), a retrograde nuclear exonuclease 2 (XRN2) inhibitor. Thus, tocopherols and PAPs are positive regulators of miRNA biogenesis because their accumulation represses XRN2, which negatively regulates mRNA and pri-miRNA levels by the degradation of uncapped 5’ mRNA. Since tocopherols (together with carotenoids) are reported to protect PSII against photoinhibition and lipid peroxidation by ROS, questions about the role of oxidative stress (including 1O2) in this regulation have been raised [61, 93, 94]. This is especially important because PAP accumulates during drought- and HL-induced oxidative stress, and its level is regulated by the SAL1 enzyme, which dephosphorylates PAP to AMP [95]. Interestingly, the alx8 mutant, which accumulates PAP, has a significantly higher level of PXMT1, which encodes a 1,7-paraxanthine methyltransferase, a target gene of miR163 [95]. Since PAP effectively suppresses not only nuclear XRN2 and XRN3 but also cytoplasmic XRN4, which functions in the degradation of miRNA target cleavage products, this result is reasonable [96]. In addition, alx8 also exhibited an increase in the pri-miR163 level, with a simultaneous decrease in mature miRNA forms (Aditional file 1: Fig. S14).Considering that the increased pri-miRNA level in HL can be controlled not only by inhibition of nuclear XRNs but also by a higher pri-miRNA transcription rate, we performed PolII-ChIP analysis. There were no changes in PolII occupancy at the tested MIR genes between LL and HL treatments, arguing against the second possibility (Fig. 1f). However, promoter fusion with the beta-glucoronidase gene showed enhanced transcription of pri-miR163 in HL (Additional file 1: Fig. S7). For miR163, there is a report that HY5 regulates light-responsive transcription of miR163 [97]. The lack of changes in PolII occupancy seem to contradict these results; however, PolII-ChIP measures the transcription rate as a function of PolII occupancy at the time of analysis – 2 h of HL. Nonetheless, there is the possibility that the transcription rate for miR163 and miR840 is in fact changed at the beginning of HL treatment because pri-miRNA expression of those miRNAs was even higher in plants exposed to HL for 30 min (Additional file 1: Fig. S15). Since the DRP level was also more abundant after 30 min of treatment, we consequently considered 1O2 to be a putative signal that induces miRNA expression. Notably, the effective transcription rate can vary even when PolII occupancy at tested genes is stable because PolII can pause or even stall [98]. Perhaps the pri-miRNA level is determined by the interplay between pri-miRNA transcription and stability; therefore, we also verified the half-life of these pri-miRNAs using a cordycepin assay (Fig. 1g-h). The obtained results suggest that pri-miR163 is more stable in HL, while the stability of pri-miR840 is reduced (Fig. 1g-h). The differences between pri-miR163 and pri-miR840 degradation in HL consistently indicate that HYL1 may be an important determinant of the observed changes. It is possible that HL causes an increase in the transcription of some miRNAs and simultaneously disturbs the activity of HYL1 and DCL1. Collectively, our results suggest that changes in pri-miRNA levels are tuned by changes in pri-miRNA transcription and stability. Although we present the link between miRNA expression and ROS (i.e., 1O2) signaling upon HL stress, the involvement of other retrograde signals remains to be investigated, especially in the context of miRNA biogenesis regulation. Further investigation of these aspects will lead to a better understanding of how miRNA expression and processing machinery cooperate upon HL. Such research should include the biological function of miRNA targets and their evolutionary context.


In summary, we provided several lines of evidence that pinpoint the function of 1O2 in the regulation of miRNA expression during HL conditions. The negative impact of high light intensity on some mature miRNA levels and simultaneous accumulation of their primary transcripts was also observed. We postulate that microprocessor component activity may be negatively regulated by HL, similar to the de-etiolation process. Our work provides deeper insight into the cross-talk between retrograde signaling and miRNA expression, creating new perspectives for further studies.


Plant materials

Arabidopsis thaliana Columbia-0 ecotype wild-type plants seeds were obtained from the Nottingham Arabidopsis Stock Centre (NASC, NASC code: N76778) and ex1 (SALK_022735, NASC code: N522735) seeds. flu and flu/ex1 seeds were kindly provided by C. Kim [67], while mbs1 mutant seeds (SAIL_661_B05) were received from N. Shao [75]. The exact times at which plant materials were harvested for analysis are presented in Additional file 1: Figures S3, S6,S8, S10-S12.

Hydroponic conditions

A. thaliana plants were grown for 4 weeks in hydroponic conditions [99] with some previously described modifications [13]. Briefly, seeds were surface sterilized using a chlorine gas method and kept for 2 days at 4 °C on agarose in high humidity to synchronize germination. Controlled growth conditions included a short-day photoperiod (8 h light/16 h dark), 22 °C/20 °C (day/night), 70% air humidity, and low light intensity (LL; 100–120 µmol m−2 s −1).

High-light treatment procedures

HL treatment was applied by exposing LL-adapted plants to high light (HL) intensity for 2 h using LED light sources (HL; 1500 µmol photons m−2 s−1; Photon Systems Instruments, Brno, Czech Republic).

Chlorophyll a fluorescence

Chlorophyll a fluorescence parameters were determined using a PAM FluorCam 800 MF PSI device (Brno, Czech Republic). The plants were transferred to darkness for 30 min prior to measurement. The Fv/Fm parameter, which reflects PSII maximum efficiency, was measured.

Microtranscriptomic sequencing

Four-week-old A. thaliana plants grown in hydroponic conditions were subjected to HL stress. Plant material was sampled just after exposure to HL stress conditions (sample name: HL) and after recovery in LL conditions for 4 h (LLr). RNA isolation was performed using the Universal RNA/miRNA purification kit (EURX, cat. no E3599 Gdańsk, Poland) according to the manufacturer’s instructions. miRNA library preparation, miRNA sequencing and data analysis were performed by GENOMED S. A (Warsaw, Poland). Briefly, miRNA libraries were prepared using the NEBNext® Small RNA Library Prep Set for Illumina® (Multiplex Compatible) and sequenced using the Illumina HiSeq 4000 platform (Illumina Inc., San Diego, CA, USA). Bioinformatic analysis was also outsourced and conducted as previously described in Barczak-Brzyżek et al. 2019 [13]. Briefly, quality control checks of raw sequence data were performed using the FASTQ tool. Then, for trimming adapters, the Cutadapt program was used with subsequent identification of novel and known miRNAs using miRDeep2. Finally, the EdgeR Bioconductor package was applied for differential expression analysis. Microtranscriptomic screening results are presented for 3 independent experiments. In all cases, each biological replicate was pooled from six plants. The trimmed sequence data were deposited in the SRA database under accession number PRJNA650313.

Chromatin immunoprecipitation (ChIP)

4-week-old A. thaliana Col-0 plants grown in hydroponic conditions were subjected to 2 h of HL treatment. Plant material was harvested to analysis at 1330. In all cases, each biological replicate was pooled from approx. 20 plants. Chromatin immunoprecipitation was performed as described in Godoy-Herz et al. 2019 [79]. IP buffer was prepared based on Kaufmann et al. 2010 [100]. Plant material was crosslinked using formaldehyde and then grinded with liquid nitrogen. Subsequently, chromatin was isolated and then sonicated before proceeding to immunoprecipitation. Next antibodies against total Pol II (Agrisera AS11 1804) were used with Dynabeads Protein G (Invitrogen, 10003D). Chelex (Biorad, 1421253) was used for de-crosslinking as described in Nelson et al. 2009 [101]. No antibody control was used to determine nonspecific background and percentage of input was calculated for each sample using qPCR.

RNA stability assay

A cordycepin RNA stability assay was performed as described before in Fedak et al. 2016 [102]. Arabidopsis seedlings were grown for 2 weeks in LL (SD; 8 h light/16 h dark, temperature 22 °C/20 °C) on Murashige and Skoog (½ MS) medium (Duchefa Biochemie, cat. no M0222), supplemented with 1% w/v Sucrose (Duchefa Biochemie, cat. no S0809 and 0.7% phytoagar (Duchefa Biochemie, cat. no P1003), pH 5.7. Seeds were stratified at 4 °C for 2 days after sowing on Petri dishes. Seedlings were kept in LL (growth chamber) or HL conditions (1 h). Plant material was harvested to analysis from 1210 to 1330. Seedlings were collected and transferred to a flask containing incubation buffer (1 mM Pipes, pH 6.25, 1 mM trisodium citrate, 1 mM KCl, 15 mM sucrose. After 15 min of incubation, cordycepin was added to a final concentration of 150 µg/mL and vacuum-infilitrated (approx.2 × 5 min). At each time point (0, 20, 40, 60, 80 min), seedlings representing approx. 0.05 g were collected and frozen in liquid nitrogen. Samples were analyzed in triplicate. RNA extraction was performed using TRIzol method. RT-qPCR. analysis was performed with primers listed in Addifional file 1:Table S2. Calculation of pri-miRNAs half-life was performed as described in Chen et al. 2008 [103].

DCMU and DBMIB treatments

Four-week-old A. thaliana plants grown in hydroponic conditions (see above) were treated with inhibitors. DCMU and DBMIB stock solutions (30 mM) were prepared by dissolving DCMU (3-(3,4-dichlorophenyl)-1,1-dimethylurea, Sigma Aldrich cat. no D2425) and DBMIB (2,5- dibromo-6-isopropyl-3-methyl-1,4-benzoquinone, Sigma Aldrich cat. no 271993) in DMSO (dimethyl sulfoxide, Sigma Aldrich cat. no 8418). DCMU and DBMIB working solutions were prepared by diluting appropriate stock solutions with sterilized water to a final concentration of 30 µM. For control treatments, plants were sprayed with 0.1% DMSO solution. In all cases, each biological replicate was pooled from six plants.

EX1-dependent 1O2 signaling study

For the results presented in Fig. 3a-c ex1, flu and flu/ex1 plants were grown for 2 weeks on 90 mm diameter Petri dishes on half-strength Murashige and Skoog (½ MS) medium (Duchefa Biochemie, cat. no M0222) supplemented with 1% w/v sucrose (Duchefa Biochemie, cat. no S0809 and 0.7% phytoagar (Duchefa Biochemie, cat. no P1003), pH 5.7. Induction of 1O2 accumulation was achieved in flu background plants by transferring plants cultivated in CL (constant light: LL intensity 90–110 µmol m−2 s −1, temperature 20 °C humidity 70%) to the dark for 12 h followed by 2 h reillumination. In all cases, each biological replicate was pooled from at least six plants.

For the results presented in Fig. 3d, 4-week-old A. thaliana Col-0 and ex1 plants grown in hydroponic conditions were subjected to HL treatment. In all cases, each biological replicate was pooled from six plants.

β-CC-dependent 1O2 signaling study

For the results presented in Fig. 4b-d, 3- to 5-week-old A. thaliana plants were grown in pots under controlled conditions with a long photoperiod (light intensity approx. 250 µmol photons m−2 s−1, temperature 20 °C and humidity 70%). β-CC (β-cyclocitral; Santa Cruz Biotechnology, cat no sc-207467) treatment was performed as previously described [30]. Briefly, plants were placed for 4 h in a transparent plexiglass box (approx. 15 l vol.) with defined volumes (50 µl and 1 ml) of β-CC applied on a cotton wick to increase the contact area with the air. For the control conditions, the β-CC was replaced by distilled water. In all cases, each biological replicate was pooled from six plants.

For the results presented in Fig. 4e, 4-week-old A. thaliana Col-0 and mbs1 plants grown in hydroponic conditions were subjected to HL treatment. In all cases, each biological replicate was pooled from six plants.

RNA preparation

RNA extraction was performed using the Universal RNA/miRNA purification kit EURX (EURX, cat. no E3599, Gdańsk, Poland) according to the manufacturer’s instructions (for results presented in Fig. 1b-e and Fig. 2c) or using TRIzol reagent (Invitrogen, cat. no 15596062) (for results presented in Fig. 3b-d and Fig. 4b-e). The obtained RNA was treated with TurboDNase (Invitrogen, cat. no AM2238) according to the manufacturer’s recommendation. The RNA concentration was estimated using a NanoDrop 1000 (Thermo Fischer Scientific, Wilmington, MA, USA).

RT-qPCR for mRNA quantification

cDNA synthesis

A Quantitect Reverse Transcription kit (Qiagen, cat. no 205313, Hindel, Germany) was used for cDNA synthesis (for results presented in Fig. 1e and Fig. 2c) according to the manufacturer’s instructions. Oligo dT-primed cDNA was prepared using Superscript III reverse transcriptase (Invitrogen cat. no 18080093; for results presented in Fig. 3b-d and Fig. 4b,c,e).

Quantitative RT PCR

qRT-PCR was performed in triplicate using the Bio-Rad CFX96 Touch ™ Real-Time PCR Detection System (Bio-Rad, Hercules, CA United States)/Roche Light Cycler® 384 System (Basel, Switzerland) with the primers listed in Additional file 1: Table S2. Real-time PCR cycling conditions were optimized depending on the primer used in the protocol, and relative expression was calculated with reference to the UPL7 (AT3G53090) and PP2A (AT1G13320) genes. Product melting curves were generated following PCR to ensure amplification product purity.

Reverse transcription for miRNAs was performed with the qScript flex cDNA synthesis kit (Quantabio, cat. no 95047, Beverly, USA) according to [46] in a 10 µl total reaction volume. RNA was diluted in TE-LPA buffer (TE buffer with Linear Poly-Acrylamide (Ambion®, cat no AM9520) at a final working concentration of 20 µg/mL). The RT reaction mixture contained 10 ng of total RNA, 1 × RT buffer, 0.05 µM RT gene-specific primer (see primers list Additional file 1: Table S2), 1 µl GSP enhancer, and 0.5 µl RT enzyme. RT reactions were incubated in PCR tubes for 45 min at 25 °C and 5 min at 85 °C and then held at 4 °C.

Two-tailed qRT PCR (TT qRT PCR)

TT qRT PCR was performed according to Androvic et al. 2017 [46], as previously described in Barczak-Brzyżek et al. 2019 [13]. The primers used in the study were designed in collaboration with BIOCEV, Institute of Biotechnology CAS, Czech Republic. Specifically, 1 × SYBR, 0.4 µM forward and reverse primer (see primers list Additional file 1:Table S2), and the cDNA product diluted 5 × were mixed in a 10 µL total reaction volume. Reactions were performed in triplicate and incubated in 96-well plates (CFX 96 Real Time Detection System (Bio-Rad)) at 95 °C for 30 s, followed by 45 cycles of 95 °C for 5 s and 60 °C for 15 s. Reaction specificity was assessed by melting curve analysis. Expression levels were calculated relative to those of snoRNA85 (NCBI Accession Number AJ505658) and snoRNA101 (NCBI Accession Number AJ505631).

Statistical analysis

Statistical analysis was performed in R (version 3.6.0).

Availability of data and materials

The trimmed data were deposited in the SRA database under accession number PRJNA650313. All the other data supporting the results of this article are included within the paper and its supplementary file as figures or tables.


1O2 :

Singlet oxygen


Chromatin immunoprecipitation


2,5- Dibromo-6-isopropyl-3-methyl-1,4-benzoquinone




Dicer-like 1


Executer 1


Grana margins


High light


Low light


Methylene Blue Sensitivity 1


Micro RNAs


3’-Phosphoadenosine 5’ phosphate


Polymerase chain reaction


Photosynthetic electron transport


RNA polymerase II






Primary micro RNAs


Photosystem I


Photosystem II


Quantitative reverse transcription PCR


RNA induced silencing complex


Reactive oxygen species


Small nucleolar85 RNA


Small nucleolar101 RNA


Singlet oxygen responsive genes


Transcription factor


Two-tailed qRT PCR


β- Cyclocitral


  1. Karpinski S, Reynolds H, Karpinska B, Wingsle G, Creissen G, Mullineaux P. Systemic signaling and acclimation in response to excess excitation energy in Arabidopsis. Science. 1999;284(5414):654–7.

    CAS  PubMed  Google Scholar 

  2. Murata N, Takahashi S, Nishiyama Y, Allakhverdiev SI. Photoinhibition of photosystem II under environmental stress. Biochim Biophys Acta. 2007;1767(6):414–21.

    CAS  PubMed  Google Scholar 

  3. Chan KX, Phua SY, Crisp P, McQuinn R, Pogson BJ. Learning the Languages of the Chloroplast: Retrograde Signaling and Beyond. Annu Rev Plant Biol. 2016;67(1):25–53.

    CAS  PubMed  Google Scholar 

  4. Karpiński S, Szechyńska-Hebda M, Wituszyńska W, Burdiak P. Light acclimation, retrograde signalling, cell death and immune defences in plants. Plant Cell Environ. 2013;36(4):736–44.

    PubMed  Google Scholar 

  5. Leister D. Retrograde signaling in plants: from simple to complex scenarios. Front Plant Sci. 2012;3:135.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Barciszewska-Pacak M, Milanowska K, Knop K, Bielewicz D, Nuc P, Plewka P, et al. Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses. Front Plant Sci. 2015;6:410.

    PubMed  PubMed Central  Google Scholar 

  7. Achkar NP, Cambiagno DA, Manavella PA. miRNA Biogenesis: A Dynamic Pathway. Trends Plant Sci. 2016;21(12):1034–44.

    CAS  PubMed  Google Scholar 

  8. Huot B, Yao J, Montgomery BL, He SY. Growth-Defense Tradeoffs in Plants: A Balancing Act to Optimize Fitness. Mol Plant. 2014;7(8):1267–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 2013;25(7):2383–99.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Li Y, Varala K, Hudson ME. A survey of the small RNA population during far-red light-induced apical hook opening. Front Plant Sci. 2014;5:156.

    PubMed  PubMed Central  Google Scholar 

  11. Zhao H, Chen D, Peng Z, Wang L, Gao Z. Identification and characterization of microRNAs in the leaf of ma bamboo (Dendrocalamus latiflorus) by deep sequencing. PLoS ONE. 2013;8(10):e78755.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Tiwari B, Habermann K, Arif MA, Top O, Frank W. Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana. Front Plant Sci. 2021;12:1179.

    Google Scholar 

  13. Anna B-B, Grzegorz B, Marek K, Piotr G, Marcin F. Exposure to High-Intensity Light Systemically Induces Micro-Transcriptomic Changes in Arabidopsis thaliana Roots. Int J Mol Sci. 2019;20(20):5131.

    CAS  PubMed Central  Google Scholar 

  14. Casati P. Analysis of UV-B regulated miRNAs and their targets in maize leaves. Plant Signal Behav. 2013;8(10).

  15. Jia X, Ren L, Chen Q-J, Li R, Tang G. UV-B-responsive microRNAs in Populus tremula. J Plant Physiol. 2009;166(18):2046–57.

    CAS  PubMed  Google Scholar 

  16. Kim JH, Go YS, Kim JK, Chung BY. Characterization of microRNAs and their target genes associated with transcriptomic changes in gamma-irradiated Arabidopsis. Genet Mol Res. 2016;15(3).

  17. Sunitha S, Loyola R, Alcalde JA, Arce-Johnson P, Matus JT, Rock CD. The Role of UV-B light on Small RNA Activity During Grapevine Berry Development. G3 Genes Genomes Genet. 2019;9(3):769–87.

    CAS  Google Scholar 

  18. Zhou X, Wang G, Zhang W. UV-B responsive microRNA genes in Arabidopsis thaliana. Mol Syst Biol. 2007;3:103.

    PubMed  PubMed Central  Google Scholar 

  19. Megraw M, Baev V, Rusinov V, Jensen ST, Kalantidis K, Hatzigeorgiou AG. MicroRNA promoter element discovery in Arabidopsis. RNA. 2006;12(9):1612–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC. Expression of Arabidopsis MIRNA Genes. Plant Physiol. 2005;138(4):2145–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhao X, Zhang H, Li L. Identification and analysis of the proximal promoters of microRNA genes in Arabidopsis. Genomics. 2013;101(3):187–94.

    CAS  PubMed  Google Scholar 

  22. Kurihara Y, Watanabe Y. Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci U S A. 2004;101(34):12753–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Kurihara Y, Takashi Y, Watanabe Y. The interaction between DCL1 and HYL1 is important for efficient and precise processing of pri-miRNA in plant microRNA biogenesis. RNA N Y N. 2006;12(2):206–12.

    CAS  Google Scholar 

  24. 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(6):841–50.

    CAS  PubMed  Google Scholar 

  25. Yu Y, Jia T, Chen X. The, “how” and “where” of plant microRNAs. New Phytol. 2017;216(4):1002–17.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Pratt AJ, MacRae IJ. The RNA-induced silencing complex: a versatile gene-silencing machine. J Biol Chem. 2009;284(27):17897–901.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Szweykowska-Kulinska Z, Jarmolowski A. Post-transcriptional Regulation of MicroRNA Accumulation and Function: New Insights from Plants. Mol Plant. 2018;11(8):1006–7.

    CAS  PubMed  Google Scholar 

  28. Bhat SS, Bielewicz D, Gulanicz T, Bodi Z, Yu X, Anderson SJ, et al. mRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana. Proc Natl Acad Sci. 2020;117(35):21785–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Dolata J, Taube M, Bajczyk M, Jarmolowski A, Szweykowska-Kulinska Z, Bielewicz D. Regulation of Plant Microprocessor Function in Shaping microRNA Landscape. Front Plant Sci. 2018;9:753.

    PubMed  PubMed Central  Google Scholar 

  30. Foyer CH, Shigeoka S. Understanding Oxidative Stress and Antioxidant Functions to Enhance Photosynthesis. Plant Physiol. 2011;155(1):93–100.

    CAS  PubMed  Google Scholar 

  31. Wagner D, Przybyla D, Op den Camp R, Kim C, Landgraf F, Lee KP, et al. The genetic basis of singlet oxygen-induced stress responses of Arabidopsis thaliana. Science. 2004;306(5699):1183–5.

    CAS  PubMed  Google Scholar 

  32. Wang L, Leister D, Guan L, Zheng Y, Schneider K, Lehmann M, et al. The Arabidopsis SAFEGUARD1 suppresses singlet oxygen-induced stress responses by protecting grana margins. Proc Natl Acad Sci U S A. 2020;117(12):6918-6927.

  33. Woodson JD, Joens MS, Sinson AB, Gilkerson J, Salomé PA, Weigel D, et al. Ubiquitin facilitates a quality-control pathway that removes damaged chloroplasts. Science. 2015;350(6259):450–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Ramel F, Birtic S, Ginies C, Soubigou-Taconnat L, Triantaphylidès C, Havaux M. Carotenoid oxidation products are stress signals that mediate gene responses to singlet oxygen in plants. Proc Natl Acad Sci U S A. 2012;109(14):5535–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Shumbe L, Chevalier A, Legeret B, Taconnat L, Monnet F, Havaux M. Singlet Oxygen-Induced Cell Death in Arabidopsis under High-Light Stress Is Controlled by OXI1 Kinase. Plant Physiol. 2016;170(3):1757–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Jung H-S, Crisp PA, Estavillo GM, Cole B, Hong F, Mockler TC, et al. Subset of heat-shock transcription factors required for the early response of Arabidopsis to excess light. Proc Natl Acad Sci. 2013;110(35):14474–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Achkar NP, Cho SK, Poulsen C, Arce AL, Re DA, Giudicatti AJ, et al. A Quick HYL1-Dependent Reactivation of MicroRNA Production Is Required for a Proper Developmental Response after Extended Periods of Light Deprivation. Dev Cell. 2018;46(2):236-247.e6.

    CAS  PubMed  Google Scholar 

  38. Petrillo E, Herz MAG, Fuchs A, Reifer D, Fuller J, Yanovsky MJ, et al. A chloroplast retrograde signal regulates nuclear alternative splicing. Science. 2014;344(6182):427–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Górecka M, Lewandowska M, Dąbrowska-Bronk J, Białasek M, Barczak-Brzyżek A, Kulasek M, et al. Photosystem II 22kDa protein level - a prerequisite for excess light-inducible memory, cross-tolerance to UV-C and regulation of electrical signalling. Plant Cell Environ. 2020;43(3):649–61.

    PubMed  Google Scholar 

  40. Szechyńska-Hebda M, Kruk J, Górecka M, Karpińska B, Karpiński S. Evidence for Light Wavelength-Specific Photoelectrophysiological Signaling and Memory of Excess Light Episodes in Arabidopsis. Plant Cell. 2010;22(7):2201–18.

    PubMed  PubMed Central  Google Scholar 

  41. Vogel MO, Moore M, König K, Pecher P, Alsharafa K, Lee J, et al. Fast retrograde signaling in response to high light involves metabolite export, MITOGEN-ACTIVATED PROTEIN KINASE6, and AP2/ERF transcription factors in Arabidopsis. Plant Cell. 2014;26(3):1151–65.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Gordon MJ, Carmody ME, Albrecht V, Pogson B. Systemic and Local Responses to Repeated HL Stress-Induced Retrograde Signaling in Arabidopsis. Front Plant Sci. 2013;3:303.

  43. Karpinski S, Escobar C, Karpinska B, Creissen G, Mullineaux PM. Photosynthetic electron transport regulates the expression of cytosolic ascorbate peroxidase genes in Arabidopsis during excess light stress. Plant Cell. 1997;9(4):627–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. DaiZhao XPX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39(Web Server issue):W155-159.

    Google Scholar 

  45. Szczesniak MW, Makalowska I. miRNEST 2.0: a database of plant and animal microRNAs. Nucleic Acids Res. 2014;42(Database issue):D74-D77.

  46. Androvic P, Valihrach L, Elling J, Sjoback R, Kubista M. Two-tailed RT-qPCR: a novel method for highly accurate miRNA quantification. Nucleic Acids Res. 2017;45(15):e144.

    PubMed  PubMed Central  Google Scholar 

  47. Chung PJ, Park BS, Wang H, Liu J, Jang I-C, Chua N-H. Light-Inducible MiR163 Targets PXMT1 Transcripts to Promote Seed Germination and Primary Root Elongation in Arabidopsis. Plant Physiol. 2016;170(3):1772–82.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Mancini E, Sanchez SE, Romanowski A, Schlaen RG, Sanchez-Lamas M, Cerdán PD, et al. Acute Effects of Light on Alternative Splicing in Light-Grown Plants. Photochem Photobiol. 2016;92(1):126–33.

    CAS  PubMed  Google Scholar 

  49. Shikata H, Hanada K, Ushijima T, Nakashima M, Suzuki Y, Matsushita T. Phytochrome controls alternative splicing to mediate light responses in Arabidopsis. Proc Natl Acad Sci U S A. 2014;111(52):18781–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Song L, Han M-H, Lesicka J, Fedoroff N. Arabidopsis primary microRNA processing proteins HYL1 and DCL1 define a nuclear body distinct from the Cajal body. Proc Natl Acad Sci U S A. 2007;104(13):5437–42.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Szarzynska B, Sobkowiak L, Pant BD, Balazadeh S, Scheible W-R, Mueller-Roeber B, et al. Gene structures and processing of Arabidopsis thaliana HYL1-dependent pri-miRNAs. Nucleic Acids Res. 2009;37(9):3083–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Ré DA, Lang PLM, Yones C, Arce AL, Stegmayer G, Milone D, et al. Alternative use of miRNA-biogenesis co-factors in plants at low temperatures. Development. 2019;146(5):dev172932.

    PubMed  Google Scholar 

  53. Dietz K-J. Redox signal integration: from stimulus to networks and genes. Physiol Plant. 2008;133(3):459–68.

    CAS  PubMed  Google Scholar 

  54. Fey V, Wagner R, Braütigam K, Wirtz M, Hell R, Dietzmann A, et al. Retrograde plastid redox signals in the expression of nuclear genes for chloroplast proteins of Arabidopsis thaliana. J Biol Chem. 2005;280(7):5318–28.

    CAS  PubMed  Google Scholar 

  55. Pfannschmidt T, Schütze K, Brost M, Oelmüller R. A novel mechanism of nuclear photosynthesis gene regulation by redox signals from the chloroplast during photosystem stoichiometry adjustment. J Biol Chem. 2001;276(39):36125–30.

    CAS  PubMed  Google Scholar 

  56. Khandelwal A, Elvitigala T, Ghosh B, Quatrano RS. Arabidopsis transcriptome reveals control circuits regulating redox homeostasis and the role of an AP2 transcription factor. Plant Physiol. 2008;148(4):2050–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Schansker G, Tóth SZ, Strasser RJ. Methylviologen and dibromothymoquinone treatments of pea leaves reveal the role of photosystem I in the Chl a fluorescence rise OJIP. Biochim Biophys Acta. 2005;1706(3):250–61.

    CAS  PubMed  Google Scholar 

  58. Gawroński P, Górecka M, Bederska M, Rusaczonek A, Ślesak I, Kruk J, et al. Isochorismate synthase 1 is required for thylakoid organization, optimal plastoquinone redox status, and state transitions in Arabidopsis thaliana. J Exp Bot. 2013;64(12):3669–79.

    PubMed  PubMed Central  Google Scholar 

  59. Dogra V, Kim C. Chloroplast protein homeostasis is coupled with retrograde signaling. Plant Signal Behav. 2019;14(11):1656037.

    PubMed  PubMed Central  Google Scholar 

  60. Dogra V, Kim C. Singlet Oxygen Metabolism: From Genesis to Signaling. Front Plant Sci. 2020;10:1640.

  61. Krieger-Liszkay A. Singlet oxygen production in photosynthesis. J Exp Bot. 2005;56(411):337–46.

    CAS  PubMed  Google Scholar 

  62. Krieger-Liszkay A, Fufezan C, Trebst A. Singlet oxygen production in photosystem II and related protection mechanism. Photosynth Res. 2008;98(1–3):551–64.

    CAS  PubMed  Google Scholar 

  63. Triantaphylidès C, Havaux M. Singlet oxygen in plants: production, detoxification and signaling. Trends Plant Sci. 2009;14(4):219–28.

    PubMed  Google Scholar 

  64. Dogra V, Duan J, Lee KP, Lv S, Liu R, Kim C. FtsH2-Dependent Proteolysis of EXECUTER1 Is Essential in Mediating Singlet Oxygen-Triggered Retrograde Signaling in Arabidopsis thaliana. Front Plant Sci. 2017;8:1145.

    PubMed  PubMed Central  Google Scholar 

  65. Dogra V, Li M, Singh S, Li M, Kim C. Oxidative post-translational modification of EXECUTER1 is required for singlet oxygen sensing in plastids. Nat Commun. 2019;10(1):2834.

    PubMed  PubMed Central  Google Scholar 

  66. Lee KP, Kim C, Landgraf F, Apel K. EXECUTER1- and EXECUTER2-dependent transfer of stress-related signals from the plastid to the nucleus of Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2007;104(24):10270–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Wang L, Kim C, Xu X, Piskurewicz U, Dogra V, Singh S, et al. Singlet oxygen- and EXECUTER1-mediated signaling is initiated in grana margins and depends on the protease FtsH2. Proc Natl Acad Sci U S A. 2016;113(26):E3792-3800.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. op den Camp RGL, Przybyla D, Ochsenbein C, Laloi C, Kim C, Danon A, et al. Rapid Induction of Distinct Stress Responses after the Release of Singlet Oxygen in Arabidopsis. Plant Cell. 2003;15(10):2320–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. Meskauskiene R, Nater M, Goslings D, Kessler F, op den Camp R, Apel K. FLU: A negative regulator of chlorophyll biosynthesis in Arabidopsis thaliana. Proc Natl Acad Sci. 2001;98(22):12826–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  70. Zhang S, Apel K, Kim C. Singlet oxygen-mediated and EXECUTER-dependent signalling and acclimation of Arabidopsis thaliana exposed to light stress. Philos Trans R Soc Lond B Biol Sci. 2014;369(1640):20130227.

    PubMed  PubMed Central  Google Scholar 

  71. Kim C, Meskauskiene R, Zhang S, Lee KP, Lakshmanan Ashok M, Blajecka K, et al. Chloroplasts of Arabidopsis are the source and a primary target of a plant-specific programmed cell death signaling pathway. Plant Cell. 2012;24(7):3026–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. Gawroński P, Witoń D, Vashutina K, Bederska M, Betliński B, Rusaczonek A, et al. Mitogen-activated protein kinase 4 is a salicylic acid-independent regulator of growth but not of photosynthesis in Arabidopsis. Mol Plant. 2014;7(7):1151–66.

    PubMed  Google Scholar 

  73. Ramel F, Birtic S, Cuiné S, Triantaphylidès C, Ravanat J-L, Havaux M. Chemical Quenching of Singlet Oxygen by Carotenoids in Plants. Plant Physiol. 2012;158(3):1267–78.

    CAS  PubMed  PubMed Central  Google Scholar 

  74. Shao N, Duan GY, Bock R. A mediator of singlet oxygen responses in Chlamydomonas reinhardtii and Arabidopsis identified by a luciferase-based genetic screen in algal cells. Plant Cell. 2013;25(10):4209–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. Shumbe L, D’Alessandro S, Shao N, Chevalier A, Ksas B, Bock R, et al. METHYLENE BLUE SENSITIVITY 1 (MBS1) is required for acclimation of Arabidopsis to singlet oxygen and acts downstream of β-cyclocitral. Plant Cell Environ. 2017;40(2):216–26.

    CAS  PubMed  Google Scholar 

  76. Zhao X, Huang J, Chory J. Unraveling the Linkage between Retrograde Signaling and RNA Metabolism in Plants. Trends Plant Sci. 2020;25(2):141–7.

    CAS  PubMed  Google Scholar 

  77. Godoy MA, Herz AR, Kornblihtt, Barta A, Kalyna M, Petrillo E. Shedding light on the chloroplast as a remote control of nuclear gene expression. Plant Signal Behav. 2014;9(11):e976150.

    Google Scholar 

  78. Riegler S, Servi L, Scarpin MR, Godoy Herz MA, Kubaczka MG, Venhuizen P, et al. Light regulates alternative splicing outcomes via the TOR kinase pathway. Cell Rep. 2021;36(10):109676.

    CAS  PubMed  PubMed Central  Google Scholar 

  79. Godoy Herz MA, Kubaczka MG, Brzyżek G, Servi L, Krzyszton M, Simpson C, et al. Light Regulates Plant Alternative Splicing through the Control of Transcriptional Elongation. Mol Cell. 2019;73(5):1066-1074.e3.

    CAS  PubMed  Google Scholar 

  80. de la Mata M, Alonso CR, Kadener S, Fededa JP, Blaustein M, Pelisch F, et al. A slow RNA polymerase II affects alternative splicing in vivo. Mol Cell. 2003;12(2):525–32.

    PubMed  Google Scholar 

  81. Hernando CE, Garcia C, Mateos JL. Casting Away the Shadows: Elucidating the Role of Light-mediated Posttranscriptional Control in Plants. Photochem Photobiol. 2017;93(3):656–65.

    CAS  PubMed  Google Scholar 

  82. Manavella PA, Yang SW, Palatnik J. Keep calm and carry on: miRNA biogenesis under stress. Plant J. 2019;99(5):832–43.

    CAS  PubMed  Google Scholar 

  83. Choi SW, Ryu MY, Viczián A, Jung HJ, Kim GM, Arce AL, et al. Light Triggers the miRNA-Biogenetic Inconsistency for De-etiolated Seedling Survivability in Arabidopsis thaliana. Mol Plant. 2020;13(3):431–45.

  84. Park SJ, Choi SW, Kim GM, Møller C, Pai H-S, Yang SW. Light-stabilized FHA2 suppresses miRNA biogenesis through interactions with DCL1 and HYL1. Mol Plant. 2021;14(4):647–63.

    CAS  PubMed  Google Scholar 

  85. Gan L, Yu B. Forkhead-associated domain 2 links light signal to miRNA biogenesis. Mol Plant. 2021;14(4):541–3.

    CAS  PubMed  Google Scholar 

  86. Sun Z, Li M, Zhou Y, Guo T, Liu Y, Zhang H, et al. Coordinated regulation of Arabidopsis microRNA biogenesis and red light signaling through Dicer-like 1 and phytochrome-interacting factor 4. PLOS Genet. 2018;14(3):e1007247.

    PubMed  PubMed Central  Google Scholar 

  87. Cho SK, Ben Chaabane S, Shah P, Poulsen CP, Yang SW. COP1 E3 ligase protects HYL1 to retain microRNA biogenesis. Nat Commun. 2014;23(5):5867.

    Google Scholar 

  88. Karlsson P, Christie MD, Seymour DK, Wang H, Wang X, Hagmann J, et al. KH domain protein RCF3 is a tissue-biased regulator of the plant miRNA biogenesis cofactor HYL1. Proc Natl Acad Sci. 2015;112(45):14096–101.

    CAS  PubMed  PubMed Central  Google Scholar 

  89. Manavella PA, Hagmann J, Ott F, Laubinger S, Franz M, Macek B, et al. Fast-forward genetics identifies plant CPL phosphatases as regulators of miRNA processing factor HYL1. Cell. 2012;151(4):859–70.

    CAS  PubMed  Google Scholar 

  90. Swida-Barteczka A, Szweykowska-Kulinska Z. Micromanagement of Developmental and Stress-Induced Senescence: The Emerging Role of MicroRNAs. Genes. 2019;10(3):210.

    PubMed Central  Google Scholar 

  91. Habermann K, Tiwari B, Krantz M, Adler SO, Klipp E, Arif MA, et al. Identification of small non-coding RNAs responsive to GUN1 and GUN5 related retrograde signals in Arabidopsis thaliana. Plant J. 2020;104(1):138–55.

    CAS  PubMed  Google Scholar 

  92. Fang X, Zhao G, Zhang S, Li Y, Gu H, Li Y, et al. Chloroplast-to-Nucleus Signaling Regulates MicroRNA Biogenesis in Arabidopsis. Dev Cell. 2019;48(3):371-382.e4.

    CAS  PubMed  Google Scholar 

  93. Choudhury NK, Behera RK. Photoinhibition of Photosynthesis: Role of Carotenoids in Photoprotection of Chloroplast Constituents. Photosynthetica. 2001;39(4):481–8.

    CAS  Google Scholar 

  94. Krieger-Liszkay A, Trebst A. Tocopherol is the scavenger of singlet oxygen produced by the triplet states of chlorophyll in the PSII reaction centre. J Exp Bot. 2006;57(8):1677–84.

    CAS  PubMed  Google Scholar 

  95. Estavillo GM, Crisp PA, Pornsiriwong W, Wirtz M, Collinge D, Carrie C, et al. Evidence for a SAL1-PAP Chloroplast Retrograde Pathway That Functions in Drought and High Light Signaling in Arabidopsis. Plant Cell. 2011;23(11):3992–4012.

    CAS  PubMed  PubMed Central  Google Scholar 

  96. Gy I, Gasciolli V, Lauressergues D, Morel J-B, Gombert J, Proux F, et al. Arabidopsis FIERY1, XRN2, and XRN3 Are Endogenous RNA Silencing Suppressors. Plant Cell. 2007;19(11):3451–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  97. Li M, Yu B. Recent advances in the regulation of plant miRNA biogenesis. RNA Biol. 2021;0(0):1–10.

    Google Scholar 

  98. Price DH. Transient pausing by RNA polymerase II. Proc Natl Acad Sci. 2018;115(19):4810–2.

    CAS  PubMed  PubMed Central  Google Scholar 

  99. Conn SJ, Hocking B, Dayod M, Xu B, Athman A, Henderson S, et al. Protocol: optimising hydroponic growth systems for nutritional and physiological analysis of Arabidopsis thaliana and other plants. Plant Methods. 2013;9(1):4.

    CAS  PubMed  PubMed Central  Google Scholar 

  100. Kaufmann I, White E, Azad A, Marguerat S, Bähler J, Proudfoot NJ. Transcriptional Activation of the General Amino Acid Permease Gene per1 by the Histone Deacetylase Clr6 Is Regulated by Oca2 Kinase. Mol Cell Biol. 2010;30(13):3396–410.

    CAS  PubMed  PubMed Central  Google Scholar 

  101. Nelson J, Denisenko O, Bomsztyk K. The Fast Chromatin Immunoprecipitation Method. In: Collas P, editor. Chromatin Immunoprecipitation Assays. Methods in Molecular Biology 567:45–57. Humana Press, Totowa, NJ.

  102. Fedak H, Palusinska M, Krzyczmonik K, Brzezniak L, Yatusevich R, Pietras Z, et al. Control of seed dormancy in Arabidopsis by a cis-acting noncoding antisense transcript. Proc Natl Acad Sci U S A. 2016;113(48):E7846–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  103. Chen C-YA, Ezzeddine N, Shyu A-B. Messenger RNA Half-Life Measurements in Mammalian Cells. Methods Enzymol. 2008;448:335–57.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank C. Kim (Shanghai Center for Plant Stress Biology, Chinese Academy of Sciences) for providing flu and flu/ex1 A. thaliana seeds, N. Shao (Institute of Genetics and Developmental Biology, Chinese Academy of Sciences) for mbs1 mutant seeds, and we are also grateful to I.-C. Jang and P.J. Chung (Rockefeller University, New York) for providing the Arabidopsis line with pri-miR163 promoter::GUS::GFP construct. We also thank S. Świeżewski and A. Krieger-Liszkay for helpful discussion of the results, B. Karpińska and M. Matuszkiewicz for help in the initial phase of this work, and M. Laszczka for technical support.

Plant/ Experimental guidelines

The presented experiments did not include rare and protected plant species. The contained use of genetically modified plants was in accordance with the European and Polish regulations.


This research was supported by PRELUDIUM12 project 2016/23/N/NZ3/02237 and OPUS13 project 2017/25/B/NZ9/02574, funded by the National Science Centre Poland.

Author information

Authors and Affiliations



A.B. B formulated the hypothesis and conceived the research plan. M. F and P. G supervised the experiments. A.B. B, G. B and E. S performed the experiments. A.B. B, M. F, P. G, G. B and M. K analyzed the results. P. G and A.B. B prepared the figures. A.B. B wrote the manuscript with the contribution of all authors. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Marcin Filipecki.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Barczak-Brzyżek, A., Brzyżek, G., Koter, M. et al. Plastid retrograde regulation of miRNA expression in response to light stress. BMC Plant Biol 22, 150 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • High light
  • miRNAs
  • Chloroplast
  • Singlet oxygen
  • Plastoquinone