NATpipe: an integrative pipeline for systematical discovery of natural antisense transcripts (NATs) and phase-distributed nat-siRNAs from de novo assembled transcriptomes

Small Interfering RNAs as Critical Regulators of Plant Life Process: New Perspectives on Regulating the Transcriptomic Machinery

Small interfering RNAs (siRNAs) are a distinct class of regulatory RNAs in plants and animals. Gene silencing by small interfering RNAs is one of the fundamental mechanisms for regulating gene expression. siRNAs are critical regulators during developmental processes. siRNAs have similar structures and functions to small RNAs but are derived from double-stranded RNA and may be involved in directing DNA methylation of target sequences. siRNAs are a less well-studied class than the miRNA group, and researchers continue to identify new classes of siRNAs that appear at specific developmental stages and in particular tissues, revealing a more complex mode of siRNA action than previously thought. This review characterizes the siRNA classes and their biogenesis process and focuses on presenting their known functions in the regulation of plant development and responses to biotic and abiotic stresses. The review also highlights the exciting potential for future research in this field, proposing methods for detecting plant siRNAs and a bioinformatic pathway for identifying siRNAs and their functions.

Genome-Wide Identification and Characterization of miRNAs and Natural Antisense Transcripts Show the Complexity of Gene Regulatory Networks for Secondary Metabolism in Aristolochia contorta

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene-enzyme genes, NAT-STs, and NAT-miRNA-target gene-enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

CCIVR2 facilitates comprehensive identification of both overlapping and non-overlapping antisense transcripts within specified regions

Pairs of sense and antisense transcriptions that are adjacent at their 5' and 3' regions are called divergent and convergent transcription, respectively. However, the structural properties of divergent/convergent transcription in different species or RNA biotypes are poorly characterized. Here, we developed CCIVR2, a program that facilitates identification of both overlapping and non-overlapping antisense transcripts produced from divergent/convergent transcription whose transcription start sites (TSS) or transcript end sites (TES) are located within a specified region. We used CCIVR2 to analyze antisense transcripts starting around the sense TSS (from divergent transcription) or ending around the sense TES (from convergent transcription) in 11 different species and found species- and RNA biotype-specific features of divergent/convergent transcription. Furthermore, we confirmed that CCIVR2 enables the identification of multiple sense/antisense transcript pairs from divergent transcription, including those with known functions in processes such as embryonic stem cell differentiation and TGFβ stimulation. CCIVR2 is therefore a valuable bioinformatics tool that facilitates the characterization of divergent/convergent transcription in different species and aids the identification of functional sense/antisense transcript pairs from divergent transcription in specified biological processes.

Biogenesis, Mode of Action and the Interactions of Plant Non-Coding RNAs

The central dogma of genetics, which outlines the flow of genetic information from DNA to RNA to protein, has long been the guiding principle in molecular biology. In fact, more than three-quarters of the RNAs produced by transcription of the plant genome are not translated into proteins, and these RNAs directly serve as non-coding RNAs in the regulation of plant life activities at the molecular level. The breakthroughs in high-throughput transcriptome sequencing technology and the establishment and improvement of non-coding RNA experiments have now led to the discovery and confirmation of the biogenesis, mechanisms, and synergistic effects of non-coding RNAs. These non-coding RNAs are now predicted to play important roles in the regulation of gene expression and responses to stress and evolution. In this review, we focus on the synthesis, and mechanisms of non-coding RNAs, and we discuss their impact on gene regulation in plants.

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

CCIVR facilitates comprehensive identification of cis-natural antisense transcripts with their structural characteristics and expression profiles

Cis-natural antisense transcripts (cis-NATs) are transcribed from the same genomic locus as their partner gene but from the opposite DNA strand and overlap with the partner gene transcript. Here, we developed a simple and convenient program termed CCIVR (comprehensive cis-NATs identifier via RNA-seq data) that comprehensively identifies all kinds of cis-NATs based on genome annotation with expression data obtained from RNA-seq. Using CCIVR with genome databases, we demonstrated total cis-NAT pairs from 11 model organisms. CCIVR analysis with RNA-seq data from parthenogenetic and androgenetic embryonic stem cells identified well-known imprinted cis-NAT pair, KCNQ1/KCNQ1OT1, ensuring the availability of CCIVR. Finally, CCIVR identified cis-NAT pairs that demonstrate inversely correlated expression upon TGFβ stimulation including cis-NATs that functionally repress their partner genes by introducing epigenetic alteration in the promoters of partner genes. Thus, CCIVR facilitates the investigation of structural characteristics and functions of cis-NATs in numerous processes in various species.

Genome-wide identification of small interfering RNAs from sRNA libraries constructed from soybean cyst nematode resistant and susceptible cultivars

Plant small-RNAs regulate various biological processes by manipulating the expression of target genes at the transcriptional and post-transcriptional levels. However, little is known about the response and the functional roles of sRNAs, particularly small-interfering RNAs (siRNAs), in the soybean-soybean cyst nematode interaction. In this study, siRNA data from 24 sRNA libraries constructed from SCN-infected and non-SCN-infected resistant and susceptible soybean roots were analysed in silico. A total of 26 novel siRNAs including 17 phasiRNAs and 9 nat-siRNAs, as well as two phasiRNAs that were differentially expressed (DE) in three comparisons, were identified. Then, using qRT-PCR, the expression of majority of siRNAs was found to be downregulated after SCN infection, and the expression patterns of DE siRNAs were confirmed. Further functional annotation analyses revealed that the target genes of these siRNA were highly related to disease resistance, which included the genes coding for the NB-ARC domain, leucine-rich repeats, and Hs1pro-1 homologous proteins. Overall, the present research identified novel siRNAs and annotated their target genes, thereby laying the foundation for deciphering the roles of siRNAs in the soybean-SCN interaction.

Biogenesis, Functions, Interactions, and Resources of Non-Coding RNAs in Plants

Plant transcriptomes encompass a large number of functional non-coding RNAs (ncRNAs), only some of which have protein-coding capacity. Since their initial discovery, ncRNAs have been classified into two broad categories based on their biogenesis and mechanisms of action, housekeeping ncRNAs and regulatory ncRNAs. With advances in RNA sequencing technology and computational methods, bioinformatics resources continue to emerge and update rapidly, including workflow for in silico ncRNA analysis, up-to-date platforms, databases, and tools dedicated to ncRNA identification and functional annotation. In this review, we aim to describe the biogenesis, biological functions, and interactions with DNA, RNA, protein, and microorganism of five major regulatory ncRNAs (miRNA, siRNA, tsRNA, circRNA, lncRNA) in plants. Then, we systematically summarize tools for analysis and prediction of plant ncRNAs, as well as databases. Furthermore, we discuss the silico analysis process of these ncRNAs and present a protocol for step-by-step computational analysis of ncRNAs. In general, this review will help researchers better understand the world of ncRNAs at multiple levels.

Biogenesis, Trafficking, and Function of Small RNAs in Plants

Small RNAs (sRNAs) encoded by plant genomes have received widespread attention because they can affect multiple biological processes. Different sRNAs that are synthesized in plant cells can move throughout the plants, transport to plant pathogens via extracellular vesicles (EVs), and transfer to mammals via food. Small RNAs function at the target sites through DNA methylation, RNA interference, and translational repression. In this article, we reviewed the systematic processes of sRNA biogenesis, trafficking, and the underlying mechanisms of its functions.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional medicines. Natural antisense transcripts (NATs) are a class of long noncoding RNAs that can regulate gene expression. Here, we identified 812 NATs, including 168 cis-NATs and 644 trans-NATs from twelve root, flower, and leaf samples of S. miltiorrhiza using RNA-seq. The expression profiles for 41 of 50 NATs and their sense transcripts (STs) obtained from RNA-Seq were validated using qRT-PCR. The expression profiles of 17 NATs positively correlated with their STs. GO and KEGG pathway analyses mapped the STs for cis-NATs to pathways for biosynthesis of secondary metabolites. We characterized four NATs in detail, including NAT0001, NAT0002, NAT0004, and NAT00023. Their STs are kaurene synthase-like 1 and the homologs of UDP-glucose flavonoid 3-O-glucosyltransferase 6, UDP-glycosyltransferase 90A1, and beta-glucosidase 40, respectively. The first gene is involved in the biosynthesis of bioactive tanshinones, the next two are involved in anthocyanin biosynthesis, whereas the last is involved in phenylpropanoid biosynthesis. Besides, we found seven STs that are potential targets of miRNAs. And we found two miRNAs including miR156a and miR7208, might originate from NATs, NAT0112 and NAT0086. The results suggest that S. miltiorrhiza NATs might interact with STs, produce miRNAs, and be regulated by miRNAs. They potentially play significant regulatory roles in the biosynthesis of bioactive compounds.

NATpare: a pipeline for high-throughput prediction and functional analysis of nat-siRNAs

Natural antisense transcript-derived small interfering RNAs (nat-siRNAs) are a class of functional small RNA (sRNA) that have been found in both plant and animals kingdoms. In plants, these sRNAs have been shown to suppress the translation of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex (RISC) to their sequence-specific mRNA target(s). Current computational tools for classification of nat-siRNAs are limited in number and can be computationally infeasible to use. In addition, current methods do not provide any indication of the function of the predicted nat-siRNAs. Here, we present a new software pipeline, called NATpare, for prediction and functional analysis of nat-siRNAs using sRNA and degradome sequencing data. Based on our benchmarking in multiple plant species, NATpare substantially reduces the time required to perform prediction with minimal resource requirements allowing for comprehensive analysis of nat-siRNAs in larger and more complex organisms for the first time. We then exemplify the use of NATpare by identifying tissue and stress specific nat-siRNAs in multiple Arabidopsis thaliana datasets.

The RNA degradome: a precious resource for deciphering RNA processing and regulation codes in plants

The plant RNA degradome was defined as an aggregate of the RNA fragments degraded from various biochemical pathways, such as RNA turnover, maturation and quality surveillance. In recent years, the degradome sequencing (degradome-seq) libraries became a rich storehouse for researchers to study on RNA processing and regulation. Here, we provided a brief overview of the uses of degradome-seq data in plant RNA biology, especially on non-coding RNA processing and small RNA-guided target cleavages. Some novel applications in RNA research area, such as in vivo mapping of the endoribonucleolytic cleavage sites, identification of conserved motifs at the 5' ends of the uncapped RNA fragments, and searching for the protein-binding regions on the transcripts, were also mentioned. More importantly, we proposed a model for the biologists to deduce the contributions of transcriptional and/or post-transcriptional regulation to gene differential expression based on degradome-seq data. Finally, we hope that the degradome-based analytical methods could be widely applied for the studies on RNA biology in eukaryotes.

Computational tools for plant small RNA detection and categorization

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

Bioinformatics-assisted, integrated omics studies on medicinal plants

The immense therapeutic and economic values of medicinal plants have attracted increasing attention from the worldwide researchers. It has been recognized that production of the authentic and high-quality herbal drugs became the prerequisite for maintaining the healthy development of the traditional medicine industry. To this end, intensive research efforts have been devoted to the basic studies, in order to pave a way for standardized authentication of the plant materials, and bioengineering of the metabolic pathways in the medicinal plants. In this paper, the recent advances of omics studies on the medicinal plants were summarized from several aspects, including phenomics and taxonomics, genomics, transcriptomics, proteomics and metabolomics. We proposed a multi-omics data-based workflow for medicinal plant research. It was emphasized that integration of the omics data was important for plant authentication and mechanistic studies on plant metabolism. Additionally, the computational tools for proper storage, efficient processing and high-throughput analyses of the omics data have been introduced into the workflow. According to the workflow, authentication of the medicinal plant materials should not only be performed at the phenomics level but also be implemented by genomic and metabolomic marker-based examination. On the other hand, functional genomics studies, transcriptional regulatory networks and protein-protein interactions will contribute greatly for deciphering the secondary metabolic pathways. Finally, we hope that our work could inspire further efforts on the bioinformatics-assisted, integrated omics studies on the medicinal plants.

OverGeneDB: a database of 5' end protein coding overlapping genes in human and mouse genomes

Gene overlap plays various regulatory functions on transcriptional and post-transcriptional levels. Most current studies focus on protein-coding genes overlapping with non-protein-coding counterparts, the so called natural antisense transcripts. Considerably less is known about the role of gene overlap in the case of two protein-coding genes. Here, we provide OverGeneDB, a database of human and mouse 5′ end protein-coding overlapping genes. The database contains 582 human and 113 mouse gene pairs that are transcribed using overlapping promoters in at least one analyzed library. Gene pairs were identified based on the analysis of the transcription start site (TSS) coordinates in 73 human and 10 mouse organs, tissues and cell lines. Beside TSS data, resources for 26 human lung adenocarcinoma cell lines also contain RNA-Seq and ChIP-Seq data for seven histone modifications and RNA Polymerase II activity. The collected data revealed that the overlap region is rarely conserved between the studied species and tissues. In ∼50% of the overlapping genes, transcription started explicitly in the overlap regions. In the remaining half of overlapping genes, transcription was initiated both from overlapping and non-overlapping TSSs. OverGeneDB is accessible at http://overgenedb.amu.edu.pl.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Genome-wide researches and applications on Dendrobium

This review summarizes current knowledge of chromosome characterization, genetic mapping, genomic sequencing, quality formation, floral transition, propagation, and identification in Dendrobium. The widely distributed Dendrobium has been studied for a long history, due to its important economic values in both medicine and ornamental. In recent years, some species of Dendrobium and other orchids had been reported on genomic sequences, using the next-generation sequencing technology. And the chloroplast genomes of many Dendrobium species were also revealed. The chromosomes of most Dendrobium species belong to mini-chromosomes, and showed 2n = 38. Only a few of genetic studies were reported in Dendrobium. After revealing of genomic sequences, the techniques of transcriptomics, proteomics and metabolomics could be employed on Dendrobium easily. Some other molecular biological techniques, such as gene cloning, gene editing, genetic transformation and molecular marker developing, had also been applied on the basic research of Dendrobium, successively. As medicinal plants, insights into the biosynthesis of some medicinal components were the most important. As ornamental plants, regulation of flower related characteristics was the most important. More, knowledge of growth and development, environmental interaction, evolutionary analysis, breeding of new cultivars, propagation, and identification of species and herbs were also required for commercial usage. All of these studies were improved using genomic sequences and related technologies. To answer some key scientific issues in Dendrobium, quality formation, flowering, self-incompatibility and seed germination would be the focus of future research. And genome related technologies and studies would be helpful.

Natural Antisense Transcripts: Molecular Mechanisms and Implications in Breast Cancers

Natural antisense transcripts are RNA sequences that can be transcribed from both DNA strands at the same locus but in the opposite direction from the gene transcript. Because strand-specific high-throughput sequencing of the antisense transcriptome has only been available for less than a decade, many natural antisense transcripts were first described as long non-coding RNAs. Although the precise biological roles of natural antisense transcripts are not known yet, an increasing number of studies report their implication in gene expression regulation. Their expression levels are altered in many physiological and pathological conditions, including breast cancers. Among the potential clinical utilities of the natural antisense transcripts, the non-coding|coding transcript pairs are of high interest for treatment. Indeed, these pairs can be targeted by antisense oligonucleotides to specifically tune the expression of the coding-gene. Here, we describe the current knowledge about natural antisense transcripts, their varying molecular mechanisms as gene expression regulators, and their potential as prognostic or predictive biomarkers in breast cancers.

Bioinformatics resources for deciphering the biogenesis and action pathways of plant small RNAs

The next-generation sequencing (NGS) technology has revolutionized our previous understanding of the plant genomes, relying on its innate advantages, such as high throughput and deep sequencing depth. In addition to the protein-coding gene loci, massive transcription signals have been detected within intergenic or intragenic regions. Most of these signals belong to non-coding ones, considering their weak protein-coding potential. Generally, these transcripts could be divided into long non-coding RNAs and small non-coding RNAs (sRNAs) based on their sequence length. In addition to the well-known microRNAs (miRNAs), many plant endogenous sRNAs were collectively referred to as small interfering RNAs. However, an increasing number of unclassified sRNA species are being discovered by NGS. The high heterogeneity of these novel sRNAs greatly hampered the mechanistic studies, especially on the clear description of their biogenesis and action pathways. Fortunately, public databases, bioinformatics softwares and NGS datasets are increasingly available for plant sRNA research. Here, by summarizing these valuable resources, we proposed a general workflow to decipher the RDR (RNA-dependent RNA polymerase)- and DCL (Dicer-like)-dependent biogenesis pathways, and the Argonaute-mediated action modes (such as target cleavages and chromatin modifications) for specific sRNA species in plants. Taken together, we hope that by summarizing a list of the public resources, this work will facilitate the plant biologists to perform classification and functional characterization of the interesting sRNA species.