A gradient of silent substitution rate in the human pseudoautosomal region

It has been demonstrated that recombination in the human p-arm pseudoautosomal region (p-PAR) is at least twenty times more frequent than the genomic average of approximately 1 cM/Mb, which may affect substitution patterns and rates in this region. Here I report the analysis of substitution patterns and rates in 10 human, chimpanzee, gorilla, and orangutan genes across the p-PAR. Between species silent divergence in the p-PAR forms a gradient, increasing toward the telomere. The correlation of silent divergence with distance from the p-PAR boundary is highly significant (rho = 0.911, P < 0.001). After exclusion of the CpG dinucleotides this correlation is still significant (rho = 0.89, P < 0.01), thus the substitution rate gradient cannot be explained solely by the differences in the extent of methylation across the p-PAR. Frequent recombination in the PAR may result in a relatively strong effect of biased gene conversion (BGC), which, because of the increased probability of fixation of the G or C nucleotides at (A or T)/(G or C) segregating sites, may affect substitution rates. BGC, however, does not seem to be the factor creating the substitution rate gradient in the p-PAR, because the only gradient is still detactable if only A<-->T and G<-->C substitutions are taken into account (rho = 0.82, P < 0.01). I hypothesize that the substitution rate gradient in the p-PAR is due to the mutagenic effect of recombination, which is very frequent in the distal human p-PAR and might be lower near the p-PAR boundary.     

The pseudoautosomal regions of the human sex chromosomes

In human females, both X chromosomes are equivalent in size and genetic content, and pairing and recombination can theoretically occur anywhere along their entire length. In human males, however, only small regions of sequence identity exist between the sex chromosomes. Recombination and genetic exchange is restricted to these regions of identity, which cover 2.6 and 0.4 Mbp, respectively, and are located at the tips of the short and the long arm of the X and Y chromosome. The unique biology of these regions has attracted considerable interest, and complete long-range restriction maps as well as comprehensive physical maps of overlapping YAC clones are already available. A dense genetic linkage map has disclosed a high rate of recombination at the short arm telomere. A consequence of the obligatory recombination within the pseudoautosomal region is that genes show only partial sex linkage. Pseudoautosomal genes are also predicted to escape X-inactivation, thus guaranteeing an equal dosage of expressed sequences between the X and Y chromosomes. Gene pairs that are active on the X and Y chromosomes are suggested as candidates for the phenotypes seen in numerical X chromosome disorders, such as Klinefelter's (47,XXY) and Turner's syndrome (45,X). Several new genes have been assigned to the Xp/Yp pseudoautosomal region. Potential associations with clinical disorders such as short stature, one of the Turner features, and psychiatric diseases are discussed. Genes in the Xq/Yq pseudoautosomal region have not been identified to date.     

Evidence for heterogeneity in recombination in the human pseudoautosomal region: high resolution analysis by sperm typing and radiation-hybrid mapping

Accurate genetic and physical maps for the human pseudoautosomal region were constructed by use of sperm typing and high-resolution radiation-hybrid mapping. PCR analysis of 1,912 sperm was done with a manual, single-sperm isolation method. Data on four donors show highly significant linkage heterogeneity among individuals. The most significant difference was observed in a marker interval located in the middle of the Xp/Yp pseudoautosomal region, where one donor showed a particularly high recombination fraction. Longitudinal models were fitted to the data to test whether linkage heterogeneity among donors was significant for multiple intervals across the region. The results indicated that increased recombination in particular individuals and regions is compensated for by reduced recombination in neighboring intervals. To investigate correspondence between physical and genetic distances within the region, we constructed a high-resolution radiation-hybrid map containing 29 markers. The recombination fraction per unit of physical distance varies between regions ranging from 13- to 70-fold greater than the genome-average rate.     

Significant positive correlation between the recombination rate and GC content in the human pseudoautosomal region.

This paper establishes that recombination drives the evolution of GC content in a significant way. Because the human P-arm pseudoautosomal region (PAR1) has been shown to have a high recombination rate, at least 20-fold more frequent than the genomic average of approximately 1 cM/Mb, this region provides an ideal system to study the role of recombination in the evolution of base composition. Nine non-coding regions of PAR1 are analyzed in this study. We have observed a highly significant positive correlation between the recombination rate and GC content (rho = 0.837, p < or = 0.005). Five regions that lie in the distal part of PAR1 are shown to be significantly higher than genomic average divergence. By comparing the intra- and inter-specific AT->GC -GC->AT ratios, we have detected no fixation bias toward GC alleles except for L254915, which has excessive AT-->GC changes in the human lineage. Thus, we conclude that the high GC content of the PAR1 genes better fits the biased gene conversion (BGC) model.     

Studying the biological and technical sources of variation in telomere length of individual chromosomes.

BACKGROUND: Consistent average length differences between species and chromosome arm differences within species indicate that telomere length is genetically determined. This seems to contradict an observed large variation in lengths of the same human telomere between metaphases of the same individual. We examined the extent to which the variation in the telomeres of the human X and Y chromosomes is heritable, induced, or technical in origin. METHODS: Metaphase chromosomes were stained by fluorescence in situ hybridization with a telomere repeat-specific probe, and fluorescence intensities of the X and Y chromosomes were measured. If telomere length variation is predominantly genetically determined and a 50% probability of meiotic recombination between the pseudo-autosomal regions of Yp and Xp in the father is taken into account, one expects an equal chance that the Yp telomere of a son is derived from his father's Xp or Yp telomere. This implies that the Yp/Yq telomere ratios in fathers and sons will be identical in the absence of paternal meiotic recombination and different when recombination occurs. RESULTS: Among five father-son pairs, four showed similar Yp/Yq ratios (P > 0.05), whereas one pair exhibited a large difference in the Yp/Yq ratio that was attributable to a significantly longer Xp than Yp telomere in the father and a presumptive meiotic exchange between X and Y during paternal meiosis. Further, the Xq telomere exhibited a generally shorter telomere length than the others. CONCLUSIONS: The high variation in telomere length appeared to be intracellular (between sister chromatids) and, hence, technical in nature. We found no measurable induced variation in the cells studied, implying that, if induced variation exists, it is small compared with the technical variation.     

High levels of sequence polymorphism and linkage disequilibrium at the telomere of 12q: implications for telomere biology and human evolution

The human Xp/Yp telomere-junction region exhibits high levels of sequence polymorphism and linkage disequilibrium. To determine whether this is a general feature of human telomeres, we have undertaken sequence analysis at the 12q telomere and have extended the analysis at Xp/Yp. A total of 22 single-nucleotide polymorphisms (SNPs) and one 30-bp duplication were detected in the 1,870 bp adjacent to the 12q telomere. Twenty polymorphic positions were in almost complete linkage disequilibrium, creating three common diverged haplotypes accounting for 80% of 12q telomeres in the white population. A further 6% of 12q telomeres contained a 1,439-bp deletion in the DNA flanking the telomere. The remaining 13% of 12q telomeres did not amplify with the primers used (nulls). The distribution of telomere (TTAGGG) and variant repeats within 12q telomeres was hypervariable, but alleles with similar distribution patterns were associated with the same haplotype in the telomere-adjacent DNA. These data suggest that 12q telomeres, like Xp/Yp telomeres, exhibit low levels of homologous recombination and evolve along haploid lineages. In contrast, high levels of homologous recombination occur in the adjacent proterminal regions of human chromosomes. This suggests that there is a localized telomere-mediated suppression of recombination. In addition, the genetic characteristics of these regions may provide a source of deep lineages for the study of early human evolution, unaffected by both natural selection and recombination. To explain the presence of a few diverged haplotypes adjacent to the Xp/Yp and 12q telomeres, we propose a model that involves the hybridization of two archaic hominoid lineages ultimately giving rise to modern Homo sapiens.     

DNA variability and recombination rates at X-linked loci in humans.

We sequenced 11,365 bp from introns of seven X-linked genes in 10 humans, one chimpanzee, and one orangutan to (i) provide an average estimate of nucleotide diversity (pi) in humans, (ii) investigate whether there is variation in pi among loci, (iii) compare ratios of polymorphism to divergence among loci, and (iv) provide a preliminary test of the hypothesis that heterozygosity is positively correlated with the local rate of recombination. The average value for pi was low 0.063%, SE = 0.036%, about one order of magnitude smaller than for Drosophila melanogaster, the species for which the best data are available. Among loci, pi varied by over one order of magnitude. Statistical tests of neutrality based on ratios of polymorphism to divergence or based on the frequency spectrum of variation within humans failed to reject a neutral, equilibrium model. However, there was a positive correlation between heterozygosity and rate of recombination, suggesting that the joint effects of selection and linkage are important in shaping patterns of nucleotide variation in humans.     

Rapid evolution of human pseudoautosomal genes and their mouse homologs

Comparative studies of genes in the pseudoautosomal region (PAR) of human and mouse sex chromosomes have thus far been very limited. The only comparisons that can presently be made indicate that the PARs of humans and mice are not identical in terms of gene content. Here we describe additional comparative studies of human pseudoautosomal genes and their mouse homologs. Using a somatic cell hybrid mapping panel, we have assigned the mouse homolog of the human pseudoautosomal interleukin 3 receptor alpha subunit (IL3RA) gene to mouse Chromosome (Chr) 14. Attempts to clone the mouse homolog of the human pseudoautosomal adenine nucleotide translocase-3 (ANT3) gene resulted in the isolation of the murine homologs of the human ANT1 and ANT2 genes. The mouse Ant1 and Ant2 genes are very similar in sequence to their human homologs, and we have mapped them to mouse Chromosomes (Chrs) (8 and X respectively) that exhibit conserved synteny with the chromosomes on which the human genes are located. In contrast, the homolog of ANT3 appears to be either very divergent or absent from the mouse genome. Southern blot analysis of DNA from a variety of mammalian species shows restricted conservation of human pseudoautosomal genes, a trend that also applies to the two cloned mouse homologs of these genes and to neighboring human genes in distal Xp22.3. Our observations combined with those of other workers lead us to propose a model for the evolution of the PAR that includes both rapid sequence evolution and the incremental reduction in size of the region during mammalian evolution. Received: 4 May 1995 / Accepted: 21 August 1995     

Elevated DNA sequence diversity in the genomic region of the phosphatase PPP2R3L gene in the human pseudoautosomal region

The evolution, inheritance and recombination rate of genes located in the pseudoautosomal region 1 (PAR1) is exceptional within the human genome. Pseudoautosomal genes are identical on X and Y chromosomes and are not inherited in a sex linked manner. Due to an obligatory recombination event in male meiosis, pseudoautosomal genes are exchanged frequently between X and Y chromosomes. During the isolation, characterization and sequencing of a novel gene PPP2R3L, which was classified by sequence homology as a novel member of the protein phosphatase regulatory subunit families, it became apparent that cosmids of different origin harboring this gene are highly polymorphic between individuals, both at the nucleotide level and in the number.     

Mapping human X-linked genes in the phalangerid marsupial Trichosurus vulpecula.

We mapped 15 human X-chromosome markers in the common brush-tailed possum, Trichosurus vulpecula (Kerr), which represents the Australian marsupial family Phalangeridae. In situ hybridization was used to localize highly conserved human X-linked genes to chromosomes of T. vulpecula diploid lines. Ten genes located on the long arm of the human X (human Xq genes) all mapped to the possum X chromosome. However, all five genes located on the short arm of the human X (human Xp genes) mapped to autosomes. These findings confirm our previous work, which showed that the X chromosome in macropodid and dasyurid marsupials bears all the human Xq genes but none of the human Xp genes studied. This suggests that the marsupial X is highly conserved, but its gene content reflects that of only part of the eutherian X, a result consistent with our hypothesis that an autosomal region was added to the X early in eutherian divergence.     

Comparative mapping of ZFY in the hominoid apes

Summary Within our project of comparative mapping of candidate genes for sex-determination/testis differentiation, we used a cloned probe from the human ZFY locus for comparative hybridization studies in hominoids. As in the human, the ZFY probe detects X- and Y-specific restriction fragments in the chimpanzee, the gorilla, the orangutan, and the gibbon. Furthermore, the X-specific hybridization site in the great apes resides in Xp21.3, the same locus defining ZFX in the human. The Y-specific locus of ZFY maps closely to the early replicating pseudoautosomal segment in the telomeric or subtelomeric position of the Y chromosomes of the great apes, again as found in the human. Thus, despite cytogenetically visible structural alterations within the euchromatic parts of the Y chromosomes of the human species and the great apes, a segment of the Y chromosome defined by the pseudoautosomal region and ZFY seems to be more strongly conserved than the rest of the Y chromosome.     

Mechanisms underlying telomere repeat turnover, revealed by hypervariable variant repeat distribution patterns in the human Xp/Yp telomere.

Sequences immediately adjacent to the human Xp/Yp telomere exhibit a high frequency of base substitutional polymorphisms, together with almost complete linkage disequilibrium, to create only a few diverged haplotypes. This sequence divergence has been used to develop a PCR-based system for mapping the distribution of the telomere (TTAGGG) and variant repeats (TGAGGG and TCAGGG) at the proximal end of the telomere repeat array. The distribution of these repeats is extremely variable. Almost all Xp/Yp telomeres are different, indicating a high mutation rate. Some telomere maps associated with the same flanking haplotype show similarities, identifying subsets of telomeres that share a recent common ancestry. Mechanisms underlying the rapid turnover of repeats at the proximal end of the Xp/Yp telomere include intra-allelic processes, such as slippage during replication. Inter-allelic exchanges may occur occasionally, but telomerase activity probably plays only a minor role in the germline turnover of proximally located telomere and variant repeats.     

The minimal mammalian Y chromosome - the marsupial Y as a model system

The mammalian X and Y chromosomes are very different in size and gene content. The Y chromosome is much smaller than the X and consists largely of highly repeated non-coding DNA, containing few active genes. The 65-Mb human Y is homologous to the X over two small pseudoautosomal regions which together contain 13 active genes. The heterochromatic distal half of the human Yq is entirely composed of highly repeated non-coding DNA, and even the euchromatic portion of the differential region is largely composed of non-coding repeated sequences, amongst which about 30 active genes are located. The basic marsupial Y chromosome (about 10 Mb) is much smaller than that of humans or other eutherian mammals. It appears to include no PAR, since it does not undergo homologous pairing, synaptonemal complex formation or recombination with the X. We show here that the tiny dunnart Y chromosome does not share cytogenetically detectable sequences with any other chromosome, suggesting that it contains many fewer repetitive DNA sequences than the human or mouse Y chromosomes. However, it shares several genes with the human and/or mouse Y chromosome, including the sex determining gene SRY and the candidate spermatogenesis gene RBMY, implying that the marsupial and eutherian Y are monophyletic. This minimal mammalian Y chromosome might provide a good model Y in which to hunt for new mammalian Y specific genes.     

How strong is the mutagenicity of recombination in mammals?

It is commonly believed that a high recombination rate such as that in a pseudoautosomal region (PAR) greatly increases the mutation rate because a 170-fold increase was estimated for the mouse PAR region. However, sequencing PAR and non-PAR introns of the Fxy gene in four Mus taxa, we found an increase of only twofold to fivefold. Furthermore, analyses of sequence data from human and orangutan PAR and X-linked regions and from autosomal regions showed a weak effect of recombination on mutation rate (a slope of less than 0.2% per cM/Mb), although a much stronger effect on GC content (1% to 2% per cM/Mb). Because typical recombination rates in mammals are much lower than those in PARs, the mutagenicity of recombination is weak or, at best, moderate, although its effect on GC% is much stronger. In addition, contrary to a previous study, we found no Fxy duplicate in Mus spretus.     

Regional assignment of Y-linked DNA probes by deletion mapping and their homology with X-chromosome and autosomal sequences.

A series of Y recombinants have been isolated from Y-specific DNA libraries and regionally located on the Y chromosome using a Y deletion panel constructed from individuals carrying structural abnormalities of the Y chromosome. Of twenty recombinants examined twelve have been assigned to Yp and eight to Yq. Five of the Yp recombinants map between Yp11.2 and Ypter and one can only be assigned to Yp. Of the former, four detect homologies on the X chromosome between Xq13 and Xq24 and the latter one between Xp22.3 and Xpter. The sixth recombinant detects autosomal homologous sequences. The six remaining Yp probes are located between Ycen and Yp11.2. One of these detects a homology on the X chromosome at Xq13-Xq24 and a series of autosomal sequences, two detect uniquely Y-specific sequences and three a complex pattern of autosomal homologies. The remaining eight recombinants have been assigned to three intervals on Yq. Of three recombinants located between Ycen and Yq11.21 two detect only Y sequences and one additional autosomal homologies. Two recombinants lie in the interval Yq11.21-Yq11-22, one of which detects only Y sequences and the other an Xp homology between Xp22.3 and Xpter. Finally, the three remaining Yq recombinants all detect autosomal homologies and are located between Yq11.22 and Yq12. The divergence between homologies on different chromosomes has been examined for three recombinants by washing Southern Blots at different levels of stringency. Additionally, Southern analysis of DNA from flow sorted chromosomes has been used to identify autosomes carrying homologies to two of the Y recombinants.     

Isolation of sequences from Xp22.3 and deletion mapping using sex chromosome rearrangements from human X-Y interchange sex reversals

A repeated DNA element (STIR) interspersed in Xp22.3 and on the Y chromosome has been used as a tag to isolate seven single-copy probes from the human sex chromosomes. The seven probes detect X-specific loci located in Xp22.3. Using a panel of X-chromosomal deletions from X-Y interchange sex reversals (XX males and XY females), these X-specific loci and some additional ones were mapped to four contiguous intervals of Xp22.3, proximal to the pseudoautosomal region and distal to STS. The construction of this deletion map of the terminal part of the human X chromosome can serve as a starting point for a long-range physical map of Xp22.3 and for a more accurate mapping of genetic diseases located in Xp22.3.     

Evolutionary rate of a gene affected by chromosomal position.

Genes evolve at different rates depending on the strength of selective pressure to maintain their function. Chromosomal position can also have an influence [1] [2]. The pseudoautosomal region (PAR) of mammalian sex chromosomes is a small region of sequence identity that is the site of an obligatory pairing and recombination event between the X and Y chromosomes during male meiosis [3] [4] [5] [6]. During female meiosis, X chromosomes can pair and recombine along their entire length. Recombination in the PAR is therefore approximately 10 times greater in male meiosis compared with female meiosis [4] [5] [6]. The gene Fxy (also known as MID1 [7]) spans the pseudoautosomal boundary (PAB) in the laboratory mouse (Mus musculus domesticus, C57BL/6) such that the 5' three exons of the gene are located on the X chromosome but the seven exons encoding the carboxy-terminal two-thirds of the protein are located within the PAR and are therefore present on both the X and Y chromosomes [8]. In humans [7] [9], the rat, and the wild mouse species Mus spretus, the gene is entirely X-unique. Here, we report that the rate of sequence divergence of the 3' end of the Fxy gene is much higher (estimated at 170-fold higher for synonymous sites) when pseudoautosomal (present on both the X and Y chromosomes) than when X-unique. Thus, chromosomal position can directly affect the rate of evolution of a gene. This finding also provides support for the suggestion that regions of the genome with a high recombination frequency, such as the PAR, may have an intrinsically elevated rate of sequence divergence.     

Recombination, GC-content and the human pseudoautosomal boundary paradox

The pseudoautosomal boundary of mammalian sex chromosomes separates a low-recombination region (X- or Y-specific) from a high-recombination region (the pseudoautosomal region), providing a good opportunity to investigate the influence of recombination on molecular evolutionary processes. The mouse and human patterns of sequence variation, however, are discordant: a striking difference of GC-content and evolutionary rate was reported between the proximal and distal sides of the pseudoautosomal boundary in the mouse genome, whereas this difference was not found in the human genome. The paradox might be explained by the mirror histories of the pseudoautosomal boundary in the two species, and by the asymmetric nature of the forces driving GC-content evolution in mammalian genomes.