Original Research ARTICLE
MicroRNAs Are Intensively Regulated during Induction of Somatic Embryogenesis in Arabidopsis
- 1Department of Genetics, Faculty of Biology and Environmental Protection, University of Silesia, Katowice, Poland
- 2Department of Gene Expression, Faculty of Biology, Institute of Molecular Biology and Biotechnology, Adam Mickiewicz University, Poznan, Poland
Several genes encoding transcription factors (TFs) were indicated to have a key role in the induction of somatic embryogenesis (SE), which is triggered in the somatic cells of plants. In order to further explore the genetic regulatory network that is involved in the embryogenic transition induced in plant somatic cells, micro-RNA (miRNAs) molecules, the products of MIRNA (MIR) genes and the common regulators of TF transcripts, were analyzed in an embryogenic culture of Arabidopsis thaliana. In total, the expression of 190 genes of the 114 MIRNA families was monitored during SE induction and the levels of the primary (pri-miRNAs) transcripts vs. the mature miRNAs were investigated. The results revealed that the majority (98%) of the MIR genes were active and that most of them (64%) were differentially expressed during SE. A distinct attribute of the MIR expression in SE was the strong repression of MIR transcripts at the early stage of SE followed by their significant up-regulation in the advanced stage of SE. Comparison of the mature miRNAs vs. pri-miRNAs suggested that the extensive post-transcriptional regulation of miRNA is associated with SE induction. Candidate miRNA molecules of the assumed function in the embryogenic response were identified among the mature miRNAs that had a differential expression in SE, including miR156, miR157, miR159, miR160, miR164, miR166, miR169, miR319, miR390, miR393, miR396, and miR398. Consistent with the central role of phytohormones and stress factors in SE induction, the functions of the candidate miRNAs were annotated to phytohormone and stress responses. To confirm the functions of the candidate miRNAs in SE, the expression patterns of the mature miRNAs and their presumed targets were compared and regulatory relation during SE was indicated for most of the analyzed miRNA-target pairs. The results of the study contribute to the refinement of the miRNA-controlled regulatory pathways that operate during embryogenic induction in plants and provide a valuable platform for the identification of the genes that are targeted by the candidate miRNAs in SE induction.
Somatic embryogenesis (SE) reflects the unique developmental potential of plant somatic cells, which results in the transition of the differentiated somatic cells that are cultured in vitro into the embryogenic ones that form the somatic embryos. Thus, studies on SE provide basic knowledge about the molecular and genetic mechanisms that govern the developmental plasticity in plants. It is believed that genes that have a regulatory function activated by plant growth regulators and stress that is imposed in vitro play a key role in the mechanism of embryogenic transition (Jiménez, 2005; Karami and Saidi, 2010). In line with this assumption, numerous genes encoding transcription factors (TFs) were indicated as being involved in the regulatory pathway that operates in SE induction, including LEAFY COTYLEDON2 (LEC2) (Gaj et al., 2005; Ledwoń and Gaj, 2009; Wójcikowska et al., 2013), BABY BOOM (BBM) (Boutilier et al., 2002), WUSCHEL (WUS) (Zuo et al., 2002), and AGAMOUS-LIKE15 (AGL15) (Harding et al., 2003; Zheng et al., 2013).
In the regulation of the TF expression micro-RNA molecules (miRNAs), the products of MIRNA (MIR) genes have an essential function. miRNAs are single-stranded RNA molecules of 21–24 nucleotides that regulate the expression of the genes that are involved in plant development (for review, Bartel, 2009; Rubio-Somoza and Weigel, 2011). The biogenesis of mature miRNAs, which are the functional products of the MIR genes, is a multi-stage process that involves numerous interacting proteins. The primary MIR transcripts (pri-miRNA) are processed by DCL1 (DICER LIKE 1) RNase III, that is accompanied by the double-stranded RNA binding protein HYPONASTIC LEAVES 1 (HYL 1), the C2H2-zinc finger protein SERRATE (SE), and two cap binding proteins, CBP20 and CBP80/ABH1 (for review, Voinnet, 2009). In addition, the DDL (DAWDLE) protein was proposed to stabilize pri-miRNAs and facilitate the maturation of miRNA (Yu et al., 2008). As a result, the miRNA/miRNA* duplex that is produced in the nucleus of a plant cell is transported to the cytoplasm where the miRNA strand is bound by the protein of the ARGONAUTE (AGO) family to form the RNA-Induced Silencing Complex (RISC) engaged in the recognition of the target transcripts that are complementary to the miRNA sequence (Baumberger and Baulcombe, 2005). Then, the miRNA-loaded RISC directs the post-transcriptional silencing of the targeted mRNA via its cleavage or translation repression (Tang et al., 2003; Brodersen et al., 2008).
The transcripts that are produced by members of the MIR gene family are processed to the identical or almost identical mature miRNA molecules. Different members of the MIR gene family are expressed in a developmental and tissue-specific manner and in response to various biotic and abiotic stimuli (Zhao et al., 2007, 2011; Moldovan et al., 2010; Kruszka et al., 2014).
Similar to the widely documented involvement of miRNA molecules in plant development in vivo (Jin et al., 2013), the expression of miRNAs was reported during in vitro induced SE in several plant species including Citrus sinensis, Dimocarpus longan, Gossypium hirsutum, Larix kaempferi, Larix leptolepis, Liriodendron tulipifera×L. chinense, Manihot esculenta, and Zea mays (Zhang et al., 2012, 2014; Li et al., 2013; Lin and Lai, 2013; Yang et al., 2013; Chávez-Hernández et al., 2015; Wu et al., 2015; Lin et al., 2015a,b; Khatabi et al., 2016). Thus, the engagement of miRNAs in the embryogenic transition that is induced in vitro is assumed, although knowledge about the function of the specific miRNA in SE induction is very limited.
In Arabidopsis, which is a model plant that has greatly contributed to the present knowledge on the genetic regulation of SE (Wójcikowska and Gaj, 2016), analysis of the MIRs/miRNAs that are associated with embryogenic induction has not yet been conducted. Thus, in the present study the expression profiles of 190 MIR genes that represented 114 MIR gene families was monitored during SE induction in an embryogenic culture of Arabidopsis. The analysis of the primary MIR transcripts was followed by the identification of mature miRNAs that were differentially accumulated during the embryogenic transition. A comparison of the pri-miRNA and the cognate mature miRNA level implied that an extensive differential processing of the primary MIR transcripts precedes the production of the functional miRNA molecules that are engaged in SE induction. The identified set of candidate miRNAs provides a valuable platform for further analysis that is aimed at deciphering the miRNA-mediated regulatory network that controls the embryogenic transition in plants.
A Vast Number of MIR Genes Is Transcribed during SE Induction
Our analysis indicated that a great majority (98%) of the analyzed MIR genes were expressed in the Col-0 explants and in the derived embryogenic culture. In total, primary transcripts of 187 MIR genes were detected at different time points of the culture (Figure 1; Table S1). A significant fraction (160; 86%) of the analyzed MIRs was transcriptionally active at all of the monitored time points (0, 5, and 10 d). SE induction resulted in the activation of a relatively low number (18) of MIR genes that were not expressed in freshly isolated explants (0 d) and a similar number (12–15) of the transcripts was induced at the early and the advanced stage of SE induction. We observed that the majority of the detected MIR transcripts (185; 99%) were expressed at both of the stages of the culture that were analyzed (5 and 10 d) and only three and six of the MIR genes were transcribed exclusively at early or advanced SE, respectively.
Figure 1. Venn diagram of the MIRNA genes transcribed in the explants (0 d) and the derived embryogenic culture (5 and 10 d).
SE Induction Is Associated with the Extensive Modulation of MIR Transcripts
Profiling of the MIR genes expression revealed that the majority (120, 64%) of the transcripts displayed a significantly modulated expression (Fold change, FC ≥ 2.0) in response to SE induction (Table S2). A closer inspection of the MIR transcript expression across the SE culture indicated a similar number of differentially expressed MIR genes in the stages of SE induction that were compared: 67 (56%) vs. 70 (58%) of the MIR genes were found to be significantly modulated during the early (5–0 d) vs. the advanced (10–5 d) SE stage, respectively (Table 1). In contrast to the similar number of the modulated MIRs, the SE stages differed distinctly in the expression patterns of these genes. During early SE induction, a majority (42; 63%) of MIRs were significantly (FC ≥ 2.0) down-regulated and a large subset (30; 71%) of these genes was found to be highly repressed (FC ≥ 10.0) (Figure 2). In contrast to early SE induction, in advanced SE, the MIR genes were predominantly (63; 90%) up-regulated (FC ≥ 2.0) and almost half (28/63; 44%) of them displayed a highly stimulated transcription (FC ≥ 10.0).
Table 1. MIRNA genes differentially expressed in the early (5–0 d) and advanced (10–5 d) stages of SE induction.
Figure 2. Venn diagrams with the number of differentially expressed MIRNA genes at FC ≥ 2 (A,B) and FC ≥ 10 (C,D) in the early (5–0 d) and advanced (10–5 d) stages of SE induction.
A set of 120 MIR genes that had a significantly modulated expression was subjected to hierarchical clustering and five distinct gene expression patterns were observed (Figure 3; Table S3). The analysis showed that numerous (61) genes that had been down-regulated at early SE were up-regulated at advanced SE induction (clusters II and III). In majority of these genes (42; 69%), the expression level at 10 d was found to be similar to 0 d of the culture (cluster II), which suggests a transient modulation of these genes. The opposite expression patterns, i.e., up-regulation in early SE followed by down-regulation in advanced SE, were noticed for a small number (7) of genes (cluster IV). Hierarchical clustering analysis also indicated numerous genes (38) that were consistently up-regulated during both of the monitored SE stages (cluster I) and a limited number of genes (14) were consistently down-regulated at both of the culture stages (cluster V).
Figure 3. Hierarchical clustering of the expression profiles of 120 MIRNA genes differentially expressed in SE culture. (A) Heat map displaying the changes in the MIR expression in the early (5–0 d) and advanced (10–5 d) stages of SE induction, (B) Five patterns (I–V) of gene expression and the relevant number of MIR genes.
In conclusion, numerous MIR genes were found to be differentially expressed during SE induction, and, for the majority of these genes, a distinct down-regulation at the early culture was followed by their up-regulation in the advanced culture.
The analyzed transcripts represented 114 MIR gene families of different sizes that ranged from 1 to 14 member genes. A closer inspection of the pri-miRNAs that were produced within the gene family indicated profound differences in the expression level in SE of the gene family members. Divergent transcript profiles of the member genes were found within the majority (79%) of the analyzed MIR families. Thus, a diverse contribution of the MIR genes of the same family to the regulation of SE induction is assumed.
Mature miRNAs of SE-Modulated Accumulation Level
To evaluate the regulatory impact of MIR genes on SE induction, the accumulation of the mature miRNAs that constitute the functional products of MIR transcripts was examined during the time course of SE. The mature miRNAs that were selected for the analysis represented the MIR genes that were differentially expressed in SE. A total of 19 mature miRNA molecules, which represented 60 pri-miRNAs encoded by 14 MIR gene families (MIR156, MIR157, MIR159, MIR160, MIR164, MIR166, MIR168, MIR169, MIR172, MIR319, MIR390, MIR393, MIR396, and MIR398), were subjected to stem-loop RT-qPCR analysis.
It was found that a great majority (85%) of the mature miRNAs that were analyzed (miR156a–f, miR156g, miR156h, miR157, miR159, miR160, miR164, miR166, miR169a–c, miR169d–g, miR169h–n, miR319a–b, miR319c, miR390, miR393, miR396, and miR398) were differentially accumulated during SE induction (Table S4). Two mature miRNAs, miR168 and miR172 that are encoded by the multigene MIR families, displayed a steady expression although the level of their pri-miRNA precursors was modulated in SE.
It was observed that the expression profiles of mature miRNAs differed distinctly in the early and advanced stage of SE induction. In general, the expression profile of the mature miRNAs that were analyzed corresponded with the global transcription pattern of the MIR genes observed during SE induction. Accordingly, in early SE the majority (69%) of mature miRNAs was found to be down-regulated whereas in the advanced SE stage, most of them were up-regulated (11; 79%) (Table 2). Corresponding expression profiles of pri- and mature miRNA were found for MIR157a–c and MIR396a in early SE and for MIR169c, MIR169d, MIR169g, MIR319b, and MIR396b in advanced SE. In contrast, some inconsistencies in the expression of pri- and mature miRNA were found for MIR157d and MIR396b in early and for MIR169a, MIR169b, MIR319a, and MIR396a in advanced SE.
Table 2. Number of differentially expressed mature miRNAs in the early (5–0 d) and advanced (10–5 d) stages of SE induction.
A closer look into the expression of individual members of the MIR families indicated that only a small subset of the family transcripts seem to contribute to the production of mature miRNAs (Figure S1; Table 3). The MIR genes that had a rather dissimilar contribution to the production of the mature miRNAs were found within the majority (11; 79%) of the MIR gene families. The majority (38; 63%) of the analyzed MIR genes was transcribed at a significantly higher level than the relevant mature miRNAs, thus suggesting that numerous pri-miRNAs were not processed into the functional miRNA molecules. In support of this assumption we observed a significantly decreased level of the transcripts that encode the key enzymes of the pri-miRNA processing machinery including DCL1 and HEN1 (Figure 4; Table S5).
Table 3. Accumulation (up-regulation in red and down-regulation in gray) of the primary MIRNA transcripts (pri-miRNA) and the relevant mature miRNAs in early (5–0 d) and advanced (10–5 d) SE induction.
Figure 4. Expression levels of AGO1, DCL1, HEN, and SE genes in the early (5 d) and the advanced stage (10 d) of SE induction. A value significantly different from 0 d is indicated with an asterisk (p ≤ 0.05).
In some instances, a delayed accumulation of mature miRNA was found in respect to the expression of the relevant pri-miRNA. Accordingly, the increased expression of MIR169h–n in early SE was followed by a large accumulation of mature-miRNA in advanced SE. A gradual accumulation of these molecules in the tissues undergoing SE induction cannot be excluded when attempting to explain a delayed increase of the mature miRNA.
In conclusion, in contrast to the robustly modulated transcription of the MIR genes accumulation of the mature miRNAs seem to be rather confined during embryogenic response. Hence, the extensive post-transcriptional regulation is assumed to be associated with the production of the functional miRNAs that control the reprogramming of the somatic cells into embryogenic cells.
Functional Annotation of the Candidate miRNAs
Mature miRNA molecules with a differential accumulation in SE induction (miR156/miR157, miR159, miR160, miR164, miR166, miR169, miR319, miR390, miR393 miR396, and miR398) were selected as the candidate regulators of the embryogenic response. In support of the involvement of the candidate miRNA in the regulation of SE, all of these molecules have been reported to control hormone and/or stress responses. Moreover, among the candidate miRNAs, those that had a documented impact on the development of zygotic embryos (miR156/157, miR164, miR166, and miR169), seeds (miR159), leaves (miR159, miR164, mir319, miR390, and miR396), roots (miR160, miR169, miR390, miR393, and miR396), and flowers (miR159, miR164), as well as the control of flowering (miR156/miR157, miR159, and miR169) and vegetative phase transition (miR156/157, miR169, and miR390) were found (Table S6).
To further assess the potential pathways that are controlled by the candidate miRNAs during SE induction, the functions of their predicted targets were annotated (Table S7). A total of 59 target genes were subjected to Gene Ontology (GO) analysis and 52 of them were significantly (at p < 0.01) enriched for 235 GO terms over two main functional categories–molecular function (MF) (16) and biological processes (BP) (219) (Figure 5; Tables S8, S9). The highly enriched GO terms in the MF category are related to transcription factors (TF) (GO:0001071; 31 genes); binding (GO:0005488; 48 genes) including nucleic acids (32 genes), ions (22 genes), and hormones (3 genes), and oxidoreductase catalytic activity (GO:0016721; 2 genes) (Figure 5A). TFs targeted by the candidate miRNA were annotated to nine gene families that are referred to as SBP-box (9), NAC (6), Homeodomain-like (6), TCP (5), MYB (4), F-box (4), CCAAT-binding (3), ARF (3), and Homeobox (2).
Figure 5. Functional categories: (A) molecular function and (B) biological process of the target genes annotated to the candidate miRNAs that are assumed to be involved in SE induction. ABA, abscisic acid; AU, auxin; ET, ethylene; GA, gibberellic acid.
The annotated to the candidate miRNA GO terms infer that a wide range of biological processes is involved in regulation of SE induction (Figure 5B). The main, highly enriched BP categories are related to metabolic (GO:0008152; 50 genes) and developmental (GO:0032502; 39 genes) processes and the responses to stimuli (GO:0051716; 19 genes) (Table S9).
Among the targets of the candidate miRNAs, those that are related to metabolic processes were found to be overrepresented in the BP category. Within this category, numerous target genes that are involved in the regulation of the macromolecule metabolic process (GO:0060255), the regulation of transcription (GO:0006355) and RNA biosynthesis (GO:2001141) were identified.
Targets of the great majority (92%) of the candidate miRNAs were annotated to plant development, a general category that covers diverse developmental processes including phyllome (GO:0048827), shoot (GO:0048367), and root development (GO:0048364), the transition of the meristem from the vegetative to the reproductive phase (GO:0010228) and the regulation of the development of the reproductive structures (GO:0048608), such as flowers (GO:0009908), embryos (GO:0009793), and seeds (GO:0048316). Among the genes that are related to zygotic embryogenesis, which is a process that corresponds to somatic embryo development, the targets of miR164, miR166, miR169, miR319, and miR396 were identified.
Importantly to the mechanism of SE induction, another functional category that was highly enriched in genes was found to be related to the response to stimuli (GO:0051716) and targets of 92% of the candidate miRNA were annotated to this functional category, including numerous genes that are involved in plant responses to hormones, particularly auxin (GO:0009733), abscisic acid (GO:0009737), ethylene (GO:0009723), and gibberellic acid (GO:0071370), osmotic stress (GO:0071470) and radiation (GO:0071478). The genes that are related to hormone signaling were identified within the targets of miR159, miR160, miR164, miR319, miR393, and miR396.
Verification of the Function of the Candidate miRNAs in SE–Target Analysis
To confirm the involvement of the candidate miRNAs in SE, the expression patterns of the presumed miRNA targets were analyzed in an embryogenic culture. In total, the expression levels of 21 genes that are targeted by seven candidate miRNAs were examined (Figure 6; Table S5). An analysis of seven SPL (SPL2,3,9,10,11,12, and 13) genes, which are possible targets of miR156/miR157, revealed that all of them displayed a significant accumulation of the transcripts during SE except for SPL2. The SPLs with an up-regulated expression (SPL3,9,10,11,1, and 13) may be controlled by miR156 due to the decreased expression of this miRNA during SE. In addition, the results of the expression profiling infer a regulatory interaction of SPL2 and miR157 in SE due to the up-regulation of miR157 and down-regulation of SPL2 that was observed in embryogenic culture (Figure 6A).
Figure 6. Expression levels of mature miRNAs (solid lines) and the relevant target genes (dotted lines) at day 0, 5, and 10 of SE induction. (A) miR156/157, (B) miR164, (C) miR169, (D) miR319, (E) miR390, and (F) miR398. Lower case letters indicate a significant difference between 5 and 0 d (a); 10–0 d (b); and 10–5 d (c) (p ≤ 0.05).
The elevated accumulation of both miR164 targets, CUC1 and CUC2 genes, throughout SE induction was found to be opposite to the down-regulation of miR164, thus suggesting that an miR164-CUC regulatory node might operate in SE induction (Figure 6B). Within the miR169 targets, three genes of NF-YA family, NF-YA1, NF-YA8 and NF-YA10, displayed a significant decrease in the early stage of SE induction, which contrasted with the strongly increased level of miR169h–n. These results suggest that a biological function of miR169 in SE seems to be related to the repression of NF-YA genes (Figure 6C). In SE, a regulatory relation between miR319 and the genes of the TCP family cannot be excluded as the substantially increased level of miR319 that was observed in the advanced stage of the embryogenic culture was found to be accompanied by a significantly decreased expression of the TCP4 and TCP10 transcripts (Figure 6D).
The expression profiling indicated that among the TF genes that are regulated by the candidate miRNAs, key regulators of the auxin response also need to be considered. Accordingly, miR390 might control SE induction via contribution to production of the tasiARFs that repress ARF2, ARF3, and ARF4 transcripts. In support for this assumption a distinct down-regulation of ARF2 and ARF3 transcripts was observed in the early and advanced embryogenic culture coupled with a significant accumulation of miR390 during SE (Figure 6E). The target analysis indicated that besides the regulation of the genes that encode TFs, the candidate miRNA might also control the key enzymes that are involved in a stress response. In support of this assumption, the CSD1 gene seems to be controlled by miR398 during both stages of SE because inverse expression patterns of the target transcripts vs. mature miRNA molecules was observed (Figure 6F).
Materials and Methods
The expression profiles of 190 MIR genes belonging to 114 gene families were monitored at the level of the primary MIR transcripts using mirEX, a high throughput RT-qPCR platform (Bielewicz et al., 2012; Zielezinski et al., 2015, http://comgen.pl/mirex2/). The analysis encompassed the tissue at different stages of embryogenic culture derived from immature zygotic embryo (IZE) explants that had been induced on an auxin medium. The experimental design followed the analysis of SE-associated TF transcripts (Gliwicka et al., 2013) and included freshly isolated explants (0 d) and the embryogenic induction at the early (5 d) and advanced (10 d) stages. To identify the MIR genes that were differentially expressed during SE, the pri-miRNA transcript levels at 5 vs. 0 d, 10 vs. 0 d, and 10 vs. 5 d were analyzed. The level of mature miRNAs products of the selected MIR that displayed a significantly modulated expression was evaluated using stem-loop RT-qPCR.
The Arabidopsis thaliana (L.) Heynh. Col-0 genotype was analyzed in this study. Seeds of Col-0 were supplied by NASC (The Nottingham Arabidopsis Stock Centre). Plants were grown in soil pots in a “walk-in” type phytotron under controlled conditions (20–22°C, 16/8 h L/D photoperiod, light intensity of 100 μE/m2s).
Somatic Embryogenesis Induced In Vitro
Somatic embryogenesis was induced following the standard protocol (Gaj, 2001). IZE at the mid-late cotyledonary stage of development (10–12 days after pollination) were used as explants. The SE induction medium (E5) was based on B5 basal micro- and macro-salts (Gamborg et al., 1976) and supplemented with 20 g/L sucrose, 8 g/L Oxoid agar, and 5 μM 2,4-D. Cultures were maintained in a growth chamber under controlled conditions: 22°C, 16/8 h (light/dark), light intensity 50 mE/m2s.
RNA Isolation and cDNA Synthesis for the Analysis of MIRs
SE cultures were sampled at three time points for the isolation of total RNA: freshly isolated IZEs (0 d) and explants at the early (5 d) and the late (10 d) stage of SE induction. At each time point, three biological replicates were used for the analysis.
To analyse the transcripts of MIR genes (pri-miRNA), total RNA was isolated using an RNAqueous Kit (Ambion by Life Technologies) according to the manufacturer's instructions. The RNA concentration was measured using Nano-Drop ND-1000 (NanoDrop Technologies, Wilmington, Delaware, USA) and RNA integrity was tested on 1% agarose gel. Reverse transcriptase reactions were performed using 3 μg of RNA, Oligo(dT)20 Primer (Invitrogen by Life Technologies) and SuperScript III Reverse Transcriptase (Invitrogen by Life Technologies) according to the manufacturers' instructions.
RNA Isolation and cDNA Synthesis for Mature miRNAs Analysis
To evaluate the accumulation level of mature miRNAs, total RNA was isolated using a miRVana Kit (Ambion by Life Technologies) according to the manufacturer's instructions. The design of the oligonucleotides and stem-loop reverse transcriptase reactions were performed according to Speth and Laubinger (2014).
RNA Isolation and cDNA Synthesis for the Analysis of Target Genes
Total RNA was isolated using an RNAqueous Kit (Ambion by Life Technologies) according to the manufacturer's instructions. Reverse transcriptase reactions were performed using a RevertAid First Strand cDNA Synthesis Kit (Fermentas) according to the manufacturer's instructions.
Quantitative PCR for Profiling the MIRs
To monitor the accumulation of pri-miRs during SE, the RT-qPCR reaction was used. Gene-specific primers for 190 microRNAs from A. thaliana genes were designed as described in Szarzynska et al. (2009) and Bielewicz et al. (2012). The primer sequences are available on the miREX website (http://www.comgen.pl/mirex1/). RT-qPCR reactions were performed using a 7900HT Fast Real-Time PCR System (Applied Biosystems by Life Technologies) and PowerSYBR Green to monitor the dsDNA synthesis. The reaction mixture (10 μL) contained 5 μL of 2x PowerSYBR Green PCR Master Mix (Applied Biosystems by Life Technologies), cDNA and gene-specific primers (200 nM each). The following thermal profile was used for all of the qPCRs: 95°C for 10 min; 40 cycles of 95°C for 15 s; and 60°C for 1 min. After each RT-qPCR run, dissociation curve analyses were performed. The results were analyzed using SDS 2.2.1 software (Applied Biosystems by Life Technologies) (Szarzynska et al., 2009). Ct values for all of the MIR transcripts were normalized to the PP2AA3 (AT1G13320) and ELONGATION FACTOR 1-ALFA (EF-1) (AT1G07930) (Czechowski et al., 2005).
Quantitative PCR for Profiling Mature miRNAs
qPCR analyses were performed using a LightCycler 480 (Roche) to monitor the accumulation of mature miRNAs. The following RT-qPCR reaction conditions were used: Denaturation one repeat of 10 min at 95°C followed by 45 repeats of 10 s at 95°C, 8 s at the specific temperature for each of the primers, 12 s at 72°C and 5 s at 80°C. Denaturation for the melt curve analysis was conducted at 95°C followed by 15 s at 65°C and 95°C (0.1°C/s for fluorescence measurement). The stem-loop primers for the reverse transcription and the primers for RT-qPCR are listed in Table S10. The universal qPCR reverse primer sequence was designed according to Wu et al. (2007). Ct values were normalized to the EF-1 (AT1G07930).
Quantitative PCR for Profiling the Target Genes
qPCR analyses were performed using a LightCycler 480 (Roche) to monitor the accumulation of target gene transcripts. The following RT-qPCR reaction conditions were used: Denaturation–one repeat of 10 min at 95°C, followed by 45 repeats of 10 s at 95°C, 8 s at the temperature that is specific for the primer pairs, 12 s at 72°C and 5 s at 80°C. Denaturation for the melt curve analysis was conducted at 95°C followed by 15 s at 65°C and 95°C (0.1°C/s for the fluorescence measurement). Ct values were normalized to the EF-1 (AT1G07930). The primers used for the RT-qPCR reactions are listed in Table S11.
Gene and Mature miRNA Expression Level Analysis
The fold change (FC) of the SE-modulated MIRs, mature miRNAs and their presumed targets was calculated using the comparative 2(−ΔΔCt) method. In all of the analyzed culture tissue samples, the control genes displayed a constant expression pattern with Ct = 17 ± 1 and Ct = 18 ± 1 for PP2AA3 and EF-1, respectively. Candidate genes were identified using the thresholds of 2- and 10-fold changes. The FC were calculated as the ratio of the transcript levels at different SE time-points (SE-modulated genes). The reactions and calculations were performed in biological triplicate. P-value was calculated by comparisons of dCt values. As the number of data was limited it was not possible to conduct normality tests, thus non-parametric Analysis of Variance (Kruskal Wallis test) was applied to calculate the significant differences between the comparisons (corrected p ≤ 0.05).
Hierarchical Clustering of SE-Modulated MIRs
The expression profiles were clustered hierarchically using the average linkage method and Euclidean distances. The experiment was repeated for the number of clusters k = 2, 3, …, 15. For each k, the results were visualized and manually verified. A number of clusters, k = 5, was selected for further analysis as the one that was characterized by the highest within-cluster coherence and between-cluster separation. The clustering procedure and the visualizations were performed using MATLAB R2014a software.
Target Prediction and Functional Annotation
The miRNA targets that were predicted using the psRNATarget tool (Dai and Zhao, 2011) were functionally annotated with using the PLAZA Dicots database version 3.0 (http://bioinformatics.psb.ugent.be/plaza/versions/plaza_v3_dicots/). The significance of the over-representation was determined using the hypergeometric distribution followed by the Bonferroni method for multiple testing corrections (corrected p ≤ 0.01).
Extensive and Specific to SE-Stage Modulation of MIRNA Genes Is Associated with the Embryogenic Transition Induced in Arabidopsis
The study indicated that numerous MIR genes are active in the embryogenically induced somatic tissue of Arabidopsis and that the majority of them are significantly modulated during SE induction. Interestingly, the expression level of MIRs seems to be specific to the stage of SE induction and a substantial repression vs. stimulation of MIR genes was found to be characteristic to the early vs. advanced stage of SE induction, respectively. Like the MIR transcripts, the majority of mature miRNAs were found to be down-regulated during early SE and up-regulated in the advanced SE induction stage. Inverse transcription profile i.e., substantial transcript accumulation in the early SE followed by transcript down-regulation in the advanced SE stage, displayed TF genes expressed in embryogenic culture of Arabidopsis (Gliwicka et al., 2013). This observation was a reason to focus on TFs as miRNA targets. Thus, it seems that in SE, similar to ZE, miRNAs might contribute to the cellular differentiation during embryonic development via the regulation of the TF genes (Nodine and Bartel, 2010). In support of this assumption, we found genes encoding TFs to be over-represented among the targets that were annotated to the differentially expressed miRNAs. Thus, the intense miRNA-mediated regulation of TF genes could be employed in the mechanism of SE induction, similar to other developmental processes (Chen and Rajewsky, 2007; reviewed in Hobert, 2008). Moreover, a regulatory feedback loop between the TF and MIR genes might be expected during SE considering that binding sites for various TFs were identified within the MIR promoters (Megraw et al., 2006). In support of this assumption, the GRF1 (GROWTH RESPONSE FACTOR1) and GRF3, which are the targets of miR396 that were found to be differentially expressed in SE, were reported to repress the expression of MIR396a and MIR396b in Arabidopsis (Hewezi and Baum, 2012). Likewise, ARFs (AUXIN RESPONSE FACTOR) might control the SE-modulated expression of MIR160, MIR167, and MIR390 due to the AUXIN RESPONSE ELEMENTS (AREs) that were detected in the promoters of these genes (Gutierrez et al., 2009; Yoon et al., 2009; Marin et al., 2010).
In conclusion, the study provides comprehensive evidence that the regulatory interactions between TFs and miRNA play a pivotal role in the re-programming of somatic cells into embryogenic cells.
Intense Post-transcriptional Regulation of miRNA is Associated to SE Induction
The present comparative analysis of pri- and mature-miRNAs levels in an embryogenic culture of Arabidopsis indicated a global similarity in the expression profiles of these molecules. However, the up-regulation of the individual pri-miRNAs that were produced within a gene family did not always result in the accumulation of the functional product, i.e., mature miRNA. A growing number of reports have highlighted the importance of the post-transcriptional regulation of miRNA biogenesis in plant and animal development (Lee et al., 2008; Nogueira et al., 2009; Bielewicz et al., 2013; Barciszewska-Pacak et al., 2015). The extensive differential processing of the primary miRNA transcripts that was inferred in the present study appears to reflect the response of the cultured tissue to the stress conditions that were applied in vitro to induce SE. Various environmental stresses have been shown to trigger the differential processing of the primary miRNA transcripts, which is relevant to this assumption (Yan et al., 2012; Jia and Rock, 2013). A fundamental role of post-transcriptional regulation of miRNA expression in the responses to various abiotic stresses was recently postulated in Arabidopsis seedlings and, similar to the present results, a broad modulation of pri-miRNA was found to distinctly contrast to the rather confined response of mature miRNA (Barciszewska-Pacak et al., 2015).
Among the various mechanisms that control the biogenesis of miRNA, splicing efficiency, alternative splicing and polyA site selection have been postulated to be frequent (Bielewicz et al., 2013; Szweykowska-Kulińska et al., 2013). Alternative splicing may result in a low level of mature miRNA in spite of the up-regulation of the relevant pri-miRNA as was indicated for miR400 and miR846 under heat stress and ABA treatment, respectively (Yan et al., 2012; Jia and Rock, 2013). The discrepancy between the accumulation of pri- and mature miRNA might result from the differential susceptibility of pre-miRNA to the processing machinery and the diverse stability of mature miRNA (Ramachandran and Chen, 2008; Köster et al., 2014; Dolata et al., 2016). Recently, the increased levels of miRNA161 and miRNA173 coupled with down-regulated expression of the relevant pri-miRNAs were described in Arabidopsis seedlings subjected to salinity stress and ARGONAUTE 1 (AGO1) was proposed to be involved in the co-transcriptional regulation of MIR gene expression (Dolata et al., 2016).
Among the post-transcriptional processes that control the accumulation of mature miRNA, the level of the proteins that are involved in microRNA biogenesis has been postulated to regulate the efficiency of pri-miRNA processing into mature miRNAs (Rogers and Chen, 2013; Wang et al., 2013). In accordance with this postulate, in the present study, a lower than expected accumulation of the functional miRNA products in SE was coupled with the down-regulation of the genes encoding key enzymes in miRNA processing (AGO1, DCL1, HEN1, and SE).
The differential stability of miRNAs has mostly been demonstrated in animal cells and the active degradation of mature miRNAs has been identified as an important mechanism in miRNA homeostasis (Bail et al., 2010; Rüegger and Großhans, 2012). Thus, it cannot be ruled out that the high degree of stability of miRNA may account for the delayed accumulation (in the advance SE stage) of the mature products of the MIR169h–n genes of up-regulated transcription in the early stage of SE.
An inconsistency in the expression of pri- and mature miRNA was also observed for miR168. A steady level of mature miR168 was detected throughout SE in contrast to the differential expression of both members of the family (MIR168a and MIR168b). miR168 was found to control miRNA processing by targeting a key gene in this pathway, AGO1 (Vaucheret et al., 2004, 2006). However, in SE, other regulatory elements besides miR168 seem to control AGO1 as we observed a down-regulated transcription of the AGO1 transcripts in the culture that had a steady miR168 level. In support of this suggestion, complex regulatory loops were postulated to be involved in the control of AGO1 in order to ensure the correct function of the miRNA and siRNA pathways (Mallory and Vaucheret, 2009).
The results show that in the majority (79%) of the analyzed MIR gene families, the members displayed distinctly divergent expression profiles during SE. Thus, the diverse contribution of different members of the MIR family to the production of mature miRNAs, and hence, the regulation of SE is assumed. In support of this supposition, the functional diversification within MIR gene families was indicated in Arabidopsis. Accordingly, it was shown that MIR393a contributes to bacterial resistance (Navarro et al., 2006), MIR164a controls leaf differentiation (Koyama et al., 2010) and the expression of MIR165a/MIR166a is specific to the abaxial epidermis (Yao et al., 2009).
Stress- and Hormone-Related Functions of the Candidate miRNAs during SE Induction
Stress factors together with hormone treatments are widely accepted to play a pivotal role in the mechanism of SE induction (Jiménez, 2005; Zavattieri et al., 2010). In support of this belief, the promoters of MIR genes, including those that are differentially expressed in the SE of Arabidopsis (present study), have been found to be highly enriched in the cis regulatory elements that control stress- and hormone-responses (Megraw et al., 2006; Zhao and Li, 2013).
The miRNAs that have a differential accumulation in SE that were reported to control plant responses to stress include miR398. It can be assumed that miR398 contributes to SE induction via the activation of a stress protective reaction (Sunkar et al., 2006). In support of this postulate, we observed that the down-regulated expression of miR398 in early SE was accompanied by a significant up-regulation of the CSD1 (Cu/Zn superoxide dismutase 1) gene encoding a key enzyme that is involved in the responses to oxidative stress (Sunkar et al., 2006). The decreased accumulation of miR398 linked with an increased transcription of the CSD genes was also indicated in embryogenic cultures of other plants (Zhang et al., 2012; Lin and Lai, 2013). Another stress-related candidate for a possible regulatory role in SE, miR169, was reported to be highly produced in response to different stresses in Arabidopsis, tomato and rice (Zhao et al., 2009, 2011; Zhang et al., 2011). The present results suggest that during SE, miR169 might target the NF-YA (NF-YA1, NF-YA8, and NF-YA10) genes encoding the HAP2-type transcription factors, which are components of the CCAAT-box binding factor complex (CBF/NF-Y/HAP) (Testa et al., 2005). miR169-NFY regulatory interactions have been indicated as operating in various plant development processes including embryogenesis and seed development (Mu et al., 2013) as well as responses to stresses (Liu and Howell, 2010; Zhao et al., 2011; Luan et al., 2014). In SE, a stress-related function of the mir169-NF-YA10 regulatory interaction can be postulated due to the observations that stress factors affect the TaNF-YA10 expression in wheat and that the overexpression of this gene in Arabidopsis resulted in enhanced stress tolerance (Ma et al., 2015). The stress conditions that are inevitably associated with in vitro cultures may also account for the differential expression of miR319 that was indicated during SE as this molecule has been indicated as controlling the general stress-responses in Arabidopsis (Barciszewska-Pacak et al., 2015). It is possible that a mechanism of the miR319-mediated regulation of SE induction is related to auxin as miR319 was found to indirectly repress of the auxin response inhibitor, SHY2 (AtIAA3) (Koyama et al., 2010). miR319 might also exert its function in SE by targeting the TCP (TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR) genes encoding TFs that are involved in the organ-specific regulation of cell growth and differentiation (Palatnik et al., 2003; Crawford et al., 2004; Nag et al., 2009). The present results showed that during the advanced stage of SE induction, which is connected with somatic embryo differentiation, miR319 appears to control two of the TCP genes, TCP4 and TCP10.
A significantly modulated expression of numerous genes encoding auxin-responsive TFs was reported in Arabidopsis (Gliwicka et al., 2013; Wickramasuriya and Dunwell, 2015) and other plant species (Legrand et al., 2007; Sharma et al., 2008; Chakrabarty et al., 2010). Consistent with these observations, numerous miRNAs that are involved in auxin responses were identified among the candidates that are engaged in SE induction. Among them, miR393, which plays a key role in auxin signaling during plant development, was identified (Navarro et al., 2006; Si-Ammour et al., 2011). A recent report on SE in Arabidopsis confirmed that miR393 contributes to embryogenic transition by targeting the auxin receptors, TIR1 and AFB2, and modulating tissue sensitivity to auxin treatment (Wójcik and Gaj, 2016). Another auxin signaling-related miRNA candidate, miR160, was indicated as controlling the development of various organs in Arabidopsis, particularly zygotic embryos, by targeting ARF10, ARF16, and ARF17 (Liu et al., 2007; Liu and Chen, 2010). The down-regulation of miR160 that was observed in SE of Arabidopsis in this study was also documented in the embryogenic cultures of other plants (Zhang et al., 2012; Lin and Lai, 2013). In support of the regulatory relation between miR160 and ARFs during SE induction is the observation about the increased accumulation of ARF10, ARF16, and ARF17 transcripts in an embryogenic culture of Arabidopsis (B. Wójcikowska and MDG., submitted for publication). In addition, the involvement of the miR160-mediated regulation of ARF10 in the regeneration of shoots in a callus culture of Arabidopsis was also reported (Qiao et al., 2012).
Auxin/ARF-related functions can also be postulated for miR390, which had a significantly modulated expression in SE of Arabidopsis (present study) and of other plants (Zhang et al., 2012; Lin and Lai, 2013; Wu et al., 2015). miR390 has been documented as controlling the auxin signaling pathway by triggering the production of tasiARFs which down-regulate expression of ARF2, ARF3, and ARF4 genes (Allen et al., 2005; Williams et al., 2005). The present results imply that the miR390-TAS3-ARFs regulatory interaction seems to operate during early SE induction and miR390-mediated regulation of the ARF2 and ARF3 genes might contribute to auxin signaling involved in embryogenic transition induced in somatic cells.
Similarity of miRNA-Mediated Control in ZE and SE
The essential role of numerous miRNAs in the control of ZE was indicated in Arabidopsis and among them miR156/miR157, which are expressed in the early morphogenic stage of ZE, were identified (Nodine and Bartel, 2010; Willmann et al., 2011). In the present study, an accumulation of miR156h and miR157 was detected in the advanced stage of SE induction, which is relevant to ZE. Thus, it might be expected that, similar to ZE, these molecules control the morphogenesis of somatic embryos possibly through targeting the SPL (SQUAMOSA PROMOTER BINDING PROTEIN LIKE) genes. The present analysis of the SPLs vs. miR156/157 expression profiles during SE, which suggests the involvement of miR156-SPL3/9/10/11/12/13 and miR157-SPL2 regulatory modules in the control of somatic embryo development supports that assumption. The miRNA-mediated regulation of the SPL transcripts might be a common mechanism that operates during the formation of the somatic embryo as an inverse expression pattern of miR156/miR157 and SPL genes was also found in the SE of citrus and cotton (Wu et al., 2011; Yang et al., 2013).
The present results infer that miR164 might contribute to the SE mechanism via the regulation of the CUC1 and CUC2 genes encoding the CUP-SHAPED COTYLEDON transcription factors of the NAC family. Thus, similar to zygotic embryos, the miR164-CUC1/CUC2 regulatory module appears to control establishment of the shoot apical meristem in somatic embryos (Aida et al., 1999).
The promotion of the seed maturation programme during ZE in Arabidopsis requires the miR166-mediated repression of PHB and PHV (Tang et al., 2012). Importantly for the SE induction mechanism, the PHB and PHV were reported to positively control the LEC2, which is the master regulator of zygotic (Stone et al., 2001) and somatic (Gaj et al., 2005; Wójcikowska et al., 2013) embryogenesis. Some evidence suggests that miR165/166 and PHB/PHV are involved in the LEC2-controlled pathway of SE induction since the up-regulation of the PHB/PHV transcript was associated with efficient SE induction and that the silencing of the MIR166/165 genes resulted in impaired embryogenic response (A. M. W. and M. D. G., unpublished). Similar to Arabidopsis, an inhibited expression of miR166 was attributed to the early stages of an embryogenic culture of C. sinensis (Wu et al., 2011), which suggests a common function of miR166 in SE induction in plants.
The enrichment of the SE-related transcriptome in MIR transcripts that were indicated in the present study together with the extensive modulation of the TF genes that have been reported in embryogenic cultures of Arabidopsis and other plants confirm that a robust regulatory burst is associated with the reprogramming of plant somatic cells toward embryogenic development. The extensive modulation of MIR gene expression that is associated with the embryogenic transition appears to be distinctly controlled at the post-transcriptional level and as a result, the final level of mature miRNA, causative for SE induction, is adjusted.
The functions annotated to the SE-involved miRNA candidates reflect a general belief about the prevalent role of stress- and hormone-related responses in the genetic mechanism that governs SE induction. In addition, notable similarities in the miRNA-mediated regulatory pathways that operate in SE to the developmental processes in ZE are evident. The results of the study provide a valuable platform for further analysis that is aimed at the identification of the miRNA-controlled regulatory pathways that contribute to embryogenic induction in plant somatic cells. Further experiments are needed to verify the involvement of the candidate miRNAs and their postulated targets in the embryogenic transition.
MG, ZS, AJ, and KS conceived and designed research. KS, DB, JD, KN, AS, and AW conducted the experiments. MG and KS analyzed the data and wrote the manuscript. All the authors read and approved the manuscript.
This work was supported by a research grant from the National Science Centre in Poland (OPUS5 2013/09/B/NZ2/03233).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors wish to thank Adam Gudyś for the hierarchical clustering analyses, Joanna Szymanowska-Pułka for the valuable assistance in the statistical analyses, Barbara Wójcikowska for setting up the embryogenic cultures, Sylwia Kowalik and Katarzyna Krowicka for their assistance with the RT-qPCR analysis of the gene expression.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00018/full#supplementary-material
Adenot, X., Elmayan, T., Lauressergues, D., Boutet, S., Bouché, N., Gasciolli, V., et al. (2006). DRB4-Dependent TAS3 trans-acting siRNAs control leaf morphology through AGO7. Curr. Biol. 16, 927–932. doi: 10.1016/j.cub.2006.03.035
Aida, M., Ishida, T., and Tasaka, M. (1999). Shoot apical meristem and cotyledon formation during Arabidopsis embryogenesis: interaction among the CUP-SHAPED COTYLEDON and SHOOT MERISTEMLESS genes. Development 126, 1563–1570.
Aida, M., Ishida, T., Fukaki, H., Fujisawa, H., and Tasaka, M. (1997). Genes involved in organ separation in Arabidopsis: an analysis of the cup-shaped cotyledon mutant. Plant Cell 9, 841–857. doi: 10.1105/tpc.9.6.841
Barciszewska-Pacak, M., Milanowska, K., Knop, K., Bielewicz, D., Nuc, P., Plewka, P., et al. (2015). Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses. Front. Plant Sci. 6:410. doi: 10.3389/fpls.2015.00410
Baumberger, N., and Baulcombe, D. C. (2005). Arabidopsis ARGONAUTE1 is an RNA slicer that selectively recruits microRNAs and short interfering RNAs. Proc. Natl. Acad. Sci.U.S.A. 102, 11928–11933. doi: 10.1073/pnas.0505461102
Bian, H., Xie, Y., Guo, F., Han, N., Ma, S., Zeng, Z., et al. (2012). Distinctive expression patterns and roles of the miRNA393/TIR1 homolog module in regulating flag leaf inclination and primary and crown root growth in rice (Oryza sativa). New Phytol. 196, 149–161. doi: 10.1111/j.1469-8137.2012.04248.x
Bielewicz, D., Dolata, J., Zielezinski, A., Alaba, S., Szarzynska, B., Szczesniak, M. W., et al. (2012). mirEX: a platform for comparative exploration of plant pri-miRNA expression data. Nucleic Acids Res. 40, 191–197. doi: 10.1093/nar/gkr878
Boutilier, K., Offringa, R., Sharma, V. K., Kieft, H., Ouellet, T., Zhang, L., et al. (2002). Ectopic expression of BABY BOOM triggers a conversion from vegetative to embryonic growth. Plant Cell Online 14, 1737–1749. doi: 10.1105/tpc.001941.tissue
Brodersen, P., Sakvarelidze-Achard, L., Bruun-Rasmussen, M., Dunoyer, P., Yamamoto, Y. Y., Sieburth, L., et al. (2008). Widespread translational inhibition by plant miRNAs and siRNAs. Science 320, 1185–1190. doi: 10.1126/science.1159151
Cai, X., Ballif, J., Endo, S., Davis, E., Liang, M., Chen, D., et al. (2007). A putative CCAAT-binding transcription factor is a regulator of flowering timing in Arabidopsis. Plant Physiol. 145, 98–105. doi: 10.1104/pp.107.102079
Chakrabarty, D., Trivedi, P. K., Shri, M., Misra, P., Asif, M. H., Dubey, S., et al. (2010). Differential transcriptional expression following thidiazuron-induced callus differentiation developmental shifts in rice. Plant Biol. 12, 46–59. doi: 10.1111/j.1438-8677.2009.00213.x
Chávez-Hernández, E. C., Alejandri-Ramírez, N. D., Juárez-González, V. T., and Dinkova, T. D. (2015). Maize miRNA and target regulation in response to hormone depletion and light exposure during somatic embryogenesis. Front. Plant Sci. 6:555. doi: 10.3389/fpls.2015.00555
Chen, L., Wang, T., Zhao, M., and Zhang, W. (2012). Ethylene-responsive miRNAs in roots of Medicago truncatula identified by high-throughput sequencing at whole genome level. Plant Sci. 184, 14–19. doi: 10.1016/j.plantsci.2011.11.007
Cheng, H., Qin, L., Lee, S., Fu, X., Richards, D. E., Cao, D., et al. (2004). Gibberellin regulates Arabidopsis floral development via suppression of DELLA protein function. Development 131, 1055–1064. doi: 10.1242/dev.00992
Combier, J., Frugier, F., de Billy, F., Boualem, A., El-yahyaoui, F., Moreau, S., et al. (2006). MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev. 20, 3084–3088. doi: 10.1101/gad.402806
Crawford, B. C., Nath, U., Carpenter, R., and Coen, E. S. (2004). CINCINNATA controls both cell differentiation and growth in petal lobes and leaves of Antirrhinum. Plant Physiol. 135, 244–253. doi: 10.1104/pp.103.036368
Czechowski, T., Stitt, M., Altmann, T., and Udvardi, M. K. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization. Plant Physiol. 139, 5–17. doi: 10.1104/pp.105.063743.1
Dello Ioio, R., Nakamura, K., Moubayidin, L., Perilli, S., Taniguchi, M., Morita, M. T., et al. (2008). A genetic framework for the control of cell division and differentiation in the root meristem. Science 322, 1380–1384. doi: 10.1126/science.1164147
Dolata, J., Bajczyk, M., Bielewicz, D., Niedojadlo, K., Niedojadlo, J., Pietrykowska, H., et al. (2016). Salt stress reveals a new role for ARGONAUTE 1 in miRNA biogenesis at the transcriptional and post-transcriptional levels. Plant Physiol. 172, 297–312. doi: 10.1104/pp.16.00830
Gaj, M. D. (2001). Direct somatic embryogenesis as a rapid and efficient system for in vitro regeneration of Arabidopsis thaliana. Plant Cell Tissue Organ Cult. 64, 39–46. doi: 10.1023/A:1010679614721
Gaj, M. D., Zhang, S., Harada, J. J., and Lemaux, P. G. (2005). LEAFY COTYLEDON genes are essential for induction of somatic embryogenesis of Arabidopsis. Planta 222, 977–988. doi: 10.1007/s00425-005-0041-y
Gliwicka, M., Nowak, K., Balazadeh, S., Mueller-Roeber, B., and Gaj, M. D. (2013). Extensive modulation of the transcription factor transcriptome during somatic embryogenesis in Arabidopsis thaliana. PLoS ONE 8:e69261. doi: 10.1371/journal.pone.0069261
Guo, H.-S., Fei, J.-F., and Chua, N.-H. (2005). MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell Online 17, 1376–1386. doi: 10.1105/tpc.105.030841
Gutierrez, L., Bussell, J. D., Pacurar, D. I., Schwambach, J., Pacurar, M., and Bellini, C. (2009). Phenotypic plasticity of adventitious rooting in Arabidopsis is controlled by complex regulation of AUXIN RESPONSE FACTOR transcripts and microRNA abundance. Plant Cell 21, 3119–3132. doi: 10.1105/tpc.108.064758
Gutierrez, L., Mongelard, G., Floková, K., Pacurar, D. I., Novák, O., Staswick, P., et al. (2012). Auxin controls Arabidopsis adventitious root initiation by regulating jasmonic acid homeostasis. Plant Cell 24, 2515–2527. doi: 10.1105/tpc.112.099119
Harding, E. W., Tang, W., Nichols, K. W., Fernandez, D. E., and Perry, S. E. (2003). Expression and maintenance of embryogenic potential is enhanced through constitutive expression of AGAMOUS-Like 15. Plant Physiol. 133, 653–663. doi: 10.1104/pp.103.023499
Hunter, C., Willmann, M. R., Wu, G., Yoshikawa, M., de la Luz Gutiérrez-Nava, M., and Poethig, S. R. (2006). Trans-acting siRNA-mediated repression of ETTIN and ARF4 regulates heteroblasty in Arabidopsis. Development 133, 2973–2981. doi: 10.1242/dev.02491
Jia, F., and Rock, C. D. (2013). MIR846 and MIR842 comprise a cistronic MIRNA pair that is regulated by abscisic acid by alternative splicing in roots of Arabidopsis. Plant Mol. Biol. 81, 447–460. doi: 10.1007/s11103-013-0015-6
Kant, S., Bi, Y.-M., Zhu, T., and Rothstein, S. J. (2009). SAUR39, a small auxin-up RNA gene, acts as a negative regulator of auxin synthesis and transport in rice. Plant Physiol. 151, 691–701. doi: 10.1104/pp.109.143875
Kasschau, K. D., Xie, Z., Allen, E., Llave, C., Chapman, E. J., Krizan, K. A., et al. (2003). P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA function. Dev. Cell 4, 205–217. doi: 10.1016/S1534-5807(03)00025-X
Khatabi, B., Arikit, S., Xia, R., Winter, S., Oumar, D., Mongomake, K., et al. (2016). High-resolution identification and abundance profiling of cassava (Manihot esculenta Crantz) microRNAs. BMC Genomics 17:85. doi: 10.1186/s12864-016-2391-1
Kim, J. H., Woo, H. R., Kim, J., Lim, P. O., Lee, I. C., Choi, S. H., Nam, H. G., et al. (2009). Trifurcate feed-forward regulation of age-dependent cell death involving miR164 in Arabidopsis. Science 323, 1053–1057. doi: 10.1126/science.1166386
Köster, T., Meyer, K., Weinholdt, C., Smith, L. M., Lummer, M., Speth, C., et al. (2014). Regulation of pri-miRNA processing by the hnRNP-like protein AtGRP7 in Arabidopsis. Nucleic Acids Res. 42, 1–12. doi: 10.1093/nar/gku716
Koyama, T., Mitsuda, N., Seki, M., Shinozaki, K., and Ohme-Takagi, M. (2010). TCP transcription factors regulate the activities of ASYMMETRIC LEAVES1 and miR164, as well as the auxin response, during differentiation of leaves in Arabidopsis. Plant Cell 22, 3574–3588. doi: 10.1105/tpc.110.075598
Kruszka, K., Pacak, A., Swida-Barteczka, A., Nuc, P., Alaba, S., Wroblewska, Z., et al. (2014). Transcriptionally and post-transcriptionally regulated microRNAs in heat stress response in barley. J. Exp. Bot. 65, 6123–6135. doi: 10.1093/jxb/eru353
Kwong, R. W., Bui, A. Q., Lee, H., Kwong, L. W., Fischer, R. L., Goldberg, R. B., et al. (2003). LEAFY COTYLEDON1-LIKE defines a class of regulators essential for embryo development. Plant Cell 15, 5–18. doi: 10.1105/tpc.006973
Ledwoń, A., and Gaj, M. D. (2009). LEAFY COTYLEDON2 gene expression and auxin treatment in relation to embryogenic capacity of Arabidopsis somatic cells. Plant Cell Rep. 28, 1677–1688. doi: 10.1007/s00299-009-0767-2
Lee, E. J., Baek, M., Gusev, Y., Brackett, D. J., Nuovo, G. J., and Schmittgen, T. D. (2008). Systematic evaluation of microRNA processing patterns in tissues, cell lines, and tumors. RNA 14, 35–42. doi: 10.1261/rna.804508
Legrand, S., Hendriks, T., Hilbert, J.-L., and Quillet, M.-C. (2007). Characterization of expressed sequence tags obtained by SSH during somatic embryogenesis in Cichorium intybus L. BMC Plant Biol. 7:27. doi: 10.1186/1471-2229-7-27
Li, H., Johnson, P., Stepanova, A., Alonso, J. M., and Ecker, J. R. (2004). Convergence of signaling pathways in the control of differential cell growth in Arabidopsis. Dev. Cell 7, 193–204. doi: 10.1016/j.devcel.2004.07.002
Li, W. F., Zhang, S. G., Han, S. Y., Wu, T., Zhang, J. H., and Qi, L. W. (2013). Regulation of LaMYB33 by miR159 during maintenance of embryogenic potential and somatic embryo maturation in Larix kaempferi (Lamb.) Carr. Plant Cell Tissure Organ Cult. 113, 131–136. doi: 10.1007/s11240-012-0233-7
Li, W. X., Oono, Y., Zhu, J., He, X. J., Wu, J. M., Iida, K., et al. (2008). The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 20, 2238–2251. doi: 10.1105/tpc.108.059444
Lin, Y., and Lai, Z. (2013). Comparative analysis reveals dynamic changes in miRNAs and their targets and expression during somatic embryogenesis in longan (Dimocarpus longan Lour.). PLoS ONE 8:e60337. doi: 10.1371/journal.pone.0060337
Lin, Y., Lai, Z., Tian, Q., Lin, L., Lai, R., Yang, M., et al. (2015a). Endogenous target mimics down-regulate miR160 mediation of ARF10, -16, and -17 cleavage during somatic embryogenesis in Dimocarpus longan Lour. Mol. Breed. 35, 1–16. doi: 10.3389/fpls.2015.00956
Lin, Y., Lin, L., Lai, R., Liu, W., Chen, Y., Zhang, Z., et al. (2015b). MicroRNA390-directed TAS3 cleavage leads to the production of tasiRNA-ARF3/4 during somatic embryogenesis in Dimocarpus longan Lour. Front. Plant Sci. 6:1119. doi: 10.3389/fpls.2015.01119
Liu, J.-X., and Howell, S. H. (2010). bZIP28 and NF-Y transcription factors are activated by ER stress and assemble into a transcriptional complex to regulate stress response genes in Arabidopsis. Plant Cell 22, 782–796. doi: 10.1105/tpc.109.072173
Liu, P. P., Montgomery, T. A., Fahlgren, N., Kasschau, K. D., Nonogaki, H., and Carrington, J. C. (2007). Repression of AUXIN RESPONSE FACTOR10 by microRNA160 is critical for seed germination and post-germination stages. Plant J. 52, 133–146. doi: 10.1111/j.1365-313X.2007.03218.x
Liu, Q., and Chen, Y. Q. (2010). A new mechanism in plant engineering: the potential roles of microRNAs in molecular breeding for crop improvement. Biotechnol. Adv. 28, 301–307. doi: 10.1016/j.biotechadv.2010.01.002
Liu, Q., Zhang, Y. C., Wang, C. Y., Luo, Y. C., Huang, Q. J., Chen, S. Y., et al. (2009). Expression analysis of phytohormone-regulated microRNAs in rice, implying their regulation roles in plant hormone signaling. FEBS Lett. 583, 723–728. doi: 10.1016/j.febslet.2009.01.020
Luan, M., Xu, M., Lu, Y., Zhang, Q., Zhang, L., Zhang, C., et al. (2014). Family-wide survey of miR169s and NF-YAs and their expression profiles response to abiotic stress in maize roots. PLoS ONE 9:e91369. doi: 10.1371/journal.pone.0091369
Ma, X., Zhu, X., Li, C., Song, Y., Zhang, W., Xia, G., et al. (2015). Overexpression of wheat NF-YA10 gene regulates the salinity stress response in Arabidopsis thaliana. Plant Physiol. Biochem. 86, 34–43. doi: 10.1016/j.plaphy.2014.11.011
Mallory, A. C., Bartel, D. P., and Bartel, B. (2005). MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Development 17, 1–16. doi: 10.1105/tpc.105.031716.1
Marin, E., Jouannet, V., Herz, A., Lokerse, A. S., Weijers, D., Vaucheret, H., Maizel, A., et al. (2010). miR390, Arabidopsis TAS3 tasiRNAs, and their AUXIN RESPONSE FACTOR targets define an autoregulatory network quantitatively regulating lateral root growth. Plant Cell 22, 1104–1117. doi: 10.1105/tpc.109.072553
Millar, A. A., and Gubler, F. (2005). The Arabidopsis GAMYB-Like Genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell 17, 705–721. doi: 10.1105/tpc.104.027920
Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2012). AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 1819, 86–96. doi: 10.1016/j.bbagrm.2011.08.004
Moldovan, D., Spriggs, A., Yang, J., Pogson, B. J., Dennis, E. S., and Wilson, I. W. (2010). Hypoxia-responsive microRNAs and trans-acting small interfering RNAs in Arabidopsis. J. Exp. Bot. 61, 165–177. doi: 10.1093/jxb/erp296
Mu, J., Tan, H., Hong, S., Liang, Y., and Zuo, J. (2013). Arabidopsis transcription factor genes NF-YA1, 5, 6, and 9 play redundant roles in male gametogenesis, embryogenesis, and seed development. Mol. Plant 6, 188–201. doi: 10.1093/mp/sss061
Nag, A., King, S., and Jack, T. (2009). miR319a targeting of TCP4 is critical for petal growth and development in Arabidopsis. Proc. Natl. Acad. Sci.U.S.A. 106, 22534–22539. doi: 10.1073/pnas.0908718106
Navarro, L., Dunoyer, P., Jay, F., Arnold, B., Dharmasiri, N., Estelle, M., et al. (2006). A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science 312:436–439. doi: 10.1126/science.1126088
Nogueira, F. T., Chitwood, D. H., Madi, S., Ohtsu, K., Schnable, P. S., Scanlon, M. J., et al. (2009). Regulation of small RNA accumulation in the maize shoot apex. PLoS Genet. 5:e1000320. doi: 10.1371/journal.pgen.1000320
Ori, N., Cohen, A. R., Etzioni, A., Brand, A., Yanai, O., Shleizer, S., et al. (2007). Regulation of LANCEOLATE by miR319 is required for compound-leaf development in tomato. Nat. Genet 39, 787–791. doi: 10.1038/ng2036
Prigge, M. J., Otsuga, D., Alonso, J. M., Ecker, J. R., Drews, G. N., and Clark, S. E. (2005). Class III homeodomain-leucine zipper gene family members have overlapping, antagonistic, and distinct roles in Arabidopsis development. Plant Cell 17, 61–76. doi: 10.1105/tpc.104.026161.1
Qiao, M., Zhao, Z., Song, Y., Liu, Z., Cao, L., Yu, Y., et al. (2012). Proper regeneration from in vitro cultured Arabidopsis thaliana requires the microRNA-directed action of an auxin response factor. Plant J. 71, 14–22. doi: 10.1111/j.1365-313X.2012.04944.x
Reyes, J. L., and Chua, N. H. (2007). ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 49, 592–606. doi: 10.1111/j.1365-313X.2006.02980.x
Rodriguez, R. E., Mecchia, M. A., Debernardi, J. M., Schommer, C., Weigel, D., and Palatnik, J. F. (2010). Control of cell proliferation in Arabidopsis thaliana by microRNA miR396. Development 137, 103–112. doi: 10.1242/dev.043067
Sharma, S. K., Millam, S., Hedley, P. E., McNicol, J., and Bryan, G. J. (2008). Molecular regulation of somatic embryogenesis in potato: an auxin led perspective. Plant Mol. Biol. 68, 185–201. doi: 10.1007/s11103-008-9360-2
Si-Ammour, A., Windels, D., Arn-Bouldoires, E., Kutter, C., Ailhas, J., Meins, F., et al. (2011). miR393 and secondary siRNAs regulate expression of the TIR1/AFB2 auxin receptor clade and auxin-related development of Arabidopsis leaves. Plant Physiol. 157, 683–691. doi: 10.1104/pp.111.180083
Speth, C., and Laubinger, S. (2014). “Rapid and parallel quantification of small and large RNA species,” in Plant Circadian Networks Methods and Protocols, ed D. Staiger (New York, NY: Springer), 93–106.
Stone, S. L., Kwong, L. W., Yee, K. M., Pelletier, J., Lepiniec, L., Fischer, R. L., et al. (2001). LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development. Proc. Natl. Acad. Sci.U.S.A. 98, 11806–11811. doi: 10.1073/pnas.201413498
Sunkar, R., Kapoor, A., and Zhu, J.-K. (2006). Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell 18, 2051–2065. doi: 10.1105/tpc.106.041673
Szarzynska, B., Sobkowiak, L., Pant, B. D., Balazadeh, S., Scheible, W. R., Mueller-Roeber, B., et al. (2009). Gene structures and processing of Arabidopsis thaliana HYL1-dependent pri-miRNAs. Nucleic Acids Res. 37, 3083–3093. doi: 10.1093/nar/gkp189
Tang, X., Bian, S., Tang, M., Lu, Q., Li, S., Liu, X., et al. (2012). MicroRNA-mediated repression of the seed maturation program during vegetative development in Arabidopsis. PLoS Genet. 8:e1003091. doi: 10.1371/journal.pgen.1003091
Testa, A., Donati, G., Yan, P., Romani, F., Huang, T. H., Viganò, M. A., et al. (2005). Chromatin immunoprecipitation (ChIP) on chip experiments uncover a widespread distribution of NF-Y binding CCAAT sites outside of core promoters. J. Biol. Chem. 280, 13606–13615. doi: 10.1074/jbc.M414039200
Tsuji, H., Aya, K., Ueguchi-Tanaka, M., Shimada, Y., Nakazono, M., Watanabe, R., et al. (2006). GAMYB controls different sets of genes and is differentially regulated by microRNA in aleurone cells and anthers. Plant J. 47, 427–444. doi: 10.1111/j.1365-313X.2006.02795.x
Vaucheret, H., Mallory, A. C., and Bartel, D. P. (2006). AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol. Cell 22, 129–136. doi: 10.1016/j.molcel.2006.03.011
Vaucheret, H., Vazquez, F., Crété, P., and Bartel, D. P. (2004). The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 18, 1187–1197. doi: 10.1101/gad.1201404
Wang, L., Song, X., Gu, L., Li, X., Cao, S., Chu, C., et al. (2013). NOT2 proteins promote polymerase II-dependent transcription and interact with multiple microRNA biogenesis factors in Arabidopsis. Plant Cell 25, 715–727. doi: 10.1105/tpc.112.105882
Warpeha, K. M., Upadhyay, S., Yeh, J., Adamiak, J., Hawkins, S. I., Lapik, Y. R., et al. (2007). The GCR1, GPA1, PRN1, NF-Y signal chain mediates both blue light and abscisic acid responses in Arabidopsis. Plant Physiol. 143, 1590–1600. doi: 10.1104/pp.106.089904
Williams, L., Carles, C. C., Osmont, K. S., and Fletcher, J. C. (2005). A database analysis method identifies an endogenous trans-acting short-interfering RNA that targets the Arabidopsis ARF2, ARF3, and ARF4 genes. Proc. Natl. Acad. Sci.U.S.A. 102, 9703–9708. doi: 10.1073/pnas.0504029102
Wójcik, A. M., and Gaj, M. D. (2016). miR393 contributes to the embryogenic transition induced in vitro in Arabidopsis via the modification of the tissue sensitivity to auxin treatment. Planta 244, 231–243. doi: 10.1007/s00425-016-2505-7
Wójcikowska, B., and Gaj, M. D. (2016). “Embryogenesis in Arabidopsis,” in Somatic Embryogenesis. Fundamental Aspects and Applications, eds V. M. Loyola Vargas and N. Ochoa-Alejo (Switzeland: Springer), 185–199.
Wójcikowska, B., Jaskóła, K., Gasiorek, P., Meus, M., Nowak, K., and Gaj, M. D. (2013). LEAFY COTYLEDON2 (LEC2) promotes embryogenic induction in somatic tissues of Arabidopsis, via YUCCA-mediated auxin biosynthesis. Planta 238, 425–440. doi: 10.1007/s00425-013-1892-2
Woldemariam, M. G., Onkokesung, N., Baldwin, I. T., and Galis, I. (2012). Jasmonoyl-L-isoleucine hydrolase 1 (JIH1) regulates jasmonoyl-L-isoleucine levels and attenuates plant defenses against herbivores. Plant J. 72, 758–767. doi: 10.1111/j.1365-313X.2012.05117.x
Wu, X. M., Kou, S. J., Liu, Y. L., Fang, Y. N., Xu, Q., and Guo, W. W. (2015). Genomewide analysis of small RNAs in nonembryogenic and embryogenic tissues of citrus: microRNA- and siRNA-mediated transcript cleavage involved in somatic embryogenesis. Plant Biotechnol. J. 13, 383–394. doi: 10.1111/pbi.12317
Wu, X. M., Liu, M. Y., Ge, X. X., Xu, Q., and Guo, W. W. (2011). Stage and tissue-specific modulation of ten conserved miRNAs and their targets during somatic embryogenesis of Valencia sweet orange. Planta 233, 495–505. doi: 10.1007/s00425-010-1312-9
Xie, Q., Frugis, G., Colgan, D., and Chua, N. H. (2000). Arabidopsis NAC1 transduces auxin signal downstream of TIR1 to promote lateral root development. Genes Dev. 14, 3024–3036. doi: 10.1101/gad.852200
Xin, M., Wang, Y., Yao, Y., Xie, C., Peng, H., Ni, Z., et al. (2010). Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.). BMC Plant Biol. 10:123. doi: 10.1186/1471-2229-10-123
Yan, K., Liu, P., Wu, C. A., Yang, G.-D., Xu, R., Guo, Q.-H., et al. (2012). Stress-induced alternative splicing provides a mechanism for the regulation of microRNA processing in Arabidopsis thaliana. Mol. Cell 48, 521–531. doi: 10.1016/j.molcel.2012.08.032
Yang, X., Wang, L., Yuan, D., Lindsey, K., and Zhang, X. (2013). Small RNA and degradome sequencing reveal complex miRNA regulation during cotton somatic embryogenesis. J. Exp. Bot. 64, 1521–1536. doi: 10.1093/jxb/ert013
Yao, X., Wang, H., Li, H., Yuan, Z., Li, F., Yang, L., et al. (2009). Two types of cis-acting elements control the abaxial epidermis-specific transcription of the MIR165a and MIR166a genes. FEBS Lett. 583, 3711–3717. doi: 10.1016/j.febslet.2009.10.076
Yoon, E. K., Yang, J. H., Lim, J., Kim, S. H., Kim, S. K., and Lee, W. S. (2009). Auxin regulation of the microRNA390-dependent transacting small interfering RNA pathway in Arabidopsis lateral root development. Nucleic Acids Res. 38, 1382–1391. doi: 10.1093/nar/gkp1128
Yu, B., Bi, L., Zheng, B., Ji, L., Chevalier, D., Agarwal, M., et al. (2008). The FHA domain proteins DAWDLE in Arabidopsis and SNIP1 in humans act in small RNA biogenesis. Proc. Natl. Acad. Sci.U.S.A. 105, 10073–10078. doi: 10.1073/pnas.0804218105
Yu, N., Cai, W.-J., Wang, S., Shan, C.-M., Wang, L.-J., and Chen, X.-Y. (2010). Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana. Plant Cell 22, 2322–2335. doi: 10.1105/tpc.109.072579
Zavattieri, M. A., Frederico, A. M., Lima, M., Sabino, R., and Arnholdt-Schmitt, B. (2010). Induction of somatic embryogenesis as an example of stress-related plant reactions. EJB 13, 1–9. doi: 10.2225/vol13-issue1-fulltext-4
Zhang, J., Zhang, S., Han, S., Wu, T., Li, X., Li, W., et al. (2012). Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis. Planta 236, 647–657. doi: 10.1007/s00425-012-1643-9
Zhang, L.-F., Li, W.-F., Xu, H.-Y., Qi, L.-W., and Han, S.-Y. (2014). Cloning and characterization of four differentially expressed cDNAs encoding NFYA homologs involved in responses to ABA during somatic embryogenesis in Japanese larch (Larix leptolepis). Plant Cell Tissue Organ Cult. 117, 293–304. doi: 10.1007/s11240-014-0440-5
Zhang, X., Zou, Z., Gong, P., Zhang, J., Ziaf, K., Li, H., et al. (2011). Over-expression of microRNA169 confers enhanced drought tolerance to tomato. Biotechnol. Lett. 33, 403–409. doi: 10.1007/s10529-010-0436-0
Zhao, B., Ge, L., Liang, R., Li, W., Ruan, K., Lin, H., et al. (2009). Members of miR-169 family are induced by high salinity and transiently inhibit the NF-YA transcription factor. BMC Mol. Biol. 10:29. doi: 10.1186/1471-2199-10-29
Zhao, B., Liang, R., Ge, L., Li, W., Xiao, H., Lin, H., et al. (2007). Identification of drought-induced microRNAs in rice. Biochem. Biophys. Res. Commun. 354, 585–590. doi: 10.1016/j.bbrc.2007.01.022
Zhao, M., Ding, H., Zhu, J. K., Zhang, F., and Li, W. X. (2011). Involvement of miR169 in the nitrogen-starvation responses in Arabidopsis. New Phytol. 190, 906–915. doi: 10.1111/j.1469-8137.2011.03647.x.
Zheng, Q., Zheng, Y., and Perry, S. E. (2013). AGAMOUS-Like15 promotes somatic embryogenesis in Arabidopsis and soybean in part by the control of ethylene biosynthesis and response. Plant Physiol. 161, 2113–2127. doi: 10.1104/pp.113.216275
Zhou, G. K., Kubo, M., Zhong, R., Demura, T., and Ye, Z. H. (2007). Overexpression of miR165 affects apical meristem formation, organ polarity establishment and vascular development in Arabidopsis. Plant Cell Physiol. 48, 391–404. doi: 10.1093/pcp/pcm00
Zielezinski, A., Dolata, J., Alaba, S., Kruszka, K., Pacak, A., Swida-Barteczka, A., et al. (2015). mirEX 2.0 - an integrated environment for expression profiling of plant microRNAs. BMC Plant Biol. 15:144. doi: 10.1186/s12870-015-0533-2
Keywords: Arabidopsis, gene expression, mature miRNA, MIRNA genes, pri-miRNA, somatic embryogenesis
Citation: Szyrajew K, Bielewicz D, Dolata J, Wójcik AM, Nowak K, Szczygieł-Sommer A, Szweykowska-Kulinska Z, Jarmolowski A and Gaj MD (2017) MicroRNAs Are Intensively Regulated during Induction of Somatic Embryogenesis in Arabidopsis. Front. Plant Sci. 8:18. doi: 10.3389/fpls.2017.00018
Received: 12 August 2016; Accepted: 04 January 2017;
Published: 23 January 2017.
Edited by:Keqiang Wu, National Taiwan University, Taiwan
Reviewed by:Marcela Dotto, Universidad Nacional del Litoral and CEFOBI (CONICET), Argentina
Jia-Wei Wang, CAS Center for Excellence in Molecular Plant Sciences, China
Copyright © 2017 Szyrajew, Bielewicz, Dolata, Wójcik, Nowak, Szczygieł-Sommer, Szweykowska-Kulinska, Jarmolowski and Gaj. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Małgorzata D. Gaj, firstname.lastname@example.org