The Divergence of Neandertal and Modern Human Y Chromosomes

Fernando L Mendez; G David Poznik; Sergi Castellano; Carlos D Bustamante

doi:10.1016/j.ajhg.2016.02.023

. 2016 Apr 7;98(4):728–734. doi: 10.1016/j.ajhg.2016.02.023

The Divergence of Neandertal and Modern Human Y Chromosomes

Fernando L Mendez ^1,^∗, G David Poznik ^1,², Sergi Castellano ³, Carlos D Bustamante ^1,^4,^∗∗

PMCID: PMC4833433 PMID: 27058445

Abstract

Sequencing the genomes of extinct hominids has reshaped our understanding of modern human origins. Here, we analyze ∼120 kb of exome-captured Y-chromosome DNA from a Neandertal individual from El Sidrón, Spain. We investigate its divergence from orthologous chimpanzee and modern human sequences and find strong support for a model that places the Neandertal lineage as an outgroup to modern human Y chromosomes—including A00, the highly divergent basal haplogroup. We estimate that the time to the most recent common ancestor (TMRCA) of Neandertal and modern human Y chromosomes is ∼588 thousand years ago (kya) (95% confidence interval [CI]: 447–806 kya). This is ∼2.1 (95% CI: 1.7–2.9) times longer than the TMRCA of A00 and other extant modern human Y-chromosome lineages. This estimate suggests that the Y-chromosome divergence mirrors the population divergence of Neandertals and modern human ancestors, and it refutes alternative scenarios of a relatively recent or super-archaic origin of Neandertal Y chromosomes. The fact that the Neandertal Y we describe has never been observed in modern humans suggests that the lineage is most likely extinct. We identify protein-coding differences between Neandertal and modern human Y chromosomes, including potentially damaging changes to PCDH11Y, TMSB4Y, USP9Y, and KDM5D. Three of these changes are missense mutations in genes that produce male-specific minor histocompatibility (H-Y) antigens. Antigens derived from KDM5D, for example, are thought to elicit a maternal immune response during gestation. It is possible that incompatibilities at one or more of these genes played a role in the reproductive isolation of the two groups.

Introduction

A central goal of human population genetics and paleoanthropology is to elucidate the relationships among ancient populations. Before the emergence of anatomically modern humans in the Middle Pleistocene ∼200 thousand years ago (kya),¹ archaic humans lived across Africa, Europe, and Asia in highly differentiated populations. Modern human populations that expanded out of Africa in the Upper Pleistocene received a modest genetic contribution from at least two archaic hominin groups, the Neandertals and Denisovans.²^,³^,⁴^,⁵ Especially in light of hypothesized genetic incompatibilities between Neandertals and modern humans,⁶ it is important to characterize differentiation between their ancestral populations and to investigate potential barriers to gene flow.

When populations diverge from one another, each retains a subset of the variation that existed in the ancestral population. Consequently, sequence divergence times usually exceed population divergence times, and this effect is more pronounced when the ancestral effective population size was large. In humans, a large fraction of genetic diversity is due to ancient polymorphisms that arose long before the emergence of anatomically modern traits. As a result, Neandertal and modern haplotypes are often no more diverged than modern human sequences are among themselves.² This fact complicates the search for introgressed genomic segments, but two features facilitate their detection.⁶^,⁷ First, due to low levels of polymorphism among Neandertals,⁵ introgressed sequences are often quite similar to those of the Neandertal reference. Second, these regions have elevated linkage disequilibrium due to the relatively recent date of admixture, ∼50 kya.⁸^,⁹^,¹⁰ Although introgressed Neandertal sequences have been identified in modern human autosomes and X chromosomes, no mitochondrial genome (mtDNA) sequences of Neandertal origin have been reported in modern humans, and Neandertal Y-chromosome sequences have not yet been characterized.

Because uniparentally inherited loci have much smaller effective population sizes than autosomal or X-linked loci, the expected differences between sequence and population divergence times are smaller. Therefore, studying these loci can help to delineate an upper bound for the time at which populations last exchanged genetic material. To date, five Neandertal individuals have been whole-genome sequenced to 0.1× coverage or higher,²^,⁵ but all were female. Full mtDNA sequences are also available for eight individuals from Spain, Germany, Croatia, and Russia,¹¹^,¹² but the relationship between Neandertal and modern human Y chromosomes remains unknown.

In this work, we analyzed ∼120 kb of exome-captured Y-chromosome sequence from an ∼49,000-year-old (uncalibrated ¹⁴C)¹³ Neandertal male from El Sidrón, Spain.¹⁴ We compare it to the human and chimpanzee reference sequences and to the sequences of two Mbo individuals¹⁵ who carry the A00 haplogroup, the most deeply branching group known.¹⁶ We identify the relationship between the Neandertal and modern human Y chromosomes and estimate the time to their most recent common ancestor (TMRCA). We also examine coding differences and explore their potential significance for reproductive isolation.

Material and Methods

Sequence Data and Processing

We used the Y-chromosome sequences from the exome capture of a Neandertal from El Sidrón, Spain,¹⁴ and we downloaded the complete sequences of two A00 Y chromosomes.¹⁵ The Neandertal data included coding, non-coding, and off-target sequences, and all three sequences were mapped against the GRCh37 reference.¹⁴ Given that the A00 sequences were closely related,¹⁵^,¹⁶ we merged them to increase coverage. We called bases for both the Neandertal and A00 sequences by using SAMtools mpileup (v.1.1),¹⁷ specifying input options to count anomalous read pairs (-A), recalculate base qualities (-E), and filter out poor-quality bases (-Q 17) and poorly mapping reads (-q 20).

We then identified overlapping regions and excluded coordinates with unusually high coverage, filtering out sites with coverage greater than the mean plus five times its square root (Figure S1). Under a Poisson model, this cutoff would elicit the loss of less than one genuine site per 10,000. Finally, we removed sites with inconsistent base calls, discarding those with more than two reads differing from the consensus allele and those for which more than one third of the observed bases did not match the consensus. This filter should minimize the effects of postmortem DNA damage and of modern contamination.

Using the blastz file chrY.hg19.panTro4.net.axt.gz,¹⁸ we identified the subset of regions within which the human sequences align to the chimpanzee reference. This yielded a total of 118,643 base pairs (bp). In what follows, we refer to this set of sites as “filter 1.” We also identified a second, more restrictive, set of regions totaling 100,324 bp, “filter 2,” by further requiring that the alignment correspond to the chimpanzee Y chromosome rather than to another chimpanzee chromosome (Tables S1A and S1B).

For each position within these regions, we determined whether the Neandertal, A00, or both differed from the human reference sequence. We then used the corresponding chimpanzee allele as a proxy for the ancestral state in order to assign the mutation to the appropriate branch of the tree relating the four sequences (Figure 1A). In doing so, we discarded five sites: two at which the chimpanzee carries a third allele, one for which the chimpanzee carries a deletion, and two that were specific to A00 but only supported by a single read. Excluding these sites had little impact on our analyses.

Tree Inference

(A) A priori, three trees could feasibly have related the Y chromosomes of the chimpanzee (Chimp), the Neandertal (Neander), haplogroup A00, and the human reference (Ref). Mutations on branch a support topology i, with the Neandertal lineage as the outgroup to those of modern humans, whereas mutations on branches b and c support topologies ii and iii, respectively. Branches d, e, and f correspond to mutations private to individual lineages.

(B) Counts of SNVs consistent with each branch. Columns refer to sets of coordinates considered (see Materials and Methods). Incompatible sites are those that cannot be explained by a single mutation on any of the three trees.

Estimating TMRCA

To estimate the TMRCA of the Neandertal and modern human Y chromsomes (T_NR), we decomposed this quantity (Figure 2) into the sum of the TMRCA of modern humans (T_AR) and the time separating the most recent common ancestor of modern humans from its common ancestor with the Neandertal lineage (T_NM):

T_{NR} = T_{AR} + T_{NM} = α T_{AR} α \equiv (1 + \frac{T_{NM}}{T_{AR}}) .

We then estimated T_AR and used two methods to estimate α.

Estimating the TMRCA of Neandertal and Modern Y Chromosomes

The quantity of primary interest is T_NR = T_NM+ T_AR. Branches are labeled as in Figure 1, and “M” denotes the most recent common ancestor of modern human lineages.

To estimate T_AR, we used sequence data from the ancient Ust’-Ishim sample,⁹ first applying the filters described for the A00 sequences. To reduce the potential impact of postmortem DNA damage, we restricted this analysis to coordinates covered by at least three sequencing reads. We further restricted to the subset of Poznik et al.¹⁹ regions in which the human reference sequence is based on bacterial artificial chromosome clones derived from the RP-11 individual,²⁰ a known carrier of haplogroup R1b. This left ∼7.83 Mb of sequence within which to assign variants to the appropriate branches (Figure S2, Appendix A). Using the known age of the Ust’-Ishim individual and the constrained optimization procedure described in Rasmussen et al.,²¹ we obtained parametric bootstrap estimates for T_AR as well as for the mutation rate and the TMRCA of haplogroup K-M526 (Appendix A). Briefly, we sampled from the process that generated the observed tree (Figure S2) by simulating the number of single nucleotide variants (SNVs) on each branch as a Poisson draw with mean equal to the observed number of mutations. To obtain bootstrap samples of the three parameters, we maximized their joint likelihood for each tree replicate.

In our first approach to estimate α, we used the relative numbers of mutations assigned to branches a, d, and e (Figure 1), assigning the four sites that did not fit the consensus topology to the A00 or reference lineages, as appropriate (Appendix A). The proportion of time represented by branch a is:

\frac{T_{a}}{T_{a} + T_{d} + T_{e}} = \frac{T_{NM}}{T_{NM} + 2 T_{AR}} = \frac{(α - 1) T_{AR}}{(α - 1) T_{AR} + 2 T_{AR}} = \frac{α - 1}{α + 1} .

Therefore, assuming a time-homogeneous mutation rate, the number of branch-a mutations is binomially distributed with parameters p = (α – 1) ∕ (α + 1) and n equal to the total number of mutations. Estimating p from the data leads directly to a point estimate and confidence interval (CI) for α. This first method has the appealing property that it is independent of both the mutation rate and the absolute values of the times. However, the estimation error might be suboptimal due to uncertainty in both the numerator and the denominator.

In the second method, we estimated α via the ratio $T_{NM} / T_{AR}$ , making use of the fact that we can estimate T_AR with greater certainty than we can T_NM. To estimate T_NM, we restricted our attention to sequences overlapping the ∼8.8 Mb of sequence analyzed by Karmin et al.,¹⁵ leaving 80,420 bp (“filter 3”) or 75,596 bp (“filter 4”) when confining our analysis to those sites that passed filter 2 (Figure 1, Tables S2a and S2b). Let l equal the total length of sequence under consideration (e.g., 80.42 kb), let μ equal the mutation rate over the full 8.8 Mb, let r equal the ratio of the mutation rate within the smaller region to that of the larger, and let s equal the number of mutations shared by A00 and the reference sequence within the smaller region. With these, we constructed the estimator ${\hat{T}}_{NM} = s / (l r μ)$ . Similarly, let L equal the subset of the 8.8 Mb for which the A00 sequence had 3× or greater coverage (also ∼8.8 Mb), and let S equal the number of mutations unique to either the reference sequence or to A00 over the entire 8.8 Mb. We can estimate T_AR with ${\hat{T}}_{AR} = S / (2 L μ)$ and α with:

\hat{α} = (1 + \frac{{\hat{T}}_{NM}}{{\hat{T}}_{AR}}) = (1 + \frac{2 s L}{\hat{r} S l}) .

We estimated r by comparing the number of mutations unique to a single branch of the Y-chromosome tree of Karmin et al.,¹⁵ both within the full 8.8-Mb region and within the ∼80-kb subset. These numbers, 32,853 and 279 (238 under filter 4), respectively, correspond to a relative mutation rate of 0.93 (95% CI: 0.82–1.04) (0.84 under filter 4 [95% CI: 0.74–0.95]). Because selection has the strongest effect on lower frequency mutations, we also estimated r by using only shared variants, and this yielded nearly identical point estimates.

Finally, to construct a CI for α, we sampled values of s and S from Poisson distributions with means equal to the observed numbers of mutations, and we sampled rl/L as the ratio of two Poisson random variables with means equal to 279 (or 238) and 32,853, respectively.

Functional Variation

We determined whether each mutation overlaps with annotated RefSeq genes and whether it overlaps with coding sequence (Figure 1, Data S5). For each coding SNV, we determined whether the mutation results in silent, missense, or nonsense mutations, but we did not consider frameshift mutations. For each coding non-synonymous mutation, we used the HumDiv model of PolyPhen-2 to evaluate ancestral-to-derived changes and MutationTaster to evaluate reference-to-alternative changes. We report findings from all sites for which these programs were able to make predictions.

Results

With the chimpanzee Y chromosome as the outgroup, three tree topologies could have related the lineages of the Neandertal, haplogroup A00, and the human reference (Figure 1A). To identify which of the three was consistent with the data, the key question was which of the three possible pairs of sequences is the most closely related. Of 118,643 sites (Figure 1B, filter 1) for which we had Neandertal data and human-chimpanzee reference alignments,¹⁸ we identified 24 biallelic SNVs for which the Neandertal sequence shared the chimpanzee allele and differed from both A00 and the human reference. In contrast, the chimpanzee and A00 sequences shared just four SNVs not present in the other sequences, and the chimpanzee and human reference sequences shared zero. Taken together, these data strongly support the tree that places the Neandertal Y as the most distantly related to the others (Figure 1A, tree i). Two of the four variants that are inconsistent with this topology are known to segregate within modern humans and are therefore the result of recurrent mutations or contamination (Appendix A).

Upon elucidating the topology of the tree relating the Neandertal Y chromosome to those of modern humans, our next goal was to estimate the divergence time. We decomposed the TMRCA, T_NR, into the sum of two intervals (Figure 2): the TMRCA of A00 and the reference, T_AR, and the time between their common ancestor and the common ancestor of the Neandertal lineage and that of modern humans, T_NM. To estimate T_NR, we estimated T_AR and the ratio α ≡ T_NR ∕ T_AR, taking care to consider uncertainty both in the mutation rate and in the expected number of mutations to construct a CI. Because the numbers of mutations that accumulate on the branches of the tree are conditionally independent of one another and are nearly uncorrelated with the estimator of T_AR, we estimated α and T_AR independently (see Materials and Methods).

Leveraging data from an ∼45,000-year-old Siberian (Ust’-Ishim),⁹ we estimated that T_AR = 275 kya (95% CI: 241–305 kya), and we estimated α by using two approaches that yielded similar results. In our first approach, we simply used the number of mutations shared by A00 and the reference (branch a of Figure 2) and the number of mutations unique to each (branches d and e) to estimate the relative times between splits. This method is insensitive to mutation-rate variability across the chromosome and led us to estimate α = 2.14 (95% CI: 1.64–2.89). In the second approach, we made use of the greater amount of data available for the denominator of the ratio and adjusted for mutation rate heterogeneity across the chromosome to estimate α = 1.82 (95% CI: 1.40–2.32). Because the main source of uncertainty is the limited sequence coverage for the Neandertal lineage, the CIs from the two approaches overlap substantially, but we prefer the first method, as it is simpler and potentially less biased. In both cases, we disregarded the number of variants unique to the Neandertal sequence (branch f) because this branch is enriched for false positives as a result of low coverage, DNA damage, and sequencing errors.

Combining the parametric bootstrap CIs of α and T_AR, we estimated T_NR = 588 kya (95% CI: 447–806 kya) with the first α estimate and T_NR = 499 kya (95% CI: 375–656 kya) with the second.

Finally, we examined the potential functional relevance of the 146 mutations that differed among the Neandertal, A00, and reference sequences (Data S5). These included 11 non-synonymous changes and one nonsense mutation (Table 1). PolyPhen-2²² predicted most missense mutations to have a benign effect, but it predicted possibly or probably damaging effects for Neandertal mutations in PCDH11Y (MIM: 400022) and USP9Y (MIM: 400005), for an A00 mutation in ZFY (MIM: 490000), and for a modern human mutation in KDM5D (MIM: 426000). The Neandertal nonsense mutation at codon 16 of TMSB4Y (MIM: 400017) might render its product non-functional, and MutationTaster²³ predicts that it is probably deleterious.

Table 1.

Protein-Changing Mutations

Coordinate	Gene	Lineage¹	Substitution²	Effect	Tool	Function	MIM No.
2,844,774	ZFY	N	p.Val140Ala	B	P2	potential transcription factor	MIM: 490000
p.Val331Ala	B	P2
2,847,322	ZFY	A	p.Ile374Thr	B	P2	potential transcription factor	MIM: 490000
p.Ile488Thr	PrD	P2
p.Ile565Thr	B	P2
4,967,724	PCDH11Y³	N	p.Lys702Thr	B, B	P2, MT	protocadherin	MIM: 400022
5,605,569	PCDH11Y³	N	p.Ser1203Arg	PrD	P2	protocadherin	MIM: 400022
6,932,032	TBL1Y	N	p.Gly100Ala	B, B	P2, MT	–	MIM: 400033
14,832,610	USP9Y	N	p.Glu62Gly	PrD	P2	peptidase	MIM: 400005
14,832,620	USP9Y	R	p.Glu65Asp	B	P2	peptidase	MIM: 400005
14,838,553	USP9Y	N	p.Ala162Thr	B	P2	peptidase	MIM: 400005
15,816,262	TMSB4Y	N	p.Ser16^∗	PrD	MT	actin sequestration	MIM: 400017
21,868,167	KDM5D	R, A	p.Arg1445Gln	B, B	P2, MT	demethylase	MIM: 426000
p.Arg1388Gln	B, B	P2, MT
p.Arg1476Gln	B, B	P2, MT
21,905,071	KDM5D	R, A	p.Ile69Val	PoD, B	P2, MT	demethylase	MIM: 426000
23,545,399	PRORY	A	p.Arg125Cys	–	–	–	–

Open in a new tab

Please see Data S5. Mutations variable among Neanderthal, A00, and the reference sequence, Data S6. Alignment (SAM) file with Neanderthal reads overlapping mutations in Table S3, Data S7. Neanderthal pileup lines for positions indicated in Table S3 for additional information on all mutations. Abbreviations are as follows: N, Neandertal; A, A00; R, reference; B, benign; PoD, possibly damaging; PrD, probably damaging; P2, PolyPhen-2 (ancestral to derived); MT, MutationTaster (reference to alternative).

Lineage(s) bearing the derived allele.

Multiple listings for a single coordinate reflect substitutions in different transcripts of the gene.

See Appendix B.

Discussion

We have estimated that the Neandertal Y chromosome from El Sidrón diverged from those of modern humans ∼590 kya, a value similar to TMRCA estimates for mtDNA sequences: 400 kya to 800 kya.¹¹^,¹² This time estimate and the genealogy we have inferred strongly support the notion that the most recent common ancestor of these Y chromosomes belonged to the population from which Neandertals and modern humans diverged, thereby refuting three alternative hypotheses. A priori, the Neandertal Y could have introgressed from a super-archaic population⁵ (Figure 3, scenario a), but this would have led to a far greater TMRCA estimate. Alternatively, it could have introgressed from the ancestors of modern humans after their divergence from Neandertals and prior to the most recent common ancestor of present-day Y chromosomes (scenario b) or from modern human populations subsequent to their migrations out of Africa (scenario c). We can also reject these hypotheses, as each requires a more recent split time.

Relationship of Neandertal Y Chromosome to Those of Modern Humans

The genealogy (red tree) can be parsimoniously explained as mirroring the population divergence (gray tree). We find no evidence for (a) a highly divergent super-archaic origin of the Neandertal Y chromosome, (b) ancient gene flow post-dating the population split, or (c) relatively recent introgression of a modern human Y chromosome into the Neandertal population.

The fact that the Neandertal Y-chromosome lineage we describe has never been observed in modern humans suggests that the lineage is most likely extinct. Although the Neandertal Y chromosome (and mtDNA) might have simply drifted out of the modern human gene pool,²⁴ it is also possible that genetic incompatibilities contributed to their loss. In comparing the Neandertal lineage to those of modern humans, we identified four coding differences with predicted functional impacts, three missense and one nonsense (Table 1). Three mutations—within PCDH11Y, USP9Y, and TMSB4Y—are unique to the Neandertal lineage, and one, within KMD5D, is fixed in modern human sequences. The first gene, PCDH11Y, resides in the X-transposed region of the Y chromosome. Together with its X-chromosome homolog PCDH11X, it might play a role in brain lateralization and language development.²⁵ The second gene, USP9Y, has been linked to ubiquitin-specific protease activity²⁶ and might influence spermatogenesis.²⁷ Expression of the third gene, TMSB4Y, might reduce cell proliferation in tumor cells, suggesting tumor suppressor function.²⁸ Finally, the fourth gene, KDM5D, encodes a lysine-specific demethylase whose activity suppresses the invasiveness of some cancers.²⁹

Polypeptides from several Y-chromosome genes act as male-specific minor histocompatibility (H-Y) antigens that can elicit a maternal immune response during gestation. Such effects could be important drivers of secondary recurrent miscarriages³⁰ and might play a role in the fraternal birth order effect of male sexual orientation.³¹ Interestingly, all three genes with potentially functional missense differences between the Neandertal and modern humans sequences are H-Y genes, including KDM5D, the first H-Y gene characterized.³² It is tempting to speculate that some of these mutations might have led to genetic incompatibilities between modern humans and Neandertals and to the consequent loss of Neandertal Y chromosomes in modern human populations. Indeed, reduced fertility or viability of hybrid offspring with Neandertal Y chromosomes is fully consistent with Haldane’s rule, which states that “when in the [first generation] offspring of two different animal races one sex is absent, rare, or sterile, that sex is the [heterogametic] sex.”³³

Acknowledgments

F.L.M. was supported by the Stanford Center for Computational, Evolutionary, and Human Genomics (CEHG). This work was supported by the National Science Foundation grant no. 1201234. G.D.P. was supported by the National Science Foundation Graduate Research Fellowship under grant no. DGE-1147470 and by the National Library of Medicine training grant LM-007033. S.C. was supported by the Max Planck Society. We thank Bence Viola for lending expertise and Chris Gignoux for helpful comments and suggestions. G.D.P. is an employee of 23andMe. C.D.B is on the scientific advisory boards (SABs) of AncestryDNA, BigDataBio, Etalon DX, Liberty Biosecurity, and Personalis. He is also a founder and SAB chair of IdentifyGenomics. None of these entities played a role in the design, execution, interpretation, or presentation of this study.

Published: April 7, 2016

Footnotes

Supplemental Data include two figures, seven data files, and a Spanish translation of this article and can be found with this article online at http://dx.doi.org/10.1016/j.ajhg.2016.02.023.

Contributor Information

Fernando L. Mendez, Email: [email protected].

Carlos D. Bustamante, Email: [email protected].

Appendix A

Incompatible and Recurrent Variants and T_AR

In estimating T_AR, we initially removed 46 sites with no chimpanzee alignment, 70 sites at which the chimpanzee base differed from both human bases, 13 sites at which the chimpanzee and reference sequences agreed (to the exclusion of the other two lineages), and 4 sites at which the chimpanzee and Ust’-Ishim lineages agreed (to the exclusion of the others). For the first two sets, A00 differs from the reference, so we could partition the mutations as either specific to the reference (branch d of Figure S2) or to the union of branches a and f (af). In both the set of 46 and the set of 70, the relative numbers of mutations assigned to branches d and af are consistent with those of the sites for which the chimpanzee data were conclusive (Fisher’s exact test p values: 0.82 and 0.13).

The 17 sites that are incompatible with the tree are principally due to recurrent and back mutations. Because the reference has accumulated more mutations than Ust’-Ishim since their common ancestor, it is expected that more incompatible sites unite A00 and Ust’-Ishim than unite the reference and A00. Indeed, 10 of the 13 mutations map to branches ancestral to the reference sequence (but not Ust’-Ishim) in the 1000 Genomes Project (see Web Resources). Likewise, one of the other four mutations could have recurred in the reference and A00. Our approach cannot detect mutations that occurred on both the lineage leading to A00 and on the lineage leading to K-M526; however, the expected number of such mutations is quite small.

Including the 116 mutations from the first two sets lowers the T_AR estimate from 287 kya (95% CI: 252–321 kya) to 284 kya (95% CI: 249–316 kya) and including the additional 11 mutations from the third and fourth sets lowers it further, to 275 kya (95% CI: 241–305 kya). However, this last estimate is most likely biased slightly downward due to the impossibility of observing mutations that recurred in the ancestors of A00 and of K-M526.

Mutation Rate and TMRCA of K-M526

With the corrections described above, and assuming that the age of Ust’-Ishim is 45,000 years,^⁹ we estimate the mutation rate in the analyzed region to be 0.78 × 10⁻⁹ mutations per bp per year (95% CI: 0.71–0.89 × 10⁻⁹ mutations per bp per year), and we estimate the TMRCA of K-M526 to be 48.1 kya (95% CI: 46.4–49.6 kya). The effective correction due to including the 127 mutations described above was small.

Recurrent Mutations

Four mutations were inconsistent with tree ii in Figure 1A. Non-A0 lineages in the 1000 Genomes panel^³⁴ share the reference allele at coordinate 2,710,154, and individuals in haplogroups B through T share the reference allele at 23,558,260. The two others were at coordinates 9,386,241 and 15,024,530.

Appendix B

The X-transposed region of the Y chromosome arose from the transposition of an ∼3.5-Mb stretch of the X chromosome at some point subsequent to the divergence of human and chimpanzee lineages.^²⁰ Due to sequence similarity of ∼99%, short-read mapping is often ambiguous in this region, but we were able to use accumulated sequence divergence to manually assess reads that mapped poorly to PCDH11Y.

The probably damaging functional mutation at GRCh37 coordinate 5,605,569 is flanked by two bases that differ between the X and Y chromosomes, at positions 5,605,520 and 5,605,622. Seven sequencing reads overlapped the putative functional mutation and at least one of these two Y-specific bases, and each supported a derived allele call for the Neandertal lineage. Thus, despite the fact that just one of the seven reads mapped with high quality, we had sufficient evidence to call the derived genotype. Furthermore, the only read that carried the ancestral allele at the functional site, and overlapped one of the two diagnostic sites, bore the X-specific base.

Web Resources

The URLs for data presented herein are as follows:

1000 Genomes Project Y Chromosome Supporting Data, ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/supporting/chrY/
MutationTaster, http://www.mutationtaster.org/
OMIM, http://www.omim.org/
PolyPhen-2, http://genetics.bwh.harvard.edu/pph2/
UCSC Genome Browser, http://genome.ucsc.edu

Supplemental Data

Document S1. Figures S1 and S2

mmc1.pdf^{(109.7KB, pdf)}

Data S1. Regions passing Filter 1

mmc2.zip^{(13.2KB, zip)}

Data S2. Regions passing Filter 2

mmc3.zip^{(10.8KB, zip)}

Data S3. Regions passing Filter 3

mmc4.zip^{(9.1KB, zip)}

Data S4. Regions passing Filter 4

mmc5.zip^{(8.4KB, zip)}

Data S5. Mutations variable among Neanderthal, A00, and the reference sequence

mmc6.zip^{(1.7KB, zip)}

Data S6. Alignment (SAM) file with Neanderthal reads overlapping mutations in Table S3

mmc7.zip^{(30.3KB, zip)}

Data S7. Neanderthal pileup lines for positions indicated in Table S3

mmc8.zip^{(1.4KB, zip)}

Document S2. Spanish Translation: La Divergencia de Los Cromosomas Y de Neandertal y de los Humanos Modernos

mmc9.pdf^{(1MB, pdf)}

Document S3. Article plus Supplemental Data

mmc10.pdf^{(433.9KB, pdf)}

References

1.McDougall I., Brown F.H., Fleagle J.G. Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature. 2005;433:733–736. doi: 10.1038/nature03258. [DOI] [PubMed] [Google Scholar]
2.Green R.E., Krause J., Briggs A.W., Maricic T., Stenzel U., Kircher M., Patterson N., Li H., Zhai W., Fritz M.H.-Y. A draft sequence of the Neandertal genome. Science. 2010;328:710–722. doi: 10.1126/science.1188021. [DOI] [PMC free article] [PubMed] [Google Scholar]
3.Reich D., Green R.E., Kircher M., Krause J., Patterson N., Durand E.Y., Viola B., Briggs A.W., Stenzel U., Johnson P.L.F. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature. 2010;468:1053–1060. doi: 10.1038/nature09710. [DOI] [PMC free article] [PubMed] [Google Scholar]
4.Meyer M., Kircher M., Gansauge M.-T., Li H., Racimo F., Mallick S., Schraiber J.G., Jay F., Prüfer K., de Filippo C. A high-coverage genome sequence from an archaic Denisovan individual. Science. 2012;338:222–226. doi: 10.1126/science.1224344. [DOI] [PMC free article] [PubMed] [Google Scholar]
5.Prüfer K., Racimo F., Patterson N., Jay F., Sankararaman S., Sawyer S., Heinze A., Renaud G., Sudmant P.H., de Filippo C. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014;505:43–49. doi: 10.1038/nature12886. [DOI] [PMC free article] [PubMed] [Google Scholar]
6.Sankararaman S., Mallick S., Dannemann M., Prüfer K., Kelso J., Pääbo S., Patterson N., Reich D. The genomic landscape of Neanderthal ancestry in present-day humans. Nature. 2014;507:354–357. doi: 10.1038/nature12961. [DOI] [PMC free article] [PubMed] [Google Scholar]
7.Vernot B., Akey J.M. Resurrecting surviving Neandertal lineages from modern human genomes. Science. 2014;343:1017–1021. doi: 10.1126/science.1245938. [DOI] [PubMed] [Google Scholar]
8.Sankararaman S., Patterson N., Li H., Pääbo S., Reich D. The date of interbreeding between Neandertals and modern humans. PLoS Genet. 2012;8:e1002947. doi: 10.1371/journal.pgen.1002947. [DOI] [PMC free article] [PubMed] [Google Scholar]
9.Fu Q., Li H., Moorjani P., Jay F., Slepchenko S.M., Bondarev A.A., Johnson P.L.F., Aximu-Petri A., Prüfer K., de Filippo C. Genome sequence of a 45,000-year-old modern human from western Siberia. Nature. 2014;514:445–449. doi: 10.1038/nature13810. [DOI] [PMC free article] [PubMed] [Google Scholar]
10.Seguin-Orlando A., Korneliussen T.S., Sikora M., Malaspinas A.-S., Manica A., Moltke I., Albrechtsen A., Ko A., Margaryan A., Moiseyev V. Paleogenomics. Genomic structure in Europeans dating back at least 36,200 years. Science. 2014;346:1113–1118. doi: 10.1126/science.aaa0114. [DOI] [PubMed] [Google Scholar]
11.Green R.E., Malaspinas A.S., Krause J., Briggs A.W., Johnson P.L.F., Uhler C., Meyer M., Good J.M., Maricic T., Stenzel U. A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing. Cell. 2008;134:416–426. doi: 10.1016/j.cell.2008.06.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
12.Briggs A.W., Good J.M., Green R.E., Krause J., Maricic T., Stenzel U., Lalueza-Fox C., Rudan P., Brajković D., Kućan Z. Targeted retrieval and analysis of five Neandertal mtDNA genomes. Science. 2009;325:318–321. doi: 10.1126/science.1174462. [DOI] [PubMed] [Google Scholar]
13.Wood R.E., Higham T.F.G., De Torres T., Tisnérat-Laborde N., Valladas H., Ortiz J.E., Lalueza-Fox C., Sánchez-Moral S., Cañaveras J.C., Rosas A. A new date for the neanderthals from el Sidrón cave (Asturias, northern Spain) Archaeometry. 2013;55:148–158. [Google Scholar]
14.Castellano S., Parra G., Sánchez-Quinto F.A., Racimo F., Kuhlwilm M., Kircher M., Sawyer S., Fu Q., Heinze A., Nickel B. Patterns of coding variation in the complete exomes of three Neandertals. Proc. Natl. Acad. Sci. USA. 2014;111:6666–6671. doi: 10.1073/pnas.1405138111. [DOI] [PMC free article] [PubMed] [Google Scholar]
15.Karmin M., Saag L., Vicente M., Wilson Sayres M.A., Järve M., Talas U.G., Rootsi S., Ilumäe A.-M., Mägi R., Mitt M. A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 2015;25:459–466. doi: 10.1101/gr.186684.114. [DOI] [PMC free article] [PubMed] [Google Scholar]
16.Mendez F.L., Krahn T., Schrack B., Krahn A.M., Veeramah K.R., Woerner A.E., Fomine F.L.M., Bradman N., Thomas M.G., Karafet T.M., Hammer M.F. An African American paternal lineage adds an extremely ancient root to the human Y chromosome phylogenetic tree. Am. J. Hum. Genet. 2013;92:454–459. doi: 10.1016/j.ajhg.2013.02.002. [DOI] [PMC free article] [PubMed] [Google Scholar]
17.Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27:2987–2993. doi: 10.1093/bioinformatics/btr509. [DOI] [PMC free article] [PubMed] [Google Scholar]
18.Kent, W.J. (2013). http://hgdownload.soe.ucsc.edu/goldenPath/hg19/vsPanTro4/axtNet.
19.Poznik G.D., Henn B.M., Yee M.-C., Sliwerska E., Euskirchen G.M., Lin A.A., Snyder M., Quintana-Murci L., Kidd J.M., Underhill P.A., Bustamante C.D. Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females. Science. 2013;341:562–565. doi: 10.1126/science.1237619. [DOI] [PMC free article] [PubMed] [Google Scholar]
20.Skaletsky H., Kuroda-Kawaguchi T., Minx P.J., Cordum H.S., Hillier L., Brown L.G., Repping S., Pyntikova T., Ali J., Bieri T. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003;423:825–837. doi: 10.1038/nature01722. [DOI] [PubMed] [Google Scholar]
21.Rasmussen M., Anzick S.L., Waters M.R., Skoglund P., DeGiorgio M., Stafford T.W., Jr., Rasmussen S., Moltke I., Albrechtsen A., Doyle S.M. The genome of a Late Pleistocene human from a Clovis burial site in western Montana. Nature. 2014;506:225–229. doi: 10.1038/nature13025. [DOI] [PMC free article] [PubMed] [Google Scholar]
22.Adzhubei I., Jordan D.M., Sunyaev S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013;Chapter 7:20. doi: 10.1002/0471142905.hg0720s76. [DOI] [PMC free article] [PubMed] [Google Scholar]
23.Schwarz J.M., Cooper D.N., Schuelke M., Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat. Methods. 2014;11:361–362. doi: 10.1038/nmeth.2890. [DOI] [PubMed] [Google Scholar]
24.Nordborg M. On the probability of Neanderthal ancestry. Am. J. Hum. Genet. 1998;63:1237–1240. doi: 10.1086/302052. [DOI] [PMC free article] [PubMed] [Google Scholar]
25.Williams N.A., Close J.P., Giouzeli M., Crow T.J. Accelerated evolution of Protocadherin11X/Y: a candidate gene-pair for cerebral asymmetry and language. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 2006;141B:623–633. doi: 10.1002/ajmg.b.30357. [DOI] [PubMed] [Google Scholar]
26.Lee K.H., Song G.J., Kang I.S., Kim S.W., Paick J.S., Chung C.H., Rhee K. Ubiquitin-specific protease activity of USP9Y, a male infertility gene on the Y chromosome. Reprod. Fertil. Dev. 2003;15:129–133. doi: 10.1071/rd03002. [DOI] [PubMed] [Google Scholar]
27.Tyler-Smith C., Krausz C. The will-o’-the-wisp of genetics--hunting for the azoospermia factor gene. N. Engl. J. Med. 2009;360:925–927. doi: 10.1056/NEJMe0900301. [DOI] [PMC free article] [PubMed] [Google Scholar]
28.Wong H.Y., Wang G.M., Croessmann S., Zabransky D.J., Chu D., Garay J.P., Cidado J., Cochran R.L., Beaver J.A., Aggarwal A. TMSB4Y is a candidate tumor suppressor on the Y chromosome and is deleted in male breast cancer. Oncotarget. 2015;6:44927–44940. doi: 10.18632/oncotarget.6743. [DOI] [PMC free article] [PubMed] [Google Scholar]
29.Li N., Dhar S.S., Chen T.-Y., Kan P.-Y., Wei Y., Kim J.-H., Chan C.-H., Lin H.-K., Hung M.-C., Lee M.G. JARID1D is a suppressor and prognostic marker of prostate cancer invasion and metastasis. Cancer Res. 2016;76:831–843. doi: 10.1158/0008-5472.CAN-15-0906. [DOI] [PMC free article] [PubMed] [Google Scholar]
30.Nielsen H.S. Secondary recurrent miscarriage and H-Y immunity. Hum. Reprod. Update. 2011;17:558–574. doi: 10.1093/humupd/dmr005. [DOI] [PubMed] [Google Scholar]
31.Bogaert A.F., Skorska M. Sexual orientation, fraternal birth order, and the maternal immune hypothesis: a review. Front. Neuroendocrinol. 2011;32:247–254. doi: 10.1016/j.yfrne.2011.02.004. [DOI] [PubMed] [Google Scholar]
32.Wang W., Meadows L.R., den Haan J.M., Sherman N.E., Chen Y., Blokland E., Shabanowitz J., Agulnik A.I., Hendrickson R.C., Bishop C.E. Human H-Y: a male-specific histocompatibility antigen derived from the SMCY protein. Science. 1995;269:1588–1590. doi: 10.1126/science.7667640. [DOI] [PubMed] [Google Scholar]
33.Haldane J.B.S. Sex ratio and unisexual sterility in hybrid animals. J. Genet. 1922;12:101–109. [Google Scholar]
34.Auton A., Brooks L.D., Durbin R.M., Garrison E.P., Kang H.M., Korbel J.O., Marchini J.L., McCarthy S., McVean G.A., Abecasis G.R., The 1000 Genomes Project Consortium A global reference for human genetic variation. Nature. 2015;526:68–74. doi: 10.1038/nature15393. [DOI] [PMC free article] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Document S1. Figures S1 and S2

mmc1.pdf^{(109.7KB, pdf)}

Data S1. Regions passing Filter 1

mmc2.zip^{(13.2KB, zip)}

Data S2. Regions passing Filter 2

mmc3.zip^{(10.8KB, zip)}

Data S3. Regions passing Filter 3

mmc4.zip^{(9.1KB, zip)}

Data S4. Regions passing Filter 4

mmc5.zip^{(8.4KB, zip)}

Data S5. Mutations variable among Neanderthal, A00, and the reference sequence

mmc6.zip^{(1.7KB, zip)}

Data S6. Alignment (SAM) file with Neanderthal reads overlapping mutations in Table S3

mmc7.zip^{(30.3KB, zip)}

Data S7. Neanderthal pileup lines for positions indicated in Table S3

mmc8.zip^{(1.4KB, zip)}

Document S2. Spanish Translation: La Divergencia de Los Cromosomas Y de Neandertal y de los Humanos Modernos

mmc9.pdf^{(1MB, pdf)}

Document S3. Article plus Supplemental Data

mmc10.pdf^{(433.9KB, pdf)}

[bib1] 1.McDougall I., Brown F.H., Fleagle J.G. Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature. 2005;433:733–736. doi: 10.1038/nature03258. [DOI] [PubMed] [Google Scholar]

[bib2] 2.Green R.E., Krause J., Briggs A.W., Maricic T., Stenzel U., Kircher M., Patterson N., Li H., Zhai W., Fritz M.H.-Y. A draft sequence of the Neandertal genome. Science. 2010;328:710–722. doi: 10.1126/science.1188021. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib3] 3.Reich D., Green R.E., Kircher M., Krause J., Patterson N., Durand E.Y., Viola B., Briggs A.W., Stenzel U., Johnson P.L.F. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature. 2010;468:1053–1060. doi: 10.1038/nature09710. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib4] 4.Meyer M., Kircher M., Gansauge M.-T., Li H., Racimo F., Mallick S., Schraiber J.G., Jay F., Prüfer K., de Filippo C. A high-coverage genome sequence from an archaic Denisovan individual. Science. 2012;338:222–226. doi: 10.1126/science.1224344. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib5] 5.Prüfer K., Racimo F., Patterson N., Jay F., Sankararaman S., Sawyer S., Heinze A., Renaud G., Sudmant P.H., de Filippo C. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014;505:43–49. doi: 10.1038/nature12886. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib6] 6.Sankararaman S., Mallick S., Dannemann M., Prüfer K., Kelso J., Pääbo S., Patterson N., Reich D. The genomic landscape of Neanderthal ancestry in present-day humans. Nature. 2014;507:354–357. doi: 10.1038/nature12961. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib7] 7.Vernot B., Akey J.M. Resurrecting surviving Neandertal lineages from modern human genomes. Science. 2014;343:1017–1021. doi: 10.1126/science.1245938. [DOI] [PubMed] [Google Scholar]

[bib8] 8.Sankararaman S., Patterson N., Li H., Pääbo S., Reich D. The date of interbreeding between Neandertals and modern humans. PLoS Genet. 2012;8:e1002947. doi: 10.1371/journal.pgen.1002947. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib9] 9.Fu Q., Li H., Moorjani P., Jay F., Slepchenko S.M., Bondarev A.A., Johnson P.L.F., Aximu-Petri A., Prüfer K., de Filippo C. Genome sequence of a 45,000-year-old modern human from western Siberia. Nature. 2014;514:445–449. doi: 10.1038/nature13810. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib10] 10.Seguin-Orlando A., Korneliussen T.S., Sikora M., Malaspinas A.-S., Manica A., Moltke I., Albrechtsen A., Ko A., Margaryan A., Moiseyev V. Paleogenomics. Genomic structure in Europeans dating back at least 36,200 years. Science. 2014;346:1113–1118. doi: 10.1126/science.aaa0114. [DOI] [PubMed] [Google Scholar]

[bib11] 11.Green R.E., Malaspinas A.S., Krause J., Briggs A.W., Johnson P.L.F., Uhler C., Meyer M., Good J.M., Maricic T., Stenzel U. A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing. Cell. 2008;134:416–426. doi: 10.1016/j.cell.2008.06.021. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib12] 12.Briggs A.W., Good J.M., Green R.E., Krause J., Maricic T., Stenzel U., Lalueza-Fox C., Rudan P., Brajković D., Kućan Z. Targeted retrieval and analysis of five Neandertal mtDNA genomes. Science. 2009;325:318–321. doi: 10.1126/science.1174462. [DOI] [PubMed] [Google Scholar]

[bib13] 13.Wood R.E., Higham T.F.G., De Torres T., Tisnérat-Laborde N., Valladas H., Ortiz J.E., Lalueza-Fox C., Sánchez-Moral S., Cañaveras J.C., Rosas A. A new date for the neanderthals from el Sidrón cave (Asturias, northern Spain) Archaeometry. 2013;55:148–158. [Google Scholar]

[bib14] 14.Castellano S., Parra G., Sánchez-Quinto F.A., Racimo F., Kuhlwilm M., Kircher M., Sawyer S., Fu Q., Heinze A., Nickel B. Patterns of coding variation in the complete exomes of three Neandertals. Proc. Natl. Acad. Sci. USA. 2014;111:6666–6671. doi: 10.1073/pnas.1405138111. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib15] 15.Karmin M., Saag L., Vicente M., Wilson Sayres M.A., Järve M., Talas U.G., Rootsi S., Ilumäe A.-M., Mägi R., Mitt M. A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 2015;25:459–466. doi: 10.1101/gr.186684.114. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib16] 16.Mendez F.L., Krahn T., Schrack B., Krahn A.M., Veeramah K.R., Woerner A.E., Fomine F.L.M., Bradman N., Thomas M.G., Karafet T.M., Hammer M.F. An African American paternal lineage adds an extremely ancient root to the human Y chromosome phylogenetic tree. Am. J. Hum. Genet. 2013;92:454–459. doi: 10.1016/j.ajhg.2013.02.002. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib17] 17.Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27:2987–2993. doi: 10.1093/bioinformatics/btr509. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib18] 18.Kent, W.J. (2013). http://hgdownload.soe.ucsc.edu/goldenPath/hg19/vsPanTro4/axtNet.

[bib19] 19.Poznik G.D., Henn B.M., Yee M.-C., Sliwerska E., Euskirchen G.M., Lin A.A., Snyder M., Quintana-Murci L., Kidd J.M., Underhill P.A., Bustamante C.D. Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females. Science. 2013;341:562–565. doi: 10.1126/science.1237619. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib20] 20.Skaletsky H., Kuroda-Kawaguchi T., Minx P.J., Cordum H.S., Hillier L., Brown L.G., Repping S., Pyntikova T., Ali J., Bieri T. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003;423:825–837. doi: 10.1038/nature01722. [DOI] [PubMed] [Google Scholar]

[bib21] 21.Rasmussen M., Anzick S.L., Waters M.R., Skoglund P., DeGiorgio M., Stafford T.W., Jr., Rasmussen S., Moltke I., Albrechtsen A., Doyle S.M. The genome of a Late Pleistocene human from a Clovis burial site in western Montana. Nature. 2014;506:225–229. doi: 10.1038/nature13025. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib22] 22.Adzhubei I., Jordan D.M., Sunyaev S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013;Chapter 7:20. doi: 10.1002/0471142905.hg0720s76. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib23] 23.Schwarz J.M., Cooper D.N., Schuelke M., Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat. Methods. 2014;11:361–362. doi: 10.1038/nmeth.2890. [DOI] [PubMed] [Google Scholar]

[bib24] 24.Nordborg M. On the probability of Neanderthal ancestry. Am. J. Hum. Genet. 1998;63:1237–1240. doi: 10.1086/302052. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib25] 25.Williams N.A., Close J.P., Giouzeli M., Crow T.J. Accelerated evolution of Protocadherin11X/Y: a candidate gene-pair for cerebral asymmetry and language. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 2006;141B:623–633. doi: 10.1002/ajmg.b.30357. [DOI] [PubMed] [Google Scholar]

[bib26] 26.Lee K.H., Song G.J., Kang I.S., Kim S.W., Paick J.S., Chung C.H., Rhee K. Ubiquitin-specific protease activity of USP9Y, a male infertility gene on the Y chromosome. Reprod. Fertil. Dev. 2003;15:129–133. doi: 10.1071/rd03002. [DOI] [PubMed] [Google Scholar]

[bib27] 27.Tyler-Smith C., Krausz C. The will-o’-the-wisp of genetics--hunting for the azoospermia factor gene. N. Engl. J. Med. 2009;360:925–927. doi: 10.1056/NEJMe0900301. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib28] 28.Wong H.Y., Wang G.M., Croessmann S., Zabransky D.J., Chu D., Garay J.P., Cidado J., Cochran R.L., Beaver J.A., Aggarwal A. TMSB4Y is a candidate tumor suppressor on the Y chromosome and is deleted in male breast cancer. Oncotarget. 2015;6:44927–44940. doi: 10.18632/oncotarget.6743. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib29] 29.Li N., Dhar S.S., Chen T.-Y., Kan P.-Y., Wei Y., Kim J.-H., Chan C.-H., Lin H.-K., Hung M.-C., Lee M.G. JARID1D is a suppressor and prognostic marker of prostate cancer invasion and metastasis. Cancer Res. 2016;76:831–843. doi: 10.1158/0008-5472.CAN-15-0906. [DOI] [PMC free article] [PubMed] [Google Scholar]

[bib30] 30.Nielsen H.S. Secondary recurrent miscarriage and H-Y immunity. Hum. Reprod. Update. 2011;17:558–574. doi: 10.1093/humupd/dmr005. [DOI] [PubMed] [Google Scholar]

[bib31] 31.Bogaert A.F., Skorska M. Sexual orientation, fraternal birth order, and the maternal immune hypothesis: a review. Front. Neuroendocrinol. 2011;32:247–254. doi: 10.1016/j.yfrne.2011.02.004. [DOI] [PubMed] [Google Scholar]

[bib32] 32.Wang W., Meadows L.R., den Haan J.M., Sherman N.E., Chen Y., Blokland E., Shabanowitz J., Agulnik A.I., Hendrickson R.C., Bishop C.E. Human H-Y: a male-specific histocompatibility antigen derived from the SMCY protein. Science. 1995;269:1588–1590. doi: 10.1126/science.7667640. [DOI] [PubMed] [Google Scholar]

[bib33] 33.Haldane J.B.S. Sex ratio and unisexual sterility in hybrid animals. J. Genet. 1922;12:101–109. [Google Scholar]

[bib34] 34.Auton A., Brooks L.D., Durbin R.M., Garrison E.P., Kang H.M., Korbel J.O., Marchini J.L., McCarthy S., McVean G.A., Abecasis G.R., The 1000 Genomes Project Consortium A global reference for human genetic variation. Nature. 2015;526:68–74. doi: 10.1038/nature15393. [DOI] [PMC free article] [PubMed] [Google Scholar]

PERMALINK

The Divergence of Neandertal and Modern Human Y Chromosomes

Fernando L Mendez

G David Poznik

Sergi Castellano

Carlos D Bustamante

Abstract

Introduction