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

Full-length transcriptome analysis of maize root tips reveals the molecular mechanism of cold stress during the seedling stage

BACKGROUND

As maize originated in tropical or subtropical zones, most maize germplasm is extremely sensitive to low temperatures during the seedling stage. Clarifying the molecular mechanism of cold acclimation would facilitate the breeding of cold tolerant maize varieties, which is one of the major sustainability factors for crop production. To meet this goal, we investigated two maize inbred lines with contrasting levels of cold tolerance at the seedling stage (IL85, a cold tolerant line; B73, a cold sensitive line), and performed full-length transcriptome sequencing on the root tips of seedlings before and after 24 h of cold treatment.

RESULTS

We identified 152,263 transcripts, including 20,993 novel transcripts, and determined per-transcript expression levels. A total of 1,475 transcripts were specifically up-regulated in the cold tolerant line IL85 under cold stress. GO enrichment analysis revealed that 25 transcripts were involved in reactive oxygen species (ROS) metabolic processes and 15 transcripts were related to the response to heat. Eight genes showed specific differential alternative splicing (DAS) in IL85 under cold stress, and were mainly involved in amine metabolism. A total of 1,111 lncRNAs were further identified, 62 of which were up-regulated in IL85 or B73 under cold stress, and their corresponding target genes were enriched in protein phosphorylation.

CONCLUSIONS

These results provide new insights into the molecular mechanism of cold acclimation during the seedling stage in maize, and will facilitate the development of cultivars with improved cold stress tolerance.

Methylation hallmarks on the histone tail as a linker of osmotic stress and gene transcription

Plants dynamically manipulate their gene expression in acclimation to the challenging environment. Hereinto, the histone methylation tunes the gene transcription via modulation of the chromatin accessibility to transcription machinery. Osmotic stress, which is caused by water deprivation or high concentration of ions, can trigger remarkable changes in histone methylation landscape and genome-wide reprogramming of transcription. However, the dynamic regulation of genes, especially how stress-inducible genes are timely epi-regulated by histone methylation remains largely unclear. In this review, recent findings on the interaction between histone (de)methylation and osmotic stress were summarized, with emphasis on the effects on histone methylation profiles imposed by stress and how histone methylation works to optimize the performance of plants under stress.

The Arabidopsis DREAM complex antagonizes WDR5A to modulate histone H3K4me2/3 deposition for a subset of genome repression

The master transcriptional repressor DREAM (dimerization partner, RB-like, E2F and multivulval class B) complex regulates the cell cycle in eukaryotes, but much remains unknown about how it transmits repressive signals on chromatin to the primary transcriptional machinery (e.g., RNA polymerase II [Pol II]). Through a forward genetic screen, we identified BTE1 (barrier of transcription elongation 1), a plant-specific component of the DREAM complex. The subsequent characterization demonstrated that DREAM complex containing BTE1 antagonizes the activity of Complex Proteins Associated with Set1 (COMPASS)-like complex to repress H3K4me3 occupancy and inhibits Pol II elongation at DREAM target genes. We showed that BTE1 is recruited to chromatin at the promoter-proximal regions of target genes by E2F transcription factors. DREAM target genes exhibit characteristic enrichment of H2A.Z and H3K4me2 modification on chromatin. We further showed that BTE1 directly interacts with WDR5A, a core component of COMPASS-like complex, repressing WDR5A chromatin binding and the elongation of transcription on DREAM target genes. H3K4me3 is known to correlate with the Pol II transcription activation and promotes efficient elongation. Thus, our study illustrates a transcriptional repression mechanism by which the DREAM complex dampens H3K4me3 deposition at a set of genes through its interaction with WDR5A.

Full-Length Transcriptome Sequencing Reveals Alternative Splicing and lncRNA Regulation during Nodule Development in Glycine max

Alternative splicing (AS) is a ubiquitous phenomenon among eukaryotic intron-containing genes, which greatly contributes to transcriptome and proteome diversity. Here we performed the isoform sequencing (Iso-Seq) of soybean underground tissues inoculated and uninoculated with Rhizobium and obtained 200,681 full-length transcripts covering 26,183 gene loci. It was found that 80.78% of the multi-exon loci produced more than one splicing variant. Comprehensive analysis of these identified 7874 differentially splicing events with highly diverse splicing patterns during nodule development, especially in defense and transport-related processes. We further profiled genes with differential isoform usage and revealed that 2008 multi-isoform loci underwent stage-specific or simultaneous major isoform switches after Rhizobium inoculation, indicating that AS is a vital way to regulate nodule development. Moreover, we took the lead in identifying 1563 high-confidence long non-coding RNAs (lncRNAs) in soybean, and 157 of them are differentially expressed during nodule development. Therefore, our study uncovers the landscape of AS during the soybean-Rhizobium interaction and provides systematic transcriptomic data for future study of multiple novel directions in soybean.

Genome-wide identification and characterization of long noncoding RNAs during peach (Prunus persica) fruit development and ripening

LncRNAs represent a class of RNA transcripts of more than 200 nucleotides (nt) in length without discernible protein-coding potential. The expression levels of lncRNAs are significantly affected by stress or developmental cues. Recent studies have shown that lncRNAs participate in fruit development and ripening processes in tomato and strawberry; however, in other fleshy fruits, the association between lncRNAs and fruit ripening remains largely elusive. Here, we constructed 9 ssRNA-Seq libraries from three different peach (Prunus persica) fruit developmental stages comprising the first and second exponential stages and the fruit-ripening stage. In total, 1500 confident lncRNAs from 887 loci were obtained according to the bioinformatics analysis. The lncRNAs identified in peach fruits showed distinct characteristics compared with protein-coding mRNAs, including lower expression levels, lower complexity of alternative splicing, shorter isoforms and smaller numbers of exons. Expression analysis identified 575 differentially expressed lncRNAs (DELs) classified into 6 clusters, among which members of Clusters 1, 2, 4 and 5 were putatively associated with fruit development and ripening processes. Quantitative real-time PCR revealed that the DELs indeed had stage-specific expression patterns in peach fruits. GO and KEGG enrichment analysis revealed that DELs might be associated with fruit-ripening-related physiological and metabolic changes, such as flavonoid biosynthesis, fruit texture softening, chlorophyll breakdown and aroma compound accumulation. Finally, the similarity analysis of lncRNAs within different plant species indicated the low sequence conservation of lncRNAs. Our study reports a large number of fruit-expressed lncRNAs and identifies fruit development phase-specific expressed lncRNA members, which highlights their potential functions in fruit development and ripening processes and lays the foundations for future functional research.

LncRNA TMPO-AS1 promotes esophageal squamous cell carcinoma progression by forming biomolecular condensates with FUS and p300 to regulate TMPO transcription

Esophageal squamous cell carcinoma (ESCC) is one of the most life- and health-threatening malignant diseases worldwide, especially in China. Long noncoding RNAs (lncRNAs) have emerged as important regulators of tumorigenesis and tumor progression. However, the roles and mechanisms of lncRNAs in ESCC require further exploration. Here, in combination with a small interfering RNA (siRNA) library targeting specific lncRNAs, we performed MTS and Transwell assays to screen functional lncRNAs that were overexpressed in ESCC. TMPO-AS1 expression was significantly upregulated in ESCC tumor samples, with higher TMPO-AS1 expression positively correlated with shorter overall survival times. In vitro and in vivo functional experiments revealed that TMPO-AS1 promotes the proliferation and metastasis of ESCC cells. Mechanistically, TMPO-AS1 bound to fused in sarcoma (FUS) and recruited p300 to the TMPO promoter, forming biomolecular condensates in situ to activate TMPO transcription in cis by increasing the acetylation of histone H3 lysine 27 (H3K27ac). Targeting TMPO-AS1 led to impaired ESCC tumor growth in a patient-derived xenograft (PDX) model. We found that TMPO-AS1 is required for cell proliferation and metastasis in ESCC by promoting the expression of TMPO, and both TMPO-AS1 and TMPO might be potential biomarkers and therapeutic targets in ESCC.

The role of a regulatory RNA in promoting esophageal squamous cell carcinoma (ESCC) has been clarified, revealing molecular details that might help in cancer diagnosis and treatment. Xiao-Jing Luo and colleagues at Sun Yat-sen University in China found that overproduction of an RNA molecule called thymopoietin-antisense RNA 1 (TMPO-AS1) in ESCC tissue samples from cancer patients was associated with shorter survival times. Overproduction of this RNA promoted proliferation and spread (metastasis) of the cancer cells. Research on details of the molecular mechanisms involved showed that the RNA ultimately activated the gene that codes for the protein hormone thymopoietin, which has previously been linked with various cancers. The authors suggest that TMPO-AS1 and thymopoietin could serve as diagnostic biomarkers of cancer and become targets for anti-cancer drugs.

Revealing the novel complexity of plant long non-coding RNA by strand-specific and whole transcriptome sequencing for evolutionarily representative plant species

Background

Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts’ direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique.

Results

A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes.

Conclusion

The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08602-9.

Nitrilases NIT1/2/3 Positively Regulate Flowering by Inhibiting MAF4 Expression in Arabidopsis

Three of the nitrilases (NITs), NIT1, NIT2, and NIT3, are ubiquitously existing in plant kingdom, which catalyze indole-3-acetonitrile into the most important auxin indole-3-acetic acid. Auxin is an indispensable hormone, which plays the important roles in almost all processes of plant growth and development. However, there are few reports on the regulation of flowering-time mediated by auxin. Here, we found that in Arabidopsis, nit1/2/3 showed a late flowering phenotype in short days. To explore the molecular mechanism by which NIT1/2/3 regulate flowering time, we performed transcriptome sequencing of nit1/2/3. The results showed that the expression of a MADS-box transcription factor gene MADS AFFECTING FLOWERING4 (MAF4) was dramatically increased in nit1/2/3 comparing to wild type (WT). MAF4 is one of the paralogs of the potent flowering inhibitor FLOWERING LOCUS C (FLC). There are four other paralogs in FLC clade in Arabidopsis, including FLOWERING LOCUS M (FLM/MAF1), MAF2, MAF3, and MAF5. The late flowering phenotype of nit1/2/3 could not be observed in the maf4 background, indicating that the phenotype was specifically dependent on MAF4 rather than other FLC clade members. Interestingly, the expression of a lncRNA gene MAS, which is transcribed in the opposite direction of MAF4, was found significantly increased in nit1/2/3. Also, MAS has been reported to activate MAF4 transcription by promoting histone 3 lysine 4 trimethylation (H3K4me3). As expected, H3K4me3 deposition at MAF4 locus in nit1/2/3 was highly enriched and significantly higher than that of WT. In summary, we show that NITs, NIT1/2/3, positively regulate flowering by repressing MAF4 through manipulating H3K4me3 modification. Further study needs to be performed to explore the largely unknown mechanisms behind it.

Long non-coding RNA transcriptome landscape of anthers at different developmental stages in response to drought stress in tomato

Pollen development is particularly susceptible to drought stress. Long non-coding RNAs (lncRNAs) are known to play a role in plant development and responses to drought stress. However, the expression profile and putative function of lncRNAs in drought-induced male sterility remain largely unknown. In this study, we investigated the lncRNA transcriptome landscape of tomato anthers at early and late stages-tetrad-vacuolated microspore (TED-VUM) and binucleate-mature pollen (BIN-MP) anthers, respectively-in response to drought stress using RNA-sequencing. In total, we identified 67,770 lncRNAs, of which 3053 lncRNAs were drought responsive. Interestingly, there were more differentially expressed (DE) lncRNAs in TED-VUM (2879) than in BIN-MP (174) anthers, which was consistent with the TED-VUM anthers being more drought sensitive. Functional enrichment analysis revealed that the target genes of DE lncRNAs were significantly enriched in diverse metabolic processes, including in carbohydrate metabolism and hormone synthesis. Co-expression analysis also identified 1407 lncRNAs that strongly co-expressed with 8 target genes that are involved in hormone (abscisic acid and jasmonic acid) and carbohydrate (sucrose and starch) metabolisms and tapetum development, highlighting the potential of lncRNA-target-gene modulation of anther development under drought stress. Our results serve as a baseline for future investigations of the potential function of lncRNAs in plant reproductive development under drought stress.

A long noncoding RNA functions in high-light-induced anthocyanin accumulation in apple by activating ethylene synthesis

Anthocyanin production in apple (Malus domestica) fruit and their consequent coloration can be induced by high-light treatment. The hormone ethylene is also essential for this coloration, but the regulatory relationships that link ethylene and light with anthocyanin-associated coloration are not well defined. In this study, we observed that high-light treatment of apple fruit increased anthocyanin accumulation more than moderate-light treatment did and was the main contributor of induced ethylene production and activation of anthocyanin biosynthesis. A transcriptome study of light-treated apple fruit suggested that a long noncoding RNA (lncRNA), MdLNC610, the corresponding gene of which is physically located downstream from the 1-aminocyclopropane-1-carboxylate oxygenase (ACO) ethylene biosynthesis gene MdACO1, likely affects anthocyanin biosynthesis under high-light treatment. Expression and promoter β-glucuronidase reporter analyses further showed that MdLNC610 upregulates expression of MdACO1 and so likely participates in high-light-induced ethylene biosynthesis. Overexpression of MdACO1 and MdLNC610 in apple fruit and calli indicated that a major increase in MdLNC610 expression activates MdACO1 expression, thereby causing an increase in ethylene production and anthocyanin levels. These results suggest that MdLNC610 participates in the regulation of high-light-induced anthocyanin production by functioning as a positive regulator to promote MdACO1 expression and ethylene biosynthesis. Our study provides insights into the relationship between mRNA and lncRNA networks in the ethylene biosynthetic pathway and anthocyanin accumulation in apple fruit.

Role of lncRNAs in cis- and trans-regulatory responses to salt in Populus trichocarpa

Long non-coding RNAs (lncRNAs) are emerging as versatile regulators in diverse biological processes. However, little is known about their cis- and trans-regulatory contributions in gene expression under salt stress. Using 27 RNA-seq data sets from Populus trichocarpa leaves, stems and roots, we identified 2988 high-confidence lncRNAs, including 1183 salt-induced differentially expressed lncRNAs. Among them, 301 lncRNAs have potential for positively affecting their neighboring genes, predominantly in a cis-regulatory manner rather than by co-transcription. Additionally, a co-expression network identified six striking salt-associated modules with a total of 5639 genes, including 426 lncRNAs, and in these lncRNA sequences, the DNA/RNA binding motifs are enriched. This suggests that lncRNAs might contribute to distant gene expression of the salt-associated modules in a trans-regulatory manner. Moreover, we found 30 lncRNAs that have potential to simultaneously cis- and trans-regulate salt-responsive homologous genes, and Ptlinc-NAC72, significantly induced under long-term salt stress, was selected for validating its regulation of the expression and functional roles of the homologs PtNAC72.A and PtNAC72.B (PtNAC72.A/B). The transient transformation of Ptlinc-NAC72 and a dual-luciferase assay of Ptlinc-NAC72 and PtNAC72.A/B promoters confirmed that Ptlinc-NAC72 can directly upregulate PtNAC72.A/B expression, and a presence/absence assay was further conducted to show that the regulation is probably mediated by Ptlinc-NAC72 recognizing the tandem elements (GAAAAA) in the PtNAC72.A/B 5' untranslated region (5'-UTR). Finally, the overexpression of Ptlinc-NAC72 produces a hypersensitive phenotype under salt stress. Altogether, our results shed light on the cis- and trans-regulation of gene expression by lncRNAs in Populus and provides an example of long-term salt-induced Ptlinc-NAC72 that could be used to mitigate growth costs by conferring plant resilience to salt stress.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Reshuffling of the ancestral core-eudicot genome shaped chromatin topology and epigenetic modification in Panax

All extant core-eudicot plants share a common ancestral genome that has experienced cyclic polyploidizations and (re)diploidizations. Reshuffling of the ancestral core-eudicot genome generates abundant genomic diversity, but the role of this diversity in shaping the hierarchical genome architecture, such as chromatin topology and gene expression, remains poorly understood. Here, we assemble chromosome-level genomes of one diploid and three tetraploid Panax species and conduct in-depth comparative genomic and epigenomic analyses. We show that chromosomal interactions within each duplicated ancestral chromosome largely maintain in extant Panax species, albeit experiencing ca. 100–150 million years of evolution from a shared ancestor. Biased genetic fractionation and epigenetic regulation divergence during polyploidization/(re)diploidization processes generate remarkable biochemical diversity of secondary metabolites in the Panax genus. Our study provides a paleo-polyploidization perspective of how reshuffling of the ancestral core-eudicot genome leads to a highly dynamic genome and to the metabolic diversification of extant eudicot plants.

The role of polyploidization generated genomic diversity in shaping the hierarchical genome architecture remains unclear. Here, the authors show that repatterning of the ancestral eudicot genome has resulted in multi-dimensional genome plasticity and secondary metabolite diversification via comparisons of Panax genomes.

Transcriptomic identification of long noncoding RNAs and their hormone-associated nearby coding genes involved in the differential development of caryopses localized on different branches in rice

Long noncoding RNAs (lncRNAs) play important regulatory roles in caryopsis development and grain size in rice. However, whether there exist differences in lncRNA expression between caryopses located on primary branches (CPB) and caryopses located on secondary branches (CSB) that contribute to their differential development remains elusive. Here, we performed transcriptome-wide analysis to identify 2,273 lncRNAs expressed in CPB and CSB at 0, 5, 12, and 20 days after flowering (DAF). Although these lncRNAs were widely distributed, the majority were located in intergenic regions of the 12 rice chromosomes. Based on gene expression cluster analysis, lncRNAs expressed in CPB and CSB were clustered into two subtypes in a position-independent manner: one includes 0- and 5-DAF CPB and CSB, and 12-DAF CSB; the second includes 12-DAF CPB and 20-DAF CPB and CSB. Furthermore, according to the expression value of each lncRNA, K-means cluster analysis revealed 135 early-stage, 116 middle-stage, and 114 late-stage expression-delayed lncRNAs in CSB. Then, we analyzed the expression values of the expression-delayed lncRNAs and nearby coding genes (100 kb upstream and downstream of the lncRNAs), and found 631 lncRNA-mRNA pairs, including 258 lncRNAs and 571 nearby coding genes, some of which are related to hormone-regulated grain development. These results suggested that expression-delayed lncRNAs in CSB may regulate the development of CPB and CSB, providing insight into the mechanism underlying the developmental differences between CPB and CSB, and the differences in grain yield.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Global identification of full-length cassava lncRNAs unveils the role of cold-responsive intergenic lncRNA 1 in cold stress response

Long noncoding RNAs (lncRNAs) have been considered to be important regulators of gene expression in a range of biological processes in plants. A large number of lncRNAs have been identified in plants. However, most of their biological functions still remain to be determined. Here, we identified a total of 3004 lncRNAs in cassava under normal or cold-treated conditions from Iso-seq data. We further characterized a cold-responsive intergenic lncRNA 1 (CRIR1) as a novel positive regulator of the plant response to cold stress. CRIR1 can be significantly induced by cold treatment. Ectopic expression of CRIR1 in cassava enhanced the cold tolerance of transgenic plants. Transcriptome analysis demonstrated that CRIR1 regulated a range of cold stress-related genes in a CBF-independent pathway. We further found that CRIR1 RNA can interact with cassava cold shock protein 5 (MeCSP5), which acts as an RNA chaperone, indicating that CRIR1 may recruit MeCSP5 to improve the translation efficiency of messenger RNA. In summary, our study extends the repertoire of lncRNAs in plants as well as their role in cold stress responses. Moreover, it reveals a mechanism by which CRIR1 affected cold stress response by modulating the expression of stress-responsive genes and increasing their translational yield.

The Identification of MATE Antisense Transcripts in Soybean Using Strand-Specific RNA-Seq Datasets

Natural antisense transcripts (NATs) have been generally reported as negative regulators of their sense counterparts. Multidrug and toxic compound extrusion (MATE) proteins mediate the transport of various substrates. Although MATEs have been identified genome-wide in various plant species, their transcript regulators remain unclear. Here, using the publicly available strand-specific RNA-seq datasets of Glycine soja (wild soybean) which have the data from various tissues including developing pods, developing seeds, embryos, cotyledons and hypocotyls, roots, apical buds, stems, and flowers, we identified 35 antisense transcripts of MATEs from 28 gene loci after transcriptome assembly. Spearman correlation coefficients suggested the positive expression correlations of eight MATE antisense and sense transcript pairs. By aligning the identified transcripts with the reference genome of Glycine max (cultivated soybean), the MATE antisense and sense transcript pairs were identified. Using soybean C08 (Glycine max), in developing pods and seeds, the positive correlations between MATE antisense and sense transcript pairs were shown by RT-qPCR. These findings suggest that soybean antisense transcripts are not necessarily negative transcription regulators of their sense counterparts. This study enhances the existing knowledge on the transcription regulation of MATE transporters by uncovering the previously unknown MATE antisense transcripts and their potential synergetic effects on sense transcripts.

Characterization of XR_311113.2 as a MicroRNA Sponge for Pre-ovulatory Ovarian Follicles of Goats via Long Noncoding RNA Profile and Bioinformatics Analysis

Long noncoding RNAs (lncRNAs) were identified recently as a large class of noncoding RNAs (ncRNAs) with a length ≥200 base pairs (bp). The function and mechanism of lncRNAs have been reported in a growing number of species and tissues. In contrast, the regulatory mechanism of lncRNAs in the goat reproductive system has rarely been reported. In the present study, we sequenced and analyzed the lncRNAs using bioinformatics to identify their expression profiles. As a result, 895 lncRNAs were predicted in the pre-ovulatory ovarian follicles of goats. Eighty-eight lncRNAs were differentially expressed in the Macheng black goat when compared with Boer goat. In addition, the lncRNA XR_311113.2 acted as a sponge of chi-miR-424-5p, as assessed via a luciferase activity assay. Taken together, our findings demonstrate that lncRNAs have potential effects in the ovarian follicles of goats and may represent a promising new research field to understand follicular development.

Genome-wide analysis and functional annotation of chromatin-enriched noncoding RNAs in rice during somatic cell regeneration

BACKGROUND

Plants have the remarkable ability to generate callus, a pluripotent cell mass that acquires competence for subsequent tissue regeneration. Global chromatin remodeling is required for this cell fate transition, but how the process is regulated is not fully understood. Chromatin-enriched noncoding RNAs (cheRNAs) are thought to play important roles in maintaining chromatin state. However, whether cheRNAs participate in somatic cell regeneration in plants has not yet been clarified.

RESULTS

To uncover the characteristics and functions of cheRNAs during somatic cell reprogramming in plants, we systematically investigate cheRNAs during callus induction, proliferation and regeneration in rice. We identify 2284 cheRNAs, most of which are novel long non-coding RNAs or small nucleolar RNAs. These cheRNAs, which are highly conserved across plant species, shuttle between chromatin and the nucleoplasm during somatic cell regeneration. They positively regulate the expression of neighboring genes via specific RNA motifs, which may interact with DNA motifs around cheRNA loci. Large-scale mutant analysis shows that cheRNAs are associated with plant size and seed morphology. Further detailed functional investigation of two che-lncRNAs demonstrates that their loss of function impairs cell dedifferentiation and plant regeneration, highlighting the functions of cheRNAs in regulating the expression of neighboring genes via specific motifs. These findings support cis- regulatory roles of cheRNAs in influencing a variety of rice traits.

CONCLUSIONS

cheRNAs are a distinct subclass of regulatory non-coding RNAs that are required for somatic cell regeneration and regulate rice traits. Targeting cheRNAs has great potential for crop trait improvement and breeding in future.

Novel insights into water-deficit-responsive mRNAs and lncRNAs during fiber development in Gossypium hirsutum

BACKGROUND

The fiber yield and quality of cotton are greatly and periodically affected by water deficit. However, the molecular mechanism of the water deficit response in cotton fiber cells has not been fully elucidated.

RESULTS

In this study, water deficit caused a significant reduction in fiber length, strength, and elongation rate but a dramatic increase in micronaire value. To explore genome-wide transcriptional changes, fibers from cotton plants subjected to water deficit (WD) and normal irrigation (NI) during fiber development were analyzed by transcriptome sequencing. Analysis showed that 3427 mRNAs and 1021 long noncoding RNAs (lncRNAs) from fibers were differentially expressed between WD and NI plants. The maximum number of differentially expressed genes (DEGs) and lncRNAs (DERs) was identified in fibers at the secondary cell wall biosynthesis stage, suggesting that this is a critical period in response to water deficit. Twelve genes in cotton fiber were differentially and persistently expressed at ≥ five time points, suggesting that these genes are involved in both fiber development and the water-deficit response and could potentially be used in breeding to improve cotton resistance to drought stress. A total of 540 DEGs were predicted to be potentially regulated by DERs by analysis of coexpression and genomic colocation, accounting for approximately 15.76% of all DEGs. Four DERs, potentially acting as target mimics for microRNAs (miRNAs), indirectly regulated their corresponding DEGs in response to water deficit.

CONCLUSIONS

This work provides a comprehensive transcriptome analysis of fiber cells and a set of protein-coding genes and lncRNAs implicated in the cotton response to water deficit, significantly affecting fiber quality during the fiber development stage.

Transcriptional Association between mRNAs and Their Paired Natural Antisense Transcripts Following Fusarium oxysporum Inoculation in Brassica rapa L

Long noncoding RNAs (lncRNAs) play important roles in abiotic and biotic stress responses; however, studies on the mechanism of regulation of lncRNA expression are limited in plants. The present study examined the relationship between lncRNA expression level and two active histone modifications (H3K4me3 and H3K36me3) in Brassica rapa. Both histone marks were enriched in the chromatin regions encoding lncRNAs, especially around the transcription start site. The transcription level of long intergenic noncoding RNAs was positively associated with the level of H3K4me3 and H3K36me3, while this association was not observed in natural antisense RNAs (NATs) and intronic noncoding RNAs. As coordinate expression of mRNAs and paired NATs under biotic stress treatment has been identified, the transcriptional relationship between mRNAs and their paired NATs following Fusarium oxysporum f. sp. conglutinans (Foc) inoculation was examined. A positive association of expression levels between mRNAs and their paired NATs following Foc inoculation was observed. This association held for several defense-response-related genes and their NAT pairs. These results suggest that coordinate expression of mRNAs and paired NATs plays a role in the defense response against Foc.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Identification of key gene networks controlling vernalization development characteristics of Isatis indigotica by full-length transcriptomes and gene expression profiles

Isatis indigotica Fort., as a common Chinese medicinal raw material, will lose its medicinal value if it blooms early, so it is highly valuable to clarify the induction mechanism of the vernalization of I. indigotica at low temperature. In this study, the concentrations of soluble sugar, proline, glutathione and zeatin in two germplasms of I. indigotica with different degrees of low temperature tolerance (Y1 and Y2) were determined at 10 days, 20 days and 30 days of low-temperature treatment, and the full-length transcriptome of 24 samples was sequenced by Nanopore sequencing with Oxford Nanopore Technologies (ONT). After that, the data of transcripts involved in the vernalization of I. indigotica at low temperature were obtained, and these transcripts were identified using weighted gene co-expression network analysis (WGCNA). The results revealed the massive accumulation of soluble sugar and proline in Y1 and Y2 after low temperature induction. A total of 18,385 new transcripts, 6168 transcription factors and 470 lncRNAs were obtained. Differential expression analysis showed that gibberellin, flavonoids, fatty acids and some processes related to low temperature response were significantly enriched. Eight key transcripts were identified by WGCNA, among which ONT.14640.1, ONT.9119.1, ONT.13080.2 and ONT.16007.1 encodes a flavonoid transporter, 9-cis-epoxycarotenoid dioxygenase 3 (NCED3), growth factor gene and L-aspartate oxidase in plants, respectively. It indicated that secondary metabolites such as hormones and flavonoids play an important role in the vernalization of I. indigotica. qRT-PCR proved the reliability of transcriptome results. These results provide important insights on the low-temperature vernalization of I. indigotica, and provide a research basis for analyzing the vernalization mechanism of I. indigotica.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01110-2.

A Comparative Transcriptional Landscape of Two Castor Cultivars Obtained by Single-Molecule Sequencing Comparative Analysis

Background and Objectives: Castor (Ricinus communis L.) is an important non-edible oilseed crop. Lm-type female strains and normal amphiprotic strains are important castor cultivars, and are mainly different in their inflorescence structures and leaf shapes. To better understand the mechanisms underlying these differences at the molecular level, we performed a comparative transcriptional analysis.

Materials and Methods: Full-length transcriptome sequencing and short-read RNA sequencing were employed.

Results: A total of 76,068 and 44,223 non-redundant transcripts were obtained from high-quality transcripts of Lm-type female strains and normal amphiprotic strains, respectively. In Lm-type female strains and normal amphiprotic strains, 51,613 and 20,152 alternative splicing events were found, respectively. There were 13,239 transcription factors identified from the full-length transcriptomes. Comparative analysis showed a great variety of gene expression of common and unique transcription factors between the two cultivars. Meanwhile, a functional analysis of the isoforms was conducted. The full-length sequences were used as a reference genome, and a short-read RNA sequencing analysis was performed to conduct differential gene analysis. Furthermore, the function of DEGs were performed to annotation analysis.

Conclusion: The results revealed considerable differences and expression diversity between the two cultivars, well beyond what was reported in previous studies and likely reflecting the differences in architecture between these two cultivars.

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.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Global profiling of RNA-chromatin interactions reveals co-regulatory gene expression networks in Arabidopsis

It is increasingly evident that various RNAs can bind chromatin to regulate gene expression and genome organization. Here we adapted a sequencing-based technique to profile RNA-chromatin interactions at a genome-wide scale in Arabidopsis seedlings. We identified more than 10,000 RNA-chromatin interactions mediated by protein-coding RNAs and non-coding RNAs. Cis and intra-chromosomal interactions are mainly mediated by protein-coding RNAs, whereas inter-chromosomal interactions are primarily mediated by non-coding RNAs. Many RNA-chromatin interactions tend to positively correlate with DNA-DNA interactions, suggesting their mutual influence and reinforcement. We further show that some RNA-chromatin interactions undergo alterations in response to biotic and abiotic stresses and that altered RNA-chromatin interactions form co-regulatory networks. Our study provides a global view on RNA-chromatin interactions in Arabidopsis and a rich resource for future investigations of regulatory roles of RNAs in gene expression and genome organization.

Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton

Long non-coding RNAs (lncRNAs) play important roles in response to biotic and abiotic stress through acting as competing endogenous RNAs (ceRNAs) to decoy mature miRNAs. However, whether this mechanism is involved in cotton salt stress response remains unknown. We report the characterization of an endogenous lncRNA, lncRNA354, whose expression was reduced in salt-treated cotton and was localized at the nucleus and cytoplasm. Using endogenous target mimic (eTM) analysis, we predicted that lncRNA354 had a potential binding site for miR160b. Transient expression in tobacco demonstrated that lncRNA354 was a miR160b eTM and attenuated miR160b suppression of its target genes, including auxin response factors (ARFs). Silencing or overexpressing lncRNA354 affected the expression of miR160b and target ARFs. Silencing lncRNA354 and targets GhARF17/18 resulted in taller cotton plants and enhanced the resistant to salt stress. Overexpression of lncRNA354 and targets GhARF17/18 in Arabidopsis led to dwarf plants, decreased root dry weight and reduced salt tolerance. Our results show that the lncRNA354-miR160b effect on GhARF17/18 expression may modulate auxin signalling and thus affect growth. These results also shed new light on a mechanism of lncRNA-associated responses to salt stress.

Integrated Analysis of mRNA and Non-coding RNA Transcriptome in Pepper (Capsicum chinense) Hybrid at Seedling and Flowering Stages

Pepper is an important vegetable in the world. In this work, mRNA and ncRNA transcriptome profiles were applied to understand the heterosis effect on the alteration in the gene expression at the seedling and flowering stages between the hybrid and its parents in Capsicum chinense. Our phenotypic data indicated that the hybrid has dominance in leaf area, plant scope, plant height, and fruit-related traits. Kyoto Encyclopedia of Genes and Genomes analysis showed that nine members of the plant hormone signal transduction pathway were upregulated in the seedling and flowering stages of the hybrid, which was supported by weighted gene coexpression network analysis and that BC332_23046 (auxin response factor 8), BC332_18317 (auxin-responsive protein IAA20), BC332_13398 (ethylene-responsive transcription factor), and BC332_27606 (ethylene-responsive transcription factor WIN1) were candidate hub genes, suggesting the important potential role of the plant hormone signal transduction in pepper heterosis. Furthermore, some transcription factor families, including bHLH, MYB, and HSF were greatly over-dominant. We also identified 2,525 long ncRNAs (lncRNAs), 47 micro RNAs (miRNAs), and 71 circle RNAs (circRNAs) in the hybrid. In particular, downregulation of miR156, miR169, and miR369 in the hybrid suggested their relationship with pepper growth vigor. Moreover, we constructed some lncRNA–miRNA–mRNA regulatory networks that showed a multi-dimension to understand the ncRNA relationship with heterosis. These results will provide guidance for a better understanding of the molecular mechanism involved in pepper heterosis.

Integrative lncRNA landscape reveals lncRNA-coding gene networks in the secondary cell wall biosynthesis pathway of moso bamboo (Phyllostachys edulis)

Background

LncRNAs are extensively involved in plant biological processes. However, the lack of a comprehensive lncRNA landscape in moso bamboo has hindered the molecular study of lncRNAs. Moreover, the role of lncRNAs in secondary cell wall (SCW) biosynthesis of moso bamboo is elusive.

Results

For comprehensively identifying lncRNA throughout moso bamboo genome, we collected 231 RNA-Seq datasets, 1 Iso-Seq dataset, and 1 full-length cDNA dataset. We used a machine learning approach to improve the pipeline of lncRNA identification and functional annotation based on previous studies and identified 37,009 lncRNAs in moso bamboo. Then, we established a network of potential lncRNA-coding gene for SCW biosynthesis and identified SCW-related lncRNAs. We also proposed that a mechanism exists in bamboo to direct phenylpropanoid intermediates to lignin or flavonoids biosynthesis through the PAL/4CL/C4H genes. In addition, we identified 4 flavonoids and 1 lignin-preferred genes in the PAL/4CL/C4H gene families, which gained implications in molecular breeding.

Conclusions

We provided a comprehensive landscape of lncRNAs in moso bamboo. Through analyses, we identified SCW-related lncRNAs and improved our understanding of lignin and flavonoids biosynthesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07953-z.

SMAD4-induced knockdown of the antisense long noncoding RNA BRE-AS contributes to granulosa cell apoptosis

Antisense long noncoding RNAs (AS-lncRNAs), a sub-class of lncRNAs, are transcribed in the opposite direction from their overlapping protein-coding genes and are implicated in various physiological and pathological processes. However, their role in female reproduction remains largely unknown. Here, we report that BRE-AS, an AS-lncRNA transcript from intron 10 of the protein-coding gene BRE, is involved in granulosa cell (GC) apoptosis. Based on our previous RNA sequencing data, we identified 28 AS-lncRNAs as important in the initiation of porcine follicular atresia, with BRE-AS showing the most significant upregulation in early atretic follicles. In this study, gain- and loss-of-function assays demonstrated that BRE-AS induces early apoptosis in GCs. Mechanistically, BRE-AS acts in cis to suppress the expression of BRE, an anti-apoptotic factor, via direct interaction with the pre-mRNA transcript of the latter, inducing increased GC apoptosis. Notably, we also found that BRE-AS was upregulated in SMAD4-silenced GCs. SMAD4 was identified as a transcriptional repressor of BRE-AS because it inhibits BRE-AS expression and BRE-AS-mediated GC apoptosis. In conclusion, we not only identified a novel AS-lncRNA related to the early apoptosis of GCs and initiation of follicular atresia but also described a novel regulatory pathway, SMAD4/BRE-AS/BRE, coordinating GC function and female fertility.

Functions of long non-coding RNA in Arabidopsis thaliana

A major part of the eukaryotic genome is transcribed into non-coding RNAs (ncRNAs) having no protein coding potential. ncRNAs which are longer than 200 nucleotides are categorized as long non coding RNAs (lncRNAs). Most lncRNAs are induced as a consequence of various environmental and developmental cues. Among plants, the functions of lncRNAs are best studied in Arabidopsis thaliana. In this review, we highlight the important functional roles of various lncRNAs during different stages of Arabidopsis life cycle and their response to environmental changes. These lncRNAs primarily govern processes such as flowering, seed germination, stress response, light- and auxin-regulated development, and RNA-dependent DNA methylation (RdDM). Major challenge is to differentiate between functional and cryptic transcripts. Genome editing, large scale RNAi and computational approaches may help to identify and characterize novel functional lncRNAs in Arabidopsis.

The Natural Antisense Transcript DONE40 Derived from the lncRNA ENOD40 Locus Interacts with SET Domain Protein ASHR3 During Inception of Symbiosis in Arachis hypogaea

The long noncoding RNA ENOD40 is required for cortical cell division during root nodule symbiosis (RNS) of legumes, though it is not essential for actinorhizal RNS. Our objective was to understand whether ENOD40 was required for aeschynomenoid nodule formation in Arachis hypogaea. AhENOD40 express from chromosome 5 (chr5) (AhENOD40-1) and chr15 (AhENOD40-2) during symbiosis, and RNA interference of these transcripts drastically affected nodulation, indicating the importance of ENOD40 in A. hypogaea. Furthermore, we demonstrated several distinct characteristics of ENOD40. (i) Natural antisense transcript (NAT) of ENOD40 was detected from the AhENOD40-1 locus (designated as NAT-AhDONE40). (ii) Both AhENOD40-1 and AhENOD40-2 had two exons, whereas NAT-AhDONE40 was monoexonic. Reverse-transcription quantitative PCR analysis indicated both sense and antisense transcripts to be present in both cytoplasm and nucleus, and their expression increased with the progress of symbiosis. (iii) RNA pull-down from whole cell extracts of infected roots at 4 days postinfection indicated NAT-AhDONE40 to interact with the SET (Su(var)3-9, enhancer of Zeste and Trithorax) domain containing absent small homeotic disc (ASH) family protein AhASHR3 and this interaction was further validated using RNA immunoprecipitation and electrophoretic mobility shift assay. (iv) Chromatin immunoprecipitation assays indicate deposition of ASHR3-specific histone marks H3K36me3 and H3K4me3 in both of the ENOD40 loci during the progress of symbiosis. ASHR3 is known for its role in optimizing cell proliferation and reprogramming. Because both ASHR3 and ENOD40 from legumes cluster away from those in actinorhizal plants and other nonlegumes in phylogenetic distance trees, we hypothesize that the interaction of DONE40 with ASHR3 could have evolved for adapting the nodule organogenesis program for legumes.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Neofunctionalization of a polyploidization-activated cotton long intergenic non-coding RNA DAN1 during drought stress regulation

The genomic shock of whole-genome duplication (WGD) and hybridization introduces great variation into transcriptomes, for both coding and noncoding genes. An altered transcriptome provides a molecular basis for improving adaptation during the evolution of new species. The allotetraploid cotton, together with the putative diploid ancestor species compose a fine model for study the rapid gene neofunctionalization over the genome shock. Here we report on Drought-Associated Non-coding gene 1 (DAN1), a long intergenic noncoding RNA (lincRNA) that arose from the cotton progenitor A-diploid genome after hybridization and WGD events during cotton evolution. DAN1 in allotetraploid upland cotton (Gossypium hirsutum) is a drought-responsive lincRNA predominantly expressed in the nucleoplasm. Chromatin isolation by RNA purification profiling and electrophoretic mobility shift assay analysis demonstrated that GhDAN1 RNA can bind with DNA fragments containing AAAG motifs, similar to DNA binding with one zinc finger transcription factor binding sequences. The suppression of GhDAN1 mainly regulates genes with AAAG motifs in auxin-response pathways, which are associated with drought stress regulation. As a result, GhDAN1-silenced plants exhibit improved tolerance to drought stress. This phenotype resembles the drought-tolerant phenotype of the A-diploid cotton ancestor species, which has an undetectable expression of DAN1. The role of DAN1 in cotton evolution and drought tolerance regulation suggests that the genomic shock of interspecific hybridization and WGD stimulated neofunctionalization of non-coding genes during the natural evolutionary process.

Plant long non-coding RNAs in the regulation of transcription

Eukaryotic genomes are pervasively transcribed, producing large numbers of non-coding RNAs (ncRNAs), including tens of thousands of long ncRNAs (lncRNAs), defined as ncRNAs longer than 200 nucleotides. Recent studies have revealed the important roles lncRNAs play in the regulation of gene expression at various levels in all eukaryotes; moreover, emerging research in plants has identified roles for lncRNAs in key processes such as flowering time control, root organogenesis, reproduction, and adaptation to environmental changes. LncRNAs participate in regulating most steps of gene expression, including reshaping nuclear organization and chromatin structure; governing multiple steps of transcription, splicing, mRNA stability, and translation; and affecting post-translational protein modifications. In this review, I present the latest progress on the lncRNA-mediated regulatory mechanisms modulating transcription in Arabidopsis thaliana, focusing on their functions in regulation of gene expression via chromatin structure and interactions with the transcriptional machinery.

Long Noncoding RNAs in Plants

Plants have an extraordinary diversity of transcription machineries, including five nuclear DNA-dependent RNA polymerases. Four of these enzymes are dedicated to the production of long noncoding RNAs (lncRNAs), which are ribonucleic acids with functions independent of their protein-coding potential. lncRNAs display a broad range of lengths and structures, but they are distinct from the small RNA guides of RNA interference (RNAi) pathways. lncRNAs frequently serve as structural, catalytic, or regulatory molecules for gene expression. They can affect all elements of genes, including promoters, untranslated regions, exons, introns, and terminators, controlling gene expression at various levels, including modifying chromatin accessibility, transcription, splicing, and translation. Certain lncRNAs protect genome integrity, while others respond to environmental cues like temperature, drought, nutrients, and pathogens. In this review, we explain the challenge of defining lncRNAs, introduce the machineries responsible for their production, and organize this knowledge by viewing the functions of lncRNAs throughout the structure of a typical plant gene.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting soybean genotypes subjected to phosphate starvation

BACKGROUND

Phosphorus (P) is essential for plant growth and development, and low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Long noncoding RNAs (lncRNAs) have recently been reported to be key regulators in the responses of plants to stress conditions, but the mechanism through which LP stress mediates the biogenesis of lncRNAs in soybean remains unclear.

RESULTS

In this study, to explore the response mechanisms of lncRNAs to LP stress, we used the roots of two representative soybean genotypes that present opposite responses to P deficiency, namely, a P-sensitive genotype (Bogao) and a P-tolerant genotype (NN94156), for the construction of RNA sequencing (RNA-seq) libraries. In total, 4,166 novel lncRNAs, including 525 differentially expressed (DE) lncRNAs, were identified from the two genotypes at different P levels. GO and KEGG analyses indicated that numerous DE lncRNAs might be involved in diverse biological processes related to phosphate, such as lipid metabolic processes, catalytic activity, cell membrane formation, signal transduction, and nitrogen fixation. Moreover, lncRNA-mRNA-miRNA and lncRNA-mRNA networks were constructed, and the results identified several promising lncRNAs that might be highly valuable for further analysis of the mechanism underlying the response of soybean to LP stress.

CONCLUSIONS

These results revealed that LP stress can significantly alter the genome-wide profiles of lncRNAs, particularly those of the P-sensitive genotype Bogao. Our findings increase the understanding of and provide new insights into the function of lncRNAs in the responses of soybean to P stress.

A Survey of Current Resources to Study lncRNA-Protein Interactions

Phenotypes are driven by regulated gene expression, which in turn are mediated by complex interactions between diverse biological molecules. Protein-DNA interactions such as histone and transcription factor binding are well studied, along with RNA-RNA interactions in short RNA silencing of genes. In contrast, lncRNA-protein interaction (LPI) mechanisms are comparatively unknown, likely directed by the difficulties in studying LPI. However, LPI are emerging as key interactions in epigenetic mechanisms, playing a role in development and disease. Their importance is further highlighted by their conservation across kingdoms. Hence, interest in LPI research is increasing. We therefore review the current state of the art in lncRNA-protein interactions. We specifically surveyed recent computational methods and databases which researchers can exploit for LPI investigation. We discovered that algorithm development is heavily reliant on a few generic databases containing curated LPI information. Additionally, these databases house information at gene-level as opposed to transcript-level annotations. We show that early methods predict LPI using molecular docking, have limited scope and are slow, creating a data processing bottleneck. Recently, machine learning has become the strategy of choice in LPI prediction, likely due to the rapid growth in machine learning infrastructure and expertise. While many of these methods have notable limitations, machine learning is expected to be the basis of modern LPI prediction algorithms.

From "Dark Matter" to "Star": Insight Into the Regulation Mechanisms of Plant Functional Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) play a vital role in a variety of biological functions in plant growth and development. In this study, we provided an overview of the molecular mechanisms of lncRNAs in interacting with other biomolecules with an emphasis on those lncRNAs validated only by low-throughput experiments. LncRNAs function through playing multiple roles, including sponger for sequestering RNA or DNA, guider or decoy for recruiting or hijacking transcription factors or peptides, and scaffold for binding with chromatin modification complexes, as well as precursor of microRNAs or small interfering RNAs. These regulatory roles have been validated in several plant species with a comprehensive list of 73 lncRNA-molecule interaction pairs in 16 plant species found so far, suggesting their commonality in the plant kingdom. Such initial findings of a small number of functional plant lncRNAs represent the beginning of what is to come as lncRNAs with unknown functions were found in orders of magnitude more than proteins.

Bud endodormancy in deciduous fruit trees: advances and prospects

Bud endodormancy is a complex physiological process that is indispensable for the survival, growth, and development of deciduous perennial plants. The timely release of endodormancy is essential for flowering and fruit production of deciduous fruit trees. A better understanding of the mechanism of endodormancy will be of great help in the artificial regulation of endodormancy to cope with climate change and in creating new cultivars with different chilling requirements. Studies in poplar have clarified the mechanism of vegetative bud endodormancy, but the endodormancy of floral buds in fruit trees needs further study. In this review, we focus on the molecular regulation of endodormancy induction, maintenance and release in floral buds of deciduous fruit trees. We also describe recent advances in quantitative trait loci analysis of chilling requirements in fruit trees. We discuss phytohormones, epigenetic regulation, and the detailed molecular network controlling endodormancy, centered on SHORT VEGETATIVE PHASE (SVP) and Dormancy-associated MADS-box (DAM) genes during endodormancy maintenance and release. Combining previous studies and our observations, we propose a regulatory model for bud endodormancy and offer some perspectives for the future.

Interplay of DNA and RNA N 6 -methyladenine with R-loops in regulating gene transcription in Arabidopsis

R-loops and covalent modifications of N ⁶ -methyladenine on DNA (D-6 mA) or RNA (R-m⁶A) have been documented to function in various cellular processes in eukaryotes. However, the relationships between R-loops and both covalent modifications are still elusive in plants. Here, we integrated existing ssDRIP-seq with D-6 mA and R-m⁶A data from Arabidopsis thaliana. We found that the presence of either of both modifications facilitates R-loop formation and transcription of overlapping genes. Interestingly, our study suggests that the presence of R-m⁶A is key to affect R-loop intensity and positively regulate gene transcription. Moreover, the presence of D-6 mA plays an additive role to facilitate the effect of R-m⁶A on R-loop intensity, however, D-6 mA may negatively regulate gene transcription when coexisted with R-m⁶A. Our analyses indicate that D-6 mA, R-m⁶A, or histone marks may act individually and cooperatively with R-loops in regulating gene transcription. Our study is the first to link R-loops with D-6 mA and R-m⁶A in plants, thereby providing new insights into interactions between R-loops with D-6 mA, R-m⁶A, and histone marks for regulating gene transcription.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01010-5.

TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes

BACKGROUND

The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data.

RESULTS

We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: (i) full-length RNA-seq for detection of splicing patterns and (ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts. We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings and Saccharomyces cerevisiae cells as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the most commonly used community gene models, TAIR10 and Araport11 for A.thaliana and SacCer3 for S.cerevisiae. In particular, we identify multiple transient transcripts missing from the existing annotations. Our new annotations promise to improve the quality of A.thaliana and S.cerevisiae genome research.

CONCLUSIONS

Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Genome-wide analysis of long noncoding RNAs affecting floral bud dormancy in pears in response to cold stress

The versatile role of long noncoding RNAs (lncRNAs) in plant growth and development has been established, but a systematic identification and analysis of lncRNAs in the pear has not been reported. Bud dormancy is a crucial and complicated protective mechanism for plants in winter. The roles of lncRNAs in the dormancy process remain largely unclear. In this study, we induced pear floral buds to enter into different dormant statuses by simulating four different chilling accumulation conditions. Then, a time series of RNA-seq analysis was performed and we identified 7594 lncRNAs in Pyrus pyrifolia (Burm. F.) Nakai that have not been identified. The sequence and expression of the lncRNAs were confirmed by PCR analysis. In total, 6253 lncRNAs were predicted to target protein-coding genes including 692 cis-regulated pairs (596 lncRNAs) and 13,158 trans-regulated pairs (6181 lncRNAs). Gene Ontology analysis revealed that most of lncRNAs' target genes were involved in catalytic activity, metabolic processes and cellular processes. In the trend analysis, 124 long-term cold response lncRNAs and 80 short-term cold response lncRNAs were predicted. Regarding the lncRNA-miRNA regulatory networks, 59 lncRNAs were identified as potential precursors for miRNA members of 20 families, 586 lncRNAs were targets of 261 pear miRNAs and 53 lncRNAs were endogenous target mimics for 26 miRNAs. In addition, three cold response lncRNAs, two miRNAs and their target genes were selected for expression confirmed. The trend of their expression was consistent with the predicted relationships among them and suggested possible roles of lncRNAs in ABA metabolic pathway. Our findings not only suggest the potential roles of lncRNAs in regulating the dormancy of pear floral buds but also provide new insights into the lncRNA-miRNA-mRNA regulatory network in plants.

Role of non-coding RNAs in plant immunity

Crops are exposed to attacks by various pathogens that cause substantial yield losses and severely threaten food security. To cope with pathogenic infection, crops have elaborated strategies to enhance resistance against pathogens. In addition to the role of protein-coding genes as key regulators in plant immunity, accumulating evidence has demonstrated the importance of non-coding RNAs (ncRNAs) in the plant immune response. Here, we summarize the roles and molecular mechanisms of endogenous ncRNAs, especially microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs), in plant immunity. We discuss the coordination between miRNAs and small interfering RNAs (siRNAs), between lncRNAs and miRNAs or siRNAs, and between circRNAs and miRNAs in the regulation of plant immune responses. We also address the role of cross-kingdom mobile small RNAs in plant-pathogen interactions. These insights improve our understanding of the mechanisms by which ncRNAs regulate plant immunity and can promote the development of better approaches for breeding disease-resistant crops.

Interplay between miRNAs and lncRNAs: Mode of action and biological roles in plant development and stress adaptation

Plants employ sophisticated mechanisms to control developmental processes and to cope with environmental changes at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs), two classes of endogenous noncoding RNAs, are key regulators of gene expression in plants. Recent studies have identified the interplay between miRNAs and lncRNAs as a novel regulatory layer of gene expression in plants. On one hand, miRNAs target lncRNAs for the production of phased small interfering RNAs (phasiRNAs). On the other hand, lncRNAs serve as origin of miRNAs or regulate the accumulation or activity of miRNAs at transcription and post-transcriptional levels. Theses lncRNA-miRNA interplays are crucial for plant development, physiology and responses to biotic and abiotic stresses. In this review, we summarize recent advances in the biological roles, interaction mechanisms and computational predication methods of the interplay between miRNAs and lncRNAs in plants.

Mutation of a major CG methylase alters genome-wide lncRNA expression in rice

Plant long non-coding RNAs (lncRNAs) function in diverse biological processes, and lncRNA expression is under epigenetic regulation, including by cytosine DNA methylation. However, it remains unclear whether 5-methylcytosine (5mC) plays a similar role in different sequence contexts (CG, CHG, and CHH). In this study, we characterized and compared the profiles of genome-wide lncRNA profiles (including long intergenic non-coding RNAs [lincRNAs] and long noncoding natural antisense transcripts [lncNATs]) of a null mutant of the rice DNA methyltransferase 1, OsMET1-2 (designated OsMET1-2-/-) and its isogenic wild type (OsMET1-2+/+). The En/Spm transposable element (TE) family, which was heavily methylated in OsMET1-2+/+, was transcriptionally de-repressed in OsMET1-2-/- due to genome-wide erasure of CG methylation, and this led to abundant production of specific lncRNAs. In addition, RdDM-mediated CHH hypermethylation was increased in the 5'-upstream genomic regions of lncRNAs in OsMET1-2-/-. The positive correlation between the expression of lincRNAs and that of their proximal protein-coding genes was also analyzed. Our study shows that CG methylation negatively regulates the TE-related expression of lncRNA and demonstrates that CHH methylation is also involved in the regulation of lncRNA expression.