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.     

Assignment of SRY, ANT3, and CSF2RA to the bovine Y chromosome by FISH and RH mapping

Three genes, SRY, ANT3, and CSF2RA, were mapped to the bovine Y chromosome (BTAY) by fluorescence in situ hybridization (FISH) and/or radiation hybrid (RH) mapping. FISH analysis indicated that the bovine SRY gene maps to the distal region of BTAYq, while ANT3 and CSF2RA are located in the pseudoautosomal region (PAR) of BTAYp and BTAXq. RH mapping with a 7000-rad cattle hamster whole-genome radiation hybrid panel further defined the ANT3 and CSF2RA position in relationship to previously mapped 12 PAR markers, and resulted in a relatively high resolution RH map for the PAR of BTAY.     

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.     

Chromosomal location and gene paucity of the male specific region on papaya Y chromosome

Sex chromosomes in flowering plants evolved recently and many of them remain homomorphic, including those in papaya. We investigated the chromosomal location of papaya’s small male specific region of the hermaphrodite Y (Y^h) chromosome (MSY) and its genomic features. We conducted chromosome fluorescence in situ hybridization mapping of Y^h-specific bacterial artificial chromosomes (BACs) and placed the MSY near the centromere of the papaya Y chromosome. Then we sequenced five MSY BACs to examine the genomic features of this specialized region, which resulted in the largest collection of contiguous genomic DNA sequences of a Y chromosome in flowering plants. Extreme gene paucity was observed in the papaya MSY with no functional gene identified in 715 kb MSY sequences. A high density of retroelements and local sequence duplications were detected in the MSY that is suppressed for recombination. Location of the papaya MSY near the centromere might have provided recombination suppression and fostered paucity of genes in the male specific region of the Y chromosome. Our findings provide critical information for deciphering the sex chromosomes in papaya and reference information for comparative studies of other sex chromosomes in animals and plants. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Y Chromosomal Variation Tracks the Evolution of Mating Systems in Chimpanzee and Bonobo

The male-specific regions of the Y chromosome (MSY) of the human and the chimpanzee (Pan troglodytes) are fully sequenced. The most striking difference is the dramatic rearrangement of large parts of their respective MSYs. These non-recombining regions include ampliconic gene families that are known to be important for male reproduction,and are consequently under significant selective pressure. However, whether the published Y-chromosomal pattern of ampliconic fertility genes is invariable within P. troglodytes is an open but fundamental question pertinent to discussions of the evolutionary fate of the Y chromosome in different primate mating systems. To solve this question we applied fluorescence in situ hybridisation (FISH) of testis-specific expressed ampliconic fertility genes to metaphase Y chromosomes of 17 chimpanzees derived from 11 wild-born males and 16 bonobos representing seven wild-born males. We show that of eleven P. troglodytes Y-chromosomal lines, ten Y-chromosomal variants were detected based on the number and arrangement of the ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y)—a so-far never-described variation of a species'' Y chromosome. In marked contrast, no variation was evident among seven Y-chromosomal lines of the bonobo, P. paniscus, the chimpanzee''s closest living relative. Although, loss of variation of the Y chromosome in the bonobo by a founder effect or genetic drift cannot be excluded, these contrasting patterns might be explained in the context of the species'' markedly different social and mating behaviour. In chimpanzees, multiple males copulate with a receptive female during a short period of visible anogenital swelling, and this may place significant selection on fertility genes. In bonobos, however, female mate choice may make sperm competition redundant (leading to monomorphism of fertility genes), since ovulation in this species is concealed by the prolonged anogenital swelling, and because female bonobos can occupy high-ranking positions in the group and are thus able to determine mate choice more freely.     

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

Key message

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement.

Abstract

To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1U-6U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.     

Testing for the Footprint of Sexually Antagonistic Polymorphisms in the Pseudoautosomal Region of a Plant Sex Chromosome Pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.     

Extreme population structure and high interspecific divergence of the Silene Y chromosome

Previous studies have demonstrated that the diversity of Y-linked genes is substantially lower than that of their X-linked homologs in the plant Silene latifolia. This difference has been attributed to selective sweeps, Muller's ratchet, and background selection, processes that are predicted to severely affect the evolution of the nonrecombining Y chromosome. We studied the DNA diversity of a noncoding region of the homologous genes DD44Y and DD44X, sampling S. latifolia populations from a wide geographical area and also including the closely related species S. dioica, S. diclinis, and S. heuffelii. On the Y chromosome of S. latifolia, we found substantial DNA diversity. Geographical population structure was far higher than on the X chromosome and differentiation between the species was also higher for the Y than for the X chromosome. Our findings indicate that the loss of genetic diversity on the Y chromosome in Silene occurs within local populations rather than within entire species. These results are compatible with background selection, Muller's ratchet, and local selective sweeps, but not with species-wide selective sweeps. The higher interspecific divergence of DD44Y, compared to DD44X, supports the hypothesis that Y chromosome differentiation between incipient species precedes reproductive isolation of the entire genome, forming an early stage in the process of speciation.     

Gene Finding on the Y: Fruitful Strategy in Drosophila does not Deliver in Anopheles

The Anopheles gambiae genome project yielded almost complete sequences for the autosomes and for a large part of the X chromosome, however, no information for the Y chromosome was obtained. Yet, by design, fragmented Y chromosome sequences should be present in the resulting assembly. Here we report the search for Anopheles Y chromosome genes using a strategy successfully applied for identification of Y genes in Drosophila. A complete set of the unmapped scaffolds was targeted in a broad TBLASTN search using both A. gambiae predicted genes and all proteins from nr database as query sequences. After filtering of the BLAST report, we selected 181 scaffolds possibly containing fragments of Y chromosome genes to experimentally test their Y-linkage. Surprisingly, none of the tested sequences appeared to originate from the Y chromosome. Several factors could account for the failure to detect Y genes, including their different organization in A. gambiae compared to Drosophila and the suboptimal quality of the assembly and annotation of the Anopheles genome. Regardless of the cause, our results illuminate problems associated with the genome analysis of outbred organisms.     

Assembly of the Genome of the Disease Vector Aedes aegypti onto a Genetic Linkage Map Allows Mapping of Genes Affecting Disease Transmission

The mosquito Aedes aegypti transmits some of the most important human arboviruses, including dengue, yellow fever and chikungunya viruses. It has a large genome containing many repetitive sequences, which has resulted in the genome being poorly assembled — there are 4,758 scaffolds, few of which have been assigned to a chromosome. To allow the mapping of genes affecting disease transmission, we have improved the genome assembly by scoring a large number of SNPs in recombinant progeny from a cross between two strains of Ae. aegypti, and used these to generate a genetic map. This revealed a high rate of misassemblies in the current genome, where, for example, sequences from different chromosomes were found on the same scaffold. Once these were corrected, we were able to assign 60% of the genome sequence to chromosomes and approximately order the scaffolds along the chromosome. We found that there are very large regions of suppressed recombination around the centromeres, which can extend to as much as 47% of the chromosome. To illustrate the utility of this new genome assembly, we mapped a gene that makes Ae. aegypti resistant to the human parasite Brugia malayi, and generated a list of candidate genes that could be affecting the trait.     

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.     

Identification in chromosome 8q11 of a region of homology with the g1 amplicon of the Y chromosome and functional analysis of the BEYLA gene

The male-specific region (MSY) of the Y chromosome contains genes involved mainly in male sex determination and in spermatogenesis. The majority of genes involved in male fertility are localized in multiple copies in the long arm of the Y chromosome, within specific regions defined as "ampliconic regions." It has been suggested that these genes derived from X-linked or autosomal ancestors during evolution, providing a benefit for male fertility when transposed onto the Y chromosome. So far, the autosomal origin has been demonstrated only for two MSY genes, DAZ and CDY. In the present study we report on the identification within chromosome 8q11.2 of a region homologous to the g amplicon, containing the VCY2 (approved gene symbol BPY2), TTTY4, and TTTY17 genes. A search for ancestor genes within the 8q11.2 region allowed us to identify a gene named BEYLA and to characterize the genomic organization and the expression patterns of this gene.     

A Distinct Type of Heterochromatin at the Telomeric Region of the Drosophila melanogaster Y Chromosome

Heterochromatin assembly and its associated phenotype, position effect variegation (PEV), provide an informative system to study chromatin structure and genome packaging. In the fruit fly Drosophila melanogaster, the Y chromosome is entirely heterochromatic in all cell types except the male germline; as such, Y chromosome dosage is a potent modifier of PEV. However, neither Y heterochromatin composition, nor its assembly, has been carefully studied. Here, we report the mapping and characterization of eight reporter lines that show male-specific PEV. In all eight cases, the reporter insertion sites lie in the telomeric transposon array (HeT-A and TART-B2 homologous repeats) of the Y chromosome short arm (Ys). Investigations of the impact on the PEV phenotype of mutations in known heterochromatin proteins (i.e., modifiers of PEV) show that this Ys telomeric region is a unique heterochromatin domain: it displays sensitivity to mutations in HP1a, EGG and SU(VAR)3-9, but no sensitivity to Su(z)2 mutations. It appears that the endo-siRNA pathway plays a major targeting role for this domain. Interestingly, an ectopic copy of 1360 is sufficient to induce a piRNA targeting mechanism to further enhance silencing of a reporter cytologically localized to the Ys telomere. These results demonstrate the diversity of heterochromatin domains, and the corresponding variation in potential targeting mechanisms.     

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.