Comparative methylomics and chromatin accessibility analysis in Osmanthus fragrans uncovers regulation of genic transcription and mechanisms of key floral scent production

Abstract

Linalool and ionone are two important aromatic components in sweet osmanthus petals, and the regulatory mechanisms that produce these two components remain unclear. In this study, we employed whole-genome methylation sequencing and ATAC-seq technology to analyze the genomic DNA methylation status and chromatin accessibility of the sweet osmanthus cultivars ‘Zaohuang’ and ‘Chenghong Dangui’. Results showed that the promoter region of TPS2, a key gene in the linalool synthesis pathway, was less methylated in ‘Chenghong Dangui’ than in ‘Zaohuang’. The chromatin was more accessible in ‘Chenghong Dangui’ than in ‘Zaohuang’, which resulted in a much stronger expression of this gene in ‘Chenghong Dangui’ than in ‘Zaohuang’. This eventually led to a high quantity of linalool and its oxides in the petals of ‘Chenghong Dangui’, but there were lower levels present in the petals of ‘Zaohuang’. These results suggest that DNA methylation and chromatin accessibility play major roles in linalool synthesis in sweet osmanthus. The methylation level of the promoter region of CCD4, a key gene for ionone synthesis, was higher in ‘Zaohuang’ than in ‘Chenghong Dangui’. The chromatin accessibility was lower in ‘Zaohuang’ than in ‘Chenghong Dangui’, although the expression of this gene was significantly higher in ‘Zaohuang’ than in ‘Chenghong Dangui’. ChIP-seq analysis and a series of experiments showed that the differential expression of CCD4 and CCD1 in the two cultivars may predominantly be the result of regulation by ERF2 and other transcription factors. However, a 183-bp deletion involving the CCD4 promoter region in ‘Chenghong Dangui’ may be the main reason for the low expression of this gene in its petals. This study provides an important theoretical basis for improving selective breeding of key floral fragrance components in sweet osmanthus.

References

[1] Xiang QB, Liu YL . An illustrated monograph of the sweet Osmanthus variety in China . Zhejiang Science & Technology Press, 2007, 93-260.

[2] Zeng XL, Liu C, Zheng RR et al. Emission and accumulation of monoterpene and the key terpene synthase (TPS) associated with monoterpene biosynthesis in Osmanthus fragrans lour. Front Plant Sci. 2016; 6: 1232.

[3] Yang X, Yue Y, Li H et al. The chromosome-level quality genome provides insights into the evolution of the biosynthesis genes for aroma compounds of Osmanthus fragrans . Hortic Res. 2018; 5: 72.

[4] Han YJ, Wang HY, Wang XD et al. Mechanism of floral scent production in Osmanthus fragrans and the production and regulation of its key floral constituents, β-ionone and linalool. Hortic Res. 2019; 6: 106.

[5] Sun BJ, Li L, Han YJ et al. HS-SPME-GC-MS analysis of different Osmanthus fragrans cultivars from Guilin Garden in Shanghai . J Fujian Coll Forest. 2012; 32: 39-42.

[6] Frusciante S, Diretto G, Bruno M et al. Novel carotenoid cleavage dioxygenase catalyzes the first dedicated step in saffron crocin biosynthesis. PNAS Nexus. 2014; 111: 12246-51.

[7] Taylor IB, Sonneveld T, Bugg TDH et al. Regulation and manipulation of the biosynthesis of abscisic acid, including the supply of xanthophyll precursors. J Plant Growth Regul. 2005; 24: 253-73.

[8] Huang FC, Molnár P, Schwab W . Cloning and functional characterization of carotenoid cleavage dioxygenase 4 genes. J Exp Bot. 2009; 60: 3011-22.

[9] Baldermann S, Kato M, Kurosawa M et al. Functional characterization of a carotenoid cleavage dioxygenase 1 and its relation to the carotenoid accumulation and volatile emission during the floral development of Osmanthus fragrans lour. J Exp Bot. 2010; 61: 2967-77.

[10] Susanne B. Biosynthesis of α- and β-ionone, prominent scent compounds, in flowers of Osmanthus fragrans . Acta Biochim Pol. 2012; 59: 79-81.

[11] Zhang X, Pei J, Zhao L et al. Overexpression and characterization of CCD4 from Osmanthus fragrans and β-ionone biosynthesis from β-carotene in vitro . J Mol Catal B Enzym. 2016; 134: 105-14.

[12] Han YJ, Wu M, Cao LY et al. Characterization of OfWRKY3, a transcription factor that positively regulates the carotenoid cleavage dioxygenase gene OfCCD4 in Osmanthus fragrans . Plant Mol Biol. 2016; 91: 485-96.

[13] Han YJ, Liu LX, Dong MF et al. cDNA cloning of the phytoene synthase (PSY) and expression analysis of PSY and carotenoid cleavage dioxygenase genes in Osmanthus fragrans . Biologia. 2013; 68: 258-63.

[14] Han YJ, Wang XH, Chen WC et al. Differential expression of carotenoid-related genes determines diversified carotenoid coloration in flower petal of Osmanthus fragrans . Tree Genet Genomes. 2014; 10: 329-38.

[15] Tholl D, Chen F, Petri J et al. Two sesquiterpene synthases are responsible for the complex mixture of sesquiterpenes emitted from Arabidopsis flowers. Plant J. 2005; 42: 757-71.

[16] Garms S, Kollner TG, Boland W . A multiproduct terpene synthase from Medicago truncatula generates cadalane sesquiterpenes via two different mechanisms. J Org Chem. 2010; 75: 5590-600.

[17] Dudareva N, Martin D, Kish CM et al. (E)-β-ocimene and myrcene synthase genes of floral scent biosynthesis in snapdragon: function and expression of three terpene synthase genes of a new terpene synthase subfamily. Plant Cell. 2003; 15: 1227-41.

[18] Nagegowda DA, Gutensohn M, Wilkerson CG et al. Two nearly identical terpene synthases catalyze the formation of nerolidol and linalool in snapdragon flowers. Plant J. 2008; 55: 224-39.

[19] Nieuwenhuizen NJ, Green SA, Chen X et al. Functional genomics reveals that a compact terpene synthase gene family can account for terpene volatile production in apple. Plant Physiol. 2013; 161: 787-804.

[20] Green SA, Chen X, Nieuwenhuizen NJ et al. Identification, functional characterization, and regulation of the enzyme responsible for floral (E)-nerolidol biosynthesis in kiwifruit (Actinidia chinensis) . J Exp Bot. 2012; 63: 1951-67.

[21] Martin DM, Aubourg S, Schouwey M et al. Functional annotation, genome organization and phylogeny of the grapevine (Vitis vinifera) terpene synthase gene family based on genome assembly, FLcDNA cloning, and enzyme assays . BMC Plant Biol. 2010; 10: 226.

[22] Sos BC, Fung HL, Gao DR et al. Characterization of chromatin accessibility with a transposome hypersensitive sites sequencing (THS-seq) assay. Genome Biol. 2016; 17: 20.

[23] Buenrostro JD, Giresi PG, Zaba LC et al. Transposition of native chromatin for multimodal regulatory analysis and personal epigenomics. Nat Methods. 2013; 10: 1213-8.

[24] Liu GF, Wang W, Hu SG et al. Inherited DNA methylation primes the establishment of accessible chromatin during genome activation. Genome Res. 2018; 28: 998-1007.

[25] Sijacic P, Bajic M, McKinney EC et al. Changes in chromatin accessibility between Arabidopsis stem cells and mesophyll cells illuminate cell type-specific transcription factor networks. Plant J. 2018; 94: 215-31.

[26] Domcke S, Bardet AF, Ginno PA et al. Competition between DNA methylation and transcription factors determines binding of NRF1. Nature. 2015; 528: 575-9.

[27] Zhu H, Wang GH, Qian J . Transcription factors as readers and effectors of DNA methylation. Nat Rev Genet. 2016; 17: 551-65.

[28] Zhou P, Zhang XD, Fatima M et al. DNA methylome and transcriptome landscapes revealed differential characteristics of dioecious flowers in papaya. Horti Res. 2020; 7: 81.

[29] Lang ZB, Wang YH, Tang K et al. Critical roles of DNA demethylation in the activation of ripening-induced genes and inhibition of ripening-repressed genes in tomato fruit. PNAS Nexus. 2017; 114: E4511-9.

[30] Yang XL, Yue YZ, Li HY et al. The chromosome-level quality genome provides insights into the evolution of the biosynthesis genes for aroma compounds of Osmanthus fragrans . Hortic. Res. 2018; 5: 72.

[31] Lister R, O’Malley RC, Tonti-Filippini J et al. Highly integrated single-base resolution maps of the epigenome in Arabidopsis . Cell. 2008; 133: 523-36.

[32] Cheng JF, Niu QF, Zhang B et al. Downregulation of RdDM during strawberry fruit ripening. Genome Biol. 2018; 19: 212.

[33] Regulski M, Lu Z, Kendall J et al. The maize methylome influences mrna splice sites and reveals widespread paramutation-like switches guided by small rna. Genome Res. 2013; 23: 1651-62.

[34] Bouyer D, Kramdi A, Kassam M et al. DNA methylation dynamics during early plant life. Genome Biol. 2017; 18: 179.

[35] Lu ZF, Hofmeister BT, Vollmers C et al. Combining ATAC-seq with nuclei sorting for discovery of cis-regulatory regions in plant genomes. Nucleic Acids Res. 2016; 45: e41.

[36] Yu ZX, Li JX, Yang CQ et al. The jasmonate-responsive AP2/ERF transcription factors AaERF1 and AaERF2 positively regulate artemisinin biosynthesis in Artemisia annua L. Mol Plant. 2012; 5: 353-65.

[37] Min T, Fang F, Ge H et al. Two novel anoxia-induced ethylene response factors that interact with promoters of deastringency-related genes from persimmon. PLoS One. 2014; 9: e97043.

[38] Chen HG, Zeng XL, Yang J et al. Whole-genome resequencing of Osmanthus fragrans provides insights into flower color evolution. Hortic. Res. 2021; 8: 98.

[39] Ohme-Takagi M, Shinshi H . Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995; 7: 173-82.

[40] Fujimoto SY, Ohta M, Usui A et al. Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC-box-mediated gene expression . Plant Cell. 2000; 12: 393-404.

[41] Wang LM, Li MT, Jin WW et al. Variations in the components of Osmanthus fragrans lour. Essential oil at different stages of flowering. Food Chem. 2009; 114: 233-6.

[42] Xin H, Wu B, Zhang H et al. Characterization of volatile compounds in flowers from four groups of sweet osmanthus (Osmanthus fragrans) cultivars . Can J Plant Sci. 2013; 93: 923-31.

[43] Bentsen M, Goymann P, Schultheis H et al. ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation. Nat Commun. 2020; 11: 4267.

[44] Albanus RD, Kyono Y, Hensley J et al. Chromatin information content landscapes inform transcription factor and DNA interactions. Nat Commun. 2021; 12: 1307.

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

[46] Jullien PE, Susaki D, Yelagandula R et al. DNA methylation dynamics during sexual reproduction in Arabidopsis thaliana . Curr Biol. 2012; 22: 1825-30.

[47] Park K, Kim MY, Vickers M et al. DNA demethylation is initiated in the central cells of Arabidopsis and rice. PNAS Nexus. 2016; 113: 15138-43.

[48] Lei M, Zhang H, Julian R et al. Regulatory link between DNA methylation and active demethylation in Arabidopsis . Proc Natl Acad Sci USA. 2015; 112: 3553-7.

[49] Ligterink W, Kroj T, Zur Nieden U et al. Receptor-mediated activation of a MAP kinase in pathogen defense of plants. Science. 1997; 276: 2054-7.

[50] Mikolajczyk M, Awotunde OS, Muszynska G et al. Osmotic stress induces rapid activation of a salicylic acid-induced protein kinase and a homolog of protein kinase ASK1 in tobacco cells. Plant Cell. 2000; 12: 165-78.

[51] Ren DT, Liu YD, Yang KY et al. A fungal-responsive MAPK cascade regulates phytoalexin biosynthesis in Arabidopsis . Proc Natl Acad Sci USA. 2008; 105: 5638-43.

[52] Sangwan V, Orvar BL, Beyerly J et al. Opposite changes in membrane fluidity mimic cold and heat stress activation of distinct plant MAP kinase pathways. Plant J. 2002; 31: 629-38.

[53] Mao G, Meng X, Liu Y et al. Phosphorylation of a WRKY transcription factor by two pathogen-responsive MAPKs drives phytoalexin biosynthesis in Arabidopsis . Plant Cell. 2011; 23: 1639-53.

[54] Zhang M, Su J, Yan Z et al. Conveying endogenous and exogenous signals: MAPK cascades in plant growth and defense. Curr Opin Plant Biol. 2018; 45: 1-10.

[55] Shoji T, Hashimoto T . Tobacco MYC2 regulates jasmonate inducible nicotine biosynthesis genes directly and by way of the NIC2-locus ERF genes. Plant Cell Physiol. 2011; 52: 1117-30.

[56] Zhang HT, Hedhili S, Montiel G et al. The basic helix-loop-helix transcription factor CrMYC2 controls the jasmonate responsive expression of the ORCA genes that regulate alkaloid biosynthesis in Catharanthus roseus . Plant J. 2011; 67: 61-71.

[57] Endo T, Fujii H, Sugiyama A et al. Overexpression of a citrus basic helix-loop-helix transcription factor (CubHLH1), which is homologous to Arabidopsis activation-tagged bri1 suppressor 1 interacting factor genes, modulates carotenoid metabolism in transgenic tomato. Plant Sci. 2016; 243: 35-48.

[58] Hong GJ, Xue XY, Mao YB et al. Arabidopsis MYC2 interacts with DELLA proteins in regulating sesquiterpene synthase gene expression . Plant Cell. 2012; 24: 2635-48.

[59] Schwab W, Fischer T, Giri A et al. Potential applications of glucosyltransferases in terpene glucoside production: impacts on the use of aroma and fragrance. European journal of applied microbiology and biotechnology. 2015; 99: 165-74.

[60] Zheng RR, Zhu ZY, Wang YL et al. Ugt85a84 catalyzes the glycosylation of aromatic monoterpenes in Osmanthus fragrans lour. Flowers. Front Plant Sci. 2019; 10: 1376-6.

[61] Felix K, Simon A . Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011; 27: 1571-2.

[62] Yu G, Wang LG, Han YY et al. ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7.

[63] Bajic M, Maher KA, Deal RB . Identification of open chromatin regions in plant genomes using ATAC-Seq. Methods Mol Biol. 2018; 1675: 183-201.

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

[65] Love MI, Huber W, Anders S . Moderated estimation of fold-change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550.

[66] Kaufmann K, Muino JM, Osteras M et al. Chromatin immunoprecipitation (ChIP) of plant transcription factors followed by sequencing (ChIP-SEQ) or hybridization to whole genome arrays (ChIP-CHIP). Protocols. 2010; 5: 457-72.

[67] Shamimuzzaman M, Vodkin L . Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-Seq and RNA-Seq. BMC Genomics. 2013; 14: 477.

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