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

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.

ldi-miR396-LdPMaT1 enhances reactive oxygen species scavenging capacity and promotes drought tolerance in Lilium distichum Nakai autotetraploids

Drought is a key environmental stress that inhibits plant growth, development, yield and quality. Whole-genome replication is an effective method for breeding drought resistant cultivars. Here, we evaluated the tolerance of Lilium distichum Nakai diploids (2n = 2× = 24) and artificially induced autotetraploids (2n = 4× = 48) to drought simulated by polyethylene glycol (PEG) stress. Autotetraploids showed stronger drought tolerance than diploids, and high-throughput sequencing during PEG stress identified five differentially expressed miRNAs. Transcriptome analysis revealed significantly different reactive oxygen species (ROS)-scavenger expression levels between diploids and autotetraploids, which increased the drought tolerance of autotetraploids. Specifically, we identified ldi-miR396b and its only target gene (LdPMaT1) for further study based on its expression level and ROS-scavenging ability in response to drought stress (DS). Autotetraploids showed higher expression of LdPMaT1 and significantly downregulated expression of ldi-miR396b under DS compared with diploids. Through a short tandem target mimic (STTM) in transgenic lilies, functional studies revealed that miR396b silencing promotes LdPMaT1 expression and the DS response. Under PEG stress, STTM393 transgenic lines showed improved drought resistance mediated by lowered MDA content but exhibited high antioxidant enzyme activity, consistent with the autotetraploid results. Collectively, these findings suggest that ldi-miR396b-LdPMaT1 potentially enhances ROS-scavenging ability, which contributes to improved stress adaptation in autotetraploid lilies.

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

Identification and functional analysis of LpNAC37 associated with somatic embryogenesis in Lilium pumilum DC. Fisch. based on transcriptome analysis

Somatic embryogenesis (SE) is important for Lilium bulb propagation, germplasm conservation, and genetic transformation. The transition of somatic cells to embryonic cells is a critical step in SE, but the associated regulatory mechanisms have not been fully elucidated. Lilium pumilum DC. Fisch has a high regenerative capacity, and this study clarifies the critical timing of embryonic cell appearance in Lilium SE. Transcriptome sequencing using RNA-seq technology was performed on 5 representative samples from the early stage of Lilium SE. The 15 established cDNA libraries yielded 91.47 GB of valid data, and a total of 11,155 genes were consistently differentially expressed in the early stages of Lilium SE. GO annotation and KEGG pathway analysis of differentially expressed genes (DEGs) suggested that transcriptional regulation, hormone signaling, and stress response pathways play essential roles in the early stages of Lilium SE. WOX8, WOX11, SHR2, NAC37, AHP2, ANT, PIN1C, LAX2, LBD4, ACS12, YUC4, NFYB3, WRKY28, SAUR50, PYL9, and WRKY39 may be candidate genes for regulating early SE in Lilium. We further cloned LpNAC37, one of the key DEGs obtained from WGCNA and screening. LpNAC37 encodes a protein of 303 amino acids with a conserved NAM structural domain. The protein is a nuclear transcription factor with the highest homology to carrot DcNAC37. Overexpression of LpNAC37 suggested that LpNAC37 promotes embryonic callus formation in Arabidopsis. These results will help reveal the molecular mechanisms of the early stages of Lilium SE and advance the application of SE in Lilium propagation and genetic transformation.

Regulatory roles of noncoding RNAs in callus induction and plant cell dedifferentiation

Plant regulatory noncoding RNAs (ncRNAs) have emerged as key modulators of gene expression during callus induction. Their further study may promote the design of innovative plant tissue culture protocols. The use of plants by humans has recently taken on a new and expanding insight due to the advent of genetic engineering technologies. In this context, callus cultures have shown remarkable potential for synthesizing valuable biomolecules, crop improvement, plant micropropagation, and biodiversity preservation. A crucial stage in callus production is the conversion of somatic cells into totipotent cells; compelling evidence indicates that stress factors, transcriptional regulators, and plant hormones can trigger this biological event. Besides, posttranscriptional regulators of gene expression might be essential participants in callus induction. However, research related to the analysis of noncoding RNAs (ncRNAs) that modulate callogenesis and plant cell dedifferentiation in vitro is still at an early stage. During the last decade, some relevant studies have enlightened the fact that different classes of ncRNAs, such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and long noncoding RNAs (lncRNAs) are implicated in plant cell dedifferentiation through regulating the expression levels of diverse gene targets. Hence, understanding the molecular relevance of these ncRNAs in the aforesaid biological processes might represent a promising source of new biotechnological approaches for callus culture and plant improvement. In this current work, we review the experimental evidence regarding the prospective roles of ncRNAs in callus induction and plant cell dedifferentiation to promote this field of study.

Uncovering a Phenomenon of Active Hormone Transcriptional Regulation during Early Somatic Embryogenesis in Medicago sativa

Somatic embryogenesis (SE) is a developmental process in which somatic cells undergo dedifferentiation to become plant stem cells, and redifferentiation to become a whole embryo. SE is a prerequisite for molecular breeding and is an excellent platform to study cell development in the majority of plant species. However, the molecular mechanism involved in M. sativa somatic embryonic induction, embryonic and maturation is unclear. This study was designed to examine the differentially expressed genes (DEGs) and miRNA roles during somatic embryonic induction, embryonic and maturation. The cut cotyledon (ICE), non-embryogenic callus (NEC), embryogenic callus (EC) and cotyledon embryo (CE) were selected for transcriptome and small RNA sequencing. The results showed that 17,251 DEGs, and 177 known and 110 novel miRNAs families were involved in embryonic induction (ICE to NEC), embryonic (NEC to EC), and maturation (EC to CE). Expression patterns and functional classification analysis showed several novel genes and miRNAs involved in SE. Moreover, embryonic induction is an active process of molecular regulation, and hormonal signal transduction related to pathways involved in the whole SE. Finally, a miRNA–target interaction network was proposed during M. sativa SE. This study provides novel perspectives to comprehend the molecular mechanisms in M. sativa SE.

Integrated Analysis of Physiological, mRNA Sequencing, and miRNA Sequencing Data Reveals a Specific Mechanism for the Response to Continuous Cropping Obstacles in Pogostemon cablin Roots

Pogostemon cablin (patchouli) is a commercially important medicinal and industrial crop grown worldwide for its medicinal and aromatic properties. Patchoulol and pogostone, derived from the essential oil of patchouli, are considered valuable components in the cosmetic and pharmaceutical industries. Due to its high application value in the clinic and industry, the demand for patchouli is constantly growing. Unfortunately, patchouli cultivation has suffered due to severe continuous cropping obstacles, resulting in a significant decline in yield and quality. Moreover, the physiological and transcriptional changes in patchouli in response to continuous cropping obstacles remain unclear. This has greatly restricted the development of the patchouli industry. To explore the mechanism underlying the rapid response of patchouli roots to continuous cropping stress, integrated analysis of the transcriptome and miRNA profiles of patchouli roots under continuous and noncontinuous cropping conditions in different growth periods was conducted using RNA sequencing (RNA-seq) and miRNA-seq and complemented with physiological data. The physiological and biochemical results showed that continuous cropping significantly inhibited root growth, decreased root activity, and increased the activity of antioxidant enzymes (superoxide dismutase, peroxidase, and catalase) and the levels of osmoregulators (malondialdehyde, soluble protein, soluble sugar, and proline). Subsequently, we found 4,238, 3,494, and 7,290 upregulated and 4,176, 3,202, and 8,599 downregulated differentially expressed genes (DEGs) in the three growth periods of continuously cropped patchouli, many of which were associated with primary carbon and nitrogen metabolism, defense responses, secondary metabolite biosynthesis, and transcription factors. Based on miRNA-seq, 927 known miRNAs and 130 novel miRNAs were identified, among which 67 differentially expressed miRNAs (DEMIs) belonging to 24 miRNA families were induced or repressed by continuous cropping. By combining transcriptome and miRNA profiling, we obtained 47 miRNA-target gene pairs, consisting of 18 DEMIs and 43 DEGs, that likely play important roles in the continuous cropping response of patchouli. The information provided in this study will contribute to clarifying the intricate mechanism underlying the patchouli response to continuous cropping obstacles. In addition, the candidate miRNAs and genes can provide a new strategy for breeding continuous cropping-tolerant patchouli.

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.

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.

Selection of reference genes for gene expression analysis in Liriodendron hybrids' somatic embryogenesis and germinative tissues

The differential expression of genes is crucial for plant somatic embryogenesis (SE), and the accurate quantification of gene expression levels relies on choosing appropriate reference genes. To select the most suitable reference genes for SE studies, 10 commonly used reference genes were examined in synchronized somatic embryogenic and subsequent germinative cultures of Liriodendron hybrids by using quantitative real-time reverse transcription PCR. Four popular normalization algorithms: geNorm, NormFinder, Bestkeeper and Delta-Ct were used to select and validate the suitable reference genes. The results showed that elongation factor 1-gamma, histone H1 linker protein, glyceraldehyde-3-phosphate dehydrogenase and α-tubulin were suitable for SE tissues, while elongation factor 1-gamma and actin were best for the germinative organ tissues. Our work will benefit future studies of gene expression and functional analyses of SE in Liriodendron hybrids. It is also serves as a guide of reference gene selection in early embryonic gene expression analyses for other woody plant species.

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.

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.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

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.

Genome-wide identification of miRNAs and their targets during early somatic embryogenesis in Dimocarpus longan Lour

miRNAs are endogenous regulatory factors that play pivotal roles in post-transcriptional regulation. However, their specific roles in early somatic embryogenesis (SE) remain unclear. Study of the SE system is fundamental for clarifying the molecular mechanisms in Dimocarpus longan. We identified 289 known miRNAs from 106 different miRNA families and 1087 novel miRNAs during early longan SE, including embryogenic callus (EC), incomplete pro-embryogenic culture (ICpEC), globular embryo (GE), and non-embryogenic callus (NEC). The abundances of known miRNAs were concentrated in GE. The differentially expression (DE) miRNAs showed five expression patterns during early SE. Largely miRNAs were expressed highly and specially in EC, ICpEC, and GE, respectively. Some miRNAs and putative target genes were enriched in lignin metabolism. Most potential targets were related to the pathways of plant hormone signal transduction, alternative splicing, tyrosine metabolism and sulfur metabolism in early longan SE. The regulatory relationships between dlo-miR166a-3p and DlHD-zip8, dlo-miR397a and DlLAC7, dlo-miR408-3p and DlLAC12 were confirmed by RNA ligase-mediated rapid amplification of cDNA ends. The expression patterns of eight DE miRNAs detected by qRT-PCR were consistent with RNA-seq. Finally, the miRNA regulatory network in early SE was constructed, which provided new insight into molecular mechanism of early SE in longan.

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.

Full-Length Transcriptome Analysis of the ABCB, PIN/PIN-LIKES, and AUX/LAX Families Involved in Somatic Embryogenesis of Lilium pumilum DC. Fisch

Plant cell totipotency is one of the 25 major topics in current scientific research, and somatic embryos are good experimental material for studying cell totipotency. Polar auxin transport plays an important regulatory role in somatic embryogenesis (SE). However, little is known about the auxin transport genes and their regulatory mechanisms in Lilium SE. In this study, we applied single-molecule real-time (SMRT) sequencing to Lilium pumilum DC. Fisch. for the first time and obtained a total of 119,649 transcripts, of which 14 encoded auxin transport genes. Correlation analyses between somatic embryo induction and gene expression under different treatments revealed that auxin transport genes, especially ATP-binding cassette (ABC) transporter B family member 21 (ABCB21) and PIN-FORMED (PIN) LIKES 7 (PILS7), may be key players in SE, and the necessary duration of picloram (PIC) treatment to induce SE is as short as 3 days. Our research provides valuable genetic information on Lilium pumilum, elucidating the candidate auxin transport genes involved in SE and their influencing factors. This study lays a foundation for elucidating the regulatory mechanism of auxin transport in SE.

Global scale transcriptome analysis reveals differentially expressed genes involve in early somatic embryogenesis in Dimocarpus longan Lour

BACKGROUND

Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Longan SE has been established and wildly used as model system for studying embryogenesis in woody plants, SE-related genes had been characterized. In spite of that, a comprehensive overview of SE at a molecular level is still absent. To understand the molecular mechanisms during longan SE, we examined the transcriptome changes by using Illumina HiSeq from the four distinct developmental stages, including non-embryogenic callus (NEC), embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), globular embryos (GE).

RESULTS

RNA-seq of the four samples generated a total of 243.78 million high quality reads, approximately 81.5% of the data were mapped to longan genome. The cDNA libraries of NEC, EC, ICpEC and GE, generated 22,743, 19,745, 21,144, 21,102 expressed transcripts, 1935, 1710, 1816, 1732 novel transcripts, 2645, 366, 505, 588 unique genes, respectively. Comparative transcriptome analysis showed that a total of 10,642, 4180, 5846 and 1785 genes were differentially expressed in the pairwise comparisons of NEC_vs_EC, EC_vs_ICpEC, EC_vs_GE, ICpEC_vs_GE, respectively. Among them, plant hormones signalling related genes were significantly enriched, especially the auxin and cytokinin signalling components. The transcripts of flavonoid biosynthesis related genes were mainly expressed in NEC, while fatty acid biosynthesis related genes mainly accumulated in early SE. In addition, the extracelluar protein encoding genes LTP, CHI, GLP, AGP, EP1 were related to longan SE. Combined with the FPKM value of longan nine tissues transcription, 27 SE specific or preferential genes (LEC1, LEC1-like, PDF1.3, GH3.6, AGL80, PIN1, BBM, WOX9, WOX2, ABI3, et al.) and 28 NEC preferential genes (LEA5, CNOT3, DC2.15, PR1-1, NsLTP2, DIR1, PIP1, PIP2.1, TIP2-1, POD-P7 and POD5 et al.) were characterized as molecular markers for longan early SE. qRT-PCR validation of SE-related genes showed a high correlation between RNA-seq and qRT-PCR data.

CONCLUSION

This study provides new insights into the role of the transcriptome during early SE in longan. Differentially expressed genes reveal that plant hormones signalling, flavonoid and fatty acid biosynthesis, and extracelluar protein related genes were involved in longan early SE. It could serve as a valuable platform resource for further functional studies addressing embryogenesis in woody plants.

Transcriptome-wide miRNA identification of Bacopa monnieri: a cross-kingdom approach

Bacopa monnieri known as 'Brahmi' is a well-known medicinal plant belonging to Scrophulariaceae family for its nootropic properties. To the best of our knowledge, no characterization data is available on the potential role of micro RNAs (miRNAs) from this plant till date. We present here the first report of computational characterizations of miRNAs from B. monnieri. Owing to the high conservation of miRNAs in nature, new and potential miRNAs can be identified in plants using in silico techniques. Using the plant miRNA sequences present in the miRBase repository, a total of 12 miRNAs were identified from B. monnieri which pertained to 11 miRNA families from the shoot and root transcriptome data. Furthermore, gene ontology analysis of the identified 68 human target genes exhibited significance in various biological processes. These human target genes were associated with signaling pathways like NF-kB and MAPK with TRAF2, CBX1, IL1B, ITGA4 and ITGB1BP1 as the top five hub nodes. This cross-kingdom study provides initial insights about the potential of miRNA-mediated cross-kingdom regulation and unravels the essential target genes of human with implications in numerous human diseases including cancer.

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 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.

Unravelling miRNA regulation in yield of rice (Oryza sativa) based on differential network model

Rice (Oryza sativa L.) is one of the essential staple food crops and tillering, panicle branching and grain filling are three important traits determining the grain yield. Although miRNAs have been reported being regulating yield, no study has systematically investigated how miRNAs differentially function in high and low yield rice, in particular at a network level. This abundance of data from high-throughput sequencing provides an effective solution for systematic identification of regulatory miRNAs using developed algorithms in plants. We here present a novel algorithm, Gene Co-expression Network differential edge-like transformation (GRN-DET), which can identify key regulatory miRNAs in plant development. Based on the small RNA and RNA-seq data, miRNA-gene-TF co-regulation networks were constructed for yield of rice. Using GRN-DET, the key regulatory miRNAs for rice yield were characterized by the differential expression variances of miRNAs and co-variances of miRNA-mRNA, including osa-miR171 and osa-miR1432. Phytohormone cross-talks (auxin and brassinosteroid) were also revealed by these co-expression networks for the yield of rice.

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.