Sex biased expression of hormone related genes at early stage of sex differentiation in papaya flowers

Abstract

Sex types of papaya are controlled by a pair of nascent sex chromosomes, but molecular genetic mechanisms of sex determination and sex differentiation in papaya are still unclear. We performed comparative analysis of transcriptomic profiles of male and female floral buds at the early development stage before the initiation of reproductive organ primordia at which there is no morphological difference between male and female flowers. A total of 1734 differentially expressed genes (DEGs) were identified, of which 923 showed female-biased expression and 811 showed male-biased expression. Functional annotation revealed that genes related to plant hormone biosynthesis and signaling pathways, especially in abscisic acid and auxin pathways, were overrepresented in the DEGs. Transcription factor binding motifs, such as MYB2, GAMYB, and AP2/EREBP, were enriched in the promoters of the hormone-related DEGs, and transcription factors with those motifs also exhibited differential expression between sex types. Among these DEGs, we also identified 11 genes in the non-recombining region of the papaya sex chromosomes and 9 genes involved in stamen and carpel development. Our results suggested that sex differentiation in papaya may be regulated by multiple layers of regulation and coordination and involved transcriptional, epigenetic, and phytohormone regulation. Hormones, especially ABA and auxin, transcription factors, and genes in the non-recombination region of the sex chromosome could be involved in this process. Our findings may facilitate the elucidation of signal transduction and gene interaction in sex differentiation of unisexual flowers in papaya.

References

[1] Mitchell, C. H. & Diggle, P. K. The evolution of unisexual flowers: morphological and functional convergence results from diverse developmental transitions. Am. J. Bot. 92, 1068-76 (2005).

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

[3] VanBuren, R. et al. Origin and domestication of papaya Yh chromosome. Genome Res. 25, 524-33 (2015).

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

[5] Ming, R. et al. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature 452, 991-6 (2008).

[6] Decraene, L. R. & Smets, E. The floral development and anatomy of Carica papaya (Caricaceae). Can. J. Bot. 77, 582-98 (1999).

[7] Ó'Maoiléidigh, D. S., Graciet, E. & Wellmer, F. Gene networks controlling Arabidopsis thaliana flower development. N. Phytol. 201, 16-30 (2014).

[8] Larue, N. C., Sullivan, A. M., Di & Stilio, V. S. Functional recapitulation of transitions in sexual systems by homeosis during the evolution of dioecy in Thalictrum. Front. Plant Sci. 4, 487 (2013).

[9] Pfent, C ., Pobursky, K. J., Sather, D. N. & Golenberg, E. M. Characterization of SpAPETALA3 and SpPISTILLATA, B class floral identity genes in Spinacia oleracea, and their relationship to sexual dimorphism. Dev. Genes Evol. 215, 132-42 (2005).

[10] Heikrujam, M., Sharma, K., Prasad, M. & Agrawal, V . Review on different mechanisms of sex determination and sex-linked molecular markers in dioecious crops: a current update. Euphytica 201, 161-94 (2015).

[11] Switzenberg, J. A., Little, H. A., Hammar, S. A. & Grumet, R. Floral primordia-targeted ACS (1-aminocyclopropane-1-carboxylate synthase) expression in transgenic Cucumis melo implicates fine tuning of ethylene production mediating unisexual flower development. Planta 240, 797-808 (2014).

[12] Veit, B., Schmidt, R. J., Hake, S. & Yanofsky, M. F. Maize floral development: new genes and old mutants. Plant Cell 5, 1205-15 (1993).

[13] Acosta, I. F. et al. tasselseed1 is a lipoxygenase affecting jasmonic acid signaling in sex determination of maize. Science 323, 262-5 (2009).

[14] Delong, A., Calderon-Urrea, A. & Dellaporta, S. L. Sex determination gene TASSELSEED2 of maize encodes a short-chain alcohol dehydrogenase required for stage-specific floral organ abortion. Cell 74, 757-68 (1993).

[15] Boualem, A. et al. A conserved mutation in an ethylene biosynthesis enzyme leads to andromonoecy in melons. Science 321, 836-8 (2008).

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

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

[18] Rocheta, M. et al. Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber. Front. Plant Sci. 5, 599 (2014).

[19] Guo, S. et al. Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types. BMC Genomics 11, 1-13 (2010).

[20] Sobral, R., Silva, H. G., Morais-Cecilio, L. & Costa, M. M. The quest for molecular regulation underlying unisexual flower development. Front. Plant Sci. 7, 160 (2016).

[21] 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-5 (2012).

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

[23] Ao, C. Q. Developmental origins of the conjoined twin mature embryo sacs in Smilax davidiana, with special notes on the formation of their embryos and endosperms. Am. J. Bot. 100, 2509-15 (2013).

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

[25] Du, Z., Zhou, X., Ling, Y., Zhang, Z. & Su, Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 38, W64-W70 (2010).

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

[27] Jiang, Z. et al. AHD2.0: an update version of Arabidopsis Hormone Database for plant systematic studies. Nucleic Acids Res. 39, D1123-D1129 (2011).

[28] Jin, J., Zhang, H., Kong, L., Gao, G. & Luo, J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 42, D1182-D1187 (2014).

[29] Bailey, T. L. et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, W202-W208 (2009).

[30] Farré, D. et al. Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 31, 3651-3 (2003).

[31] Ackerman, C. M. B-class MADS-box genes in trioecious papaya: two paleoAP3 paralogs, CpTM6-1 and CpTM6-2, and a PI ortholog CpPI. Planta 227, 741 (2007).

[32] Yu, Q., Steiger, D., Kramer, E. M., Moore, P. H. & Ming, R. Floral MADS-box genes in trioecious papaya: characterization of AG and AP1 subfamily genes revealed a sex-type-specific gene. Trop. Plant Biol. 1, 97-107 (2008).

[33] Westergaard, M. The mechanism of sex determination in dioecious flowering plants. Adv. Genet. 9, 217-81 (1958).

[34] Sather, D. N., Jovanovic, M. & Golenberg, E. M. Functional analysis of B and C class floral organ genes in spinach demonstrates their role in sexual dimorphism. BMC Plant Biol. 10, 46 (2010).

[35] Alvarez, J. & Smyth, D. R. CRABS CLAW and SPATULA, two Arabidopsis genes that control carpel development in parallel with AGAMOUS. Development 126, 2377-86 (1999).

[36] Fourquin, C., Primo, A., Martinez-Fernandez, I., Huet-Trujillo, E. & Ferrandiz, C. The CRC orthologue from Pisum sativum shows conserved functions in carpel morphogenesis and vascular development. Ann. Bot. 114, 1535-44 (2014).

[37] Hartwig, T. & Schulz, B. Brassinosteroid control of sex determination in maize. Proc. Natl Acad. Sci. USA 108, 19814-9 (2011).

[38] Bensen, R. J. et al. Cloning and characterization of the maize An1 gene. Plant Cell 7, 75-84 (1995).

[39] Kumar, A. & Jaiswal, V. S. Sex reversal and fruit formation on male plants of Carica Papaya L by ethrel and chlorflurenol. Proc. Indian Acad. Sci. 93, 635-41 (1984).

[40] Peterson, C. E. & Anhder, L. D. Induction of staminate flowers on gynoecious cucumbers with gibberellin A3. Science 131, 1673-4 (1960).

[41] Rudich, J. & Halevy, A. H. Involvement of abscisic acid in the regulation of sex expression in the cucumber. Plant Cell Physiol. 15, 635-42 (1974).

[42] Han, J., Murray, J. E., Yu, Q., Moore, P. H., & Ming, R. The effects of gibberellic acid on sex expression and secondary sexual characteristics in papaya. HortScience 49, 378-83 (2014).

[43] Gerashchenkov, G. A. & Rozhnova, N. A. The involvement of phytohormones in the plant sex regulation. Russian J. Plant Physiol. 60, 597-610 (2013).

[44] Tivendale, N. D., Ross, J. J. & Cohen, J. D. The shifting paradigms of auxin biosynthesis. Trends Plant Sci. 19, 44-51 (2014).

[45] Cheng, Y., Dai, X. & Zhao, Y. Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis. Genes Dev. 20, 1790-9 (2006).

[46] van Mourik, S. Simulation of organ patterning on the floral meristem using a polar auxin transport model. PloS ONE 7, e28762 (2012).

[47] Zeng, J. et al. Arabidopsis cryptochrome-1 restrains lateral roots growth by inhibiting auxin transport. J. Plant Physiol. 167, 670-3 (2010).

[48] Staswick, P. E. et al. Characterization of an Arabidopsis enzyme family that conjugates amino acids to indole-3-acetic acid. Plant Cell 17, 616-27 (2005).

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

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

[51] Costa-Nunes, J. A. D. et al. The AtRAD21.1 and AtRAD21.3 Arabidopsis cohesins play a synergistic role in somatic DNA double strand break damage repair. BMC Plant Biol. 14, 1-14 (2014).

[52] Qin, S. & Min, J. Structure and function of the nucleosome-binding PWWP domain. Trends Biochem. Sci. 39, 536-47 (2014).

[53] Chen, T., Tsujimoto, N. & Li, E. The PWWP domain of Dnmt3a and Dnmt3b is required for directing DNA methylation to the major satellite repeats at pericentric heterochromatin. Mol. Cell. Biol. 24, 9048-58 (2004).

[54] Aryal, R. Sex specific expression and distribution of small RNAs in papaya. BMC Genomics 15, 20 (2014).