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

Abstract

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

Acknowledgements

This work was supported by the CAMS Innovation Fund for Medical Sciences (CIFMS) (2021-I2M-2-001).

Author contributions

S.L. and J.L. conceived and designed the experiments. J.L., C.L., and Y.D. performed the experiments. J.L. and H.W analysed the data. J.L. and S.L. wrote the manuscript. All authors read and agreed to the final article.

Data availability

All sequencing data generated are available at the NCBI SRA database under accession number PRJNA975272.

Conflict of interest statement

The authors declare no conflicts of interest.

References

[1] Chan SW, Henderson IR, Jacobsen SE . Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005; 6: 351-60.

[2] Law JA, Jacobsen SE . Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010; 113: 204-20.

[3] Zhang H, Lang Z, Zhu JK . Dynamics and function of DNA methylation in plants. Nat Rev Mol Cell Biol. 2018; 19: 489-506.

[4] Arora H, Singh RK, Sharma S et al. DNA methylation dynamics in response to abiotic and pathogen stress in plants. Plant Cell Rep. 2022; 41: 1931-44.

[5] Hou Q, Zhang T, Qi Y et al. Epigenetic dynamics and regulation of plant male reproduction. Int J Mol Sci. 2022; 23: 10420.

[6] Lucibelli F, Valoroso MC, Aceto S . Plant DNA methylation: an epigenetic mark in development, environmental interactions, and evolution. Int J Mol Sci. 2022; 23: 8299.

[7] Nguyen V, Gutzat R . Epigenetic regulation in the shoot apical meristem. Curr Opin Plant Biol. 2022; 69: 102267.

[8] Gehring M, Henikoff S . DNA methylation dynamics in plant genomes. Biochim Biophys Acta. 2007; 1769: 276-86.

[9] To TK, Kakutani T . Crosstalk among pathways to generate DNA methylome. Curr Opin Plant Biol. 2022; 68: 102248.

[10] Cuerda-Gil D, Slotkin RK . Non-canonical RNA-directed DNA methylation. Nat Plants. 2016; 2: 16163.

[11] Wendte JM, Zhang Y, Ji L et al. Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation. elife. 2019; 8: e47891.

[12] Papareddy RK, Paldi K, Smolka AD et al. Repression of CHROMOMETHYLASE 3 prevents epigenetic collateral damage in Arabidopsis. elife. 2021; 10: e69396.

[13] Stroud H, Do T, Du J et al. Non-CG methylation patterns shape the epigenetic landscape in Arabidopsis. Nat Struct Mol Biol. 2014; 21: 64-72.

[14] Zhu JK . Active DNA demethylation mediated by DNA glycosylases. Annu Rev Genet. 2009; 43: 143-66.

[15] Cokus SJ, Feng S, Zhang X et al. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning . Nature. 2008; 452: 215-9.

[16] Li X, Zhu J, Hu F et al. Single-base resolution maps of cultivated and wild rice methylomes and regulatory roles of DNA methylation in plant gene expression. BMC Genomics. 2012; 13: 300.

[17] Gent JI, Ellis NA, Guo L et al. CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize. Genome Res. 2013; 23: 628-37.

[18] Song QX, Lu X, Li QT et al. Genome-wide analysis of DNA methylation in soybean. Mol Plant. 2013; 6: 1961-74.

[19] Gardiner LJ, Quinton-Tulloch M, Olohan L et al. A genome-wide survey of DNA methylation in hexaploid wheat. Genome Biol. 2015; 16: 273.

[20] Xu W, Yang T, Dong X et al. Genomic DNA methylation analyses reveal the distinct profiles in castor bean seeds with persistent endosperms. Plant Physiol. 2016; 171: 1242-58.

[21] Zakrzewski F, Schmidt M, Van Lijsebettens M et al. DNA methylation of retrotransposons, DNA transposons and genes in sugar beet (Beta vulgaris L.). Plant J. 2017; 90: 1156-75.

[22] Zemach A, McDaniel IE, Silva P et al. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010; 328: 916-9.

[23] Zemach A, Kim MY, Silva P et al. Local DNA hypomethylation activates genes in rice endosperm. Proc Natl Acad Sci U S A. 2010; 107: 18729-34.

[24] Niederhuth CE, Bewick AJ, Ji L et al. Widespread natural variation of DNA methylation within angiosperms. Genome Biol. 2016; 17: 194.

[25] Xing MQ, Zhang YJ, Zhou SR et al. Global analysis reveals the crucial roles of DNA methylation during rice seed development. Plant Physiol. 2015; 168: 1417-32.

[26] Candaele J, Demuynck K, Mosoti D et al. Differential methylation during maize leaf growth targets developmentally regulated genes. Plant Physiol. 2014; 164: 1350-64.

[27] Niu LZ, Xu W, Ma PF et al. Single-base methylome analysis reveals dynamic changes of genome-wide DNA methylation associated with rapid stem growth of woody bamboos. Planta. 2022; 256: 53.

[28] Groszmann M, Greaves IK, Albert N et al. Epigenetics in plants-vernalisation and hybrid vigour. Biochim Biophys Acta. 2011; 1809: 427-37.

[29] Zhong S, Fei Z, Chen YR et al. Single-base resolution methylomes of tomato fruit development reveal epigenome modifications associated with ripening. Nat Biotechnol. 2013; 31: 154-9.

[30] Liu R, How-Kit A, Stammitti L et al. A DEMETER-like DNA demethylase governs tomato fruit ripening. Proc Natl Acad Sci U S A. 2015; 112: 10804-9.

[31] Yuan X, Ma K, Zhang M et al. Integration of transcriptome and methylome analyses provides insight into the pathway of floral scent biosynthesis in Prunus mume. Front Genet. 2021; 12: 779557.

[32] Han Y, Lu M, Yue S et al. Comparative methylomics and chromatin accessibility analysis in Osmanthus fragrans uncovers regulation of genic transcription and mechanisms of key floral scent production . Hortic Res. 2022; 9: uhac096.

[33] Liu HN, Shu Q, Lin-Wang K et al. DNA methylation reprogramming provides insights into light-induced anthocyanin biosynthesis in red pear. Plant Sci. 2022; 326: 111499.

[34] Lu S. Biosynthesis and regulatory mechanisms of bioactive compounds in Salvia miltiorrhiza, a model system for medicinal plant biology . Crit Rev Plant Sci. 2021; 40: 243-83.

[35] Zhang G, Tian Y, Zhang J et al. Hybrid de novo genome assembly of the Chinese herbal plant danshen (Salvia miltiorrhiza Bunge). Gigascience. 2015; 4: 62.

[36] Xu H, Song J, Luo H et al. Analysis of the genome sequence of the medicinal plant Salvia miltiorrhiza. Mol Plant. 2016; 9: 949-52.

[37] Song Z, Lin C, Xing P et al. A high-quality reference genome sequence of Salvia miltiorrhiza provides insights into tanshinone synthesis in its red rhizomes . Plant Genome. 2020; 13: e20041.

[38] Ma Y, Cui G, Chen T et al. Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis in Salvia miltiorrhiza. Nat Commun. 2021; 12: 685.

[39] Li S, Zhu N, Tang C et al. Differential distribution of characteristic constituents in root, stem and leaf tissues of Salvia miltiorrhiza using MALDI mass spectrometry imaging . Fitoterapia. 2020; 146: 104679.

[40] Ma Y, Yuan L, Wu B et al. Genome-wide identification and characterization of novel genes involved in terpenoid biosynthesis in Salvia miltiorrhiza. J Exp Bot. 2012; 63: 2809-23.

[41] Guo J, Zhou YJ, Hillwig ML et al. CYP76AH1 catalyzes turnover of miltiradiene in tanshinones biosynthesis and enables heterologous production of ferruginol in yeasts. Proc Natl Acad Sci U S A. 2013; 110: 12108-13.

[42] Guo J, Ma X, Cai Y et al. Cytochrome P450 promiscuity leads to a bifurcating biosynthetic pathway for tanshinones. New Phytol. 2016; 210: 525-34.

[43] Song JJ, Fang X, Li CY et al. A 2-oxoglutarate-dependent dioxygenase converts dihydrofuran to furan in Salvia diterpenoids . Plant Physiol. 2022; 188: 1496-506.

[44] Li J, Li C, Lu S . Identification and characterization of the cytosine-5 DNA methyltransferase gene family in Salvia miltiorrhiza. PeerJ. 2018; 6: e4461.

[45] Li J, Li C, Lu S . Systematic analysis of DEMETER-like DNA glycosylase genes shows lineage-specific Smi-miR7972 involved in SmDML1 regulation in Salvia miltiorrhiza. Sci Rep. 2018; 8: 7143.

[46] Yang BC, Lee MS, Lin MK et al. 5-Azacytidine increases tanshinone production in Salvia miltiorrhiza hairy roots through epigenetic modulation . Sci Rep. 2022; 12: 9349.

[47] Liang Z, Liu W, Feng X et al. Biological characters of danshen. In: Yan X, ed. Danshen (Salvia Miltiorrhiza) in Medicine. Dordrecht: Springer, 2015, 19-47.

[48] Wang H, Beyene G, Zhai J et al. CG gene body DNA methylation changes and evolution of duplicated genes in cassava. Proc Natl Acad Sci U S A. 2015; 112: 13729-34.

[49] Zhang H, Jin W, Zhu X et al. Identification and characterization of Salvia miltiorrhiza miRNAs in response to replanting disease . PLoS One. 2016; 11: e0159905.

[50] Wang L, Ma R, Liu C et al. Salvia miltiorrhiza: a potential red light to the development of cardiovascular diseases . Curr Pharm Des. 2017; 23: 1077-97.

[51] Meim XD, Cao YF, Che YY et al. Danshen: a phytochemical and pharmacological overview. Chin J Nat Med. 2019; 17: 59-80.

[52] Weinhofer I, Hehenberger E, Roszak P et al. H3K27me3 profiling of the endosperm implies exclusion of polycomb group protein targeting by DNA methylation. PLoS Genet. 2010; 6: e1001152.

[53] Lu S. Recent advances in biosynthesis and regulation of bioactive compounds in medicinal plants. Curr Pharm Biotechnol. 2021; 22: 720-1.

[54] Wang L, Shi Y, Chang X et al. DNA methylome analysis provides evidence that the expansion of the tea genome is linked to TE bursts. Plant Biotechnol J. 2019; 17: 826-35.

[55] Huang H, Liu R, Niu Q et al. Global increase in DNA methylation during orange fruit development and ripening. Proc Natl Acad Sci U S A. 2019; 116: 1430-6.

[56] Shao F, Lu S . Genome-wide identification, molecular cloning, expression profiling and posttranscriptional regulation analysis of the Argonaute gene family in salvia miltiorrhiza, an emerging model medicinal plant . BMC Genomics. 2013; 14: 512.

[57] Shao F, Lu S . Identification, molecular cloning and expression analysis of five RNA-dependent RNA polymerase genes in Salvia miltiorrhiza. PLoS One. 2014; 9: e95117.

[58] Shao F, Qiu D, Lu S . Comparative analysis of the Dicer-like gene family reveals loss of miR162 target site in SmDCL1 from Salvia miltiorrhiza. Sci Rep. 2015; 5: 9891.

[59] Zhao P, Ma B, Cai C et al. Transcriptome and methylome changes in two contrasting mungbean genotypes in response to drought stress. BMC Genomics. 2022; 23: 80.

[60] Ng DW, Miller M, Yu HH et al. A role for CHH methylation in the parent-of-origin effect on altered circadian rhythms and biomass heterosis in Arabidopsis intraspecific hybrids . Plant Cell. 2014; 26: 2430-40.

[61] Lu X, Wang X, Chen X et al. Single-base resolution methylomes of upland cotton (Gossypium hirsutum L.) reveal epigenome modifications in response to drought stress . BMC Genomics. 2017; 18: 297.

[62] Chiou CY, Pan HA, Chuang YN et al. Differential expression of carotenoid-related genes determines diversified carotenoid coloration in floral tissues of Oncidium cultivars . Planta. 2010; 232: 937-48.

[63] Shi M, Luo X, Ju G et al. Enhanced diterpene tanshinone accumulation and bioactivity of transgenic Salvia miltiorrhiza hairy roots by pathway engineering . J Agric Food Chem. 2016; 64: 2523-30.

[64] Shi M, Luo X, Ju G et al. Increased accumulation of the cardio-cerebrovascular disease treatment drug tanshinone in Salvia miltiorrhiza hairy roots by the enzymes 3-hydroxy-3-methylglutaryl CoA reductase and 1-deoxy-D-xylulose 5-phosphate reductoisomerase . Funct Integr Genomics. 2014; 14: 603-15.

[65] Dai Z, Cui G, Zhou SF et al. Cloning and characterization of a novel 3-hydroxy-3-methylglutaryl coenzyme a reductase gene from Salvia miltiorrhiza involved in diterpenoid tanshinone accumulation . J Plant Physiol. 2011; 168: 148-57.

[66] Kai G, Xu H, Zhou C et al. Metabolic engineering tanshinone biosynthetic pathway in Salvia miltiorrhiza hairy root cultures . Metab Eng. 2011; 13: 319-27.

[67] Zhang X, Guan H, Dai Z et al. Functional analysis of the isopentenyl diphosphate isomerase of Salvia miltiorrhiza via color complementation and RNA interference . Molecules. 2015; 20: 20206-18.

[68] Cheng Q, Su P, Hu Y et al. RNA interference-mediated repression of SmCPS (copalyldiphosphate synthase) expression in hairy roots of Salvia miltiorrhiza causes a decrease of tanshinones and sheds light on the functional role of SmCPS. Biotechnol Lett. 2014; 36: 363-9.

[69] Cui G, Duan L, Jin B et al. Functional divergence of diterpene syntheses in the medicinal plant Salvia miltiorrhiza. Plant Physiol. 2015; 169: 1607-18.

[70] Hu Z, Ren L, Bu J et al. Functional characterization of a 2OGD involved in abietane-type diterpenoids biosynthetic pathway in Salvia miltiorrhiza. Front Plant Sci. 2022; 13: 947674.

[71] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[72] Krueger F, Andrews SR . Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011; 27: 1571-2.

[73] Pertea M, Kim D, Pertea GM et al. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016; 11: 1650-67.

[74] Langmead B . Aligning short sequencing reads with bowtie. Curr Protoc Bioinformatics. 2010; 32: Unit 11.7.