Long noncoding RNAs in the model species Brachypodium distachyon

Genome-Wide Identification and Characterization of Long Non-Coding RNAs Associated with Floral Scent Formation in Jasmine (Jasminum sambac)

Long non-coding RNAs (lncRNAs) have emerged as curial regulators of diverse biological processes in plants. Jasmine (Jasminum sambac) is a world-renowned ornamental plant for its attractive and exceptional flower fragrance. However, to date, no systematic screening of lncRNAs and their regulatory roles in the production of the floral fragrance of jasmine flowers has been reported. In this study, we identified a total of 31,079 novel lncRNAs based on an analysis of strand-specific RNA-Seq data from J. sambac flowers at different stages. The lncRNAs identified in jasmine flowers exhibited distinct characteristics compared with protein-coding genes (PCGs), including lower expression levels, shorter transcript lengths, and fewer exons. Certain jasmine lncRNAs possess detectable sequence conservation with other species. Expression analysis identified 2752 differentially expressed lncRNAs (DE_lncRNAs) and 8002 DE_PCGs in flowers at the full-blooming stage. DE_lncRNAs could potentially cis- and trans-regulate PCGs, among which DE_lincRNAs and their targets showed significant opposite expression patterns. The flowers at the full-blooming stage are specifically enriched with abundant phenylpropanoids and terpenoids potentially contributed by DE_lncRNA cis-regulated PCGs. Notably, we found that many cis-regulated DE_lncRNAs may be involved in terpenoid and phenylpropanoid/benzenoid biosynthesis pathways, which potentially contribute to the production of jasmine floral scents. Our study reports numerous jasmine lncRNAs and identifies floral-scent-biosynthesis-related lncRNAs, which highlights their potential functions in regulating the floral scent formation of jasmine and lays the foundations for future molecular breeding.

ICAnnoLncRNA: A Snakemake Pipeline for a Long Non-Coding-RNA Search and Annotation in Transcriptomic Sequences

Long non-coding RNAs (lncRNAs) are RNA molecules longer than 200 nucleotides that do not encode proteins. Experimental studies have shown the diversity and importance of lncRNA functions in plants. To expand knowledge about lncRNAs in other species, computational pipelines that allow for standardised data-processing steps in a mode that does not require user control up until the final result were actively developed recently. These advancements enable wider functionality for lncRNA data identification and analysis. In the present work, we propose the ICAnnoLncRNA pipeline for the automatic identification, classification and annotation of plant lncRNAs in assembled transcriptomic sequences. It uses the LncFinder software for the identification of lncRNAs and allows the adjustment of recognition parameters using genomic data for which lncRNA annotation is available. The pipeline allows the prediction of lncRNA candidates, alignment of lncRNA sequences to the reference genome, filtering of erroneous/noise transcripts and probable transposable elements, lncRNA classification by genome location, comparison with sequences from external databases and analysis of lncRNA structural features and expression. We used transcriptomic sequences from 15 maize libraries assembled by Trinity and Hisat2/StringTie to demonstrate the application of the ICAnnoLncRNA pipeline.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Non-Coding RNAs in Response to Drought Stress

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

Genome wide identification and characterization of abiotic stress responsive lncRNAs in Capsicum annuum

Long non-coding RNAs (lncRNAs) are a type of non-coding transcripts having length of more than 200 nucleotides lacking protein-coding ability. In the present study, 12807 lncRNAs were identified in Capsicum annuum tissues exposed to abiotic stress conditions viz. heat, cold, osmotic and salinity stress. Expression analysis of lncRNAs in different treatment conditions demonstrates their stress-specific expression. Thirty lncRNAs were found to act as precursors for 10 microRNAs (miRNAs) of C. annuum. Additionally, a total of 1807 lncRNAs were found to interact with 194 miRNAs which targeted 621 mRNAs of C. annuum. Among these, 344 lncRNAs were found to act as target mimics for 621 genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that out of those 621 gene sequences, 546 were tagged with GO terms, 105 Enzyme Code (EC) numbers were assigned to 246 genes and 223 genes are found to be involved in 63 biological pathways. In this report, we have highlighted the prospective role of lncRNAs in different abiotic stress conditions by interacting with miRNAs and regulating stress responsive transcription factors (TFs) such as DoF, WRKY, MYB, bZIP and ERF in C. annuum.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

The regulatory landscape of early maize inflorescence development

BACKGROUND

The functional genome of agronomically important plant species remains largely unexplored, yet presents a virtually untapped resource for targeted crop improvement. Functional elements of regulatory DNA revealed through profiles of chromatin accessibility can be harnessed for fine-tuning gene expression to optimal phenotypes in specific environments.

RESULT

Here, we investigate the non-coding regulatory space in the maize (Zea mays) genome during early reproductive development of pollen- and grain-bearing inflorescences. Using an assay for differential sensitivity of chromatin to micrococcal nuclease (MNase) digestion, we profile accessible chromatin and nucleosome occupancy in these largely undifferentiated tissues and classify at least 1.6% of the genome as accessible, with the majority of MNase hypersensitive sites marking proximal promoters, but also 3' ends of maize genes. This approach maps regulatory elements to footprint-level resolution. Integration of complementary transcriptome profiles and transcription factor occupancy data are used to annotate regulatory factors, such as combinatorial transcription factor binding motifs and long non-coding RNAs, that potentially contribute to organogenesis, including tissue-specific regulation between male and female inflorescence structures. Finally, genome-wide association studies for inflorescence architecture traits based solely on functional regions delineated by MNase hypersensitivity reveals new SNP-trait associations in known regulators of inflorescence development as well as new candidates.

CONCLUSIONS

These analyses provide a comprehensive look into the cis-regulatory landscape during inflorescence differentiation in a major cereal crop, which ultimately shapes architecture and influences yield potential.

Long Non-coding RNA in Plants in the Era of Reference Sequences

The discovery of non-coding RNAs (ncRNAs), and the subsequent elucidation of their functional roles, was largely delayed due to the misidentification of non-protein-coding parts of DNA as "junk DNA," which forced ncRNAs into the shadows of their protein-coding counterparts. However, over the past decade, insight into the important regulatory roles of ncRNAs has led to rapid progress in their identification and characterization. Of the different types of ncRNAs, long non-coding RNAs (lncRNAs), has attracted considerable attention due to their mRNA-like structures and gene regulatory functions in plant stress responses. While RNA sequencing has been commonly used for mining lncRNAs, a lack of widespread conservation at the sequence level in addition to relatively low and highly tissue-specific expression patterns challenges high-throughput in silico identification approaches. The complex folding characteristics of lncRNA molecules also complicate target predictions, as the knowledge about the interaction interfaces between lncRNAs and potential targets is insufficient. Progress in characterizing lncRNAs and their targets from different species may hold the key to efficient identification of this class of ncRNAs from transcriptomic and potentially genomic resources. In wheat and barley, two of the most important crops, the knowledge about lncRNAs is very limited. However, recently published high-quality genomes of these crops are considered as promising resources for the identification of not only lncRNAs, but any class of molecules. Considering the increasing demand for food, these resources should be used efficiently to discover molecular mechanisms lying behind development and a/biotic stress responses. As our understanding of lncRNAs expands, interactions among ncRNA classes, as well as interactions with the coding sequences, will likely define novel functional networks that may be modulated for crop improvement.

Transcriptome and Phytochemical Analyses Provide New Insights Into Long Non-Coding RNAs Modulating Characteristic Secondary Metabolites of Oolong Tea (Camellia sinensis) in Solar-Withering

Oolong tea is a popular and semi-fermented beverage. During the processing of tea leaves, withering is the first indispensable process for improving flavor. However, the roles of long non-coding RNAs (lncRNAs) and the characteristic secondary metabolites during the withering of oolong tea leaves remain unknown. In this study, phytochemical analyses indicated that total polyphenols, flavonoids, catechins, epigallocatechin (EGC), catechin gallate (CG), gallocatechin gallate (GCG), epicatechin gallate (ECG), and epigallocatechin gallate (EGCG) were all less abundant in the solar-withered leaves (SW) than in the fresh leaves (FL) and indoor-withered leaves (IW). In contrast, terpenoid, jasmonic acid (JA), and methyl jasmonate (MeJA) contents were higher in the SW than in the FL and IW. By analyzing the transcriptome data, we detected 32,036 lncRNAs. On the basis of the Kyoto Encyclopedia of Genes and Genomes analysis, the flavonoid metabolic pathway, the terpenoid metabolic pathway, and the JA/MeJA biosynthesis and signal transduction pathway were enriched pathways. Additionally, 63 differentially expressed lncRNAs (DE-lncRNAs) and 23 target genes were identified related to the three pathways. A comparison of the expression profiles of the DE-lncRNAs and their target genes between the SW and IW revealed four up-regulated genes (FLS, CCR, CAD, and HCT), seven up-regulated lncRNAs, four down-regulated genes (4CL, CHI, F3H, and F3'H), and three down-regulated lncRNAs related to flavonoid metabolism; nine up-regulated genes (DXS, CMK, HDS, HDR, AACT, MVK, PMK, GGPPS, and TPS), three up-regulated lncRNAs, and six down-regulated lncRNAs related to terpenoid metabolism; as well as six up-regulated genes (LOX, AOS, AOC, OPR, ACX, and MFP2), four up-regulated lncRNAs, and three down-regulated lncRNAs related to JA/MeJA biosynthesis and signal transduction. These results suggested that the expression of DE-lncRNAs and their targets involved in the three pathways may be related to the low abundance of the total polyphenols, flavonoids, and catechins (EGC, CG, GCG, ECG, and EGCG) and the high abundance of terpenoids in the SW. Moreover, solar irradiation, high JA and MeJA contents, and the endogenous target mimic (eTM)-related regulatory mechanism in the SW were also crucial for increasing the terpenoid levels. These findings provide new insights into the greater contribution of solar-withering to the high-quality flavor of oolong tea compared with the effects of indoor-withering.

Screening and characterization of long noncoding RNAs involved in the albinism of Ananas comosus var. bracteatus leaves

Long noncoding RNAs (lncRNAs) have been reported to play key regulatory roles in plant growth, development, and biotic and abiotic stress physiology. Revealing the mechanism of lncRNA regulation in the albino portions of leaves is important for understanding the development of chimeric leaves in Ananas comosus var. bracteatus. In this study, a total of 3,543 candidate lncRNAs were identified, among which 1,451 were differentially expressed between completely green (CGr) and completely white (CWh) leaves. LncRNAs tend to have shorter transcripts, lower expression levels, and greater expression specificity than protein-coding genes. Predicted lncRNA targets were functionally annotated by the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A lncRNA-mRNA interaction network was constructed, and 36 target mRNAs related to chlorophyll metabolism were predicted to interact with 86 lncRNAs. Among these, 25 significantly differentially expressed lncRNAs putatively interacted with 16 target mRNAs. Based on an expression pattern analysis of the lncRNAs and their target mRNAs, the lncRNAs targeting magnesium chelatase subunit H (ChlH), protochlorophyllide oxidoreductase (POR), and heme o synthase (COX10) were suggested as key regulators of chlorophyll metabolism. This study provides the first lncRNA database for A. comosus var. bracteatus and contributes greatly to understanding the mechanism of epigenetic regulation of leaf albinism.

Prediction of regulatory long intergenic non-coding RNAs acting in trans through base-pairing interactions

Background

Long intergenic non-coding RNAs (lincRNAs) can act as regulators of expression of protein-coding genes. Trans-natural antisense transcripts (trans-NATs) are a type of lincRNAs that contain sequence complementary to mRNA from other loci. The regulatory potential of trans-NATs has been poorly studied in eukaryotes and no example of trans-NATs regulating gene expression in plants are reported. The goal of this study was to identify lincRNAs, and particularly trans-NATs, in Arabidopsis thaliana that have a potential to regulate expression of target genes in trans at the transcriptional or translational level.

Results

We identified 1001 lincRNAs using an RNAseq dataset from total polyA⁺ and polysome-associated RNA of seedlings grown under high and low phosphate, or shoots and roots treated with different phytohormones, of which 550 were differentially regulated. Approximately 30% of lincRNAs showed conservation amongst Brassicaceae and 25% harbored transposon element (TE) sequences. Gene co-expression network analysis highlighted a group of lincRNAs associated with the response of roots to low phosphate. A total of 129 trans-NATs were predicted, of which 88 were significantly differentially expressed under at least one pairwise comparison. Five trans-NATs showed a positive correlation between their expression and target mRNA steady-state levels, and three showed a negative correlation. Expression of four trans-NATs positively correlated with a change in target mRNA polysome association. The regulatory potential of these trans-NATs did not implicate miRNA mimics nor siRNAs. We also looked for lincRNAs that could regulate gene expression in trans by Watson-Crick DNA:RNA base pairing with target protein-encoding loci. We identified 100 and 81 with a positive or negative correlation, respectively, with steady-state level of their predicted target. The regulatory potential of one such candidate lincRNA harboring a SINE TE sequence was validated in a protoplast assay on three distinct genes containing homologous TE sequence in their promoters. Construction of networks highlighted other putative lincRNAs with multiple predicted target loci for which expression was positively correlated with target gene expression.

Conclusions

This study identified lincRNAs in Arabidopsis with potential in regulating target gene expression in trans by both RNA:RNA and RNA:DNA base pairing and highlights lincRNAs harboring TE sequences in such activity.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5946-0) contains supplementary material, which is available to authorized users.

Transcriptome Analysis Suggests the Roles of Long Intergenic Non-coding RNAs in the Growth Performance of Weaned Piglets

Long intergenic non-coding RNAs (lincRNAs) have been considered to play a key regulatory role in various biological processes. An increasing number of studies have utilized transcriptome analysis to obtain lincRNAs with functions related to cancer, but lincRNAs affecting growth rates in weaned piglets are rarely described. Although lincRNAs have been systematically identified in various mouse tissues and cell lines, studies of lincRNA in pigs remain rare. Therefore, identifying and characterizing novel lincRNAs affecting the growth performance of weaned piglets is of great importance. Here, we reconstructed 101,988 lincRNA transcripts and identified 1,078 lincRNAs in two groups of longissimus dorsi muscle (LDM) and subcutaneous fat (SF) based on published RNA-seq datasets. These lincRNAs exhibit typical characteristics, such as shorter lengths and lower expression relative to protein-encoding genes. Gene ontology analysis revealed that some lincRNAs could be involved in weaned piglet related processes, such as insulin resistance and the AMPK signaling pathway. We also compared the positional relationship between differentially expressed lincRNAs (DELs) and quantitative trait loci (QTL) and found that some of DELs may play an important role in piglet growth and development. Our work details part of the lincRNAs that may affect the growth performance of weaned piglets and promotes future studies of lincRNAs for molecular-assisted development in weaned piglets.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans)

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis–antenna proteins, circadian rhythm–plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.

Transposon-Derived Non-coding RNAs and Their Function in Plants

Transposable elements (TEs) are often regarded as harmful genomic factors and indeed they are strongly suppressed by the epigenetic silencing mechanisms. On the other hand, the mobilization of TEs brings about variability of genome and transcriptome which are essential in the survival and evolution of the host species. The vast majority of such controlling TEs influence the neighboring genes in cis by either promoting or repressing the transcriptional activities. Although TEs are highly repetitive in the genomes and transcribed in specific stress conditions or developmental stages, the trans-acting regulatory roles of TE-derived RNAs have been rarely studied. It was only recently that TEs were investigated for their regulatory roles as a form of RNA. Particularly in plants, TEs are ample source of small RNAs such as small interfering (si) RNAs and micro (mi) RNAs. Those TE-derived small RNAs have potentials to affect non-TE transcripts by sequence complementarity, thereby generating novel gene regulatory networks including stress resistance and hybridization barrier. Apart from the small RNAs, a number of long non-coding RNAs (lncRNAs) are originated from TEs in plants. For example, a retrotransposon-derived lncRNA expressed in rice root acts as a decoy RNA or miRNA target mimic which negatively controls miRNA171. The post-transcriptional suppression of miRNA171 in roots ensures the stabilization of the target transcripts encoding SCARECROW-LIKE transcription factors, the key regulators of root development. In this review article, the recent discoveries of the regulatory roles of TE-derived RNAs in plants will be highlighted.