Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues

The Emerging Role of Non-Coding RNAs (ncRNAs) in Plant Growth, Development, and Stress Response Signaling

Plant species utilize a variety of regulatory mechanisms to ensure sustainable productivity. Within this intricate framework, numerous non-coding RNAs (ncRNAs) play a crucial regulatory role in plant biology, surpassing the essential functions of RNA molecules as messengers, ribosomal, and transfer RNAs. ncRNAs represent an emerging class of regulators, operating directly in the form of small interfering RNAs (siRNAs), microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). These ncRNAs exert control at various levels, including transcription, post-transcription, translation, and epigenetic. Furthermore, they interact with each other, contributing to a variety of biological processes and mechanisms associated with stress resilience. This review primarily concentrates on the recent advancements in plant ncRNAs, delineating their functions in growth and development across various organs such as root, leaf, seed/endosperm, and seed nutrient development. Additionally, this review broadens its scope by examining the role of ncRNAs in response to environmental stresses such as drought, salt, flood, heat, and cold in plants. This compilation offers updated information and insights to guide the characterization of the potential functions of ncRNAs in plant growth, development, and stress resilience in future research.

Identification of antisense and sense RNAs of intracrine fibroblast growth factor components as novel biomarkers in colorectal cancer and in silico studies for drug and nanodrug repurposing

INTRODUCTION

Colorectal cancer (CRC) is one of the most malignant tumors and in which various efforts for screening is inconclusive.The intracrine FGF panel, the non-tyrosine kinase receptors (NTKR) FGFs and affiliated antisenses play a pivotal role in FGF signaling.The expression levels of coding and non-coding intracrine FGFs were assessed in CRC donors.Also, substantial costs and slow pace of drug discovery give high attraction to repurpose of previously discovered drugs to new opportunities.

OBJECTIVES

The aim of present study was to evaluate the potential role of the coding and non-coding intracrine FGFs as a new biomarkers for CRC cases and defining drug repurposing to alleviate FGF down regulation.

METHODS

RNA-seq data of colon adenocarcinomas (COAD) was downloaded using TCGA biolinks package in R.The DrugBank database (https://go.drugbank.com/) was used to extract interactions between drugs and candidate genes. A total of 200 CRC patients with detailed criteria were enrolled.RNAs were extracted with TRIzol-based protocol and amplified via LightCycler® instrument.FGF11 and FGF13 proteins validation was performed by used of immunohistochemistry technique in tumor and non-tumoral samples.Pearson's correlation analysis and ROC curve plotted by Prism 8.0 software.

RESULTS

RNA-seq data from TCGA was analyzed by normalizing with edgeR.Differentially expressed gene (DEG) analysis was generated. WCC algorithm extracted the most significant genes with a total of 47 genes. Expression elevation of iFGF antisenses (12AS,13As,14AS) compared with the normal colon tissue were observed (P = 0.0003,P = 0.042,P = 0.026, respectively). Moreover,a significant decrease in expression of the corresponding sense iFGF genes was detected (P < 0.0001).Plotted receiver operating characteristic (ROC) curves for iFGF components' expression showed an area of over 0.70 (FGF11-13: 0.71% and FGF12-14: 0.78%, P < 0.001) for sense mRNA expression, with the highest sensitivity for FGF12 (92.8%) and lowest for FGF11 (61.41%).The artificial intelligence (AI) revealed the valproic acid as a repurposing drug to relief the down regulation of FGF12 and 13 in CRC patients.

CONCLUSION

Intracrine FGFs panel was down regulated versus up regulation of dependent antisenses. Thus, developing novel biomarkers based on iFGF can be considered as a promising strategy for CRC screening.In advanced, valporic acid detected by AI as a repurposing drug which may be applied in clinical trials for CRC treatment.

Altered expression levels of long non-coding natural antisense transcripts overlapping the UGT73C6 gene affect rosette size in Arabidopsis thaliana

Natural antisense long non-coding RNAs (lncNATs) are involved in the regulation of gene expression in plants, modulating different relevant developmental processes and responses to various stimuli. We have identified and characterized two lncNATs (NAT1_UGT73C6 and NAT2_UGT73C6 , collectively NATs_UGT73C6 ) from Arabidopsis thaliana that are transcribed from a gene fully overlapping UGT73C6, a member of the UGT73C subfamily of genes encoding UDP-glycosyltransferases (UGTs). Expression of both NATs_UGT73C6 is developmentally controlled and occurs independently of the transcription of UGT73C6 in cis. Downregulation of NATs_UGT73C6 levels through artificial microRNAs results in a reduction of the rosette area, while constitutive overexpression of NAT1_UGT73C6 or NAT2_UGT73C6 leads to the opposite phenotype, an increase in rosette size. This activity of NATs_UGT73C6 relies on its RNA sequence and, although modulation of UGT73C6 in cis cannot be excluded, the observed phenotypes are not a consequence of the regulation of UGT73C6 in trans. The NATs_UGT73C6 levels were shown to affect cell proliferation and thus individual leaf size. Consistent with this concept, our data suggest that the NATs_UGT73C6 influence the expression levels of key transcription factors involved in regulating leaf growth by modulating cell proliferation. These findings thus reveal an additional regulatory layer on the process of leaf growth. In this work, we characterized at the molecular level two long non-coding RNAs (NATs_UGT73C6 ) that are transcribed in the opposite direction to UGT73C6, a gene encoding a glucosyltransferase involved in brassinosteroid homeostasis in A. thaliana. Our results indicate that NATs_UGT73C6 expression influences leaf growth by acting in trans and by modulating the levels of transcription factors that are involved in the regulation of cell proliferation.

The non-coding RNA SVALKA locus produces a cis-natural antisense transcript that negatively regulates the expression of CBF1 and biomass production at normal temperatures

Non-coding transcription is present in all eukaryotic genomes, but we lack fundamental knowledge about its importance for an organism's ability to develop properly. In plants, emerging evidence highlights the essential biological role of non-coding transcription in the regulation of coding transcription. However, we have few molecular insights into this regulation. Here, we show that a long isoform of the long non-coding RNA SVALKA-L (SVK-L) forms a natural antisense transcript to the host gene CBF1 and negatively regulates CBF1 mRNA levels at normal temperatures in the model plant Arabidopsis thaliana. Furthermore, we show detailed evidence for the specific mode of action of SVK-L. This pathway includes the formation of double-stranded RNA that is recognized by the DICER proteins and subsequent downregulation of CBF1 mRNA levels. Thus, the CBF1-SVK regulatory circuit is not only important for its previously known role in cold temperature acclimation but also for biomass production at normal temperatures. Our study characterizes the developmental role of SVK-L and offers mechanistic insight into how biologically important overlapping natural antisense transcripts can act on and fine-tune the steady-state levels of their host gene's mRNA.

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Non-coding RNAs fine-tune the balance between plant growth and abiotic stress tolerance

To survive in adverse environmental conditions, plants have evolved sophisticated genetic and epigenetic regulatory mechanisms to balance their growth and abiotic stress tolerance. An increasing number of non-coding RNAs (ncRNAs), including small RNAs (sRNAs) and long non-coding RNAs (lncRNAs) have been identified as essential regulators which enable plants to coordinate multiple aspects of growth and responses to environmental stresses through modulating the expression of target genes at both the transcriptional and posttranscriptional levels. In this review, we summarize recent advances in understanding ncRNAs-mediated prioritization towards plant growth or tolerance to abiotic stresses, especially to cold, heat, drought and salt stresses. We highlight the diverse roles of evolutionally conserved microRNAs (miRNAs) and small interfering RNAs (siRNAs), and the underlying phytohormone-based signaling crosstalk in regulating the balance between plant growth and abiotic stress tolerance. We also review current discoveries regarding the potential roles of ncRNAs in stress memory in plants, which offer their descendants the potential for better fitness. Future ncRNAs-based breeding strategies are proposed to optimize the balance between growth and stress tolerance to maximize crop yield under the changing climate.

Tissue-specific transcriptomic analysis uncovers potential roles of natural antisense transcripts in Arabidopsis heat stress response

Natural antisense transcripts (NATs) are an important class of non-coding ribonucleic acids (RNAs) that have been shown to regulate gene expression. Using strand-specific RNA sequencing, 36,317 NAT pairs were identified, and 5,536 were specifically expressed under heat stress. We found distinct expression patterns between vegetative and reproductive tissues for both coding genes and genes encoding NATs. Genes for heat-responsive NATs are associated with relatively high levels of H3K4me3 and low levels of H3K27me2/3. On the other hand, small RNAs are significantly enriched in sequence overlapping regions of NAT pairs, and a large number of heat-responsive NATs pairs serve as potential precursors of nat-siRNAs. Collectively, our results suggest epigenetic modifications and small RNAs play important roles in the regulation of NAT expression, and highlight the potential significance of heat-inducible NATs.

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.

Regulation of small RNA-mediated high temperature stress responses in crop plants

Small RNAs have emerged as key players of gene expression regulation. Several lines of evidences highlight their role in modulating high temperature stress responsiveness in plants. Throughout their life cycle, plants have to regulate their gene expression at various developmental phases, physiological changes, and in response to biotic or environmental stress. High temperature is one the most common abiotic stress for crop plants, that results in impaired morphology, physiology, and yield. However, plants have certain mechanisms that enable them to withstand such conditions by modulating the expression of stress-related genes. Small RNA (sRNA)-regulated gene expression is one such mechanism which is ubiquitous in all eukaryotes. The sRNAs mainly include micro RNAs (miRNAs) and small interfering RNAs (siRNAs). They are primarily associated with the gene silencing either through translation inhibition, mRNA degradation, or DNA methylation. During high temperature stress the increased or decreased level of miRNAs altered the protein accumulation of target transcripts and, therefore, regulate stress responses. Several reports are available in plants which are genetically engineered through expressing artificial miRNAs resulted in thermotolerance. sRNAs have also been reported to bring the epigenetic changes on chromatin region through RNA-dependent DNA methylation (RdDM). The present article draws a brief illustration of sRNA origin, their functional mechanisms, role in high temperature stress, and possible application for developing stress tolerant crop plants.

Small RNAs Participate in Plant-Virus Interaction and Their Application in Plant Viral Defense

Small RNAs are significant regulators of gene expression, which play multiple roles in plant development, growth, reproductive and stress response. It is generally believed that the regulation of plants’ endogenous genes by small RNAs has evolved from a cellular defense mechanism for RNA viruses and transposons. Most small RNAs have well-established roles in the defense response, such as viral response. During viral infection, plant endogenous small RNAs can direct virus resistance by regulating the gene expression in the host defense pathway, while the small RNAs derived from viruses are the core of the conserved and effective RNAi resistance mechanism. As a counter strategy, viruses evolve suppressors of the RNAi pathway to disrupt host plant silencing against viruses. Currently, several studies have been published elucidating the mechanisms by which small RNAs regulate viral defense in different crops. This paper reviews the distinct pathways of small RNAs biogenesis and the molecular mechanisms of small RNAs mediating antiviral immunity in plants, as well as summarizes the coping strategies used by viruses to override this immune response. Finally, we discuss the current development state of the new applications in virus defense based on small RNA silencing.

A natural antisense transcript acts as a negative regulator for the maize drought stress response gene ZmNAC48

Although plant-specific NAC transcription factors play crucial roles in response to abiotic stress, few reports describe the regulation of NAC genes in maize (Zea mays) by the cis-natural antisense transcripts (cis-NATs). In this study, 521 NAC genes from Gramineae were classified, of which 51 NAC genes contained cis-NATs. ZmNAC48 and cis-NATZmNAC48 co-localized to the same cell nucleus, and both transcripts responded to drought stress. Arabidopsis plants overexpressing ZmNAC48 had improved drought tolerance, lower rate of water loss, enhanced stomatal closure, and higher rates of survival. Transient expression in both maize protoplasts and tobacco leaves indicated that cis-NATZmNAC48 reduced ZmNAC48 expression. Western blotting and ribosome profiling analyses confirmed that cis-NATZmNAC48 lacked protein coding potential. Furthermore, the cis-NAT-derived small-interfering RNAs (nat-siRNAs) generated from the overlapping regions of ZmNAC48 and cis-NATZmNAC48 were detected in maize and transgenic Arabidopsis. Cis-NATZmNAC48 overexpressing maize showed higher water loss rate, increased stomatal opening, and had more dead leaves. Expression of ZmNAC48 and nat-siRNA was decreased in these plants. Taken together, our study indicates that both ZmNAC48 and cis-NATZmNAC48 are involved in plant drought stress responses, and that the double-stranded RNA-dependent mechanism is involved in the interaction between cis-NATZmNAC48 and ZmNAC48. Additionally, cis-NATZmNAC48 may negatively regulate ZmNAC48 to affect stomatal closure of maize.

Secondary siRNAs in Plants: Biosynthesis, Various Functions, and Applications in Virology

The major components of RNA silencing include both transitive and systemic small RNAs, which are technically called secondary sRNAs. Double-stranded RNAs trigger systemic silencing pathways to negatively regulate gene expression. The secondary siRNAs generated as a result of transitive silencing also play a substantial role in gene silencing especially in antiviral defense. In this review, we first describe the discovery and pathways of transitivity with emphasis on RNA-dependent RNA polymerases followed by description on the short range and systemic spread of silencing. We also provide an in-depth view on the various size classes of secondary siRNAs and their different roles in RNA silencing including their categorization based on their biogenesis. The other regulatory roles of secondary siRNAs in transgene silencing, virus-induced gene silencing, transitivity, and trans-species transfer have also been detailed. The possible implications and applications of systemic silencing and the different gene silencing tools developed are also described. The details on mobility and roles of secondary siRNAs derived from viral genome in plant defense against the respective viruses are presented. This entails the description of other compatible plant-virus interactions and the corresponding small RNAs that determine recovery from disease symptoms, exclusion of viruses from shoot meristems, and natural resistance. The last section presents an overview on the usefulness of RNA silencing for management of viral infections in crop plants.

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.

Global identification of a marine diatom long noncoding natural antisense transcripts (NATs) and their response to phosphate fluctuations

Often ignored and regarded as mere transcriptional noise, long noncoding RNAs (lncRNAs) are starting to be considered key regulators of gene expression across the Eukarya domain of life. In the model diatom Phaeodactylum tricornutum, we have previously reported the occurrence of 1,510 intergenic lncRNAs (lincRNAs), many of which displaying specific patterns of expression under phosphate fluctuation (Pi). Using strand-specific RNA-sequencing data we now expand the repertoire of P. tricornutum lncRNAs by identifying 2,628 novel natural antisense transcripts (NATs) that cover 21.5% of the annotated genomic loci. We found that NAT expression is tightly regulated by phosphate depletion and other naturally occurring environmental stresses. Furthermore, we identified 121 phosphate stress responsive NAT-mRNA pairs, the great majority of which showing a positive correlation (concordant pairs) and a small fraction with negative correlation (discordant pairs). Taken together our results show that NATs are highly abundant transcripts in P. tricornutum and that their expression is under tight regulation by nutrient and environmental stresses. Furthermore, our results suggest that in P. tricornutum Pi stress response NAT pairs predominantly regulate positively the expression of their cognate sense genes, the latter being involved in several biological processes underlying the control of cellular homeostasis under stress.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Full-Length Transcript-Based Proteogenomics of Rice Improves Its Genome and Proteome Annotation

Rice (Oryza sativa) molecular breeding has gained considerable attention in recent years, but inaccurate genome annotation hampers its progress and functional studies of the rice genome. In this study, we applied single-molecule long-read RNA sequencing (lrRNA_seq)-based proteogenomics to reveal the complexity of the rice transcriptome and its coding abilities. Surprisingly, approximately 60% of loci identified by lrRNA_seq are associated with natural antisense transcripts (NATs). The high-density genomic arrangement of NAT genes suggests their potential roles in the multifaceted control of gene expression. In addition, a large number of fusion and intergenic transcripts have been observed. Furthermore, 906,456 transcript isoforms were identified, and 72.9% of the genes can generate splicing isoforms. A total of 706,075 posttranscriptional events were subsequently categorized into 10 subtypes, demonstrating the interdependence of posttranscriptional mechanisms that contribute to transcriptome diversity. Parallel short-read RNA sequencing indicated that lrRNA_seq has a superior capacity for the identification of longer transcripts. In addition, over 190,000 unique peptides belonging to 9,706 proteoforms/protein groups were identified, expanding the diversity of the rice proteome. Our findings indicate that the genome organization, transcriptome diversity, and coding potential of the rice transcriptome are far more complex than previously anticipated.

Transcriptome-wide-scale-predicted dsRNAs potentially involved in RNA homoeostasis are remarkably excluded from genes with no/very low expression in all developmental stages

RNA interference (RNAi) refers to a conserved posttranscriptional mechanism for the degradation of RNA by short dsRNAs. A genome-wide analysis of mRNAs that are complementary to RNAs of variable length that are transcribed from the full transcriptome and susceptible to being loaded onto Argonaute type 2 was performed through computational searches in the Drosophila model. We report the segments of RNAs that are complementary to mRNAs originating from introns, the exons of mRNAs and lncRNAs as a potential source of siRNAs. A full catalogue of the mRNAs that fulfill these criteria is presented, along with the quantification of multiple annealing. The catalogue was assessed for biological validation using three published lists: two for Ago2-associated RNAs and one for dsRNAs isolated from a crude extract. A broad spectrum of mRNAs were found to theoretically form intermolecular segmental dsRNAs, which should qualify them as Dicer/Ago2 substrates if they exist in vivo. These results suggest a genome-wide scale of mRNA homoeostasis via RNAi metabolism and could extend the known roles of canonical miRNAs and hairpin RNAs. The distribution of the genes for which transcripts are engaged in intermolecular segmental pairing is largely lacking in the gene collections defined as showing no expression in each individual developmental stage from early embryos to adulthood. This trend was also observed for the genes showing very low expression from the 8-12-hour embryonic to larval stage 2. This situation was also suggested by the 3 lists generated with minimal 20-, 25- and 30-base pairing lengths.

Deep-sequencing of Solanum commersonii small RNA libraries reveals riboregulators involved in cold stress response

Among wild species used in potato breeding, Solanum commersonii displays the highest tolerance to low temperatures under both acclimated (ACC) and non-acclimated (NACC) conditions. It is also the first wild potato relative with a known whole genome sequence. Recent studies have shown that abiotic stresses induce changes in the expression of many small non-coding RNA (sncRNA). We determined the small non-coding RNA (sncRNAome) of two clones of S. commersonii contrasting in their cold response phenotypes via smRNAseq. Differential analysis provided evidence that expression of several miRNAs changed in response to cold stress conditions. Conserved miR408a and miR408b changed their expression under NACC conditions, whereas miR156 and miR169 were differentially expressed only under ACC conditions. We also report changes in tasiRNA and secondary siRNA expression under both stress conditions. Our results reveal possible roles of sncRNA in the regulatory networks associated with tolerance to low temperatures and provide useful information for a more strategic use of genomic resources in potato breeding.

Full-length mRNA sequencing and gene expression profiling reveal broad involvement of natural antisense transcript gene pairs in pepper development and response to stresses

Pepper is an important vegetable with great economic value and unique biological features. In the past few years, significant development has been made toward understanding the huge complex pepper genome; however, pepper functional genomics has not been well studied. To better understand the pepper gene structure and pepper gene regulation, we conducted full-length mRNA sequencing by PacBio sequencing and obtained 57 862 high-quality full-length mRNA sequences derived from 18 362 previously annotated and 5769 newly detected genes. New gene models were built that combined the full-length mRNA sequences and corrected approximately 500 fragmented gene models from previous annotations. Based on the full-length mRNA, we identified 4114 and 5880 pepper genes forming natural antisense transcript (NAT) genes in-cis and in-trans, respectively. Most of these genes accumulate small RNAs in their overlapping regions. By analyzing these NAT gene expression patterns in our transcriptome data, we identified many NAT pairs responsive to a variety of biological processes in pepper. Pepper formate dehydrogenase 1 (FDH1), which is required for R-gene-mediated disease resistance, may be regulated by nat-siRNAs and participate in a positive feedback loop in salicylic acid biosynthesis during resistance responses. Several cis-NAT pairs and subgroups of trans-NAT genes were responsive to pepper pericarp and placenta development, which may play roles in capsanthin and capsaicin biosynthesis. Using a comparative genomics approach, the evolutionary mechanisms of cis-NATs were investigated, and we found that an increase in intergenic sequences accounted for the loss of most cis-NATs, while transposon insertion contributed to the formation of most new cis-NATs. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at http://bigd.big.ac.cn/gsa Accession number, CRA001412.

The Involvement of Long Noncoding RNAs in Response to Plant Stress

Plant growth and productivity are greatly impacted by environmental stresses. Therefore, plants have evolved mechanisms which allow them to adapt to abiotic stresses through alterations in gene expression and metabolism. In recent years, studies have investigated the role of long noncoding RNA (lncRNA) in regulating gene expression in plants and characterized their involvement in various biological functions through their regulation of DNA methylation, DNA structural modifications, histone modifications, and RNA-RNA interactions. Genome-wide transcriptome analyses have identified various types of noncoding RNAs (ncRNAs) that respond to abiotic stress. These ncRNAs are in addition to the well-known housekeeping ncRNAs, such as rRNAs, tRNAs, snoRNAs, and snRNAs. In this review, recent research pertaining to the role of lncRNAs in the response of plants to abiotic stress is summarized and discussed.

Evolution and Diversification of Small RNA Pathways in Flowering Plants

Small regulatory RNAs guide gene silencing at the DNA or RNA level through repression of complementary sequences. The two main forms of small RNAs are microRNA (miRNA) and small interfering RNA (siRNAs), which are generated from the processing of different forms of double-stranded RNA (dsRNA) precursors. These two forms of small regulatory RNAs function in distinct but overlapping gene silencing pathways in plants. Gene silencing pathways in eukaryotes evolved from an ancient prokaryotic mechanism involved in genome defense against invasive genetic elements, but has since diversified to also play a crucial role in regulation of endogenous gene expression. Here, we review the biogenesis of the different forms of small RNAs in plants, including miRNAs, phased, secondary siRNAs (phasiRNAs) and heterochromatic siRNAs (hetsiRNAs), with a focus on their functions in genome defense, transcriptional and post-transcriptional gene silencing, RNA-directed DNA methylation, trans-chromosomal methylation and paramutation. We also discuss the important role that gene duplication has played in the functional diversification of gene silencing pathways in plants, and we highlight recently discovered components of gene silencing pathways in plants.

Increased expression of antisense lncRNA SPINT1-AS1 predicts a poor prognosis in colorectal cancer and is negatively correlated with its sense transcript

Purpose

Colorectal cancer (CRC) is a leading cause of cancer-associated mortality worldwide. Natural antisense transcripts (NATs) are pervasively expressed in human genome and have been confirmed to contribute to cancer progression. In our study, we aimed to investigate the expression and clinical pertinence of serine peptidase inhibitor, Kunitz type 1 antisense RNA1 (SPINT1-AS1) in CRC.

Materials and methods

The expression levels of SPINT1-AS1 and the corresponding sense transcript SPINT1 mRNA were analyzed in 150 pairs of CRC tissues and adjacent normal (AN) tissues, along with 45 pairs of preoperative and postoperative serum exosome samples by the strand-specific real-time quantitative polymerase chain reaction.

Results

Compared with AN tissues, the expression of SPINT1-AS1 was increased (P<0.001, 3.771 vs 0.980) in CRC tissues, while SPINT1 mRNA expression was decreased in CRC (P<0.001, 0.927 vs 1.165), and there was an obviously negative correlation between SPINT1-AS1 expression and its sense transcript (r=-0.701, P<0.001). SPINT1-AS1 yielded an area under the receiver operating characteristic curve value of 0.865 (95% confidence interval, 0.821-0.902) for discriminating CRC tissues from AN tissues. Moreover, high SPINT1-AS1 expression was correlated with regional lymph node metastasis (P<0.001), distant metastasis (P<0.001), and shorter relapse-free survival (RFS) time (P<0.001), and Cox regression analysis indicated that SPINT1-AS1 was an independent prognostic factor for RFS. Meanwhile, significant reduction of SPINT1-AS1 expression level (P=0.001) was observed in CRC serum exosomes after surgical resection.

Conclusion

SPINT1-AS1 is upregulated in CRC tissues and plays an essential role in CRC progression and prognosis. Thereby, SPINT1-AS1 may serve as a candidate prognostic biomarker and molecular therapy target for CRC.

DAWDLE Interacts with DICER-LIKE Proteins to Mediate Small RNA Biogenesis

DAWDLE (DDL) is a conserved forkhead-associated (FHA) domain-containing protein with essential roles in plant development and immunity. It acts in the biogenesis of microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs), which regulate gene expression at the transcriptional and/or posttranscriptional levels. However, the functional mechanism of DDL and its impact on exogenous siRNAs remain elusive. Here, we report that DDL is required for the biogenesis of siRNAs derived from sense transgenes and inverted-repeat transgenes. Furthermore, we show that a mutation in the FHA domain of DDL disrupts the interaction of DDL with DICER-LIKE1 (DCL1), which is the enzyme that catalyzes miRNA maturation from primary miRNA transcripts (pri-miRNAs), resulting in impaired pri-miRNA processing. Moreover, we demonstrate that DDL interacts with DCL3, which is a DCL1 homolog responsible for siRNA production, and this interaction is crucial for optimal DCL3 activity. These results reveal that the interaction of DDL with DCLs is required for the biogenesis of miRNAs and siRNAs in Arabidopsis (Arabidopsis thaliana).

Analysis of RDR1/RDR2/RDR6-independent small RNAs in Arabidopsis thaliana improves MIRNA annotations and reveals unexplained types of short interfering RNA loci

Plant small RNAs (sRNAs) modulate key physiological mechanisms through post-transcriptional and transcriptional silencing of gene expression. Small RNAs fall into two major categories: those are reliant on RNA-dependent RNA polymerases (RDRs) for biogenesis and those that are not. Known RDR1/2/6-dependent sRNAs include phased and repeat-associated short interfering RNAs, while known RDR1/2/6-independent sRNAs are primarily microRNAs (miRNA) and other hairpin-derived sRNAs. In this study we produced and analyzed sRNA-seq libraries from rdr1/rdr2/rdr6 triple mutant plants. We found 58 previously annotated miRNA loci that were reliant on RDR1, -2, or -6 function, casting doubt on their classification. We also found 38 RDR1/2/6-independent sRNA loci that are not MIRNAs or otherwise hairpin-derived, and did not fit into other known paradigms for sRNA biogenesis. These 38 sRNA-producing loci have as-yet-undescribed biogenesis mechanisms, and are frequently located in the vicinity of protein-coding genes. Altogether, our analysis suggests that these 38 loci represent one or more undescribed types of sRNA in Arabidopsis thaliana.

Cell-specific cis-natural antisense transcripts (cis-NATs) in the sperm and the pollen vegetative cells of Arabidopsis thaliana

Background: cis-NATs (cis-natural antisense transcripts ) are transcribed from opposite strands of adjacent genes and have been shown to regulate gene expression by generating small RNAs from the overlapping region. cis-NATs are important for plant development and resistance to pathogens and stress. Several genome-wide investigations identified a number of cis-NAT pairs, but these investigations predicted cis-NATS using expression data from bulk samples that included lots of cell types. Some cis-NAT pairs identified from those investigations might not be functional, because both transcripts of cis-NAT pairs need to be co-expressed in the same cell. Pollen only contains two cell types, two sperm and one vegetative cell, which makes cell-specific investigation of cis-NATs possible.

Methods: We investigated potential protein-coding cis-NATs in pollen and sperm using pollen RNA-seq data and TAIR10 gene models using the Integrated Genome Browser.  We then used sperm microarray data and sRNAs in sperm and pollen to determine possibly functional cis-NATs in the sperm or vegetative cell, respectively.

Results: We identified 1471 potential protein-coding cis-NAT pairs, including 131 novel pairs that were not present in TAIR10 gene models. In pollen, 872 possibly functional pairs were identified. 72 and 56 pairs were potentially functional in sperm and vegetative cells, respectively. sRNAs were detected at 794 genes, belonging to 739 pairs.

Conclusion: These potential candidates in sperm and the vegetative cell are tools for understanding gene expression mechanisms in pollen.

Biogenesis and regulatory hierarchy of phased small interfering RNAs in plants

Several varieties of small RNAs including microRNAs (miRNAs) and small interfering RNAs (siRNAs) are generated in plants to regulate development, genome stability and response to adverse environments. Phased siRNA (phasiRNA) is a type of secondary siRNA that is processed from a miRNA-mediated cleavage of RNA transcripts, increasing silencing efficiency or simultaneously suppressing multiple target genes. Trans-acting siRNAs (ta-siRNAs) are a particular class of phasiRNA produced from noncoding transcripts that silence targets in trans. It was originally thought that 'one-hit' and 'two-hit' models were essential for processing distinct TAS precursors; however, a single hit event was recently shown to be sufficient at triggering all types of ta-siRNAs. This review discusses the findings about biogenesis, targeting modes and regulatory networks of plant ta-siRNAs. We also summarize recent advances in the generation of other phasiRNAs and their possible biological benefits to plants.

The antiphasic regulatory module comprising CDF5 and its antisense RNA FLORE links the circadian clock to photoperiodic flowering

Circadian rhythms of gene expression are generated by the combinatorial action of transcriptional and translational feedback loops as well as chromatin remodelling events. Recently, long noncoding RNAs (lncRNAs) that are natural antisense transcripts (NATs) to transcripts encoding central oscillator components were proposed as modulators of core clock function in mammals (Per) and fungi (frq/qrf). Although oscillating lncRNAs exist in plants, their functional characterization is at an initial stage. By screening an Arabidopsis thaliana lncRNA custom-made array we identified CDF5 LONG NONCODING RNA (FLORE), a circadian-regulated lncRNA that is a NAT of CDF5. Quantitative real-time RT-PCR confirmed the circadian regulation of FLORE, whereas GUS-staining and flowering time evaluation were used to determine its biological function. FLORE and CDF5 antiphasic expression reflects mutual inhibition in a similar way to frq/qrf. Moreover, whereas the CDF5 protein delays flowering by directly repressing FT transcription, FLORE promotes it by repressing several CDFs (CDF1, CDF3, CDF5) and increasing FT transcript levels, indicating both cis and trans function. We propose that the CDF5/FLORE NAT pair constitutes an additional circadian regulatory module with conserved (mutual inhibition) and unique (function in trans) features, able to fine-tune its own circadian oscillation, and consequently, adjust the onset of flowering to favourable environmental conditions.

Non-coding RNAs and Their Roles in Stress Response in Plants

Eukaryotic genomes encode thousands of non-coding RNAs (ncRNAs), which play crucial roles in transcriptional and post-transcriptional regulation of gene expression. Accumulating evidence indicates that ncRNAs, especially microRNAs (miRNAs) and long ncRNAs (lncRNAs), have emerged as key regulatory molecules in plant stress responses. In this review, we have summarized the current progress on the understanding of plant miRNA and lncRNA identification, characteristics, bioinformatics tools, and resources, and provided examples of mechanisms of miRNA- and lncRNA-mediated plant stress tolerance.

Natural antisense transcripts are significantly involved in regulation of drought stress in maize

Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency

Long non-coding RNAs (lncRNAs) have been identified as important regulatory factors of gene expression in eukaryotic species, such as Homo sapiens, Arabidopsis thaliana, and Oryza sativa. However, the systematic identification of potential lncRNAs in trees is comparatively rare. In particular, the characteristics, expression, and regulatory roles of lncRNAs in trees under nutrient stress remain largely unknown. A genome-wide strategy was used in this investigation to identify and characterize novel and low-nitrogen (N)-responsive lncRNAs in Populus tomentosa; 388 unique lncRNA candidates belonging to 380 gene loci were detected and only seven lncRNAs were found to belong to seven conserved non-coding RNA families indicating the majority of P. tomentosa lncRNAs are species-specific. In total, 126 lncRNAs were significantly altered under low-N stress; 8 were repressed, and 118 were induced. Furthermore, 9 and 5 lncRNAs were detected as precursors of 11 known and 14 novel Populus miRNAs, respectively, whereas 4 lncRNAs were targeted by 29 miRNAs belonging to 5 families, including 22 conserved and 7 non-conserved miRNAs. In addition, 15 antisense lncRNAs were identified to be generated from opposite strands of 14 corresponding protein-coding genes. In total, 111 protein-coding genes with regions complementary to 38 lncRNAs were also predicted with some lncRNAs corresponding to multiple genes and vice versa, and their functions were annotated, which further demonstrated the complex regulatory relationship between lncRNAs and protein-coding genes in plants. Moreover, an interaction network among lncRNAs, miRNAs, and mRNAs was investigated. These findings enrich our understanding of lncRNAs in Populus, expand the methods of miRNA identification. Our results present the first global characterization of lncRNAs and their potential target genes in response to nitrogen stress in trees, which provides more information on low-nutrition adaptation mechanisms in woody plants.

Small RNAs: essential regulators of gene expression and defenses against environmental stresses in plants

Eukaryotic genomes produce thousands of diverse small RNAs (smRNAs), which play vital roles in regulating gene expression in all conditions, including in survival of biotic and abiotic environmental stresses. SmRNA pathways intersect with most of the pathways regulating different steps in the life of a messenger RNA (mRNA), starting from transcription and ending at mRNA decay. SmRNAs function in both nuclear and cytoplasmic compartments; the regulation of mRNA stability and translation in the cytoplasm and the epigenetic regulation of gene expression in the nucleus are the main and best-known modes of smRNA action. However, recent evidence from animal systems indicates that smRNAs and RNA interference (RNAi) also participate in the regulation of alternative pre-mRNA splicing, one of the most crucial steps in the fast, efficient global reprogramming of gene expression required for survival under stress. Emerging evidence from bioinformatics studies indicates that a specific class of plant smRNAs, induced by various abiotic stresses, the sutr-siRNAs, has the potential to target regulatory regions within introns and thus may act in the regulation of splicing in response to stresses. This review summarizes the major types of plant smRNAs in the context of their mechanisms of action and also provides examples of their involvement in regulation of gene expression in response to environmental cues and developmental stresses. In addition, we describe current advances in our understanding of how smRNAs function in the regulation of pre-mRNA splicing. WIREs RNA 2016, 7:356-381. doi: 10.1002/wrna.1340 For further resources related to this article, please visit the WIREs website.

The Pokeweed Leaf mRNA Transcriptome and Its Regulation by Jasmonic Acid

The American pokeweed plant, Phytolacca americana, is recognized for synthesizing pokeweed antiviral protein (PAP), a ribosome inactivating protein (RIP) that inhibits the replication of several plant and animal viruses. The plant is also a heavy metal accumulator with applications in soil remediation. However, little is known about pokeweed stress responses, as large-scale sequencing projects have not been performed for this species. Here, we sequenced the mRNA transcriptome of pokeweed in the presence and absence of jasmonic acid (JA), a hormone mediating plant defense. Trinity-based de novo assembly of mRNA from leaf tissue and BLASTx homology searches against public sequence databases resulted in the annotation of 59 096 transcripts. Differential expression analysis identified JA-responsive genes that may be involved in defense against pathogen infection and herbivory. We confirmed the existence of several PAP isoforms and cloned a potentially novel isoform of PAP. Expression analysis indicated that PAP isoforms are differentially responsive to JA, perhaps indicating specialized roles within the plant. Finally, we identified 52 305 natural antisense transcript pairs, four of which comprised PAP isoforms, suggesting a novel form of RIP gene regulation. This transcriptome-wide study of a Phytolaccaceae family member provides a source of new genes that may be involved in stress tolerance in this plant. The sequences generated in our study have been deposited in the SRA database under project # SRP069141.

Role of the ACL2 locus in flower stalk elongation in Arabidopsis thaliana

The acaulis2 (acl2) mutant of Arabidopsis thaliana shows a defect in flower stalk elongation. We identified the mutation point of acl2 by map-based cloning. The ACL2 locus is located within an approximately 320-kb region at around 100 map units on chromosome 1. One nucleotide substitution was detected in this region in the acl2 mutant, but no significant open reading frames were found around this mutation point. When wild-type DNA fragments containing the mutation point were introduced into acl2 mutant plants, some transgenic plants partially or almost completely recovered from the defect in flower stalk elongation. 3'-RACE experiments showed that bidirectional transcripts containing the acl2 mutation point were expressed, and the Plant MPSS database revealed that several small RNAs were produced from this region. Microarray analysis showed that transcription of many genes is activated in flower stalks of acl2 mutant plants. Overexpression of some of these genes caused a dwarf phenotype in wild-type plants. These results suggest the following novel mechanism for control of the elongation of flower stalks. Bidirectional non-coding RNAs are transcribed from the ACL2 locus, and small RNAs are generated from them in flower stalks. These small RNAs repress the transcription of a set of genes whose expression represses flower stalk elongation, and flower stalks are therefore fully elongated.

Salt-induced transcription factor MYB74 is regulated by the RNA-directed DNA methylation pathway in Arabidopsis

Salt stress is one of the major abiotic stresses in agriculture worldwide that causes crop failure by interfering with the profile of gene expression and cell metabolism. Transcription factors and RNA-directed DNA methylation (RdDM) play an important role in the regulation of gene activation under abiotic stress in plants. This work characterized AtMYB74, a member of the R2R3-MYB gene family, which is transcriptionally regulated mainly by RdDM as a response in salt stress in Arabidopsis. Bisulphite sequencing indicated that 24-nt siRNAs target a region approximately 500bp upstream of the transcription initiation site of AtMYB74, which is heavily methylated. Levels of DNA methylation in this region were significantly reduced in wild type plants under salt stress, whereas no changes were found in RdDM mutants. Northern blot and quantitative real-time reverse transcription PCR analysis showed that the accumulation of 24-nt siRNAs was decreased in WT plants under salt stress. Further promoter deletion analysis revealed that the siRNA target region is essential for maintaining AtMYB74 expression patterns. In addition, transgenic plants overexpressing AtMYB74 displayed hypersensitivity to NaCl during seed germination. These results suggest that changes in the levels of the five 24-nt siRNAs regulate the AtMYB74 transcription factor via RdDM in response to salt stress.

Post-transcriptional gene silencing triggered by sense transgenes involves uncapped antisense RNA and differs from silencing intentionally triggered by antisense transgenes

Although post-transcriptional gene silencing (PTGS) has been studied for more than a decade, there is still a gap in our understanding of how de novo silencing is initiated against genetic elements that are not supposed to produce double-stranded (ds)RNA. Given the pervasive transcription occurring throughout eukaryote genomes, we tested the hypothesis that unintended transcription could produce antisense (as)RNA molecules that participate to the initiation of PTGS triggered by sense transgenes (S-PTGS). Our results reveal a higher level of asRNA in Arabidopsis thaliana lines that spontaneously trigger S-PTGS than in lines that do not. However, PTGS triggered by antisense transgenes (AS-PTGS) differs from S-PTGS. In particular, a hypomorphic ago1 mutation that suppresses S-PTGS prevents the degradation of asRNA but not sense RNA during AS-PTGS, suggesting a different treatment of coding and non-coding RNA by AGO1, likely because of AGO1 association to polysomes. Moreover, the intended asRNA produced during AS-PTGS is capped whereas the asRNA produced during S-PTGS derives from 3′ maturation of a read-through transcript and is uncapped. Thus, we propose that uncapped asRNA corresponds to the aberrant RNA molecule that is converted to dsRNA by RNA-DEPENDENT RNA POLYMERASE 6 in siRNA-bodies to initiate S-PTGS, whereas capped asRNA must anneal with sense RNA to produce dsRNA that initiate AS-PTGS.

Functions of plants long non-coding RNAs

Long non-coding RNAs (lncRNAs) have been emerged as important players for various biological pathways, including dosage compensation, genomic imprinting, chromatin regulation, alternative splicing and nuclear organization. A large number of lncRNAs had already been identified by different approaches in plants, while the functions of only a few of them have been investigated. This review will summarize our current understanding of a wide range of plant lncRNAs functions, and highlight their roles in the regulation of diverse pathways in plants. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Ancient endo-siRNA pathways reveal new tricks

Endogenously produced small interfering RNAs (endo-siRNAs, 18-30 nucleotides) play a key role in gene regulatory pathways, guiding Argonaute effector proteins as a part of a functional ribonucleoprotein complex called the RISC (RNA induced silencing complex) to complementarily target nucleic acid. Enabled by the advent of high throughput sequencing, there has been an explosion in the identification of endo-siRNAs in all three kingdoms of life since the discovery of the first microRNA in 1993. Concurrently, our knowledge of the variety of cellular processes in which small RNA pathways related to RNA interference (RNAi) play key regulatory roles has also expanded dramatically. Building on the strong foundation of RNAi established over the past fifteen years, this review uses a historical context to highlight exciting recent developments in endo-siRNA pathways. Specifically, my focus will be on recent insights regarding the Argonaute effectors, their endo-siRNA guides and the functional outputs of these pathways in several model systems that have been longstanding champions of small RNA research. I will also touch on newly discovered roles for bacterial Argonautes, which have been integral in deciphering Argonaute structure and demonstrate key functions of these conserved pathways in genome defense.

Bioinformatics analysis of small RNAs in pima (Gossypium barbadense L.)

Small RNAs (sRNAs) are ~20 to 24 nucleotide single-stranded RNAs that play crucial roles in regulation of gene expression. In plants, sRNAs are classified into microRNAs (miRNAs), repeat-associated siRNAs (ra-siRNAs), phased siRNAs (pha-siRNAs), cis and trans natural antisense transcript siRNAs (cis- and trans-nat siRNAs). Pima (Gossypium barbadense L.) is one of the most economically important fiber crops, producing the best and longest spinnable fiber. Although some miRNAs are profiled in Pima, little is known about siRNAs, the largest subclass of plant sRNAs. In order to profile these gene regulators in Pima, a comprehensive analysis of sRNAs was conducted by mining publicly available sRNA data, leading to identification of 678 miRNAs, 3,559,126 ra-siRNAs, 627 pha-siRNAs, 136,600 cis-nat siRNAs and 79,994 trans-nat siRNAs. The 678 miRNAs, belonging to 98 conserved and 402 lineage-specific families, were produced from 2,138 precursors, of which 297 arose from introns, exons, or intron/UTR-exon junctions of protein-coding genes. Ra-siRNAs were produced from various repeat loci, while most (97%) were yielded from retrotransposons, especially LTRs (long terminal repeats). The genes encoding auxin-signaling-related proteins, NBS-LRRs and transcription factors were major sources of pha-siRNAs, while two conserved TAS3 homologs were found as well. Most cis-NATs in Pima overlapped in enclosed and convergent orientations, while a few hybridized in divergent and coincided orientations. Most cis- and trans-nat siRNAs were produced from overlapping regions. Additionally, characteristics of length and the 5’-first nucleotide of each sRNA class were analyzed as well. Results in this study created a valuable molecular resource that would facilitate studies on mechanism of controlling gene expression.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0601-9) contains supplementary material, which is available to authorized users.

Genome-wide determination of poly(A) sites in Medicago truncatula: evolutionary conservation of alternative poly(A) site choice

Background

Alternative polyadenylation (APA) plays an important role in the post-transcriptional regulation of gene expression. Little is known about how APA sites may evolve in homologous genes in different plant species. To this end, comparative studies of APA sites in different organisms are needed. In this study, a collection of poly(A) sites in Medicago truncatula, a model system for legume plants, has been generated and compared with APA sites in Arabidopsis thaliana.

Results

The poly(A) tags from a deep-sequencing protocol were mapped to the annotated M. truncatula genome, and the identified poly(A) sites used to update the annotations of 14,203 genes. The results show that 64% of M. truncatula genes possess more than one poly(A) site, comparable to the percentages reported for Arabidopsis and rice. In addition, the poly(A) signals associated with M. truncatula genes were similar to those seen in Arabidopsis and other plants. The 3′-UTR lengths are correlated in pairs of orthologous genes between M. truncatula and Arabidopsis. Very little conservation of intronic poly(A) sites was found between Arabidopsis and M. truncatula, which suggests that such sites are likely to be species-specific in plants. In contrast, there is a greater conservation of CDS-localized poly(A) sites in these two species. A sizeable number of M. truncatula antisense poly(A) sites were found. A high percentage of the associated target genes possess Arabidopsis orthologs that are also associated with antisense sites. This is suggestive of important roles for antisense regulation of these target genes.

Conclusions

Our results reveal some distinct patterns of sense and antisense poly(A) sites in Arabidopsis and M. truncatula. In so doing, this study lends insight into general evolutionary trends of alternative polyadenylation in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-615) contains supplementary material, which is available to authorized users.

The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis

In eukaryotic RNA silencing, RNase-III classes of enzymes in the Dicer family process double-stranded RNA of cellular or exogenous origin into small-RNA (sRNA) molecules. sRNAs are then loaded into effector proteins known as ARGONAUTEs (AGOs), which, as part of RNA-induced silencing complexes, target complementary RNA or DNA for silencing. Plants have evolved a large variety of pathways over the Dicer-AGO consortium, which most likely underpins part of their phenotypic plasticity. Dicer-like proteins produce all known classes of plant silencing sRNAs, which are invariably stabilized via 2'-O-methylation mediated by HUA ENHANCER 1 (HEN1), potentially amplified by the action of several RNA-dependent RNA polymerases, and function through a variety of AGO proteins. Here, we review the known characteristics and biochemical properties of the core silencing factors found in the model plant Arabidopsis thaliana. We also describe how interactions between these core factors and more specialized proteins allow the production of a plethora of silencing sRNAs involved in a large array of biological functions. We emphasize in particular the biogenesis and activities of silencing sRNAs of endogenous origin.

Global analysis of cis-natural antisense transcripts and their heat-responsive nat-siRNAs in Brassica rapa

BACKGROUND

Brassica rapa includes several important leaf vegetable crops whose production is often damaged by high temperature. Cis-natural antisense transcripts (cis-NATs) and cis-NATs-derived small interfering RNAs (nat-siRNAs) play important roles in plant development and stress responses. However, genome-wide cis-NATs in B. rapa are not known. The NATs and nat-siRNAs that respond to heat stress have never been well studied in B. rapa. Here, we took advantage of RNA-seq and small RNA (sRNA) deep sequencing technology to identify cis-NATs and heat responsive nat-siRNAs in B. rapa.

RESULTS

Analyses of four RNA sequencing datasets revealed 1031 cis-NATs B. rapa ssp. chinensis cv Wut and B. rapa ssp. pekinensis cv. Bre. Based on sequence homology between Arabidopsis thaliana and B. rapa, 303 conserved cis-NATs in B. rapa were found to correspond to 280 cis-NATs in Arabidopsis; the remaining 728 novel cis-NATs were identified as Brassica-specific ones. Using six sRNA libraries, 4846 nat-siRNAs derived from 150 cis-NATs were detected. Differential expression analysis revealed that nat-siRNAs derived from 12 cis-NATs were responsive to heat stress, and most of them showed strand bias. Real-time PCR indicated that most of the transcripts generating heat-responsive nat-siRNAs were upregulated under heat stress, while the transcripts from the opposite strands of the same loci were downregulated.

CONCLUSIONS

Our results provide the first subsets of genome-wide cis-NATs and heat-responsive nat-siRNAs in B. rapa; these sRNAs are potentially useful for the genetic improvement of heat tolerance in B. rapa and other crops.

Arabidopsis non-coding RNA regulation in abiotic stress responses

Plant growth and productivity are largely affected by environmental stresses. Therefore, plants have evolved unique adaptation mechanisms to abiotic stresses through fine-tuned adjustment of gene expression and metabolism. Recent advanced technologies, such as genome-wide transcriptome analysis, have revealed that a vast amount of non-coding RNAs (ncRNAs) apart from the well-known housekeeping ncRNAs such as rRNAs, tRNAs, small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs) are expressed under abiotic stress conditions. These various types of ncRNAs are involved in chromatin regulation, modulation of RNA stability and translational repression during abiotic stress response. In this review, we summarize recent progress that has been made on ncRNA research in plant abiotic stress response.

Transcription termination and chimeric RNA formation controlled by Arabidopsis thaliana FPA

Alternative cleavage and polyadenylation influence the coding and regulatory potential of mRNAs and where transcription termination occurs. Although widespread, few regulators of this process are known. The Arabidopsis thaliana protein FPA is a rare example of a trans-acting regulator of poly(A) site choice. Analysing fpa mutants therefore provides an opportunity to reveal generic consequences of disrupting this process. We used direct RNA sequencing to quantify shifts in RNA 3′ formation in fpa mutants. Here we show that specific chimeric RNAs formed between the exons of otherwise separate genes are a striking consequence of loss of FPA function. We define intergenic read-through transcripts resulting from defective RNA 3′ end formation in fpa mutants and detail cryptic splicing and antisense transcription associated with these read-through RNAs. We identify alternative polyadenylation within introns that is sensitive to FPA and show FPA-dependent shifts in IBM1 poly(A) site selection that differ from those recently defined in mutants defective in intragenic heterochromatin and DNA methylation. Finally, we show that defective termination at specific loci in fpa mutants is shared with dicer-like 1 (dcl1) or dcl4 mutants, leading us to develop alternative explanations for some silencing roles of these proteins. We relate our findings to the impact that altered patterns of 3′ end formation can have on gene and genome organisation.

The ends of almost all eukaryotic protein-coding genes are defined by a poly(A) signal. When genes are transcribed into mRNA by RNA polymerase II, the poly(A) signal guides cleavage of the precursor mRNA at a particular site; this is accompanied by the addition of a poly(A) tail to the mRNA and termination of transcription. Many genes have more than one poly(A) signal and the regulated choice of which to select can effectively determine what the gene will code for, how the gene can be regulated and where transcription termination occurs. We discovered a rare example of a regulator of poly(A) site choice, called FPA, while studying flower development in the model plant Arabidopsis thaliana. Studying FPA therefore provides an opportunity to understand not only its roles in plant biology but also the generic consequences of disrupting alternative polyadenylation. In this study, we use a technique called direct RNA sequencing to quantify genome-wide shifts in poly(A) site selection in plants that lack FPA function. One of our most striking findings is that in the absence of FPA we detect chimeric RNAs formed between two otherwise separate and well-characterised genes.