Genome-wide prediction and identification of cis-natural antisense transcripts in Arabidopsis thaliana

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

Recent insights into the functions and mechanisms of antisense RNA: emerging applications in cancer therapy and precision medicine

The antisense RNA molecule is a unique DNA transcript consisting of 19-23 nucleotides, characterized by its complementary nature to mRNA. These antisense RNAs play a crucial role in regulating gene expression at various stages, including replication, transcription, and translation. Additionally, artificial antisense RNAs have demonstrated their ability to effectively modulate gene expression in host cells. Consequently, there has been a substantial increase in research dedicated to investigating the roles of antisense RNAs. These molecules have been found to be influential in various cellular processes, such as X-chromosome inactivation and imprinted silencing in healthy cells. However, it is important to recognize that in cancer cells; aberrantly expressed antisense RNAs can trigger the epigenetic silencing of tumor suppressor genes. Moreover, the presence of deletion-induced aberrant antisense RNAs can lead to the development of diseases through epigenetic silencing. One area of drug development worth mentioning is antisense oligonucleotides (ASOs), and a prime example of an oncogenic trans-acting long noncoding RNA (lncRNA) is HOTAIR (HOX transcript antisense RNA). NATs (noncoding antisense transcripts) are dysregulated in many cancers, and researchers are just beginning to unravel their roles as crucial regulators of cancer's hallmarks, as well as their potential for cancer therapy. In this review, we summarize the emerging roles and mechanisms of antisense RNA and explore their application in cancer therapy.

Global identification of natural antisense transcripts in Gossypium hirsutum and Gossypium barbadense under chilling stress

Natural antisense transcripts (NATs) in model plants have been recognized as important regulators of gene expression under abiotic stresses. However, the functional roles of NATs in crops under low temperature are still unclear. Here, we identified 815 and 689 NATs from leaves of Gossypium hirsutum and G. barbadense under chilling stress. Among those, 224 NATs were identified as interspecific homologs between the two species. The correlation coefficients for expression of NATs and their cognate sense genes (CSG) were 0.43 and 0.37 in G. hirsutum and G. barbadense, respectively. Furthermore, expression of interspecific NATs and CSGs alike was highly consistent under chilling stress with correlation coefficients of 0.90-0.91. Four cold-associated NATs were selected for functional validation using virus-induced gene silencing (VIGS). Our results suggest that CAN1 engage in the molecular regulation of chilling stress by regulating SnRK2.8 expression. This highly conserved NAT have valuable potential for applications in breeding cold-tolerant cotton.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Comprehensive investigation of long non-coding RNAs in an endophytic fungus Calcarisporium arbuscula NRRL 3705

Long non-coding RNAs (lncRNAs) play an important role in eukaryotic cells. However, there is no report of lncRNAs in endophytic fungi Calcarisporium arbuscula. Here, in Calcarisporium arbuscula NRRL 3705, an endophytic fungus predominantly producing mycotoxins aurovertins, the genome-wide identification of lncRNAs was carried out based on RNA-Seq. Totally, 1332 lncRNAs were identified, including 1082 long intergenic noncoding RNAs, 64 long intronic noncoding RNAs and 186 long noncoding natural antisense transcripts. The average length of lncRNA and mRNA were 254 and 1102 bp, respectively. LncRNAs were shorter, with fewer exons and lower expression levels. Moreover, there were 39 up-regulated lncRNAs and 10 down-regulated lncRNAs in the ΔaurA mutant, which lacks the aurovertin biosynthetic enzyme AurA. Interestingly, expression of genes related to the metabolism of linoleic acid and methane were significantly down regulated in the ΔaurA mutant. This study enriches the endophytic fungal lncRNA database and provide a basis for further research.

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.

Genome-wide identification of antisense lncRNAs and their association with susceptibility to Flavobacterium psychrophilum in rainbow trout

Eukaryotic genomes encode long noncoding natural antisense transcripts (lncNATs) that have been increasingly recognized as regulatory members of gene expression. Recently, we identified a few antisense transcripts correlating in expression with immune-related genes. However, a systematic genome-wide analysis of lncNATs in rainbow trout is lacking. This study used 134 RNA-Seq datasets from five different projects to identify antisense transcripts. A total of 13,503 lncNATs were identified genome-wide. About 75% of lncNATs showed multiple exons compared to 36.5% of the intergenic lncRNAs. RNA-Seq datasets from resistant, control, and susceptible rainbow trout genetic lines with significant differences in survival rate following Flavobacterium psychrophilum (Fp) infection were analyzed to investigate the potential role of the lncNATs during infection. Twenty-four pairwise comparisons between the different genetic lines, infectious status, and time points revealed 581 differentially expressed (DE) lncNATs and 179 differentially used exons (DUEs). Most of the DE lncNATs strongly and positively correlated in expression with their corresponding sense transcripts across 24 RNA-Seq datasets. LncNATs complementary to genes related to immunity, muscle contraction, proteolysis, and iron/heme metabolism were DE following infection. LncNATs complementary to hemolysis-related genes were DE in the resistant fish compared to susceptible fish on day 5 post-infection, suggesting enhanced clearance of free hemoglobin (Hb) and heme and increased erythropoiesis. LncNATs complementary to hepcidin, a master negative regulator of the plasma iron concentration, were the most downregulated lncNATs on day 5 of bacterial infection in the resistant fish. Ninety-four DE lncNAT, including five complementary to hepcidin, are located within 26 QTL regions previously identified in association with bacterial cold water disease (BCWD) in rainbow trout. Collectively, lncNATs are involved in the molecular architecture of fish immunity and should be further investigated for potential applications in genomic selection and genetic manipulation in aquaculture.

Application of RNAi technology: a novel approach to navigate abiotic stresses

BACKGROUND

Due to the rising population globally, and the demand for food, it is critical to significantly increase crop production by 2050. However, climate change estimates show that droughts and heatwaves will become more prevalent in many parts of the world, posing a severe danger to food output.

METHODS

Selective breeding based on genetic diversity is falling short of meeting the expanding need for food and feed. However, the advent of modern plant genetic engineering, genome editing, and synthetic biology provides precise techniques for producing crops capable of sustaining yield under stress situations.

RESULTS

As a result, crop varieties with built-in genetic tolerance to environmental challenges are desperately needed. In the recent years, small RNA (sRNA) data has progressed to become one of the most effective approaches for the improvement of crops. So many sRNAs (18-30nt) have been found with the use of hi-tech bioinformatics and sequencing techniques which are involved in the regulation of sequence specific gene noncoding RNAs (short ncRNAs) i.e., microRNA (miRNA) and small interfering RNA (siRNA). Such research outcomes may advance our understanding of the genetic basis of adaptability of plants to various environmental challenges and the genetic variation of plant's tolerance to a number of abiotic stresses.

CONCLUSION

The review article highlights current trends and advances in sRNAs' critical role in responses of plants to drought, heat, cold, and salinity, and also the potential technology that identifies the abiotic stress-regulated sRNAs, and techniques for analyzing and validating the target genes.

LncPheDB: a genome-wide lncRNAs regulated phenotypes database in plants

LncPheDB (https://www.lncphedb.com/) is a systematic resource of genome-wide long non-coding RNAs (lncRNAs)-phenotypes associations for multiple species. It was established to display the genome-wide lncRNA annotations, target genes prediction, variant-trait associations, gene-phenotype correlations, lncRNA-phenotype correlations, and the similar non-coding regions of the queried sequence in multiple species. LncPheDB sorted out a total of 203,391 lncRNA sequences, 2000 phenotypes, and 120,271 variants of nine species (Zea mays L., Gossypium barbadense L., Triticum aestivum L., Lycopersicon esculentum Mille, Oryza sativa L., Hordeum vulgare L., Sorghum bicolor L., Glycine max L., and Cucumis sativus L.). By exploring the relationship between lncRNAs and the genomic position of variants in genome-wide association analysis, a total of 68,862 lncRNAs were found to be related to the diversity of agronomic traits. More importantly, to facilitate the study of the functions of lncRNAs, we analyzed the possible target genes of lncRNAs, constructed a blast tool for performing similar fragmentation studies in all species, linked the pages of phenotypic studies related to lncRNAs that possess similar fragments and constructed their regulatory networks. In addition, LncPheDB also provides a user-friendly interface, a genome visualization platform, and multi-level and multi-modal convenient data search engine. We believe that LncPheDB plays a crucial role in mining lncRNA-related plant data.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00084-3.

miR398b and AtC2GnT form a negative feedback loop to regulate Arabidopsis thaliana resistance against Phytophthora parasitica

Oomycetes are diploid eukaryotic microorganisms that seriously threaten sustainable crop production. MicroRNAs (miRNAs) and corresponding natural antisense transcripts (NATs) are important regulators of multiple biological processes. However, little is known about their roles in plant immunity against oomycete pathogens. In this study, we report the identification and functional characterization of miR398b and its cis-NAT, the core-2/I-branching beta-1,6-N-acetylglucosaminyltransferase gene (AtC2GnT), in plant immunity. Gain- and loss-of-function assays revealed that miR398b mediates Arabidopsis thaliana susceptibility to Phytophthora parasitica by targeting Cu/Zn-Superoxidase Dismutase1 (CSD1) and CSD2, leading to suppressed expression of CSD1 and CSD2 and decreased plant disease resistance. We further showed that AtC2GnT transcripts could inhibit the miR398b-CSDs module via inhibition of pri-miR398b expression, leading to elevated plant resistance to P. parasitica. Furthermore, quantitative reverse transcription PCR, RNA ligase-mediated 5'-amplification of cDNA ends (RLM-5' RACE), and transient expression assays indicated that miR398b suppresses the expression of AtC2GnT. We generated AtC2GnT-silenced A. thaliana plants by CRISPR/Cas9 or RNA interference methods, and the Nicotiana benthamiana NbC2GnT-silenced plants by virus-induced gene silencing. Pathogenicity assays showed that the C2GnT-silenced plants were more susceptible, while AtC2GnT-overexpressing plants exhibited elevated resistance to P. parasitica. AtC2GnT encodes a Golgi-localized protein, and transient expression of AtC2GnT enhanced N. benthamiana resistance to Phytophthora pathogens. Taken together, our results revealed a positive role of AtC2GnT and a negative regulatory loop formed by miR398b and AtC2GnT in regulating plant resistance to P. parasitica.

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.

An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development

Background

Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies.

Result

In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses.

Conclusions

This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02620-2.

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.

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.

Making sense of the natural antisense transcript puzzle

In plants, thousands of genes are associated with antisense transcription, which often produces noncoding RNAs. Although widespread, sense-antisense pairs have been implicated in a limited variety of functions in plants and are often thought to form extensive dsRNA stretches triggering gene silencing. In this opinion, we show that evidence does not support gene silencing as a major role for antisense transcription. In fact, it is more likely that antisense transcripts play diverse functions in gene regulation. We propose a general framework for the initial functional dissection of antisense transcripts, suggesting testable hypotheses relying on an experiment-based decision tree. By moving beyond the gene silencing paradigm, we argue that a broad and diverse role for natural antisense transcription will emerge.

Tuning Beforehand: A Foresight on RNA Interference (RNAi) and In Vitro-Derived dsRNAs to Enhance Crop Resilience to Biotic and Abiotic Stresses

Crop yield is severely affected by biotic and abiotic stresses. Plants adapt to these stresses mainly through gene expression reprogramming at the transcriptional and post-transcriptional levels. Recently, the exogenous application of double-stranded RNAs (dsRNAs) and RNA interference (RNAi) technology has emerged as a sustainable and publicly acceptable alternative to genetic transformation, hence, small RNAs (micro-RNAs and small interfering RNAs) have an important role in combating biotic and abiotic stresses in plants. RNAi limits the transcript level by either suppressing transcription (transcriptional gene silencing) or activating sequence-specific RNA degradation (post-transcriptional gene silencing). Using RNAi tools and their respective targets in abiotic stress responses in many crops is well documented. Many miRNAs families are reported in plant tolerance response or adaptation to drought, salinity, and temperature stresses. In biotic stress, the spray-induced gene silencing (SIGS) provides an intelligent method of using dsRNA as a trigger to silence target genes in pests and pathogens without producing side effects such as those caused by chemical pesticides. In this review, we focus on the potential of SIGS as the most recent application of RNAi in agriculture and point out the trends, challenges, and risks of production technologies. Additionally, we provide insights into the potential applications of exogenous RNAi against biotic stresses. We also review the current status of RNAi/miRNA tools and their respective targets on abiotic stress and the most common responsive miRNA families triggered by stress conditions in different crop species.

An antisense noncoding RNA enhances translation via localized structural rearrangements of its cognate mRNA

A large portion of eukaryotic genes are associated with noncoding, natural antisense transcripts (NATs). Despite sharing extensive sequence complementarity with their sense mRNAs, mRNA-NAT pairs elusively often evade dsRNA-cleavage and siRNA-triggered silencing. More surprisingly, some NATs enhance translation of their sense mRNAs by yet unknown mechanism(s). Here, we show that translation enhancement of the rice (Oryza sativa) PHOSPHATE1.2 (PHO1.2) mRNA is enabled by specific structural rearrangements guided by its noncoding antisense RNA (cis-NATpho1.2). Their interaction in vitro revealed no evidence of widespread intermolecular dsRNA formation, but rather specific local changes in nucleotide base pairing, leading to higher flexibility of PHO1.2 mRNA at a key high guanine-cytosine�(GC) regulatory region inhibiting translation, ∼350-nt downstream of the start codon. Sense-antisense RNA interaction increased formation of the 80S complex in PHO1.2, possibly by inducing structural rearrangement within this inhibitory region, thus making this mRNA more accessible to 60S. This work presents a framework for nucleotide resolution studies of functional mRNA-antisense pairs.

Investigation of the putative role of antisense transcripts as regulators of sense transcripts by correlation analysis of sense-antisense pairs in colorectal cancers

Antisense transcription occurs widely more expected than when it was first identified in bacteria in the 1980s. However, the functional relevance of antisense transcripts in transcription remains controversial. Here, we investigated the putative role of antisense transcripts in regulating their corresponding sense transcripts by analyzing changes in correlative relationships between sense-antisense pairs under tumor and normal conditions. A total of 3469 sense-antisense gene pairs (SAGPs) downloaded from BioMart mapped to a list of sense and antisense genes in RNA-seq data derived from 80 paired colorectal cancer (CRC) samples were analyzed. As a result, cancer-related genes were significantly enriched in the significantly correlated SAGPs (SCPs). Differentially expressed genes estimated between normal and tumor conditions were also significantly more enriched in SCPs than in non-SCPs. Interestingly, using differential correlation analysis, we found that tumor samples had a significantly larger density of genes with higher correlation coefficients than normal samples, as verified by various cancer transcriptomes from The Cancer Genome Atlas (TCGA). Moreover, we found that the magnitude of the correlation between SAGPs could distinguish poor prognostic CRCs from good prognostic CRCs, showing that correlation coefficients between the SAGPs of CRCs with a poor prognosis were significantly stronger than CRCs with a good prognosis. Consistent with this finding, the Cox proportion hazards model showed that the survival rates were significantly different between patients with high and low expression of genes in the SCPs. All these results strongly support the idea that antisense transcripts are important regulators of their corresponding sense transcripts.

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.

Plant Long Noncoding RNAs: New Players in the Field of Post-Transcriptional Regulations

The first reference to the "C-value paradox" reported an apparent imbalance between organismal genome size and morphological complexity. Since then, next-generation sequencing has revolutionized genomic research and revealed that eukaryotic transcriptomes contain a large fraction of non-protein-coding components. Eukaryotic genomes are pervasively transcribed and noncoding regions give rise to a plethora of noncoding RNAs with undeniable biological functions. Among them, long noncoding RNAs (lncRNAs) seem to represent a new layer of gene expression regulation, participating in a wide range of molecular mechanisms at the transcriptional and post-transcriptional levels. In addition to their role in epigenetic regulation, plant lncRNAs have been associated with the degradation of complementary RNAs, the regulation of alternative splicing, protein sub-cellular localization, the promotion of translation and protein post-translational modifications. In this review, we report and integrate numerous and complex mechanisms through which long noncoding transcripts regulate post-transcriptional gene expression in plants.

A novel antisense long non-coding RNA participates in asexual and sexual reproduction by regulating the expression of GzmetE in Fusarium graminearum

Fusarium graminearum is an important worldwide pathogen that causes Fusarium head blight in wheat, barley, maize and other grains. LncRNAs play important roles in many biological processes, but little is known about their functions and mechanisms in filamentous fungi. Here, we report that a natural antisense RNA, GzmetE-AS, is transcribed from the opposite strand of GzmetE. GzmetE encodes a homoserine O-acetyltransferase, which is important for sexual development and plant infection. The expression of GzmetE-AS was increased significantly during the conidiation stage, while GzmetE was upregulated in the late stage of sexual reproduction. Overexpression of GzmetE-AS inhibited the transcription of GzmetE. In contrast, the expression of GzmetE was significantly increased in GzmetE-AS transcription termination strain GzmetE-AS-T. Furthermore, GzmetE-AS-T produced more perithecia and facilitated the ascospore discharge, resembling the phenotype of GzmetE overexpressing strains. However, overexpression of GzmetE-AS in ∆dcl1/2 strain cannot inhibit the expression of GzmetE, and the GzmetE nat-siRNA is also significantly reduced in ∆dcl1/2 mutant. Taken together, we have identified a novel antisense lncRNA GzmetE-AS, which is involved in asexual and sexual reproduction by regulating its antisense gene GzmetE through RNAi pathway. Our findings reveal that the lncRNA plays critical roles in the development of F. graminearum.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana

The biological significance of non-coding RNAs (ncRNAs) has been firmly established to be important for the regulation of genes involved in stress acclimation. Light plays an important role for the growth of plants providing the energy for photosynthesis; however, excessive light conditions can also cause substantial defects. Small RNAs (sRNAs) are a class of non-coding RNAs that regulate transcript levels of protein-coding genes and mediate epigenetic silencing. Next generation sequencing facilitates the identification of small non-coding RNA classes such as miRNAs (microRNAs) and small-interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), but changes in the ncRNA transcriptome in response to high light are poorly understood. We subjected Arabidopsis plants to high light conditions and performed a temporal in-depth study of the transcriptome data after 3 h, 6 h, and 2 days of high light treatment. We identified a large number of high light responsive miRNAs and sRNAs derived from NAT gene pairs, lncRNAs and TAS transcripts. We performed target predictions for differentially expressed miRNAs and correlated their expression levels through mRNA sequencing data. GO analysis of the targets revealed an overrepresentation of genes involved in transcriptional regulation. In A. thaliana, sRNA-mediated regulation of gene expression in response to high light treatment is mainly carried out by miRNAs and sRNAs derived from NAT gene pairs, and from lncRNAs. This study provides a deeper understanding of sRNA-dependent regulatory networks in high light acclimation.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Natural antisense transcripts of MIR398 genes suppress microR398 processing and attenuate plant thermotolerance

MicroRNAs (miRNAs) and natural antisense transcripts (NATs) control many biological processes and have been broadly applied for genetic manipulation of eukaryotic gene expression. Still unclear, however, are whether and how NATs regulate miRNA production. Here, we report that the cis-NATs of MIR398 genes repress the processing of their pri-miRNAs. Through genome-wide analysis of RNA sequencing data, we identify cis-NATs of MIRNA genes in Arabidopsis and Brassica. In Arabidopsis, MIR398b and MIR398c are coexpressed in vascular tissues with their antisense genes NAT398b and NAT398c, respectively. Knock down of NAT398b and NAT398c promotes miR398 processing, resulting in stronger plant thermotolerance owing to silencing of miR398-targeted genes; in contrast, their overexpression activates NAT398b and NAT398c, causing poorer thermotolerance due to the upregulation of miR398-targeted genes. Unexpectedly, overexpression of MIR398b and MIR398c activates NAT398b and NAT398c. Taken together, these results suggest that NAT398b/c repress miR398 biogenesis and attenuate plant thermotolerance via a regulatory loop.

MiRNAs and natural antisense transcripts can both regulate gene expression and plant development. Here, the authors show that cis-NATs to MIR398 repress processing of pri-miR398 and that cis-NAT expression is downregulated at high temperatures, contributing to miR398 mediated thermotolerance responses.

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.

Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome

BACKGROUND

Phased small interfering RNA (phasiRNA) is primarily derived from the 22-nt miRNA targeting loci. GhMYB2, a gene with potential roles in cotton fiber cell fate determination, is a target gene of miR828 and miR858 in the generation of phasiRNAs.

RESULTS

In the presented work, through the evaluation of phasing scores and phasiRNA distribution pattern, we found that phasiRNAs from GhMYB2 were derived from the 3' cleavage fragments of 22-nt miR828 and 21-nt miR858 respectively. These two miRNA targeting sites initiated two phasing frames on transcripts of one locus. By means of RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE), we further demonstrated that phasiRNAs derived from the two phasing frames played a role in cis-regulation of GhMYB2. The phasiRNAs derived from GhMYB2 were expressed in the somatic tissues, especially in anther and hypocotyl. We further employed our previous small RNA sequencing data as well as the degradome data of cotton fiber bearing ovules, anthers, hypocotyls and embryogenic calli tissues published in public databases, to validate the expression, phasing pattern and functions of phasiRNAs.

CONCLUSIONS

The presenting research provide insights of the molecular mechanism of phasiRNAs in regulation of GhMYB2 loci.

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.

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.

Control of Cognate Sense mRNA Translation by cis-Natural Antisense RNAs

Cis-Natural Antisense Transcripts (cis-NATs), which overlap protein coding genes and are transcribed from the opposite DNA strand, constitute an important group of noncoding RNAs. Whereas several examples of cis-NATs regulating the expression of their cognate sense gene are known, most cis-NATs function by altering the steady-state level or structure of mRNA via changes in transcription, mRNA stability, or splicing, and very few cases involve the regulation of sense mRNA translation. This study was designed to systematically search for cis-NATs influencing cognate sense mRNA translation in Arabidopsis (Arabidopsis thaliana). Establishment of a pipeline relying on sequencing of total polyA⁺ and polysomal RNA from Arabidopsis grown under various conditions (i.e. nutrient deprivation and phytohormone treatments) allowed the identification of 14 cis-NATs whose expression correlated either positively or negatively with cognate sense mRNA translation. With use of a combination of cis-NAT stable over-expression in transgenic plants and transient expression in protoplasts, the impact of cis-NAT expression on mRNA translation was confirmed for 4 out of 5 tested cis-NAT:sense mRNA pairs. These results expand the number of cis-NATs known to regulate cognate sense mRNA translation and provide a foundation for future studies of their mode of action. Moreover, this study highlights the role of this class of noncoding RNAs in translation regulation.

Revealing the dominant long noncoding RNAs responding to the infection with Colletotrichum gloeosporioides in Hevea brasiliensis

Background

Rubber tree (Hevea brasiliensis) acts as an important tropic economic crop and rubber tree anthracnose, mainly caused by Colletotrichum gloeosporioides, is one of the most common fungal disease, which leads to serious loss of rubber production. Therefore, the investigation on disease resistance is of great worldwide significance. In the past decades, substantial progress has been made on coding gene families related with plant disease resistance. However, in rubber tree, whether the disease resistance mechanism involves noncoding RNAs, especially long noncoding RNAs (lncRNAs), still remains poorly understood.

Results

Here, we modeled the development of H. brasiliensis leaf samples inoculated with C. gloeosporioides at divergent stages, explored to identify the expressed ncRNAs by RNA-seq, and investigated the dominant lncRNAs responding to the infection, through constructing a co-expressed network systematically. On the dominant lncRNAs, we explored the potential functional role of lncRNA11254 recruiting the transcription factor, and that lncRNA11041 and lncRNA11205 probably stimulate the accumulation of corresponding disease responsive miRNAs, and further modulate the expressions of target genes, accompanying with experimental examination.

Conclusions

Take together, computational analyses in silico and experimental evidences in our research collectively revealed the responsive roles of dominant lncRNAs to the pathogen. The results will provide new perspectives to unveil the plant disease resistance mechanisms, and will presumably provide a new theoretical basis and candidate prognostic markers for the optimization and innovation of genetic breeding for rubber tree.

Reviewers

This article was reviewed by Ryan McGinty and Roland Huber.

Electronic supplementary material

The online version of this article (10.1186/s13062-019-0235-z) contains supplementary material, which is available to authorized users.

Recent advances in small RNA mediated plant-virus interactions

Small RNAs (sRNA) are reported to play pivotal roles in the epigenetic and post-transcriptional regulation of gene expression during growth, development, and stress response in plants. Recently, the involvement of two different classes of sRNAs namely, miRNAs (microRNAs), and siRNAs (small interfering RNAs) in biotic stress response has been underlined. Notably, during virus infection, these sRNAs deploy antiviral defense by regulating the gene expression of the modulators of host defense pathways. As a counter defense, viruses have evolved strategic pathways involving the production of suppressors that interfere with the host silencing machinery. This molecular arms race between the sophisticated gene regulatory mechanism of host plants fine-tuned by sRNAs and the defense response exhibited by the virus has gained much attention among the researchers. So far, several reports have been published showing the mechanistic insights on sRNA-regulated defense mechanism in response to virus infection in several crop plants. In this context, our review enumerates the molecular mechanisms underlying host immunity against viruses mediated by sRNAs, the counter defense strategies employed by viruses to surpass this immunogenic response and the advances made in our understanding of plant-virus interactions. Altogether, the report would be insightful for the researchers working to decode the sRNA-mediated defense response in crop plants challenged with virus infection.

Long Non-Coding RNAs: Rising Regulators of Plant Reproductive Development

Long non-coding RNAs (lncRNAs) exert a multitude of functions in regulating numerous biological processes. Recent studies have uncovered a growing number of lncRNAs within the plant genome. These molecules show striking tissue-specific expression patterns, suggesting that they exert regulatory functions in the growth and development processes of plants. Plant reproductive development is tightly regulated by both environmental and endogenous factors. As plant reproductive development is a crucial aspect of crop breeding, lncRNAs that modulate reproductive development are now particularly worth regarding. Here, we summarize findings that implicate lncRNAs in the control of plant reproductive development, especially in flowering control. Additionally, we elaborate on the regulation mechanisms of lncRNAs, tools for research on their function and mechanism, and potential directions of future research.

Mechanisms of Antisense Transcription Initiation with Implications in Gene Expression, Genomic Integrity and Disease Pathogenesis

Non-coding antisense transcripts arise from the strand opposite the sense strand. Over 70% of the human genome generates non-coding antisense transcripts while less than 2% of the genome codes for proteins. Antisense transcripts and/or the act of antisense transcription regulate gene expression and genome integrity by interfering with sense transcription and modulating histone modifications or DNA methylation. Hence, they have significant pathological and physiological relevance. Indeed, antisense transcripts were found to be associated with various diseases including cancer, diabetes, cardiac and neurodegenerative disorders, and, thus, have promising potentials for prognostic and diagnostic markers and therapeutic development. However, it is not clearly understood how antisense transcription is initiated and epigenetically regulated. Such knowledge would provide new insights into the regulation of antisense transcription, and hence disease pathogenesis with therapeutic development. The recent studies on antisense transcription initiation and its epigenetic regulation, which are limited, are discussed here. Furthermore, we concisely describe how antisense transcription/transcripts regulate gene expression and genome integrity with implications in disease pathogenesis and therapeutic development.

LncRNA expression profile and ceRNA analysis in tomato during flowering

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of lncRNAs during the flowering period of tomato are largely unknown. To explore the lncRNA profiles and functions during flowering in tomato, we performed strand-specific paired-end RNA sequencing of tomato leaves, flowers and roots, with three biological replicates. We identified 10919 lncRNAs including 248 novel lncRNAs, of which 65 novel lncRNAs were significantly differentially expressed (DE) in the flowers, leaves, and roots. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out to identify the cis target gene of DE lncRNAs. The results showed that the lncRNAs might play an important role in the growth, development, and apoptosis of flowering tomato plant by regulating the formation of intima in flower tissues, binding to various molecules, influencing metabolic pathways, and inducing apoptosis. Moreover, we identified the interaction between 32, 78, and 397 kinds of miRNAs, lncRNAs, and mRNAs. The results suggest that the lncRNAs can regulate the expression of mRNA during flowering period in tomato by forming competitive endogenous RNA, and further regulate various biological metabolism pathways in tomato.

Biotic Stress-Induced Priming and De-Priming of Transcriptional Memory in Arabidopsis and Apple

Under natural growth conditions, plants experience various and repetitive biotic and abiotic stresses. Salicylic acid (SA) is a key phytohormone involved in the response to biotic challenges. Application of synthetic SA analogues can efficiently prime defense responses, and leads to improved pathogen resistance. Because SA analogues can result in long-term priming and memory, we identified genes for which expression was affected by the SA analogue and explored the role of DNA methylation in this memorization process. We show that treatments with an SA analogue can lead to long-term transcriptional memory of particular genes in Arabidopsis. We found that subsequent challenging of such plants with a bacterial elicitor reverted this transcriptional memory, bringing their expression back to the original pre-treatment level. We also made very similar observations in apple (Malus domestica), suggesting that this expression pattern is highly conserved in plants. Finally, we found a potential role for DNA methylation in the observed transcriptional memory behavior. We show that plants defective in DNA methylation pathways displayed a different memory behavior. Our work improves our understanding of the role of transcriptional memory in priming, and has important implication concerning the application of SA analogues in agricultural settings.

Differential co-expression networks of long non-coding RNAs and mRNAs in Cleistogenes songorica under water stress and during recovery

BACKGROUND

Water stress seriously constrains plant growth and yield. Long non-coding RNAs (lncRNAs) serve as versatile regulators in various biological regulatory processes. To date, the systematic screening and potential functions of lncRNA have not yet been characterized in Cleistogenes songorica, especially under water stress conditions.

RESULTS

In this study, we obtained the root and shoot transcriptomes of young C. songorica plants subjected to different degrees of water stress and recovery treatments by Illumina-based RNA-seq. A total of 3397 lncRNAs were identified through bioinformatics analysis. LncRNA differential expression analysis indicated that the higher response of roots compared to shoots during water stress and recovery. We further identified the 1644 transcription factors, 189 of which were corresponded to 163 lncRNAs in C. songorica. Though comparative analyses with major Poaceae species based on blast, 81 water stress-related orthologues regulated to lncRNAs were identified as a core of evolutionary conserved genes important to regulate water stress responses in the family. Among these target genes, two genes were found to be involved in the abscisic acid (ABA) signalling pathway, and four genes were enriched for starch and sucrose metabolism. Additionally, the 52 lncRNAs were predicted as target mimics for microRNAs (miRNAs) in C. songorica. RT-qPCR results suggested that MSTRG.43964.1 and MSTRG.4400.2 may regulate the expression of miRNA397 and miRNA166, respectively, as target mimics under water stress and during recovery. Finally, a co-expression network was constructed based on the lncRNAs, miRNAs, protein-coding genes (PCgenes) and transcription factors under water stress and during recovery in C. songorica.

CONCLUSIONS

In C. songorica, lncRNAs, miRNAs, PCgenes and transcription factors constitute a complex transcriptional regulatory network which lncRNAs can regulate PCgenes and miRNAs under water stress and recovery. This study provides fundamental resources to deepen our knowledge on lncRNAs during ubiquitous water stress.

Global identification of Arabidopsis lncRNAs reveals the regulation of MAF4 by a natural antisense RNA

Long non-coding RNAs (lncRNAs) have emerged as important regulators of gene expression and plant development. Here, we identified 6,510 lncRNAs in Arabidopsis under normal or stress conditions. We found that the expression of natural antisense transcripts (NATs) that are transcribed in the opposite direction of protein-coding genes often positively correlates with and is required for the expression of their cognate sense genes. We further characterized MAS, a NAT-lncRNA produced from the MADS AFFECTING FLOWERING4 (MAF4) locus. MAS is induced by cold and indispensable for the activation of MAF4 transcription and suppression of precocious flowering. MAS activates MAF4 by interacting with WDR5a, one core component of the COMPASS-like complexes, and recruiting WDR5a to MAF4 to enhance histone 3 lysine 4 trimethylation (H3K4me3). Our study greatly extends the repertoire of lncRNAs in Arabidopsis and reveals a role for NAT-lncRNAs in regulating gene expression in vernalization response and likely in other biological processes.

Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus

Background

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis.

Results

A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels.

Conclusions

We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5117-8) contains supplementary material, which is available to authorized users.

Transcriptome characterization of moso bamboo (Phyllostachys edulis) seedlings in response to exogenous gibberellin applications

BACKGROUND

Moso bamboo (Phyllostachys edulis) is a well-known bamboo species of high economic value in the textile industry due to its rapid growth. Phytohormones, which are master regulators of growth and development, serve as important endogenous signals. However, the mechanisms through which phytohormones regulate growth in moso bamboo remain unknown to date.

RESULTS

Here, we reported that exogenous gibberellins (GA) applications resulted in a significantly increased internode length and lignin condensation. Transcriptome sequencing revealed that photosynthesis-related genes were enriched in the GA-repressed gene class, which was consistent with the decrease in leaf chlorophyll concentrations and the lower rate of photosynthesis following GA treatment. Exogenous GA applications on seedlings are relatively easy to perform, thus we used 4-week-old whole seedlings of bamboo for GA- treatment followed by high throughput sequencing. In this study, we identified 932 cis-nature antisense transcripts (cis-NATs), and 22,196 alternative splicing (AS) events in total. Among them, 42 cis-nature antisense transcripts (cis-NATs) and 442 AS events were differentially expressed upon exposure to exogenous GA3, suggesting that post-transcriptional regulation might be also involved in the GA3 response. Targets of differential expression of cis-NATs included genes involved in hormone receptor, photosynthesis and cell wall biogenesis. For example, LAC4 and its corresponding cis-NATs were GA3-induced, and may be involved in the accumulation of lignin, thus affecting cell wall composition.

CONCLUSIONS

This study provides novel insights illustrating how GA alters post-transcriptional regulation and will shed light on the underlying mechanism of growth modulated by GA in moso bamboo.

Widespread antisense transcription of Populus genome under drought

Antisense transcription is widespread in many genomes and plays important regulatory roles in gene expression. The objective of our study was to investigate the extent and functional relevance of antisense transcription in forest trees. We employed Populus, a model tree species, to probe the antisense transcriptional response of tree genome under drought, through stranded RNA-seq analysis. We detected nearly 48% of annotated Populus gene loci with antisense transcripts and 44% of them with co-transcription from both DNA strands. Global distribution of reads pattern across annotated gene regions uncovered that antisense transcription was enriched in untranslated regions while sense reads were predominantly mapped in coding exons. We further detected 1185 drought-responsive sense and antisense gene loci and identified a strong positive correlation between the expression of antisense and sense transcripts. Additionally, we assessed the antisense expression in introns and found a strong correlation between intronic expression and exonic expression, confirming antisense transcription of introns contributes to transcriptional activity of Populus genome under drought. Finally, we functionally characterized drought-responsive sense-antisense transcript pairs through gene ontology analysis and discovered that functional groups including transcription factors and histones were concordantly regulated at both sense and antisense transcriptional level. Overall, our study demonstrated the extensive occurrence of antisense transcripts of Populus genes under drought and provided insights into genome structure, regulation pattern and functional significance of drought-responsive antisense genes in forest trees. Datasets generated in this study serve as a foundation for future genetic analysis to improve our understanding of gene regulation by antisense transcription.

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.

Knockdown of collagen α-1(III) inhibits glioma cell proliferation and migration and is regulated by miR128-3p

As a member of the collagen family, collagen α-1(III) (COL3A1) is an important protein in the development and progression of several tumors. However, the role of COL3A1 in glioma is not yet clear. The present study examined the expression and function of COL3A1 in glioma cell behavior and identified microRNA (miRNA) regulators. It was demonstrated that COL3A1 expression was upregulated in glioma and directly correlated with the tumor grade. Analysis of the GSE4290 and GSE7696 profiles acquired from the Gene Expression Omnibus database also revealed an increased COL3A1 expression in malignant gliomas compared with the lower grade gliomas and non-tumor brain tissue, which was directly correlated with glioma grade. To explore the functional role of COL3A1 in glioma cell growth, small interfering RNA interference was applied to inhibit COL3A1 expression in Hs683 and U251 cells. The relative COL3A1 mRNA and protein expression levels were significantly reduced in the knockdown cells as determined by western blot analysis. In addition, decreased COL3A1 expression in Hs683 and U251 glioma cells resulted in a delay in cell growth and colony disruption as determined by MTS and colony formation assays. Wound healing analysis indicated that cells with suppressed expression of COL3A1 had a reduced ability to migrate. COL3A1 mRNA levels were inversely correlated with the miR128-3p level in glioma, suggesting that miR128-3p expression is associated with COL3A1 inhibition as verified by reverse transcription-quantified polymerase chain reaction. These results suggest that COL3A1 may be a novel regulator of glioblastoma cell behavior and may represent a novel target for gene therapies against glioma.

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.

A novel antisense long noncoding RNA, TWISTED LEAF, maintains leaf blade flattening by regulating its associated sense R2R3-MYB gene in rice

Natural antisense long noncoding RNAs (lncRNAs) are widespread in many organisms. However, their biological functions remain largely unknown, particularly in plants. We report the identification and characterization of an endogenous lncRNA, TWISTED LEAF (TL), which is transcribed from the opposite strand of the R2R3 MYB transcription factor gene locus, OsMYB60, in rice (Oryza sativa). TL and OsMYB60 were found to be coexpressed in many different tissues, and the expression level of TL was higher than that of OsMYB60. Downregulation of TL by RNA interference (RNAi) and overexpression of OsMYB60 resulted in twisted leaf blades in transgenic rice. The expression level of OsMYB60 was significantly increased in TL-RNAi transgenic plants. This suggests that TL may play a cis-regulatory role on OsMYB60 in leaf morphological development. We also determined that the antisense transcription suppressed the sense gene expression by mediating chromatin modifications. We further discovered that a C2H2 transcription factor, OsZFP7, is an OsMYB60 binding partner and involved in leaf development. Taken together, these findings reveal that the cis-natural antisense lncRNA plays a critical role in maintaining leaf blade flattening in rice. Our study uncovers a regulatory mechanism of lncRNA in plant leaf development.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

Molecular cloning, expression, purification and functional characterization of an antifungal cyclophilin protein from Panax ginseng

Cyclophilins (CyPs), a member of peptidyl-prolyl cis-trans isomerases (PPIases), are ubiquitously distributed in organisms such as bacteria, yeast, plants and animals. CyPs have diverse biological functions, with some exhibiting antifungal and antiviral activities. In this study, Panax ginseng cyclophilin (pgCyP), a novel gene encoding an antifungal protein from Panax ginseng, was cloned, and its protein product was expressed in Escherichia coli, and then fractionated by affinity chromatography. The open reading frame of the pgCyP full-length coding sequence was found to encode a single-domain CyP-like protein of 174 amino residues with a calculated molecular weight of 18.7 kDa. The pGEX system was used to express pgCyP fused to glutathione S-transferase. After affinity purification, the protein showed a strong fungal resistance effect on Phytophthora cactorum. In addition, pgCyP showed high PPIase activity. To the best of our knowledge, the present study is the first successful effort to clone and characterize a CyP-like protein gene from Panax ginseng.

Novel Stress-Inducible Antisense RNAs of Protein-Coding Loci Are Synthesized by RNA-Dependent RNA Polymerase

Our previous study identified approximately 6,000 abiotic stress-responsive noncoding transcripts existing on the antisense strand of protein-coding genes and implied that a type of antisense RNA was synthesized from a sense RNA template by RNA-dependent RNA polymerase (RDR). Expression analyses revealed that the expression of novel abiotic stress-induced antisense RNA on 1,136 gene loci was reduced in the rdr1/2/6 mutants. RNase protection indicated that the RD29A antisense RNA and other RDR1/2/6-dependent antisense RNAs are involved in the formation of dsRNA. The accumulation of stress-inducible antisense RNA was decreased and increased in dcp5 and xrn4, respectively, but not changed in dcl2/3/4, nrpd1a and nrpd1b RNA-seq analyses revealed that the majority of the RDR1/2/6-dependent antisense RNA loci did not overlap with RDR1/2/6-dependent 20-30 nt RNA loci. Additionally, rdr1/2/6 mutants decreased the degradation rate of the sense RNA and exhibited arrested root growth during the recovery stage following a drought stress, whereas dcl2/3/4 mutants did not. Collectively, these results indicate that RDRs have stress-inducible antisense RNA synthesis activity and a novel biological function that is different from the known endogenous small RNA pathways from protein-coding genes. These data reveal a novel mechanism of RNA regulation during abiotic stress response that involves complex RNA degradation pathways.