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

An Antisense Circular RNA Regulates Expression of RuBisCO Small Subunit Genes in Arabidopsis

Characterization and Function of Circular RNAs in Plants

Identification, Characterization, and Functional Prediction of Circular RNAs in Maize

Genome-wide discovery of circular rnas in the leaf and seedling tissues of Arabidopsis thaliana

Identifying and Characterizing the Circular RNAs during the Lifespan of Arabidopsis Leaves

A circRNA from SEPALLATA3 regulates splicing of its cognate mRNA through R-loop formation

Circrna Expression Pattern And Cerna And Mirna-Mrna Networks Involved In Anther Development In The Cms Line Of Brassica Campestris

Integrative analysis of Arabidopsis thaliana transcriptomics reveals intuitive splicing mechanism for circular RNA

Genome-Wide Identification of Dicer-Like, Argonaute, and Rna-Dependent Rna Polymerase Gene Families in Brassica Species and Functional Analyses of Their Arabidopsis Homologs in Resistance to Sclerotinia Sclerotiorum

Circular RNA ath‐circ032768, a competing endogenous RNA, response the drought stress by targeting miR472‐RPS5 module

The antisense CircRNA VvcircABH controls salt tolerance and the brassinosteroid signaling response by suppressing cognate mRNA splicing in grape

Small RNA-mediated chromatin silencing directed to the 3′ region of the Arabidopsis gene encoding the developmental regulator, FLC

A lariat-derived circular RNA is required for plant development in Arabidopsis

The RNA-binding protein RBP45D of Arabidopis plays a role in epigenetic control of flowering time and DCL3-independent RNA-directed DNA methylation

Camelina CircRNA Landscape: Implications for Gene Regulation and Fatty Acid Metabolism

Vesicular and non-vesicular extracellular small RNAs direct gene silencing in a plant-interacting bacterium

CMA33/XCT Regulates Small RNA Production Through Modulating the Transcription of Dicer-Like Genes in Arabidopsis

A 5′ Trna-Ala-derived Small RNA Regulates Anti-Fungal Defense in Plants

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

Antisense inhibition of Rubisco activase increases Rubisco content and alters the proportion of Rubisco activase in stroma and thylakoids in chloroplasts of rice leaves.

The untranslated mRNA leader modulates carotenoid feedback and post-transcriptional regulation