A chromosome-scale genome assembly and epigenomic profiling reveal temperature-dependent histone methylation in iridoid biosynthesis regulation in Scrophularia ningpoensis

Abstract

Understanding how medicinal plants adapt to global warming, particularly through epigenetic mechanisms that modify phenotypes without changing DNA sequences is crucial. Scrophularia ningpoensis Hemsl., a traditional Chinese Medicine (TCM), produces bioactive compounds that are influenced by environmental temperatures, making it an ideal model for studying the biological basis of TCM geoherbalism. However, the adaptive potential of epigenetic marks in S. ningpoensis under varying temperatures remains understudied, partly due to the absence of a reference genome. Here, it was demonstrated that mild warm temperatures contribute to the metabolic accumulation and the cultivated migration of S. ningpoensis using a global dataset. A high-quality chromosome-level genome was assembled, and an atlas of epigenetic, metabolic, and transcriptomic profiles across different tissues. Transcriptome analysis identified 3401 allele-specific expressed genes (ASEGs) across nine tissues by comparing two haplotypes. ChIP-seq and BS-seq data from leaf and root tissues revealed that ASEGs are associated with distinct epigenetic patterns, particularly the active mark H3K36me3, which functions differently in these tissues. Notably, genes marked with H3K36me3 in iridoid synthesis pathway predominantly expressed in roots. Additionally, the histone methyltransferase SnSDG8 was identified to regulate ectopic H3K36me3 in iridoid biosynthesis in response to warming temperatures. Our results highlight the epigenetic mechanisms of global warming on herb-derived products, significant for medicinal plant breeding under temperature stress.

Acknowledgements

The financial assistance for this research was provided by the grants from the Open Project of Hubei Key Laboratory of Wudang Local Chinese Medicine Research (WDCM2023014) to C.Z. and Principal Investigator Program (HBMUPI202104) to Y. Z. at Hubei University of Medicine, and Open Foundation of Hubei Province Key Laboratory of Tumor Microenvironment and Immunotherapy (2023KZL026) to Q.X.

Author Contributions

C.Z. and Y.Z. figured the original idea; Q.X., B.L. and K.T. analyzed the data; C.L., L.X., H.W. and L.Y. conducted experiments; W.Z. and M.Z. kept and offered the seed of S. ningpoensis; C.Z. wrote the paper with input from other authors.

Data availability

The raw sequence data included in this manuscript have been deposited in the Genome Sequence Archive at the BIG Data Center of the Beijing Institute of Genomics (BIG), Chinese Academy of Sciences. Access to this data is available to the public under the accession number PRJCA026437, which can be found at the following URL: https://bigd.big.ac.cn/gsa/.

Conflict of interest statement

The authors declare no conflicts of interest.

Supplementary data

Supplementary data is available at HORTRE Journal online.

References

[1] Huang L, Guo L, Ma C. et al. Top-geoherbs of traditional Chinese medicine: common traits, quality characteristics and forma-tion. Front Med. 2011;5: 185-94

[2] Li Y, Kong D, Fu Y. et al. The effect of developmental and environmental factors on secondary metabolites in medicinal plants. Plant Physiol Biochem. 2020;148: 80-9

[3] Zandalinas SI, Balfagon D, Gomez-Cadenas A. et al. Plant responses to climate change: metabolic changes under com-bined abiotic stresses. JExp Bot. 2022;73: 3339-54

[4] Sun Y, Alseekh S, Fernie AR. Plant secondary metabolic responses to global climate change: a meta-analysis in medic-inal and aromatic plants. Glob Chang Biol. 2023;29: 477-504

[5] Thuiller W. Biodiversity: climate change and the ecologist. Nature. 2007;448: 550-2

[6] Stocker TF. Climate change. The closing door of climate targets. Science. 2013;339: 280-2

[7] Zhao Q, Mi ZY, Lu C. et al. Predicting potential distribution of Ziziphus spinosa (Bunge) H.H. Hu ex F.H. Chen in China under climate change scenarios. Ecol Evol. 2022;12:e8629

[8] Yang Y, He J, Liu Y. et al. Assessment of Chinese suitable habitats of Zanthoxylum nitidum in different climatic conditions by Maxent model, HPLC, and chemometric methods. Ind Crop Prod. 2023;196:116515

[9] Zhao Z, Guo P, Brand E. The formation of daodi medicinal materials. J Ethnopharmacol. 2012;140: 476-81

[10] Kim NH, Jayakodi M, Lee SC. et al. Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol J. 2018;16: 1904-17

[11] Liu T, Yang J, Liu S. et al. Regulation of chlorogenic acid, flavonoid, and iridoid biosynthesis by histone H3K4 and H3K9 methylation in Lonicera japonica. Mol Biol Rep. 2020;47: 9301-11

[12] Charron JB, He H, Elling AA. et al. Dynamic landscapes of four histone modifications during deetiolation in Arabidopsis. Plant Cell. 2009;21: 3732-48

[13] Pikaard CS, Mittelsten Scheid O. Epigenetic regulation in plants. Cold Spring Harb Perspect Biol. 2014;6:a019315

[14] Malabarba J, Windels D, Xu W. et al. Regulation of DNA (de)methylation positively impacts seed germination during seed development under heat stress. Genes (Basel). 2021;12:457

[15] Li J, Huang Q, Sun M. et al. Global DNA methylation variations after short-term heat shock treatment in cultured microspores of Brassica napus cv. Topas. Sci Rep. 2016;6:38401

[16] Hossain MS, Kawakatsu T, Kim KD. et al. Divergent cytosine DNA methylation patterns in single-cell, soybean root hairs. New Phytol. 2017;214: 808-19

[17] Ma Y, Min L, Wang M. et al. Disrupted genome methyla-tion in response to high temperature has distinct affects on microspore abortion and anther indehiscence. Plant Cell. 2018;30: 1387-403

[18] Pecinka A, Dinh HQ, Baubec T. et al. Epigenetic regulation of repetitive elements is attenuated by prolonged heat stress in Arabidopsis. Plant Cell. 2010;22: 3118-29

[19] Correia B, Valledor L, Meijon M. et al. Is the interplay between epigenetic markers related to the acclimation of cork oak plants to high temperatures? PLoS One. 2013;8:e53543

[20] Hou H, Zhao L, Zheng X. et al. Dynamic changes in histone modification are associated with upregulation of Hsf and rRNA genes during heat stress in maize seedlings. Protoplasma. 2019;256: 1245-56

[21] Song ZT, Zhang LL, Han JJ. et al. Histone H3K4 methyltrans-ferases SDG25 and ATX1 maintain heat-stress gene expression during recovery in Arabidopsis. Plant J. 2021;105: 1326-38

[22] Ueda M, Matsui A, Nakamura T. et al. Versatility of HDA19-deficiency in increasing the tolerance of Arabidopsis to different environmental stresses. Plant Signal Behav. 2018; 13:1-4

[23] Pajoro A, Severing E, Angenent GC. et al. Histone H 3 lysine 36 methylation affects temperature-induced alternative splicing and flowering in plants. Genome Biol. 2017;18:102

[24] Zhao Z, Yu Y, Meyer D. et al. Prevention of early flowering by expression of FLOWERING LOCUS C requires methylation of histone H3 K36. Nat Cell Biol. 2005;7: 1256-60

[25] Adrienne BN, Deborah LS, Gemma LH. et al. Adaptive plasticity and epigenetic variation in response to warming in an Alpine plant. Ecol Evol. 2015;5: 634-47

[26] Yang S, Li J, Zhao Y. et al. Harpagoside variation is positively correlated with temperature in Scrophularia ningpoensis Hemsl. J Agric Food Chem. 2011;59: 1612-21

[27] Diaz AM, Abad MJ, Fernandez L. et al. Phenylpropanoid glyco-sides from Scrophularia scorodonia: in vitro anti-inflammatory activity. Life Sci. 2004;74: 2515-26

[28] Lee EJ, Kim SR, Kim J. et al. Hepatoprotective phenylpropanoids from Scrophularia buergeriana roots against CCl4-induced tox-icity: action mechanism and structure-activity relationship. Planta Med. 2002;68: 407-11

[29] Wu Q, Wen XD, Qi LW. et al. An in vivo microdialysis mea-surement of harpagoside in rat blood and bile for predicting hepatobiliary excretion and its interaction with cyclosporin a and verapamil. J Chromatogr B Analyt Technol Biomed Life Sci. 2009;877: 751-6

[30] Mariana C, Sílvia C, João L. Genome size variation and incidence of polyploidy in Scrophulariaceae sensu lato from the Iberian Peninsula. AoB Plants. 2012;2012:pls037

[31] Ai HL, Qin CL, Ye K. et al. Characterization of the complete chloroplast genome of a well-known Chinese medicinal herb, Scrophularia ningpoensis. Mitochondrial DNA B Resour. 2020;5: 484-5

[32] SunJ LiX,QuZ. et al. Comparative proteomic analysis reveals novel insights into the continuous cropping induced response in Scrophularia ningpoensis. J Sci Food Agric. 2023;103: 1832-45

[33] Garcia RA, Cabeza M, Rahbek C. et al. Multiple dimensions of climate change and their implications for biodiversity. Science. 2014;344:1247579

[34] Wu J, Li X, Huang L. et al. A new GIS model for ecologically suitable distributions of medicinal plants. Chin Med. 2019;14:4

[35] Scheunert A, Heubl G. Diversification of Scrophularia (Scro-phulariaceae) in the Western Mediterranean and Macaronesia-phylogenetic relationships, reticulate evolution and biogeo-graphic patterns. Mol Phylogenet Evol. 2014;70: 296-313

[36] Zhang X, Zhang S, Zhao Q. et al. Assembly of allele-aware, chromosomal-scale autopolyploid genomes based on hi-C data. Nat Plants. 2019;5: 833-45

[37] Wang Y, Tang H, Debarry JD. et al. MCScanX: a toolkit for detec-tion and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40:e49

[38] Goel M, Sun H, Jiao WB. et al. SyRI: finding genomic rearrange-ments and local sequence differences from whole-genome assemblies. Genome Biol. 2019;20:277

[39] Waters AJ, Makarevitch I, Noshay J. et al. Natural variation for gene expression responses to abiotic stress in maize. Plant J. 2017;89: 706-17

[40] Cooper RD, Shaffer HB. Allele-specific expression and gene regulation help explain transgressive thermal tolerance in non-native hybrids of the endangered California tiger salamander (Ambystoma californiense). Mol Ecol. 2021;30: 987-1004

[41] Zhang C, Yang Z, Tang D. et al. Genome design of hybrid potato. Cell. 2021;184: 3873-3883.e12

[42] Zhou C, Liu X, Li X. et al. A genome doubling event reshapes Rice morphology and products by modulating chromatin signatures and gene expression profiling. Rice (NY). 2021;14:72

[43] Tan F, Zhou C, Zhou Q. et al. Analysis of chromatin regulators reveals specific features of rice DNA methylation pathways. Plant Physiol. 2016;171: 2041-54

[44] Zhang M, Xu C, Yan H. et al. Limited tissue culture-induced mutations and linked epigenetic modifications in F hybrids of sorghum pure lines are accompanied by increased transcrip-tion of DNA methyltransferases and 5-methylcytosine glyco-sylases. Plant J. 2009;57: 666-79

[45] Feng S, Cokus SJ, Zhang X. et al. Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci. 2010;107: 8689-94

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

[47] Li J, Ahn JH, Wang GG. Understanding histone H3 lysine 36 methylation and its deregulation in disease. Cell Mol Life Sci. 2019;76: 2899-916

[48] Zhou C, Wang C, Liu H. et al. Identification and analysis of adenine N(6)-methylation sites in the rice genome. Nat Plants. 2018;4: 554-63

[49] Zhou C, Zhou H, Ma X. et al. Genome-wide identification and characterization of Main histone modifications in sorghum decipher regulatory mechanisms involved by mRNA and long noncoding RNA genes. J Agric Food Chem. 2021;69: 2337-47

[50] He G, Zhu X, Elling AA. et al. Global epigenetic and transcrip-tional trends among two rice subspecies and their reciprocal hybrids. Plant Cell. 2010;22: 17-33

[51] Du Z, Li H, Wei Q. et al. Genome-wide analysis of histone modifications: H3K4me2, H3K4me3, H3K9ac, and H3K27ac in Oryza sativa L. Japonica. Mol Plant. 2013;6: 1463-72

[52] Zhou Q, Tang D, Huang W. et al. Haplotype-resolved genome analyses of a heterozygous diploid potato. Nat Genet. 2020;52: 1018-23

[53] Ernst J, Kellis M. Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat Biotechnol. 2010;28: 817-25

[54] Lee JS, Smith E, Shilatifard A. The language of histone crosstalk. Cell. 2010;142: 682-5

[55] Sequeira-Mendes J, Gutierrez C. Links between genome repli-cation and chromatin landscapes. Plant J. 2015;83: 38-51

[56] Ren D, Shen ZY, Qin LP. et al. Pharmacology, phytochem-istry, and traditional uses of Scrophularia ningpoensis Hemsl. J Ethnopharmacol. 2021;269:113688

[57] Jensen SR, Franzyk H, Wallander E. Chemotaxonomy of the Oleaceae: iridoids as taxonomic markers. Phytochemistry. 2002;60: 213-31

[58] Vranova E, Coman D, Gruissem W. Structure and dynamics of the isoprenoid pathway network. Mol Plant. 2012;5: 318-33

[59] Pu X, Dong X, Li Q. et al. An update on the function and regula-tion of methylerythritol phosphate and mevalonate pathways and their evolutionary dynamics. J Integr Plant Biol. 2021;63: 1211-26

[60] Vranova E, Coman D, Gruissem W. Network analysis of the MVA and MEP pathways for isoprenoid synthesis. Annu Rev Plant Biol. 2013;64: 665-700

[61] Jin H, Song Z, Nikolau BJ. Reverse genetic characterization of two paralogous acetoacetyl CoA thiolase genes in Arabidopsis reveals their importance in plant growth and development. Plant J. 2012;70: 1015-32

[62] KanY, MuXR, GaoJ. et al. The molecular basis of heat stress responses in plants. Mol Plant. 2023;16: 1612-34

[63] Jagadish SVK, Way DA, Sharkey TD. Plant heat stress: concepts directing future research. Plant Cell Environ. 2021;44: 1992-2005

[64] Zhang X, Chen S, Shi L. et al. Haplotype-resolved genome assembly provides insights into evolutionary history of the tea plant Camellia sinensis. Nat Genet. 2021;53: 1250-9

[65] Gil N, Ulitsky I. Regulation of gene expression by cis-acting long non-coding RNAs. Nat Rev Genet. 2020;21: 102-17

[66] Samo N, Ebert A, Kopka J. et al. Plant chromatin, metabolism and development - an intricate crosstalk. Curr Opin Plant Biol. 2021;61:102002

[67] Lindermayr C, Rudolf EE, Durner J. et al. Interactions between metabolism and chromatin in plant models. Mol Metab. 2020;38:100951

[68] Wang X, Zhang L, Yao G. et al. De novo chromosome-level genome assembly of Chinese motherwort (Leonurus japonicus). Sci Data. 2024;11:55

[69] Marcais G, Kingsford C. A fast, lock-free approach for effi-cient parallel counting of occurrences of k-mers. Bioinformatics. 2011;27: 764-70

[70] Ranallo-Benavidez TR, Jaron KS, Schatz MC. GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat Commun. 2020;11:1432

[71] Wang C, Liu C, Roqueiro D. et al. Genome-wide analysis of local chromatin packing in Arabidopsis thaliana. Genome Res. 2015;25: 246-56

[72] Cheng H, Concepcion GT, Feng X. et al. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021;18: 170-5

[73] Durand NC, Shamim MS, Machol I. et al. Juicer provides a one-click system for analyzing loop-resolution hi-C experiments. Cell Syst. 2016;3: 95-8

[74] Durand NC, Robinson JT, Shamim MS. et al. Juicebox provides a visualization system for hi-C contact maps with unlimited zoom. Cell Syst. 2016;3: 99-101

[75] Manni M, Berkeley MR, Seppey M. et al. BUSCO: assessing genomic data quality and beyond. Curr Protoc. 2021;1:e323

[76] Chen N. Using RepeatMasker to identify repetitive ele-ments in genomic sequences. Curr Protoc Bioinformatics. 2004;5: 4.10.1-14

[77] Bao W, Kojima KK, Kohany O. Repbase update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015;6:11

[78] Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27: 573-80

[79] Xu Z, Wang H. LTR_FINDER: an efficient tool for the pre-diction of full-length LTR retrotransposons. Nucleic Acids Res. 2007;35: W265-8

[80] Flynn JM, Hubley R, Goubert C. et al. RepeatModeler2 for auto-mated genomic discovery of transposable element families. Proc Natl Acad Sci. 2020;117: 9451-7

[81] Griffiths-Jones S, Moxon S, Marshall M. et al. Rfam: annotat-ing non-coding RNAs in complete genomes. Nucleic Acids Res. 2005;33: D121-4

[82] Chan PP, Lin BY, Mak AJ. et al. tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes. Nucleic Acids Res. 2021;49: 9077-96

[83] Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12: 357-60

[84] Pertea M, Pertea GM, Antonescu CM. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33: 290-5

[85] Slater GS, Birney E. Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005; 6:31

[86] Stanke M, Keller O, Gunduz I. et al. AUGUSTUS: ab initio predic-tion of alternative transcripts. Nucleic Acids Res. 2006;34: W435-9

[87] Majoros WH, Pertea M, Salzberg SL. TigrScan and Glim-merHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004;20: 2878-9

[88] Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12:491

[89] Boeckmann B, Bairoch A, Apweiler R. et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003;31: 365-70

[90] Apweiler R, Attwood TK, Bairoch A. et al. The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res. 2001;29: 37-40

[91] Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28: 27-30

[92] Emms DM, Kelly S. OrthoFinder: phylogenetic orthology infer-ence for comparative genomics. Genome Biol. 2019;20:238

[93] Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usabil-ity. Mol Biol Evol. 2013;30: 772-80

[94] Nguyen LT, Schmidt HA, von Haeseler A. et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32: 268-74

[95] Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24: 1586-91

[96] Mendes FK, Vanderpool D, Fulton B. et al. CAFE 5 models vari-ation in evolutionary rates among gene families. Bioinformatics. 2021;36: 5516-8

[97] Zhang Z, Xiao J, Wu J. et al. ParaAT: a parallel tool for construct-ing multiple protein-coding DNA alignments. Biochem Biophys Res Commun. 2012;419: 779-81

[98] Wang D, Zhang Y, Zhang Z. et al. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding win-dow strategies. Genomics Proteomics Bioinformatics. 2010;8: 77-80

[99] Chen S, Zhou Y, Chen Y. et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34: i884-90

[100] Liao Y, Smyth GK, Shi W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019;47:e47

[101] Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559

[102] Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9: 357-9

[103] Li H, Handsaker B, Wysoker A. et al. The sequence alignmen-t/map format and SAMtools. Bioinformatics. 2009;25: 2078-9

[104] Ramirez F, Dundar F, Diehl S. et al. deepTools: a flexible plat-form for exploring deep-sequencing data. Nucleic Acids Res. 2014;42: W187-91

[105] Zhang Y, Liu T, Meyer CA. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137

[106] Yu G, Wang LG, He QY. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioin-formatics. 2015;31: 2382-3

[107] Robinson JT, Thorvaldsdottir H, Winckler W. et al. Integrative genomics viewer. Nat Biotechnol. 2011;29: 24-6

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

[109] Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Methods. 2012;9: 215-6

[110] Chong J, Yamamoto M, Xia J. MetaboAnalystR 2.0: from raw spectra to biological insights. Meta. 2019;9:57

[111] Zhou H, Yan R, He H. et al. Transcriptional profiling of long noncoding RNAs associated with flower color formation in Ipomoea nil. Planta. 2023;258:6