Plant Circular RNAs (circRNAs): Transcriptional Regulation Beyond miRNAs in Plants

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.

Multiple novel functions of circular RNAs in diabetes mellitus

Circular RNAs (circRNAs), as an emerging group of non-coding RNAs (ncRNAs), have received the attention given evidence indicating that these novel ncRNAs are implicated in various biological processes. Due to the absence of 5' and 3' ends in circ-RNAs, their two ends are covalently bonded together, and they are synthesised from pre-mRNAs in a process called back-splicing, which makes them more stable than linear RNAs. There is accumulating evidence showing that circRNAs play a critical role in the pathogenesis of diabetes mellitus (DM). Moreover, it has been indicated that dysregulation of circRNAs has made them promising diagnostic biomarkers for the detection of DM. Recently, increasing attention has been paid to investigate the mechanisms underlying the DM process. It has been demonstrated that there is a strong correlation between the expression of circRNAs and DM. Hence, our aim is to discuss the crosstalk between circRNAs and DM and its complications.

Identification and validation of coding and non-coding RNAs involved in high-temperature-mediated seed dormancy in common wheat

INTRODUCTION

Seed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusiveSeed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusive.

METHODS

Here, the wheat landrace 'Waitoubai' with strong SD and PHS resistance was treated with HT from 21 to 35 days post anthesis (DPA). Then, the seeds under HT and normal temperature (NT) environments were collected at 21 DPA, 28 DPA, and 35 DPA and subjected to whole-transcriptome sequencing.

RESULTS

The phenotypic data showed that the seed germination percentage significantly increased, whereas SD decreased after HT treatment compared with NT, consistent with the results of previous studies. In total, 5128 mRNAs, 136 microRNAs (miRNAs), 273 long non-coding RNAs (lncRNAs), and 21 circularRNAs were found to be responsive to HT, and some of them were further verified through qRT-PCR. In particular, the known gibberellin (GA) biosynthesis gene TaGA20ox1 (TraesCS3D02G393900) was proved to be involved in HT-mediated dormancy by using the EMS-mutagenized wheat cultivar Jimai 22. Similarly, a novel gene TaCDPK21 (TraesCS7A02G267000) involved in the calcium signaling pathway was validated to be associated with HT-mediated dormancy by using the EMS mutant. Moreover, TaCDPK21 overexpression in Arabidopsis and functional complementarity tests supported the negative role of TaCDPK21 in SD. We also constructed a co-expression regulatory network based on differentially expressed mRNAs, miRNAs, and lncRNAs and found that a novel miR27319 was located at a key node of this regulatory network. Subsequently, using Arabidopsis and rice lines overexpressing miR27319 precursor or lacking miR27319 expression, we validated the positive role of miR27319 in SD and further preliminarily dissected the molecular mechanism of miR27319 underlying SD regulation through phytohormone abscisic acid and GA biosynthesis, catabolism, and signaling pathways.

DISCUSSION

These findings not only broaden our understanding of the complex regulatory network of HT-mediated dormancy but also provide new gene resources for improving wheat PHS resistance to minimize PHS damage by using the molecular pyramiding approach.

Hypoxia-induced circADAMTS6 in a TDP43-dependent manner accelerates glioblastoma progression via ANXA2/ NF-κB pathway

Circular RNAs (circRNAs) play important roles in the malignant progression of tumours. Herein, we identified an unreported circRNA (hsa-circ-0072688, also named circADAMTS6) that is specifically upregulated in the hypoxic microenvironment of glioblastoma and closely correlated with poor prognosis of gliblastoma patients. We found that circADAMTS6 promotes the malignant progression of glioblastoma by promoting cell proliferation and inhibiting apoptosis. Mechanistically, the hypoxic tumour microenvironment upregulates circADAMTS6 expression through transcription factor activator protein 1 (AP-1) and RNA-binding protein TAR DNA-binding protein 43 (TDP43). Moreover, circADAMTS6 accelerates glioblastoma progression by recruiting and stabilising annexin A2 (ANXA2) in a proteasomes-dependent manner. Furthermore, we found T-5224 (AP-1 inhibitor) treatment induces downregulation of circADAMTS6 and then inhibits tumour growth. In conclusion, our findings highlight the important role of the circADAMTS6/ANXA2 axis based on hypoxic microenvironment in glioblastoma progression, as well as its regulation in NF-κB pathway. Targeting circADAMTS6 is thus expected to become a novel therapeutic strategy for glioblastoma.

Expression profile of circular RNAs in patients with unstable angina

BACKGROUND

Circular RNAs (circRNAs) are a new class of noncoding RNAs, which interfere with gene transcription by absorbing microRNAs (miRNAs).

OBJECTIVE

The expression profile and roles of circRNAs in unstable angina (UA) patients remains unclear.

METHODS

An initial screening of circRNA expression by microarray analysis was performed using blood samples from three pairs of UA patients and matched healthy individuals. The differential expression of the chosen six circRNAs from the results of the microarray analysis was validated by quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

The microarray results demonstrated that some circRNAs were markedly different in UA patients, when compared with matched healthy individuals. In these UA patients, 22 circRNAs were upregulated and six circRNAs were downregulated when a P-value of < 0.05 was considered as a cut-off level and the fold change was > 1.5. Among the six circRNAs chosen for further analysis, qRT-PCR identified that five of these were upregulated, and the remaining circRNA was downregulated. By comparing the outcome of the six candidate circRNAs between the circRNAs microarray assay and RT-PCR validation, it was found that four circRNAs (hsa_circ_0002229, hsa_circ_0005580, hsa_circ_0046667, and hsa_circ_0001451) had the same variation trend.

CONCLUSION

The present study provided the expression profile of circRNAs in UA patients. Moreover, some circRNAs have the potential to be biomarkers for the detection of UA patients. Further studies with a larger population will focus on hsa_circ_0002229, hsa_circ_0005580, hsa_circ_0046667 and hsa_circ_0001451.

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.

Characterization of the circRNA–miRNA–mRNA Network to Reveal the Potential Functional ceRNAs Associated With Dynamic Changes in the Meat Quality of the Longissimus Thoracis Muscle in Tibetan Sheep at Different Growth Stages

Circular RNAs (circRNAs) have a regulatory role in animal skeletal muscle development. In this study, RNA sequencing was performed to reveal the temporal regularity of circRNA expression and the effect of the circRNA–miRNA–mRNA ceRNA regulatory network on the meat quality of longissimus thoracis (LT) muscle in Tibetan sheep at different growth stages (4 months old, 4 m; 1.5 years old, 1.5 y; 3.5 years old, 3.5 y; 6 years old, 6 y). There were differences in the carcass performance and meat quality of Tibetan sheep at different ages. Especially, the meat tenderness significantly decreased (p &lt; 0.05) with the increase of age. GO functional enrichment indicated that the source genes of the DE circRNAs were mainly involved in the protein binding, and myofibril and organelle assembly. Moreover, there was a significant KEGG enrichment in the adenosine 5′-monophosphate (AMP)-activated protein kinase (AMPK) signaling pathway, as well as the calcium signaling pathway, regulating the pluripotency of the stem cells. The circRNA–miRNA–mRNA ceRNA interaction network analysis indicated that circRNAs such as circ_000631, circ_000281, and circ_003400 combined with miR-29-3p and miR-185-5p regulate the expression of LEP, SCD, and FASN related to the transformation of muscle fiber types in the AMPK signaling pathway. The oxidized muscle fibers were transformed into the glycolytic muscle fibers with the increase of age, the content of intramuscular fat (IMF) was lowered, and the diameter of the muscle fiber was larger in the glycolytic muscle fibers, ultimately increasing the meat tenderness. The study revealed the role of the circRNAs in the transformation of skeletal muscle fiber types in Tibetan sheep and its influence on meat quality. It improves our understanding of the role of circRNAs in Tibetan sheep muscle development.

Biogenesis, Functions, Interactions, and Resources of Non-Coding RNAs in Plants

Plant transcriptomes encompass a large number of functional non-coding RNAs (ncRNAs), only some of which have protein-coding capacity. Since their initial discovery, ncRNAs have been classified into two broad categories based on their biogenesis and mechanisms of action, housekeeping ncRNAs and regulatory ncRNAs. With advances in RNA sequencing technology and computational methods, bioinformatics resources continue to emerge and update rapidly, including workflow for in silico ncRNA analysis, up-to-date platforms, databases, and tools dedicated to ncRNA identification and functional annotation. In this review, we aim to describe the biogenesis, biological functions, and interactions with DNA, RNA, protein, and microorganism of five major regulatory ncRNAs (miRNA, siRNA, tsRNA, circRNA, lncRNA) in plants. Then, we systematically summarize tools for analysis and prediction of plant ncRNAs, as well as databases. Furthermore, we discuss the silico analysis process of these ncRNAs and present a protocol for step-by-step computational analysis of ncRNAs. In general, this review will help researchers better understand the world of ncRNAs at multiple levels.

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 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.

IMS-CDA: Prediction of CircRNA-Disease Associations From the Integration of Multisource Similarity Information With Deep Stacked Autoencoder Model

Emerging evidence indicates that circular RNA (circRNA) has been an indispensable role in the pathogenesis of human complex diseases and many critical biological processes. Using circRNA as a molecular marker or therapeutic target opens up a new avenue for our treatment and detection of human complex diseases. The traditional biological experiments, however, are usually limited to small scale and are time consuming, so the development of an effective and feasible computational-based approach for predicting circRNA-disease associations is increasingly favored. In this study, we propose a new computational-based method, called IMS-CDA, to predict potential circRNA-disease associations based on multisource biological information. More specifically, IMS-CDA combines the information from the disease semantic similarity, the Jaccard and Gaussian interaction profile kernel similarity of disease and circRNA, and extracts the hidden features using the stacked autoencoder (SAE) algorithm of deep learning. After training in the rotation forest (RF) classifier, IMS-CDA achieves 88.08% area under the ROC curve with 88.36% accuracy at the sensitivity of 91.38% on the CIRCR2Disease dataset. Compared with the state-of-the-art support vector machine and K -nearest neighbor models and different descriptor models, IMS-CDA achieves the best overall performance. In the case studies, eight of the top 15 circRNA-disease associations with the highest prediction score were confirmed by recent literature. These results indicated that IMS-CDA has an outstanding ability to predict new circRNA-disease associations and can provide reliable candidates for biological experiments.

Advance in circular RNA modulation effects of heart failure

CircRNA (circular RNA) is a kind of closed circular structure of noncoding RNA molecules without 5′ hat structure and 3′ polyA, mainly located in the cytoplasm or stored in exosomes. It is not affected by RNA exonuclease, so it's stable and hard to be degraded. Proved to be widespread in a variety of eukaryotes, most circRNAs are cyclized by exons, some are lasso structures formed by intron cyclization. Recently, circRNAs have been demonstrated to play crucial roles in cardiomyocyte hypertrophy, fibrosis, autophagy and apoptosis, participating in the development of heart failure. There is increasing evidence that circRNAs may be a novel target for the treatment of heart failure.

Regulatory non-coding RNAs: a new frontier in regulation of plant biology

Beyond the most crucial roles of RNA molecules as a messenger, ribosomal, and transfer RNAs, the regulatory role of many non-coding RNAs (ncRNAs) in plant biology has been recognized. ncRNAs act as riboregulators by recognizing specific nucleic acid targets through homologous sequence interactions to regulate plant growth, development, and stress responses. Regulatory ncRNAs, ranging from small to long ncRNAs (lncRNAs), exert their control over a vast array of biological processes. Based on the mode of biogenesis and their function, ncRNAs evolved into different forms that include microRNAs (miRNAs), small interfering RNAs (siRNAs), miRNA variants (isomiRs), lncRNAs, circular RNAs (circRNAs), and derived ncRNAs. This article explains the different classes of ncRNAs and their role in plant development and stress responses. Furthermore, the applications of regulatory ncRNAs in crop improvement, targeting agriculturally important traits, have been discussed.

Circad: a comprehensive manually curated resource of circular RNA associated with diseases

Circular RNAs (circRNAs) are unique transcript isoforms characterized by back splicing of exon ends to form a covalently closed loop or circular conformation. These transcript isoforms are now known to be expressed in a variety of organisms across the kingdoms of life. Recent studies have shown the role of circRNAs in a number of diseases and increasing evidence points to their potential application as biomarkers in these diseases. We have created a comprehensive manually curated database of circular RNAs associated with diseases. This database is available at URL http://clingen.igib.res.in/circad/. The Database lists more than 1300 circRNAs associated with 150 diseases and mapping to 113 International Statistical Classification of Diseases (ICD) codes with evidence of association linked to published literature. The database is unique in many ways. Firstly, it provides ready-to-use primers to work with, in order to use circRNAs as biomarkers or to perform functional studies. It additionally lists the assay and PCR primer details including experimentally validated ones as a ready reference to researchers along with fold change and statistical significance. It also provides standard disease nomenclature as per the ICD codes. To the best of our knowledge, circad is the most comprehensive and updated database of disease associated circular RNAs.

Availability: http://clingen.igib.res.in/circad/

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

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

The role of circular RNAs in viral infection and related diseases

•

CircRNAs are formed by pre-mRNA through “back-splicing”.

•

CircRNAs regulate host immune response and virus replication.

•

CircRNAs have potential as diagnostic markers or treatment targets for viral infection.

Circular RNAs (circRNAs) are a class of non-coding RNAs with a special covalently closed circular structure, which is formed by precursor mRNA (pre-mRNA) through “back-splicing”. CircRNAs are more stable than linear RNAs because they are resistant to exoribonucleases. Viral infections often cause abnormal expression of circRNAs, which could serve as novel biomarkers for the diagnosis of viral infections by detecting specific circRNAs in cells, body fluids, or tissues. CircRNAs also play a critical role in regulating host immune response and virus replication. Here, we reviewed the production and function of circRNAs, mainly focusing on their regulation on virus infection, to provide novel insights into the potential role of circRNAs as diagnostic marker or treatment targets for viral infection.

Transcriptome-wide identification of novel circular RNAs in soybean in response to low-phosphorus stress

Low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Circular RNAs (circRNAs) are novel noncoding RNAs that play a crucial role in plant responses to abiotic stress. However, how LP stress mediates the biogenesis of circRNAs in soybean remains unclear. Here, to explore the response mechanisms of circRNAs to LP stress, the roots of two representative soybean genotypes with different P-use efficiency, Bogao (a LP-sensitive genotype) and Nannong 94156 (a LP-tolerant genotype), were used for the construction of RNA sequencing (RNA-seq) libraries and circRNA identification. In total, 371 novel circRNA candidates, including 120 significantly differentially expressed (DE) circRNAs, were identified across different P levels and genotypes. More DE circRNAs were significantly regulated by LP stress in Bogao than in NN94156, suggesting that the tolerant genotype was less affected by LP stress than the sensitive genotype was; in other words, NN94156 may have a better ability to maintain P homeostasis under LP stress. Moreover, a positive correlation was observed between the expression patterns of P stress-induced circRNAs and their circRNA-host genes. Gene Ontology (GO) enrichment analysis of these circRNA-host genes and microRNA (miRNA)-targeted genes indicated that these DE circRNAs were involved mainly in defense responses, ADP binding, nucleoside binding, organic substance catabolic processes, oxidoreductase activity, and signal transduction. Together, our results revealed that LP stress can significantly alter the genome-wide profiles of circRNAs and indicated that the regulation of circRNAs was both genotype and environment specific in response to LP stress. LP-induced circRNAs might provide a rich resource for LP-responsive circRNA candidates for future studies.

Circular RNAs as biomarkers and therapeutic targets in environmental chemical exposure-related diseases

Chemical contamination in the environment is known to cause abnormal circular RNA (circRNA) expression through multiple exposure routes; yet, the underlying molecular mechanisms remain unclear. Non-coding RNAs (ncRNAs), especially circRNAs, play important roles in epigenetic regulation and disease pathogenesis; however, few studies have examined the function of circRNAs in chemical contamination-induced diseases. CircRNAs are covalently closed continuous loops that do not possess 5' and 3' ends, increasing their structural stability and limiting degradation by exoribonucleases. In addition, environmental chemical exposure-related diseases are often accompanied by aberrant expression of specific circRNAs and those circRNAs are often detected in tissues and body fluids. Based on these characteristics, circRNAs may serve as candidate biomarkers for the diagnosis of diseases related to environmental chemical exposure. Here, we review the generation and function of circRNAs, and the possible molecular mechanisms underlying the regulation of environmental chemical exposure-related disorders by circRNAs. This is the first comprehensive review of the relationship between environmental chemical exposure and circRNAs in chemical exposure-induced diseases.

Cereal Circular RNAs (circRNAs): An Overview of the Computational Resources for Identification and Analysis

Circular RNAs (circRNAs) are a widespread class of endogenous noncoding RNAs and they have been studied in the past few years, implying important biological functions in all kingdoms of life. Recently, circRNAs have been identified in many plant species, including cereal crops, showing differential expression during stress response and developmental programs, which suggests their role in these process. In the following years, it is expected that insights into the functional roles of circRNAs can be used by cereal scientists and molecular breeders with the aim to develop new strategies for crop improvement. Here, we briefly outline the current knowledge about circRNAs in plants and we also outline available computational resources for their validation and analysis in cereal species.

The Uroboros Theory of Life's Origin: 22-Nucleotide Theoretical Minimal RNA Rings Reflect Evolution of Genetic Code and tRNA-rRNA Translation Machineries

Theoretical minimal RNA rings attempt to mimick life's primitive RNAs. At most 25 22-nucleotide-long RNA rings code once for each biotic amino acid, a start and a stop codon and form a stem-loop hairpin, resembling consensus tRNAs. We calculated, for each RNA ring's 22 potential splicing positions, similarities of predicted secondary structures with tRNA vs. rRNA secondary structures. Assuming rRNAs partly derived from tRNA accretions, we predict positive associations between relative secondary structure similarities with rRNAs over tRNAs and genetic code integration orders of RNA ring anticodon cognate amino acids. Analyses consider for each secondary structure all nucleotide triplets as potential anticodon. Anticodons for ancient, chemically inert cognate amino acids are most frequent in the 25 RNA rings. For RNA rings with primordial cognate amino acids according to tRNA-homology-derived anticodons, tRNA-homology and coding sequences coincide, these are separate for predicted cognate amino acids that presumably integrated late the genetic code. RNA ring secondary structure similarity with rRNA over tRNA secondary structures associates best with genetic code integration orders of anticodon cognate amino acids when assuming split anticodons (one and two nucleotides at the spliced RNA ring 5' and 3' extremities, respectively), and at predicted anticodon location in the spliced RNA ring's midst. Results confirm RNA ring homologies with tRNAs and CDs, ancestral status of tRNA half genes split at anticodons, the tRNA-rRNA axis of RNA evolution, and that single theoretical minimal RNA rings potentially produce near-complete proto-tRNA sets. Hence genetic code pre-existence determines 25 short circular gene- and tRNA-like RNAs. Accounting for each potential splicing position, each RNA ring potentially translates most amino acids, realistically mimicks evolution of the tRNA-rRNA translation machinery. These RNA rings 'of creation' remind the uroboros' (snake biting its tail) symbolism for creative regeneration.

Transcriptome-Wide Identification and Characterization of Circular RNAs in Leaves of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) in Response to Calcium Deficiency-Induced Tip-burn

Circular RNA (circRNA) is a newly discovered non-coding RNA, which play significant roles in the function and transcriptional regulation of microRNA. To date, in Chinese cabbage, the functional characteristic of circRNAs in response to calcium deficiency-induced tip-burn have not been reported. In this study, 730 circRNAs were isolated from Chinese cabbage leaves, of which 23 and 22 were differentially expressed in different calcium deficiency stages compared with the control. Forty-six host genes of the differentially expressed circRNAs were identified, and one circRNA was found to act as miRNAs sponges. Based on the functional analysis of host genes and target mRNAs of the corresponding miRNAs, the identified circRNAs might participated in response to stimulus, electron carrier activity, ATPase activity, cell wall metabolism, transcription factors and plant hormone signal transduction. ABF2, a positive regulator of the abiotic stress response in the abscisic acid (ABA) pathway, may play a role in calcium deficiency tolerance through a circRNA regulatory pathway. Correspondingly, the concentration of ABA is also increased during the Ca²⁺ deficiency stress. Our results suggest that circRNAs participate in a broad range of biological processes and physiological functions in the response to calcium deficiency-induced tip-burn and provide a basis for further studies of the biological roles that circRNAs play in the plant stress response.

Circular and Micro RNAs from Arabidopsis thaliana Flowers Are Simultaneously Isolated from AGO-IP Libraries

Competing endogenous RNAs (ceRNAs) are natural transcripts that can act as endogenous sponges of microRNAs (miRNAs), modulating miRNA action upon target mRNAs. Circular RNAs (circRNAs) are one among the various classes of ceRNAs. They are produced from a process called back-splicing and have been identified in many eukaryotes. In plants, their effective action as a miRNA sponge was not yet demonstrated. To address this question, public mRNAseq data from Argonaute-immunoprecipitation libraries (AGO-IP) of Arabidopsis thaliana flowers were used in association with a bioinformatics comparative multi-method to identify putative circular RNAs. A total of 27,812 circRNAs, with at least two reads at the back-splicing junction, were identified. Further analyses were used to select those circRNAs with potential miRNAs binding sites. As AGO forms a ternary complex with miRNA and target mRNA, targets count in AGO-IP and input libraries were compared, demonstrating that mRNA targets of these miRNAs are enriched in AGO-IP libraries. Through this work, five circRNAs that may function as miRNA sponges were identified and one of them were validated by PCR and sequencing. Our findings indicate that this post-transcriptional regulation can also occur in plants.

NAMS: Noncoding Assessment of long RNAs in Magnoliophyta Species

Regulation of plant transcriptional machinery has been recently demonstrated to be widely regulated by a class of long noncoding RNAs (lncRNAs) with size larger than 200 nt. The lncRNAs have been demonstrated to play key roles in abiotic stress. Taking into account the rapid pace in the development of the sequencing technologies, accelerated identification of lncRNAs with potential involvement in regulating the gene expression has been witnessed. Although progress has been made to identify the lncRNAs, however, accurate classification of the lncRNAs particularly in the case of plants is still challenging. In this protocol chapter, we present NAMS, which provides large-scale noncoding assessment of the lncRNAs specifically designed for Magnoliophyta species. We describe the approach and the usage of NAMS with potential applications for the lncRNA discovery.

Identification and characterization of CircRNAs involved in the regulation of wheat root length

BACKGROUND

Recent studies indicate that circular RNAs (circRNAs) may play important roles in the regulation of plant growth and development. Plant roots are the main organs of nutrient and water uptake. However, whether circRNAs involved in the regulation of plant root growth remains to be elucidated.

METHODS

LH9, XN979 and YN29 are three Chinese wheat varieties with contrasting root lengths. Here, the root circRNA expression profiles of LH9, XN979 and YN29 were examined by using high-throughput sequencing technology.

RESULTS

Thirty-three and twenty-two differentially expressed circRNAs (DECs) were identified in the YN29-LH9 comparison and YN29-XN979 comparison, respectively. Among them, ten DECs coexisted in both comparisons. As the roots of both LH9 and XN979 were significantly larger and deeper than YN29, the ten DECs coexisting in the two comparisons were highly likely to be involved in the regulation of wheat root length. Moreover, three of the ten DECs have potential miRNA binding sites. Real-time PCR analysis showed that the expression levels of the potential binding miRNAs exhibited significant differences between the long root plants and the short root plants.

CONCLUSIONS

The expression levels of some circRNAs exhibited significant differences in wheat varieties with contrasting root phenotypes. Ten DECs involved in the regulation of wheat root length were successfully identified in which three of them have potential miRNAs binding sites. The expression levels of putative circRNA-binding miRNAs were correlated with their corresponding circRNAs. Our results provide new clues for studying the potential roles of circRNAs in the regulation of wheat root length.

Diff isomiRs: Large-scale detection of differential isomiRs for understanding non-coding regulated stress omics in plants

Plants have an amazing ability to cope with wide variety of stresses by regulating the expression of genes and thus by altering the physiological status. In the past few years, canonical microRNA variants (isomiRs) have been shown to play pivotal roles by acting as regulators of the transcriptional machinery. In the present research, we present Diff isomiRs, a web-based exploratory repository of differential isomiRs across 16 sequenced plant species representing a total of 433 datasets across 21 different stresses and 158 experimental states. Diff isomiRs provides the high-throughput detection of differential isomiRs using mapping-based and model-based differential analysis revealing a total of 16,157 and 2,028 differential isomiRs, respectively. Easy-to-use and web-based exploration of differential isomiRs provides several features such as browsing of the differential isomiRs according to stress or species, as well as association of the differential isomiRs to targets and plant endogenous target mimics (PeTMs). Diff isomiRs also provides the relationship between the canonical miRNAs, isomiRs and the miRNA-target interactions. This is the first web-based large-scale repository for browsing differential isomiRs and will facilitate better understanding of the regulatory role of the isomiRs with respect to the canonical microRNAs. Diff isomiRs can be accessed at: www.mcr.org.in/diffisomirs.

Identification and characterization of circRNAs involved in the regulation of low nitrogen-promoted root growth in hexaploid wheat

BACKGROUND

CircRNAs are widespread in plants and play important roles in response to abiotic stresses. Low nitrogen (LN) promotes the growth of plant root system, allowing it to explore more nitrogen. However, whether circRNAs involved in the response to LN stress and the regulation of LN-promoted root growth in wheat remains unclear.

METHODS

Two wheat varieties (LH9 and XN979) with contrasting root phenotypes to LN stress were used as materials to identify circRNAs under control and LN conditions by using high-throughput sequencing technology.

RESULTS

Six differentially expressed circRNAs (DECs) involved in the common response to LN stress and 23 DECs involved in the regulation of LN-promoted root growth were successfully identified. GO analysis of the DEC-host genes involved in the regulation of LN-promoted root growth showed that GO terms related to biological regulation, responses to stimuli and signalling were significantly enriched. Moreover, seven DECs were predicted to have miRNA binding sites and may serve as miRNA sponges to capture miRNAs from their target genes.

CONCLUSIONS

LN stress altered the expression profiles of circRNAs in wheat. This is the first report of LN stress responsive circRNAs in plants. Our results provided new clues for investigating the functions of circRNAs in response to LN stress and in the regulation of LN-promoted wheat root growth.

Detection and profiling of circular RNAs in uninfected and maize Iranian mosaic virus-infected maize

Circular RNAs (circRNAs) are covalently closed non-coding RNAs that are usually derived from exonic regions of genes, but can also arise from intronic and intergenic regions. Studies of circRNAs in humans, animals and several plant species have shown an altered population of circRNAs in response to abiotic and biotic stress. Recently it was shown that circRNAs also occur in maize, but it is unknown if maize circRNAs are responsive to stress. Maize Iranian mosaic virus (MIMV, genus Nucleorhabdovirus, family Rhabdoviridae) causes an economically important disease in maize and other gramineous crops in Iran. In this study, we used data from RNA-Seq of MIMV-infected maize and uninfected controls to identify differentially expressed circRNAs. Such circRNAs were confirmed by two-dimensional polyacrylamide gel electrophoresis, northern blot, RT-qPCR and sequencing. A total of 1443 circRNAs were identified in MIMV-infected maize and 1165 circRNAs in uninfected maize. Two hundred and one circRNAs were in common between MIMV-infected and uninfected samples. Of these, 155 circRNAs were up-regulated and 5 down-regulated in MIMV infected plants, compared to the uninfected control. This study for the first time identified and profiled circRNA expression in maize in response to virus infection. Moreover, we predict that 33 circRNAs may bind 23 maize miRNAs, possibly affecting plant metabolism and development. Our data suggest a role for circRNAs in plant cell regulation and response to biotic stress such as virus infection, and give new insights into the complexity of plant-microbe interactions.

Genome-wide Identification and characterization of circular RNAs in the rice blast fungus Magnaporthe oryzae

Numerous circRNAs have been identified in different organisms, but little attention has been addressed on fungal circRNAs. Here, we identified a total of 8,848 circRNAs from the model plant pathogenic fungus M. oryzae. 5,840 circRNAs were identified from mycelium, 2,721 circRNAs from conidium, while only 287 circRNAs from both tissues. This indicated that most of the M. oryzae circRNAs were specifically expressed in mycelium or in conidium. Parental genes of circRNAs in mycelium were enriched in basic metabolisms required for normal growth, while in conidium, they were enriched in biogenesis of storages potentially used for infection. M. oryzae circRNAs could also bind to miRNAs, suggesting they may also function as sponges in fungi. This study suggested M. oryzae circRNAs could play important roles in regulation of growth and development.

Genome-wide identification of circRNAs involved in tomato fruit coloration

Circular RNAs (circRNAs) are an emerging class of non-coding RNAs in plants. Our aim is to identify the circRNAs with different expression levels between tomato fruits from two cultivars ('Jinling Fenyu' and 'Jinling Moyu') at two different developmental stages. Following high-throughput sequencing, 3796 circRNAs were identified, and 243 circRNAs were shared in the four samples. As compared with the fruit at mature green stage, the expression levels of 273 and 89 circRNAs were significantly altered in the fruit at turning stage from 'Jinling Fenyu' and 'Jinling Moyu', respectively. Moreover, the parental genes of the circRNAs with significantly different expression level were mainly involved in metabolic, cellular and single-organism process and played roles in catalytic activity and binding based on GO (Gene Ontology) analysis. The results suggested that circRNAs were widespread in tomato and were generated from different chromosomes and diverse genomic regions. Some circRNAs were specific in tomato fruits at different developmental stages, which enriches the number of circRNAs in plants involved in fruit coloration and ripening. This study provides the first genome-wide profile of circRNAs involved in tomato fruit coloration and lays a foundation for studying the potential biological functions of circRNAs involved in fruit ripening.

Genome-wide identification and characterization of circular RNAs by high throughput sequencing in soybean

Circular RNAs (circRNAs) arise during pre-mRNA splicing, in which the 3' and 5' ends are linked to each other by a covalent bond. Soybean is an ancient tetraploid, which underwent two whole genome duplications. Most of soybean genes are paralogous genes with multiple copies. Although many circRNAs have been identified in animals and plants, little is known about soybean circRNAs, especially about circRNAs derived from paralogous genes. Here, we used deep sequencing technology coupled with RNase R enrichment strategy and bioinformatic approach to uncover circRNAs in soybean. A total of 5,372 circRNAs were identified, approximately 80% of which were paralogous circRNAs generated from paralogous genes. Despite high sequence homology, the paralogous genes could produce different paralogous circRNAs with different expression patterns. Two thousand and one hundred thirty four circRNAs were predicted to be 92 miRNAs target mimicry. CircRNAs and circRNA isoforms exhibited tissue-specific expression patterns in soybean. Based on the function of circRNA-host genes, the soybean circRNAs may participate in many biological processes such as developmental process, multi-organism process, and metabolic process. Our study not only provided a basis for research into the function of circRNAs in soybean but also new insights into the plant circRNA kingdom.

isomiR2Function: An Integrated Workflow for Identifying MicroRNA Variants in Plants

In plants, post transcriptional regulation by non-coding RNAs (ncRNAs), in particular miRNAs (19-24 nt) has been involved in modulating the transcriptional landscape in developmental, biotic and abiotic interactions. In past few years, considerable focus has been leveraged on delineating and deciphering the role of miRNAs and their canonical isomiRs in plants. However, proper classification and accurate prediction of plant isomiRs taking into account the relative features by which we define isomiRs, such as templated or non-templated is still lacking. In the present research, we present isomiR2Function, a standalone easily deployable tool that allows for the robust and high-throughput discovery of templated and non-templated isomiRs. Additionally, isomiR2Function allows for identification of differentially expressed isomiRs and in parallel target prediction based on both transcripts or PARE-Seq either using Targetfinder or Cleaveland. isomiR2Function allows for the functional enrichment of the detected targets using TopGO package. Benchmarking of isomiR2Function revealed highly accurate prediction and classification of isomiRs as compared to the previously developed isomiR prediction tools. Additionally, the downstream implementation of additional features allows isomiR2Function to be classified as a single standalone tool for isomiR profiling from discovery to functional roles. All in all, isomiR2Function allows the streamline processing of the miRNA-seq for the identification and characterization of isomiRs with minimal efforts. isomiR2Function can be accessed through: https://github.com/347033139/isomiR2Function.

Genome-Wide Identification of Circular RNAs in Arabidopsis thaliana

Circular RNAs (circRNAs) are a family of transcripts with covalently closed circular structures and still largely unknown functions. Large numbers of circRNAs have been found in various biological processes in humans and animals, but fewer circRNAs have been identified in plants. We performed a genome-wide analysis of circRNAs in Arabidopsis thaliana via deep sequencing. We constructed 14 strand-specific libraries from 13 samples of plants from four developmental stages, four stress treatments, and five organs and a mixed sample across the lifespan. In total, we identified 5861 circRNAs, including 1275 novel ones, using the strict threshold of at least two unique back-spliced supporting reads. The circRNAs were non-randomly distributed in all chromosomes; most were exonic. Sequence similarity analysis of circRNAs between A. thaliana and four other species showed that some circRNAs are conserved in plants. Functional annotation indicated that many parent genes of circRNAs are involved in many fundamental processes including plant development, reproduction, and response to stimulus. In addition, a small proportion of circRNAs was shown to be potential targets of miRNAs, indicating that the circRNAs could interact with miRNAs to regulate gene expression. qRT-PCR analysis revealed that circRNAs displayed diverse expression patterns at different growth stages. Our results provide an important resource for continuing circRNA research in A. thaliana, and should enhance our understanding of circRNAs in plants.

Identification of Circular RNAs in Kiwifruit and Their Species-Specific Response to Bacterial Canker Pathogen Invasion

Research studies have recently focused on circle RNAs (circRNAs) in relation to their regulatory functions in animals. However, the systematic identification of circRNAs in plants, especially non-model plants, is limited. In addition, raw report on the prediction of the potential role of circRNAs in plant response to pathogen invasion is currently available. We conducted the systematic identification of circRNAs from four materials originating from three species belonging to genus Actinidia under different situations using ribosomal RNA (rRNA) depleted RNA-Seq data. A total of 3,582 circRNAs were identified in Actinidia, of which 64.01, 21.44, and 14.55% were intergenic circRNAs, exonic circRNAs, and intronic circRNAs, respectively. Tissue-specific expression of circRNAs was observed in kiwifruit, and a species-specific response was detected when infected with Pseudomonas syringae pv. actinidiae (Psa), which is the causative agent of kiwifruit bacterial canker disease. Furthermore, we found that both exonic and intronic circRNAs were significantly positively correlated to parent protein-coding genes, and intronic circRNAs are a class of highly remarkable regulators the parent genes comparing to that of exonic circRNAs. Expression and weighted gene co-expression network analysis (WGCNA) identified a set of circRNAs that were closely associated with plant defense response. The findings of the presents study suggest that circRNAs exhibit tissue- and species-specific expression, as well as play an important role in plant immune response.

Circular RNAs are miRNA sponges and can be used as a new class of biomarker

Circular RNAs (circRNAs) are a class of non-coding RNAs (ncRNAs) that are involved in transcriptional and posttranscriptional gene expression regulation. The development of deep sequencing of ribosomal RNA (rRNA)-depleted RNA libraries, associated with improved computational tools, has provided the identification of several new circRNAs in all sorts of organisms, from protists, plants and fungi to animals. Recently, it was discovered that endogenous circRNAs can work as microRNA (miRNA) sponges. This means that the circRNAs bind to miRNAs and consequently repress their function, providing a new model of action for this class of ncRNA, as well as indicating another mechanism that regulates miRNA activity. As miRNAs control a large set of biological processes, circRNA sponge activity will also affect these pathways. Several studies have associated miRNA sponges with human diseases, including osteoarthritis, diabetes, neurodegenerative pathologies and several types of cancer. Additionally, high stability, abundance and tissue-specific expression patterns make circRNA sponges very attractive for clinical research. Herein, we review the biogenesis, properties and function of endogenous circRNA sponges, with a special focus on those related to human cancer. A list of web tools available for the study of circRNAs is also given. Additionally, we discuss the possibility of using circRNAs as molecular markers for the diagnosis of diseases.

Identification of Circular RNAs and Their Targets in Leaves of Triticum aestivum L. under Dehydration Stress

Circular RNAs (circRNAs) are a type of newly identified non-coding RNAs through high-throughput deep sequencing, which play important roles in miRNA function and transcriptional controlling in human, animals, and plants. To date, there is no report in wheat seedlings regarding the circRNAs identification and roles in the dehydration stress response. In present study, the total RNA was extracted from leaves of wheat seedlings under dehydration-stressed and well-watered conditions, respectively. Then, the circRNAs enriched library based deep sequencing was performed and the circRNAs were identified using bioinformatics tools. Around 88 circRNAs candidates were isolated in wheat seedlings leaves while 62 were differentially expressed in dehydration-stressed seedlings compared to well-watered control. Among the dehydration responsive circRNAs, six were found to act as 26 corresponding miRNAs sponges in wheat. Sixteen circRNAs including the 6 miRNAs sponges and other 10 randomly selected ones were further validated to be circular by real-time PCR assay, and 14 displayed consistent regulation patterns with the transcriptome sequencing results. After Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the targeted mRNAs functions, the circRNAs were predicted to be involved in dehydration responsive process, such as photosynthesis, porphyrin, and chlorophyll metabolism, oxidative phosphorylation, amino acid biosynthesis, and metabolism, as well as plant hormone signal transduction, involving auxin, brassinosteroid, and salicylic acid. Herein, we revealed a possible connection between the regulations of circRNAs with the expressions of functional genes in wheat leaves associated with dehydration resistance.