DNA methylation is involved in sexual differentiation and sex chromosome evolution in the dioecious plant garden asparagus

Abstract

DNA methylation is a crucial regulatory mechanism in many biological processes. However, limited studies have dissected the contribution of DNA methylation to sexual differentiation in dioecious plants. In this study, we investigated the variances in methylation and transcriptional patterns of male and female flowers of garden asparagus. Compared with male flowers, female flowers at the same stages showed higher levels of DNA methylation. Both male and female flowers gained DNA methylation globally from the premeiotic to meiotic stages. Detailed analysis revealed that the increased DNA methylation was largely due to increased CHH methylation. Correlation analysis of differentially expressed genes and differentially methylated regions suggested that DNA methylation might not have contributed to the expression variation of the sex-determining genes SOFF and TDF1 but probably played important roles in sexual differentiation and flower development of garden asparagus. The upregulated genes AoMS1, AoLAP3, AoAMS, and AoLAP5 with varied methylated CHH regions might have been involved in sexual differentiation and flower development of garden asparagus. Plant hormone signaling genes and transcription factor genes also participated in sexual differentiation and flower development with potential epigenetic regulation. In addition, the CG and CHG methylation levels in the Y chromosome were notably higher than those in the X chromosome, implying that DNA methylation might have been involved in Y chromosome evolution. These data provide insights into the epigenetic modification of sexual differentiation and flower development and improve our understanding of sex chromosome evolution in garden asparagus.

References

[1] Panda, K. et al. Full-length autonomous transposable elements are preferentially targeted by expression-dependent forms of RNA-directed DNA methylation. Genome Biol. 17, 170 (2016).

[2] Meng, H. et al. DNA methylation, its mediators and genome integrity. Int. J. Biol. Sci. 11, 604-617 (2015).

[3] Li, X. et al. High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between DNA methylation, histone methylation, and gene expression. Plant Cell 20, 259-276 (2008).

[4] Wang, H. et al. CG gene body DNA methylation changes and evolution of duplicated genes in cassava. Proc. Natl Acad. Sci. USA 112, 13729-13734 (2015).

[5] Zemach, A., McDaniel, I. E., Silva, P. & Zilberman, D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science 328, 916-919 (2010).

[6] Niederhuth, C. E., Bewick, A. J., Ji, L., Alabady, M. S. & Schmitz, R. J. Widespread natural variation of DNA methylation within angiosperms. Genome Biol. 17, 194 (2016).

[7] Cheng, J. et al. Downregulation of RdDM during strawberry fruit ripening. Genome Biol. 19, 212 (2018).

[8] An, Y. C. et al. Dynamic changes of genome-wide DNA methylation during soybean seed development. Sci. Rep. 7, 12263 (2017).

[9] Huang, H. et al. Global increase in DNA methylation during orange fruit development and ripening. Proc. Natl Acad. Sci. USA 116, 1430-1436 (2019).

[10] Yaish, M. W., Al-Lawati, A., Al-Harrasi, I. & Patankar, H. V. Genome-wide DNA methylation analysis in response to salinity in the model plant caliph medic (Medicage truncatula). BMC Genomics 19, 78 (2018).

[11] Renner, S. S. The relative and absolute frequencies of angiosperm sexual systems: dioecy, monocy, gynodioecy, and un updated online database. Am. J. Bot. 101, 1588-1596 (2014).

[12] Bachtrog, D. et al. Sex determination: why so many ways of doing it? PLoS One 12, e1001899 (2014).

[13] Ming, R., Wang, J., Moore, P. H. & Paterson, A. H. Sex chromosomes in flowering plants. Am. J. Bot. 94, 141-150 (2007).

[14] Li, S. F., Zhang, G. J., Yuan, J. H., Deng, C. L. & Gao, W. J. Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes. Planta 243, 1083-1095 (2016).

[15] Lorenzo, J. L. R., Hobza, R. & Vyskot, B. DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia. BMC Genomics 19, 540 (2018).

[16] Chawla, A. et al. Sex-biased temporal gene expression in male and female floral buds of seabuckthorn (Hippophae rhamnoides). PLoS One 10, e0124890 (2015).

[17] Diggle, P. K. et al. Multiple developmental processes underlie sex differentiation in angiosperms. Trends Genet. 27, 368-376 (2011).

[18] Boualem, A. et al. A cucurbit androecy gene reveals how unisexual flowers develop and dioecy emerges. Science 350, 688-691 (2015).

[19] Akagi, T., Henry, I., Tao, R. & Comai, L. A Y-chromosome-encoded small RNA acts as a sex determinant in persimmons. Science 346, 646-650 (2014).

[20] Harkess, A. et al. The asparagus genome sheds light on the origin and evolution of a young Y chromosome. Nat. Commun. 8, 1279 (2017).

[21] Torres, M. F. et al. Genus-wide sequencing supports a two-locus model for sex-determination in Phoenix. Nat. Commun. 9, 3969 (2018).

[22] Akagi, T. et al. A Y-encoded suppressor of feminization arose via lineage-specific duplication of a cytokinin response regulator in kiwifruit. Plant Cell 30, 780-795 (2018).

[23] Urasaki, N. et al. Digital transcriptome analysis of putative sex-determination genes in papaya (Carica papaya). PLoS One 7, e40904 (2012).

[24] Harkess, A. et al. Sex-biased gene expression in dioecious garden asparagus (Asparagus officinalis). N. Phytol. 207, 883-892 (2015).

[25] Li, S. F. et al. Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis). BMC Plant Biol. 17, 143 (2017).

[26] Arumuganathan, K. & Earle, E. D. Nuclear DNA content of some important plant species. Plant Mol. Biol. Rep. 9, 208-218 (1991).

[27] Štajner, N., Bohanec, B. & Javornik, B. Genetic variability of economically important Asparagus species as revealed by genome size analysis and rDNA ITS polymorphisms. Plant Sci. 162, 931-937 (2002).

[28] Harkess, A. et al. Sex determination by two Y-linked genes in garden asparagus. Plant Cell 32, 1790-1796 (2020).

[29] Caporali, E., Carboni, A. & Galli, M. Development of male and female flowers in Asparagus officinalis. Search for point of transition from hermaphroditic to unisexual developmental pathway. Sex. Plant Reprod. 7, 239-249 (1994).

[30] Xi, Y. & Li, W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinforma. 10, 232 (2009).

[31] Lister, R. et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 462, 315 (2009).

[32] Akalin, A. et al. MethyKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 13, R87 (2012).

[33] Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[34] Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[35] Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 12, 1 (2011).

[36] Zhang, J. et al. Autotetraploid rice methylome analysis reveals methylation variation of transposable elements and their effects on gene expression. Proc. Natl Acad. Sci. USA 112, E7022-E7029 (2015).

[37] Zhou, P. et al. DNA methylome and transcriptome landscapes revealed differential characteristics of dioecious flowers in papaya. Hortic. Res. 7, 81 (2020).

[38] Rajkumar, M. S., Shankar, R., Garg, R. & Jain, M. Bisulphite sequencing reveals dynamic DNA methylation under desiccation and salinity stresses in rice cultivars. Genomics 112, 3537-3548 (2020).

[39] Shen, Y. et al. DNA methylation footprints during soybean domestication and improvement. Genome Biol. 19, 128 (2018).

[40] Janoršek, B., Široký, J. & Vyskot, B. Epigenetic control of sexual phenotype in a dioecious plant, Melandrium album. Mol. Gen. Genet. 250, 483-490 (1996).

[41] Akagi, T., Henry, I. M., Kawai, T., Comai, L. & Tao, R. Epigenetic regulation of the sex determination gene MeGI in polyploid persimmon. Plant Cell 28, 2905-2915 (2016).

[42] Abbott, J. K., Nordén, A. K. & Hansson, B. Sex chromosome evolution: historical insights and future perspectives. Proc. R. Soc. B 284, 20162806 (2017).

[43] Gorelick, R. Evolution of dioecy and sex chromosomes via methylation driving Muller′ ratchet. Biol. J. Linn. Soc. 80, 353-368 (2003).

[44] Metzger, D. C. H. & Schulte, P. M. The DNA methylation landscape of stickleback reveals patterns of sex chromosome evolution and effects of environmental salinity. Genome Biol. Evol. 10, 775-785 (2018).

[45] Charlesworth, D. Plant contributions to our understanding of sex chromosome evolution. N. Phytol. 208, 52-65 (2015).

[46] de Bello Cioffi, M. et al. The key role of repeated DNAs in sex chromosome evolution in two fish species with ZW sex chromosome system. Mol. Cytogenet. 5, 28 (2012).

[47] Noshay, J. M. et al. Monitoring the interplay between transposable element families and DNA methylation in maize. PLoS Genet. 15, e1008291 (2019).

[48] Zhou, W., Liang, G., Molloy, P. L. & Jones, P. A. DNA methylation enables transposable element-driven genome expansion. Proc. Natl Acad. Sci. USA 117, 19359-19366 (2020).

[49] Mirouze, M. et al. Loss of DNA methylation affects the recombination landscape in Arabidopsis. Proc. Natl Acad. Sci. USA 109, 5880-5885 (2012).

[50] Yelina, N. E. et al. DNA methylation epigenetically silences crossover hot spots and controls chromosomal domains of meiotic recombination in Arabidopsis. Genes Dev. 29, 2183-2202 (2015).

[51] Zhang, W., Wang, X., Yu, Q., Ming, R. & Jiang, J. DNA methylation and heterochromatinization in the male-specific region of the primitive Y chromosome of papaya. Genome Res. 18, 1938-1943 (2008).

[52] Li, H. et al. Time-series transcriptome provides insights into the gene regulation network involved in the volatile terpenoid metabolism during the flower development of lavender. BMC Plant Biol. 19, 313 (2019).

[53] Wellmer, F. & Riechmann, J. L. Gene networks controlling the initiation of flower development. Trends Genet. 26, 519-527 (2010).

[54] Li, N. et al. Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L. BMC Genomics 21, 850 (2020).

[55] Wilson, Z. A., Morroll, S. M., Dawson, J., Swarup, R. & Tighe, P. J. The Arabidopsis MALE STERILITY1 (MS1) gene is a transcriptional regulator of male gametogenesis, with homology to the PHD-finger family of transcription factors. Plant J. 28, 27-39 (2001).

[56] Xu, J. et al. The ABORTED MICROSPORES regulatory network is required for postmeiotic male reproductive development in Arabidopsis thaliana. Plant Cell 22, 91-107 (2010).

[57] Dobritsa, A. A. et al. LAP3, a novel plant protein required for pollen development, is essential for proper exine formation. Sex. Plant Reprod. 22, 167 (2009).

[58] Dobritsa, A. A. et al. LAP5 and LAP6 encode anther-specific proteins with similarity to chalcone synthase essential for pollen exine development in Arabidopsis. Plant Physiol. 153, 937-955 (2010).

[59] Yang, H. et al. Whole-genome DNA methylation patterns and complex associations with gene structure and expression during flower development in Arabidopsis. Plant J. 81, 268-281 (2015).

[60] Jacobsen, S. E., Sakai, H., Finnegan, E. J., Cao, X. & Meyerowitz, E. M. Ectopic hypermethylation of flower-specific genes in Arabidopsis. Curr. Biol. 10, 179-186 (2000).

[61] Gray, W. M. Hormonal regulation of plant growth and development. PLoS Biol. 2, e311 (2004).

[62] Bräutigam, K. et al. Sexual epigenetics: gender-specific methylation of a gene in the sex determining region of Populus balsamifera. Sci. Rep. 7, 45388 (2017).

[63] Xie, M. et al. A B-ARR-mediated cytokinin transcriptional network directs hormone cross-regulation and shoot development. Nat. Commun. 9, 1604 (2018).

[64] Han, Y., Yang, H. & Jiao, Y. Regulation of infloresence architecture by cytokinins. Front. Plant Sci. 5, 669 (2014).

[65] Reyes-Olaldem, J. I. et al. The bHLH transcription factor SPATULA enables cytokinin signaling, and both activate auxin biosysthesis and transport genes at the medial domain of the gynoecium. PLoS Genet. 13, e1006726 (2017).

[66] Nishino, K., Hattori, N., Tanaka, S. & Shiota, K. DNA methylation-mediated control of Sry gene expression in mouse gonadal development. J. Biol. Chem. 279, 22306-22313 (2004).

[67] Han, J. et al. High temperature induced masculinization of zebrafish by down-regulation of sox9b and esr1 via DNA methylation. J. Environ. Sci. 107, 160-170 (2021).

[68] Guo, S. et al. Co-ordination of flower development through epigenetic regulation in two model species: rice and Arabidopsis. Plant Cell Physiol. 56, 830-842 (2015).

[69] Pelayo, M. A., Yamaguchi, N. & Ito, T. One factor, many systems: the floral homeotic protein AGAMOUS and its epigenetic regulatory mechanisms. Curr. Opin. Plant Biol. 61, 102009 (2021).

[70] Martin, A. et al. A transposon-induced epigenetic change leads to sex determination in melon. Nature 461, 1135-1138 (2009).

[71] Sheldon, C. C. et al. The FLF MADS box gene: a repressor of flowering in Arabidopsis regulated by vernalization and methylation. Plant Cell 11, 445-458 (1999).