Deciphering the roles of circRNAs on chilling injury in tomato

AScirRNA: A novel computational approach to discover abiotic stress-responsive circular RNAs in plant genome

In the realm of plant biology, understanding the intricate regulatory mechanisms governing stress responses stands as a pivotal pursuit. Circular RNAs (circRNAs), emerging as critical players in gene regulation, have garnered attention in recent days for their potential roles in abiotic stress adaptation. A comprehensive grasp of circRNAs' functions in stress response offers avenues for breeders to manipulating plants to develop abiotic stress resistant crop cultivars to thrive in challenging climates. This study pioneers a machine learning-based model for predicting abiotic stress-responsive circRNAs. The K-tuple nucleotide composition (KNC) and Pseudo KNC (PKNC) features were utilized to numerically represent circRNAs. Three different feature selection strategies were employed to select relevant and non-redundant features. Eight shallow and four deep learning algorithms were evaluated to build the final predictive model. Following five-fold cross-validation process, XGBoost learning algorithm demonstrated superior performance with LightGBM-chosen 260 KNC features (Accuracy: 74.55 %, auROC: 81.23 %, auPRC: 76.52 %) and 160 PKNC features (Accuracy: 74.32 %, auROC: 81.04 %, auPRC: 76.43 %), over other combinations of learning algorithms and feature selection techniques. Further, the robustness of the developed models were evaluated using an independent test dataset, where the overall accuracy, auROC and auPRC were found to be 73.13 %, 72.34 % and 72.68 % for KNC feature set and 73.52 %, 79.53 % and 73.09 % for PKNC feature set, respectively. This computational approach was also integrated into an online prediction tool, AScirRNA (https://iasri-sg.icar.gov.in/ascirna/) for easy prediction by the users. Both the proposed model and the developed tool are poised to augment ongoing efforts in identifying stress-responsive circRNAs in plants.

Drought-induced circular RNAs in maize roots: Separating signal from noise

Circular RNAs (circRNAs) play an important role in diverse biological processes; however, their origin and functions, especially in plants, remain largely unclear. Here, we used 2 maize (Zea mays) inbred lines, as well as 14 of their derivative recombination inbred lines with different drought sensitivity, to systematically characterize 8,790 circRNAs in maize roots under well-watered (WW) and water-stress (WS) conditions. We found that a diverse set of circRNAs expressed at significantly higher levels under WS. Enhanced expression of circRNAs was associated with longer flanking introns and an enrichment of long interspersed nuclear element retrotransposable elements. The epigenetic marks found at the back-splicing junctions of circRNA-producing genes were markedly different from canonical splicing, characterized by increased levels of H3K36me3/H3K4me1, as well as decreased levels of H3K9Ac/H3K27Ac. We found that genes expressing circRNAs are subject to relaxed selection. The significant enrichment of trait-associated sites along their genic regions suggested that genes giving rise to circRNAs were associated with plant survival rate under drought stress, implying that circRNAs play roles in plant drought responses. Furthermore, we found that overexpression of circMED16, one of the drought-responsive circRNAs, enhances drought tolerance in Arabidopsis (Arabidopsis thaliana). Our results provide a framework for understanding the intricate interplay of epigenetic modifications and how they contribute to the fine-tuning of circRNA expression under drought stress.

PlantCircRNA: a comprehensive database for plant circular RNAs

Circular RNAs (circRNAs) represent recently discovered novel regulatory non-coding RNAs. While they are present in many eukaryotes, there has been limited research on plant circRNAs. We developed PlantCircRNA (https://plant.deepbiology.cn/PlantCircRNA/) to fill this gap. The two most important features of PlantCircRNA are (i) it incorporates circRNAs from 94 plant species based on 39 245 RNA-sequencing samples and (ii) it imports the original AtCircDB and CropCircDB databases. We manually curated all circRNAs from published articles, and imported them into the database. Furthermore, we added detailed information of tissue as well as abiotic stresses to the database. To help users understand these circRNAs, the database includes a detection score to measure their consistency and a naming system following the guidelines recently proposed for eukaryotes. Finally, we developed a comprehensive platform for users to visualize, analyze, and download data regarding specific circRNAs. This resource will serve as a home for plant circRNAs and provide the community with unprecedented insights into these mysterious molecule.

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.

Comparative genome-wide analysis of circular RNAs in Brassica napus L.: target-site versus non-target-site resistance to herbicide stress

Circular RNAs (circRNAs), a class of non-coding RNA molecules, are recognized for their unique functions; however, their responses to herbicide stress in Brassica napus remain unclear. In this study, the role of circRNAs in response to herbicide treatment was investigated in two rapeseed cultivars: MH33, which confers non-target-site resistance (NTSR), and EM28, which exhibits target-site resistance (TSR). The genome-wide circRNA profiles of herbicide-stressed and non-stressed seedlings were analyzed. The findings indicate that NTSR seedlings exhibited a greater abundance of circRNAs, shorter lengths of circRNAs and their parent genes, and more diverse functions of parent genes compared with TSR seedlings. Compared to normal-growth plants, the herbicide-stressed group exhibited similar trends in the number of circRNAs, functions of parent genes, and differentially expressed circRNAs as observed in NTSR seedlings. In addition, a greater number of circRNAs that function as competing microRNA (miRNA) sponges were identified in the herbicide stress and NTSR groups compared to the normal-growth and TSR groups, respectively. The differentially expressed circRNAs were validated by qPCR. The differntially expressed circRNA-miRNA networks were predicted, and the mRNAs targeted by these miRNAs were annotated. Our results suggest that circRNAs play a crucial role in responding to herbicide stress, exhibiting distinct responses between NTSR and TSR in rapeseed. These findings offer valuable insights into the mechanisms underlying herbicide resistance in rapeseed.

CircRNA: An emerging star in plant research: A review

CircRNAs are a class of covalently closed non-coding RNA formed by linking the 5' terminus and the 3' terminus after reverse splicing. CircRNAs are widely found in eukaryotes, and they are highly conserved, with spatio-temporal expression specificity and stability. CircRNAs can act as miRNA sponges to regulate the expression of downstream target genes, regulating the transcription of parental genes and some can even be translated into peptides or proteins. Research on circRNAs in plants is still in its infancy compared to that in animals. With the deepening of research, the results of a variety of plant circRNAs suggest that they play an important role in growth and development, and tolerance towards abiotic stresses such as salt, drought, low temperature, high temperature and other adverse environments. In this review paper, we elaborated the molecular characteristics, mechanism of action, function and bioinformatics databases of plant circRNAs, combined with the progress of circRNA research in animals, discussed the potential mechanism of action of plant circRNAs, and proposed the unsolved problems and prospects for future application of plant circRNAs.

Identification and Characterisation of the CircRNAs Involved in the Regulation of Leaf Colour in Quercus mongolica

Circular RNAs (circRNAs) are important regulatory molecules involved in various biological processes. However, the potential function of circRNAs in the turning red process of Quercus mongolica leaves is unclear. This study used RNA-seq data to identify 6228 circRNAs in leaf samples from four different developmental stages and showed that 88 circRNAs were differentially expressed. A correlation analysis was performed between anthocyanins and the circRNAs. A total of 16 circRNAs that may be involved in regulating the colour of Mongolian oak leaves were identified. CircRNAs may affect the colour of Q. mongolica leaves by regulating auxin, cytokinin, gibberellin, ethylene, and abscisic acid. This study revealed the potential role of circRNAs in the colour change of Q. mongolica leaves.

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

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

Non-coding RNAs and leaf senescence: Small molecules with important roles

Non-coding RNAs (ncRNAs) are a special class of functional RNA molecules that are not translated into proteins. ncRNAs have emerged as pivotal regulators of diverse developmental processes in plants. Recent investigations have revealed the association of ncRNAs with the regulation of leaf senescence, a complex and tightly regulated developmental process. However, a comprehensive review of the involvement of ncRNAs in the regulation of leaf senescence is still lacking. This manuscript aims to summarize the molecular mechanisms underlying ncRNAs-mediated leaf senescence and the potential applications of ncRNAs to manipulate the onset and progression of leaf senescence. Various classes of ncRNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), are discussed in terms of their regulatory mechanisms in leaf senescence. Furthermore, we explore the interactions between ncRNA and the key regulators of senescence, including transcription factors as well as core components in phytohormone signaling pathways. We also discuss the possible challenges and approaches related to ncRNA-mediated leaf senescence. This review contributes to a further understanding of the intricate regulatory network involving ncRNAs in leaf senescence.

An accretive detection method for in silico identification and validation of circular RNAs in wheat (Triticum aestivum L.) using RT-qPCR

BACKGROUND

Circular RNAs (circRNAs) are novel class of non-coding RNAs, which are involved in various functions at the transcriptional and post-transcriptional level in response to a fungal pathogen (Puccinia triticina), including microRNA (miRNA) sponge, RNA binding proteins sponge, regulation of parental gene and biomarkers. Detailed analysis of wheat circRNAs is essential to accelerate the regulated expression of fungal miRNAs. Therefore, we suggest a protocol to aid circRNA identification through RNA-Seq data using various algorithms based on perl script followed by validation through divergent primer designing, standard PCR, and RT-qPCR assays.

METHODS AND RESULT

The divergent primer has been widely used to detect, validate, and quantify back-spliced junction (BSJ) of circRNAs. The procedure covers index file formation, circRNA identification and BSJ detections. However, the laboratory validation of circRNA includes wheat genomic DNA isolation, RNA isolation and its cDNA conversion upto validation. In this study, we identified 28 circRNAs from RNA-Seq of S0 and R0, wherein six circRNAs are commonly present and 75% of the identified circRNAs were belongs to inter-genic, 14% were exonic and intronic category were 11%. Divergent primer designing method successfully validated the two circRNAs via RT-qPCR assay, where circRNA_2 showed less relative expression pattern than circRNA_1 in contrast with housekeeping genes.

CONCLUSION

Thus, our results of identified and validated circRNAs showed that, this protocol is quite helpful, relatively easy, reliable, and accurate for large datasets as other algorithms need various dependencies and have complex scripts with high chances of error occurrence. Additionally, analysis time will vary depending on the expertise level and the number of RNA-Seq data. This proposed protocol can also be used for a wide range of monocotyledons belonging to the Poaceae plant family.

Comprehensive non-coding RNA analysis reveals specific lncRNA/circRNA-miRNA-mRNA regulatory networks in the cotton response to drought stress

Root and leaf are essential organs of plants in sensing and responding to drought stress. However, comparative knowledge of non-coding RNAs (ncRNAs) of root and leaf tissues in the regulation of drought response in cotton is limited. Here, we used deep sequencing data of leaf and root tissues of drought-resistant and drought-sensitive cotton varieties for identifying miRNAs, lncRNAs and circRNAs. A total of 1531 differentially expressed (DE) ncRNAs was identified, including 77 DE miRNAs, 1393 DE lncRNAs and 61 DE circRNAs. The tissue-specific and variety-specific competing endogenous RNA (ceRNA) networks of DE lncRNA-miRNA-mRNA response to drought were constructed. Furthermore, the novel drought-responsive lncRNA 1 (DRL1), specifically and differentially expressed in root, was verified to positively affect phenotypes of cotton seedlings under drought stress, competitively binding to miR477b with GhNAC1 and GhSCL3. In addition, we also constructed another ceRNA network consisting of 18 DE circRNAs, 26 DE miRNAs and 368 DE mRNAs. Fourteen circRNA were characterized, and a novel molecular regulatory system of circ125- miR7484b/miR7450b was proposed under drought stress. Our findings revealed the specificity of ncRNA expression in tissue- and variety-specific patterns involved in the response to drought stress, and uncovered novel regulatory pathways and potentially effective molecules in genetic improvement for crop drought resistance.

LncRNAs exert indispensable roles in orchestrating the interaction among diverse noncoding RNAs and enrich the regulatory network of plant growth and its adaptive environmental stress response

With the advent of advanced sequencing technologies, non-coding RNAs (ncRNAs) are increasingly pivotal and play highly regulated roles in the modulation of diverse aspects of plant growth and stress response. This includes a spectrum of ncRNA classes, ranging from small RNAs to long non-coding RNAs (lncRNAs). Notably, among these, lncRNAs emerge as significant and intricate components within the broader ncRNA regulatory networks. Here, we categorize ncRNAs based on their length and structure into small RNAs, medium-sized ncRNAs, lncRNAs, and circle RNAs. Furthermore, the review delves into the detailed biosynthesis and origin of these ncRNAs. Subsequently, we emphasize the diverse regulatory mechanisms employed by lncRNAs that are located at various gene regions of coding genes, embodying promoters, 5'UTRs, introns, exons, and 3'UTR regions. Furthermore, we elucidate these regulatory modes through one or two concrete examples. Besides, lncRNAs have emerged as novel central components that participate in phase separation processes. Moreover, we illustrate the coordinated regulatory mechanisms among lncRNAs, miRNAs, and siRNAs with a particular emphasis on the central role of lncRNAs in serving as sponges, precursors, spliceosome, stabilization, scaffolds, or interaction factors to bridge interactions with other ncRNAs. The review also sheds light on the intriguing possibility that some ncRNAs may encode functional micropeptides. Therefore, the review underscores the emergent roles of ncRNAs as potent regulatory factors that significantly enrich the regulatory network governing plant growth, development, and responses to environmental stimuli. There are yet-to-be-discovered roles of ncRNAs waiting for us to explore.

CircRNAs responsive to winter dormancy and spring flushing conditions of tea leaf buds

Circular RNAs (circRNAs) are important regulators of diverse biological processes of plants. However, the evolution and potential functions of circRNAs during winter dormancy and spring bud flushing of tea plant is largely unknown. Using RNA-seq data, a total of 1184 circRNAs were identified in the winter dormant and spring bud flushing leaf samples of tea plants in two different cultivars exhibiting different duration of winter dormancy. A total of 156 circRNAs are found to be differentially expressed and the weighted gene co-expression network (WGCNA) analysis revealed that 22 and 20 differentially expressed-circRNAs (DE-circRNAs) positively correlated with the flushing and dormant leaf traits, respectively, in both the tea cultivars used. Some transcription factors (TFs) viz. MYB, WRKY, ERF, bHLH and several genes related to secondary metabolite biosynthetic pathways are found to co-express with circRNAs. DE-circRNAs also predicted to interact with miRNAs and can regulate phytohormone biosynthesis and various signalling pathways in tea plant. This study uncovers the potential roles of circRNAs to determine winter dormancy and spring bud flushing conditions in tea plants.

Is winter coming? Impact of the changing climate on plant responses to cold temperature

Climate change is causing alterations in annual temperature regimes worldwide. Important aspects of this include the reduction of winter chilling temperatures as well as the occurrence of unpredicted frosts, both significantly affecting plant growth and yields. Recent studies advanced the knowledge of the mechanisms underlying cold responses and tolerance in the model plant Arabidopsis thaliana. However, how these cold-responsive pathways will readjust to ongoing seasonal temperature variation caused by global warming remains an open question. In this review, we highlight the plant developmental programmes that depend on cold temperature. We focus on the molecular mechanisms that plants have evolved to adjust their development and stress responses upon exposure to cold. Covering both genetic and epigenetic aspects, we present the latest insights into how alternative splicing, noncoding RNAs and the formation of biomolecular condensates play key roles in the regulation of cold responses. We conclude by commenting on attractive targets to accelerate the breeding of increased cold tolerance, bringing up biotechnological tools that might assist in overcoming current limitations. Our aim is to guide the reflection on the current agricultural challenges imposed by a changing climate and to provide useful information for improving plant resilience to unpredictable cold regimes.

Evaluation of efficacy of non-coding RNA in abiotic stress management of field crops: Current status and future prospective

Abiotic stresses are responsible for the major losses in crop yield all over the world. Stresses generate harmful ROS which can impair cellular processes in plants. Therefore, plants have evolved antioxidant systems in defence against the stress-induced damages. The frequency of occurrence of abiotic stressors has increased several-fold due to the climate change experienced in recent times and projected for the future. This had particularly aggravated the risk of yield losses and threatened global food security. Non-coding RNAs are the part of eukaryotic genome that does not code for any proteins. However, they have been recently found to have a crucial role in the responses of plants to both abiotic and biotic stresses. There are different types of ncRNAs, for example, miRNAs and lncRNAs, which have the potential to regulate the expression of stress-related genes at the levels of transcription, post-transcription, and translation of proteins. The lncRNAs are also able to impart their epigenetic effects on the target genes through the alteration of the status of histone modification and organization of the chromatins. The current review attempts to deliver a comprehensive account of the role of ncRNAs in the regulation of plants' abiotic stress responses through ROS homeostasis. The potential applications ncRNAs in amelioration of abiotic stresses in field crops also have been evaluated.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Genome-wide identification and characterization of circRNAs in wheat tiller

KEY MESSAGE

This study found that the intergenic circRNAs of wheat were more abundant than those of other plants. More importantly, a circRNA-mediated network associated with tillering was constructed for the first time. Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs with covalently closed circular structures, which play an important role in transcriptional and post-transcriptional regulation. Tiller is an important agronomic trait that determines plant morphological architecture and affects spike number in wheat. However, no studies on the characteristics and functions of circRNAs involved in the regulation of wheat tiller. Here, we performed a genome-wide identification of circRNAs using ribosomal-depleted RNA-seq from wheat tiller of two pairs near-isogenic lines. A total of 686 circRNAs were identified and distributed on 21 chromosomes of wheat, of which 537 were novel circRNAs. Unlike other plants, the majority of these circRNAs (61.8%) were derived from intergenic regions. One circRNA-mediated network associated with tillering was constructed through weighted gene co-expression network analysis, including 323 circRNAs, 117 miRNAs, and 968 mRNAs. GO and pathway enrichment analysis of mRNAs suggested that these circRNAs are involved in cell cycle, ncRNA export from nucleus, developmental process, plant hormone signal transduction, MAPK signaling pathway, RNA degradation. Of these circRNAs, ten circRNAs are associated with known tillering/branching genes in rice or Arabidopsis thaliana, including OsCesA7, EBR1, DTE1, CRD1, LPA1, PAY1, LRK1, OsNR2, OsCCA1, OsBZR1. In summary, we present the first study of the identification and characterization of circRNAs in wheat tiller, and the results suggest these circRNAs associated with tillering could play an important role in wheat tiller formation and development.

Identification and characterization of CircRNA-associated CeRNA networks in moso bamboo under nitrogen stress

BACKGROUND

Nitrogen is a macronutrient element for plant growth and development. Circular RNAs (circRNAs) serve as pivotal regulators for the coordination between nutrient supply and plant demand. Moso bamboo (Phyllostachys edulis) is an excellent plant with fast growth, and the mechanism of the circRNA-target module in response to nitrogen remains unclear.

RESULTS

Deep small RNA sequencing results of moso bamboo seedlings under different concentrations of KNO₃ (N0 = 0 mM, N6 = 6 mM, N18 = 18 mM) were used to identify circRNAs. A total of 549 circRNAs were obtained, of which 309 were generated from corresponding parental coding genes including 66 new ones. A total of 536 circRNA-parent genes were unevenly distributed in 24 scaffolds and were associated with root growth and development. Furthermore, 52 differentially expressed circRNAs (DECs) were obtained, including 24, 33 and 15 DECs from three comparisons of N0 vs. N6, N0 vs. N18 and N6 vs. N18, respectively. Based on integrative analyses of the identified DECs, differentially expressed mRNAs (DEGs), and miRNAs (DEMs), a competitive endogenous RNA (ceRNA) network was constructed, including five DECs, eight DEMs and 32 DEGs. A regulatory module of PeSca_6:12,316,320|12,372,905-novel_miR156-PH02Gene35622 was further verified by qPCR and dual-luciferase reporter assays.

CONCLUSION

The results indicated that circRNAs could participate in multiple biological processes as miRNA sponges, including organ nitrogen compound biosynthesis and metabolic process regulation in moso bamboo. Our results provide valuable information for further study of circRNAs in moso bamboo under fluctuating nitrogen conditions.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Identification, biogenesis, function, and mechanism of action of circular RNAs in plants

Circular RNAs (circRNAs) are a class of single-stranded, closed RNA molecules with unique functions that are ubiquitously expressed in all eukaryotes. The biogenesis of circRNAs is regulated by specific cis-acting elements and trans-acting factors in humans and animals. circRNAs mainly exert their biological functions by acting as microRNA sponges, forming R-loops, interacting with RNA-binding proteins, or being translated into polypeptides or proteins in human and animal cells. Genome-wide identification of circRNAs has been performed in multiple plant species, and the results suggest that circRNAs are abundant and ubiquitously expressed in plants. There is emerging compelling evidence to suggest that circRNAs play essential roles during plant growth and development as well as in the responses to biotic and abiotic stress. However, compared with recent advances in human and animal systems, the roles of most circRNAs in plants are unclear at present. Here we review the identification, biogenesis, function, and mechanism of action of plant circRNAs, which will provide a fundamental understanding of the characteristics and complexity of circRNAs in plants.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

CircRNA: a rising star in plant biology

Circular RNAs (circRNAs) are covalently closed single-stranded RNA molecules, which are widespread in eukaryotic cells. As regulatory molecules, circRNAs have various functions, such as regulating gene expression, binding miRNAs or proteins, and being translated into proteins, which are important for cell proliferation and cell differentiation, individual growth and development, as well as many other biological processes. However, compared with that in animal models, studies of circRNAs in plants lags behind and, particularly, the regulatory mechanisms of biogenesis and molecular functions of plant circRNAs remain elusive. Recent studies have shown that circRNAs are wide spread in plants with tissue- or development-specific expression patterns and are responsive to a variety of environmental stresses. In this review, we summarize these advances, focusing on the regulatory mechanisms of biogenesis, molecular and biological functions of circRNAs, and the methods for investigating circRNAs. We also discuss the challenges and the prospects of plant circRNA studies.

Regulatory non-coding RNA: The core defense mechanism against plant pathogens

Plant pathogens damage crops and threaten global food security. Plants have evolved complex defense networks against pathogens, using crosstalk among various signaling pathways. Key regulators conferring plant immunity through signaling pathways include protein-coding genes and non-coding RNAs (ncRNAs). The discovery of ncRNAs in plant transcriptomes was first considered "transcriptional noise". Recent reviews have highlighted the importance of non-coding RNAs. However, understanding interactions among different types of noncoding RNAs requires additional research. This review attempts to consider how long-ncRNAs, small-ncRNAs and circular RNAs interact in response to pathogenic diseases within different plant species. Developments within genomics and bioinformatics could lead to the further discovery of plant ncRNAs, knowledge of their biological roles, as well as an understanding of their importance in exploiting the recent molecular-based technologies for crop protection.

Nitric oxide, salicylic acid and oxidative stress: Is it a perfect equilateral triangle?

Nitric oxide (NO) is an endogenous free radical involved in the regulation of a wide array of physio-biochemical phenomena in plants. The biological activity of NO directly depend on its cellular concentration which usually changes under stress conditions, it participates in maintaining cellular redox equilibrium and regulating target checkpoints which control switches among development and stress. It is one of the key players in plant signalling and a plethora of evidence supports its crosstalk with other phytohormones. NO and salicylic acid (SA) cooperation is also of great physiological relevance, where NO modulates the immune response by regulating SA linked target proteins i.e., non-expressor of pathogenesis-related genes (NPR-1 and NPR-2) and Group D bZIP (basic leucine zipper domain transcription factor). Many experimental data suggest a functional cooperative role between NO and SA in mitigating the plant oxidative stress which suggests that these relationships could constitute a metabolic "equilateral triangle".

The Intersection of Non-Coding RNAs Contributes to Forest Trees' Response to Abiotic Stress

Non-coding RNAs (ncRNAs) play essential roles in plants by modulating the expression of genes at the transcriptional or post-transcriptional level. In recent years, ncRNAs have been recognized as crucial regulators for growth and development in forest trees, and ncRNAs that respond to various abiotic stresses are now under intense study. In this review, we summarized recent advances in the understanding of abiotic stress-responsive microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) in forest trees. Furthermore, we analyzed the intersection of miRNAs, and epigenetic modified ncRNAs of forest trees in response to abiotic stress. In particular, the abiotic stress-related lncRNA/circRNA-miRNA-mRNA regulatory network of forest trees was explored.

Identification and Characterization of Circular RNAs in Brassica rapa in Response to Plasmodiophora brassicae

Plasmodiophora brassicae is a soil-borne pathogen that attacks the roots of cruciferous plants and causes clubroot disease. CircRNAs are noncoding RNAs, widely existing in plant and animal species. Although knowledge of circRNAs has been updated continuously and rapidly, information about circRNAs in the regulation of clubroot disease resistance is extremely limited in Brassica rapa. Here, Chinese cabbage (BJN 222) containing clubroot resistance genes (CRa) against P. brassicae Pb4 was susceptible to PbE. To investigate the mechanism of cicRNAs responsible for clubroot disease resistance in B. rapa, circRNA-seq was performed with roots of 'BJN 222' at 0, 8, and 23 days post-inoculated (dpi) with Pb4 and PbE. A total of 231 differentially expressed circRNAs were identified between the groups. Based on the differentially expressed circRNAs, the circRNA-miRNA-mRNA network was constructed using the target genes directly or indirectly related to plant resistance. Upregulated novel_circ_000495 suppressed the expression of miR5656-y, leading to the upregulation of Bra026508, which might cause plant resistance. Our results provide new insights into clubroot resistance mechanisms and lay a foundation for further studies exploring complex gene regulation networks in B. rapa.

Low Light Conditions Alter Genome-Wide Profiles of Circular RNAs in Rice Grains during Grain Filling

In animals and plants, circRNAs regulate gene expression and act as sponges that inhibit the activity of microRNAs. This study aimed to determine how specific circRNAs are expressed in rice grains at different stages of grain filling, under normal and low light conditions. We extracted total RNA from rice grains under low and sufficient light conditions. Deep sequencing was performed using circRNA libraries, and bioinformatics tools were used to identify the circRNAs. In addition, we analyzed targeted messenger RNA functions using two databases to predict the processes involved in rice grain development, and we conducted real-time PCR on 15 of the circRNAs as well as Sanger sequencing. During the grain development process, 8015 candidate circRNAs were isolated, among which the number of known circRNAs was 1661. We also found that the number of circRNAs changed with the time of development. Among them, six circRNAs acted as sponges that targeted more than two microRNAs at different stages of development, and these circRNAs showed a regulatory pattern consistent with the transcriptome sequencing results. More circRNA diversity was found under low light treatment compared to normal light. These findings reveal a possible link between circRNA regulation and the expression of the functional genes associated with photosignal-mediated rice grain development.

Expression Characteristics in Roots, Phloem, Leaves, Flowers and Fruits of Apple circRNA

Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in different tissues of apple are currently unknown. A total of 6495 unique circRNAs were identified from roots, phloem, leaves, flowers and fruits; 65.99% of them were intergenic circRNAs. Similar to other plants, tissue-specific expression was also observed for apple circRNAs; only 175 (2.69%) circRNAs were prevalently expressed in all five different tissues, while 1256, 1064, 912, 904 and 1080 circRNAs were expressed only in roots, phloem, leaves, flowers and fruit, respectively. The hosting-genes of circRNAs showed significant differences enriched in COG, GO terms or KEGG pathways in five tissues, suggesting the special functions of circRNAs in different tissues. Potential binding interactions between circRNAs and miRNAs were investigated using TargetFinder; 2989 interactions between 647 circRNAs and 192 miRNA were predicated in the present study. It also predicted that Chr00:18744403|18744580-mdm-miR160 might play an important role in the formation of flowers or in regulating the coloration of flowers, Chr10:6857496|6858910-mdm-miR168 might be involved in response to drought stress in roots, and Chr03:1226434|1277176 may absorb mdm-miR482a-3p and play a major role in disease resistance. Two circRNAs were experimentally analyzed by qRT-PCR with divergent primers, the expression levels were consistent with RNA-seq, which indicates that the RNA-seq datasets were reliable.

Identification and Characterization of circRNAs under Drought Stress in Moso Bamboo (Phyllostachys edulis)

Circular RNAs (circRNAs) are a class of endogenous noncoding RNAs formed by 3′-5′ ligation during splicing. They play an important role in the regulation of transcription and miRNA in eukaryotes. Drought is one of the detrimental abiotic stresses that limit plant growth and productivity. How circRNAs influence the response to drought stress in moso bamboo (Phyllostachys edulis) remains elusive. In this study, we investigate the expression pattern of circRNAs in moso bamboo at 6 h, 12 h, 24 h and 48 h after drought treatment by deep sequencing and bioinformatics analysis and identify 4931 circRNAs, 52 of which are differentially expressed (DEcircRNAs) in drought-treated and untreated moso bamboo. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the host genes that generate the DEcircRNAs indcate that these DEcircRNAs are predicted to be involved in biochemical processes in response to drought, such as ubiquitin-mediated proteolysis, calcium-dependent protein kinase phosphorylation, amino acid biosynthesis and plant hormone signal transduction including abscisic acid. In addition, some circRNAs are shown to act as sponges for 291 miRNAs. Taken together, our results characterize the transcriptome profiles of circRNAs in drought responses and provide new insights into resistance breeding of moso bamboo.

Exonic Circular RNAs Are Involved in Arabidopsis Immune Response Against Bacterial and Fungal Pathogens and Function Synergistically with Corresponding Linear RNAs

Circular RNAs (circRNAs) are a group of covalently closed RNAs, and their biological function is largely unknown. In this study, we focused on circRNAs that are generated from exon back-splicing (exonic circRNAs). The linear RNA counterparts encode functional proteins so that we can compare and investigate the relationship between circular and linear RNAs. We compared circRNA expression profiles between untreated and Pseudomonas syringae-infected Arabidopsis and identified and experimentally validated differentially expressed exonic circRNAs by multiple approaches. We found that exonic circRNAs are preferentially enriched in biological processes that associate with biotic and abiotic stress responses. We discovered that circR194 and circR4022 are involved in plant response against P. syringae infection, whereas circR11208 is involved in response against Botrytis cinerea infection. Intriguingly, our results indicate that these exonic circRNAs function synergistically with their corresponding linear RNAs. Furthermore, circR4022 and circR11208 also play substantial roles in Arabidopsis tolerance to salt stress. This study extends our understanding of the molecular functions of plant circRNAs.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

Identification and Characterization of Rice Circular RNAs Responding to Xanthomonas oryzae pv. oryzae Invasion

Emerging roles of circular RNAs (circRNAs) in various biological processes have advanced our knowledge of transcriptional and posttranscriptional gene regulation. To date, no research has been conducted to explore their roles in the rice-Xanthomonas oryzae pv. oryzae interaction. Therefore, we identified 3,517 circRNAs from rice leaves infected with the highly virulent X. oryzae pv. oryzae strain PXO99^A by using rRNA depleted RNA sequencing technique coupled with the CIRI2 and CIRCexplorer2 pipeline. Characterization analyses showed that these circRNAs were distributed across the whole genome of rice, and most circRNAs arose from exons (85.13%), ranged from 200 to 1,000 bp, and were with a noncanonical GT/AG (including CT/AC equivalent) splicing signal. Functional annotation and enrichment analysis of the host genes that produced the differentially expressed circRNAs (DEcircRNAs) suggested that these identified circRNAs might play an important role in reprogramming rice responses to PXO99^A invasion, mainly by mediating photorespiration and chloroplast, peroxisome, and diterpenoid biosynthesis. Moreover, 31 DEcircRNAs were predicted to act as microRNA decoys in rice. The expression profile of four DEcircRNAs were validated by quantitative real-time PCR with divergent primers, and the back-splicing sites of seven DEcircRNAs were verified by PCR analysis and Sanger sequencing. Collectively, these results inferred a potential functional role of circRNAs in the regulation of rice immunity and provide novel clues about the molecular mechanisms of rice-PXO99^A interaction.

Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza

Circular RNAs (circRNAs) play various roles in cellular functions. However, no studies have been reported on the potential involvement of circRNAs in the biosynthesis of secondary metabolites in plants. Here, we performed a genome-wide discovery of circRNAs from root, stem and leaf samples of Salvia miltiorrhiza using RNA-Seq. We predicted a total of 2,476 circRNAs with at least two junction reads using circRNA_finder and CIRI, of which 2,096, 151 and 229 were exonic, intronic and intergenic circRNAs, respectively. Sequence similarity analysis showed that 294 out of 2,476 circRNAs were conserved amongst multiple plants. Of the 55 predicted circRNAs, 31 (56%) were validated successfully by PCR and Sanger sequencing using convergent and divergent primer pairs. Alternative circularisation analysis showed that most parental genes produced two circRNAs. Functional enrichment analyses of the parental genes showed that the primary metabolism pathways were significantly enriched, particularly the carbon metabolism. Differential expression analysis showed that the expression profiles of circRNAs were tissue-specific. Co-expression analysis showed 275 circRNAs, and their parental genes had significantly positive correlations. However, 14 had significantly negative correlations. Weighted gene co-expression network analysis showed that nine circRNAs were co-expressed with four modules of protein-coding genes. Next, we found 416 exonic circRNAs with miRNA-binding sites, suggesting possible interactions between circRNAs and miRNAs. Lastly, we found six validated circRNAs, namely, SMscf2473-46693-46978, SMscf3091-29256-29724, SMscf16-111773-112193, SMscf432-13232-13866, SMscf7007-10563-10888 and SMscf1730-1749-2013, which were originated from the genes involved in the biosynthesis of secondary metabolites. Their parental genes were acetyl-CoA C-acetyltransferase 1 (SmAACT1), 1-deoxy-d-xylulose-5-phosphate synthase 2 (SmDXS2), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1 (SmHDR1), kaurene synthase-like 2 (SmKSL2), DWF4 and CYP88A3, respectively. In particular, the correlation coefficient of SMscf2473-46693-46978 and SmDXS2 gene was 0.86 (p = 0.003), indicating a potential interaction between this pair of circRNA and its parent gene. Our results provided the first comprehensive catalogue of circRNAs in S. miltiorrhiza and identified one circRNA that might play important roles in the biosynthesis of secondary metabolites.

Genome-Wide Characterization of Salt-Responsive miRNAs, circRNAs and Associated ceRNA Networks in Tomatoes

Soil salinization is a major environmental stress that causes crop yield reductions worldwide. Therefore, the cultivation of salt-tolerant crops is an effective way to sustain crop yield. Tomatoes are one of the vegetable crops that are moderately sensitive to salt stress. Global market demand for tomatoes is huge and growing. In recent years, the mechanisms of salt tolerance in tomatoes have been extensively investigated; however, the molecular mechanism through which non-coding RNAs (ncRNAs) respond to salt stress is not well understood. In this study, we utilized small RNA sequencing and whole transcriptome sequencing technology to identify salt-responsive microRNAs (miRNAs), messenger RNAs (mRNAs), and circular RNAs (circRNAs) in roots of M82 cultivated tomato and Solanum pennellii (S. pennellii) wild tomato under salt stress. Based on the theory of competitive endogenous RNA (ceRNA), we also established several salt-responsive ceRNA networks. The results showed that circRNAs could act as miRNA sponges in the regulation of target mRNAs of miRNAs, thus participating in the response to salt stress. This study provides insights into the mechanisms of salt tolerance in tomatoes and serves as an effective reference for improving the salt tolerance of salt-sensitive cultivars.

Expression Patterns and Regulation of Non-Coding RNAs during Synthesis of Cellulose in Eucalyptus grandis Hill

Cellulose, an essential structural component in the plant cell wall and a renewable biomass resource, plays a significant role in nature. Eucalyptus’s excellent timber tree species (including Eucalyptus grandis Hill) provide many raw materials for the paper and wood industries. The synthesis of cellulose is a very complex process involving multiple genes and regulated by various biological networks. However, research on regulating associated genes and non-coding RNAs during cellulose synthesis in E. grandis remains lacking. In this study, the wood anatomical characteristics and chemical indexes of E. grandis were analyzed by taking three different parts (diameter at breast height (DBH), middle and upper part of the trunk) from the main stem of E. grandis as raw materials. The role of non-coding RNAs (Long non-coding RNA, lncRNA; Micro RNA, miRNA; Circle RNA, circRNA) on regulating candidate genes was presented, and the network map of ceRNA (Competing endogenous RNA) regulation during wood cellulose biosynthesis of E. grandis was constructed. The transcriptome sequencing of nine samples obtained from the trunk of the immature xylem in E. grandis at DBH, middle and upper parts had a 95.81 G clean reading, 57,480 transcripts, 7365 lncRNAs, and 5180 circRNAs. Each sample had 172–306 known miRNAs and 1644–3508 new miRNAs. A total of 190 DE-lncRNAs (Differentially expressed long non-coding RNAs), 174 DE-miRNAs (Differentially expressed micro RNAs), and 270 DE-circRNAs (Differentially expressed circle RNAs) were obtained by comparing transcript expression levels. Four lncRNAs and nine miRNAs were screened out, and the ceRNA regulatory network was constructed. LncRNA1 and lncRNA4 regulated the genes responsible for cellulose synthesis in E. grandis, which were overexpressed in 84K (Populus Alba × Populus glandulosa) poplar. The cellulose and lignin content in lncRNA4-oe were significantly higher than wild type 84K poplar and lncRNA1-oe. The average plant height, middle and basal part of the stem diameter in lncRNA4-oe were significantly higher than the wild type. However, there was no significant difference between the growth of lncRNA1-oe and the wild type. Further studies are warranted to explore the molecular regulatory mechanism of cellulose biosynthesis in Eucalyptus species.

GreenCircRNA: a database for plant circRNAs that act as miRNA decoys

Circular RNAs (circRNAs) are endogenous non-coding RNAs that form a covalently closed continuous loop, are widely distributed and play important roles in a series of developmental processes. In plants, an increasing number of studies have found that circRNAs can regulate plant metabolism and are involved in plant responses to biotic or abiotic stress. Acting as miRNA decoys is a critical way for circRNAs to perform their functions. Therefore, we developed GreenCircRNA—a database for plant circRNAs acting as miRNA decoys that is dedicated to providing a plant-based platform for detailed exploration of plant circRNAs and their potential decoy functions. This database includes over 210 000 circRNAs from 69 species of plants; the main data sources of circRNAs in this database are NCBI, EMBL-EBI and Phytozome. To investigate the function of circRNAs as competitive endogenous RNAs, the possibility of circRNAs from 38 plants to act as miRNA decoys was predicted. Moreover, we provide basic information for the circRNAs in the database, including their locations, host genes and relative expression levels, as well as full-length sequences, host gene GO (Gene Ontology) numbers and circRNA visualization. GreenCircRNA is the first database for the prediction of circRNAs that act as miRNA decoys and contains the largest number of plant species.

Database URL: http://greencirc.cn

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.

Efficient deletion of multiple circle RNA loci by CRISPR-Cas9 reveals Os06circ02797 as a putative sponge for OsMIR408 in rice

CRISPR-Cas9 is an emerging genome editing tool for reverse genetics in plants. However, its application for functional study of non-coding RNAs in plants is still at its infancy. Despite being a major class of non-coding RNAs, the biological roles of circle RNAs (circRNAs) remain largely unknown in plants. Previous plant circRNA studies have focused on identification and annotation of putative circRNAs, with their functions largely uninvestigated by genetic approaches. Here, we applied a multiplexed CRISPR-Cas9 strategy to efficiently acquire individual null mutants for four circRNAs in rice. We showed each of these rice circRNA loci (Os02circ25329, Os06circ02797, Os03circ00204 and Os05circ02465) can be deleted at 10% or higher efficiency in both protoplasts and stable transgenic T0 lines. Such high efficiency deletion enabled the generation of circRNA null allele plants without the CRISPR-Cas9 transgene in the T1 generation. Characterization of the mutants reveals these circRNAs' participation in salt stress response during seed germination and in particular the Os05circ02465 null mutant showed high salt tolerance. Notably, the seedlings of the Os06circ02797 mutant showed rapid growth phenotype after seed germination with the seedlings containing higher chlorophyll A/B content. Further molecular and computational analyses suggested a circRNA-miRNA-mRNA regulatory network where Os06circ02797 functions to bind and sequester OsMIR408, an important and conserved microRNA in plants. This study not only presents genetic evidence for the first time in plants that certain circRNAs may serve as sponges to negatively regulate miRNAs, a phenomenon previously demonstrated in mammalian cells, but also provides important insights for improving agronomic traits through gene editing of circRNA loci in crops.

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.

Molecular and Hormonal Mechanisms Regulating Fleshy Fruit Ripening

This article focuses on the molecular and hormonal mechanisms underlying the control of fleshy fruit ripening and quality. Recent research on tomato shows that ethylene, acting through transcription factors, is responsible for the initiation of tomato ripening. Several other hormones, including abscisic acid (ABA), jasmonic acid (JA) and brassinosteroids (BR), promote ripening by upregulating ethylene biosynthesis genes in different fruits. Changes to histone marks and DNA methylation are associated with the activation of ripening genes and are necessary for ripening initiation. Light, detected by different photoreceptors and operating through ELONGATED HYPOCOTYL 5(HY5), also modulates ripening. Re-evaluation of the roles of 'master regulators' indicates that MADS-RIN, NAC-NOR, Nor-like1 and other MADS and NAC genes, together with ethylene, promote the full expression of genes required for further ethylene synthesis and change in colour, flavour, texture and progression of ripening. Several different types of non-coding RNAs are involved in regulating expression of ripening genes, but further clarification of their diverse mechanisms of action is required. We discuss a model that integrates the main hormonal and genetic regulatory interactions governing the ripening of tomato fruit and consider variations in ripening regulatory circuits that operate in other fruits.

Emerging roles of centromeric RNAs in centromere formation and function

BACKGROUND

Centromeres are specialized chromosomal domains involved in kinetochore formation and faithful chromosome segregation. Despite a high level of functional conservation, centromeres are not identified by DNA sequences, but by epigenetic means. Universally, centromeres are typically formed on highly repetitive DNA, which were previously considered to be silent. However, recent studies have shown that transcription occurs in this region, known as centromeric-derived RNAs (cenRNAs). CenRNAs that contribute to fundamental aspects of centromere function have been recently investigated in detail. However, the distribution, behavior and contributions of centromeric transcripts are still poorly understood.

OBJECTIVE

The aim of this article is to provide an overview of the roles of cenRNAs in centromere formation and function.

METHODS

We describe the structure and DNA sequence of centromere from yeast to human. In addition, we briefly introduce the roles of cenRNAs in centromere formation and function, kinetochore structure, accurate chromosome segregation, and pericentromeric heterochromatin assembly. Centromeric circular RNAs (circRNAs) and R-loops are rising stars in centromere function. CircRNAs have been successfully identified in various species with the assistance of high-throughput sequencing and novel computational approaches for non-polyadenylated RNA transcripts. Centromeric R-loops can be identified by the single-strand DNA ligation-based library preparation technique. But the molecular features and function of these centromeric R-loops and circRNAs are still being investigated.

CONCLUSION

In this review, we summarize recent findings on the epigenetic regulation of cenRNAs across species, which would provide useful information about cenRNAs and interesting hints for further studies.

PCirc: random forest-based plant circRNA identification software

Background

Circular RNA (circRNA) is a novel type of RNA with a closed-loop structure. Increasing numbers of circRNAs are being identified in plants and animals, and recent studies have shown that circRNAs play an important role in gene regulation. Therefore, identifying circRNAs from increasing amounts of RNA-seq data is very important. However, traditional circRNA recognition methods have limitations. In recent years, emerging machine learning techniques have provided a good approach for the identification of circRNAs in animals. However, using these features to identify plant circRNAs is infeasible because the characteristics of plant circRNA sequences are different from those of animal circRNAs. For example, plants are extremely rich in splicing signals and transposable elements, and their sequence conservation in rice, for example is far less than that in mammals. To solve these problems and better identify circRNAs in plants, it is urgent to develop circRNA recognition software using machine learning based on the characteristics of plant circRNAs.

Results

In this study, we built a software program named PCirc using a machine learning method to predict plant circRNAs from RNA-seq data. First, we extracted different features, including open reading frames, numbers of k-mers, and splicing junction sequence coding, from rice circRNA and lncRNA data. Second, we trained a machine learning model by the random forest algorithm with tenfold cross-validation in the training set. Third, we evaluated our classification according to accuracy, precision, and F1 score, and all scores on the model test data were above 0.99. Fourth, we tested our model by other plant tests, and obtained good results, with accuracy scores above 0.8. Finally, we packaged the machine learning model built and the programming script used into a locally run circular RNA prediction software, Pcirc (https://github.com/Lilab-SNNU/Pcirc).

Conclusion

Based on rice circRNA and lncRNA data, a machine learning model for plant circRNA recognition was constructed in this study using random forest algorithm, and the model can also be applied to plant circRNA recognition such as Arabidopsis thaliana and maize. At the same time, after the completion of model construction, the machine learning model constructed and the programming scripts used in this study are packaged into a localized circRNA prediction software Pcirc, which is convenient for plant circRNA researchers to use.

Generation of Transgenic Rice Expressing CircRNA and Its Functional Characterization

Circular RNA (CircRNA) is yet another vital addition to the noncoding RNA family. They are mainly derived by fusion of downstream 3' splice donor with upstream 5' splice acceptor by a noncanonical form of alternative splicing mechanism called backsplicing. An array of functional aspects of these circRNAs has been reported in animal systems. However, functional investigation of circRNA in plants is very limited. In this chapter, we described a methodological outline to study the circRNA biogenesis and to characterize its function(s). Sequence of a newly identified Oryza sativa Indica circRNA flanked by complementary repeat sequences of a rice intron was assembled to yield a circRNA expression cassette. This cassette can be cloned into any plant expression vector which has a suitable promoter (CaMV 35S or ubiquitin promoter) and terminator, and can be used for any circRNA-mediated functional studies. Subsequent agroinfection of rice calli with this cassette yielded circRNA expressing transgenic plants. These transgenic plants were used to establish a correlation between the expressing circRNA, parental gene, and interacting miRNAs. Moreover, effect of circRNA overexpression on plant phenotype under various stress conditions can be studied using these transgenic plants. Also, RNA pull-down assay can be performed to identify the circRNA interacting proteins and the expression of these RBPs can also be studied from these transgenic plants.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

Identification, characterization and expression profiling of circular RNAs in the early cotton fiber developmental stages

Circular RNA is one of the endogenous non-coding RNAs with a covalently closed loop structure and largely involved in regulating gene expression. However, the abundance of circular RNAs and their regulatory functions during the early stages of fiber development are still not known. In this work, we conducted high-throughput sequencing of the Ligonlintless-1 and its wild-type at 0 DPA, 8 DPA and stem. A total of 2811 circular RNAs were identified and unevenly distributed across cotton chromosomes. We found 34, 142 and 27 circular RNAs were differentially expressed between Ligonlintless-1 and wild-type at 0 DPA, 8 DPA and stem, respectively. Both circular RNA-microRNA-mRNA network and MeJA treatment results in Ligonlintless-1 and wild-type might provide a strong indication of four circular RNAs and ghr_miR169b being important biological molecular associating with fiber development. The results provide new insight into the putative molecular function of circular RNAs in the regulation of fiber development.

Genome-Wide Identification of Circular RNAs in Response to Low-Temperature Stress in Tomato Leaves

Abiotic stress adversely inhibits the growth and development of plants, by changing the expression of multiple genes. Circular RNAs (circRNAs), as a class of non-coding RNAs, function in transcriptional and posttranscriptional regulation. Yet, the involvement of circRNAs in abiotic stress response is rarely reported. In this study, the participation and function of circRNAs in low-temperature (LT)-induced stress response were investigated in tomato leaves. We generated genome-wide profiles of circRNAs and mRNAs in tomato leaves grown at 25°C room temperature (RT) and 12°C LT. Our results show that 1,830 circRNAs were identified in tomato leaves in both RT and LT treatments, among which 1,759 were differentially induced by the LT treatment. We find that the identified circRNAs are mainly located at exons of genes, but less distributed at introns of genes or intergenic regions. Our results suggest that there are 383 differentially expressed circRNAs predicted to function as putative sponges of 266 miRNAs to target 4,476 mRNAs in total. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis assays indicate that multiple pathways were enriched in both differentially expressed genes induced by LT and parental genes of differentially expressed circRNAs induced by LT, revealing the key functions of circRNAs and the corresponding targeted genes in response to LT stress. Our results suggest that circRNAs may be involved in regulating metabolism (i.e., carbohydrate, amino acid, lipid, and energy), signal transduction, and environmental adaptation-related pathways and that these circRNAs were predicted to regulate the expression of transcription factors, genes in signal transduction pathways, and genes related to the Ca²⁺ channel through targeting the corresponding proteins, such as WRKY, NAC, cytochrome P450, and calmodulin binding protein. Taken together, our study uncovers that multiple circRNAs are isolated and differently regulated in response to LT stress and provides the resource and potential networks of circRNA–miRNA–mRNA under LT stress for further investigations in tomato leaves.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.