Methylation related genes affect sex differentiation in dioecious and gynodioecious papaya

Abstract

Morphological, genic and epigenetic differences often exist in separate sexes of dioecious and trioecious plants. However, the connections and relationships among them in different breeding systems are still unclear. Papaya has three sex types, which is genetically determined and epigenetically regulated, and was chosen as a model to study sex differentiation. Bisulfite sequencing of genomic DNA extracted from early-stage flowers revealed sex-specific genomic methylation landscapes and seasonally methylome reprogramming processes in dioecious and gynodioecious papaya grown in spring and summer. Extensive methylation of sex-determining region (SDR) was the distinguishing epigenetic characteristics of nascent XY sex chromosomes in papaya. Seasonal methylome reprogramming of early-stage flowers in both dioecy and gynodioecy systems were detected, resulting from transcriptional expression pattern alterations of methylation-modification-related and chromatin-remodeling-related genes, particularly from those genes involved in active demethylation. Genes involved in phytohormone signal transduction pathway in male flowers have played an important role in the formation of male-specific characteristics. These findings enhanced the understanding of the genetic and epigenetic contributions to sex differentiation and the complexity of sex chromosome evolution in trioecious plants.

References

[1] Wils CR, Kaufmann K . Gene-regulatory networks controlling inflorescence and flower development in Arabidopsis thaliana . Biochim Biophys Acta Gene Regul Mech. 2017; 1860: 95-105.

[2] Thomson B, Wellmer F . Molecular regulation of flower development. Curr Top Dev Biol. 2019; 131: 185-210.

[3] Ming R, Bendahmane A, Renner SS . Sex chromosomes in land plants. Annu Rev Plant Biol. 2011; 62: 485-514.

[4] Lebel-Hardenack S, Grant SR . Genetics of sex determination in flowering plants. Trends Plant Sci. 1997; 2: 130-6.

[5] Aryal R, Ming R . Sex determination in flowering plants: papaya as a model system. Plant Sci. 2014; 217-218: 56-62.

[6] Vyskot B, Hobza R . The genomics of plant sex chromosomes. Plant Sci. 2015; 236: 126-35.

[7] Lin H, Liao Z, Zhang L et al. Transcriptome analysis of the male-to-hermaphrodite sex reversal induced by low temperature in papaya. Tree Genet Genomes. 2016; 12: 1-14.

[8] Zhou P, Fatima M, Ma X et al. Auxin regulation involved in gynoecium morphogenesis of papaya flowers. Hortic Res. 2019; 6: 119-9.

[9] Wang J, Na JK, Yu Q et al. Sequencing papaya X and Yh chromosomes reveals molecular basis of incipient sex chromosome evolution. Proc Natl Acad Sci. 2012; 109: 13710.

[10] VanBuren R, Zeng F, Chen C et al. Origin and domestication of papaya Yh chromosome. Genome Biol. 2015; 25: 524-33.

[11] Ramos HCC, Pereira MG, Silva FF et al. Seasonal and genetic influences on sex expression in a backcrossed segregating papaya population. Crop Breed Applied Biotechnol. 2011; 11: 97-105.

[12] Damasceno Junior P, Pereira T, Pereira M . Estimation of genetic parameters for flower anomalies in papaya. Crop Breed Applied Biotechnol. 2018; 18: 9-15.

[13] Moreira N et al. Quantification of floral abnormalities in a population generated from sexual polymorphism aiming at recurrent selection in papaya. Bragantia. 2019; 78.

[14] Song Y, Ma K, Bo W et al. Sex-specific DNA methylation and gene expression in andromonoecious poplar. Plant Cell Rep. 2012; 31: 1393-405.

[15] Bräutigam K, Soolanayakanahally R, Champigny M et al. Sexual epigenetics: gender-specific methylation of a gene in the sex determining region of Populus balsamifera. Sci Rep. 2017; 7: 45388.

[16] Müller NA, Kersten B, Leite Montalvão AP et al. A single gene underlies the dynamic evolution of poplar sex determination. Nature Plants. 2020; 6: 630-7.

[17] Yang W, Wang D, Li Y et al. A general model to explain repeated turnovers of sex determination in the Salicaceae. Mol Biol Evol. 2020; 38: 968-80.

[18] Xue L, Wu H, Chen Y et al. Evidences for a role of two Y-specific genes in sex determination in Populus deltoides. Nat Commun. 2020; 11: 5893.

[19] Rodríguez Lorenzo JL, Hobza R, Vyskot B . DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia. BMC Genomics. 2018; 19: 540.

[20] Bačovský V, Houben A, Kumke K et al. The distribution of epigenetic histone marks differs between the X and Y chromosomes in Silene latifolia. Planta. 2019; 250: 487-494.

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

[22] Chen J, Hu Q, Zhang Y et al. P-MITE: a database for plant miniature inverted-repeat transposable elements. Nucleic Acids Res. 2014; 42: D1176-81.

[23] Li SF, Zhang GJ, Zhang XJ et al. DPTEdb, an integrative database of transposable elements in dioecious plants. Database. 2016; 2016: 1-10.

[24] Parrilla-Doblas JT, Roldán-Arjona T, Ariza RR et al. Active DNA demethylation in plants. Int J Mol Sci. 2019; 20: 4683.

[25] Elhamamsy AR . DNA methylation dynamics in plants and mammals: overview of regulation and dysregulation. Cell Biochem Funct. 2016; 34: 289-98.

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

[27] Kenchanmane Raju SK, Ritter EJ, Niederhuth CE et al. Establishment, maintenance, and biological roles of non-CG methylation in plants. Essays Biochem. 2019; 63: 743-55.

[28] Movahedi A, Sang M, Zhang J et al. Functional analyses of PtROS1-RNAi in poplars and evaluation of its effect on DNA methylation. J Plant Biol. 2018; 61: 227-40.

[29] Liu J, Wu X, Yao X et al. Mutations in the DNA demethylase OsROS1 result in a thickened aleurone and improved nutritional value in rice grains. Proc Natl Acad Sci U S A. 2018; 115: 11327-32.

[30] Lang Z, Lei M, Wang X et al. The methyl-CpG-binding protein MBD7 facilitates active DNA demethylation to limit DNA hypermethylation and transcriptional gene silencing. Mol Cell. 2015; 57: 971-83.

[31] Liang X, Hou X, Li J et al. High-resolution DNA methylome reveals that demethylation enhances adaptability to continuous cropping comprehensive stress in soybean. BMC Plant Biol. 2019; 19: 79-9.

[32] Song Y, Ma K, Ci D et al. Sexual dimorphic floral development in dioecious plants revealed by transcriptome, phytohormone, and DNA methylation analysis in Populus tomentosa. Plant Mol Biol. 2013;83:559-76.

[33] Nie WF, Lei M, Zhang M et al. Histone acetylation recruits the SWR1 complex to regulate active DNA demethylation in Arabidopsis. Proc Natl Acad Sci. 2019;116:16641.

[34] Hasegawa J, Sakamoto T, Fujimoto S et al. Auxin decreases chromatin accessibility through the TIR1/AFBs auxin signaling pathway in proliferative cells. Sci Rep. 2018; 8: 7773.

[35] Hobza R, Hudzieczek V, Kubat Z et al. Sex and the flower - developmental aspects of sex chromosome evolution. Ann Bot. 2018; 122: 1085-101.

[36] Zhang W, Wang X, Yu Q et al. DNA methylation and heterochromatinization in the male-specific region of the primitive Y chromosome of papaya. Genome Res. 2008; 18: 1938-43.

[37] Yu Q, Hou S, Hobza R et al. Chromosomal location and gene paucity of the male specific region on papaya Y chromosome. Mol Gen Genomics. 2007; 278: 177-85.

[38] Iovene M, Yu Q, Ming R et al. Evidence for emergence of sex-determining gene(s) in a Centromeric region in Vasconcellea parviflora. Genetics. 2015; 199: 413-21.

[39] Gremski K, Ditta G, Yanofsky MF . The HECATE genes regulate female reproductive tract development in Arabidopsis thaliana. Development. 2007; 134: 3593-601.

[40] Zuñiga-Mayo VM, Baños-Bayardo CR, Díaz-Ramírez D et al. Conserved and novel responses to cytokinin treatments during flower and fruit development in Brassica napus and Arabidopsis thaliana. Sci Rep. 2018; 8: 6836.

[41] Reeves P, Ellis CM, Ploense SE et al. A regulatory network for coordinated flower maturation. PLoS Genet. 2012; 8: e1002506.

[42] Liu K, Yuan C, Li H et al. Genome-wide identification and characterization of auxin response factor (ARF) family genes related to flower and fruit development in papaya (Carica papaya L.). BMC Genomics. 2015; 16: 901-1.

[43] Wang L, Wang C, Liu X et al. Peroxisomal β-oxidation regulates histone acetylation and DNA methylation in Arabidopsis . Proc Natl Acad Sci. 2019; 116: 10576.

[44] Li SF, Lv CC, Lan LN et al. DNA methylation is involved in sexual differentiation and sex chromosome evolution in the dioecious plant garden asparagus. Hortic Res. 2021; 8: 198.

[45] Gonsalves D. Transgenic papaya in Hawaii and beyond. AgBioforum. 2004; 7: 36-40.

[46] Gonsalves D . In: Ming R, Moore PH, eds. Hawaii’s Transgenic Papaya Story 1978-2012: A Personal Account in Genetics and Genomics of Papaya . Springer: New York, 2014, 115-42.

[47] Tuo D, Shen W, Yang Y et al. Development and validation of a multiplex reverse transcription PCR assay for simultaneous detection of three papaya viruses. Viruses. 2014; 6: 3893-906.

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

[49] Akalin A, Kormaksson M, Li S et al. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012; 13: R87.

[50] Dobin A, Davis CA, Schlesinger F et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2012; 29: 15-21.

[51] Trapnell C, Hendrickson DG, Sauvageau M et al. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013; 31: 46-53.

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