Low X/Y divergence in four pairs of papaya sex-linked genes

Sex chromosomes in flowering plants, in contrast to those in animals, evolved relatively recently and only a few are heteromorphic. The homomorphic sex chromosomes of papaya show features of incipient sex chromosome evolution. We investigated the features of paired X- and Y-specific bacterial artificial chromosomes (BACs), and estimated the time of divergence in four pairs of sex-linked genes. We report the results of a comparative analysis of long contiguous genomic DNA sequences between the X and hermaphrodite Y (Y(h)) chromosomes. Numerous chromosomal rearrangements were detected in the male-specific region of the Y chromosome (MSY), including inversions, deletions, insertions, duplications and translocations, showing the dynamic evolutionary process on the MSY after recombination ceased. DNA sequence expansion was documented in the two regions of the MSY, demonstrating that the cytologically homomorphic sex chromosomes are heteromorphic at the molecular level. Analysis of sequence divergence between four X and Y(h) gene pairs resulted in a estimated age of divergence of between 0.5 and 2.2 million years, supporting a recent origin of the papaya sex chromosomes. Our findings indicate that sex chromosomes did not evolve at the family level in Caricaceae, and reinforce the theory that sex chromosomes evolve at the species level in some lineages.     

Recent Origin of Dioecious and Gynodioecious Y Chromosomes in Papaya

Sex of dioecious and gynodioecious papayas is controlled by two slightly different Y chromosomes, Y for males and Y^h for hermaphrodites. All combinations of the Y and/or Y^h chromosomes are lethal. We investigated the features of paired dioecious X- and Y-specific bacterial artificial chromosomes (BACs) and compared their sequences to corresponding gynodioecious X- and Y-specific BACs. Numerous chromosomal rearrangements were detected between the X- and Y-specific BACs, including inversions, deletions, insertions, and duplications. DNA sequence expansion was documented on the Y BAC. Dioecious and gynodioecious X-specific BACs were virtually identical. The Y- and Y^h-specific BACs shared high degree of DNA sequence identity, but local chromosomal rearrangements were detected, as the consequence of suppression of recombination in the male specific region and the isolation of Y and Y^h chromosomes enforced by the lethal effect. Analysis of sequence divergence between three dioecious X and Y gene pairs resulted in the estimated ages of divergence from 0.6 to 2.5 million years, reinforcing the hypothesis of a recent origin of the papaya sex chromosomes. The estimated age of divergence between Y and Y^h chromosomes was approximately 73,000 years for Gene 5. Our findings indicate that Y and Y^h chromosomes evolved from a common ancestral Y chromosome, possibly prior to the origin of agriculture. The existence of a hermaphrodite Y^h chromosome is less likely to have resulted from human selection as once suggested.     

Organelle DNA accumulation in the recently evolved papaya sex chromosomes

Sex chromosomes are a pair of specialized chromosomes containing a sex determination region that is suppressed for recombination. Without recombination, Y chromosomes are thought to accumulate repetitive DNA sequences which contribute to their degeneration. A pair of primitive sex chromosomes controls sex type in papaya with male and hermaphrodite determined by the slightly different male-specific region of the Y (MSY) and hermaphrodite-specific region of Y^h (HSY) chromosomes, respectively. Here, we show that the papaya HSY and MSY in the absence of recombination have accumulated nearly 12 times the amount of chloroplast-derived DNA than the corresponding region of the X chromosome and 4 times the papaya genome-wide average. Furthermore, a chloroplast genome fragment containing the rsp15 gene has been amplified 23 times in the HSY, evidence of retrotransposon-mediated duplication. Surprisingly, mitochondria-derived sequences are less abundant in the X and HSY compared to the whole genome. Shared organelle integrations are sparse between X and HSY, with only 11 % of chloroplast and 12 % of mitochondria fragments conserved, respectively, suggesting that the accelerated accumulation of organelle DNA occurred after the HSY was suppressed for recombination. Most of the organelle-derived sequences have divergence times of <7 MYA, reinforcing this notion. The accumulated chloroplast DNA is evidence of the slow degeneration of the HSY.     

Evidence for Emergence of Sex-Determining Gene(s) in a Centromeric Region in Vasconcellea parviflora

Sex chromosomes have been studied in many plant and animal species. However, few species are suitable as models to study the evolutionary histories of sex chromosomes. We previously demonstrated that papaya (Carica papaya) (2n = 2x = 18), a fruit tree in the family Caricaceae, contains recently emerged but cytologically heteromorphic X/Y chromosomes. We have been intrigued by the possible presence and evolution of sex chromosomes in other dioecious Caricaceae species. We selected a set of 22 bacterial artificial chromosome (BAC) clones that are distributed along the papaya X/Y chromosomes. These BACs were mapped to the meiotic pachytene chromosomes of Vasconcellea parviflora (2n = 2x = 18), a species that diverged from papaya ∼27 million years ago. We demonstrate that V. parviflora contains a pair of heteromorphic X/Y chromosomes that are homologous to the papaya X/Y chromosomes. The comparative mapping results revealed that the male-specific regions of the Y chromosomes (MSYs) probably initiated near the centromere of the Y chromosomes in both species. The two MSYs, however, shared only a small chromosomal domain near the centromere in otherwise rearranged chromosomes. The V. parviflora MSY expanded toward the short arm of the chromosome, whereas the papaya MSY expanded in the opposite direction. Most BACs mapped to papaya MSY were not located in V. parviflora MSY, revealing different DNA compositions in the two MSYs. These results suggest that mutation of gene(s) in the centromeric region may have triggered sex chromosome evolution in these plant species.     

Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-20) contains supplementary material, which is available to authorized users.     

Accumulation of interspersed and sex-specific repeats in the non-recombining region of papaya sex chromosomes

Background

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 Y^h 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.

Results

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.

Conclusion

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-335) contains supplementary material, which is available to authorized users.     

Construction of a sequence-tagged high-density genetic map of papaya for comparative structural and evolutionary genomics in brassicales

A high-density genetic map of papaya (Carica papaya L.) was constructed using microsatellite markers derived from BAC end sequences and whole-genome shot gun sequences. Fifty-four F(2) plants derived from varieties AU9 and SunUp were used for linkage mapping. A total of 707 markers, including 706 microsatellite loci and the morphological marker fruit flesh color, were mapped into nine major and three minor linkage groups. The resulting map spanned 1069.9 cM with an average distance of 1.5 cM between adjacent markers. This sequence-based microsatellite map resolved the very large linkage group 2 (LG 2) of the previous high-density map using amplified fragment length polymorphism markers. The nine major LGs of our map represent papaya's haploid nine chromosomes with LG 1 of the sex chromosome being the largest. This map validates the suppression of recombination at the male-specific region of the Y chromosome (MSY) mapped on LG 1 and at potential centromeric regions of other LGs. Segregation distortion was detected in a large region on LG 1 surrounding the MSY region due to the abortion of the YY genotype and in a region of LG6 due to an unknown cause. This high-density sequence-tagged genetic map is being used to integrate genetic and physical maps and to assign genome sequence scaffolds to papaya chromosomes. It provides a framework for comparative structural and evolutional genomic research in the order Brassicales.     

Genomics of sex chromosomes

Sex chromosomes in plants and animals are distinctive, not only because of their gender-determining role but also for genomic features that reflect their evolutionary history. The genomic sequences in the ancient sex chromosomes of humans and in the incipient sex chromosomes of medaka, stickleback, papaya, and poplar exhibit unusual features as consequences of their evolution. These include the enormous palindrome structure in human MSY, a duplicated genomic fragment that evolved into a Y chromosome in medaka, and a 700 kb extra telomeric sequence of the W chromosome in poplar. Comparative genomic analysis of ancient and incipient sex chromosomes highlights common features that implicate the selection forces that shaped them, even though evolutionary origin, pace, and fate vary widely among individual sex-determining systems.     

Gender in plants: sex chromosomes are emerging from the fog

Although most plants have flowers with both male and female sex organs, there are several thousands of plant species where male or female flowers form on different individuals. Surprisingly, the presence of well-established sex chromosomes in these dioecious plants is rare. The best-described example is white campion, for which large sex chromosomes have been identified and mapped partially. A recent study presented a comprehensive genetic and physical mapping of the genome of dioecious papaya. It revealed a short male specific region on the Y chromosome (MSY) that does not recombine with the X chromosome, providing strong evidence that the sex chromosomes originated from a regular pair of autosomes. The primitive papaya Y chromosome thus represents an early event in sex chromosome evolution. In this article, we review the current status of plant sex-chromosome research and discuss the advantages of different dioecious models.     

Plant sex chromosomes: molecular structure and function

Recent molecular and genomic studies carried out in a number of model dioecious plant species, including Asparagus officinalis, Carica papaya, Silene latifolia, Rumex acetosa and Marchantia polymorpha, have shed light on the molecular structure of both homomorphic and heteromorphic sex chromosomes, and also on the gene functions they have maintained since their evolution from a pair of autosomes. The molecular structure of sex chromosomes in species from different plant families represents the evolutionary pathway followed by sex chromosomes during their evolution. The degree of Y chromosome degeneration that accompanies the suppression of recombination between the Xs and Ys differs among species. The primitive Ys of A. officinalis and C. papaya have only diverged from their homomorphic Xs in a short male-specific and non-recombining region (MSY), while the heteromorphic Ys of S. latifolia, R. acetosa and M. polymorpha have diverged from their respective Xs. As in the Y chromosomes of mammals and Drosophila, the accumulation of repetitive DNA, including both transposable elements and satellite DNA, has played an important role in the divergence and size enlargement of plant Ys, and consequently in reducing gene density. Nevertheless, the degeneration process in plants does not appear to have reached the Y-linked genes. Although a low gene density has been found in the sequenced Y chromosome of M. polymorpha, most of its genes are essential and are expressed in the vegetative and reproductive organs in both male and females. Similarly, most of the Y-linked genes that have been isolated and characterized up to now in S. latifolia are housekeeping genes that have X-linked homologues, and are therefore expressed in both males and females. Only one of them seems to be degenerate with respect to its homologous region in the X. Sequence analysis of larger regions in the homomorphic X and Y chromosomes of papaya and asparagus, and also in the heteromorphic sex chromosomes of S. latifolia and R. acetosa, will reveal the degenerative changes that the Y-linked gene functions have experienced during sex chromosome evolution.     

Genetic and physical maps around the sex-determining <Emphasis Type="Italic">M</Emphasis>-locus of the dioecious plant asparagus

Asparagus officinalis L. is a dioecious plant. A region called the M-locus located on a pair of homomorphic sex chromosomes controls the sexual dimorphism in asparagus. The aim of this work was to clone the region determining sex in asparagus from its position in the genome. The structure of the region encompassing M should be investigated and compared to the sex-determining regions in other dioecious model species. To establish an improved basis for physical mapping, a high-resolution genetic map was enriched with AFLP markers closely linked to the target locus by carrying out a bulked segregant analysis. By screening a BAC library with AFLP- and STS-markers followed by chromosome walking, a physical map with eight contigs could be established. However, the gaps between the contigs could not be closed due to a plethora of repetitive elements. Surprisingly, two of the contigs on one side of the M-locus did not overlap although they have been established with two markers, which mapped in a distance as low as 0.25 cM flanking the sex locus. Thus, the clustering of the markers indicates a reduced recombination frequency within the M-region. On the opposite side of the M-locus, a contig was mapped in a distance of 0.38 cM. Four closely linked BAC clones were partially sequenced and 64 putative ORFs were identified. Interestingly, only 25% of the ORFs showed sequence similarity to known proteins and ESTs. In addition, an accumulation of repetitive sequences and a low gene density was revealed in the sex-determining region of asparagus. Molecular cytogenetic and sequence analysis of BACs flanking the M-locus indicate that the BACs contain highly repetitive sequences that localize to centromeric and pericentromeric locations on all asparagus chromosomes, which hindered the localization of the M-locus to the single pair of sex chromosomes. We speculate that dioecious Silene, papaya and Asparagus species may represent three stages in the evolution of XX, XY sex determination systems. Given that asparagus still rarely produces hermaphroditic flowers and has homomorphic sex chromosomes, this species may be an ideal system to further investigates early sex chromosome evolution and the origins of dioecy.     

Sex determination in papaya

Sex determination is an intriguing system in trioecious papaya. Over the past seven decades various hypotheses, based on the knowledge and information available at the time, have been proposed to explain the genetics of the papaya's sex determination. These include a single gene with three alleles, a group of closely linked genes, a genic balance of sex chromosome over autosomes, classical XY chromosomes, and regulatory elements of the flower development pathway. Recent advancements in genomic technology make it possible to characterize the genomic region involved in sex determination at the molecular level. High density linkage mapping validated the hypothesis that predicted recombination suppression at the sex determination locus. Physical mapping and sample sequencing of the non-recombination region led to the conclusion that sex determination is controlled by a pair of primitive sex chromosomes with a small male-specific region (MSY) of the Y chromosome. We now postulate that two sex determination genes control the sex determination pathway. One, a feminizing or stamen suppressor gene, causes stamen abortion before or at flower inception while the other, a masculinizing or carpel suppressor gene, causes carpel abortion at a later flower developmental stage. Detailed physical mapping is beginning to reveal structural details about the sex determination region and sequencing is expected to uncover candidate sex determining genes. Cloning of the sex determination genes and understanding the sex determination process could have profound application in papaya production.     

Integration of Genetic and Cytological Maps and Development of a Pachytene Chromosome-based Karyotype in Papaya

A significant amount of genetic and genomic resources have been developed in papaya (Carica papaya, $ {\hbox{2n = 2}} \times { = 18} $ ), including genetic linkage maps consisting of nine major and three minor linkage groups. However, the 12 genetic linkage groups have not been integrated with the nine chromosomes of papaya. Bacterial artificial chromosome (BAC) clones associated with each linkage group were recently isolated. These linkage group-specific BACs were mapped to meiotic pachytene chromosomes of papaya using fluorescence in situ hybridization (FISH). The FISH mapping results integrated the 12 linkage groups into the nine papaya chromosomes. We developed a pachytene chromosome-based high resolution karyotype for the hermaphrodite plant genome of papaya cultivar SunUp. The chromosomal distribution of heterochromatin in the papaya genome is provided in the karyotype with the X chromosome representing the most euchromatic chromosome in the papaya genome. FISH mapping also revealed a significant amplification of sequences related to the 5S ribosomal RNA genes, which was detected in the male-specific region of the Y chromosome, but not in the corresponding region in the X chromosome.     

Reduced Rates of Sequence Evolution of Y-Linked Satellite DNA in Rumex (Polygonaceae)

One characteristic of sex chromosomes is the accumulation of a set of different types of repetitive DNA sequences in the Y chromosomes. However, little is known about how this occurs or about how the absence of recombination affects the subsequent evolutionary fate of the repetitive sequences in the Y chromosome. Here we compare the evolutionary pathways leading to the appearance of three different families of satellite-DNA sequences within the genomes of Rumex acetosa and R. papillaris, two dioecious plant species with a complex XX/XY₁Y₂ sex-chromosome system. We have found that two of these families, one autosomic (the RAE730 family) and one Y-linked (the RAYSI family), arose independently from the ancestral duplication of the same 120-bp repeat unit. Conversely, a comparative analysis of the three satellite-DNA families reveals no evolutionary relationships between these two and the third, RAE180, also located in the Y chromosomes. However, we have demonstrated that, regardless of the mechanisms that gave rise to these families, satellite-DNA sequences have different evolutionary fates according to their location in different types of chromosomes. Specifically, those in the Y chromosomes have evolved at half the rate of those in the autosomes, our results supporting the hypothesis that satellite DNAs in nonrecombining Y chromosomes undergo lower rates of sequence evolution and homogenization than do satellite DNAs in autosomes.[Reviewing Editor: DR. Jerzy Jurka]     

Molecular cytogenetic characterization of <Emphasis Type="Italic">Rumex papillaris</Emphasis>, a dioecious plant with an XX/XY<Subscript>1</Subscript>Y<Subscript>2</Subscript> sex chromosome system

Rumex papillaris Boiss, & Reut., an Iberian endemic, belongs to the section Acetosa of the genus Rumex whose main representative is R. acetosa L., a species intensively studied in relation to sex-chromosome evolution. Here, we characterize cytogenetically the chromosomal complement of R. papillaris in an effort to enhance future comparative genomic approaches and to better our understanding of sex chromosome structure in plants. Rumex papillaris, as is common in this group, is a dioecious species characterized by the presence of a multiple sex chromosome system (with females 2n = 12 + XX and males 2n = 12 + XY₁Y₂). Except for the X chromosome both Y chromosomes are the longest in the karyotype and appear heterochromatic due to the accumulation of at least two satellite DNA families, RAE180 and RAYSI. Each chromosome of pair VI has an additional major heterochromatin block at the distal region of the short arm. These supernumerary heterochromatic blocks are occupied by RAE730 satellite DNA family. The Y-related RAE180 family is also present in an additional minor autosomal locus. Our comparative study of the chromosomal organization of the different satellite-DNA sequences in XX/XY and XX/XY₁Y₂ Rumex species demonstrates that of active mechanisms of heterochromatin amplification occurred and were accompanied by chromosomal rearrangements giving rise to the multiple XX/XY₁Y₂ chromosome systems observed in Rumex. Additionally, Y₁ and Y₂ chromosomes have undergone further rearrangements leading to differential patterns of Y-heterochromatin distribution between Rumex species with multiple sex chromosome systems.