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

Descripción completa

Guardado en:
Detalles Bibliográficos
Autores principales: Mei Jiang, Haimei Chen, Qing Du, Liqiang Wang, Xinyue Liu, Chang Liu
Formato: article
Lenguaje:EN
Publicado: Frontiers Media S.A. 2021
Materias:
Acceso en línea:https://doaj.org/article/ce083f0ba15e482a983ad08f5e73b97d
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
id oai:doaj.org-article:ce083f0ba15e482a983ad08f5e73b97d
record_format dspace
spelling oai:doaj.org-article:ce083f0ba15e482a983ad08f5e73b97d2021-11-05T07:36:13ZGenome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza1664-802110.3389/fgene.2021.645115https://doaj.org/article/ce083f0ba15e482a983ad08f5e73b97d2021-11-01T00:00:00Zhttps://www.frontiersin.org/articles/10.3389/fgene.2021.645115/fullhttps://doaj.org/toc/1664-8021Circular 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.Mei JiangMei JiangMei JiangHaimei ChenQing DuQing DuLiqiang WangXinyue LiuChang LiuFrontiers Media S.A.articleSalvia miltiorrhizacircular RNAsSmAACT1SmDXS2SmHDR1SmKSL2GeneticsQH426-470ENFrontiers in Genetics, Vol 12 (2021)
institution DOAJ
collection DOAJ
language EN
topic Salvia miltiorrhiza
circular RNAs
SmAACT1
SmDXS2
SmHDR1
SmKSL2
Genetics
QH426-470
spellingShingle Salvia miltiorrhiza
circular RNAs
SmAACT1
SmDXS2
SmHDR1
SmKSL2
Genetics
QH426-470
Mei Jiang
Mei Jiang
Mei Jiang
Haimei Chen
Qing Du
Qing Du
Liqiang Wang
Xinyue Liu
Chang Liu
Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
description 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.
format article
author Mei Jiang
Mei Jiang
Mei Jiang
Haimei Chen
Qing Du
Qing Du
Liqiang Wang
Xinyue Liu
Chang Liu
author_facet Mei Jiang
Mei Jiang
Mei Jiang
Haimei Chen
Qing Du
Qing Du
Liqiang Wang
Xinyue Liu
Chang Liu
author_sort Mei Jiang
title Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
title_short Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
title_full Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
title_fullStr Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
title_full_unstemmed Genome-Wide Identification of Circular RNAs Potentially Involved in the Biosynthesis of Secondary Metabolites in Salvia miltiorrhiza
title_sort genome-wide identification of circular rnas potentially involved in the biosynthesis of secondary metabolites in salvia miltiorrhiza
publisher Frontiers Media S.A.
publishDate 2021
url https://doaj.org/article/ce083f0ba15e482a983ad08f5e73b97d
work_keys_str_mv AT meijiang genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT meijiang genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT meijiang genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT haimeichen genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT qingdu genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT qingdu genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT liqiangwang genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT xinyueliu genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
AT changliu genomewideidentificationofcircularrnaspotentiallyinvolvedinthebiosynthesisofsecondarymetabolitesinsalviamiltiorrhiza
_version_ 1718444498925649920