Background: Papaya is a major fruit crop in tropical and subtropical regions worldwide. It is trioecious with three sex forms: male, female, and hermaphrodite. Sex determination is controlled by a pair of nascent sex chromosomes with two slightly different Y chromosomes, Y for male and Yh for hermaphrodite. The sex chromosome genotypes are XY (male), XYh (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Yh chromosome region (HSY) is pericentromeric and heterochromatic. Physical mapping of HSY and its X counterpart is essential for sequencing these regions and uncovering the early events of sex chromosome evolution and to identify the sex determination genes for crop improvement. Results: A reiterate chromosome walking strategy was applied to construct the two physical maps with three bacterial artificial chromosome (BAC) libraries. The HSY physical map consists of 68 overlapped BACs on the minimum tiling path, and covers all four HSY-specific Knobs. One gap remained in the region of Knob 1, the only knob structure shared between HSY and X, due to the lack of HSY-specific sequences. This gap was filled on the physical map of the HSY corresponding region in the X chromosome. The X physical map consists of 44 BACs on the minimum tiling path with one gap remaining in the middle, due to the nature of highly repetitive sequences. This gap was filled on the HSY physical map. The borders of the non-recombining HSY were defined genetically by fine mapping using 1460 F2 individuals. The genetically defined HSY spanned approximately 8.5 Mb, whereas its X counterpart extended about 5.4 Mb including a 900 Kb region containing the Knob 1 shared by the HSY and X. The 8.5 Mb HSY corresponds to 4.5 Mb of its X counterpart, showing 4 Mb (89%) DNA sequence expansion. Conclusion: The 89% increase of DNA sequence in HSY indicates rapid expansion of the Yh chromosome after genetic recombination was suppressed 2–3 million years ago. The genetically defined borders coincide with the common BACs on the minimum tiling paths of HSY and X. The minimum tiling paths of HSY and its X counterpart are being used for sequencing these X and Yh-specific regions.
The papaya Y chromosome has undergone a degenerative expansion from its ancestral autosome, as a consequence of recombination suppression in the sex determining region of the sex chromosomes. The non-recombining feature led to the accumulation of repetitive sequences in the male- or hermaphrodite-specific regions of the Y or the Yh chromosome (MSY or HSY). Therefore, repeat composition and distribution in the sex determining region of papaya sex chromosomes would be informative to understand how these repetitive sequences might be involved in the early stages of sex chromosome evolution. Detailed composition of interspersed, sex-specific, and tandem repeats was analyzed from 8.1 megabases (Mb) HSY and 5.3 Mb corresponding X chromosomal regions. Approximately 77% of the HSY and 64% of the corresponding X region were occupied by repetitive sequences. Ty3-gypsy retrotransposons were the most abundant interspersed repeats in both regions. Comparative analysis of repetitive sequences between the sex determining region of papaya X chromosome and orthologous autosomal sequences of Vasconcellea monoica, a close relative of papaya lacking sex chromosomes, revealed distinctive differences in the accumulation of Ty3-Gypsy, suggesting that the evolution of the papaya sex determining region may accompany Ty3-Gypsy element accumulation. In total, 21 sex-specific repeats were identified from the sex determining region; 20 from the HSY and one from the X. Interestingly, most HSY-specific repeats were detected in two regions where the HSY expansion occurred, suggesting that the HSY expansion may result in the accumulation of sex-specific repeats or that HSY-specific repeats might play an important role in the HSY expansion. The analysis of simple sequence repeats (SSRs) revealed that longer SSRs were less abundant in the papaya sex determining region than the other chromosomal regions. Major repetitive elements were Ty3-gypsy retrotransposons in both the HSY and the corresponding X. Accumulation of Ty3-Gypsy retrotransposons in the sex determining region of papaya X chromosome was significantly higher than that in the corresponding region of V. monoica, suggesting that Ty3-Gypsy could be crucial for the expansion and evolution of the sex determining region in papaya. Most sex-specific repeats were located in the two HSY expansion regions.