Ancestral alleles and population origins: inferences depend on mutation rate.

Previous studies have found that at most human loci, ancestral alleles are "African," in the sense that they reach their highest frequency there. Conventional wisdom holds that this reflects a recent African origin of modern humans. This paper challenges that view by showing that the empirical pattern (of elevated allele frequencies within Africa) is not as pervasive as has been thought. We confirm this African bias in a set of mainly protein-coding loci, but find a smaller bias in Alu insertion polymorphisms, and an even smaller bias in noncoding loci. Thus, the strong bias that was originally observed must reflect some factor that varies among data sets--something other than population history. This factor may be the per-locus mutation rate: the African bias is most pronounced in loci where this rate is high. The distribution of ancestral alleles among populations has been studied using 2 methods. One of these involves comparing the fractions of loci that reach maximal frequency in each population. The other compares the average frequencies of ancestral alleles. The first of these methods reflects history in a manner that depends on the mutation rate. When that rate is high, ancestral alleles at most loci reach their highest frequency in the ancestral population. When that rate is low, the reverse is true. The other method--comparing averages--is unresponsive. Average ancestral allele frequencies are affected neither by mutation rate nor by the history of population size and migration. In the absence of selection and ascertainment bias, they should be the same everywhere. This is true of one data set, but not of 2 others. This also suggests the action of some factor, such as selection or ascertainment bias, that varies among data sets.     

Contrasting signatures of population growth for mitochondrial DNA and Y chromosomes among human populations in Africa

A history of Pleistocene population expansion has been inferred from the frequency spectrum of polymorphism in the mitochondrial DNA (mtDNA) of many human populations. Similar patterns are not typically observed for autosomal and X-linked loci. One explanation for this discrepancy is a recent population bottleneck, with different rates of recovery for haploid and autosomal loci as a result of their different effective population sizes. This hypothesis predicts that mitochondrial and Y chromosomal DNA will show a similar skew in the frequency spectrum in populations that have experienced a recent increase in effective population size. We test this hypothesis by resequencing 6.6 kb of noncoding Y chromosomal DNA and 780 basepairs of the mtDNA cytochrome c oxidase subunit III (COIII) gene in 172 males from 5 African populations. Four tests of population expansion are employed for each locus in each population: Fu's Fs statistic, the R(2) statistic, coalescent simulations, and the mismatch distribution. Consistent with previous results, patterns of mtDNA polymorphism better fit a model of constant population size for food-gathering populations and a model of population expansion for food-producing populations. In contrast, none of the tests reveal evidence of Y chromosome growth for either food-gatherers or food-producers. The distinct mtDNA and Y chromosome polymorphism patterns most likely reflect sex-biased demographic processes in the recent history of African populations. We hypothesize that males experienced smaller effective population sizes and/or lower rates of migration during the Bantu expansion, which occurred over the last 5,000 years.     

The distribution of the ancestral haplotype in finite stepping-stone models with population expansion

Recent extensive analyses of human DNA polymorphism reveal that the ancestral haplotype at various genetic loci occurs almost exclusively in African samples. We develop a coalescence-based simulation method in stepping-stone models with population expansion and examine the probability (P(A)) that the ancestral haplotype is found in African samples and the probability (Q(A)) that the most recent common ancestor of sampled genes occurs in Africa. These probabilities and other summary statistics are used to infer the human demographic history. It is shown that the high observed P(A) value cannot be explained simply by sampling bias. Rather, it suggests that the African population has been more strongly subdivided and isolated from each other than the non-African population and that there must have been some African populations which were not directly involved in the Out-of-Africa expansion in the late Pleistocene.     

Genetic similarities within and between human populations

The proportion of human genetic variation due to differences between populations is modest, and individuals from different populations can be genetically more similar than individuals from the same population. Yet sufficient genetic data can permit accurate classification of individuals into populations. Both findings can be obtained from the same data set, using the same number of polymorphic loci. This article explains why. Our analysis focuses on the frequency, omega, with which a pair of random individuals from two different populations is genetically more similar than a pair of individuals randomly selected from any single population. We compare omega to the error rates of several classification methods, using data sets that vary in number of loci, average allele frequency, populations sampled, and polymorphism ascertainment strategy. We demonstrate that classification methods achieve higher discriminatory power than omega because of their use of aggregate properties of populations. The number of loci analyzed is the most critical variable: with 100 polymorphisms, accurate classification is possible, but omega remains sizable, even when using populations as distinct as sub-Saharan Africans and Europeans. Phenotypes controlled by a dozen or fewer loci can therefore be expected to show substantial overlap between human populations. This provides empirical justification for caution when using population labels in biomedical settings, with broad implications for personalized medicine, pharmacogenetics, and the meaning of race.     

Is the doubly deleted alpha-thalassemia gene a "fugitive" allele?

The alpha-thalassemias (particularly in Asia) can be considered as a three-allele system, with one "normal" allele (N) consisting of a pair of closely linked alpha-chain loci, a second "single" allele (S) in which one of the loci has been lost by deletion, and a third "double" allele (D) in which both have been lost. Representatives of all the sets of fitnesses leading to the maintenance of this condition by selection for malaria resistance have been explored, and after the discarding of unlikely sets of fitnesses, it is found that there will be three outcomes: (1) the S chromosome is fixed, (2) the S and N chromosomes form a stable polymorphism, and (3) the N and D chromosomes form a stable polymorphism, but this can be lost and the population forced to fixation by the introduction of sufficiently large number of S chromosomes. Some Melanesian populations appear to have reached outcome (1), while frequencies in African, Mediterranean, and Middle Eastern populations are not incompatible with outcome (2). Southeast Asian populations, however, which carry S and D chromosomes in high frequency, may be in a state of flux. The D chromosome may form a polymorphism with N locally, but it can be driven from the local population by the introduction of large numbers of S chromosomes. The D chromosome would thus be somewhat analogous to a fugitive species, which can only exist in certain transient environments and is displaced as the environment changes. The possibility that N, S, and D are coexisting as a stable polymorphism can almost certainly be ruled out by a consideration of fitness sets required.     

Approximate Bayesian analysis of Drosophila melanogaster polymorphism data reveals a recent colonization of Southeast Asia

Southeast Asian populations of the fruit fly Drosophila melanogaster differ from ancestral African and derived European populations by several morphological characteristics. It has been argued that this morphological differentiation could be the result of an early colonization of Southeast Asia that predated the migration of D. melanogaster to Europe after the last glacial period (around 10,000 years ago). To investigate the colonization process of Southeast Asia, we collected nucleotide polymorphism data for more than 200 X-linked fragments and 50 autosomal loci from a population of Malaysia. We analyzed this new single nucleotide polymorphism data set jointly with already existing data from an African and a European population by employing an Approximate Bayesian Computation approach. By contrasting different demographic models of these three populations, we do not find any evidence for an early divergence between the African and the Asian populations. Rather, we show that Asian and European populations of D. melanogaster share a non-African most recent common ancestor that existed about 2,500 years ago.     

SINEs of speciation: tracking lineages with retroposons

The value of short interspersed elements (SINEs) for diagnosing common ancestry is being expanded to examine the differential sorting of lineages through the course of speciation events. Because most SINEs are neutral markers of identical descent, are not precisely excised from the genome and have a known ancestral condition, they are advantageous for reconciling gene trees and species trees with minimal phylogenetic error. A population perspective on SINE evolution combined with coalescence theory provides a context for investigating the phenomenon of ancestral polymorphism and its role in producing incongruent SINE insertion patterns among multiple loci. Studies of human Alu repeats demonstrate the value of young polymorphic SINEs for assessing human genomic diversity and tracking ancient demographics of human populations, whereas incongruent insertion patterns revealed by older fixed SINE loci, such as those in African cichlid fishes, contain information that might help identify ancient radiations that are otherwise obscured by accumulated mutations in sequence data. Here, we review the utility of retroposons for inferring common ancestry, discuss limits to the method, and clarify confusion by providing examples from the literature that illustrate how discordant multi-locus insertion patterns of retroelements can indicate lineage-sorting events that should not be misinterpreted as phylogenetic noise.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.     

A genetic affinity analysis of human populations

Genetic affinity of human populations based on allele frequency data was studied from two viewpoints. (1) The effect of the number of polymorphic loci on the reconstruction of a phylogenetic tree of human populations was empirically investigated. Genetic affinity trees were constructed based on data for 1–12 polymorphic loci, by using the neighbor-joining method. Geographical clustering of populations gradually appeared when the number of loci was increased. A new classification and terminology of higher order human population clusters is proposed based on these and other studies. (2) A new method of estimating the absolute divergence time of two populations is proposed, which is based on a diffusion equation that describes random genetic drift.     

Non-African populations of Drosophila melanogaster have a unique origin

Drosophila melanogaster is widely used as a model in DNA variation studies. Patterns of polymorphism have, however, been affected by the history of this species, which is thought to have recently spread out of Africa to the rest of the world. We analyzed DNA sequence variation in 11 populations, including four continental African and seven non-African samples (including Madagascar), at four independent X-linked loci. Variation patterns at all four loci followed neutral expectations in all African populations, but departed from it in all non-African ones due to a marked haplotype dimorphism at three out of four loci. We also found that all non-African populations show the same major haplotypes, though in various frequencies. A parsimonious explanation for these observations is that all non-African populations are derived from a single ancestral population having undergone a substantial reduction of polymorphism, probably through a bottleneck. Less likely alternatives involve either selection at all four loci simultaneously (including balancing selection at three of them), or admixture between two divergent populations. Small but significant structure was observed among African populations, and there were indications of differentiation across Eurasia for non-African ones. Since population history may result in non-equilibrium variation patterns, our study confirms that the search for footprints of selection in the D. melanogaster genome must include a sufficient understanding of its history.     

mtDNA variation in the South African Kung and Khwe-and their genetic relationships to other African populations

The mtDNA variation of 74 Khoisan-speaking individuals (Kung and Khwe) from Schmidtsdrift, in the Northern Cape Province of South Africa, was examined by high-resolution RFLP analysis and control region (CR) sequencing. The resulting data were combined with published RFLP haplotype and CR sequence data from sub-Saharan African populations and then were subjected to phylogenetic analysis to deduce the evolutionary relationships among them. More than 77% of the Kung and Khwe mtDNA samples were found to belong to the major mtDNA lineage, macrohaplogroup L* (defined by a HpaI site at nucleotide position 3592), which is prevalent in sub-Saharan African populations. Additional sets of RFLPs subdivided macrohaplogroup L* into two extended haplogroups-L1 and L2-both of which appeared in the Kung and Khwe. Besides revealing the significant substructure of macrohaplogroup L* in African populations, these data showed that the Biaka Pygmies have one of the most ancient RFLP sublineages observed in African mtDNA and, thus, that they could represent one of the oldest human populations. In addition, the Kung exhibited a set of related haplotypes that were positioned closest to the root of the human mtDNA phylogeny, suggesting that they, too, represent one of the most ancient African populations. Comparison of Kung and Khwe CR sequences with those from other African populations confirmed the genetic association of the Kung with other Khoisan-speaking peoples, whereas the Khwe were more closely linked to non-Khoisan-speaking (Bantu) populations. Finally, the overall sequence divergence of 214 African RFLP haplotypes defined in both this and an earlier study was 0.364%, giving an estimated age, for all African mtDNAs, of 125,500-165,500 years before the present, a date that is concordant with all previous estimates derived from mtDNA and other genetic data, for the time of origin of modern humans in Africa.     

The Episode of Genetic Drift Defining the Migration of Humans out of Africa Is Derived from a Large East African Population Size

Human genetic variation particularly in Africa is still poorly understood. This is despite a consensus on the large African effective population size compared to populations from other continents. Based on sequencing of the mitochondrial Cytochrome C Oxidase subunit II (MT-CO2), and genome wide microsatellite data we observe evidence suggesting the effective size (Ne) of humans to be larger than the current estimates, with a foci of increased genetic diversity in east Africa, and a population size of east Africans being at least 2-6 fold larger than other populations. Both phylogenetic and network analysis indicate that east Africans possess more ancestral lineages in comparison to various continental populations placing them at the root of the human evolutionary tree. Our results also affirm east Africa as the likely spot from which migration towards Asia has taken place. The study reflects the spectacular level of sequence variation within east Africans in comparison to the global sample, and appeals for further studies that may contribute towards filling the existing gaps in the database. The implication of these data to current genomic research, as well as the need to carry out defined studies of human genetic variation that includes more African populations; particularly east Africans is paramount.     

Historicity and the Population Genetics of Drosophila Melanogaster and D. Simulans

We summarize data showing that there is population structure in African populations of Drosophila from the melanogaster-simulans complex. In D. melanogaster, population structuring is found at individual loci, but is obscured by population structuring for large inversions that simultaneously affect several loci. In D. simulans, molecular polymorphism at the X-linked vermilion locus suggests that different groups of populations have been geographically isolated for some time. Invading populations are probably derived from different areas in Africa. European populations originate from an east African population that was probably not at a demographic equilibrium. The origin of the Antilles population is apparently different and is as yet unknown. In south-western France, populations from these two species undergo different population structuring at the scale of a few kilometres: D. melanogaster makes up a large panmictic population, whereas D. simulans forms a metapopulation that is divided into smaller demes.     

Accurate gene diversity estimates from amplified fragment length polymorphism (AFLP) markers

Three procedures for the estimation of null allele frequencies and gene diversity from dominant multilocus data were empirically tested in natural populations of the outcrossing angiosperm Persoonia mollis (Proteaceae). The three procedures were the square root transform of the null homozygote frequency, the Lynch & Milligan procedure, and the Bayesian method. Genotypes for each of 116 polymorphic loci generated by amplified fragment length polymorphism (AFLP) were inferred from segregation patterns in progeny arrays. Therefore, for the plus phenotype (band present), heterozygotes were distinguished from homozygotes. In contrast to previous studies, all three procedures produced very similar mean estimates of heterozygosity, which were in turn accurate estimators of the direct value (HO = 0.28). A second population of P. mollis displayed markedly lower levels of heterozygosity (HO = 0.20) but approximately twice as many polymorphic loci (284). These AFLP results show that biases in estimates of average null allele frequency and heterozygosity are largely eliminated in highly polymorphic dominant marker data sets displaying a J-shaped beta distribution with a high percentage of loci containing more than three null homozygotes and relatively few loci with no null homozygotes. This distribution may be typical of outcrossing angiosperms.     

Evolutionary relationships of human populations on a global scale

Using gene frequency data for 29 polymorphic loci (121 alleles), we conducted a phylogenetic analysis of 26 representative populations from around the world by using the neighbor-joining (NJ) method. We also conducted a separate analysis of 15 populations by using data for 33 polymorphic loci. These analyses have shown that the first major split of the phylogenetic tree separates Africans from non-Africans and that this split occurs with a 100% bootstrap probability. The second split separates Caucasian populations from all other non-African populations, and this split is also supported by bootstrap tests. The third major split occurs between Native American populations and the Greater Asians that include East Asians (mongoloids), Pacific Islanders, and Australopapuans (native Australians and Papua New Guineans), but Australopapuans are genetically quite different from the rest of the Greater Asians. The second and third levels of population splitting are quite different from those of the phylogenetic tree obtained by Cavalli- Sforza et al. (1988), where Caucasians, Northeast Asians, and Ameridians from the Northeurasian supercluster and the rest of non- Africans form the Southeast Asian supercluster. One of the major factors that caused the difference between the two trees is that Cavalli-Sforza et al. used unweighted pair-group method with arithmetic mean (UPGMA) in phylogenetic inference, whereas we used the NJ method in which evolutionary rate is allowed to vary among different populations. Bootstrap tests have shown that the UPGMA tree receives poor statistical support whereas the NJ tree is well supported. Implications that the phylogenetic tree obtained has on the current controversy over the out-of-Africa and the multiregional theories of human origins are discussed.      

Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data

We estimate parameters of a general isolation-with-migration model using resequence data from mitochondrial DNA (mtDNA), the Y chromosome, and two loci on the X chromosome in samples of 25-50 individuals from each of 10 human populations. Application of a coalescent-based Markov chain Monte Carlo technique allows simultaneous inference of divergence times, rates of gene flow, as well as changes in effective population size. Results from comparisons between sub-Saharan African and Eurasian populations estimate that 1500 individuals founded the ancestral Eurasian population approximately 40 thousand years ago (KYA). Furthermore, these small Eurasian founding populations appear to have grown much more dramatically than either African or Oceanian populations. Analyses of sub-Saharan African populations provide little evidence for a history of population bottlenecks and suggest that they began diverging from one another upward of 50 KYA. We surmise that ancestral African populations had already been geographically structured prior to the founding of ancestral Eurasian populations. African populations are shown to experience low levels of mitochondrial DNA gene flow, but high levels of Y chromosome gene flow. In particular, Y chromosome gene flow appears to be asymmetric, i.e., from the Bantu-speaking population into other African populations. Conversely, mitochondrial gene flow is more extensive between non-African populations, but appears to be absent between European and Asian populations.     

When did the human population size start increasing?

We analyze the frequency spectra of all available human nuclear sequence data sets by using a model of constant population size followed by exponential growth. Parameters of growth (more extreme than or) comparable to what has been suggested from mtDNA data can be rejected for 6 out of the 10 largest data sets. When the data are separated into African and non-African samples, a constant size no-growth model can be rejected for 4 out of 8 non-African samples. Long-term growth (i.e., starting 50-100 kya) can be rejected for 2 out of 8 African samples and 5 out of 8 non-African ones. Under more complex demographic models, including a bottleneck or population subdivision, more of the data are compatible with long-term growth. One problem with the data used here is that a subset of loci may reflect the action of natural selection as well as of demography. It remains possible that the correct demographic model is one of constant population size followed by long-term growth but that at several loci the demographic signature has been obscured by balancing or diversifying selection. However, it is not clear that the data at these loci are consistent with a simple model of balancing selection; more complicated selective alternatives cannot be tested unless they are made explicit. An alternative explanation is that population size growth is more recent (e.g., upper Paleolithic) and that some of the loci have experienced recent directional selection. Given the available data, the latter hypothesis seems more likely.     

Distinctly different sex ratios in African and European populations of Drosophila melanogaster inferred from chromosomewide single nucleotide polymorphism data

It has been hypothesized that the ratio of X-linked to autosomal sequence diversity is influenced by unequal sex ratios in Drosophila melanogaster populations. We conducted a genome scan of single nucleotide polymorphism (SNP) of 378 autosomal loci in a derived European population and of a subset of 53 loci in an ancestral African population. On the basis of these data and our already available X-linked data, we used a coalescent-based maximum-likelihood method to estimate sex ratios and demographic histories simultaneously for both populations. We confirm our previous findings that the African population experienced a population size expansion while the European population suffered a population size bottleneck. Our analysis also indicates that the female population size in Africa is larger than or equal to the male population size. In contrast, the European population shows a huge excess of males. This unequal sex ratio and the bottleneck alone, however, cannot account for the overly strong decrease of X-linked diversity in the European population (compared to the reduction on the autosome). The patterns of the frequency spectrum and the levels of linkage disequilibrium observed in Europe suggest that, in addition, positive selection must have acted in the derived population.     

Empirical tests of the reliability of phylogenetic trees constructed with microsatellite DNA

Microsatellite DNA loci or short tandem repeats (STRs) are abundant in eukaryotic genomes and are often used for constructing phylogenetic trees of closely related populations or species. These phylogenetic trees are usually constructed by using some genetic distance measure based on allele frequency data, and there are many distance measures that have been proposed for this purpose. In the past the efficiencies of these distance measures in constructing phylogenetic trees have been studied mathematically or by computer simulations. Recently, however, allele frequencies of 783 STR loci have been compiled from various human populations. We have therefore used these empirical data to investigate the relative efficiencies of different distance measures in constructing phylogenetic trees. The results showed that (1) the probability of obtaining the correct branching pattern of a tree (PC) is generally highest for DA distance; (2) FST*, standard genetic distance (DS), and FST/(1-FST) give similar PC-values, FST* being slightly better than the other two; and (3) (deltamu)2 shows PC-values much lower than the other distance measures. To have reasonably high PC-values for trees similar to ours, at least 30 loci with a minimum of 15 individuals are required when DA distance is used.     

Effects of ascertainment bias on recovering human demographic history.

In recent years multilocus data sets have been used to study the demographic history of human populations. In this paper (1) analyses previously done on 60 short tandem repeat (STR) loci are repeated on 30 restriction site polymorphism (RSP) markers; (2) relative population weights are estimated from the RSP data set and compared to previously published estimates from STR and craniometric data sets; and (3) computer simulations are performed to show the effects of ascertainment bias on relative population weight estimates. Not surprisingly, given that the RSP markers were originally identified in a small panel of Caucasians, estimates of relative population weights are biased and the European population weight is artificially inflated. However, the effects of ascertainment bias are not apparent in a principal components plot or estimates of FST. Ascertainment bias can have a large effect in other genetic systems with inherently low heterozygosity such as Alus or single nucleotide polymorphisms (SNPs), and care must be taken to have prior knowledge of how polymorphic markers in a given data set were originally identified. Otherwise, results can be skewed and interpretations faulty.