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.

Transcriptomes Analysis Reveals the Regulatory Roles of Noncoding RNA in Tanshinones Synthesis Pathway of Salvia Miltiorrhiza

Identification of the Genome-Wide Expression Patterns of Non-Coding RNAs Associated with Tanshinones Synthesis Pathway in Salvia miltiorrhiza

Transcriptome Analysis of Medicinal Plant Salvia Miltiorrhiza and Identification of Genes Related to Tanshinone Biosynthesis

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Deep Sequencing Identifies Tissue-Specific Micrornas and Their Target Genes Involving in the Biosynthesis of Tanshinones in Salvia Miltiorrhiza

Comparative transcriptome analyses of different Salvia miltiorrhiza varieties during the accumulation of tanshinones

Spatiotemporal and Transcriptional Characterization on Tanshinone Initial Synthesis in <i>Salvia miltiorrhiza</i> Roots

Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues of Salvia miltiorrhiza and tanshinone biosynthesis

[Screening of long non-coding RNA related to CYP450s involved in biosynthesis of tanshinones]

Characterization of NAC family genes in<i> Salvia</i><i> miltiorrhiza</i> and NAC2 potentially involved in the biosynthesis of tanshinones

Integrated Transcriptomics and Proteomics to Reveal Regulation Mechanism and Evolution of SmWRKY61 on Tanshinone Biosynthesis in Salvia miltiorrhiza and Salvia castanea

Metabolic profiles and cDNA-AFLP analysis of Salvia miltiorrhiza and Salvia castanea Diel f. tomentosa Stib

Comprehensive transcriptome profiling of Salvia miltiorrhiza for discovery of genes associated with the biosynthesis of tanshinones and phenolic acids

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis

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

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

RNA interference targeting CYP76AH1 in hairy roots of Salvia miltiorrhiza reveals its key role in the biosynthetic pathway of tanshinones

Coding single nucleotide polymorphisms and SmCPS1 and SmKSL1 subcellular localization are associated with tanshinone biosynthesis in Salvia miltiorrhiza Bunge roots

Transcriptome sequencing of Salvia miltiorrhiza after infection by its endophytic fungi and identification of genes related to tanshinone biosynthesis

Functional Divergence of Diterpene Syntheses in the Medicinal Plant Salvia Miltiorrhiza

Transcriptional Data Mining of Salvia Miltiorrhiza in Response to Methyl Jasmonate to Examine the Mechanism of Bioactive Compound Biosynthesis and Regulation.