# Recombination Dynamics of a Human Y-Chromosomal Palindrome: Rapid GC-Biased Gene Conversion, Multi-kilobase Conversion Tracts, and Rare Inversions

Pille Hallast, Patricia Balaresque, Georgina R Bowden, Stéphane Ballereau, Mark A Jobling

<figure><img src="image_1.png" /></figure>

To cite this version:

Pille Hallast, Patricia Balaresque, Georgina R Bowden, Stéphane Ballereau, Mark A Jobling. Recombination Dynamics of a Human Y-Chromosomal Palindrome: Rapid GC-Biased Gene Conversion, Multi-kilobase Conversion Tracts, and Rare Inversions. PLoS Genetics, 2013, 9 (7), pp.e1003666. 10.1371/journal.pgen.1003666. [https://hal.science/hal-05203136v1](https://hal.science/hal-05203136v1)

HAL Id: hal-05203136

Submitted on 7 Aug 2025

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Recombination Dynamics of a Human Y-Chromosomal Palindrome: Rapid GC-Biased Gene Conversion, Multi-kilobase Conversion Tracts, and Rare Inversions

Pille Hallast, Patricia Balaresque<sup>a,a</sup>, Georgina R. Bowden<sup>a,b</sup>, Stéphane Ballereau<sup>a,c</sup>, Mark A. Jobling<sup>*</sup>

Department of Genetics, University of Leicester, Leicester, United Kingdom

## Abstract

The male-specific region of the human Y chromosome (MSY) includes eight large inverted repeats (palindromes) in which arm-to-arm similarity exceeds 99.9%, due to gene conversion activity. Here, we studied one of these palindromes, P6, in order to illuminate the dynamics of the gene conversion process. We genotyped ten paralogous sequence variants (PSVs) within the arms of P6 in 378 Y chromosomes whose evolutionary relationships within the SNP-defined Y phylogeny are known. This allowed the identification of 146 historical gene conversion events involving individual PSVs, occurring at a rate of 2.9–8.4×10<sup>-4</sup> events per generation. A consideration of the nature of nucleotide change and the ancestral state of each PSV showed that the conversion process was significantly biased towards the fixation of G or C nucleotides (GC-biased), and also towards the ancestral state. Determination of haplotypes by long-PCR allowed likely co-conversion of PSVs to be identified, and suggested that conversion tract lengths are large, with a mean of 2068 bp, and a maximum in excess of 9 kb. Despite the frequent formation of recombination intermediates implied by the rapid observed gene conversion activity, resolution via crossover is rare: only three inversions within P6 were detected in the sample. An analysis of chimpanzee and gorilla P6 orthologs showed that the ancestral state bias has existed in all three species, and comparison of human and chimpanzee sequences with the gorilla outgroup confirmed that GC bias of the conversion process has apparently been active in both the human and chimpanzee lineages.

**Citation:** Hallast P, Balaresque P, Bowden GR, Ballereau S, Jobling MA (2013) Recombination Dynamics of a Human Y-Chromosomal Palindrome: Rapid GC-Biased Gene Conversion, Multi-kilobase Conversion Tracts, and Rare Inversions. PLoS Genet 9(7): e1003666. doi:10.1371/journal.pgen.1003666

**Editor:** Mikkel H. Schierup, Aarhus University, Denmark

**Received** February 7, 2013; **Accepted** June 7, 2013; **Published** July 25, 2013

**Copyright:** © 2013 Hallast et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** MAJ was supported by a Wellcome Trust Senior Fellowship in Basic Biomedical Science (grant no. 087576); PH, PB, GRB and SB were supported by the Wellcome Trust. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: maj4@le.ac.uk

## Introduction

The male-specific region of the human Y chromosome (MSY) is constitutively haploid, yet contains a high proportion (~35%) of pseudo-diploid duplicated regions, eight of which are arranged as large inverted repeats ('palindromes', known as P1 - P8; Figure 1a), with arms in most cases separated by non-duplicated spacers [1]. The arms of each palindrome are >99.9% similar in sequence due to the homogenising effect of gene conversion. Human-chimpanzee sequence divergence within palindrome arms is significantly lower than that within spacers, and compared to the MSY (non-palindrome, non-spacer) average [2]. This suggests that gene conversion since speciation may have been directional, tending to return new mutations that arise within arms to their ancestral states. Most palindrome arms are enriched in testis-specific genes, important in spermatogenesis, and the suggestion has been made that directional gene conversion between pseudo-diploid copies may protect these genes against evolutionary decay [2].

It is becoming increasingly recognised that such palindromic structures are far from being a peculiarity of great ape Y chromosomes, but have more general biological significance as a feature of independently arising constitutively haploid sex chro-

mosomes in other mammals [3–5], birds [6,7] and insects [8], as well as of the mammalian X chromosome [9,10], which is haploid in males. Yet despite this general importance, and despite some theoretical analyses of palindrome evolution [11,12], little is known about the dynamics of conversion within these remarkable structures.

Large stretches of sequence identity between palindrome arms represent compelling evidence for rapid gene conversion, and yet, paradoxically, provide a barrier to understanding the dynamics of the conversion process. Conversion rate, tract length, and directionality cannot be examined when there are no sequence differences (paralogous sequence variants; PSVs) between arms that might allow specific conversion events to be recognised. However, when a PSV does exist (e.g. the 'pseudoheterozygous' state G/A), then the observation in other chromosomes of the two other possible genotypes, the 'pseudohomozygous' G/G and A/A, indicates that gene conversion must have occurred within the history of the examined sequences (Figure 2a), assuming that recurrent mutation can be neglected. Such an observation tells us nothing about how many independent conversion events underlie the three genotypes. But the availability of a detailed and robust Y phylogeny, defined by stable single nucleotide polymorphisms

<footnote><sup>a,a</sup> Current address: Molecular Anthropology and Image Synthesis Laboratory, Université Paul Sabatier, Toulouse, France.</footnote>

<footnote><sup>b</sup> Current address: Illumina, Chesterford Research Park, Little Chesterford, Nr Saffron Walden, Essex, United Kingdom.</footnote>

<footnote><sup>c,c</sup> Current address: ENS-CNRS USR3010, Campus Charles Mérieux - Université Claude Bernard Lyon 1, Lyon, France.</footnote>
## Author Summary

The sex-determining role of the human Y chromosome makes it male-specific, and always present in only a single copy. This solo lifestyle has endowed it with some bizarre features, among which are eight large DNA units constituting about a quarter of the chromosome's length, and containing many genes important for sperm production. These units are called palindromes, since, taking into account the polarity of the DNA strands, the sequence is the same read from either end of the unit. We investigated the details of a process (gene conversion) that transfers sequence variants in one half of a palindrome into the other, thereby maintaining >99.9% similarity between the halves. We analysed patterns of sequence variants within one palindrome in a set of Y chromosomes whose evolutionary relationships are known. This allowed us to identify past gene conversion events, and to demonstrate a bias towards events that eliminate new variants, and retain old ones. Gene conversion has therefore acted during human evolution to retard sequence change in these regions. Analysis of the chimpanzee and gorilla versions of the palindrome shows that the dynamic processes we see in human Y chromosomes have a deep evolutionary history.

(SNPs) outside the palindromic regions, allows the evolutionary relationships of palindrome sequences to be known, and genotyping within this phylogenetic context can then provide an estimate of the minimum number of conversion events (Figure 2b–d).

Genotyping PSVs within a phylogenetic context provides evidence for past gene conversion events, but the resulting genotypes (pursuing the analogy of diploidy) are ‘unphased’ – we do not know which allele of a PSV lies on which palindrome arm. Because of the high degree of sequence identity and the scarcity of PSVs within palindromes, phasing is technically challenging, but nonetheless important if we are to gain an understanding of the lengths of conversion tracts, suggested by sets of co-converted adjacent PSVs (Figure 2e). If phased PSV data for palindromes were available, it would also be possible to address an additional important aspect of the dynamics of these structures: the ratio of

non-reciprocal exchanges (conversions) to reciprocal exchanges (inversions; Figure 2f).

Here, we analyse paralogous sequence variants (PSVs) within the arms of human palindrome P6, taking the approaches outlined above. We demonstrate through a phylogenetic analysis of conversion events five cardinal features of the palindrome conversion process during human evolution: (i) the conversion process has been rapid throughout the evolution of modern human Y-chromosomal lineages; (ii) it shows significant bias to the fixation of GC base pairs; (iii) it is biased towards the retention of ancestral states of PSVs; (iv) conversion tracts can encompass several kilobases; and (v) despite the high frequency of recombination events within palindrome arms, these resolve overwhelmingly via non-reciprocal exchange (conversions) rather than reciprocal exchange (inversions). We then extend our findings to deeper evolutionary time by determining the sequence of gorilla P6, showing that ancestral state bias has existed in the gorilla lineage as well as in humans and chimpanzees, and allowing us to ascertain the direction of evolutionary changes in the human and chimpanzee lineages, revealing a possible long-term GC conversion bias.

## Results

### Palindrome P6 as a study system

To study gene conversion dynamics we first sought a segment of a palindrome carrying a suitable number and density of PSVs. Arm-to-arm alignment of the reference sequence (belonging to haplogroup R1b1b2* [13]) for palindrome P6 (Figure 1b) demonstrated a 99.97% sequence similarity between its 110-kb arms, but revealed a total of 49 discrete sequence differences, which we supplemented with two additional single-nucleotide PSVs identified from the sequencing of a flow-sorted Y chromosome from a different source, belonging to haplogroup O3a [14]. Twenty-nine of these represent simple single-nucleotide PSVs (SN-PSVs) that are unlikely to undergo mutational reversion or recurrence. Furthermore, 16 SN-PSVs lie within 20 kb of the outer arm boundaries, potentially allowing arm-specific PCR anchored in flanking single-copy DNA to determine in which arm a particular variant lies (‘phasing’). Two additional factors favour P6: chimpanzee and gorilla orthologs exist that allow the ancestral

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location and structure of palindrome P6, showing positions of PSVs analysed. a) Idiogram of Y chromosome, showing positions of the 8 palindromes, with structure and coordinates (in GRCh37) of P6 below. b) Position and nature of the differences between the arms of P6, indicating the 10 SN-PSVs analysed, and the positions of PCR primers used in arm-specific amplifications. STSs marking the arm boundaries are also shown (with ‘sY’ prefixes). Asterisks indicate the two SN-PSVs identified from a haplogroup O3a chromosome.
doi:10.1371/journal.pgen.1003666.g001</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Recognition of gene conversion, co-conversion and inversion events. a) Existence of three genotypes at a hypothetical PSV indicates that gene conversion has taken place, if recurrent mutation is neglected. Genotyping the PSV in a phylogenetic context, and applying the principle of maximum parsimony, allows the recognition of: b) Haplogroup descending from an ancestor in which the PSV mutation has not yet arisen (G/G), and is therefore uninformative; c) Haplogroup descending from an ancestor in which the PSV mutation has arisen (G/A), but shows no variation, and therefore no evidence for gene conversion; d) Haplogroup descending from an ancestor in which the PSV mutation has arisen, and shows evidence of at least two bidirectional conversion events (G/G and A/A); e) Recognition of co-conversion of more than one PSV requires 'phase' information, as does f) recognition of inversions.
doi:10.1371/journal.pgen.1003666.g002</figcaption></figure>

state of its PSVs to be determined; and P6 lacks protein-coding genes, meaning that direct effects of natural selection are less likely than for other palindromes.

We sought to design reliable typing assays for all SN-PSVs, and this was successfully accomplished (see Materials & Methods) for ten, indicated in Figure 1b.

In order to identify gene conversion events between the arms of P6, we required a set of Y chromosomes for which detailed phylogenetic relationships were well established. We exploited the availability of the CEPH-Human Genome Diversity Project (HGDP) panel of DNA samples [15], which has good global coverage and for which data were available for 184 Y-
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Phylogenetic analysis of gene conversion events within P6. Adjacent to the Y phylogeny based on binary markers is shown a schematic representation of allelic states of nine informative PSVs (PSV10 is omitted because it is invariant). Note that the number of haplotypes shown within each haplogroup is less than the number of samples genotyped. Circles represent PSVs with inter-PSV distances not to scale, and the colours of circles indicate uninformative, pseudoheterozygous and pseudohomozygous PSVs, as shown in the key.
doi:10.1371/journal.pgen.1003666.g003</figcaption></figure>

chromosomal binary markers ([16–18]; www.cephb.fr/en/hgdp/),
supplementing this by typing an additional 23 binary markers, to
define a total of 63 haplogroups. The tree thus defined, the
markers, the haplogroup nomenclature and the sources of data are
shown in Figure S1. A simplified version of the tree is shown in
Figure 3.

PSV genotypes and gene conversion bias within P6

The ten PSVs were analysed in a subset of 378 of the 684
HGDP male samples, chosen to cover the haplogroup diversity of
the sample set. Each PSV genotype was recorded as pseudoheter-
ozygous (e.g. G/A) or pseudohomozygous (e.g. G/G or A/A), and,
by comparison to the orthologous sequences in chimpanzee [19]
and gorilla ([20], and our own gorilla sequence – see below) each
PSV allele coded as ancestral (0) or derived (1).

Figure 3 illustrates the patterns of variation observed in the
sample, and full details are given in Table S1. Some PSVs (e.g.
PSV6) are variable across all haplogroups, suggesting that the
variant arose at the root of the Y phylogeny. Others show
variability only in specific haplogroups (hg), suggesting (assuming
maximum parsimony, and no recurrent mutation) that they arose
in their founders (e.g. PSV2 in hgF, PSV5 in hgP, and PSV9 in
hgO3a). PSV10 was monomorphic in all 378 cases tested,
suggesting that it represents a recently arising variant. For any
haplogroup, we can deduce whether the founder was pseudohe-
terozygous (0/1); when this is so, the finding of pseudohomozygous
states (0/0 or 1/1) among chromosomes within the haplogroup
indicates that conversion must have occurred (Figure 2d). Treating
each PSV as an independent site of gene conversion, it is thus
possible to both count the total number of conversion events, and
to ask what proportion of these are conversions to the ancestral state (i.e. 0/1 to 0/0), or the derived state (i.e. 0/1 to 1/1).

This analysis identified a total of 146 converted SN-PSVs, of which 86 represent conversion to the ancestral, and 60 to the derived state (Tables S1 and S2). This difference is statistically significant (p = 0.0314; Chi-square test), which is consistent with published observations based on human-chimpanzee comparisons [2].

We can also ask if there is a tendency towards the fixation of GC base-pairs rather than AT base-pairs: this is so-called GC-biased gene conversion, and results from a bias in the repair of AC and GT mismatches that form in heteroduplex recombination intermediates [21]. Of the 146 converted SN-PSVs, some are uninformative because they involve transversions (from CG to GC, or AT to TA), but among the 79 informative cases 59 involve the fixation of GC, and 20 of AT (p = 1.1×10⁻⁵; Chi-square test). From these observations, gene conversion among the studied PSVs appears to be strongly GC-biased, and slightly but significantly biased towards the retention of ancestral state.

Having counted the number of observed gene conversion events in our dataset, we can estimate an average rate of gene conversion by dividing by the number of generations encompassed in the phylogeny that relates the studied Y chromosomes (Materials & Methods). For a 25-year generation time, this yields a rate of 2.9–8.4×10⁻⁴ events per generation.

## Infrequent inversions within P6

The above analysis provides evidence of a highly active gene conversion process within P6: but does the frequent formation of recombination intermediates that this implies also give rise to frequent inversions of the palindrome arms? As explained in the Introduction (Figure 2f), identification of inversion events requires the palindrome arms to be ‘phased’ at pseudoheterozygous sites. In order to do this, an arm-specific long-range PCR approach was developed, using one universal primer binding within the arm, and another binding to a distal-arm-specific region outside the outer palindrome boundary. This generated a product of ~18.9 kb incorporating seven of the studied PSVs (PSV1–7) that could then be typed in an arm-specific manner, thus determining their phase.

Arm-specific haplotypes from 83 selected DNA samples representing all of the haplogroups were compared to the Y-chromosome reference sequence, whose phase is known from BAC clone sequencing [1]. All but five samples were found to have identical phase to the reference sequence (Figure 4) at informative (pseudoheterozygous) sites; this corresponds to just three independent inversion events, in haplogroups A3b2*, B2a, and D2. Where phase information is available for several chromosomes within a haplogroup, these are always concordant – in other words, inversions are rare. This strong preponderance of conversion over inversion allows us to infer the phase of unphased chromosomes within haplogroups.

Among the 83 phased chromosomes, the three inversion events compare with 56 gene conversion events (assuming each converted PSV represents a single event). In the same set of chromosomes, and under the same assumptions, the per-generation rate of inversion is 1.36–1.72×10⁻⁵, compared to 2.54–3.21×10⁻⁴ for conversion. The latter rate differs from that given in the section above due to the smaller number of chromosomes phased and analysed here.

## Co-conversion of PSVs

All of the analysis above assumes that PSVs are independently converted, but from simple inspection of the behaviour of the adjacent PSVs 1 and 2, separated by only 81 bp, it is evident that

### Reference sequence

<figure><img src="image_6.png" /></figure>

### Haplogroups showing inversions

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Schematic representation of the inversions identified within three haplogroups. The approximate position of each inversion is indicated by a red cross.
doi:10.1371/journal.pgen.1003666.g004</figcaption></figure>

co-conversions must be occurring: for example, of the 14 instances where conversion affecting PSV1 and PSV2 is informative, 11 involve apparent co-conversion of both markers (e.g. in hgM1*; Figure 3). We therefore wished to examine co-conversion more systematically, and the phasing information allows us to do this (as shown in Figure 2e). The true number of co-conversion events is impossible to estimate, because the apparent co-conversion of adjacent variants could actually reflect the sum of two independent events. However, we can estimate the minimum number of co-conversions that explain the observed data: first, we identify adjacent pairs of pseudohomozygous PSVs within a haplogroup whose founder is deduced to be pseudoheterozygous for the same PSVs; and second, to exclude independent opposite conversions as an explanation (Figure 2e), we count only those PSV pairs that match a single arm-specific haplotype of the reference sequence. We then assume that these reflect a single conversion tract. On this
<table><caption>**Table 1.** Inter-specific sequence divergence in arms and spacers of palindrome P6.</caption><thead><tr><th>Species comparison</th><th>Region of P6</th><th>Ungapped length (bp)</th><th>No of nt substitutions</th><th>Divergence (%)</th><th>P-value (arm vs spacer)<sup>a</sup></th></tr></thead><tbody><tr><td rowspan="2">Human vs chimpanzee</td><td>Arm</td><td>104230</td><td>1497</td><td>1.44</td><td rowspan="2">1.52×10<sup>-15</sup></td></tr><tr><td>Spacer</td><td>45959</td><td>919</td><td>2.00</td></tr><tr><td rowspan="2">Human vs gorilla</td><td>Arm</td><td>88031</td><td>1726</td><td>1.96</td><td rowspan="2">5.04×10<sup>-8</sup></td></tr><tr><td>Spacer</td><td>31206</td><td>773</td><td>2.48</td></tr><tr><td rowspan="2">Chimpanzee vs gorilla</td><td>Arm</td><td>84096</td><td>1820</td><td>2.16</td><td rowspan="2">6.04×10<sup>-7</sup></td></tr><tr><td>Spacer</td><td>31097</td><td>828</td><td>2.66</td></tr></tbody></table>

<sup>a</sup>2×2 contingency table, Chi-square test with Yates correction.
doi:10.1371/journal.pgen.1003666.t001

basis, 49 of the 107 (45.8%) individual conversion events among ‘phased’ PSVs (1–7) can be explained by a minimum of 18 co-conversion tracts.

We cannot arrive at a useful estimate of maximum co-conversion tract length, because most tracts are not flanked by informative genotyped markers that would indicate their outer limits. However, we can estimate minimum lengths by considering the distance between the outer converted markers within each tract. The mean value of these minimum estimates is 2068 bp: this is much longer than most recorded gene conversion events, which are typically a few hundred bp in length, and rarely exceed 1 kb [22].

Some apparent co-conversion tracts are very long indeed. For example, within hgQ1a* we observe PSVs 1–8 in the pseudoherozygous state, but also a case where the first seven of these variants are pseudohomozygous. This case seems unlikely to have arisen as a result of a series of consecutive small-scale conversion events, because the allelic state of the variants matches a single arm-specific haplotype in the same haplogroup. An alternative trivial explanation is that one arm in this chromosome has been lost by deletion, and that the PSVs are being observed in a pseudohemizygous, rather than pseudohomozygous state. To eliminate this possibility we confirmed that two arms were present, and were of the expected length, using two methods: a paralog-ratio test (PRT) to measure the copy-number of the palindrome arm with respect to a reference sequence on the X chromosome; and a long-PCR assay specific for each arm in turn. The most parsimonious explanation for the observed genotype in this chromosome is therefore a massive conversion event that spans at least 9023 bp (the distance between PSVs 1 and 7 on the proximal arm).

The analysis carried out above, to detect biases in gene conversion towards retention of the ancestral state and fixation of GC base-pairs, treated each converted nucleotide as an independent replicate in a statistical test. However, since we have inferred that co-conversions occur, some variants are not independent; we therefore repeated both tests after removing the putative co-conversion events. In both cases, the statistical significance of the bias is retained (Table S2).

## Deeper evolutionary history of conversion bias

In order to study the deeper history of gene conversion activity and its impact on palindrome evolution, we required an outgroup to human and chimpanzee P6. A high-quality MSY sequence is available for rhesus macaque that contains three palindromes, but

a P6 ortholog is not among them [3]. A gorilla Y-chromosome reference sequence is not yet available, but this species is known to carry both P6 arm-spacer boundaries with almost identical sequence to human and chimpanzee [2]. We constructed a partial sequence of gorilla P6 by merging Illumina paired-end sequencing data from two whole-genome-sequenced male gorillas [20] and from an independent male analysed in a sequence capture experiment. A total of 88,031 bp of merged gorilla P6 arm and 31,206 bp of gorilla spacer were assembled using the human Y-chromosome sequence as a reference. These data represent 80% of the human proximal arm and 67.5% of human spacer. The presence of both P6 arms in gorilla is confirmed by the fact that the mean coverage of proximal arm for all three gorillas is approximately twice that of the spacer (Protocol S1).

Pairwise alignments between human, chimpanzee and gorilla show that nucleotide divergence in all three comparisons is highly significantly reduced in the arms of P6 compared to spacer (Table 1). This is consistent with previous results [2] showing a similar pattern when comparing segments of Y-chromosome palindromes between human and chimpanzee. Our findings therefore confirm that the processes influencing palindrome evolution are active in both human and chimpanzee lineages, and also probably active in gorilla.

Availability of an outgroup sequence also allows possible long-term GC-bias to be examined in human and chimpanzee lineages. We used a phylogenetic approach to study nucleotide replacements in palindrome arms and spacer. Since the universally low (~0.02%) arm-to-arm divergence suggests that conversion is highly active within each species, all replacements found in arms can be assumed to be due to mutation followed by gene conversion; in spacers the divergence is expected to arise solely from mutational processes. From the alignment of human, chimp and gorilla P6 sequences (Dataset S1), we determined the types of all fixed differences, noting G or C (S) nucleotides that changed to A or T (W) nucleotides, and vice versa. We also determined the evolutionary direction of each of these fixed differences: if a nucleotide was identical between chimpanzee and gorilla but divergent in human, a replacement on the human lineage was assumed; if human and gorilla were identical, a replacement in chimpanzee was assumed.

Table 2 summarises the numbers and types of nucleotide replacements in both the human and chimpanzee arms and spacers. In the arm, the proportion of W to S changes slightly exceeds that in the spacer, but the proportion of S to W changes is significantly lower than that in the spacer. Furthermore, in human
<table><caption>**Table 2.** Patterns of P6 nucleotide replacements in the human and chimpanzee lineages.</caption><thead><tr><th></th><th></th><th>Total no. replacements</th><th>W to S changes/total AT nt (%)</th><th>Ratio of changes to W/to S<sup>a</sup></th><th>S to W changes/total GC nt (%)</th></tr></thead><tbody><tr><td rowspan="3">Before dinucleotide removal</td><td>Human</td><td>487</td><td>235/51201 (0.46)</td><td>1.07</td><td>161/32892 (0.49)</td></tr><tr><td>Spacer</td><td>238</td><td>72/19687 (0.37)</td><td>2.97**</td><td>124/11404 (1.09)</td></tr><tr><td>Arm: spacer ratio<sup>a</sup></td><td></td><td>1.25</td><td></td><td>0.45**</td></tr><tr><td rowspan="3">Chimpanzee</td><td>Human</td><td>659</td><td>366/51201 (0.71)</td><td>0.77**</td><td>181/32892 (0.55)</td></tr><tr><td>Spacer</td><td>297</td><td>110/19687 (0.56)</td><td>2.15**</td><td>137/11404 (1.20)</td></tr><tr><td>Arm: spacer ratio<sup>a</sup></td><td></td><td>1.28*</td><td></td><td>0.46**</td></tr><tr><td rowspan="3">After dinucleotide removal</td><td>Human</td><td>180</td><td>80/38071 (0.21)</td><td>1.83**</td><td>74/19215 (0.39)</td></tr><tr><td>Spacer</td><td>100</td><td>24/14935 (0.16)</td><td>5.28**</td><td>56/6596 (0.85)</td></tr><tr><td>Arm: spacer ratio<sup>a</sup></td><td></td><td>1.31</td><td></td><td>0.45**</td></tr><tr><td rowspan="3">Chimpanzee</td><td>Human</td><td>270</td><td>135/38071 (0.35)</td><td>1.29</td><td>88/19215 (0.46)</td></tr><tr><td>Spacer</td><td>134</td><td>43/14935 (0.29)</td><td>3.95**</td><td>75/6596 (1.14)</td></tr><tr><td>Arm: spacer ratio<sup>a</sup></td><td></td><td>1.23</td><td></td><td>0.40**</td></tr></tbody></table>

<sup>a</sup>2×2 contingency table, Chi-square test with Yates correction.

*p-value<0.05.

**p-value<0.01.

doi:10.1371/journal.pgen.1003666.t002

the proportion of W to S and S to W changes in the arm are approximately equal, while in the spacer S to W changes significantly predominate (as has been observed previously for substitutions not associated with gene conversion [23]). These observations indicate a relative bias towards W to S changes in arms. In chimpanzee P6, the proportion of W to S changes in the arm is significantly higher than that of S to W. In order to eliminate the potential influence of hypermutable CpG dinucleotides, all sites in CpG, TpG, or CpA sequences were removed from the raw sequence alignment, and the comparisons repeated. In both human and chimpanzee P6, the differences between arm and spacer remain. These striking differences in substitution patterns in arms and spacer seems likely to reflect the preferential fixation of GC base-pairs in arms due to the action of GC-biased gene conversion.

## Discussion

In this study we have used a phylogenetic approach to the diversity of sequences within a Y-chromosomal palindrome, P6, to illuminate the dynamic processes of recombination that distinguish these remarkable structures.

Analysis of a set of ten PSVs in 378 chromosomes has revealed 146 individual PSV conversion events in the Y phylogeny, and confirms that gene conversion is an ongoing and rapid process. Our findings add to the body of evidence showing that, despite its exemption from the otherwise ubiquitous process of meiotic crossing over, the MSY is highly active in gene conversion, involving not only palindromes [2,24], but also widely separated direct repeats [25] and minisatellite arrays [26]. As well as intrachromosomal conversions, gametologous regions on the Y have been shown to exchange sequences with the X chromosome in humans [27,28] as well as other mammals [29–32].

## Basic parameters of gene conversion

We observe a conversion rate of 2.9–8.4×10⁻⁴ events per generation among the 10 surveyed PSVs. This equates to a per-

PSV rate of 2.9–8.4×10⁻⁵ events per generation, though this represents a minimum estimate, since not all of the PSVs are informative in all studied chromosomes. Based on a measure of inter-arm divergence and an estimate of the base-substitution rate, Rozen et al. [2] estimated a conversion rate per generation per nucleotide in Y palindromes of 2.2×10⁻⁴.

Although gene conversion tracts several kilobases in length occur frequently in yeast [33], in mammals tracts are short, typically ranging from a few tens of base pairs [34] to 1 kb [22]. In palindrome P6, we infer minimum gene conversion tract lengths up to 9023 bp with mean minimum length of 2068 bp. These lengths do not represent direct measurements, and it remains possible that the inferred patterns of long conversion tracts could be created by multiple independent events. However, the longest inferred tract, including 7 PSVs, would require several independent events all in the same direction (from proximal to distal arm), so the most parsimonious explanation is a single event. It is possible that long conversion tracts are a typical characteristic of palindromes, but this remains to be tested by future studies.

## Resolution of recombination intermediates as inversions

Recombination is initiated by double-strand breaks (DSBs), and their repair can result in either reciprocal crossover, or non-reciprocal conversion. In considering the effects of these different pathways in P6, we need to differentiate between inter- and intramolecular events since, while conversion between or within chromatids will have the same molecular outcomes, this is not the case for crossover. Inter-chromatid crossover is expected to lead to an isodicentric chromosome and an acentric fragment, both of which are selected against. For example, 7/8 human palindromes are involved in crossover events between sister chromatids resulting in large-scale rearrangements in patients with disorders ranging from spermatogenic failure to sex reversal and Turner syndrome [35]. By contrast, intra-chromatid crossover will lead to simple inversion of palindrome arms, which seems unlikely to have strong effects on fitness. As an example, crossover between IR3 inverted repeats on Yp, resulting in apparently asymptomatic
inversion, has occurred at least twelve times in the history of the Y
phylogeny [36]. This different consequences of the two pathways
means that while observed conversions reflect both inter- and
intramolecular events, observed inversions are the result of intra-
chromatid events only, and this complicates the interpretation of
conversion: crossover ratios.

Phylogenetic detection of intra-chromatid crossovers leading to
palindrome inversions is possible if the phase of the PSVs is
known. Phasing of seven of the studied PSVs, located within the
first ~19 kb from the outer palindrome boundaries provides
evidence of only three independent inversions among the studied
chromosomes (Figure 4). The deduced rate of inversion, 1.36–
1.72×10⁻⁵ per generation, compares to a published rate of
2.3×10⁻⁴ for the IR3 inverted repeats [36]. Notably, we have
ascertained only those inversions with breakpoints occurring in the
outer ~16% of the arms of P6, whereas the published study was
able to ascertain all intra-chromatid inversions by determining the
orientation of markers between the IR3 repeats. Our finding of 56
conversion events in the same chromosome set indicates that
observed recombination events in P6 are strongly biased towards
conversions rather than crossovers. Among the studied chromo-
somes, intra-chromatid inversions are comparatively well ascer-
tained, because a crossover in the interval between any pair of
informative PSVs will be detected reliably. Conversion, however,
is under-ascertained because it is only observed when it transfers a
particular informative PSV. The scarcity of PSVs means that the
observed conversion: intra-chromatid crossover bias is actually an
underestimate of the true value. Additional uncertainty is
introduced by our inability to accurately identify co-conversion.

A bias towards non-crossovers is commonly observed in
recombination analysis. According to cytological studies the repair
of only 10% of DSBs in mammals results in crossovers, while the
remainder are assumed to be repaired as non-crossovers [37].
Most mammalian data on conversion: crossover ratios come from
studies of meiotic recombination hotspots in humans and mice.
The ratio varies significantly between different human hotspots
(from 2.7:1 at hotspot DNA3 to <1:12 at the β-globin hotspot);
there are also considerable differences among individuals, driven
in part by variation in trans-acting factors [38–41].

In comparing MSY gene conversion with conversion affecting
other chromosomes, its singular status as a constitutively haploid
chromosome must be remembered. As discussed above, both
inter- and intra-chromatid conversion can occur, but neither of
these processes is linked with the highly regulated ‘normal’
processes of synapsis and meiotic crossing over. Many questions
therefore remain about the timing and mechanism of MSY
conversion processes.

GC-bias in palindrome gene conversion

In a number of organisms recombination has been associated
with GC-bias arising from biased repair of mismatches in
heteroduplex DNA [42]. Consistent with this, we found evidence
of highly statistically significant GC-bias among the P6 gene
conversion events within the Y phylogeny.

We also asked whether GC-bias in gene conversion had a
deeper evolutionary history by comparing the patterns of
nucleotide replacements among human, chimpanzee and gorilla
P6 sequences. Spacers show a statistically significantly greater
proportion of replacements of S nucleotides by W nucleotides than
arms do (Table 2). This is true for both human- and chimpanzee-
specific nucleotide replacements. It is possible that these differ-
ences could be due to regional variation in GC-content, repeat
content, mutation rates or some other factors, but the observed
replacement patterns in palindrome arms are nonetheless consis-

tent with the action of GC-biased gene conversion. We might
expect the long-term action of such bias to lead to elevated GC-
content in arms compared to spacers. For P6, this is the case
(Table S6): 38.8% (arms) is significantly greater (p = 2.7×10⁻¹¹;
Chi-square test) than 37.0% (spacer). We can make the same
comparisons for other palindromes, setting aside P1 and P2, which
have very large arms and very small spacers. P3 also shows a
significant elevation of GC-content in its arms (p = 1.0×10⁻⁵⁶),
while P4, P5, P7 and P8 show no significant differences; the
pattern is therefore complex, but notably none of these
palindromes shows significantly higher GC-content in spacer
compared to arms. The observed differences could in principle
reflect the enrichment of protein-coding genes in palindrome arms
compared to spacers; however, the observed pattern persists when
the genes are removed (Table S6).

Apparent bias to ancestral state in gene conversion

Our comparisons of human, chimpanzee and gorilla P6
sequences concur with previous observations [2] in revealing
significantly lower inter-specific divergence among arms than
among spacers, in all three possible comparisons (Table 1). This
suggests either that the rate of initial mutation in arms is lower
than that in spacers, or that gene conversion is acting to
preferentially return new mutations arising in one arm to the
ancestral state, via conversion from the unmutated arm. Our
observation that individual gene conversion events among human
Y chromosomes are significantly biased towards retention of the
ancestral states of PSVs tends to support the second explanation.
Natural selection acting directly on the PSV sites seems an unlikely
explanation for the bias: examination of ENCODE [43] data (as
represented in the UCSC Genome Browser; April 2013) shows P6
to be devoid of functionally significant elements, apart from a 107-
bp snRNA gene in the arms ~20 kb proximal to the inner arm
boundary. There is no evidence for functional elements overlap-
ping the variants tested. An alternative explanation is that the
ancestral state bias emerges from the GC-bias. Notably, of the six
PSVs that are informative about GC-bias acting at individual sites,
five have a G or C nucleotide as their ancestral state. Whether
GC-bias provides a more general explanation for the conservation
of palindrome sequences will require more data on a larger
number of palindrome sequence variants.

Y-chromosomal palindromes are not alone in showing apparent
ancestral-state bias in conversion: comparison of human and
chimpanzee orthologs of an X-chromosomal palindrome [44] also
display significantly reduced interspecific divergence in arms
compared to spacers. This bias in conversion may therefore be a
general property of palindromic repeats. Its consequence is that
palindromes are ‘hard wired’ for conservation; although this will
be largely beneficial because most mutations are deleterious, it
may also ultimately limit adaptability of genes in palindromes by
limiting the opportunity for fixing beneficial mutations.

Future developments

Our understanding of the molecular evolution of the Y
chromosome would be greatly improved by the availability of
additional accurate sequences both from non-human primates and
humans. In principle, next-generation sequencing technologies
offer the opportunity to generate such sequences, but in practice
the complex repetitive structure of the Y chromosome means that
sequence assembly is impossible with current methods. Successful
generation of useful Y-chromosome sequences from humans and
other species [1,3,15,19] has required shot-gun sequencing of
assembled tiling arrays of BAC clones, an expensive and laborious
process. An additional problem is that genome sequencing projects
in non-human primates focus on females, in order to provide good coverage of the X chromosome. The structures of palindromes, the phase of variants within them, and gene conversion tract lengths will be illuminated by the advance of third-generation sequencing methods that have very long read lengths, and also high-throughput haplotyping of single sperm molecules, a method that has already proved successful in identifying the longest known allelic gene conversion tract of 22 kb [45].

# Materials and Methods

## DNA samples and Y haplotyping

We analysed a total of 378 male samples chosen from the CEPH-HGDP Cell Line Panel (Table S1) [15]. Choice was motivated by existing information on haplogroup, and practicality: we wanted to ensure representation of several members of each known haplogroup in order to detect gene conversion events (Figure 2), but to avoid analysing all 684 males in the panel due to the laborious nature of PSV typing and phasing.

Y-chromosome binary polymorphism data for HGDP samples were compiled as follows: 145 SNPs from CEPH 2011 (www. cephb.fr/en/hgdp/ - data supplement 10), 37 from Shi et al. 2010 [16] and Peter de Knijff (unpublished observations), 10 from Li et al. 2008 [18] and three from Sengupta et al. 2006 [17]. In addition, 23 SNPs (M112, M119, M150, M182, M217, M223, M231, M267, M285, M287, M3, M32, M35, M38, M6, M75, M78, M8, P15, P2, P37, P45 and P312) were typed as part of a GoldenGate custom genotyping assay (see section below). Eleven samples representing haplogroup K*(xL,M,N,O,P) [16] were typed for M254 and P204 using published PCR primers [13] and Sanger sequencing.

The whole dataset is described in Figure S1 and Table S1. For the phylogeny, the total of 200 mutational events gave rise to 122 possible Y-chromosome haplogroups, of which we observed 63 among the 378 samples analysed. Haplogroup nomenclature is as described [13], with shorthand names for some haplogroups, as described in Table S1.

There were two inconsistencies between data sources: (i) The phylogenetic relationships of markers P7 and M169 within hgB2 were consistent with the data of [46] rather than the original description [13]; (ii) Four samples (HGDP numbers 541, 542, 553 and 662) are listed in the data of [16] as belonging to hgK(xL,M1,NO,P), with the hg-M1-defining marker M106 ancestral; however, these same samples are listed under CEPH 2011 (www.cephb.fr/en/hgdp/ - data supplement 10) as derived for both M106 and the phylogenetically equivalent marker M189. Given that two markers are in agreement in this dataset, we regard them here as hgM1 chromosomes.

## Ethics statement

This study uses human DNA samples from the CEPH-HGDP panel, a widely available anonymised set of lymphoblastoid cell-lines (LCLs). The original paper describing this panel [15] states that the blood specimens that served as sources of the LCLs were freely donated under conditions of informed consent and confidentiality by reviewing consent forms, institutional review board approvals, or detailed reports from those who organised collections.

## Genotyping of P6 PSVs

The ten typed PSVs were labelled PSV1 to PSV10 based on their proximal-to-distal order on the proximal palindrome arm in the reference sequence (Figure 1b, Table S4).

As a convenient medium-throughput system for typing SN-PSVs, we chose the Illumina GoldenGate Genotyping Assay (Illumina, San Diego, CA). This system does not allow assay design or reliable calling for some variants in particular sequence contexts, and was eventually used for the successful typing of seven analysed PSVs (PSV2, 3, 5, 7-10). Experiments were carried out at the Genomics Core Facility of the University of Leicester. Genotypes were called with the Illumina GenomeStudio software version 3.1.0.0 (Illumina). Results were validated by Sanger sequencing of 5% of samples (n = 19) for each PSV (133 sequencing reads in total). PSV1 and PSV4 were typed by PCR-RFLP analysis using the restriction enzyme *Tst*I (Fermentas) for the former and *Hpy166II* (NEB) for the latter. PSV6 was typed by allele-specific PCR.

## Phasing of palindrome arms

In order to phase the palindrome arms an arm-specific long-range PCR approach was developed, using one universal primer binding within the arm and an arm-specific primer binding outside the outer palindrome boundary, generating a distal-arm-specific fragment of 18,893 bp incorporating seven of the studied PSVs. This fragment was then used as a template in nested PCR followed by re-typing of the seven PSVs. Five of the PSVs were typed by PCR-RFLP analysis using the following restriction enzymes (all NEB except PSV1): PSV1 - *Tst*I, PSV2 - *Aci*I, PSV3 - *Hinf*I, PSV4 - *Hpy166II* and PSV7 - *Mnl*I. Sanger sequencing and allele-specific PCR were used for PSV5 and PSV6, respectively. All primer sequences are listed in Table S5.

Arm-specific haplotypes were compared to the known phase of the human Y-chromosome reference sequence. In total 83 samples were examined and all but three found to have identical phase to the reference sequence (Table S1).

## PCR approaches to verify the presence of both palindrome arms

In order to ascertain the presence of both palindrome arms in samples with long apparently pseudohomozygous stretches, a paralog ratio test (PRT) [47] was developed. PRT primers were designed to amplify fluorescently labelled 390-bp test fragments from both arms of P6 (Figure 1b), plus a single 387-bp reference region from chromosome X (Table S5). Products were resolved and quantified using an ABI3130xl Genetic Analyzer and GeneMapper software v4.0 (Applied Biosystems, Carlsbad, CA). A normal male is expected to have two palindrome arms and one X chromosome, resulting in a test-to-reference ratio of 2:1. In total 50 samples were tested, each at least twice, including pseudoherozygous controls known to contain both palindrome arms (Table S1). All samples showed the expected ~2:1 ratio except one (HGDP00445), which showed a ratio of ~1:1. Semi-quantitative analysis using the amelogenin sex test [48], which simultaneously amplifies different-sized X- and Y-specific fragments, showed an X:Y ratio of 2:1, consistent with this cell-line having a 47,XXY karyotype.

The presence of both palindrome arms was also checked by an additional PCR-based approach. Firstly, PCR primers were designed to specifically extend over and amplify both the inner and outer boundaries of the palindrome. Secondly, long-range PCR primers were used to amplify ~10-kb fragments arm-specifically from the outer boundary of both arms followed by gel electrophoresis to check for changes in product length. The presence of all four palindrome boundaries and expected lengths of arm-specific PCR products was confirmed for all samples tested.
## Estimation of gene conversion rate

Mean gene conversion rate (assuming each converted SN-PSV represented an independent event) was estimated by dividing the number of conversion events (n), by the number of generations (g) encompassed in the phylogeny relating the 378 tested Y chromosomes. Estimation of g was based on a study [36] in which ~80 kb of DNA were resequenced in 47 Y chromosomes covering most of the major branches of the Y phylogeny to ascertain unbiased nucleotide divergence, revealing a total of 95 base substitutions. Assuming a TMRCA of 118,000 years (supported by more recent large-scale resequencing [49]), a generation time of 25 years, and a human-chimpanzee divergence time of 6.5 million years, the 47 chromosomes encompassed 52,000 generations [36]. The 378 Y chromosomes we studied also included most haplogroups in the phylogeny, but also multiple examples in individual haplogroups. We estimated the number of additional generations contributed by these: for the lower bound we assumed that all chromosomes sharing major haplogroups contributed no additional base substitutions in excess of the haplogroup-specific branch lengths; for the upper bound we assumed that each additional chromosome in a given haplogroup contributed an additional number of base substitutions equivalent to its descending from the root of the clade independently. This led to a range of total base substitutions of 323–935, corresponding to ~175,000–505,500 generations (Table S3).

## Sequencing and assembly of gorilla P6, and evolutionary comparisons

A partial consensus sequence of gorilla P6 was constructed from Illumina paired-end sequencing reads from: (i) whole genomes of two male gorillas giving an overall ~6× Y-chromosome coverage [20]; (ii) a sequence-capture library (SureSelect, Agilent), using a repeat-masked probe-design based on the human reference sequence, of a male gorilla giving a mean coverage of targetable portions of P6 of 232× (Protocol S1). Reads from all samples were mapped against the spacer and proximal arm of P6 in the human reference (GRCh37) and a consensus sequence for a given nucleotide called where it was covered by at least 5 concordant reads and minimum base quality score 20.

Chimpanzee MSY sequence used in interspecific comparisons was taken from [19]. Sequence alignments were performed using the web-based ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2) and Stretcher implemented in the EMBOSS package (http://emboss.sourceforge.net/).

## References

1. Skaletsky H, Kuroda-Kawaguchi T, Minx PJ, Cordum HS, Hillier L, et al. (2003) The male-specific region of the human Y chromosome: a mosaic of discrete sequence classes. Nature 423: 825–837.

2. Rozen S, Skaletsky H, Marszalek JD, Minx PJ, Cordum HS, et al. (2003) Abundant gene conversion between arms of massive palindromes in human and ape Y chromosomes. Nature 423: 873–876.

3. Hughes JF, Skaletsky H, Brown LG, Pyntikova T, Graves T, et al. (2012) Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes. Nature 483: 82–86.

4. Alföldi JE (2008) Sequence of the mouse Y chromosome. Cambridge, Massachusetts: Massachusetts Institute of Technology.

5. Geraldes A, Rambo T, Wing RA, Ferrand N, Nachman MW (2010) Extensive gene conversion drives the concerted evolution of paralogous copies of the SRY gene in European rabbits. Mol Biol Evol 27: 2437–2440.

6. Backström N, Ceplitis H, Berlin S, Ellegren H (2005) Gene conversion drives the evolution of HINTW, an ampliconic gene on the female-specific avian W chromosome. Mol Biol Evol 22: 1992–1999.

7. Davis JK, Thomas PJ, Thomas JW (2010) A W-linked palindrome and gene conversion in New World sparrows and blackbirds. Chromosome Res 18: 543–553.

## Supporting Information

**Dataset S1** Sequence alignment of human, chimpanzee and gorilla P6. (DOCX)

**Figure S1** Y phylogeny, showing markers typed and data sources. (TIF)

**Table S1** Genotyping results. (XLSX)

**Table S2** Summary of identified gene conversion events. (XLSX)

**Table S3** Estimation of gene conversion rate. (XLSX)

**Table S4** Genotypes of 10 studied PSVs in different primate species and genomic coordinates in the human reference sequence. (XLSX)

**Table S5** Primer sequences. (XLSX)

**Table S6** Comparison of GC-content between human chromosome arms and spacers. (XLSX)

**Protocol S1** Sequencing and analysis of gorilla Y-chromosome palindrome, containing summary of NGS sequence data and coverage. (DOCX)

## Acknowledgments

We thank Chris Tyler-Smith, Yali Xue, and Quan Long for access to sequence data from a haplogroup O3a Y chromosome; Peter de Knijff for supplying HGDP Y-haplogroup data; Rob Hardwick for help with design of the PRT; Pierpaolo Maisano Delser for statistical advice; Celia May for assistance with long-PCR design; the staff of NUCLEUS at the University of Leicester for assistance with GoldenGate genotyping; Chiara Batini, Dan Zadik and Edinburgh GenePool for assistance with Illumina sequencing of gorilla P6; Aylwyn Scally for sharing sequence data from two male gorillas; and three anonymous reviewers for helpful comments on our manuscript.

## Author Contributions

Conceived and designed the experiments: PH PB GRB MAJ. Performed the experiments: PH PB. Analyzed the data: PH PB MAJ. Contributed reagents/materials/analysis tools: PH GRB SB. Wrote the paper: PH MAJ.

8. Mendez-Lago M, Bergman CM, de Pablos B, Tracey A, Whitehead SL, et al. (2011) A large palindrome with interchromosomal gene duplications in the pericentromeric region of the D. melanogaster Y chromosome. Mol Biol Evol 28: 1967–1971.

9. Warburton PE, Giordano J, Cheung F, Gelfand Y, Benson G (2004) Inverted repeat structure of the human genome: the X-chromosome contains a preponderance of large, highly homologous inverted repeats that contain testes genes. Genome Res 14: 1861–1869.

10. Bagnall RD, Ayres KL, Green PM, Giannelli F (2005) Gene conversion and evolution of Xq28 duplicates involved in recurring inversions causing severe hemophilia A. Genome Res 15: 214–223.

11. Connallon T, Clark AG (2010) Gene duplication, gene conversion and the evolution of the Y chromosome. Genetics 186: 277–286.

12. Marais GA, Campos PR, Gordo I (2010) Can intra-Y gene conversion oppose the degeneration of the human Y chromosome? A simulation study. Genome Biol Evol 2: 347–357.

13. Karafet TM, Mendez FL, Meilerman M, Underhill PA, Zegura SL, et al. (2008) New binary polymorphisms reshape and increase resolution of the human Y-chromosomal haplogroup tree. Genome Res 18: 830–838.
14. Xue Y, Wang Q, Long Q, Ng BL, Swerdlow H, et al. (2009) Human Y chromosome base-substitution mutation rate measured by direct sequencing in a deep-rooting pedigree. Curr Biol 19: 1453-1457.

15. Cann HM, de Toma C, Cazes L, Legrand MF, Morel V, et al. (2002) A human genome diversity cell line panel. Science 296: 261-262.

16. Shi W, Ayub Q, Vermeulen M, Shao RG, Zuniga S, et al. (2010) A worldwide survey of human male demographic history based on Y-SNP and Y-STR data from the HGDP-CEPH populations. Mol Biol Evol 27: 385-393.

17. Sengupta S, Zhivotovsky LA, King R, Mehdi SQ, Edmonds CA, et al. (2006) Polarity and temporality of high-resolution Y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists. Am J Hum Genet 78: 202-221.

18. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, et al. (2008) Worldwide human relationships inferred from genome-wide patterns of variation. Science 319: 1100-1104.

19. Hughes JF, Skaletsky H, Pyntikova T, Graves TA, van Daalen SK, et al. (2010) Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content. Nature 463: 536-539.

20. Scally A, Dutheil JY, Hillier LW, Jordan GE, Goodhead I, et al. (2012) Insights into hominid evolution from the gorilla genome sequence. Nature 483: 169-175.

21. Marais G (2003) Biased gene conversion: implications for genome and sex evolution. Trends Genet 19: 330-338.

22. Chen JM, Cooper DN, Chuzhanova N, Ferec C, Patrinos GP (2007) Gene conversion: mechanisms, evolution and human disease. Nat Rev Genet 8: 762-775.

23. Assis R, Kondrashov AS (2012) Nonallelic gene conversion is not GC-biased in Drosophila or primates. Mol Biol Evol 29: 1291-1295.

24. Adams SM, King TE, Bosch E, Jobling MA (2006) The case of the unreliable SNP: recurrent back-mutation of Y-chromosomal marker P25 through gene conversion. Forensic Sci Int 159: 14-20.

25. Bosch E, Hurles ME, Navarro A, Jobling MA (2004) Dynamics of a human interparalog gene conversion hotspot. Genome Res 14: 835-844.

26. Bouzekri N, Taylor PG, Hammer MF, Jobling MA (1998) Novel mutation processes in the evolution of a haploid minisatellite, MSY1: array homogenization without homogenization. Hum Mol Genet 7: 655-659.

27. Trombetta B, Cruciani F, Underhill PA, Sellitto D, Scozzari R (2010) Footprints of X-to-Y gene conversion in recent human evolution. Mol Biol Evol 27: 714-725.

28. Rosser ZH, Balaresque P, Jobling MA (2009) Gene conversion between the X chromosome and the male-specific region of the Y chromosome at a translocation hotspot. Am J Hum Genet 85: 130-134.

29. Iwase M, Satta Y, Hirai H, Hirai Y, Takahata N (2010) Frequent gene conversion events between the X and Y homologous chromosomal regions in primates. BMC Evol Biol 10: 225.

30. Iwase M, Satta Y, Hirai Y, Hirai H, Imai H, et al. (2003) The amelogenin loci span an ancient pseudoautosomal boundary in diverse mammalian species. Proc Natl Acad Sci U S A 100: 5258-5263.

31. Marais G, Galtier N (2003) Sex chromosomes: how X-Y recombination stops. Curr Biol 13: R641-643.

32. Pecon Slattery J, Sanner-Wachter L, O'Brien SJ (2000) Novel gene conversion between X-Y homologues located in the nonrecombining region of the Y chromosome in Felidae (Mammalia). Proc Natl Acad Sci USA 97: 5307-5312.

33. Mancera E, Bourgon R, Brozzi A, Huber W, Steinmetz LM (2008) High-resolution mapping of meiotic crossovers and non-crossovers in yeast. Nature 454: 479-485.

34. Jeffreys AJ, May CA (2004) Intense and highly localized gene conversion activity in human meiotic crossover hot spots. Nature Genet 36: 152-156.

35. Lange J, Skaletsky H, van Daalen SK, Embry SL, Korver CM, et al. (2009) Isodicentric Y chromosomes and sex disorders as byproducts of homologous recombination that maintains palindromes. Cell 138: 855-869.

36. Repping S, van Daalen SK, Brown LG, Korver CM, Lange J, et al. (2006) High mutation rates have driven extensive structural polymorphism among human Y chromosomes. Nat Genet 38: 463-467.

37. Guillon H, Baudat F, Grey C, Liskay RM, de Massy B (2005) Crossover and noncrossover pathways in mouse meiosis. Mol Cell 20: 563-573.

38. Jeffreys AJ, Holloway JK, Kauppi L, May CA, Neumann R, et al. (2004) Meiotic recombination hot spots and human DNA diversity. Philos Trans R Soc Lond B Biol Sci 359: 141-152.

39. Holloway K, Lawson VE, Jeffreys AJ (2006) Allelic recombination and de novo deletions in sperm in the human beta-globin gene region. Hum Mol Genet 15: 1099-1111.

40. Berg IL, Neumann R, Sarbajna S, Odenthal-Hesse L, Butler NJ, et al. (2011) Variants of the protein PRDM9 differentially regulate a set of human meiotic recombination hotspots highly active in African populations. Proc Natl Acad Sci U S A 108: 12378-12383.

41. Sarbajna S, Denniff M, Jeffreys AJ, Neumann R, Soler Artigas M, et al. (2012) A major recombination hotspot in the XqYq pseudoautosomal region gives new insight into processing of human gene conversion events. Hum Mol Genet 21: 2029-2038.

42. Duret L, Galtier N (2009) Biased gene conversion and the evolution of mammalian genomic landscapes. Annu Rev Genomics Hum Genet 10: 285-311.

43. ENCODE Project Consortium (2012) An integrated encyclopedia of DNA elements in the human genome. Nature 489: 57-74.

44. Losch FO, Bredenbeck A, Hollstein VM, Walden P, Wrede P (2007) Evidence for a large double-cruciform DNA structure on the X chromosome of human and chimpanzee. Hum Genet 122: 337-343.

45. Wang J, Fan HC, Behr B, Quake SR (2012) Genome-wide single-cell analysis of recombination activity and de novo mutation rates in human sperm. Cell 150: 402-412.

46. Batini C, Ferri G, Destro-Bisol G, Brisighelli F, Luiselli D, et al. (2011) Signatures of the pre-agricultural peopling processes in sub-Saharan Africa as revealed by the phylogeography of early Y chromosome lineages. Mol Biol Evol 28: 2603-2613.

47. Armour JA, Palla R, Zeeuwen PL, den Heijer M, Schalkwijk J, et al. (2007) Accurate, high-throughput typing of copy number variation using paralogue ratios from dispersed repeats. Nucleic Acids Res 35: e19.

48. Sullivan KM, Mannucci A, Kimpton CP, Gill P (1993) A rapid and quantitative DNA sex test - fluorescence-based PCR analysis of X-Y homologous gene amelogenin. Biotechniques 15: 636.

49. Wei W, Ayub Q, Chen Y, McCarthy S, Hou Y, et al. (2013) A calibrated human Y-chromosomal phylogeny based on resequencing. Genome Res 23: 388-395.