DNA methylome and transcriptome landscapes revealed differential characteristics of dioecious flowers in papaya

Abstract

Separate sexes in dioecious plants display different morphology and physiological characteristics. The differences between the two sexes lie in their highly differentiated floral characteristics and in sex-related phenotype, which is genetically determined and epigenetically modified. In dioecious papaya (Carica papaya L.), global comparisons of epigenetic DNA methylation and gene expressions were still limited. We conducted bisulfite sequencing of early-stage flowers grown in three seasons (spring, summer and winter) and compared their methylome and transcriptome profiles to investigate the differential characteristics of male and female in papaya. Methylation variances between female and male papaya were conserved among three different seasons. However, combined genome-scale transcriptomic evidence revealed that most methylation variances did not have influence on the expression profiles of neighboring genes, and the differentially expressed genes were most overrepresented in phytohormone signal transduction pathways. Further analyses showed diverse stress-responsive methylation alteration in male and female flowers. Male flower methylation was more responsive to stress whereas female flower methylation varied less under stress. Early flowering of male papaya in spring might be associated with the variation in the transcription of CpSVP and CpAP1 coinciding with their gene-specific hypomethylation. These findings provide insights into the sex-specific DNA methylation and gene expression landscapes of dioecious papaya and a foundation to investigate the correlation between differentiated floral characteristics and their candidate genes.

References

[1] Ming, R., Bendahmane, A. & Renner, S. S. Sex chromosomes in land plants. Annu. Rev. Plant Biol. 62, 485-514 (2011).

[2] Aryal, R. & Ming, R. Sex determination in flowering plants: papaya as a model system. Plant Sci. 217−218, 56-62 (2014).

[3] Juarez, C. & Banks, J. A. Sex determination in plants. Curr. Top. Dev. Biol. 49, 70-72 (1998).

[4] 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 (2016).

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

[6] Janousek, B., Siroký, J. & Vyskot, B. Epigenetic control of sexual phenotype in a dioecious plant, Melandrium album. Mol. Gen. Genet 250, 483-490 (1996).

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

[8] Liu, Z. et al. A primitive Y chromosome in papaya marks incipient sex chromosome evolution. Nature 427, 348-352 (2004).

[9] Vanburen, R. et al. Origin and domestication of papaya Yh chromosome. Genome Res. 25, 524-533 (2015).

[10] Yu, Q. et al. Recent origin of dioecious and gynodioecious Y chromosomes in papaya. Trop. Plant Biol. 1, 49-57 (2008).

[11] Lin, H., Liao, Z., Zhang, L. & Yu, Q. Transcriptome analysis of the male-to-hermaphrodite sex reversal induced by low temperature in papaya. Tree Genet. Genomes 12, 94 (2016).

[12] Hofmeyr, J. D. J. Sex reversal in carica Papaya L. South Afr. J. Sci. 36, 286-287 (1939).

[13] 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 (2008).

[14] Tuo, D. et al. Development and validation of a multiplex reverse transcription PCR assay for simultaneous detection of three papaya viruses. Viruses 6, 3893-3906 (2014).

[15] Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120 (2014).

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

[17] Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589-595 (2010).

[18] Garrison, E. & Marth, G. Haplotype-based variant detection from short-read sequencing. arXiv 1207, 3907 (2012).

[19] Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562-578 (2012).

[20] Wu, J., Mao, X., Cai, T., Luo, J. & Wei, L. KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 34, W720-W724 (2006).

[21] Xie, C. et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39, W316-W322 (2011).

[22] Chen, C., Xia, R., Chen, H. & He, Y . TBtools, a Toolkit for biologists integrating various HTS-data handling tools with a user-friendly interface. bioRxiv. https://doi.org/10.1101/289660 (2018).

[23] Elhamamsy, A. R. DNA methylation dynamics in plants and mammals: overview of regulation and dysregulation. Cell Biochem. Funct. 34, 289-298 (2016).

[24] Zhang, H., Lang, Z. & Zhu, J.-K. Dynamics and function of DNA methylation in plants. Nat. Rev. Mol. Cell Biol. 19, 489-506 (2018).

[25] Wang, J. et al. Sequencing papaya X and Yh chromosomes reveals molecular basis of incipient sex chromosome evolution. Proc. Natl. Acad. Sci. USA 109, 13710-13715 (2012).

[26] Zhou, P., Fatima, M., Ma, X., Liu, J. & Ming, R. Auxin regulation involved in gynoecium morphogenesis of papaya flowers. Hortic. Res. 6, 119 (2019).

[27] Eva, H., Nijuscha, G. & Alexander, H. The more, the merrier: cytokinin signaling beyond Arabidopsis. Plant Signal Behav. 5, 1384-1390 (2010).

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

[29] Nayelli, M. M. et al. The role of cytokinin during Arabidopsis gynoecia and fruit morphogenesis and patterning. Plant J. Cell Mol. Biol. 72, 222-234 (2012).

[30] Isabel, B., Elisabeth, O., Miroslav, S., Tomáš, W. & Thomas, S. Cytokinin regulates the activity of reproductive meristems, flower organ size, ovule formation, and thus seed yield in Arabidopsis thaliana. Plant Cell 23, 69-80 (2011).

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

[32] 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).

[33] Huang, X. et al. The antagonistic action of abscisic acid and cytokinin signaling mediates drought stress response in Arabidopsis. Mol. Plant 11, 970-982 (2018).

[34] Varas, J., Santos, J. L. & Pradillo, M. The absence of the Arabidopsis chaperone complex CAF-1 produces mitotic chromosome abnormalities and changes in the expression profiles of genes involved in DNA repair. Front. Plant Sci. 8, 525 (2017).

[35] Pavlištová, V. et al. Phenotypic reversion in fas mutants of Arabidopsis thaliana by reintroduction of FAS genes: variable recovery of telomeres with major spatial rearrangements and transcriptional reprogramming of 45S rDNA genes. Plant J. 88, 414-424 (2016).

[36] Veronika, M. et al. Homology-dependent repair is involved in 45 S rDNA loss in plant CAF-1 mutants. Plant J. Cell Mol. Biol. 81, 198-209 (2015).

[37] Iva, M., Petr, M. & Jirí, F. Dysfunction of chromatin assembly factor 1 induces shortening of telomeres and loss of 45 S rDNA in Arabidopsis thaliana. Plant Cell 22, 2768-2780 (2010).

[38] Hasegawa, J. et al. Auxin decreases chromatin accessibility through the TIR1/AFBs auxin signaling pathway in proliferative cells. Sci. Rep. 8, 7773 (2018).

[39] Aerts, N., Bruijn, S. D., Mourik, H. V., Angenent, G. C. & Dijk, A. D. J. V. Comparative analysis of binding patterns of MADS-domain proteins in Arabidopsis thaliana. Bmc Plant Biol. 18, 131 (2018).

[40] Kerstin, K. et al. Orchestration of floral initiation by APETALA1. Science 328, 85-89 (2010).

[41] Hartmann, U. et al. Molecular cloning of SVP: a negative regulator of the floral transition in Arabidopsis. Plant J. 21, 351-360 (2010).

[42] Lee, J. H., Park, S. H., Lee, J. S. & Ji, H. A. A conserved role of SHORT VEGETATIVE PHASE (SVP) in controlling flowering time of Brassica plants. Biochim. Biophys. Acta Gene Struct. Expr. 1769, 455-461 (2007).

[43] Wu, R., Wang, T., Allan, A. C., Macknight, R. C. & Varkonyi-Gasic, E. Overexpression of both AcSVP1 and AcSVP4 delays budbreak in kiwifruit A. chinensis var. deliciosa, but only AcSVP1 delays flowering in model plants. Environ. Exp. Bot. 153, 262-270 (2018).

[44] Mandel, M. A. & Yanofsky, M. F. A gene triggering flower formation in Arabidopsis. Nature 377, 522-524 (1995).

[45] Sanae, S. et al. A genetic network of flowering-time genes in wheat leaves, in which an APETALA1/FRUITFULL-like gene, VRN1, is upstream of FLOWERING LOCUS T. Plant J. 58, 668-681 (2009).

[46] Chi, Y., Fang, H., Liu, H., Yang, S. & Yu, D. An APETALA1-like gene of soybean regulates flowering time and specifies floral organs. J. Plant Physiol. 168, 2251-2259 (2011).

[47] Winterhagen, P., Tiyayon, P., Samach, A., Hegele, M. & Wünsche, J. N. Isolation and characterization of FLOWERING LOCUS T subforms and APETALA1 of the subtropical fruit tree Dimocarpus longan. Plant Physiol. Biochem. 71, 184-190 (2013).

[48] Wilschut, R. A., Oplaat, C., Snoek, L. B., Kirschner, J. & Verhoeven, K. J. F. Natural epigenetic variation contributes to heritable flowering divergence in a widespread asexual dandelion lineage. Mol. Ecol. 25, 1759-1768 (2016).

[49] Finnegan, E. J. et al. Zebularine treatment is associated with deletion of FT-B1 leading to an increase in spikelet number in bread wheat. Plant, Cell Environ. 41, 1346-1360 (2018).

[50] Iwase, Y., Shiraya, T. & Takeno, K. Flowering and dwarfism induced by DNA demethylation in Pharbitis nil. Physiologia Plant. 139, 118-127 (2010).

[51] Brown, J. C. L., De Decker, M. M. & Fieldes, M. A. A comparative analysis of developmental profiles for DNA methylation in 5-azacytidine-induced early-flowering flax lines and their control. Plant Sci. 175, 217-225 (2008).

[52] Zhang, J.-Z., Mei, L., Liu, R., Khan, M. R. G. & Hu, C.-G. Possible involvement of locus-specific methylation on expression regulation of LEAFY homologous gene (CiLFY) during precocious trifoliate orange phase change process. PLoS ONE. https://doi.org/10.1371/journal.pone.0088558 (2014).

[53] Herden, J., Eckert, S., Stift, M., Joshi, J. & van Kleunen, M . No evidence for local adaptation and an epigenetic underpinning in native and non-native ruderal plant species in Germany. Ecol. Evolution 9, 9412-9426 (2019).

[54] Baubec, T., Pecinka, A., Rozhon, W. & Mittelsten Scheid, O. Effective, homogeneous and transient interference with cytosine methylation in plant genomic DNA by zebularine. Plant J. 57, 542-554 (2009).

[55] Ramos, H. C. C., Pereira, M. G., Silva, F. F., Viana, A. P. & Ferreguetti, G. A. Seasonal and genetic influences on sex expression in a backcrossed segregating papaya population. Crop Breed. Appl. Biotechnol. 11, 97-105 (2011).

[56] Jiménez, V. M., Mora-Newcomer, E. & Gutiérrez-Soto, M. V. in Genetics and Genomics of Papaya (eds Ming, R. & Moore, P. H.) 17-33 (Springer, New York, 2014).

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

[58] Bačovský, V., Houben, A., Kumke, K. & Hobza, R. The distribution of epigenetic histone marks differs between the X and Y chromosomes in Silene latifolia. Planta. https://doi.org/10.1007/s00425-019-03182-7 (2019).