Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis

miR156-SPL and miR169-NF-YA Modules Regulate the Induction of Somatic Embryogenesis in Arabidopsis via LEC- and Auxin-Related Pathways

The embryogenic transition of plant somatic cells to produce somatic embryos requires extensive reprogramming of the cell transcriptome. The prominent role of transcription factors (TFs) and miRNAs in controlling somatic embryogenesis (SE) induction in plants was documented. The profiling of MIRNA expression in the embryogenic culture of Arabidopsis implied the contribution of the miR156 and miR169 to the embryogenic induction. In the present study, the function of miR156 and miR169 and the candidate targets, SPL and NF-YA genes, were investigated in Arabidopsis SE. The results showed that misexpression of MIRNA156 and candidate SPL target genes (SPL2, 3, 4, 5, 9, 10, 11, 13, 15) negatively affected the embryogenic potential of transgenic explants, suggesting that specific fine-tuning of the miR156 and target genes expression levels seems essential for efficient SE induction. The results revealed that SPL11 under the control of miR156 might contribute to SE induction by regulating the master regulators of SE, the LEC (LEAFY COTYLEDON) genes (LEC1, LEC2, FUS3). Moreover, the role of miR169 and its candidate NF-YA targets in SE induction was demonstrated. The results showed that several miR169 targets, including NF-YA1, 3, 5, 8, and 10, positively regulated SE. We found, that miR169 via NF-YA5 seems to modulate the expression of a master SE regulator LEC1/NF-YA and other auxin-related genes: YUCCA (YUC4, 10) and PIN1 in SE induction. The study provided new insights into miR156-SPL and miR169-NF-YA functions in the auxin-related and LEC-controlled regulatory network of SE.

mRNA-Seq and miRNA-Seq Analyses Provide Insights into the Mechanism of Pinellia ternata Bulbil Initiation Induced by Phytohormones

Pinellia ternata (Thunb.) Breit (abbreviated as P. ternata) is a plant with an important medicinal value whose yield is restricted by many factors, such as low reproductive efficiency and continuous cropping obstacles. As an essential breeding material for P. ternata growth and production, the bulbils have significant advantages such as a high survival rate and short breeding cycles. However, the location effect, influencing factors, and molecular mechanism of bulbil occurrence and formation have not been fully explored. In this study, exogenously applied phytohormones were used to induce in vitro petiole of P. ternata to produce bulbil structure. Transcriptome sequencing of mRNA and miRNA were performed in the induced petiole (TCp) and the induced bulbil (TCb). Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for the identification of key genes and pathways involved in bulbil development. A total of 58,019 differentially expressed genes (DEGs) were identified. The GO and KEGG analysis indicated that DEGs were mainly enriched in plant hormone signal transduction and the starch and sucrose metabolism pathway. The expression profiles of miR167a, miR171a, and miR156a during bulbil induction were verified by qRT-PCR, indicating that these three miRNAs and their target genes may be involved in the process of bulbil induction and play an important role. However, further molecular biological experiments are required to confirm the functions of the identified bulbil development-related miRNAs and targets.

Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae

Endogenous microRNAs (miRNAs) are small non-coding RNAs that perform post-transcriptional regulatory roles across diverse cellular processes, including defence responses to biotic stresses. Pseudocercospora musae, the causal agent of Sigatoka leaf spot disease in banana (Musa spp.), is an important fungal pathogen of the plant. Illumina HiSeq 2500 sequencing of small RNA libraries derived from leaf material in Musa acuminata subsp. burmannicoides, var. Calcutta 4 (resistant) after inoculation with fungal conidiospores and equivalent non-inoculated controls revealed 202 conserved miRNAs from 30 miR-families together with 24 predicted novel miRNAs. Conserved members included those from families miRNA156, miRNA166, miRNA171, miRNA396, miRNA167, miRNA172, miRNA160, miRNA164, miRNA168, miRNA159, miRNA169, miRNA393, miRNA535, miRNA482, miRNA2118, and miRNA397, all known to be involved in plant immune responses. Gene ontology (GO) analysis of gene targets indicated molecular activity terms related to defence responses that included nucleotide binding, oxidoreductase activity, and protein kinase activity. Biological process terms associated with defence included response to hormone and response to oxidative stress. DNA binding and transcription factor activity also indicated the involvement of miRNA target genes in the regulation of gene expression during defence responses. sRNA-seq expression data for miRNAs and RNAseq data for target genes were validated using stem-loop quantitative real-time PCR (qRT-PCR). For the 11 conserved miRNAs selected based on family abundance and known involvement in plant defence responses, the data revealed a frequent negative correlation of expression between miRNAs and target host genes. This examination provides novel information on miRNA-mediated host defence responses, applicable in genetic engineering for the control of Sigatoka leaf spot disease.

Construction of ceRNA Networks at Different Stages of Somatic Embryogenesis in Garlic

LncRNA (long non-coding RNA) and mRNA form a competitive endogenous RNA (ceRNA) network by competitively binding to common miRNAs. This network regulates various processes of plant growth and development at the post-transcriptional level. Somatic embryogenesis is an effective means of plant virus-free rapid propagation, germplasm conservation, and genetic improvement, which is also a typical process to study the ceRNA regulatory network during cell development. Garlic is a typical asexual reproductive vegetable. Somatic cell culture is an effective means of virus-free rapid propagation in garlic. However, the ceRNA regulatory network of somatic embryogenesis remains unclear in garlic. In order to clarify the regulatory role of the ceRNA network in garlic somatic embryogenesis, we constructed lncRNA and miRNA libraries of four important stages (explant stage: EX; callus stage: AC; embryogenic callus stage: EC; globular embryo stage: GE) in the somatic embryogenesis of garlic. It was found that 44 lncRNAs could be used as precursors of 34 miRNAs, 1511 lncRNAs were predicted to be potential targets of 144 miRNAs, and 45 lncRNAs could be used as eTMs of 29 miRNAs. By constructing a ceRNA network with miRNA as the core, 144 miRNAs may bind to 1511 lncRNAs and 12,208 mRNAs. In the DE lncRNA-DE miRNA-DE mRNA network of adjacent stages of somatic embryo development (EX-VS-CA, CA-VS-EC, EC-VS-GE), by KEGG enrichment of adjacent stage DE mRNA, plant hormone signal transduction, butyric acid metabolism, and C5-branched dibasic acid metabolism were significantly enriched during somatic embryogenesis. Since plant hormones play an important role in somatic embryogenesis, further analysis of plant hormone signal transduction pathways revealed that the auxin pathway-related ceRNA network (lncRNAs-miR393s-TIR) may play a role in the whole stage of somatic embryogenesis. Further verification by RT-qPCR revealed that the lncRNA125175-miR393h-TIR2 network plays a major role in the network and may affect the occurrence of somatic embryos by regulating the auxin signaling pathway and changing the sensitivity of cells to auxin. Our results lay the foundation for studying the role of the ceRNA network in the somatic embryogenesis of garlic.

Integrated transcriptomic and metabolic analyses provide insights into the maintenance of embryogenic potential and the biosynthesis of phenolic acids and flavonoids involving transcription factors in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis (SE) techniques have been established for micropropagation or basic research related to plant development in many conifer species. The frequent occurrence of non-embryogenic callus (NEC) during SE has impose constraints on the application of somatic embryogenesis SE in Larix kaempferi (Lamb.) Carr, but the potential regulatory mechanisms are poorly understood. In this study, integrated transcriptomic and metabolomic analyses were performed in embryogenic callus (EC) and NEC originating from a single immature zygotic embryo to better decipher the key molecular and metabolic mechanisms required for embryogenic potential maintenance. The results showed that a total of 13,842 differentially expressed genes (DEGs) were found in EC and NEC, among which many were enriched in plant hormone signal transduction, starch and sucrose metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis, and the biosynthesis of amino acids pathways. Metabolite profiling showed that 441 differentially accumulated metabolites (DAMs) were identified in EC and NEC. Both EC and NEC had vigorous primary metabolic activities, while most secondary metabolites were upregulated in NEC. Many totipotency-related transcription factor (TF) genes such as BBMs, WUSs, and LEC1 showed higher expression levels in EC compared with NEC, which may result in the higher accumulation of indole 3-acetic acid (IAA) in EC. NEC was characterized by upregulation of genes and metabolites associated with stress responses, such as DEGs involved in jasmonic acid (JA) and ethylene (ETH) biosynthesis and signal transduction pathways, and DEGs and DAMs related to phenylpropanoid and flavonoid biosynthesis. We predicted and analyzed TFs that could target several key co-expressed structural DEGs including two C4H genes, two CcoAOMT genes and three HCT genes involved in phenylpropanoid and flavonoid biosynthesis. Based on the targeted relationship and the co-expression network, two ERFs (Lk23436 and Lk458687), one MYB (Lk34626) and one C2C2-dof (Lk37167) may play an important role in regulating phenolic acid and flavonoid biosynthesis by transcriptionally regulating the expression of these structural genes. This study shows an approach involving integrated transcriptomic and metabolic analyses to obtain insights into molecular events underlying embryogenic potential maintenance and the biosynthesis mechanisms of key metabolites involving TF regulation, which provides valuable information for the improvement of SE efficiency in L. kaempferi.

Regulatory non-coding RNAs: Emerging roles during plant cell reprogramming and in vitro regeneration

Plant regeneration is a well-known capacity of plants occurring either in vivo or in vitro. This potential is the basis for plant micropropagation and genetic transformation as well as a useful system to analyse different aspects of plant development. Recent studies have proven that RNA species with no protein-coding capacity are key regulators of cellular function and essential for cell reprogramming. In this review, the current knowledge on the role of several ncRNAs in plant regeneration processes is summarized, with a focus on cell fate reprogramming. Moreover, the involvement/impact of microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and small-interfering RNAs (siRNAs) in the regulatory networks of cell dedifferentiation, proliferation and differentiation is also analysed. A deeper understanding of plant ncRNAs in somatic cell reprogramming will allow a better modulation of in vitro regeneration processes such as organogenesis and somatic embryogenesis.

Genome-wide identification of miRNAs and targets associated with cell wall biosynthesis: Differential roles of dlo-miR397a and dlo-miR408-3p during early somatic embryogenesis in longan

The dynamic alterations in cell wall (CW) biosynthesis play an essential role in physiological isolation during the plant somatic embryogenesis (SE). However, the mechanisms underlying the functions of cell wall-associated miRNAs (CW-miRNA) remain poorly understood in plant SE. Here, we have identified 36 distinct candidate miRNAs associated with CW biosynthesis from longan third-generation genome as well as miRNA transcriptome, and modified RLM-RACE validated four distinct miRNA, which specifically targeted four CW-related genes. More importantly, we found that the dlo-miR397a-antagomir significantly enhanced DlLAC7 expression and improved laccase activity. Interestingly, inhibition of dlo-miR397a increased CW lignin deposition and promoted the tightening of protodermal cell by miRNA-mimic technology during early SE. Moreover, overexpression of dlo-miR408-3p (dlo-miR408-3p-agomir) markedly decreased DlLAC12 expression. dlo-miR408-3p-agomir activated rapid cell division, thus promoting the globular embryo (GE) development, which might be due to high DNA synthesis activity in protoepidermal cells, rather than affecting lignin synthesis. The subcellular location also indicated that both DlLAC7 and DlLAC12 proteins were primarily localized in CW and regulated CW biosynthesis. Overall, our findings provided new insight on the molecular regulatory networks comprising various miRNAs associated with cell wall, and established that dlo-miR397a and dlo-miR408-3p played differential roles during early SE in longan. The findings also shed some light on the potential role of miRNA target DlLAC regulating in vivo embryonic development of plant.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Global DNA Methylation and mRNA-miRNA Variations Activated by Heat Shock Boost Early Microspore Embryogenesis in Cabbage (Brassica oleracea)

Microspore culture, a type of haploid breeding, is extensively used in the cultivation of cruciferous crops such as cabbage. Heat shock (HS) treatment is essential to improve the embryo rate during the culture process; however, its molecular role in boosting early microspore embryogenesis (ME) remains unknown. Here we combined DNA methylation levels, miRNAs, and transcriptome profiles in isolated microspores of cabbage '01-88' under HS (32 °C for 24 h) and normal temperature (25 °C for 24 h) to investigate the regulatory roles of DNA methylation and miRNA in early ME. Global methylation levels were significantly different in the two pre-treatments, and 508 differentially methylated regions (DMRs) were identified; 59.92% of DMRs were correlated with transcripts, and 39.43% of miRNA locus were associated with methylation levels. Significantly, the association analysis revealed that 31 differentially expressed genes (DEGs) were targeted by methylation and miRNA and were mainly involved in the reactive oxygen species (ROS) response and abscisic acid (ABA) signaling, indicating that HS induced DNA methylation, and miRNA might affect ME by influencing ROS and ABA. This study revealed that DNA methylation and miRNA interfered with ME by modulating key genes and pathways, which could broaden our understanding of the molecular regulation of ME induced by HS pre-treatment.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

Function and Regulation of microRNA171 in Plant Stem Cell Homeostasis and Developmental Programing

MicroRNA171 (miR171), a group of 21-nucleotide single-strand small RNAs, is one ancient and conserved microRNA family in land plants. This review focuses on the recent progress in understanding the role of miR171 in plant stem cell homeostasis and developmental patterning, and the regulation of miR171 by developmental cues and environmental signals. Specifically, miR171 regulates shoot meristem activity and phase transition through repressing the HAIRYMERISTEM (HAM) family genes. In the model species Arabidopsis, miR171 serves as a short-range mobile signal, which initiates in the epidermal layer of shoot meristems and moves downwards within a limited distance, to pattern the apical-basal polarity of gene expression and drive stem cell dynamics. miR171 levels are regulated by light and various abiotic stresses, suggesting miR171 may serve as a linkage between environmental factors and cell fate decisions. Furthermore, miR171 family members also demonstrate both conserved and lineage-specific functions in land plants, which are summarized and discussed here.

Functional analysis of the eTM-miR171-SCL6 module regulating somatic embryogenesis in Lilium pumilum DC. Fisch

Somatic embryogenesis (SE) is of great significance in Lilium bulb production, germplasm preservation and genetic improvement. miRNAs are important regulators of plant growth and development at the transcriptional level. Previous research by our group has shown that lpu-miR171 and its target gene SCARECROW-LIKE 6 (SCL6) play an important regulatory role in lily SE, and we predicted and identified that endogenous target mimics (eTMs) can regulate lpu-miR171. However, the associated mechanism and internal regulatory network are not yet clear. In the present study, lpu-miR171 was used as an entry point to explore the regulatory network between its upstream eTMs and its downstream target gene LpSCL6, as well as to identify the mechanism of this regulatory network in Lilium SE. Tobacco transient transformation confirmed that miRNA171 significantly inhibited the expression of LpSCL6. On this basis, the Lilium stable genetic transformation system was used to demonstrate that silencing lpu-miR171a and lpu-miR171b and overexpressing LpSCL6-II and LpSCL6-I promoted starch accumulation in calli and the expression of key cell cycle genes, thus providing energy to meet preconditions for SE and accelerate the formation and development of Lilium somatic embryos. LpSCL6-II and LpSCL6-I are nuclear proteins with self-activation activity in yeast cells. In addition, we confirmed in Lilium that lpu-eTM171 is the eTM of lpu-miR171 that binds lpu-miR171 to prevent cleavage of the target gene LpSCL6, thereby promoting SE. Therefore, the present study established a new mechanism whereby the eTM-miR171-SCL6 module regulates SE in Lilium pumilum DC. Fisch. and provided new insights clarifying the mechanism of SE.

Identification and profiling of conserved microRNAs in different developmental stages of crown imperial (Fritillaria imperialis L.) using high-throughput sequencing

BACKGROUND

Novel strategies for improvement of ornamental plants and their properties relay on miRNA control of differential plant gene expression modulation. Still, in response to the same abiotic stresses, some conserved miRNA families show different expression patterns in different plant species. In parallel, the use of deep sequencing technologies reveals new levels of complexity of regulatory networks in plants through identification of new miRNAs.

METHODS AND RESULTS

Fritillaria imperialis plants were collected from their natural habitats in Koohrang, Chaharmahal va Bakhtiari, Iran. Several tissues including stamen, pistil, petal, sepal, leaf, stem, bulb and fruit were collected during three developmental stages (stem elongation, flower development and seed head stages). Using RNAseq and qRT-PCR approach, this research revealed 21 conserved miRNAs, matching 15 miRNA families, in Fritilaria imperialis.

CONCLUSIONS

The expression of seven conserved miRNAs (Fim-miR156b, Fim-miR159, Fim-miR166a-5p, Fim-miR169d-5p, Fim-miR171c, Fim-miR393 and Fim-miR396e-3p) was further investigated in different tissues and three developmental stages, suggesting different roles for these miRNAs during growth and development of crown imperial. Gained knowledge from this research can open the door to find efficient ways to secure crown imperial survival, preservation and utilization and if proven useful may be applied in other plant species as well.

Embryogenic callus induction from immature zygotic embryos and genetic transformation of Larix kaempferi 3x Larix gmelinii 9

To date, there are few reports of the successful genetic transformation of larch and other conifers, mainly because it is difficult to transform and integrate exogenous genes. In this study, hybrid larch Larix kaempferi 3x Larix gmelinii 9 cones were collected on June 27, July 1, July 4, July 7 and July 16, 2017. Embryogenic callus induction was studied using a combination of different plant growth regulators and concentrations. The results showed that July 1 was the best stage; the highest induction rate was 10.83%, which cultured in BM medium (Button medium, which formula was listed in S1 Table) with 1.0 mg/L 2,4-D (2,4-dichlorophenoxyacetic acid) and 0.2 mg/L KT(kinetin). When cultured on a proliferation medium for 12 days, proliferation was the fastest, reaching 323.08%, which could also maintain the freshness and vitality. The suitable pre-culture medium for somatic embryogenesis was 1/4 BM medium containing 10 g/L inositol and 60 g/L sucrose. The combination of 45 mg/L ABA (abscisic acid) and 75 g/L PEG4000 (Polyethyene glycol 4000) could promote the number of somatic embryos, and reached the maximum, 210 140 per 1 g FW. The genetic transformation was carried out by the Agrobacterium-mediated transformation method with embryogenic callus cultured for 12 days. The results showed the optimal OD600 of the infection solution(suspension of A. tumefaciens) was 0.5, co-culture time was 2 days, and screening concentration of Hyg (hygromycin B) was 4 mg/L. In this study, the transformation rate of resistance callus was 32.1%. It provides a reference for low genetic transformation efficiency of larch at present. This study could be beneficial for the innovation and breeding of larch by genetic engineering and provides a certain basis for rapid propagation of excellent larch germplasm resources and genetic engineering breeding of larch and other conifers.

Small Non-Coding RNAs at the Crossroads of Regulatory Pathways Controlling Somatic Embryogenesis in Seed Plants

Small non-coding RNAs (sncRNAs) are molecules with important regulatory functions during development and environmental responses across all groups of terrestrial plants. In seed plants, the development of a mature embryo from the zygote follows a synchronized cell division sequence, and growth and differentiation events regulated by highly regulated gene expression. However, given the distinct features of the initial stages of embryogenesis in gymnosperms and angiosperms, it is relevant to investigate to what extent such differences emerge from differential regulation mediated by sncRNAs. Within these, the microRNAs (miRNAs) are the best characterized class, and while many miRNAs are conserved and significantly represented across angiosperms and other seed plants during embryogenesis, some miRNA families are specific to some plant lineages. Being a model to study zygotic embryogenesis and a relevant biotechnological tool, we systematized the current knowledge on the presence and characterization of miRNAs in somatic embryogenesis (SE) of seed plants, pinpointing the miRNAs that have been reported to be associated with SE in angiosperm and gymnosperm species. We start by conducting an overview of sncRNA expression profiles in the embryonic tissues of seed plants. We then highlight the miRNAs described as being involved in the different stages of the SE process, from its induction to the full maturation of the somatic embryos, adding references to zygotic embryogenesis when relevant, as a contribution towards a better understanding of miRNA-mediated regulation of SE.

Somatic embryogenesis in Larix: the state of art and perspectives

Clonal propagation of conifers using somatic embryogenesis is essential for the selection of tree species, and for the implementation of afforestation and reforestation. In combination with cryopreservation, somatic embryogenesis creates the basis for the development of economically valuable lines of clones and elite genotypes. The industrial use of such genetically verified clone lines in forestry can significantly increase forest productivity compared to any conventional methods for improving tree crops that are available. Larch is considered as one of the main conifer candidates for large-scale reforestation, not only due to the vastness of its habitat, but also due to the unique quality of its wood, rapid growth and high ecological plasticity. However, the vast majority of larch species are characterized by uneven yields and extremely low seed quality. In this regard, obtaining planting material for reforestation from larch seeds on seed plantations is not advisable, but can be successfully implemented in afforestation programs using somatic embryogenesis technologies. Research on the somatic embryogenesis of larch has been conducted for over 30 years, which allowed considerable experience in this field to be accumulated. To date, the conditions for the initiation and maintenance of embryogenic cultures, as well as for the formation and development of somatic embryos have been determined. Significant progress has been made in the study of both the factors affecting these processes and the molecular mechanisms that underlie the various stages of embryogenesis. Nevertheless, despite the successes achieved, knowledge available today on the somatic embryogenesis of representatives of the genus Larix is still not enough to develop technologies for producing valuable plant-breeding material in vitro. This review analyzes the current state of research on the problem of somatic embryogenesis of representatives of the genus Larix. Particular attention is paid to the choice of explants for somatic embryogenesis, the composition of the media for cultivation, the dependence of the potential of somatic embryogenesis on the duration of cultivation, and the genetic control of somatic embryogenesis.

Plant miR397 and its functions

MicroRNAs (miRNAs) are noncoding, small RNAs of 20-24 nucleotides (nt) and function critically at the post-transcriptional level to regulate gene expression through cleaving mRNA targets or interfering with translation of the target mRNAs. They are broadly involved in many biological processes in plants. The miR397 family in plants contains several conserved members either in 21-nt or in 22-nt that mainly target the laccase (LAC) genes functioning in lignin synthesis and are involved in the development of plants under various conditions. Recent findings showed that miR397b in Arabidopsis could also target to Casein Kinase II Subunit Beta 3 (CKB3) and mediate circadian regulation and plant flowering. This review aims to summarise recent updates on miR397 and provides the available basis for understanding the functional mechanisms of miR397 in plant growth and development regulation and in response to external adverse stimulation.

Molecular characterization of miRNA genes and their expression in Dimocarpus longan Lour

A genome-wide analysis of longan miRNA genes was conducted, and full-length pri-miRNA transcripts were cloned. Bioinformatics and expression analyses contributed to the functional characterization of longan miRNA genes. MicroRNAs are important for the post-transcriptional regulation of target genes. However, little is known about the transcription and regulation of miRNA genes in longan (Dimocarpus longan Lour.). In this study, 80 miRNA precursors (pre-miRNA) were predicted, and their secondary structure, size, conservation, and diversity were analyzed. Furthermore, the full-length cDNA sequences of 13 longan primary miRNAs (pri-miRNAs) were amplified by RLM-RACE and SMART-RACE and analyzed, which revealed that longan pri-miRNA transcripts have multiple transcription start sites (TSSs) and the downstream pre-miRNAs are polymorphic. Accordingly, the longan pri-miRNAs and protein-encoding genes may have similar transcriptional specificities. An analysis of the longan miRNA gene promoter elements indicated that the three most abundant cis-acting elements were light-responsive, stress-responsive, and hormone-responsive elements. A quantitative real-time PCR assay elucidated the potential spatial and temporal expression patterns of longan pre-miRNAs during the early stages of somatic embryogenesis (SE) and in different longan organs/tissues. This is the first report regarding the molecular characterization of miRNA genes and their expression profiles in longan. The generated data may serve as a foundation for future research aimed at clarifying the longan miRNA gene functions.

Metabolome and Transcriptome Association Analysis Reveals Regulation of Flavonoid Biosynthesis by Overexpression of LaMIR166a in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis is an ideal model process for studying early plant development. Embryonic cell lines of Larix kaempferi (Lamb.) Carr overexpressing LaMIR166a were obtained in our previous study. Here, a combination of de novo transcriptomics and extensively targeted metabolomics was used to study the transcriptional profiles and metabolic changes in wild-type and LaMIR166a-overexpressed embryonic cell lines. A total of 459 metabolites were found in the wild-type and transgenic cell lines. Compared to those in the wild-type cell lines, transcripts and metabolites were significantly altered in the LaMIR166a-overexpressed cell lines. Among differentially expressed genes (DEGs), phenylalanine and flavonoid synthesis genes were significantly enriched, and among differentially accumulated metabolites (DAMs), phenolic acids and flavonoids accumulated in particularly high amounts. Thus, the flavonoid biosynthetic pathway seems to be the most abundant pathway in response to LaMIR166a overexpression. Based on the Kyoto Encyclopedia of Genes and Genomes database, the association analysis of metabolome and transcriptome data showed that flavonoid biosynthesis and plant hormone signal transduction processes were significantly changed in miR166a-overexpression lines, suggesting that miR166 might be involved in these processes. The present study identified a number of potential metabolites associated with LaMIR166a overexpression, providing a significant foundation for a better understanding of the regulatory mechanisms underlying miR166.

miRNA expression profiling and zeatin dynamic changes in a new model system of in vivo indirect regeneration of tomato

Callus formation and adventitious shoot differentiation could be observed on the cut surface of completely decapitated tomato plants. We propose that this process can be used as a model system to investigate the mechanisms that regulate indirect regeneration of higher plants without the addition of exogenous hormones. This study analyzed the patterns of trans-zeatin and miRNA expression during in vivo regeneration of tomato. Analysis of trans-zeatin revealed that the hormone cytokinin played an important role in in vivo regeneration of tomato. Among 183 miRNAs and 1168 predicted target genes sequences identified, 93 miRNAs and 505 potential targets were selected based on differential expression levels for further characterization. Expression patterns of six miRNAs, including sly-miR166, sly-miR167, sly-miR396, sly-miR397, novel 156, and novel 128, were further validated by qRT-PCR. We speculate that sly-miR156, sly-miR160, sly-miR166, and sly-miR397 play major roles in callus formation of tomato during in vivo regeneration by regulating cytokinin, IAA, and laccase levels. Overall, our microRNA sequence and target analyses of callus formation during in vivo regeneration of tomato provide novel insights into the regulation of regeneration in higher plants.

miRNAs as key regulators via targeting the phytohormone signaling pathways during somatic embryogenesis of plants

Somatic embryogenesis is the regeneration of embryos from the somatic cell via dedifferentiation and redifferentiation without the occurrence of fertilization. A complex network of genes regulates the somatic embryogenesis process. Especially, microRNAs (miRNAs) have emerged as key regulators by affecting phytohormone biosynthesis, transport and signal transduction pathways. miRNAs are small, non-coding small RNA regulatory molecules involved in various developmental processes including somatic embryogenesis. Several types of miRNAs such as miR156, miR157, miR 159, miR 160, miR165, miR166, miR167, miR390, miR393 and miR396 have been reported to intricate in regulating somatic embryogenesis via targeting the phytohormone signaling pathways. Here we review current research progress on the miRNA-mediated regulation involved in somatic embryogenesis via regulating auxin, ethylene, abscisic acid and cytokinin signaling pathways. Further, we also discussed the possible role of other phytohormone signaling pathways such as gibberellins, jasmonates, nitric oxide, polyamines and brassinosteroids. Finally, we conclude by discussing the expression of miRNAs and their targets involved in somatic embryogenesis and possible regulatory mechanisms cross talk with phytohormones during somatic embryogenesis.

tasiR-ARFs Production and Target Regulation during In Vitro Maize Plant Regeneration

During in vitro maize plant regeneration somatic cells change their normal fate and undergo restructuring to generate pluripotent cells able to originate new plants. Auxins are essential to achieve such plasticity. Their physiological effects are mediated by auxin response factors (ARFs) that bind auxin responsive elements within gene promoters. Small trans-acting (ta)-siRNAs, originated from miR390-guided TAS3 primary transcript cleavage, target ARF3/4 class (tasiR-ARFs). Here we found that TAS3b precursor as well as derived tasiR-ARFbD5 and tasiR-ARFbD6 display significantly lower levels in non-embryogenic callus (NEC), while TAS3g, miR390 and tasiR-ARFg are more abundant in the same tissue. However, Argonaute (AGO7) and leafbladeless 1 (LBLl) required for tasiR-ARF biogenesis showed significantly higher transcript levels in EC suggesting limited tasiR-ARF biogenesis in NEC. The five maize ARFs targeted by tasiR-ARFs were also significantly enriched in EC and accompanied by higher auxin accumulation with punctuate patterns in this tissue. At hormone half-reduction and photoperiod implementation, plant regeneration initiated from EC with transient TAS3g, miR390 and tasiR-ARFg increase. Upon complete hormone depletion, TAS3b became abundant and derived tasiR-ARFs gradually increased at further regeneration stages. ZmARF transcripts targeted by tasiR-ARFs, as well as AGO7 and LBL1 showed significantly lower levels during regeneration than in EC. These results indicate a dynamic tasiR-ARF mediated regulation throughout maize in vitro plant regeneration.

Identification and Characterization of microRNAs in the Developing Seed of Linseed Flax (Linum usitatissimum L.)

Seed development plays an important role during the life cycle of plants. Linseed flax is an oil crop and the seed is a key organ for fatty acids synthesis and storage. So it is important to understand the molecular mechanism of fatty acid biosynthesis during seed development. In this study, four small RNA libraries from early seeds at 5, 10, 20 and 30 days after flowering (DAF) were constructed and used for high-throughput sequencing to identify microRNAs (miRNAs). A total of 235 miRNAs including 114 known conserved miRNAs and 121 novel miRNAs were identified. The expression patterns of these miRNAs in the four libraries were investigated by bioinformatics and quantitative real-time polymerase chain reaction (qPCR) analysis. It was found that several miRNAs, including Lus-miRNA156a was significantly correlated with seed development process. In order to confirm the actual biological function of Lus-miRNA156a, over-expression vector was constructed and transformed to Arabidopsis. The phenotypes of homozygous transgenic lines showed decreasing of oil content and most of the fatty acid content in seeds as well as late flowering time. The results provided a clue that miRNA156a participating the fatty acid biosynthesis pathway and the detailed molecular mechanism of how it regulates the pathway needs to be further investigated.

Epigenetic Regulation of Auxin-Induced Somatic Embryogenesis in Plants

Somatic embryogenesis (SE) that is induced in plant explants in response to auxin treatment is closely associated with an extensive genetic reprogramming of the cell transcriptome. The significant modulation of the gene transcription profiles during SE induction results from the epigenetic factors that fine-tune the gene expression towards embryogenic development. Among these factors, microRNA molecules (miRNAs) contribute to the post-transcriptional regulation of gene expression. In the past few years, several miRNAs that regulate the SE-involved transcription factors (TFs) have been identified, and most of them were involved in the auxin-related processes, including auxin metabolism and signaling. In addition to miRNAs, chemical modifications of DNA and chromatin, in particular the methylation of DNA and histones and histone acetylation, have been shown to shape the SE transcriptomes. In response to auxin, these epigenetic modifications regulate the chromatin structure, and hence essentially contribute to the control of gene expression during SE induction. In this paper, we describe the current state of knowledge with regard to the SE epigenome. The complex interactions within and between the epigenetic factors, the key SE TFs that have been revealed, and the relationships between the SE epigenome and auxin-related processes such as auxin perception, metabolism, and signaling are highlighted.

Dynamic changes in the expression pattern of miRNAs and associated target genes during coconut somatic embryogenesis

Genome-wide analysis of small RNAs identifies somatic embryogenesis- specific miRNAs and their targets and provides novel insights into the mechanisms governing somatic embryogenesis in coconut, a highly in vitro recalcitrant species. Coconut, a major plantation crop of the tropics is recalcitrant to in vitro culture with a very low rate of somatic embryo turnover. Clonal propagation to enhance the production of high yielding, disease-free planting material in coconut has remained a distant reality. To better understand the molecular basis of this recalcitrance and to throw light on the complex regulatory network involved in the transition of coconut somatic cells to embryogenic calli, genome-wide profiling of small RNAs from embryogenic (EC) and non-embryogenic calli (NEC) was undertaken using Illumina Hiseq 2000 platform. We have identified a total of 110 conserved miRNAs (representing 46 known miRNA families) in both types of calli. In addition, 97 novel miRNAs (48 specific to EC, 21 specific to NEC and 28 common to both the libraries) were also identified. Among the conserved miRNAs, 10 were found to be differentially expressed between NEC and EC libraries with a log2 fold change > 2 following RPM-based normalization. miR156f, miR167c, miR169a, miR319a, miR535a, and miR5179 are upregulated and miR160a, miR166a, miR171a, and miR319b are down-regulated in NEC. To confirm the differential expression pattern and their regulatory role in SE, the expression patterns of miRNAs and their putative targets were analyzed using qRT- PCR and most of the analyzed miRNA-target pairs showed inverse correlation during somatic embryogenesis. Selected targets were further validated by RNA ligase mediated rapid amplification of 5' cDNA ends (5'RLM-RACE). Our data suggest that a few conserved miRNAs and species-specific miRNAs act in concert to regulate the process of somatic embryogenesis in coconut. The results of this study provide the first overview into the regulatory landscape of somatic embryogenesis in coconut and possible strategies for fine-tuning or reprogramming to enhance somatic embryo turn over in coconut.

The miR396-GRF Regulatory Module Controls the Embryogenic Response in Arabidopsis via an Auxin-Related Pathway

In plants, microRNAs have been indicated to control various developmental processes, including somatic embryogenesis (SE), which is triggered in the in vitro cultured somatic cells of plants. Although a transcriptomic analysis has indicated that numerous MIRNAs are differentially expressed in the SE of different plants, the role of specific miRNAs in the embryogenic reprogramming of the somatic cell transcriptome is still poorly understood. In this study, we focused on performing a functional analysis of miR396 in SE given that the transcripts of MIR396 genes and the mature molecules of miR396 were found to be increased during an SE culture of Arabidopsis. In terms of miR396 in embryogenic induction, we observed the SE-associated expression pattern of MIR396b in explants of the β-glucuronidase (GUS) reporter line. In order to gain insight into the miR396-controlled mechanism that is involved in SE induction, the embryogenic response of mir396 mutants and the 35S:MIR396b overexpressor line to media with different 2,4-Dichlorophenoxyacetic acid (2,4-D) concentrations was evaluated. The results suggested that miR396 might contribute to SE induction by controlling the sensitivity of tissues to auxin treatment. Within the targets of miR396 that are associated with SE induction, we identified genes encoding the GROWTH-REGULATING FACTOR (GRF) transcription factors, including GRF1, GRF4, GRF7, GRF8, and GRF9. Moreover, the study suggested a regulatory relationship between miR396, GRF, and the PLETHORA (PLT1 and PLT2) genes during SE induction. A complex regulatory relationship within the miR396–GRF1/4/8/9–PLT1/2 module that involves the negative and positive control of GRFs and PLT (respectively) by miR396 might be assumed.

The explant developmental stage profoundly impacts small RNA-mediated regulation at the dedifferentiation step of maize somatic embryogenesis

Maize somatic embryogenesis (SE) requires the induction of embryogenic callus and establishment of proliferation before plant regeneration. The molecular mechanisms underlying callus embryogenic potential are not well understood. Here we explored the role of small RNAs (sRNAs) and the accumulation of their target transcripts in maize SE at the dedifferentiation step using VS-535 zygotic embryos collected at distinct developmental stages and displaying contrasting in vitro embryogenic potential and morphology. MicroRNAs (miRNAs), trans-acting siRNAs (tasiRNAs), heterochromatic siRNAs (hc-siRNAs) populations and their RNA targets were analyzed by high-throughput sequencing. Abundances of specific miRNAs, tasiRNAs and targets were validated by qRT-PCR. Unique accumulation patterns were found for immature embryo at 15 Days After Pollination (DAP) and for the callus induction from this explant, as compared to 23 DAP and mature embryos. miR156, miR164, miR166, tasiARFs and the 24 nt hc-siRNAs displayed the most strikingly different patterns between explants and during dedifferentiation. According to their role in auxin responses and developmental cues, we conclude that sRNA-target regulation operating within the 15 DAP immature embryo explant provides key molecular hints as to why this stage is relevant for callus induction with successful proliferation and plant regeneration.

Small RNA profiling in Pinus pinaster reveals the transcriptome of developing seeds and highlights differences between zygotic and somatic embryos

Regulation of seed development by small non-coding RNAs (sRNAs) is an important mechanism controlling a crucial phase of the life cycle of seed plants. In this work, sRNAs from seed tissues (zygotic embryos and megagametophytes) and from somatic embryos of Pinus pinaster were analysed to identify putative regulators of seed/embryo development in conifers. In total, sixteen sRNA libraries covering several developmental stages were sequenced. We show that embryos and megagametophytes express a large population of 21-nt sRNAs and that substantial amounts of 24-nt sRNAs were also detected, especially in somatic embryos. A total of 215 conserved miRNAs, one third of which are conifer-specific, and 212 high-confidence novel miRNAs were annotated. MIR159, MIR171 and MIR394 families were found in embryos, but were greatly reduced in megagametophytes. Other families, like MIR397 and MIR408, predominated in somatic embryos and megagametophytes, suggesting their expression in somatic embryos is associated with in vitro conditions. Analysis of the predicted miRNA targets suggests that miRNA functions are relevant in several processes including transporter activity at the cotyledon-forming stage, and sulfur metabolism across several developmental stages. An important resource for studying conifer embryogenesis is made available here, which may also provide insightful clues for improving clonal propagation via somatic embryogenesis.

Control of root meristem establishment in conifers

The evolution of terrestrial plant life was made possible by the establishment of a root system, which enabled plants to migrate from aquatic to terrestrial habitats. During evolution, root organization has gradually progressed from a very simple to a highly hierarchical architecture. Roots are initiated during embryogenesis and branch afterward through lateral root formation. Additionally, adventitious roots can be formed post-embryonically from aerial organs. Induction of adventitious roots (ARs) forms the basis of the vegetative propagation via cuttings in horticulture, agriculture and forestry. This method, together with somatic embryogenesis, is routinely used to clonally multiply conifers. In addition to being utilized as propagation techniques, adventitious rooting and somatic embryogenesis have emerged as versatile models to study cellular and molecular mechanisms of embryo formation and organogenesis of coniferous species. Both formation of the embryonic root and the AR primordia require the establishment of auxin gradients within cells that coordinate the developmental response. These processes also share key elements of the genetic regulatory networks that, e.g. are triggering cell fate. This minireview gives an overview of the molecular control mechanisms associated with root development in conifers, from initiation in the embryo to post-embryonic formation in cuttings.

Genome-wide identification and profiling of microRNAs in Paulownia tomentosa cambial tissues in response to seasonal changes

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, have been shown to play essential roles in the regulation of gene expression at the post-transcriptional level. Although Paulownia tomentosa is an ecologically and economically important timber species due to its rapid growth, few efforts have focused on small RNAs (sRNAs) in the cambial tissues during winter and summer transition. In the present study, we identified 33 known miRNA families and 29 novel miRNAs which include 20 putative novel miRNAs* in P. tomentosa cambial tissues during winter and summer transition. Through differential expression analysis, we showed that 15 known miRNAs and 8 novel miRNAs were preferentially abundant in certain stage of cambial tissues. Based on the P. tomentosa mRNA transcriptome database, 1667 and 78 potential targets were predicted for 29 known and 20 novel miRNAs, respectively and the predicted targets are mostly transcription factors and functional genes. The targets of these miRNAs were enriched in "metabolic process" and "transcription regulation" by using Gene Ontology enrichment analysis. In addition, KEGG pathway analyses revealed the involvement of miRNAs in starch and sucrose metabolism and plant-pathogen interaction metabolism pathways. Noticeably, qRT-PCR expression analysis demonstrated that 9 miRNAs and their targets were existed a negative correlation in P. tomentosa cambial tissues. This study is the first to examine known and novel miRNAs and their potential targets in P. tomentosa cambial tissues during winter and summer transition and identify several candidate genes potentially regulating cambial phase transition, and thus provide a framework for further understanding of miRNAs functions in the regulation of cambial phase transition and wood formation in trees.

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

Genome-wide identification and characterization of long non-coding RNAs involved in the early somatic embryogenesis in Dimocarpus longan Lour

BACKGROUND

Long non-coding RNAs (lncRNAs) are involved in variable cleavage, transcriptional interference, regulation of DNA methylation and protein modification. However, the regulation of lncRNAs in plant somatic embryos remains unclear. The longan (Dimocarpus longan) somatic embryogenesis (SE) system is a good system for research on longan embryo development.

RESULTS

In this study, 7643 lncRNAs obtained during early SE in D. longan were identified by high-throughput sequencing, among which 6005 lncRNAs were expressed. Of the expressed lncRNAs, 4790 were found in all samples and 160 were specifically expressed in embryogenic callus (EC), 154 in incomplete embryogenic compact structures (ICpECs), and 376 in globular embryos (GEs). We annotated the 6005 expressed lncRNAs, and 1404 lncRNAs belonged to 506 noncoding RNA (ncRNA) families and 4682 lncRNAs were predicted to target protein-coding genes. The target genes included 5051 cis-regulated target genes (5712 pairs) and 1605 trans-regulated target genes (3618 pairs). KEGG analysis revealed that most of the differentially expressed target genes (mRNAs) of the lncRNAs were enriched in the "plant-pathogen interaction" and "plant hormone signaling" pathways during early longan SE. Real-time quantitative PCR confirmed that 20 selected lncRNAs showed significant differences in expression and that five lncRNAs were related to auxin response factors. Compared with the FPKM expression trends, 16 lncRNA expression trends were the same in qPCR. In lncRNA-miRNA-mRNA relationship prediction, 40 lncRNAs were predicted to function as eTMs for 15 miRNAs and 7 lncRNAs were identified as potential miRNA precursors. In addition, we verified the lncRNA-miRNA-mRNA regulatory relationships by transient expression of miRNAs (miR172a, miR159a.1 and miR398a).

CONCLUSION

Analyses of lncRNAs during early longan SE showed that differentially expressed lncRNAs were involved in expression regulation at each SE stage, and may form a regulatory network with miRNAs and mRNAs. These findings provide new insights into lncRNAs and lay a foundation for future functional analysis of lncRNAs during early longan SE.

LaMIR166a-mediated auxin biosynthesis and signalling affect somatic embryogenesis in Larix leptolepis

Somatic embryogenesis (SE) involves complex molecular signalling pathways. Understanding molecular mechanism of SE in Larix leptolepis (L. leptolepis) can aid research on genetic improvement of gymnosperms. Previously, we obtained five LaMIR166a (miR166a precursor) -overexpression embryonic cell lines in the gymnosperm Larix leptolepis. The proliferation rates of pro-embryogenic masses in transgenic and wild-type lines were calculated. Overexpression of the miR166a precursor LaMIR166a led to slower proliferation. When pro-embryogenic masses were transferred to maturation medium, the relative expression of LaMIR166a and miR166a in the LaMIR166a-overexpression lines was higher than in the wild-type during SE, while LaHDZ31-34 expression levels also increased without negative control by miR166, suggesting that regulation of HD-ZIP III by miR166a exits stage-specific characteristics. The key indole-3-acetic acid (IAA) biosynthetic gene Nitrilase of L. leptolepis (LaNIT) was identified and the effects of miR166a on auxin biosynthesis and signalling genes were studied. During SE, LaNIT, Auxin response factor1 (LaARF1) and LaARF2 mRNA levels and IAA contents were markedly higher in LaMIR166a-overexpression lines, which revealed lower deformity rate of embryos, indicating endogenous IAA synthesis is required for somatic embryo maturation in L. leptolepis. Additionally, the IAA biosynthesis and signalling genes showed similar expression patterns to LaHDZ31-34, suggesting HD-ZIP III genes have a positive regulatory effect on LaNIT. Our results suggest miR166a and LaHDZ31-34 have important roles in auxin biosynthesis and signalling during SE, which might determine if the somatic embryo normally developed to mature in L. leptolepis.

miR156-SPL modules regulate induction of somatic embryogenesis in citrus callus

miR156 is a highly conserved plant miRNA and has been extensively studied because of its versatile roles in plant development. Here, we report a novel role of miR156 in regulating somatic embryogenesis (SE) in citrus, one of the most widely cultivated fruit crops in the world. SE is an important means of in vitro regeneration, but over the course of long-term sub-culturing there is always a decline in the SE potential of the preserved citrus embryogenic callus, and this represents a key obstacle for citrus biotechnology. In this study, the SE competence of citrus callus of wild kumquat (Fortunella hindsii) was significantly enhanced by either overexpression of csi-miR156a or by individual knock-down of the two target genes, CsSPL3 and CsSPL14, indicating that the effect of miR156-SPL modules was established during the initial phases of SE induction. Biological processes that might promote SE in response to miR156 overexpression were explored using RNA-seq, and mainly included hormone signaling pathways, stress responses, DNA methylation, and the cell cycle. CsAKIN10 was identified as interacting protein of CsSPL14. Our results provide insights into the regulatory pathway through which miR156-SPL modules enhance the SE potential of citrus callus, and provide a theoretical basis for improvement of plant SE competence.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21-22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

Whole Mount in situ Localization of miRNAs and mRNAs During Somatic Embryogenesis in Arabidopsis

Somatic embryogenesis (SE) results from the transition of differentiated plant somatic cells into embryogenic cells that requires the extensive reprogramming of the somatic cell transcriptome. Commonly, the SE-involved genes are identified by analyzing the heterogeneous population of explant cells and thus, it is necessary to validate the expression of the candidate genes in the cells that are competent for embryogenic transition. Here, we optimized and implemented the whole mount in situ hybridization (WISH) method (Bleckmann and Dresselhaus, 2016; Dastidar et al., 2016) in order to analyze the spatiotemporal localization of miRNAs (miR156, miR166, miR390, miR167) and mRNAs such as WOX5 and PHABULOSA-target of miR165/166 during the SE that is induced in Arabidopsis explants. This study presents a detailed step-by-step description of the WISH procedure in which DIG-labeled LNA and RNA probes were used to detect miRNAs and mRNAs, respectively. The usefulness of the WISH in the functional analysis of the SE-involved regulatory pathways is demonstrated and the advantages of this method are highlighted: (i) the ability to analyze intact non-sectioned plant tissue; (ii) the specificity of transcript detection; (iii) the detection of miRNA; and (iv) a semi-quantitative assessment of the RNA abundance.

miR160 and miR166/165 Contribute to the LEC2-Mediated Auxin Response Involved in the Somatic Embryogenesis Induction in Arabidopsis

MicroRNAs are non-coding small RNA molecules that are involved in the post-transcriptional regulation of the genes that control various developmental processes in plants, including zygotic embryogenesis (ZE). miRNAs are also believed to regulate somatic embryogenesis (SE), a counterpart of the ZE that is induced in vitro in plant somatic cells. However, the roles of specific miRNAs in the regulation of the genes involved in SE, in particular those encoding transcription factors (TFs) with an essential function during SE including LEAFY COTYLEDON2 (LEC2), remain mostly unknown. The aim of the study was to reveal the function of miR165/166 and miR160 in the LEC2-controlled pathway of SE that is induced in in vitro cultured Arabidopsis explants.In ZE, miR165/166 controls the PHABULOSA/PHAVOLUTA (PHB/PHV) genes, which are the positive regulators of LEC2, while miR160 targets the AUXIN RESPONSE FACTORS (ARF10, ARF16, ARF17) that control the auxin signaling pathway, which plays key role in LEC2-mediated SE. We found that a deregulated expression/function of miR165/166 and miR160 resulted in a significant accumulation of auxin in the cultured explants and the spontaneous formation of somatic embryos. Our results show that miR165/166 might contribute to SE induction via targeting PHB, a positive regulator of LEC2 that controls embryogenic induction via activation of auxin biosynthesis pathway (Wójcikowska et al., 2013). Similar to miR165/166, miR160 was indicated to control SE induction through auxin-related pathways and the negative impact of miR160 on ARF10/ARF16/ARF17 was shown in an embryogenic culture. Altogether, the results suggest that the miR165/166- and miR160-node contribute to the LEC2-mediated auxin-related pathway of embryogenic transition that is induced in the somatic cells of Arabidopsis. A model summarizing the suggested regulatory interactions between the miR165/166-PHB and miR160-ARF10/ARF16/ARF17 nodes that control SE induction in Arabidopsis was proposed.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21–22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21–22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

Expression profiling of AUXIN RESPONSE FACTOR genes during somatic embryogenesis induction in Arabidopsis

Extensive modulation of numerous ARF transcripts in the embryogenic culture of Arabidopsis indicates a substantial role of auxin signaling in the mechanism of somatic embryogenesis induction. Somatic embryogenesis (SE) is induced by auxin in plants and auxin signaling is considered to play a key role in the molecular mechanism that controls the embryogenic transition of plant somatic cells. Accordingly, the expression of AUXIN RESPONSE FACTOR (ARF) genes in embryogenic culture of Arabidopsis was analyzed. The study revealed that 14 of the 22 ARFs were transcribed during SE in Arabidopsis. RT-qPCR analysis indicated that the expression of six ARFs (ARF5, ARF6, ARF8, ARF10, ARF16, and ARF17) was significantly up-regulated, whereas five other genes (ARF1, ARF2, ARF3, ARF11, and ARF18) were substantially down-regulated in the SE-induced explants. The activity of ARFs during SE was also monitored with GFP reporter lines and the ARFs that were expressed in areas of the explants engaged in SE induction were detected. A functional test of ARFs transcribed during SE was performed and the embryogenic potential of the arf mutants and overexpressor lines was evaluated. ARFs with a significantly modulated expression during SE coupled with an impaired embryogenic response of the relevant mutant and/or overexpressor line, including ARF1, ARF2, ARF3, ARF5, ARF6, ARF8, and ARF11 were indicated as possibly being involved in SE induction. The study provides evidence that embryogenic induction strongly depends on ARFs, which are key regulators of the auxin signaling. Some clues on the possible functions of the candidate ARFs, especially ARF5, in the mechanism of embryogenic transition are discussed. The results provide guidelines for further research on the auxin-related functional genomics of SE and the developmental plasticity of somatic cells.

Identification of novel miRNAs and miRNA expression profiling in embryogenic tissues of Picea balfouriana treated by 6-benzylaminopurine

Here, we compared miRNA expression profiles in embryonic cell cultures of the conifer Picea balfouriana following application of the synthetic cytokinin 6-benzylaminopurine (6-BAP). We used next-generation sequencing to analyze three libraries of small RNAs from the treated embryogenic cell cultures and generated 24,000,000 raw reads from each of the libraries. Over 70 differentially regulated micro RNA (miRNA) families (≥2 fold change in expression) were identified between pairs of treatments. A quantitative analysis showed that miR3633 and miR1026 were upregulated in tissues with the highest embryogenic ability. These two miRNAs were predicted to target genes encoding receptor-like protein kinase and GAMYB transcription factors, respectively. In one library, miR1160, miR5638, miR1315, and miR5225 were downregulated. These four miRNAs were predicted to target genes encoding APETALA2, calmodulin-binding protein, and calcium-dependent protein kinase transcription factors. The expression patterns of the miRNAs and their targets were negatively correlated. Approximately 181 potentially novel P. balfouriana miRNAs were predicted from the three libraries, and seven were validated during the quantitative analysis. This study is the first report of differential miRNA regulation in tissues treated with 6-BAP during somatic embryogenesis. The differentially expressed miRNAs will be of value for investigating the mechanisms of embryogenic processes that are responsive to 6-BAP in P. balfouriana.

The pivotal role of small non-coding RNAs in the regulation of seed development

Seeds represent a crucial stage of the seed plants life cycle. It is during seed development that the foundations of the future plant body, and the ability to give rise to a new plant capable of growing under sometimes adverse environmental conditions, are established. Small non-coding RNAs are major regulators of gene expression both at the post-transcriptional and transcriptional levels and, not surprisingly, these elements play major roles in seed development and germination. We review here the current knowledge about small RNA expression and functions in seed development, going from the morphogenesis phase comprehending embryo development and patterning, to the several steps of the maturation phase, ending in the transition to the germination. A special focus is given to the small RNAs for which functional studies have been conducted and their participation in regulatory networks operating in seeds. Many challenges remain ahead for dissecting the complex small RNA landscape in seeds, but this is a highly relevant issue in plant biology and advances in this area will most certainly impact plant breeding.

Small RNA and Transcriptome Sequencing Reveal a Potential miRNA-Mediated Interaction Network That Functions during Somatic Embryogenesis in Lilium pumilum DC. Fisch

Plant somatic embryos are widely used in the fields of germplasm conservation, breeding for genetic engineering and artificial seed production. MicroRNAs (miRNAs) play pivotal roles in somatic embryogenesis (SE) regulation. However, their regulatory roles during various stages of SE remain unclear. In this study, six types of embryogenic samples of Lilium pumilum DC. Fisch., including organogenic callus, embryogenic callus induced for 4 weeks, embryogenic callus induced for 6 weeks, globular embryos, torpedo embryos and cotyledon embryos, were prepared for small RNA sequencing. The results revealed a total of 2,378,760 small RNA reads, among which the most common size was 24 nt. Four hundred and fifty-two known miRNAs, belonging to more than 86 families, 57 novel miRNAs and 40 miRNA*s were identified. The 86 known miRNA families were sorted according to an alignment with their homologs across 24 land plants into the following four categories: 23 highly conserved, 4 moderately conserved, 15 less conserved and 44 species-specific miRNAs. Differentially expressed known miRNAs were identified during various stages of SE. Subsequently, the expression levels of 12 differentially expressed miRNAs and 4 targets were validated using qRT-PCR. In addition, six samples were mixed in equal amounts for transcript sequencing, and the sequencing data were used as transcripts for miRNA target prediction. A total of 66,422 unigenes with an average length of 800 bp were assembled from 56,258,974 raw reads. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that 38,004 and 15,497 unigenes were successfully assigned to GO terms and KEGG pathways, respectively. Among the unigenes, 2,182 transcripts were predicted to be targets for 396 known miRNAs. The potential targets of the identified miRNAs were mostly classified into the following GO terms: cell, binding and metabolic process. Enriched KEGG analysis demonstrated that carbohydrate metabolism was the predominant pathway in Lilium SE. Thus, we performed systemic characterization, homology comparisons and profiling of miRNA expression, and we constructed an miRNA-target network during Lilium SE for the first time. Our findings establish a foundation for the further exploration of critical genes and elucidation of SE in Lilium.

Conservation and diversification of the miR166 family in soybean and potential roles of newly identified miR166s

BACKGROUND

microRNA166 (miR166) is a highly conserved family of miRNAs implicated in a wide range of cellular and physiological processes in plants. miR166 family generally comprises multiple miR166 members in plants, which might exhibit functional redundancy and specificity. The soybean miR166 family consists of 21 members according to the miRBase database. However, the evolutionary conservation and functional diversification of miR166 family members in soybean remain poorly understood.

RESULTS

We identified five novel miR166s in soybean by data mining approach, thus enlarging the size of miR166 family from 21 to 26 members. Phylogenetic analyses of the 26 miR166s and their precursors indicated that soybean miR166 family exhibited both evolutionary conservation and diversification, and ten pairs of miR166 precursors with high sequence identity were individually grouped into a discrete clade in the phylogenetic tree. The analysis of genomic organization and evolution of MIR166 gene family revealed that eight segmental duplications and four tandem duplications might occur during evolution of the miR166 family in soybean. The cis-elements in promoters of MIR166 family genes and their putative targets pointed to their possible contributions to the functional conservation and diversification. The targets of soybean miR166s were predicted, and the cleavage of ATHB14-LIKE transcript was experimentally validated by RACE PCR. Further, the expression patterns of the five newly identified MIR166s and 12 target genes were examined during seed development and in response to abiotic stresses, which provided important clues for dissecting their functions and isoform specificity.

CONCLUSION

This study enlarged the size of soybean miR166 family from 21 to 26 members, and the 26 soybean miR166s exhibited evolutionary conservation and diversification. These findings have laid a foundation for elucidating functional conservation and diversification of miR166 family members, especially during seed development or under abiotic stresses.

Identification and Evaluation of Suitable Reference Genes for Normalization of MicroRNA Expression in Helicoverpa armigera (Lepidoptera: Noctuidae) Using Quantitative Real-Time PCR

More and more studies have focused on microRNAs (miRNAs) expression in the pest Helicoverpa armigera (Lepidoptera: Noctuidae) recently. Quantitative real-time PCR (qRT-PCR) is being widely used in miRNA expression studies. Suitable reference genes are necessary for the correct analysis of results. In this study, 10 candidate genes of H. armigera were selected and analyzed for their expression stability under different biotic and abiotic conditions with 3 statistical methods, including geNorm, NormFinder, and Bestkeeper. Combination the best number of reference genes was calculated by geNorm. One target gene, let-7, was used to validate the selection of reference genes. The suitable candidate reference genes were shown as follows: miR-9 and U6 snRNA for developmental stages, miR-100 and U6 snRNA for larval tissues, miR-100 and miR-305 for adult tissues, miR-9 and miR-279 for parasitic treatment, miR-998 and U6 snRNA for nuclear polyhedrosis virus infection, miR-9 and U6 snRNA for insecticide treatment, miR-92a, miR-100, and miR-279 for temperature treatment, miR-92a, miR-305, and miR-998 for starvation treatment, miR-9 and miR-279 for light treatment, miR-305 and miR-998 for hormone treatment, and there was not one reference gene suitable for all samples. This study could promote future research on miRNAs expression in H. armigera with optimal reference genes under different experimental conditions.

MicroRNAs Are Intensively Regulated during Induction of Somatic Embryogenesis in Arabidopsis

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.