Human paternal and maternal demographic histories: insights from high-resolution Y chromosome and mtDNA sequences
© Lippold et al.; licensee BioMed Central Ltd. 2014
Received: 14 August 2014
Accepted: 22 August 2014
Published: 24 September 2014
Comparisons of maternally-inherited mitochondrial DNA (mtDNA) and paternally-inherited non-recombining Y chromosome (NRY) variation have provided important insights into the impact of sex-biased processes (such as migration, residence pattern, and so on) on human genetic variation. However, such comparisons have been limited by the different molecular methods typically used to assay mtDNA and NRY variation (for example, sequencing hypervariable segments of the control region for mtDNA vs. genotyping SNPs and/or STR loci for the NRY). Here, we report a simple capture array method to enrich Illumina sequencing libraries for approximately 500 kb of NRY sequence, which we use to generate NRY sequences from 623 males from 51 populations in the CEPH Human Genome Diversity Panel (HGDP). We also obtained complete mtDNA genome sequences from the same individuals, allowing us to compare maternal and paternal histories free of any ascertainment bias.
We identified 2,228 SNPs in the NRY sequences and 2,163 SNPs in the mtDNA sequences. Our results confirm the controversial assertion that genetic differences between human populations on a global scale are bigger for the NRY than for mtDNA, although the differences are not as large as previously suggested. More importantly, we find substantial regional variation in patterns of mtDNA versus NRY variation. Model-based simulations indicate very small ancestral effective population sizes (<100) for the out-of-Africa migration as well as for many human populations. We also find that the ratio of female effective population size to male effective population size (Nf/Nm) has been greater than one throughout the history of modern humans, and has recently increased due to faster growth in Nf than Nm.
The NRY and mtDNA sequences provide new insights into the paternal and maternal histories of human populations, and the methods we introduce here should be widely applicable for further such studies.
KeywordsPopulation genetics Population size Simulations HGDP
Comparisons of mtDNA and NRY variation have provided numerous important insights into the maternal and paternal histories of human populations [1–3]. However, such comparisons are limited by methodological differences in how mtDNA and NRY variation have been typically assayed. MtDNA variation is usually investigated by sequencing hypervariable segments of the control region, (or, increasingly, via complete mtDNA genome sequences), while human NRY variation is routinely assayed by genotyping SNPs of interest, often in combination with short tandem repeat (STR) loci. Nevertheless, NRY SNP typing has several drawbacks due to the ascertainment bias inherent in the selection of SNPs [1, 4, 5]. This ascertainment bias complicates many analyses of interest, such as dating the age of the NRY ancestor or particular divergence events in the NRY phylogeny, as well as demographic inferences such as population size changes . Moreover, the difference in molecular methods used to assay NRY versus mtDNA variation can complicate the interpretation of differences between patterns of NRY and mtDNA variation. For example, the seminal finding that NRY differences are bigger than mtDNA differences among global populations of humans, and that this is due to a higher rate of female than male migration due to patrilocality , may instead reflect methodological differences in how mtDNA versus NRY variation was assayed in that study .
Another fundamental question concerns whether or not male and female effective population sizes have been the same over time. Attempts to address this question using the ratio of X chromosome to autosomal DNA diversity have come up with conflicting answers [8, 9], which may in part reflect the use of different methods that capture information about effective population size at different times in the past . Moreover, the ratio of X to autosome diversity varies along the X chromosome, depending how far polymorphic sites are from genes [11–13], indicating a potential role for selection in distorting effective population size estimates from comparisons of X chromosome to autosomal DNA diversity. These and other fundamental aspects of human maternal and paternal demographic history remain unanswered.
Recently, analyses have been carried out of NRY sequences obtained as part of whole genome sequencing projects [14–16]. While these studies provide very detailed insights into the NRY phylogeny, they are nonetheless limited by the expense of whole genome sequencing, which precludes comprehensive global sampling. To allow for more accurate comparisons between mtDNA and NRY variation and to permit demographic inferences based on the NRY, we developed a capture-based array to enrich Illumina sequencing libraries for approximately 500 kb of NRY sequence. We used this approach to obtain NRY sequences from 623 males from 51 globally-distributed populations, and we also obtained complete mtDNA genome sequences from the same individuals, allowing us to investigate and directly compare the paternal and maternal relationships of global human populations in unprecedented detail.
Samples and sequencing library preparation
The samples consist of 623 males (Additional file 1: Table S1) from the CEPH Human Genome Diversity Panel (HGDP) . The samples were taken from the subset ‘H952’, which excludes atypical, duplicated, and closely-related samples . Approximately 200 ng of genomic DNA from each sample was sheared by sonication using a Bioruptor system (Diogenode) and used to construct an Illumina Sequencing library with a specific double-index as described previously . The libraries were then enriched separately for NRY and mtDNA sequences as described below.
Y-chromosome capture array design
We targeted unique regions on the NRY that are free of repeats and to which the typically short next-generation sequencing reads could be mapped with high confidence. We used the UCSC table browser  and the February 2009 (GRCh37/hg19) assembly and applied the following filter criteria. First, from the group ‘variation and repeats’, sequence regions annotated in the following tracks were removed: Interrupted Repeats, RepeatMasker, Simple Repeats, and Segmental Duplications. Next, we used the ‘mapability’ table ‘CRG Align 75’ from the group ‘mapping and sequencing tracks’ to identify and remove regions with mapability scores below 1. We then removed regions of less than 500 bp in order to reduce the number of fragments and thereby the number of fragment ends, which have low probe densities. We also removed 15mers that occurred more than 100 times in the hg19 genome assembly, as described previously , which resulted in splitting some target regions into sub-regions that were less than 500 bp. The final result was a total of approximately 500 kb of unique NRY sequence, distributed among 655 target regions ranging from 61 bp to 3.9 kb (Additional file 2: Table S2). These regions were then used to design a custom array (SureSelect 1 M capture array, Agilent) with 60 nt probes that were printed twice with a tiling density of 1 bp.
Up to 60 barcoded libraries were pooled in equimolar ratio. The library mix was enriched for target NRY regions by hybridization-capture on the custom designed array following the protocol described previously . After enrichment the library-pool was quantified by qPCR and then amplified to a total of approximately 1012 molecules. The final concentration and length distribution was measured on an Agilent DNA 100 microchip, and 10 nmol of the amplified library pool was used for sequencing. Each pool, consisting of 48 to 60 samples, was sequenced on a Solexa GAII lane using a paired end 75 cycle run plus two 7 nt index reads.
Up to 94 libraries were pooled in equimolar ratio and the library pool was enriched for mtDNA sequences by an in-solution hybridization capture method . The hybridization eluate was measured by qPCR and then amplified to produce a final concentration of 10 nmol. Up to 200 samples were sequenced on a Solexa GAII lane using a paired end 75 cycle run, plus two 7 nt index reads.
In each Solexa GAII lane, 1% PhiX174 phage DNA was spiked in and used as a training set to estimate base quality scores with the IBIS base-caller . Reads with more than five bases having a PHRED scaled quality score below Q15 were discarded, as were reads having a single base quality in the index read (7 nt) score below Q10. Reads with no mismatches to the expected double index sequences were assigned to each individual sample library.
For the NRY-enriched data, reads were mapped to the human reference genome (GRCh37) using default settings with BWA v0.5.10 . We mapped to the whole genome rather than just the target region, in order to identify reads that might, with equal probability, map to another position in the genome. The bam files containing the mapping information and reads were processed with samtools v0.1.18 . We used Picard 1.42 to mark duplicates, based on the start and end coordinates of the read pairs. The final SNP call was done on all samples simultaneously using the UnifiedGenotyper from the GATK v2.0-35 package  and the following options: --output_mode EMIT_ALL_CONFIDENT_SITES, --genotype_likelihoods_model SNP, --min_base_quality_score 20 and --heterozygosity 0.0000000001. The result was stored in a VCF file containing information for each callable site of the target region, and a second VCF file was created that contained only the variable positions among the 623 samples. For each sample at each variable position the PL scores were calculated with samtools ; PL scores are normalized, PHRED-scaled likelihoods for the three genotypes (0/0, 0/1, 1/1) and are based on, among other things, coverage, base quality, and mapping quality. Positions that showed a difference in the PL score of less than 30 between homozygote reference (0/0) and homozygote alternative (1/1) were called an ‘N’ in that sample, as were positions where heterozygote calls (0/1) either had a higher PL score than the most likely homozygous genotype, or differed by less than 30 from the most likely homozygous genotype. Note that a PL score of 30 between genotype 0/0 and 1/1 means that the former is 1,000 times more likely than the latter, for example, the genotype-calling error rate is expected to be less than 1 in 1,000. Sites where more than two bases were called (that is, multi-allelic sites) were also removed.
For the mtDNA-enriched data, reads were mapped to the revised mtDNA reference sequence (GenBank number: NC_012920) using the software MIA . The consensus sequences were aligned using MUSCLE v3.8.31  (cmd line: muscle -maxiters 1 -diags mt_623seq.fasta mt_623seq.aln), and haplogroups were called with the HaploGrep software .
Imputation for the NRY
After quality filtering, there were 2,276 variable sites in the NRY sequences, with a total of 2.54% of the individual genotypes at variable positions scored as ‘N’ (that is, as missing data; the number of missing sites per individual ranged from 9 to 1,173, with an average of 122 missing sites per individual). Since missing data can influence the results of some analyses, we took advantage of the fact that the NRY target regions are completely linked with no recombination to impute missing data as follows. First, all sites with no missing data (605 sites) were used as the reference set to define haplotypes and calculate the number of differences between each haplotype. Sites with missing data were then imputed, beginning with the site with the smallest amount of missing data and proceeding sequentially. For each haplotype with missing data for that site, the missing base was imputed as the allele present in the reference haplotype that had the fewest differences (based on the sites with no missing data). After imputation was finished for that site, it was added to the reference set, and the procedure continued for the next site with the smallest amount of missing data.
As a check on the accuracy of the imputation, we randomly deleted 2.54% of the known alleles, following the distribution of missing alleles in the full dataset, thereby creating an artificial dataset with a similar distribution of missing alleles as in the observed dataset. We then imputed the missing data according to the above procedure and compared the imputed alleles to the true alleles; this procedure was carried out 1,000 times. The imputed allele matched the true allele in 99.1% of the comparisons, indicating that the imputation procedure is quite accurate.
Recurrent NRY mutations
We expect the majority of the NRY SNPs to have mutated only once, as recurrent mutations in the known NRY phylogeny are quite rare [31, 32]. Therefore, as a further quality control measure, we investigated the NRY data for recurrent mutations by constructing a maximum parsimony tree for the 2,276 SNPs using programs in PHYLIP. We then estimated the number of mutations at each SNP, and removed 48 SNPs that had mutated more than twice, and only in terminal branches, as these are likely to reflect sequencing errors. The final dataset contains 2,228 SNPs.
Basic summary statistics (haplotype diversity, mean number of pairwise differences, nucleotide diversity, Tajima’s D value and theta(S)) were calculated using Arlequin v188.8.131.52 . Arlequin was further used to estimate pairwise ΦST values and for Analysis of Molecular Variance (AMOVA). The observed ratio of the mean pairwise differences (mpd) for the NRY versus mtDNA was calculated as mpdNRY/mpdmt. In order to detect group-specific deviations from the mean distribution of the mpd ratio in the dataset, we carried out a resampling approach. For each group sample size (Ngroup) we chose randomly Ngroup individuals (out of 623) and calculated the mpd ratio using the dist.dna command from the APE package  in R. This was repeated 10,000 times for each Ngroup sample size to obtain the distribution of resampled mpd ratios.
Divergence times in the NRY and mtDNA phylogenies were estimated using a Bayesian approach implemented in BEAST v1.6.2 . For the mtDNA genome sequences we divided the alignment into two partitions consisting of the coding and non-coding regions, respectively. For both partitions we estimated the best fitting substitution model using jModeltest  and the mutation rates estimated previously . These rates were calibrated by a combination of chimpanzee-human divergence and archaeological colonization dates, take into account time-dependency in the molecular clock, and are in the range of recent estimates of the mtDNA mutation rate [15, 38, 39]. For the non-coding region we used the GTR + I + G substitution model and a mutation rate of 9.883 × 10−8 substitutions/site/year, while for the coding region we used the TrN + I + G model and a mutation rate of 1.708 × 10−8 substitutions/site/year. A strict clock and a constant size coalescence model were used, and the MCMC was run for 10 million steps with sampling from the posterior every 2,000 steps. The MCMC was run on five independent chains in parallel. After careful inspection of the log files in Tracer, the tree files of the five runs were merged after discarding the first 2,500 trees (50%) of each run as burn-in. A consensus tree was built from the merged trees using TreeAnnotator, and the consensus tree showing the divergence times for each node was visualized with FigTree.
For the NRY sequences the same procedure was used, but modified as only variable sites were included in the BEAST analysis in order to reduce the computational time. The substitution model used was HKY without I + G, and the substitution rate was multiplied by the number of callable sites (501,108 sites) divided by the number of variable sites (2,228 sites). As there is uncertainty regarding the mutation rate, we ran the analysis twice, with a ‘fast’ rate  of 1.00 × 10−9 substitutions/site/year (transformed to 2.25 × 10−7) and with a ‘slow’ rate  of 6.17 × 10−10 substitutions/site/year (transformed to 1.39 × 10−7).
Bayesian skyline plots  were used to estimate population size change through time, using the same mutation rates and substitution models described above. The piecewise-linear Skyline coalescence model was chosen and the number of groups (bins) was set to half the sample size per group with a maximum of 20. A single MCMC chain was ran for 30 million steps and sampled every 3,000 steps from the posterior. The log file was inspected in Tracer for convergence of the chain and ESS values and the Bayesian Skyline Reconstruction was run.
NRY and mtDNA diversity
We obtained approximately 500 kb of NRY sequence from the 623 males in the HGDP, and complete mtDNA genome sequences from these 623 males plus an additional 329 females from the HGDP. The average coverage of the NRY sequences was 14.5X (range, 5X-37.5X, Additional file 3: Figure S1), while for the mtDNA genome sequences the average coverage was 640X (range, 46X-4123X, Additional file 3: Figure S1). After quality-filtering, imputation, and removal of sites with a high number of recurrent mutations, there remained 2,228 SNPs in the NRY sequences. The mtDNA analyses here are restricted to the 623 males for which NRY sequences were obtained, for which there were 2,163 SNPs; results based on the mtDNA genome sequences from the entire set of HGDP samples (952 individuals) did not differ from those based on the subset of 623 males (for example, Additional file 3: Figure S2). More details about the results from each individual, including mtDNA and NRY haplogroups, are provided in Additional file 1: Table S1. The mtDNA sequences have been deposited in Genbank with accession numbers KF450814-KF451871. A datafile with the alleles at each of the NRY SNPs in each sample has been provided to the CEPH-HGDP and additionally is available from the authors. The NRY raw sequencing data are in the European Nucleotide Archive with the study accession number PRJEB4417 (sample accession numbers ERS333252-ERS333873).
Basic summary statistics for the mtDNA and NRY diversity in each population are provided in Additional file 3: Table S3. As the sample sizes for many of the individual populations are quite small, for most subsequent analyses we grouped the populations into the following regions (based on analyses of genome-wide SNP data [43, 47]): Africa, America, Central Asia, East Asia, Europe, Middle East/North Africa (ME/NA), and Oceania (the regional affiliation for each population is in Additional file 1: Table S1). The Adygei, Hazara, and Uygur were excluded from these groupings as they show evidence of substantial admixture between these regional groups [43, 47]. We stress that the use of regional names is a convenience to refer to these groupings of these specific populations, and should not be taken to represent the entirety of the regions (for example, ‘Africa’ refers to the results based on the analysis of the combined African HGDP samples, not to Africa in general).
Summary statistics for regional groups
mpd ± SE
π ± SE a
mpd ± SE
π ± SE b
41.0 ± 18.0
80 ± 40
78.3 ± 34.0
47 ± 23
32.1 ± 14.1
62 ± 31
42.4 ± 18.5
26 ± 12
35.0 ± 15.3
71 ± 36
42.3 ± 18.5
26 ± 12
42.7 ± 18.7
85 ± 40
42.0 ± 18.4
25 ± 12
30.0 ± 13.2
58 ± 31
29.3 ± 12.9
18 ± 9
34.7 ± 15.9
71 ± 36
41.9 ± 19.2
25 ± 13
11.8 ± 5.5
22 ± 13
34.9 ± 15.8
21 ± 11
NRY and mtDNA population differentiation
An outstanding question is whether or not there are differences in the relative amounts of between-population versus within-population diversity for mtDNA versus the NRY, as some studies have found much larger between-population differences for the NRY than for mtDNA  while others have not . To address this question, we carried out an AMOVA; the results (Figure 2B) show that in the entire worldwide dataset, the between-population differences are indeed bigger for the NRY (approximately 36% of the variance) than for mtDNA (approximately 25% of the variance). However, there are substantial differences among the regional groups. The ME/NA, East Asia, and Europe regional groups follow the worldwide pattern in having bigger between-population differences for the NRY than for mtDNA. In contrast, Africa, Oceania, and the Americas have substantially bigger between-population differences for mtDNA than for the NRY, while for central Asia the between-population variation is virtually identical for the NRY and mtDNA. These regional differences likely reflect the influence of sex-biased migrations and admixture, as discussed in more detail below, and moreover indicate that focusing exclusively on the worldwide pattern of mtDNA versus NRY variation misses these important regional differences.
We also investigated the relationship between geography and genetic distance. Despite the small sample sizes at the population level, both mtDNA and NRY ΦST distances are significantly correlated with geographic distances between populations (Mantel tests with 1,000 replications: mtDNA, r = 0.41, P <0.001; NRY, r = 0.36, P = 0.002) as well as with each other (r = 0.23, P = 0.025). Thus, NRY and mtDNA divergence are both highly associated with geographic distances among populations.
MtDNA and NRY phylogenies
The age of the mtDNA ancestor is estimated to be about 160 thousand years ago (kya), and the ages of the non-African mtDNA lineages M and N are about 65 to 70 kya, in good agreement with previous estimates . Our estimate for the age of the NRY ancestor is 103 kya based on the fast rate, and 165 kya based on the slow rate; however these estimates do not include the recently-discovered ‘A00’ lineage , which would result in much older ages for the NRY ancestor. The close agreement between the slow NRY ancestor age (165 kya) and the mtDNA ancestor age (160 kya) might be taken as evidence in favor of the slow NRY mutation rate. However, the slow NRY mutation rate gives an estimated age for the initial out-of-Africa divergence of about 100 kya, and an age for the divergence of Amerindian-specific haplogroup Q lineages of about 20 kya, while the fast rate gives corresponding estimates of about 60 kya for out-of-Africa and about 12.5 kya for Amerindian haplogroup Q lineages, in better agreement with the mtDNA and other evidence for these events [54–57]. Given the current uncertainty over mutation rate estimates, we have chosen to use either both estimates in further analyses (for example, Bayesian skyline plots) or an average of the fast and slow rates (for example, in simulation-based analyses); in Additional file 3: Table S4 we provide divergence time estimates and associated 95% credible intervals for the branching events shown in the phylogenies in Figure 3.
NRY and mtDNA haplogroup frequencies per population are shown in Additional file 3: Table S5 and Additional file 3: Table S6, respectively. The mtDNA haplogroups were called from the sequences determined here, while the NRY haplogroups were previously determined by SNP genotyping [58, 59]. The NRY haplogroup information we provide is taken only from these published data; we did not infer haplogroups from the sequences, in order to have an independent comparison of the NRY tree with the haplogroups. The phylogenetic relationships for the NRY sequences are generally concordant with the SNP-genotyping results (with some exceptions, discussed in the legends to Figures S3 to S12 in Additional file 3). The haplogroup frequencies provide further insights into some of the different regional patterns of mtDNA versus NRY diversity noted previously. For example, the comparatively low diversity and smaller differences among populations for the NRY in Africa is due to the high frequency of NRY haplogroup E (55% to 100% in the non-Khoisan groups; Additional file 3: Table S5). This haplogroup is widespread in western Africa, and specific subhaplogroups of haplogroup E are associated with the Bantu expansion [59–61]. The comparatively low NRY diversity in the HGDP Africa regional group thus likely reflects a ‘homogenizing’ effect of the Bantu expansion. NRY haplogroup E is also of interest because it occurs in some European and ME/NA groups, at frequencies of up to 17%, as well as in a few individuals from Central Asia (Additional file 3: Table S5). Inspection of the phylogeny of haplogroup E sequences (Additional file 3: Figure S7) reveals that all of the European and most of the ME/NA haplogroup E sequences form a clade distinct from the African haplogroup E sequences, and the age of this clade is about 18 kya. Moreover, all of the European haplogroup E sequences fall into a subclade that is about 14 kya. These results may reflect a migration from North Africa to Europe suggested from analyses of genome-wide SNP data , and would thus provide a timeframe for this migration.
In Oceania, the bigger differences between populations for mtDNA than for the NRY (Figure 2B, Table 1) probably reflect the high frequency of mtDNA haplogroup B in just one of the two Oceania populations (75% in the Melanesian population vs. 0% in the Papuan population; Additional file 3: Table S6). MtDNA haplogroup B is associated with the Austronesian expansion [63–65]. By contrast, NRY haplogroups associated with the Austronesian expansion, such as haplogroup O [63, 66, 67] are absent in the HGDP Oceania populations (Additional file 3: Table S5). This contrast further testifies to the larger maternal than paternal impact of the Austronesian expansion on Oceanian populations [63, 66–69].
In the Americas, there are dramatic differences in mtDNA haplogroup frequencies among populations (the Karitiana and Surui are 100% haplogroup D, the Pima are 100% haplogroup C, the Maya are 100% haplogroup A, and the Colombians are 50% haplogroup B and 50% haplogroup C; Additional file 3: Table S6), which are at least partly due to the small sample sizes but also in keeping with previous studies . However, all NRY sequences from the Americas fall into haplogroup Q (with the exception of one Pima with a haplogroup G sequence that likely reflects recent European admixture), and overall NRY diversity is substantially reduced in the Americas, compared to mtDNA diversity (Table 1, Figure 2). While the small number of HGDP males from the Americas precludes any definitive statements, the apparently much greater mtDNA than NRY diversity in the Americas might indicate that fewer males than females were involved in the colonization of the Americas, and deserves further investigation.
We note some additional features pertaining to specific populations in the individual NRY haplogroup phylogenies provided in Figures S3 to S12 in Additional file 3, while the full mtDNA phylogeny for the HGDP samples is provided in Figure S13 in Additional file 3.
Prior estimates of divergence time (all priors uniformly distributed) and the mean, mode, and 95% HPD (highest posterior density) intervals
Coverage is the proportion of times the true value for the parameter lies within the 90% credible interval around the parameter estimate; and Factor 2 is the proportion of estimated values that are within 50% and 200% of the true value.
Current and ancestral estimates of male effective population size (N m ) based on simulations of the HGDP NRY sequences
Current and ancestral estimates of female effective population size (N f ) based on simulations of the HGDP mtDNA sequences
We report here the development and implementation of a capture-based array method to enrich Illumina sequencing libraries for NRY sequences. We then used this method to obtain approximately 500 kb of NRY sequence for 623 males from 51 populations of the CEPH-HGDP, and we also obtained complete mtDNA genome sequences from the same individuals. The molecular resolution (that is, number of SNPs) provided by the NRY and mtDNA sequences was roughly equivalent overall (2,228 NRY SNPs, vs. 2,163 mtDNA SNPs), allowing us to compare the maternal and paternal histories of human populations without the usual concerns about different methodologies (for example, mtDNA HV1 sequences vs. genotyping NRY SNPs and/or STRs) having an influence on the results. However, note that in other respects the molecular resolution still differs between the mtDNA and NRY sequences, for example, we obtained complete mtDNA genome sequences but only partial NRY sequences.
Our data provide new insights into the maternal versus paternal history of humans. First, a longstanding controversy has been whether or not genetic differences between human populations are bigger, on average, for the NRY than for mtDNA. The first comparative study of human mtDNA and NRY diversity found significantly bigger differences between populations for the NRY than for mtDNA , which was attributed to a higher female than male migration rate between populations due to patrilocality. A subsequent study found bigger differences between populations for the NRY than for mtDNA in patrilocal populations and the opposite pattern in matrilocal populations,  which was viewed as supporting an influence of residence pattern and associated male versus female migration rates on NRY versus mtDNA diversity. However, these studies used different methods to assay NRY versus mtDNA diversity; notably, a later study that used similar methods to assay NRY and mtDNA diversity (by comparing 6.7 kb of NRY sequence and 770 bp of mtDNA sequence in 389 individuals from 10 populations) concluded that genetic differences between populations were in fact similar for the NRY and mtDNA .
Our results, based on a more comprehensive sampling of worldwide human populations, indicate that genetic differences among human populations at the global scale are indeed bigger for the NRY than for mtDNA, although the differences are not as large as suggested by previous studies (between-population variance of 36% for the NRY vs. 25% for mtDNA in this study, compared to previous estimates of 65% for the NRY vs. 20% for mtDNA ). More importantly, our results indicate substantial differences among regional groups in the between-group variance for the NRY versus mtDNA (Figure 2) as well as in overall levels of NRY versus mtDNA diversity (Figure 2, Table 1). Thus, focusing on global patterns of variation misses this important regional variation, which (as discussed in more detail above in the Results) likely reflects differences in the paternal versus maternal demographic history of specific human populations (for example, the large impact of the Bantu expansion on African NRY diversity [59, 60], and of the Austronesian expansion on Oceanic mtDNA diversity [63, 65]).
Another question of interest is the extent to which the genetic contributions of males versus females have been the same or differed (as measured by their respective effective population sizes, Nm and Nf, respectively). Previous studies of Nm and Nf have largely relied on comparisons of X chromosome vs. autosomal variation, and have come to varying conclusions concerning the historical Nf/Nm ratio, for example, finding that this ratio suggests a large excess of Nf to Nm , a moderate excess of Nf to Nm , or even a decreased Nf relative to Nm . These differences variously reflect methodological differences, difficulties in accounting for differences in male versus female mutation rates, and/or the potentially greater effect of selection on the X chromosome than on the autosomes [10, 11]. Comparison of mtDNA versus NRY variation offers a more direct assessment of Nf/Nm that is free of some of the issues concerning X:autosome comparisons (albeit not all, as discussed below), but requires unbiased estimates of NRY variation, which until our study were only available from either whole genome sequencing studies [5, 14–16] or more limited targeted studies of NRY sequence variation [7, 77]. Our results indicate a consistent strong excess of Nf versus Nm starting even before the out-of-Africa migration that has been carried through almost all subsequent migrations. East Asia may be an exception, and indeed our estimates of Nf and Nm are substantially larger than previous estimates of Ne in east Asians based on autosomal diversity [78, 79]. However, these previous studies were based solely on data from Han Chinese and Japanese, whereas the HGDP includes a much more diverse sampling of east Asian populations, which may account for the higher effective population size estimates for the HGDP. The excess of Nfversus Nm become even more pronounced in recent times due to higher rates of growth in Nf than in Nm (Figures 4, 5, and 6); these results are in line with previous studies of smaller datasets that used different methods [4, 80]. These results suggest, in turn, that sex-specific processes that reduce Nm, such as polygyny and/or sex-specific migration , have characterized humans over most of our prehistory.
However, there are several reasons why this conclusion should be viewed as tentative. First, the sample sizes of some of the regional groups in the HGDP are quite low, precluding confident estimates of effective population sizes. Moreover, there are some surprising features of our results, such as the much larger effective size estimates for East and Central Asians than for Europeans. Whether these features are truly indicative of these regions, or rather specific to the particular populations sampled in the HGDP, will require further studies to elucidate. Nonetheless, given that the HGDP overall is a much more comprehensive sampling of worldwide genetic diversity than in previous studies that estimated effective population sizes for various human populations, it perhaps is not surprising that we obtain different results.
Second, while focusing on NRY versus mtDNA variation avoids some of the drawbacks of comparing X versus autosomal DNA variation in estimating Nf and Nm, the uncertainty associated with the resulting estimates is significantly larger for NRY:mtDNA than for X:autosome comparisons. This is because the X:autosome comparisons are averaged across many independent loci, whereas the NRY and mtDNA are each just a single independent locus.
Third, the model used in the simulations is obviously a very simplified version of reality, and indeed there are some clear differences between the observed values for some summary statistics and the posteriors (for example, the ΦST values in Figures S14 to S16 in Additional file 3). In particular, to reduce the computational complexity we did not consider migration between regional groups (after the initial colonization events) in the simulations to estimate Nf and Nm. There is some justification for doing so, as in general migration within the regional groups has been more important than migration between regional groups, as evidenced by genetic structure analyses [43, 47, 81] and by attempts to estimate migration rates directly from genetic data . Moreover, no mtDNA sequences are shared between regional groups, and only one NRY sequence is shared between regional groups, suggesting very limited recent migration between regional groups. Furthermore, by not including migration we are overestimating the ancestral Nf and Nm (because some of the diversity reflects later migration rather than genetic diversity that was present in the ancestral population). Thus, the effect of such migration would be even smaller estimates of Nf and Nm than those we obtained. Still, in future analyses migration and other complexities should be considered.
Fourth, we have here interpreted differences in levels of NRY versus mtDNA diversity and divergence as reflecting neutral, demographic history. However, a recent study has shown that background selection on the Y chromosome is probably also influencing levels of NRY diversity in human populations . The substantial regional variation that we see in comparisons of mtDNA versus NRY diversity does suggest that there are regional differences in the demographic history of males and females, as it seems unlikely that levels of background selection would vary so drastically across human populations. Moreover, recently-described regional variation in ratios of X:autosomal diversity also point to sex-biased demographic processes . Still, the overall differences we find in Nmversus Nf may be influenced by background selection, and hence may not be as large as inferred by the simulations (for example, Figure 6). More detailed investigations are warranted into the relative importance of background selection versus purely demographic processes in influencing regional variation in Nm vs. Nf.
We have developed a rapid and cost-effective means of obtaining unbiased, high-resolution NRY sequence information. Comparative analysis of NRY and mtDNA sequences from a large sample of individuals and populations from the HGDP provides new insights into the comparative demographic history of males and females. In particular, we find on average larger genetic differences between populations for the NRY than for mtDNA (albeit with substantial regional variation), and that the effective population size of females has been larger than that of males throughout human history. We anticipate that using this approach to investigate additional populations should provide a rich source of new information about the genetic history of our species.
Approximate Bayesian computation
Analysis of Molecular Variance
Analyses of Phylogenetics and Evolution
Binary alignment map
Bayesian Evolutionary Analysis Sampling Trees
Bayesian Skyline Plot
Effective sample size
Genome Analysis Toolkit
Human Genetic Diversity Panel
Improved Base Identification System
Thousand years ago
Markov Chain Monte Carlo
Middle East and North Africa
Mapping Iterative Assembler
Mean number of pairwise differences
Multiple Sequence Comparison by Log Expectation
Female effective population size
Group sample size
Male effective population size
Non-recombining Y chromosome
Phylogeny Inference Package
Quantitative polymerase chain reaction
Relative mean square error
Short tandem repeat
University of California Santa Cruz
Variant Call Format.
We thank: Maanasa Raghavan and Martina Codina for technical assistance; the MPI-EVA Evolutionary Genetics Sequencing Group and the Bioinformatics Group for carrying out the sequencing and the initial processing of the sequence data; and Chiara Barbieri, Peter de Knijff, Brigitte Pakendorf, Kay Prüfer, and Udo Stenzel for helpful discussion. This research was supported by the Max Planck Society; we also acknowledge financial support from the Biotechnology and Biological Sciences Research Council (to Francois Balloux). We dedicate this paper to the memory of Howard Cann, in recognition of his instrumental role in establishing and guiding the HGDP.
- Jobling MA: The impact of recent events on human genetic diversity. Philos Trans R Soc Lond B Biol Sci. 2012, 367: 793-799.PubMed CentralView ArticlePubMedGoogle Scholar
- Heyer E, Chaix R, Pavard S, Austerlitz F: Sex-specific demographic behaviours that shape human genomic variation. Mol Ecol. 2012, 21: 597-612.View ArticlePubMedGoogle Scholar
- Wilkins JF: Unraveling male and female histories from human genetic data. Curr Opin Genet Dev. 2006, 16: 611-617.View ArticlePubMedGoogle Scholar
- Dupanloup I, Pereira L, Bertorelle G, Calafell F, Prata MJ, Amorim A, Barbujani G: A recent shift from polygyny to monogamy in humans is suggested by the analysis of worldwide Y-chromosome diversity. J Mol Evol. 2003, 57: 85-97.View ArticlePubMedGoogle Scholar
- Wei W, Ayub Q, Xue Y, Tyler-Smith C: A comparison of Y-chromosomal lineage dating using either resequencing or Y-SNP plus Y-STR genotyping. Forensic Sci Int Genet. 2013, 7: 568-572.PubMed CentralView ArticlePubMedGoogle Scholar
- Seielstad MT, Minch E, Cavalli-Sforza LL: Genetic evidence for a higher female migration rate in humans. Nat Genet. 1998, 20: 278-280.View ArticlePubMedGoogle Scholar
- Wilder JA, Kingan SB, Mobasher Z, Pilkington MM, Hammer MF: Global patterns of human mitochondrial DNA and Y-chromosome structure are not influenced by higher migration rates of females versus males. Nat Genet. 2004, 36: 1122-1125.View ArticlePubMedGoogle Scholar
- Hammer MF, Mendez FL, Cox MP, Woerner AE, Wall JD: Sex-biased evolutionary forces shape genomic patterns of human diversity. PLoS Genet. 2008, 4: e1000202-PubMed CentralView ArticlePubMedGoogle Scholar
- Keinan A, Mullikin JC, Patterson N, Reich D: Accelerated genetic drift on chromosome X during the human dispersal out of Africa. Nat Genet. 2009, 41: 66-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Emery LS, Felsenstein J, Akey JM: Estimators of the human effective sex ratio detect sex biases on different timescales. Am J Hum Genet. 2010, 87: 848-856.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammer MF, Woerner AE, Mendez FL, Watkins JC, Cox MP, Wall JD: The ratio of human X chromosome to autosome diversity is positively correlated with genetic distance from genes. Nat Genet. 2010, 42: 830-831.View ArticlePubMedGoogle Scholar
- Arbiza L, Gottipati S, Siepel A, Keinan A: Contrasting X-linked and autosomal diversity across 14 human populations. Am J Hum Genet. 2014, 94: 827-844.PubMed CentralView ArticlePubMedGoogle Scholar
- Gottipati S, Arbiza L, Siepel A, Clark AG, Keinan A: Analyses of X-linked and autosomal genetic variation in population-scale whole genome sequencing. Nat Genet. 2011, 43: 741-743.PubMed CentralView ArticlePubMedGoogle Scholar
- Francalacci P, Morelli L, Angius A, Berutti R, Reinier F, Atzeni R, Pilu R, Busonero F, Maschio A, Zara I, Sanna D, Useli A, Urru MF, Marcelli M, Cusano R, Oppo M, Zoledziewska M, Pitzalis M, Deidda F, Porcu E, Poddie F, Kang HM, Lyons R, Tarrier B, Gresham JB, Li B, Tofanelli S, Alonso S, Dei M, Lai S, et al: Low-pass DNA sequencing of 1200 Sardinians reconstructs European Y-chromosome phylogeny. Science. 2013, 341: 565-569.View ArticlePubMedGoogle Scholar
- Poznik GD, Henn BM, Yee MC, Sliwerska E, Euskirchen GM, Lin AA, Snyder M, Quintana-Murci L, Kidd JM, Underhill PA, Bustamante CD: Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females. Science. 2013, 341: 562-565.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei W, Ayub Q, Chen Y, McCarthy S, Hou Y, Carbone I, Xue Y, Tyler-Smith C: A calibrated human Y-chromosomal phylogeny based on resequencing. Genome Res. 2013, 23: 388-395.PubMed CentralView ArticlePubMedGoogle Scholar
- Cann HM, de Toma C, Cazes L, Legrand MF, Morel V, Piouffre L, Bodmer J, Bodmer WF, Bonne-Tamir B, Cambon-Thomsen A, Chen Z, Chu J, Carcassi C, Contu L, Du R, Excoffier L, Ferrara GB, Friedlaender JS, Groot H, Gurwitz D, Jenkins T, Herrera RJ, Huang X, Kidd J, Kidd KK, Langaney A, Lin AA, Mehdi SQ, Parham P, Piazza A, et al: A human genome diversity cell line panel. Science. 2002, 296: 261-262.View ArticlePubMedGoogle Scholar
- Rosenberg NA: Standardized subsets of the HGDP-CEPH Human Genome Diversity Cell Line Panel, accounting for atypical and duplicated samples and pairs of close relatives. Ann Hum Genet. 2006, 70: 841-847.View ArticlePubMedGoogle Scholar
- Meyer M, Kircher M: Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010, 2010: pdb.prot5448-View ArticlePubMedGoogle Scholar
- Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, Kent WJ: The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004, 32: D493-D496.PubMed CentralView ArticlePubMedGoogle Scholar
- Hodges E, Xuan Z, Balija V, Kramer M, Molla MN, Smith SW, Middle CM, Rodesch MJ, Albert TJ, Hannon GJ, McCombie WR: Genome-wide in situ exon capture for selective resequencing. Nat Genet. 2007, 39: 1522-1527.View ArticlePubMedGoogle Scholar
- Hodges E, Rooks M, Xuan Z, Bhattacharjee A: Benjamin Gordon D, Brizuela L, Richard McCombie W, Hannon GJ: Hybrid selection of discrete genomic intervals on custom-designed microarrays for massively parallel sequencing. Nat Protoc. 2009, 4: 960-974.PubMed CentralView ArticlePubMedGoogle Scholar
- Maricic T, Whitten M, Paabo S: Multiplexed DNA sequence capture of mitochondrial genomes using PCR products. PLoS One. 2010, 5: e14004-PubMed CentralView ArticlePubMedGoogle Scholar
- Kircher M, Stenzel U, Kelso J: Improved base calling for the Illumina Genome Analyzer using machine learning strategies. Genome Biol. 2009, 10: R83-PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: Genome Project Data Processing S: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079.PubMed CentralView ArticlePubMedGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43: 491-498.PubMed CentralView ArticlePubMedGoogle Scholar
- Briggs AW, Good JM, Green RE, Krause J, Maricic T, Stenzel U, Lalueza-Fox C, Rudan P, Brajkovic D, Kucan Z, Gusic I, Schmitz R, Doronichev VB, Golovanova LV, de la Rasilla M, Fortea J, Rosas A, Paabo S: Targeted retrieval and analysis of five Neandertal mtDNA genomes. Science. 2009, 325: 318-321.View ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797.PubMed CentralView ArticlePubMedGoogle Scholar
- Kloss-Brandstatter A, Pacher D, Schonherr S, Weissensteiner H, Binna R, Specht G, Kronenberg F: HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups. Hum Mutat. 2011, 32: 25-32.View ArticlePubMedGoogle Scholar
- Karafet TM, Mendez FL, Meilerman MB, Underhill PA, Zegura SL, Hammer MF: New binary polymorphisms reshape and increase resolution of the human Y chromosomal haplogroup tree. Genome Res. 2008, 18: 830-838.PubMed CentralView ArticlePubMedGoogle Scholar
- Underhill PA, Passarino G, Lin AA, Shen P, Mirazon Lahr M, Foley RA, Oefner PJ, Cavalli-Sforza LL: The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations. Ann Hum Genet. 2001, 65: 43-62.View ArticlePubMedGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resources. 2010, 10: 564-567.View ArticleGoogle Scholar
- Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004, 20: 289-290.View ArticlePubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973.PubMed CentralView ArticlePubMedGoogle Scholar
- Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256.View ArticlePubMedGoogle Scholar
- Soares P, Ermini L, Thomson N, Mormina M, Rito T, Rohl A, Salas A, Oppenheimer S, Macaulay V, Richards MB: Correcting for purifying selection: an improved human mitochondrial molecular clock. Am J Hum Genet. 2009, 84: 740-759.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Q, Mittnik A, Johnson PL, Bos K, Lari M, Bollongino R, Sun C, Giemsch L, Schmitz R, Burger J, Ronchitelli AM, Martini F, Cremenosi RG, Svoboda J, Bauer P, Caramelli D, Castellano S, Reich D, Paabo S, Krause J: A revised timescale for human evolution based on ancient mitochondrial genomes. Curr Biol. 2013, 23: 553-559.View ArticlePubMedGoogle Scholar
- Rieux A, Eriksson A, Li M, Sobkowiak B, Weinert LA, Warmuth V, Ruiz-Linares A, Manica A, Balloux F: Improved calibration of the human mitochondrial clock using ancient genomes. Mol Biol Evol. 2014, in pressGoogle Scholar
- Xue YL, Wang QJ, Long Q, Ng BL, Swerdlow H, Burton J, Skuce C, Taylor R, Abdellah Z, Zhao YL, Asan , MacArthur DG, Quail MA, Carter NP, Yang H, Tyler-Smith C: Human Y chromosome base-substitution mutation rate measured by direct sequencing in a deep-rooting pedigree. Curr Biol. 2009, 19: 1453-1457.PubMed CentralView ArticlePubMedGoogle Scholar
- Mendez FL, Krahn T, Schrack B, Krahn AM, Veeramah KR, Woerner AE, Fomine FL, Bradman N, Thomas MG, Karafet TM, Hammer MF: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005, 22: 1185-1192.View ArticlePubMedGoogle Scholar
- Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319: 1100-1104.View ArticlePubMedGoogle Scholar
- Reich D, Patterson N, Kircher M, Delfin F, Nandineni MR, Pugach I, Ko AM, Ko YC, Jinam TA, Phipps ME, Saitou N, Wollstein A, Kayser M, Paabo S, Stoneking M: Denisova admixture and the first modern human dispersals into Southeast Asia and Oceania. Am J Hum Genet. 2011, 89: 516-528.PubMed CentralView ArticlePubMedGoogle Scholar
- Wollstein A, Lao O, Becker C, Brauer S, Trent RJ, Nurnberg P, Stoneking M, Kayser M: Demographic history of Oceania inferred from genome-wide data. Curr Biol. 2010, 20: 1983-1992.View ArticlePubMedGoogle Scholar
- Rasmussen M, Guo X, Wang Y, Lohmueller KE, Rasmussen S, Albrechtsen A, Skotte L, Lindgreen S, Metspalu M, Jombart T, Kivisild T, Zhai W, Eriksson A, Manica A, Orlando L, De La Vega FM, Tridico S, Metspalu E, Nielsen K, Avila-Arcos MC, Moreno-Mayar JV, Muller C, Dortch J, Gilbert MT, Lund O, Wesolowska A, Karmin M, Weinert LA, Wang B, Li J, et al: An Aboriginal Australian genome reveals separate human dispersals into Asia. Science. 2011, 334: 94-98.PubMed CentralView ArticlePubMedGoogle Scholar
- Lopez Herraez D, Bauchet M, Tang K, Theunert C, Pugach I, Li J, Nandineni MR, Gross A, Scholz M, Stoneking M: Genetic variation and recent positive selection in worldwide human populations: evidence from nearly 1 million SNPs. PLoS One. 2009, 4: e7888-PubMed CentralView ArticlePubMedGoogle Scholar
- Excoffier L, Foll M: fastsimcoal: a continuous-time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios. Bioinformatics. 2011, 27: 1332-1334.View ArticlePubMedGoogle Scholar
- Beaumont MA, Zhang W, Balding DJ: Approximate Bayesian computation in population genetics. Genetics. 2002, 162: 2025-2035.PubMed CentralPubMedGoogle Scholar
- Wilson Sayres MA, Lohmueller KE, Nielsen R: Natural selection reduced diversity on human y chromosomes. PLoS Genet. 2014, 10: e1004064-PubMed CentralView ArticlePubMedGoogle Scholar
- Scally A, Durbin R: Revising the human mutation rate: implications for understanding human evolution. Nat Rev Genet. 2012, 13: 745-753.View ArticlePubMedGoogle Scholar
- Elhaik E, Tatarinova TV, Klyosov AA, Graur D: The ‘extremely ancient’ chromosome that isn’t: a forensic bioinformatic investigation of Albert Perry’s X-degenerate portion of the Y chromosome. Eur J Hum Genet. 2014, 22: 1111-1116.PubMed CentralView ArticlePubMedGoogle Scholar
- van Oven M, Kayser M: Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009, 30: E386-E394.View ArticlePubMedGoogle Scholar
- Behar DM, Villems R, Soodyall H, Blue-Smith J, Pereira L, Metspalu E, Scozzari R, Makkan H, Tzur S, Comas D, Bertranpetit J, Quintana-Murci L, Tyler-Smith C, Wells RS: Rosset S; Genographic Consortium: The dawn of human matrilineal diversity. Am J Hum Genet. 2008, 82: 1130-1140.PubMed CentralView ArticlePubMedGoogle Scholar
- Zegura SL, Karafet TM, Zhivotovsky LA, Hammer MF: High-resolution SNPs and microsatellite haplotypes point to a single, recent entry of Native American Y chromosomes into the Americas. Mol Biol Evol. 2004, 21: 164-175.View ArticlePubMedGoogle Scholar
- O’Rourke DH, Raff JA: The human genetic history of the Americas: the final frontier. Curr Biol. 2010, 20: R202-R207.View ArticlePubMedGoogle Scholar
- Ray N, Wegmann D, Fagundes NJ, Wang S, Ruiz-Linares A, Excoffier L: A statistical evaluation of models for the initial settlement of the american continent emphasizes the importance of gene flow with Asia. Mol Biol Evol. 2010, 27: 337-345.View ArticlePubMedGoogle Scholar
- Shi W, Ayub Q, Vermeulen M, Shao RG, Zuniga S, van der Gaag K, de Knijff P, Kayser M, Xue Y, Tyler-Smith C: A worldwide survey of human male demographic history based on Y-SNP and Y-STR data from the HGDP-CEPH populations. Mol Biol Evol. 2010, 27: 385-393.PubMed CentralView ArticlePubMedGoogle Scholar
- de Filippo C, Barbieri C, Whitten M, Mpoloka SW, Gunnarsdottir ED, Bostoen K, Nyambe T, Beyer K, Schreiber H, de Knijff P, Luiselli D, Stoneking M, Pakendorf B: Y-chromosomal variation in sub-Saharan Africa: insights into the history of Niger-Congo groups. Mol Biol Evol. 2011, 28: 1255-1269.PubMed CentralView ArticlePubMedGoogle Scholar
- Wood ET, Stover DA, Ehret C, Destro-Bisol G, Spedini G, McLeod H, Louie L, Bamshad M, Strassmann BI, Soodyall H, Hammer MF: Contrasting patterns of Y chromosome and mtDNA variation in Africa: evidence for sex-biased demographic processes. Eur J Hum Genet. 2005, 13: 867-876.View ArticlePubMedGoogle Scholar
- Sikora MJ, Colonna V, Xue Y, Tyler-Smith C: Modeling the contrasting Neolithic male lineage expansions in Europe and Africa. Investig Genet. 2013, 4: 25-PubMed CentralView ArticlePubMedGoogle Scholar
- Botigue LR, Henn BM, Gravel S, Maples BK, Gignoux CR, Corona E, Atzmon G, Burns E, Ostrer H, Flores C, Bertranpetit J, Comas D, Bustamante CD: Gene flow from North Africa contributes to differential human genetic diversity in southern Europe. Proc Natl Acad Sci U S A. 2013, 110: 11791-11796.PubMed CentralView ArticlePubMedGoogle Scholar
- Kayser M, Brauer S, Cordaux R, Casto A, Lao O, Zhivotovsky LA, Moyse-Faurie C, Rutledge RB, Schiefenhoevel W, Gil D, Lin AA, Underhill PA, Oefner PJ, Trent RJ, Stoneking M: Melanesian and Asian origins of Polynesians: mtDNA and Y chromosome gradients across the Pacific. Mol Biol Evol. 2006, 23: 2234-2244.View ArticlePubMedGoogle Scholar
- Trejaut JA, Kivisild T, Loo JH, Lee CL, He CL, Hsu CJ, Lee ZY, Lin M: Traces of archaic mitochondrial lineages persist in Austronesian-speaking Formosan populations. PLoS Biol. 2005, 3: e247-PubMed CentralView ArticlePubMedGoogle Scholar
- Duggan AT, Evans B, Friedlaender FR, Friedlaender JS, Koki G, Merriwether DA, Kayser M, Stoneking M: Maternal history of Oceania from complete mtDNA genomes: contrasting ancient diversity with recent homogenization due to the Austronesian expansion. Am J Hum Genet. 2014, 94: 721-733.PubMed CentralView ArticlePubMedGoogle Scholar
- Delfin F, Myles S, Choi Y, Hughes D, Illek R, van Oven M, Pakendorf B, Kayser M, Stoneking M: Bridging near and remote Oceania: mtDNA and NRY variation in the Solomon Islands. Mol Biol Evol. 2012, 29: 545-564.View ArticlePubMedGoogle Scholar
- Scheinfeldt L, Friedlaender F, Friedlaender J, Latham K, Koki G, Karafet T, Hammer M, Lorenz J: Unexpected NRY chromosome variation in Northern Island Melanesia. Mol Biol Evol. 2006, 23: 1628-1641.View ArticlePubMedGoogle Scholar
- Friedlaender JS, Friedlaender FR, Hodgson JA, Stoltz M, Koki G, Horvat G, Zhadanov S, Schurr TG, Merriwether DA: Melanesian mtDNA complexity. PLoS One. 2007, 2: e248-PubMed CentralView ArticlePubMedGoogle Scholar
- Kayser M, Choi Y, van Oven M, Mona S, Brauer S, Trent RJ, Suarkia D, Schiefenhovel W, Stoneking M: The impact of the Austronesian expansion: evidence from mtDNA and Y chromosome diversity in the Admiralty Islands of Melanesia. Mol Biol Evol. 2008, 25: 1362-1374.View ArticlePubMedGoogle Scholar
- Schurr TG: The peopling of the New World: Perspectives from molecular anthropology. Ann Rev Anthropol. 2004, 33: 551-583.View ArticleGoogle Scholar
- Gunnarsdottir ED, Li M, Bauchet M, Finstermeier K, Stoneking M: High-throughput sequencing of complete human mtDNA genomes from the Philippines. Genome Res. 2011, 21: 1-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Neuenschwander S, Largiader CR, Ray N, Currat M, Vonlanthen P, Excoffier L: Colonization history of the Swiss Rhine basin by the bullhead (Cottus gobio): inference under a Bayesian spatially explicit framework. Mol Ecol. 2008, 17: 757-772.View ArticlePubMedGoogle Scholar
- Gronau I, Hubisz MJ, Gulko B, Danko CG, Siepel A: Bayesian inference of ancient human demography from individual genome sequences. Nat Genet. 2011, 43: 1031-1034.PubMed CentralView ArticlePubMedGoogle Scholar
- Laval G, Patin E, Barreiro LB, Quintana-Murci L: Formulating a historical and demographic model of recent human evolution based on resequencing data from noncoding regions. PLoS One. 2010, 5: e10284-PubMed CentralView ArticlePubMedGoogle Scholar
- Oota H, Settheetham-Ishida W, Tiwawech D, Ishida T, Stoneking M: Human mtDNA and Y-chromosome variation is correlated with matrilocal versus patrilocal residence. Nat Genet. 2001, 29: 20-21.View ArticlePubMedGoogle Scholar
- Labuda D, Lefebvre JF, Nadeau P, Roy-Gagnon MH: Female-to-male breeding ratio in modern humans-an analysis based on historical recombinations. Am J Hum Genet. 2010, 86: 353-363.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilder JA, Mobasher Z, Hammer MF: Genetic evidence for unequal effective population sizes of human females and males. Mol Biol Evol. 2004, 21: 2047-2057.View ArticlePubMedGoogle Scholar
- Keinan A, Mullikin JC, Patterson N, Reich D: Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in Europeans. Nat Genet. 2007, 39: 1251-1255.PubMed CentralView ArticlePubMedGoogle Scholar
- Vernot B, Akey JM: Resurrecting surviving Neandertal lineages from modern human genomes. Science. 2014, 343: 1017-1021.View ArticlePubMedGoogle Scholar
- Garrigan D, Kingan SB, Pilkington MM, Wilder JA, Cox MP, Soodyall H, Strassmann B, Destro-Bisol G, de Knijff P, Novelletto A, Friedlaender J, Hammer MF: Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data. Genetics. 2007, 177: 2195-2207.PubMed CentralView ArticlePubMedGoogle Scholar
- Hellenthal G, Busby GB, Band G, Wilson JF, Capelli C, Falush D, Myers S: A genetic atlas of human admixture history. Science. 2014, 343: 747-751.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.