Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla,But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

Localization of the active gene of aldolase on chromosome 16, and two aldolase A pseudogenes on chromosomes 3 and 10

Summary Southern blot analysis of human genomic DNA hybridized with a coding region aldolase A cDNA probe (600 bases) revealed four restriction fragments with EcoRI restriction enzyme: 7.8 kb, 13 kb, 17 kb and >30 kb. By human-hamster hybrid analysis (Southern technique) the principal fragments, 7.8 kb, 13 kb, >30 kb, were localized to chromosomes 10, 16 and 3 respectively. The 17-kb fragment was very weak in intensity; it co-segregated with the >30-kb fragment and is probably localized on chromosome 3 with the >30-kb fragment. Analysis of a second aldolase A labelled probe protected against S1 nuclease digestion by RNAs from different hybrid cells, indicated the presence of aldolase A mRNAs in hybrid cells containing only chromosome 16. Under the stringency conditions used, the EcoRI sequences detected by the coding region aldolase A cDNA probe did not correspond to aldolase B or C. The 7.8-kb and >30-kb EcoRI sequences, localized respectively on chromosomes 10 and 3, correspond to aldolase A pseudogenes, the 13-kb EcoRI sequence localized on chromosome 16 corresponds to the aldolase active gene. The fact that the aldolase A gene and pseudogenes are located on three different chromosomes supports the hypothesis that the pseudogenes originated from aldolase A mRNAs, copied into DNA and integrated in unrelated chromosomal loci.     

Initiation sites for human DNA replication at a putative ribulose-5-phosphate 3-epimerase gene

Replication of the human genome requires the activation of thousands of replicons distributed along each one of the chromosomes. Each replicon contains an initiation, or origin, site, at which DNA synthesis begins. However, very little information is known about the nature and positioning of these initiation sites along human chromosomes. We have recently focused our attention to a 1.1 kb region of human chromosome 2 which functioned as an episomal origin in the yeast Saccharomyces cerevisiae. This region corresponded to the largest exon of a putative ribulose-5-phosphate-3-epimerase gene (RPE). In the present study we have used a real-time PCR-based nascent strand DNA abundance assay to map initiation sites for DNA replication in in vivo human chromosomes around a 13.4 kb region encompassing the putative RPE gene. By applying this analysis to a 1-1.4 kb nascent strand DNA fraction isolated from both normal skin fibroblasts, and the breast cell line MCF10; we have identified five initiation sites within the 13.4 kb region of chromosome 2. The initiation sites appear to map to similar positions in both cell lines and occur outside the coding regions of the putative RPE gene.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Asynchronous Replication,Mono-Allelic Expression,and Long Range Cis-Effects of ASAR6

Mammalian chromosomes initiate DNA replication at multiple sites along their length during each S phase following a temporal replication program. The majority of genes on homologous chromosomes replicate synchronously. However, mono-allelically expressed genes such as imprinted genes, allelically excluded genes, and genes on female X chromosomes replicate asynchronously. We have identified a cis-acting locus on human chromosome 6 that controls this replication-timing program. This locus encodes a large intergenic non-coding RNA gene named Asynchronous replication and Autosomal RNA on chromosome 6, or ASAR6. Disruption of ASAR6 results in delayed replication, delayed mitotic chromosome condensation, and activation of the previously silent alleles of mono-allelic genes on chromosome 6. The ASAR6 gene resides within an ∼1.2 megabase domain of asynchronously replicating DNA that is coordinated with other random asynchronously replicating loci along chromosome 6. In contrast to other nearby mono-allelic genes, ASAR6 RNA is expressed from the later-replicating allele. ASAR6 RNA is synthesized by RNA Polymerase II, is not polyadenlyated, is restricted to the nucleus, and is subject to random mono-allelic expression. Disruption of ASAR6 leads to the formation of bridged chromosomes, micronuclei, and structural instability of chromosome 6. Finally, ectopic integration of cloned genomic DNA containing ASAR6 causes delayed replication of entire mouse chromosomes.     

Human U1 small nuclear RNA genes: extensive conservation of flanking sequences suggests cycles of gene amplification and transposition.

The DNA immediately flanking the 164-base-pair U1 RNA coding region is highly conserved among the approximately 30 human U1 genes. The U1 multigene family also contains many U1 pseudogenes (designated class I) with striking although imperfect flanking homology to the true U1 genes. Using cosmid vectors, we now have cloned, characterized, and partially sequenced three 35-kilobase (kb) regions of the human genome spanning U1 homologies. Two clones contain one true U1 gene each, and the third bears two class I pseudogenes 9 kb apart in the opposite orientation. We show by genomic blotting and by direct DNA sequence determination that the conserved sequences surrounding U1 genes are much more extensive than previously estimated: nearly perfect sequence homology between many true U1 genes extends for at least 24 kb upstream and at least 20 kb downstream from the U1 coding region. In addition, the sequences of the two new pseudogenes provide evidence that class I U1 pseudogenes are more closely related to each other than to true genes. Finally, it is demonstrated elsewhere (Lindgren et al., Mol. Cell. Biol. 5:2190-2196, 1985) that both true U1 genes and class I U1 pseudogenes map to chromosome 1, but in separate clusters located far apart on opposite sides of the centromere. Taken together, these results suggest a model for the evolution of the U1 multigene family. We speculate that the contemporary family of true U1 genes was derived from a more ancient family of U1 genes (now class I U1 pseudogenes) by gene amplification and transposition. Gene amplification provides the simplest explanation for the clustering of both U1 genes and class I pseudogenes and for the conservation of at least 44 kb of DNA flanking the U1 coding region in a large fraction of the 30 true U1 genes.     

Juxta-centromeric region of human chromosome 21 is enriched for pseudogenes and gene fragments

A physical map including four pseudogenes and 10 gene fragments and spanning 500 kb in the juxta-centromeric region of the long arm of human chromosome 21 is presented. cDNA fragments isolated from a selected cDNA library were characterized and mapped to the 831B6 YAC and to two BAC contigs that cover 250 kb of the region. An 85 kb genomic sequence located in the proximal region of the map was analyzed for putative exons. Four pseudogenes were found, including psiIGSF3, psiEIF3, psiGCT-rel whose functional copies map to chromosome 1p13, chromosome 2 and chromosome 22q11, respectively. The TTLL1 pseudogene corresponds to a new gene whose functional copy maps to chromosome 22q13. Ten gene fragments represent novel sequences that have related sequences on different human chromosomes and show 97-100% nucleotide identity to chromosome 21. These may correspond to pseudogenes on chromosome 21 and to functional genes in other chromosomes. The 85 kb genomic sequence was analyzed also for GC content, CpG islands, and repetitive sequence distribution. A GC-poor L isochore spanning 40 kb from satellite 1 was observed in the most centromeric region, next to a GC-rich H isochore that is a candidate region for the presence of functional genes. The pericentric duplication of a 7.8 kb region that is derived from the 22q13 chromosome band is described. We showed that the juxta-centromeric region of human chromosome 21 is enriched for retrotransposed pseudogenes and gene fragments transferred by interchromosome duplications, but we do not rule out the possibility that the region harbors functional genes also.     

Detection of single-copy genes and chromosome rearrangements in Petunia hybrida by fluorescence in situ hybridization

DNA sequences homologous to single-copy genes were labelled with biotinylated dUTP or digoxygenin-labelled dUTP and hybridized to chromosome spreads. The hybridization signals were visualized with fluorescent avidin- or antibody-conjugates. This method allowed the detection of DNA targets on metaphase chromosomes as small as 1.4 kb. The hybridization signals were identified as fluorescent spots on both sister chromatids. Using an 18S rDNA probe as marker to identify chromosomes II and III it was possible to assign single-copy genes to these chromosomes. In the line V30 the endogenous chalcone synthase gene (chsA) was mapped at the distal end of the short arm of chromosome 5. The cDNA probe for this single-copy gene was 1.4 kb. In contrast, in the lines Mitchell and V26 chsA was localized at the distal end of the long arm of chromosome 3, suggesting that a chromosomal rearrangement had taken place. In a transformed Petunia uidA, transgenes were detected using a 2.7 kb probe. One transgene was mapped on one of the homologues of chromosome II proximal to the ribosomal genes. This homologue could be distinguished from the other by having the ribosomal genes at the distal end of the long arm. Using multicolour fluorescence in situ hybridization it was shown that it is possible to detect the endogenous chsA genes and both transgenes simultaneously.     

Structure, sequence, chromosomal location, and evolution of the human ferredoxin gene family

Ferredoxin is an iron-sulfur protein that serves as an electron transport intermediate for mitochondrial cytochromes P450 involved in steroid, vitamin D, and bile acid metabolism. We cloned and characterized the human ferredoxin gene family, which includes two expressed genes and two pseudogenes. Sequence analysis of this gene family revealed that it encodes only one protein product. The expressed genes were assigned to chromosome 11 and pseudogenes to chromosomes 20 and 21 by identifying single-copy probes from each gene segment and hybridizing them to DNA from rodent-human hybrid cells. The pseudogenes lacked introns and contained numerous mutations, including insertion, deletion, and substitution which rendered them inactive. They were 96% and 85% homologous to the expressed gene, yet they were only 78% homologous with each other. The intronless nature, higher diversity among themselves, and distinct chromosomal location of the pseudogenes suggests that they arose by independent, retroposon-mediated events.     

Origin and spread of the SRY gene on the X and Y chromosomes of the rodent Microtus cabrerae: role of L1 elements

In the rodent species Microtus cabrerae, males as well as females present several copies of the SRY gene, a single-copy gene located on the Y chromosome in most mammals. Using different PCR approaches, we have characterized the sequence, structure, and organization of the SRY copies and their flanking regions distributed on the X and Y chromosomes of this species. All copies of SRY analyzed, including those from the Y chromosome, proved to be nonfunctional pseudogenes, as they have internal stop codons. In addition, we demonstrated the association of SRY pseudogenes with different fragments of L1 and LTR retroelements in both sex chromosomes of M. cabrerae. Examining the possible origin of SRY pseudogene and retroposons association, we propose that retroposons could have been involved in the mechanism of SRY gene amplification on the Y chromosome and in the transference of the Y-linked SRY copies to the X-chromosome heterochromatin.     

A revised root for the human Y chromosomal phylogenetic tree: the origin of patrilineal diversity in Africa

To shed light on the structure of the basal backbone of the human Y chromosome phylogeny, we sequenced about 200 kb of the male-specific region of the human Y chromosome (MSY) from each of seven Y chromosomes belonging to clades A1, A2, A3, and BT. We detected 146 biallelic variant sites through this analysis. We used these variants to construct a patrilineal tree, without taking into account any previously reported information regarding the phylogenetic relationships among the seven Y chromosomes here analyzed. There are several key changes at the basal nodes as compared with the most recent reference Y chromosome tree. A different position of the root was determined, with important implications for the origin of human Y chromosome diversity. An estimate of 142 KY was obtained for the coalescence time of the revised MSY tree, which is earlier than that obtained in previous studies and easier to reconcile with plausible scenarios of modern human origin. The number of deep branchings leading to African-specific clades has doubled, further strengthening the MSY-based evidence for a modern human origin in the African continent. An analysis of 2204 African DNA samples showed that the deepest clades of the revised MSY phylogeny are currently found in central and northwest Africa, opening new perspectives on early human presence in the continent.     

Structure of the major block of alphoid satellite DNA on the human Y chromosome

Alphoid DNA is a family of tandemly repeated simple sequences found mainly at the centromeres of the chromosomes of many primates. This paper describes the structure of the alphoid DNA at the centromere of the human Y chromosome. We have used pulsedfield gradient gel electrophoresis, cosmid cloning and DNA sequencing to determine the organization of the alphoid DNA on each of the Y chromosomes present in two somatic cell hybrids. In each case there is a single major block of alphoid DNA. This is approximately 470,000 bases (475 kb) long on one chromosome and approximately 575 kb long on the other. Apart from the size difference, the structures of the two blocks and the surrounding sequences are very similar. However, one restriction enzyme, AvaII, detects two clusters of sites within one block but does not cleave the other. The alphoid DNA within each block is organized into tandemly repeating units, most of which are about 5.7 kb long. A few variant units present on one chromosome are about 6.0 kb long. These variants, like the AvaII site variants, are clustered. The 5.7 kb and 6.0 kb units themselves consist of tandemly repeating 170 base-pair subunits. The 6.0 kb unit has two more of these subunits than the 5.7 kb unit. Our results provide a basis for further structural analysis of the human Y chromosome centromeric region, and suggest that long-range structural polymorphisms of tandemly repeated sequence families may be frequent.     

The mammalian Y chromosome: a new perspective

For several decades, the mammalian Y chromosome was considered a genetic “desert,” with the testis determinant being the sole survivor of the attrition that followed the chromosome's inception. Aside from the addition of a genetic factor required for spermatogenesis to the human Y chromosome in 1976, this view held sway until the mid-1980s. The ensuing molecular genetic analysis, culminating in the recent paper in Science by Lahn and Page,¹ has identified more than 20 genes or gene families on the human Y. This has led to a reappraisal of the evolution and functions of this unique chromosome. BioEssays 20 :363–366, 1998. © 1998 John Wiley & Sons Inc.     

An African American Paternal Lineage Adds an Extremely Ancient Root to the Human Y Chromosome Phylogenetic Tree

We report the discovery of an African American Y chromosome that carries the ancestral state of all SNPs that defined the basal portion of the Y chromosome phylogenetic tree. We sequenced ∼240 kb of this chromosome to identify private, derived mutations on this lineage, which we named A00. We then estimated the time to the most recent common ancestor (TMRCA) for the Y tree as 338 thousand years ago (kya) (95% confidence interval = 237–581 kya). Remarkably, this exceeds current estimates of the mtDNA TMRCA, as well as those of the age of the oldest anatomically modern human fossils. The extremely ancient age combined with the rarity of the A00 lineage, which we also find at very low frequency in central Africa, point to the importance of considering more complex models for the origin of Y chromosome diversity. These models include ancient population structure and the possibility of archaic introgression of Y chromosomes into anatomically modern humans. The A00 lineage was discovered in a large database of consumer samples of African Americans and has not been identified in traditional hunter-gatherer populations from sub-Saharan Africa. This underscores how the stochastic nature of the genealogical process can affect inference from a single locus and warrants caution during the interpretation of the geographic location of divergent branches of the Y chromosome phylogenetic tree for the elucidation of human origins.     

Phylogenetic assessment of introns and SINEs within the Y chromosome using the cat family felidae as a species tree

The cat family Felidae was used as a species tree to assess the phylogenetic performance of genes, and their embedded SINE elements, within the nonrecombining region of the Y chromosome (NRY). Genomic segments from single-copy X-Y homologs SMCY, UBE1Y, and ZFY (3,604 bp) were amplified in 36 species of cat. These genes are located within the X-degenerate region of the NRY and are thought to be molecular "fossils" that ceased conventional recombination with the X chromosome early within the placental mammal evolution. The pattern and tempo of evolution at these three genes is significant in light of the recent, rapid evolution of the family over approximately 12 Myr and provides exceptional support for each of the eight recognized felid lineages, as well as clear diagnostic substitutions identifying nearly all species. Bootstrap support and Bayesian posterior probabilities are uniformly high for defining each of the eight monophyletic lineages. Further, the preferential use of specific target-site motifs facilitating SINE insertion is empirically supported by sequence analyses of SINEs embedded within the three genes. Target-site insertion is thought to explain the contradiction between intron phylogeny and results of the SMCY SINE phylogeny that unites distantly related species. Overall, our data suggest X-degenerate genes within the NRY are singularly powerful markers and offer a valuable patrilineal perspective in species evolution.