An integrated cytogenetic and physical map reveals unevenly distributed recombination spots along the papaya sex chromosomes

Descriptive Genomic Analysis and Sequence Genotyping of the Two Papaya Species (Vasconcellea pubescens and Vasconcellea chilensis) Using GBS Tools

A genotyping by sequencing (GBS) approach was used to analyze the organization of genetic diversity in V. pubescens and V. chilensis. GBS identified 4675 and 4451 SNPs/INDELs in two papaya species. The cultivated orchards of V. pubescens exhibited scarce genetic diversity and low but significant genetic differentiation. The neutrality test yielded a negative and significant result, suggesting that V. pubescens suffered a selective sweep or a rapid expansion after a bottleneck during domestication. In contrast, V. chilensis exhibited a high level of genetic diversity. The genetic differentiation among the populations was slight, but it was possible to distinguish the two genetic groups. The neutrality test indicated no evidence that natural selection and genetic drift affect the natural population of V. chilensis. Using the Carica papaya genome as a reference, we identified critical SNPs/INDELs associated with putative genes. Most of the identified genes are related to stress responses (salt and nematode) and vegetative and reproductive development. These results will be helpful for future breeding and conservation programs of the Caricaceae family.

Some thoughts about the words we use for thinking about sex chromosome evolution

Sex chromosomes are familiar to most biologists since they first learned about genetics. However, research over the past 100 years has revealed that different organisms have evolved sex-determining systems independently. The differences in the ages of systems, and in how they evolved, both affect whether sex chromosomes have evolved. However, the diversity means that the terminology used tends to emphasize either the similarities or the differences, sometimes causing misunderstandings. In this article, I discuss some concepts where special care is needed with terminology. The following four terms regularly create problems: ‘sex chromosome’, ‘master sex-determining gene’, ‘evolutionary strata’ and ‘genetic degeneration’. There is no generally correct or wrong use of these words, but efforts are necessary to make clear how they are to be understood in specific situations. I briefly outline some widely accepted ideas about sex chromosomes, and then discuss these ‘problem terms’, highlighting some examples where careful use of the words helps bring to light current uncertainties and interesting questions for future work.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Analysis of Meiosis in Nonmodel Tropical Plants: The Case of Carica papaya Linn

To develop plants that are more tolerant to drought, marginal soil fertility, and diseases and that satisfy demands for high yield, new cultivars of the tropical fruit papaya (Carica papaya L.) are needed. Nonetheless, in many cases, these traits are available in only wild relatives found throughout Latin America. Understanding meiotic progression may facilitate the introgression of desirable traits into commercial cultivars that maintain high fertility. In this protocol, we describe a practical and simple method to effectively isolate male meiocytes in order to document the behavior of papaya meiotic chromosomes.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

The Guppy Sex Chromosome System and the Sexually Antagonistic Polymorphism Hypothesis for Y Chromosome Recombination Suppression

Sex chromosomes regularly evolve suppressed recombination, distinguishing them from other chromosomes, and the reason for this has been debated for many years. It is now clear that non-recombining sex-linked regions have arisen in different ways in different organisms. A major hypothesis is that a sex-determining gene arises on a chromosome and that sexually antagonistic (SA) selection (sometimes called intra-locus sexual conflict) acting at a linked gene has led to the evolution of recombination suppression in the region, to reduce the frequency of low fitness recombinant genotypes produced. The sex chromosome system of the guppy (Poecilia reticulata) is often cited as supporting this hypothesis because SA selection has been demonstrated to act on male coloration in natural populations of this fish, and probably contributes to maintaining polymorphisms for the genetic factors involved. I review classical genetic and new molecular genetic results from the guppy, and other fish, including approaches for identifying the genome regions carrying sex-determining loci, and suggest that the guppy may exemplify a recently proposed route to sex chromosome evolution.

Extremely low nucleotide diversity in the X-linked region of papaya caused by a strong selective sweep

Background

The papaya Y-linked region showed clear population structure, resulting in the detection of the ancestral male population that domesticated hermaphrodite papayas were selected from. The same populations were used to study nucleotide diversity and population structure in the X-linked region.

Results

Diversity is very low for all genes in the X-linked region in the wild dioecious population, with nucleotide diversity π_syn = 0.00017, tenfold lower than the autosomal region (π_syn = 0.0017) and 12-fold lower than the Y-linked region (π_syn = 0.0021). Analysis of the X-linked sequences shows an undivided population, suggesting a geographically wide diversity-reducing event, whereas two subpopulations were observed in the autosomes separating gynodioecy and dioecy and three subpopulations in the Y-linked region separating three male populations. The extremely low diversity in the papaya X-linked region was probably caused by a recent, strong selective sweep before domestication, involving either the spread of a recessive mutation in an X-linked gene that is beneficial to males or a partially dominant mutation that benefitted females or both sexes. Nucleotide diversity in the domesticated X samples is about half that in the wild Xs, probably due to the bottleneck when hermaphrodites were selected during domestication.

Conclusions

The extreme low nucleotide diversity in the papaya X-linked region is much greater than observed in humans, great apes, and the neo-X chromosome of Drosophila miranda, which show the expected pattern of Y-linked genes < X-linked genes < autosomal genes; papaya shows an unprecedented pattern of X-linked genes < autosomal genes < Y-linked genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1095-9) contains supplementary material, which is available to authorized users.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga Ectocarpus Exhibit Unusual Features

The recombining regions of sex chromosomes (pseudoautosomal regions, PARs) are predicted to exhibit unusual features due to their being genetically linked to the nonrecombining, sex-determining region. This phenomenon is expected to occur in both diploid (XY, ZW) and haploid (UV) sexual systems, with slightly different consequences for UV sexual systems because of the absence of masking during the haploid phase (when sex is expressed) and because there is no homozygous sex in these systems. Despite a considerable amount of theoretical work on PAR genetics and evolution, these genomic regions have remained poorly characterized empirically. We show here that although the PARs of the U/V sex chromosomes of the brown alga Ectocarpus recombine at a similar rate to autosomal regions of the genome, they exhibit many genomic features typical of nonrecombining regions. The PARs were enriched in clusters of genes that are preferentially, and often exclusively, expressed during the sporophyte generation of the life cycle, and many of these genes appear to have evolved since the Ectocarpales diverged from other brown algal lineages. A modeling-based approach was used to investigate possible evolutionary mechanisms underlying this enrichment in sporophyte-biased genes. Our results are consistent with the evolution of the PAR in haploid systems being influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

Construction of cytogenetic map of Gossypium herbaceum chromosome 1 and its integration with genetic maps

Background

Cytogenetic map can provide not only information of the genome structure, but also can build a solid foundation for genetic research. With the developments of molecular and cytogenetic studies in cotton (Gossypium), the construction of cytogenetic map is becoming more and more imperative.

Results

A cytogenetic map of chromosome 1 (A₁01) of Gossypium herbaceum (A₁) which includes 10 bacterial artificial chromosome (BAC) clones was constructed by using fluorescent in situ hybridization (FISH). Meanwhile, comparison and analysis were made for the cytogenetic map of chromosome 1 (A₁01) of G. herbaceum with four genetic linkage maps of chromosome 1 (A_h01) of G. hirsutum ((AD)₁) and one genetic linkage map of chromosome 1 of (A₁01) G. arboreum (A₂). The 10 BAC clones were also used to be localized on G. raimondii (D₅) chromosome 1 (D₅01), and 2 of them showed clear unique hybridized signals. Furthermore, these 2 BAC clones were also shown localized on chromosome 1 of both A sub-genome and D sub-genome of G. hirsutum.

Conclusion

The comparison of the cytogenetic map with genetic linkage maps showed that most of the identified marker-tagged BAC clones appearing same orders in different maps except three markers showing different positions, which might indicate chromosomal segmental rearrangements. The positions of the 2 BAC clones which were localized on A_h01 and D_h01 chromosomes were almost the same as that on A₁01 and D₅01 chromosomes. The corresponding anchored SSR markers of these 2 BAC clones were firstly found to be localized on chromosome D₅01 (D_h01) as they were not seen mapped like this in any genetic map reported.

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.

Fluorescence in situ hybridization and optical mapping to correct scaffold arrangement in the tomato genome

The order and orientation (arrangement) of all 91 sequenced scaffolds in the 12 pseudomolecules of the recently published tomato (Solanum lycopersicum, 2n = 2x = 24) genome sequence were positioned based on marker order in a high-density linkage map. Here, we report the arrangement of these scaffolds determined by two independent physical methods, bacterial artificial chromosome–fluorescence in situ hybridization (BAC-FISH) and optical mapping. By localizing BACs at the ends of scaffolds to spreads of tomato synaptonemal complexes (pachytene chromosomes), we showed that 45 scaffolds, representing one-third of the tomato genome, were arranged differently than predicted by the linkage map. These scaffolds occur mostly in pericentric heterochromatin where 77% of the tomato genome is located and where linkage mapping is less accurate due to reduced crossing over. Although useful for only part of the genome, optical mapping results were in complete agreement with scaffold arrangement by FISH but often disagreed with scaffold arrangement based on the linkage map. The scaffold arrangement based on FISH and optical mapping changes the positions of hundreds of markers in the linkage map, especially in heterochromatin. These results suggest that similar errors exist in pseudomolecules from other large genomes that have been assembled using only linkage maps to predict scaffold arrangement, and these errors can be corrected using FISH and/or optical mapping. Of note, BAC-FISH also permits estimates of the sizes of gaps between scaffolds, and unanchored BACs are often visualized by FISH in gaps between scaffolds and thus represent starting points for filling these gaps.

Strong accumulation of chloroplast DNA in the Y chromosomes of Rumex acetosa and Silene latifolia

Chloroplast DNA (cpDNA) sequences are often found in plant nuclear genomes, but patterns of their chromosomal distribution are not fully understood. The distribution of cpDNA on the sex chromosomes can only be studied in dioecious plant species possessing heteromorphic sex chromosomes. We reconstructed the whole chloroplast genome of Rumex acetosa (sorrel, XY1Y2 system) from next generation sequencing data. We systematically mapped the chromosomal localization of various regions of cpDNA in R. acetosa and in Silene latifolia (white campion, XY system) using fluorescence in situ hybridization. We found that cpDNA was accumulated on the Y chromosomes of both studied species. In R. acetosa, the entire Y chromosome gathered all parts of cpDNA equally. On the contrary, in S. latifolia, the majority of the cpDNA, corresponding to the single copy regions, was localized in the centromere of the Y chromosome, while the inverted repeat region was present also in other loci. We found a stronger accumulation of cpDNA on the more degenerated Y1 and Y2 chromosomes of R. acetosa than in evolutionary younger S. latifolia Y chromosome. Our data stressed the prominent role of the Y chromosome centromere in cpDNA accumulation.

A high-resolution cucumber cytogenetic map integrated with the genome assembly

Background

High-resolution cytogenetic map can provide not only important biological information on genome organization but also solid foundation for genetic and genomic research. The progress in the molecular and cytogenetic studies has created the basis for developing the cytogenetic map in cucumber (Cucumis sativus L.).

Results

Here, the cytogenetic maps of four cucumber chromosomes (chromosomes 1, 3–5) were constructed by fluorescence in situ hybridization (FISH) analysis on cucumber pachytene chromosomes. Together with our previously constructed cytogenetic maps of three cucumber chromosomes (chromosomes 2, 6–7), cucumber has a complete cytogenetic map with 76 anchoring points between the genetic, the cytogenetic and the draft genome assembly maps. To compare our pachytene FISH map directly to the genetic linkage and draft genome assembly maps, we used a standardized map unit—relative map position (RMP) to produce the comparative map alignments. The alignments allowed a global view of the relationship of genetic and physical distances along each cucumber chromosome, and accuracy and coverage of the draft genome assembly map.

Conclusions

We demonstrated a good correlation between positions of the markers in the linkage and physical maps, and essentially complete coverage of chromosome arms by the draft genome assembly. Our study not only provides essential information for the improvement of sequence assembly but also offers molecular tools for cucumber genomics research, comparative genomics and evolutionary study.

Integration of high-resolution physical and genetic map reveals differential recombination frequency between chromosomes and the genome assembling quality in cucumber

Cucumber is an important model crop and the first species sequenced in Cucurbitaceae family. Compared to the fast increasing genetic and genomics resources, the molecular cytogenetic researches in cucumber are still very limited, which results in directly the shortage of relation between plenty of physical sequences or genetic data and chromosome structure. We mapped twenty-three fosmids anchored by SSR markers from LG-3, the longest linkage group, and LG-4, the shortest linkage group on pachytene chromosomes 3 and 4, using uorescence in situ hybridization (FISH). Integrated molecular cytogenetic maps of chromosomes 3 and 4 were constructed. Except for three SSR markers located on heterochromatin region, the cytological order of markers was concordant with those on the linkage maps. Distinct structural differences between chromosomes 3 and 4 were revealed by the high resolution pachytene chromosomes. The extreme difference of genetic length between LG-3 and LG-4 was mainly attributed to the difference of overall recombination frequency. The significant differentiation of heterochromatin contents in chromosomes 3 and 4 might have a direct correlation with recombination frequency. Meanwhile, the uneven distribution of recombination frequency along chromosome 4 was observed, and recombination frequency of the long arm was nearly 3.5 times higher than that of the short arm. The severe suppression of recombination was exhibited in centromeric and heterochromatin domains of chromosome 4. Whereas a close correlation between the gene density and recombination frequency was observed in chromosome 4, no significant correlation was observed between them along chromosome 3. The comparison between cytogenetic and sequence maps revealed a large gap on the pericentromeric heterochromatin region of sequence map of chromosome 4. These results showed that integrated molecular cytogenetic maps can provide important information for the study of genetic and genomics in cucumber.