Genome-wide analysis of long non-coding RNAs at the mature stage of sea buckthorn (Hippophae rhamnoides Linn) fruit

Advances in CircRNAs in the Past Decade: Review of CircRNAs Biogenesis, Regulatory Mechanisms, and Functions in Plants

Circular RNA (circRNA) is a type of non-coding RNA with multiple biological functions. Whole circRNA genomes in plants have been identified, and circRNAs have been demonstrated to be widely present and highly expressed in various plant tissues and organs. CircRNAs are highly stable and conserved in plants, and exhibit tissue specificity and developmental stage specificity. CircRNAs often interact with other biomolecules, such as miRNAs and proteins, thereby regulating gene expression, interfering with gene function, and affecting plant growth and development or response to environmental stress. CircRNAs are less studied in plants than in animals, and their regulatory mechanisms of biogenesis and molecular functions are not fully understood. A variety of circRNAs in plants are involved in regulating growth and development and responding to environmental stress. This review focuses on the biogenesis and regulatory mechanisms of circRNAs, as well as their biological functions during growth, development, and stress responses in plants, including a discussion of plant circRNA research prospects. Understanding the generation and regulatory mechanisms of circRNAs is a challenging but important topic in the field of circRNAs in plants, as it can provide insights into plant life activities and their response mechanisms to biotic or abiotic stresses as well as new strategies for plant molecular breeding and pest control.

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Whole-transcriptome profiling and identification of cold tolerance-related ceRNA networks in japonica rice varieties

Introduction

Low-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined.

Methods

We conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress.

Results

We identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice.

Discussion and conclusion

We identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.

The long noncoding RNAs lnc10 and lnc11 regulating flavonoid biosynthesis in Ginkgo biloba

Although long non-coding RNAs have been recognized to play important roles in plant, their possible functions and potential mechanism in Ginkgo biloba flavonoid biosynthesis are poorly understood. Flavonoids are important secondary metabolites and healthy components of Ginkgo biloba. They have been widely used in food, medicine, and natural health products. Most previous studies have focused on the molecular mechanisms of structural genes and transcription factors that regulate flavonoid biosynthesis. Few reports have examined the biological functions of flavonoid biosynthesis by long non-coding RNAs in G. biloba. Long noncoding RNAs associated with flavonoid biosynthesis in G. biloba have been identified through RNA sequencing, but the function of lncRNAs has not been reported. In this study, the expression levels of lnc10 and lnc11 were identified. Quantitative real-time polymerase chain reaction analysis revealed that lnc10 and lnc11 were expressed in all detected organs, and they showed significantly higher levels in immature and mature leaves than in other organs. In addition, to fully identify the function of lnc10 and lnc11 in flavonoid biosynthesis in G. biloba, lnc10 and lnc11 were cloned from G. biloba, and were transformed into Arabidopsis and overexpressed. Compared with the wild type, the flavonoid content was increased in transgenic plants. Moreover, the RNA-sequencing analysis of wild-type, lnc10-overexpression, and lnc11-overexpression plants screened out 2019 and 2552 differentially expressed genes, and the transcript levels of structural genes and transcription factors associated with flavonoid biosynthesis were higher in transgenic Arabidopsis than in the wild type, indicating that lnc10 and lnc11 activated flavonoid biosynthesis in the transgenic lines. Overall, these results suggest that lnc10 and lnc11 positively regulate flavonoid biosynthesis in G. biloba.

15-cis-Phytoene Desaturase and 15-cis-Phytoene Synthase Can Catalyze the Synthesis of β-Carotene and Influence the Color of Apricot Pulp

Fruit color affects its commercial value. β-carotene is the pigment that provides color for many fruits and vegetables. However, the molecular mechanism of β-carotene metabolism during apricot ripening is largely unknown. Here, we investigated whether β-carotene content affects apricot fruit color. First, the differences in β-carotene content between orange apricot 'JTY' and white apricot 'X15' during nine developmental stages (S1-S9) were compared. β-carotene contents highly significantly differed between 'JTY' and 'X15' from S5 (color transition stage) onwards. Whole-transcriptome analysis showed that the β-carotene synthesis genes 15-cis-phytoene desaturase (PaPDS) and 15-cis-phytoene synthase (PaPSY) significantly differed between the two cultivars during the color transition stage. There was a 5 bp deletion in exon 11 of PaPDS in 'X15', which led to early termination of amino acid translation. Gene overexpression and virus-induced silencing analysis showed that truncated PaPDS disrupted the β-carotene biosynthesis pathway in apricot pulp, resulting in decreased β-carotene content and a white phenotype. Furthermore, virus-induced silencing analysis showed that PaPSY was also a key gene in β-carotene biosynthesis. These findings provide new insights into the molecular regulation of apricot carotenoids and provide a theoretical reference for breeding new cultivars of apricot.

Genome-wide analysis of long non-coding RNAs in sugar beet (Beta vulgaris L.) under drought stress

Drought stress is one of the most severe abiotic stresses that restrict global crop production. Long non-coding RNAs (lncRNAs) have been proved to play a key role in response to drought stress. However, genome-wide identification and characterization of drought-responsive lncRNAs in sugar beet is still lacking. Thus, the present study focused on analyzing lncRNAs in sugar beet under drought stress. We identified 32017 reliable lncRNAs in sugar beet by strand-specific high-throughput sequencing. A total of 386 differentially expressed lncRNAs (DElncRNAs) were found under drought stress. The most significantly upregulated and downregulated lncRNAs were TCONS_00055787 (upregulated by more than 6000 fold) and TCONS_00038334 (downregulated by more than 18000 fold), respectively. Quantitative real-time PCR results exhibited a high concordance with RNA sequencing data, which conformed that the expression patterns of lncRNAs based on RNA sequencing were highly reliable. In addition, we predicted 2353 and 9041 transcripts that were estimated to be the cis- and trans-target genes of the drought-responsive lncRNAs. As revealed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the target genes of DElncRNAs were significantly enriched in organelle subcompartment, thylakoid, endopeptidase activity, catalytic activity, developmental process, lipid metabolic process, RNA polymerase activity, transferase activity, flavonoid biosynthesis and several other terms associated with abiotic stress tolerance. Moreover, 42 DElncRNAs were predicted as potential miRNA target mimics. LncRNAs have important effects on plant adaptation to drought conditions through the interaction with protein-encoding genes. The present study leads to greater insights into lncRNA biology and offers candidate regulators for improving the drought tolerance of sugar beet cultivars at the genetic level.

Coordinating Diverse Functions of miRNA and lncRNA in Fleshy Fruit

Non-coding RNAs play vital roles in the diverse biological processes of plants, and they are becoming key topics in horticulture research. In particular, miRNAs and long non-coding RNAs (lncRNAs) are receiving increased attention in fruit crops. Recent studies in horticulture research provide both genetic and molecular evidence that miRNAs and lncRNAs regulate biological function and stress responses during fruit development. Here, we summarize multiple regulatory modules of miRNAs and lncRNAs and their biological roles in fruit sets and stress responses, which would guide the development of molecular breeding techniques on horticultural crops.

Exosomal long non-coding RNA TRAFD1-4:1 derived from fibroblast-like synoviocytes suppresses chondrocyte proliferation and migration by degrading cartilage extracellular matrix in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic, autoimmune and systemic inflammatory disease affecting 1% of the population worldwide. Immune suppression of the activity and progress of RA is vital to reduce the disability and mortality rate as well as improve the quality of life of RA patients. However, the immune molecular mechanism of RA has not been clarified yet. Our results indicated that exosomes derived from TNFα-stimulated RA fibroblast-like synoviocytes (RA-FLSs) suppressed chondrocyte proliferation and migration through modulating cartilage extracellular matrix (CECM) determining by MTS assay, cell cycle analysis, Transwell assay and Western blot (WB). Besides, RNA sequencing and verification by qRT-PCR revealed that exosomal long non-coding RNA (lncRNA) tumor necrosis factor-associated factor 1 (TRAF1)-4:1 derived from RA-FLSs treated with TNFα was a candidate lncRNA, which also inhibited chondrocyte proliferation and migration through degrading CECM. Moreover, RNA sequencing and bioinformatics analysis identified that C-X-C motif chemokine ligand 1 (CXCL1) was a target mRNA of miR-27a-3p while miR-27a-3p was a target miRNA of lnc-TRAF1-4:1 in chondrocytes. Mechanistically, lnc-TRAF1-4:1 upregulated CXCL1 expression through sponging miR-27a-3p as a competing endogenous RNA (ceRNA) in chondrocytes identifying by Dual-luciferase reporter gene assay. Summarily, exosomal lncRNA TRAFD1-4:1 derived from RA-FLSs suppressed chondrocyte proliferation and migration through degrading CECM by upregulating CXCL1 as a sponge of miR-27a-3p. This study uncovered a novel RA-related lncRNA and investigated the roles of RA-FLS-derived exosomes and exosomal lnc-TRAF1-4:1 in articular cartilage impairment, which might provide novel therapeutic targets for RA.

Genome-wide identification and comprehensive analysis reveal potential roles of long non-coding RNAs in fruit development of southern highbush blueberry (Vaccinium corymbosum L.)

Introduction

Blueberries have a high antioxidant content and are produced as healthy food worldwide. Long non-coding RNAs (lncRNAs) are a type of regulatory RNAs that play a variety of roles in plants. Nonetheless, information on lncRNAs and their functions during blueberry fruit development is scarce in public databases.

Methods

In the present study, we performed genome-wide identification of lncRNAs in a southern highbush blueberry using strand-specific RNA sequencing (ssRNA-Seq). Differentially expressed lncRNAs (DE-lncRNAs) and their potential target genes were analyzed at four stages of fruit development. Cis-regulatory DE-lncRNAs were predicted using co-localization analysis.

Results

These findings included a total of 25,036 lncRNAs from 17,801 loci. Blueberry lncRNAs had shorter transcript lengths, smaller open reading frame (ORF) sizes, fewer exons, and fewer isoforms than protein-coding RNAs, as well as lower expression levels and higher stage-specificity during fruit development. A total of 105 DE-lncRNAs were identified among the comparison group of PAD vs. CUP, 443 DE-lncRNAs were detected when comparing CUP with PINK fruits, and 285 DE-lncRNAs were revealed when comparing PINK and BLUE fruits. According to Kyoto Encyclopedia of Genes and Genomes annotation, target genes of DE-lncRNAs were primarily enriched in the "Autophagy-other", "DNA replication", "Endocytosis", 'photosynthesis' and 'chlorophyll metabolism' pathways, suggesting that lncRNAs may pay potential roles in fruit expansion and ripening. Moreover, several lncRNAs have been proposed as cis-regulators of the key genes involved in flavonoid biosynthesis. MSTRG.107242.6, and its putative target gene, BTB/POZ and TAZ domain-containing protein, might play critical roles in anthocyanin accumulation in blueberries.

Discussion

These findings highlight the regulatory function of lncRNAs and aid in elucidating the molecular mechanism underlying blueberry fruit growth.

Identification and Analysis of Long Non-Coding RNAs Related to UV-B-Induced Anthocyanin Biosynthesis During Blood-Fleshed Peach (Prunus persica) Ripening

Blood flesh is a key fruit trait in peaches (Prunus persica) and can be attributed to the accumulation of anthocyanins. The roles of long non-coding RNAs (lncRNAs) have been highlighted by multiple studies in regulating fruit ripening, anthocyanin accumulation, and abiotic stress responses in many flowering plants. Such regulatory functions of lncRNAs in Prunus persica, nonetheless, have not been reported. In this research, we sequenced and analyzed the complete transcriptome of C3-20 (a blood-fleshed peach) fruit at four developmental stages. Analyses of the correlated genes and differentially expressed lncRNA target genes helped to forecast lncRNAs’ possible functions. The RNA-seq data were generated using high-throughput sequencing. In total, 17,456 putative lncRNAs, including 4,800 intergenic lncRNAs, 2,199 antisense lncRNAs, and 10,439 intronic lncRNAs were discovered, of which 4,871 differentially expressed lncRNAs (DE-lncRNAs) were annotated in the fruit developmental processes. The target genes of these DE-lncRNAs and their regulatory relationship identifying 21,795 cis-regulated and 18,271 trans-regulated targets of the DE-lncRNAs were in a similar way predicted by us. The enriched GO terms for the target genes included anthocyanin biosynthesis. Flavonoid biosynthesis and plant hormone signal transduction were also included in the enriched KEGG pathways. Co-expression network construction demonstrated that the highly expressed genes might co-regulate multiple other genes associated with auxin signal transduction and take effect in equal pathways. We discovered that lncRNAs, including LNC_000987, LNC_000693, LNC_001323, LNC_003610, LNC_001263, and LNC_003380, correlated with fruit that ripened and could take part in ethylene biosynthesis and metabolism and the ABA signaling pathway. Several essential transcription factors, such as ERFs, WRKY70, NAC56, and NAC72, may in a similar way regulate fruit ripening. Three DE-lncRNAs, XLOC_011933, XLOC_001865, and XLOC_042291, are involved in UV-B-induced anthocyanin biosynthesis and positively regulating UVR8 and COP10, were identified and characterized. Our discovery and characterization of XLOC_011933, XLOC_001865, and XLOC_042291 provide a more precise understanding and preliminarily establishes a theoretical framework for UV-B-induced flesh anthocyanin biosynthesis. This phenomenon might encourage more in-depth investigations to study the molecular mechanisms underlying peach flesh coloring.

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.

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.

Transcriptome and DNA methylome provide insights into the molecular regulation of drought stress in sea buckthorn

Sea buckthorn is a typical drought-resistant tree species. However, there is a general lack of understanding of the pattern of DNA methylation linked with sea buckthorn responses to drought, and its relationship with drought tolerance mechanisms. In this study, we performed whole-transcriptome RNA sequencing and methylome sequencing in response to drought stress to explore differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs in sea buckthorn leaves. Based on predicted DE pairs, we constructed a competitive endogenous RNA network, which revealed potential transcriptional regulatory roles in response to drought stress. The results of methylome sequencing revealed that the DNA methylation level was increased in sea buckthorn leaves under drought stress. We identified 13,405 differentially methylated regions between CK and TR. We found one DMR-associated DEG (Vacuolar-sorting receptor 6) involved in the ABA accumulation pathway. In addition, two DNA methyltransferases (HrMET1 and HrDRM1) were closely associated with drought-induced hypermethylation in sea buckthorn. Together, we firstly conducted a comprehensive transcriptomic and epigenetic analysis of sea buckthorn under drought stress, providing a resource for further study of the potential functions of genes, miRNAs, lncRNAs, circRNAs and DNA methyltransferases.

HrTCP20 dramatically enhance drought tolerance of sea buckthorn (Hippophae rhamnoides L). by mediating the JA signaling pathway

Sea buckthorn, an important ecological and economical tree species, have remarkable drought and salt resistance. The plant-specific transcription factor TCPs play important roles in plant growth, development, and stress responses. However, in sea buckthorn, the molecular mechanism of TCP proteins and their involvement in drought stress are unknown. Here, we found that the expression of HrTCP20 was significantly up-regulated in sea buckthorn under drought stress. Overexpression of HrTCP20 in Arabidopsis thaliana showed that the superoxide dismutase (SOD), polyphenol oxidase (POD), and chlorophyll (SPAD) content was significantly increased by 1.37 and 1.35 times. However, the malondialdehyde (MDA) content decreased by 0.51 times. Our studies further confirmed that silencing HrTCP20 by virus-induced gene silencing (VIGS) led to a decrease in the content of defense enzymes, relative water content (RWC), and an increase of relative electrical conductivity (REC). Silencing HrTCP20 also caused the jasmonic acid (JA) content to decrease in the VIGS-treated tree. Interestingly, we found that JA accumulation content and the expression of HrLOX2, an essential enzyme for JA synthesis, was significantly inhibited in HrTCP20-silenced sea buckthorn under drought stress. Yeast two-hybrid analysis also showed that HrTCP20 is directly bound to HrLOX2. Taken together, the HrTCP20 transcription factor was a positive regulator in drought resistance of sea buckthorn. Further, our findings will provide comprehensive insights into the forest tree defence system of drought stress.

Analysis of long non-coding RNAs and mRNAs in harvested kiwifruit in response to the yeast antagonist, Wickerhamomyces anomalus

•

W. anomalus exhibits good biocontrol activity against blue and gray mold on kiwifruit.

•

LncRNAs in kiwifruit may be involved in activating plant hormone signal transduction pathways in response to the biocontrol yeast.

•

LncRNAs in kiwifruit may modulate the production of related TFs and secondary metabolites.

•

The expression of downstream defense-related genes in kiwifruit increases in response to the application of the biocontrol yeast.

Biological control utilizing antagonistic yeasts is an effective method for controlling postharvest diseases. Long non-coding RNAs (lncRNAs) have been found to be involved in a variety of plant growth and development processes, including those associated with plant disease resistance. In the present study, the yeast antagonist, Wickerhamomyces anomalus, was found to strongly inhibit postharvest blue mold (Penicillium expansum) and gray mold (Botrytis cinerea) decay of kiwifruit. Additionally, lncRNA high-throughput sequencing and bioinformatic analysis was used to identify lncRNAs in W. anomalus-treated wounds in kiwifruit and predict their function based on putative target genes. Our results indicate that lncRNAs may be involved in increasing ethylene (ET), jasmonic acid (JA), abscisic acid (ABA), and auxin (IAA) levels, as well as activating signal transduction pathways that regulate the expression of several transcription factors (WRKY72, WRKY53, JUB1AP2). These transcription factors (TFs) then mediate the expression of downstream, defense-related genes (ZAR1, PAD4, CCR4, NPR4) and the synthesis of secondary metabolites, thus, potentially enhancing disease resistance. Notably, by stimulating the accumulation of antifungal compounds, such as phenols and lignin, disease resistance in kiwifruit was enhanced. Our study provides new information on the mechanism underlying the induction of disease resistance in kiwifruit by W. anomalus, as well as a new disease resistance strategy that can be used to enhance the defense response of fruit to pathogenic fungi.

Identification of Long Non-Coding RNAs Associated with Tomato Fruit Expansion and Ripening by Strand-Specific Paired-End RNA Sequencing

As emerging essential regulators in plant development, long non-coding RNAs (lncRNAs) have been extensively investigated in multiple horticultural crops, as well as in different tissues of plants. Tomato fruits are an indispensable part of people’s diet and are consumed as fruits and vegetables. Meanwhile, tomato is widely used as a model to study the ripening mechanism in fleshy fruit. Although increasing evidence shows that lncRNAs are involved in lots of biological processes in tomato plants, the comprehensive identification of lncRNAs in tomato fruit during its expansion and ripening and their functions are partially known. Here, we performed strand-specific paired-end RNA sequencing (ssRNA-seq) of tomato Heinz1706 fruits at five different developmental stages, as well as flowers and leaves. We identified 17,674 putative lncRNAs by referencing the recently released SL4.0 and annotation ITAG4.0 in tomato plants. Many lncRNAs show different expression patterns in fleshy fruit at different developmental stages compared with leaves or flowers. Our results indicate that lncRNAs play an important role in the regulation of tomato fruit expansion and ripening, providing informative lncRNA candidates for further studies in tomato fruits. In addition, we also summarize the recent advanced progress in lncRNAs mediated regulation on horticultural fruits. Hence, our study updates the understanding of lncRNAs in horticultural plants and provides resources for future studies relating to the expansion and ripening of tomato fruits.

Identification and characterization of early Fusarium wilt responsive mRNAs and long non-coding RNAs in banana root using high-throughput sequencing

Fusarium wilt disease, caused by Fusarium oxysporum f.sp. cubense (Foc), has been recognized as the most devastating disease to banana. The regulatory role of long non-coding RNAs (lncRNAs) in plant defense has been verified in many plant species. However, the understanding of their role during early FocTR4 (Foc tropical race 4) infection stage is very limited. In this study, lncRNA sequencing was used to reveal banana root transcriptome profile changes during early FocTR4 infection stages. Quantitative real time PCR (qRT-PCR) was performed to confirm the expression of eight differentially expressed (DE) lncRNAs (DELs) and their predicted target genes (DETs), and three DE genes (DEGs). Totally, 12,109 lncRNAs, 36,519 mRNAs and 2642 novel genes were obtained, of which 1398 (including 78 DELs, 1220 DE known genes and 100 DE novel genes) were identified as FocTR4 responsive DE transcripts. Gene function analysis revealed that most DEGs were involved in biosynthesis of secondary metabolites, plant–pathogen interaction, plant hormone signal transduction, phenylalanine metabolism, phenylpropanoid biosynthesis, alpha-linolenic acid metabolism and so on. Coincidently, many DETs have been identified as DEGs in previous transcriptome studies. Moreover, many DETs were found to be involved in ribosome, oxidative phosphorylation, lipoic acid metabolism, ubiquitin mediated proteolysis, N-glycan biosynthesis, protein processing in endoplasmic reticulum and DNA damage response pathways. QRT-PCR result showed the expression patterns of the selected transcripts were mostly consistent with our lncRNA sequencing data. Our present study showed the regulatory role of lncRNAs on known biotic and abiotic stress responsive genes and some new-found FocTR4 responsive genes, which can provide new insights into FocTR4-induced changes in the banana root transcriptome during the early pathogen infection stage.

NGS Methodologies and Computational Algorithms for the Prediction and Analysis of Plant Circular RNAs

Circular RNAs (circRNAs) are a class of single-stranded RNAs derived from exonic, intronic, and intergenic regions from precursor messenger RNAs (pre-mRNA), where a noncanonical back-splicing event occurs, in which the 5' and 3' ends are attached by covalent bond. CircRNAs participate in the regulation of gene expression at the transcriptional and posttranscriptional level primarily as miRNA and RNA-binding protein (RBP) sponges, but also involved in the regulation of alternative RNA splicing and transcription. CircRNAs are widespread and abundant in plants where they have been involved in stress responses and development. Through the analysis of all publications in this field in the last five years, we can summarize that the identification of these molecules is carried out through next generation sequencing studies, where samples have been previously treated to eliminate DNA, rRNA, and linear RNAs as a means to enrich circRNAs. Once libraries are prepared, they are sequenced and subsequently studied from a bioinformatics point of view. Among the different tools for identifying circRNAs, we can highlight CIRI as the most used (in 60% of the published studies), as well as CIRCExplorer (20%) and find_circ (20%). Although it is recommended to use more than one program in combination, and preferably developed specifically to treat with plant samples, this is not always the case. It should also be noted that after identifying these circular RNAs, most of the authors validate their findings in the laboratory in order to obtain bona fide results.

Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening

Kiwifruit has gained increasing attention worldwide for its unique flavor and high nutritional value. Rapid softening after harvest greatly shortens its shelf-life and reduces the commercial value. Therefore, it is imperative and urgent to identify and clarify its softening mechanism. This study aimed to analyze and compare the long noncoding RNA (lncRNA) and mRNA expression patterns in ABA-treated (ABA) and room temperature (RT)-stored fruits with those in freshly harvested fruits (CK) as control. A total of 697 differentially expressed genes (DEGs) and 81 differentially expressed lncRNAs (DELs) were identified while comparing ABA with CK, and 458 DEGs and 143 DELs were detected while comparing RT with CK. The Kyoto Encyclopedia of Genes and Genomes analysis of the identified DEGs and the target genes of DELs revealed that genes involved in starch and sucrose metabolism, brassinosteroid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis accounted for a large part. The co-localization networks, including 38 DEGs and 31 DELs in ABA vs. CK, and 25 DEGs and 25 DELs in RT vs. CK, were also performed. Genes related to fruit ripening, such as genes encoding β-galactosidase, mannan endo-1,4-β-mannosidase, pectinesterase/pectinesterase inhibitor, and NAC transcription factor, were present in the co-localization network, suggesting that lncRNAs were involved in regulating kiwifruit ripening. Notably, several ethylene biosynthesis- and signaling-related genes, including one 1-aminocyclopropane-1-carboxylic acid oxidase gene and three ethylene response factor genes, were found in the co-localization network of ABA vs. CK, suggesting that the promoting effect of ABA on ethylene biosynthesis and fruit softening might be embodied by increasing the expression of these lncRNAs. These results may help understand the regulatory mechanism of lncRNAs in ripening and ABA-induced fruit softening of kiwifruit.

Characterization and complete genome sequences of two novel variants of the family Closteroviridae from Chinese kiwifruit

Two distinct closterovirus-like genome sequences (termed AdV-1 v1 and v2) were identified in Actinidia chinensis var. deliciosa 'Miliang-1' that had no disease symptoms using high-throughput sequencing. Using overlapping reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends, the genomic sequences of AdV-1 v1 and v2 were confirmed as 17,646 and 18,578 nucleotides in length, respectively. The two complete genomes contained 9 and 15 open reading frames, respectively, coding for proteins having domains typical of Closteroviridae, such as RNA-dependent RNA polymerase (RdRp), heat shock protein 70 homolog (HSP70h) and coat protein (CP). Sequence analysis showed that the amino acid sequences of RdRp, HSP70h, and CP of the two variants exhibited high similarity (> 80%), while their genomic organization was somewhat different. This suggested that the two viral genomes identified here are variants of the family Closteroviridae in a single kiwifruit host. Furthermore, phylogenetic relationship analysis revealed that the two variants had a closer relationship with the unclassified virus Persimmon virus B (PeVB) and Actinidia virus 1 (AcV-1) than with other members of the family Closteroviridae, as did their genomic organization. It is speculated that the two variants, together with PeVB and AcV-1 belong to a new subfamily of Closteroviridae.

Discovery, identification, and functional characterization of long noncoding RNAs in Arachis hypogaea L

BACKGROUND

Long noncoding RNAs (lncRNAs), which are typically > 200 nt in length, are involved in numerous biological processes. Studies on lncRNAs in the cultivated peanut (Arachis hypogaea L.) largely remain unknown.

RESULTS

A genome-wide scan of the peanut (Arachis hypogaea L.) transcriptome identified 1442 lncRNAs, which were encoded by loci distributed over every chromosome. Long intergenic noncoding RNAs accounted for 85.58% of these lncRNAs. Additionally, 189 lncRNAs were differentially abundant in the root, leaf, or seed. Generally, lncRNAs showed lower expression levels, tighter tissue-specific expression, and less splicing than mRNAs. Approximately 44.17% of the lncRNAs with an exon/intron structure were alternatively spliced; this rate was slightly lower than the splicing rate of mRNA. Transcription at the start site event was the alternative splicing (AS) event with the highest frequency (28.05%) in peanut lncRNAs, whereas the occurrence rate (30.19%) of intron retention event was the highest in mRNAs. AS changed the target gene profiles of lncRNAs and increased the diversity and flexibility of lncRNAs, which may be important for lncRNAs to execute their functions. Additionally, a substantial number of the peanut AS isoforms generated from protein-encoding genes appeared to be noncoding because they were truncated transcripts; such isoforms can be legitimately regarded as a class of lncRNAs. The predicted target genes of the lncRNAs were involved in a wide range of biological processes. Furthermore, expression pattern of several selected lncRNAs and their target genes were examined under salt stress, results showed that all of them could respond to salt stress in different manners.

CONCLUSIONS

This study provided a resource of candidate lncRNAs and expression patterns across tissues, and whether these lncRNAs are functional will be further investigated in our subsequent experiments.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.

RESULTS

A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.

CONCLUSION

The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Identification of Novel lncRNAs Differentially Expressed in Placentas of Chinese Ningqiang Pony and Yili Horse Breeds

As a nutrient sensor, the placenta plays a key role in regulating fetus growth and development. Long non-coding RNAs (lncRNAs) have been shown to regulate growth-related traits. However, the biological function of lncRNAs in horse placentas remains unclear. To compare the expression patterns of lncRNAs in the placentas of the Chinese Ningqiang (NQ) and Yili (YL) breeds, we performed a transcriptome analysis using RNA sequencing (RNA-seq) technology. NQ is a pony breed with an average adult height at the withers of less than 106 cm, whereas that of YL is around 148 cm. Based on 813 million high-quality reads and stringent quality control procedures, 3011 transcripts coding for 1464 placental lncRNAs were identified and mapped to the horse reference genome. We found 107 differentially expressed lncRNAs (DELs) between NQ and YL, including 68 up-regulated and 39 down-regulated DELs in YL. Six (TBX3, CACNA1F, EDN3, KAT5, ZNF281, TMED2, and TGFB1) out of the 233 genes targeted by DELs were identified as being involved in limb development, skeletal myoblast differentiation, and embryo development. Two DELs were predicted to target the TBX3 gene, which was found to be under strong selection and associated with small body size in the Chinese Debao pony breed. This finding suggests the potential functional significance of placental lncRNAs in regulating horse body size.

Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots

BACKGROUND

Environmental stimuli can activate a series of physiological and biochemical responses in plants accompanied by extensive transcriptional reprogramming. Long non-coding RNAs (lncRNAs), as versatile regulators, control gene expression in multiple ways and participate in the adaptation to biotic and abiotic stresses.

RESULTS

In this study, soybean seedlings were continuously cultured for 15 days with high salinity solutions started from seed germination. Strand-specific whole transcriptome sequencing and stringent bioinformatic analysis led to the identification of 3030 long intergenic non-coding RNAs (lincRNAs) and 275 natural antisense transcripts (lncNATs) in soybean roots. In contrast to mRNAs, newly identified lncRNAs exhibited less exons, similar AU content to UTRs, even distribution across the genome and low evolutionary conservation. Remarkably, more than 75% of discovered lncRNAs that were activated or up-regulated by continuous salt stress mainly targeted proteins with binding and catalytic activities. Furthermore, two DNA methylation maps with single-base resolution were generated by using reduced representation bisulfite sequencing, offering a genome-wide perspective and important clues for epigenetic regulation of stress-associated lncRNAs and protein-coding genes.

CONCLUSIONS

Taken together, our findings systematically demonstrated the characteristics of continuous salt stress-induced lncRNAs and extended the knowledge of corresponding methylation profiling, providing valuable evidence for a better understanding of how plants cope with long-term salt stress circumstances.

Plant Noncoding RNAs: Hidden Players in Development and Stress Responses

A large and significant portion of eukaryotic transcriptomes consists of noncoding RNAs (ncRNAs) that have minimal or no protein-coding capacity but are functional. Diverse ncRNAs, including both small RNAs and long ncRNAs (lncRNAs), play essential regulatory roles in almost all biological processes by modulating gene expression at the transcriptional and posttranscriptional levels. In this review, we summarize the current knowledge of plant small RNAs and lncRNAs, with a focus on their biogenesis, modes of action, local and systemic movement, and functions at the nexus of plant development and environmental responses. The complex connections among small RNAs, lncRNAs, and small peptides in plants are also discussed, along with the challenges of identifying and investigating new classes of ncRNAs.

Long noncoding RNAs and miRNAs regulating terpene and tartaric acid biosynthesis in rose-scented geranium

This study aimed to explore the noncoding RNAs, which have emerged as key regulatory molecules in biological processes, in rose-scented geranium. We analyzed RNA-seq data revealing 26 784 long noncoding RNAs (lncRNAs) and 871 miRNAs in rose-scented geranium. A total of 466 lncRNAs were annotated using different plant lncRNA public databases. Furthermore, 372 lncRNAs and 99 miRNAs were detected that target terpene and tartarate biosynthetic pathways. An interactome, comprising of lncRNAs, miRNAs, and mRNAs, was constructed that represents a noncoding RNA regulatory network of the target mRNAs. Real-time quantitative PCR expression validation was done for selected lncRNAs involved in the regulation of terpene and tartaric acid pathways. This study provides the first insights into the regulatory functioning of noncoding RNAs in rose-scented geranium.

Genome-wide analysis of long non-coding RNA expression profile in porcine circovirus 2-infected intestinal porcine epithelial cell line by RNA sequencing

Porcine circovirus-associated disease (PCVAD), which is induced by porcine circovirus type 2 (PCV2), is responsible for severe economic losses. Recently, the role of noncoding RNAs, and in particular microRNAs, in PCV2 infection has received great attention. However, the role of long noncoding RNA (lncRNA) in PCV2 infection is unclear. Here, for the first time, we describe the expression profiles of lncRNAs in an intestinal porcine epithelial cell line (IPEC-J2) after PCV2 infection, and analyze the features of differently expressed lncRNAs and their potential target genes. After strict filtering of approximately 150 million reads, we identified 13,520 lncRNAs, including 199 lncRNAs that were differentially expressed in non-infected and PCV2-infected cells. Furthermore, trans analysis found lncRNA-regulated target genes enriched for specific Gene Ontology terms (P < 0.05), such as DNA binding, RNA binding, and transcription factor activity, which are closely associated with PCV2 infection. In addition, we analyzed the predicted target genes of differentially expressed lncRNAs, including SOD2, TNFAIP3, and ARG1, all of which are involved in infectious diseases. Our study identifies many candidate lncRNAs involved in PCV2 infection and provides new insight into the mechanisms underlying the pathogenesis of PCVAD.

Identification and characterization of circular RNAs during the sea buckthorn fruit development

As a rising star of noncoding RNA, circular RNAs (circRNAs) have a covalently closed loop structure, which formed by 3'-5' ligation during splicing. A few circRNAs were identified and thought to be transcriptional noise due to cognitive defect over the past 40 years. Recently, with the development of high-throughput RNA sequencing techniques and specific algorithms for circRNA detection and quantification, plenty of potential circRNAs were identified in many species which play important roles in various biological processes. However, researches on circRNAs in fruit ripening process were lacking. Here, we totally identified 2616 circRNAs in sea buckthorn fruit development process, which uniformly distributed in sea buckthorn chromosome. Among them, 1721 (65.8%) circRNAs were arising from the exons of their host genes, 252 circRNAs were identified as the differentially expressed circRNAs (DEcircRNAs) between three different development stages, and 181 (71.8%) DEcircRNAs had sequence similarity with 235 identified circRNAs from five know plant species. Functional annotation revealed that host genes of DEcircRNAs were predicted to be involved in carotenoid biosynthesis, lipid synthesis and plant hormone signal transduction. Additionally, 53 DEcircRNAs were predicted as the corresponding nine miRNAs sponges in sea buckthorn. Divergent reverse-transcription PCR and RT-qPCR were used for validate the differential expression and back-splicing sites of six DEcircRNAs. These results revealed the role of circRNAs in sea buckthorn fruit ripening process and promoted the noncoding RNA researches in plants.

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

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

Transcriptomic and functional analyses unveil the role of long non-coding RNAs in anthocyanin biosynthesis during sea buckthorn fruit ripening

Fruit ripening is a developmental process regulated by a complex network of endogenous and exogenous cues. Sea buckthorn is an excellent material for fruit ripening studies due to its dramatic ripening process and high contents of nutritional and anti-oxidant compounds in berries. Here, the whole transcriptome of sea buckthorn fruit at three development stages were analysed using multiple high-throughput sequencings. We assembled and annotated 9,008 long non-coding RNAs (lncRNAs) in sea buckthorn fruits, and identified 118 differentially expressed lncRNAs (DE-lncRNAs) and 32 differentially expressed microRNAs in fruit developmental process. In addition, we predicted 1,061 cis-regulated and 782 trans-regulated targets of DE-lncRNAs, and these DE-lncRNAs are specifically enriched in the biosynthesis of ascorbic acid, carotenoids and flavonoids. Moreover, the silencing of two lncRNAs (LNC1 and LNC2) in vivo and expression analysis revealed that LNC1 and LNC2 can act as endogenous target mimics of miR156a and miR828a to reduce SPL9 and induce MYB114 expression, respectively, which lead to increased and decreased anthocyanin content as revealed by high-performance liquid chromatography analysis. Our results present the first global functional analysis of lncRNA in sea buckthorn and provide two essential regulators of anthocyanin biosynthesis, which provides new insights into the regulation of fruit quality.