Microsatellite diversity and the demographic history of modern humans

Estimating divergence time with the use of microsatellite genetic distances: impacts of population growth and gene flow

Genetic distances play an important role in estimating divergence time of bifurcated populations. However, they can be greatly affected by demographic processes, such as migration and population dynamics, which complicate their interpretation. For example, the widely used distance for microsatellite loci, (deltamu)2, assumes constant population size, no gene flow, and mutation-drift equilibrium. It is shown here that (deltamu)2 strongly underestimates divergence time if populations are growing and/or connected by gene flow. In recent publications, the average estimate of divergence time between African and non-African populations obtained by using (deltamu)2 is about 34,000 years, although archaeological data show a much earlier presence of modern humans out of Africa. I introduce a different estimator of population separation time based on microsatellite statistics, T(D), that does not assume mutation-drift equilibrium, is independent of population dynamics in the absence of gene flow, and is robust to weak migration flow for growing populations. However, it requires a knowledge of the variance in the number of repeats at the beginning of population separation, V(0). One way to overcome this problem is to find minimal and maximal bounds for the variance and thus obtain the earliest and latest bounds for divergence time (this is not a confidence interval, and it simply reflects an uncertainty about the value of V(0) in an ancestral population). Another way to avoid the uncertainty is to choose from among present populations a reference whose variation is presumably close to what it might have been in an ancestral population. A different approach for using T(D) is to estimate the time difference between adjacent nodes on a phylogenetic population tree. Using data on variation at autosomal short tandem repeat loci with di-, tri-, and tetranucleotide repeats in worldwide populations, T(D) gives an estimate of 57,000 years for the separation of the out-of-Africa branch of modern humans from Africans based on the value of V(0) in the Southern American Indian populations; the earliest bound for this event has been estimated to be about 135,000 years. The data also suggest that the Asian and European populations diverged from each other about 20,000 years, after the occurrence of the out-of-Africa branch.

Genetic predisposition and somatic diversification in tumor development and progression

Studies on human cancer predisposition syndromes have contributed significantly to our understanding on tumor initiation and progression. Work performed on hereditary colon cancer has been particularly fruitful. Much of the molecular background of the various intestinal polyposis syndromes, such as familial adenomatous polyposis (FAP), juvenile polyposis, and Peutz-Jeghers syndrome, has been revealed, pinpointing several key cancer-associated genes. Studies on hereditary nonpolyposis colorectal cancer (HNPCC) have revealed a novel mechanism of tumorigenesis; genomic instability caused by defective DNA mismatch repair (MMR). Understanding the molecular background of these diseases helps us to understand tumor initiation in the affected individuals. Relatively little is known about the details of tumor progression in hereditary and sporadic neoplasia. Certain additional gene mutations can be associated with advancing stages of the disease, but the pace and tempo of the process have remained obscure. A high mutation rate in MMR-deficient tumors has provided a new approach in the analysis of human tumor dynamics. Microsatellite (MS) sequences are frequently mutated in MMR deficient tumors. The high mutation rate allows the use of microsatellite mutations as a tool for analyzing the past patterns of tumor progression. This approach is similar to the use of MS mutations in studying human evolution and migrations. Such tumor studies have revealed progression pathways that differ from the classic adenoma-cancer sequence. The reasons why and how molecular clocks may reveal something new about a well-studied problem are discussed.

Patterns of ancestral human diversity: an analysis of Alu-insertion and restriction-site polymorphisms

We have analyzed 35 widely distributed, polymorphic Alu loci in 715 individuals from 31 world populations. The average frequency of Alu insertions (the derived state) is lowest in Africa (.42) but is higher and similar in India (.55), Europe (.56), and Asia (.57). A comparison with 30 restriction-site polymorphisms (RSPs) for which the ancestral state has been determined shows that the frequency of derived RSP alleles is also lower in Africa (.35) than it is in Asia (.45) and in Europe (.46). Neighbor-joining networks based on Alu insertions or RSPs are rooted in Africa and show African populations as separate from other populations, with high statistical support. Correlations between genetic distances based on Alu and nuclear RSPs, short tandem-repeat polymorphisms, and mtDNA, in the same individuals, are high and significant. For the 35 loci, Alu gene diversity and the diversity attributable to population subdivision is highest in Africa but is lower and similar in Europe and Asia. The distribution of ancestral alleles is consistent with an origin of early modern human populations in sub-Saharan Africa, the isolation and preservation of ancestral alleles within Africa, and an expansion out of Africa into Eurasia. This expansion is characterized by increasing frequencies of Alu inserts and by derived RSP alleles with reduced genetic diversity in non-African populations.

Worldwide genetic analysis of the CFTR region

Mutations at the cystic fibrosis transmembrane conductance regulator gene (CFTR) cause cystic fibrosis, the most prevalent severe genetic disorder in individuals of European descent. We have analyzed normal allele and haplotype variation at four short tandem repeat polymorphisms (STRPs) and two single-nucleotide polymorphisms (SNPs) in CFTR in 18 worldwide population samples, comprising a total of 1,944 chromosomes. The rooted phylogeny of the SNP haplotypes was established by typing ape samples. STRP variation within SNP haplotype backgrounds was highest in most ancestral haplotypes-although, when STRP allele sizes were taken into account, differences among haplotypes became smaller. Haplotype background determines STRP diversity to a greater extent than populations do, which indicates that haplotype backgrounds are older than populations. Heterogeneity among STRPs can be understood as the outcome of differences in mutation rate and pattern. STRP sites had higher heterozygosities in Africans, although, when whole haplotypes were considered, no significant differences remained. Linkage disequilibrium (LD) shows a complex pattern not easily related to physical distance. The analysis of the fraction of possible different haplotypes not found may circumvent some of the methodological difficulties of LD measure. LD analysis showed a positive correlation with locus polymorphism, which could partly explain the unusual pattern of similar LD between Africans and non-Africans. The low values found in non-Africans may imply that the size of the modern human population that emerged "Out of Africa" may be larger than what previous LD studies suggested.

Short tandem-repeat polymorphism/alu haplotype variation at the PLAT locus: implications for modern human origins

Two dinucleotide short tandem-repeat polymorphisms (STRPs) and a polymorphic Alu element spanning a 22-kb region of the PLAT locus on chromosome 8p12-q11.2 were typed in 1,287-1,420 individuals originating from 30 geographically diverse human populations, as well as in 29 great apes. These data were analyzed as haplotypes consisting of each of the dinucleotide repeats and the flanking Alu insertion/deletion polymorphism. The global pattern of STRP/Alu haplotype variation and linkage disequilibrium (LD) is informative for the reconstruction of human evolutionary history. Sub-Saharan African populations have high levels of haplotype diversity within and between populations, relative to non-Africans, and have highly divergent patterns of LD. Non-African populations have both a subset of the haplotype diversity present in Africa and a distinct pattern of LD. The pattern of haplotype variation and LD observed at the PLAT locus suggests a recent common ancestry of non-African populations, from a small population originating in eastern Africa. These data indicate that, throughout much of modern human history, sub-Saharan Africa has maintained both a large effective population size and a high level of population substructure. Additionally, Papua New Guinean and Micronesian populations have rare haplotypes observed otherwise only in African populations, suggesting ancient gene flow from Africa into Papua New Guinea, as well as gene flow between Melanesian and Micronesian populations.

Archaic lineages in the history of modern humans

An important question in the ongoing debate on the origin of Homo sapiens is whether modern human populations issued from a single lineage or whether several, independently evolving lineages contributed to their genetic makeup. We analyzed haplotypes composed of 35 polymorphisms from a segment of the dystrophin gene. We find that the bulk of a worldwide sample of 868 chromosomes represents haplotypes shared by different continental groups. The remaining chromosomes carry haplotypes specific for the continents or for local populations. The haplotypes specific for non-Africans can be derived from the most frequent ones through simple recombination or a mutation. In contrast, chromosomes specific for sub-Saharan Africans represent a distinct group, as shown by principal component analysis, maximum likelihood tree, structural comparison, and summary statistics. We propose that African chromosomes descend from at least two lineages that have been evolving separately for a period of time. One of them underwent range expansion colonizing different continents, including Africa, where it mixed with another, local lineage represented today by a large fraction of African-specific haplotypes. Genetic admixture involving archaic lineages appears therefore to have occurred within Africa rather than outside this continent, explaining greater diversity of sub-Saharan populations observed in a variety of genetic systems.

Association mapping in structured populations

The use, in association studies, of the forthcoming dense genomewide collection of single-nucleotide polymorphisms (SNPs) has been heralded as a potential breakthrough in the study of the genetic basis of common complex disorders. A serious problem with association mapping is that population structure can lead to spurious associations between a candidate marker and a phenotype. One common solution has been to abandon case-control studies in favor of family-based tests of association, such as the transmission/disequilibrium test (TDT), but this comes at a considerable cost in the need to collect DNA from close relatives of affected individuals. In this article we describe a novel, statistically valid, method for case-control association studies in structured populations. Our method uses a set of unlinked genetic markers to infer details of population structure, and to estimate the ancestry of sampled individuals, before using this information to test for associations within subpopulations. It provides power comparable with the TDT in many settings and may substantially outperform it if there are conflicting associations in different subpopulations.

Microsatellite variation in Central Africa: an analysis of intrapopulational and interpopulational genetic diversity

As a part of a research project on molecular variation in Central Africa, we have analyzed 10 microsatellites (CD4, CSFO, D3S1358, D18S51, D21S11, F13A1, FES, TH01, TPOX, and VWA) in the Bamileke and Ewondo from Cameroon and the Sanga and Mbenzele Pygmies from the Central African Republic (a total of 390 chromosomes). A statistically significant trend towards heterozygote deficiency was detected in the Mbenzele Pygmies. This was established through the use of powerful exact tests for the Hardy-Weinberg equilibrium. A certain degree of isolation and a small effective size may explain this finding. However, the lack of any substantial reduction in allelic diversity in the Mbenzele does not support the possibility that this group has a smaller effective size in evolutionary terms. A possible explanation based on ethnographic studies suggests that the gene flow from non-Pygmies to Pygmies could have been interrupted only in relatively recent times. The analysis of association between genotypes at pairs of independent loci indicates that the level of subheterogeneity is markedly lower in the Bamileke than in other sampled populations. This may be explained by the combined effect of larger population size, more rigid respect of clanic exogamy, and higher matrimonial mobility of the Bamileke. Finally, we have analyzed interpopulational relationships among our sampled populations and other Central African populations. The results are consistent with a previous study of protein loci (Spedini et al. 1999), which suggests the recent history of the Bamileke and Ewondo has led them to aquire a substantial genetic similarity. Furthermore, the Mbenzele Pygmies diverge from Biaka Pygmies, despite their common origin and geographical proximity. This is probably due to the differentiating effect of genetic drift, which is enhanced by the small effective size of Pygmy populations.

Population structure of the rural communities on the island of Krk (Croatia): A comparison of genetic, cultural, and geographic data

The population structure of seven rural communities from the island of Krk (Croatia) was studied on the basis of its genetic constitution and genetic structure, and different sets of biological (genetic), biocultural (migration), cultural (linguistic), and geographic data. Genetic constitution was estimated from the allelic frequencies of nine tetranucleotide STR loci, three HLA class I loci, five HLA class II loci, and five red blood cell enzyme and serum protein loci. Genetic structure was based on coefficients of genetic diversity G(ST) and genetic kinship R(ST). Population structure was assessed by the structure of mutual relationships among biological, biocultural, cultural, and geographic matrices through Mantel's test, or by relationships among biological, biocultural, cultural measures of distance and/or similarity, and geographic distance with regression analysis of the "isolation by distance" model. One group of genetic markers (red blood enzymes, serum proteins, and HLA polymorphisms) reflects heterogeneity of the rural communities (a reflection of important historical processes and migration patterns), and a second group of genetic markers (STR polymorphisms) reflects the currently observed genetic homogeneity of the rural population. The positive correlations between the matrices of genetic (HLA) and linguistic distances indicate their possible co-evolution under conditions of genetic isolation (low gene flow). Furthermore, the negative correlations between the matrices of genetic (STR) distances and the migration kinship of the second and the third cohort indicate temporal coincidence between genetic constitution and structure assessed from STRs and migration patterns in the period between 1892 and 1940. Finally, the positive correlations between the matrices of genetic (STR) and geographic distances indicate isolation by distance. Am. J. Hum. Biol. 12:509-525, 2000. Copyright 2000 Wiley-Liss, Inc.

Human population expansion and microsatellite variation

Polymorphisms at di-, tri-, and tetranucleotide microsatellite loci have been analyzed in 14 worldwide populations. A statistical index of population expansion, denoted S(k), is introduced to detect historical changes in population size using the variation at the microsatellites. The index takes the value 0 at equilibrium with constant population size and is positive or negative according to whether the population is expanding or contracting, respectively. The use of S(k) requires estimation of properties of the mutation distribution for which we use both family data of Dib et al. for dinucleotide loci and our population data on tri- and tetranucleotide loci. Statistical estimates of the expansion index, as well as their confidence intervals from bootstrap resampling, are provided. In addition, a dynamical analysis of S(k) is presented under various assumptions on population growth or decline. The studied populations are classified as having high, intermediate, or low values of S(k) and genetic variation, and we use these to interpret the data in terms of possible population dynamics. Observed values of S(k) for samples of di-, tri-, and tetranucleotide data are compatible with population expansion earlier than 60,000 years ago in Africa, Asia, and Europe if the initial population size before the expansion was on the order of 500. Larger initial population sizes force the lower bound for the time since expansion to be much earlier. We find it unlikely that bottlenecks occurred in Central African, East Asian, or European populations, and the estimated expansion times are rather similar for all of these populations. This analysis presented here suggests that modern human populations departed from Africa long before they began to expand in size. Subsequently, the major groups (the African, East Asian, and European groups) started to grow at approximately same time. Populations of South America and Oceania show almost no growth. The Mbuti population from Zaire appears to have experienced a bottleneck during its expansion.

Microsatellite structures in the context of human evolution

Six microsatellite - or short tandem repeat (STR) - systems with uniform repetitive sequences (HumTH01, HumCD4, HumFES/FPS, HumF13B, HumTPO, HumLPL) and three compound repeat systems (HumVWA, HumFIBRA, D21S11) were used, including data from the literature, to determine genetic distances among eight populations worldwide. The TH01- and VWA homologous loci in nonhuman primates (chimpanzees, gorillas, orangutans, rhesus monkeys, ring-tailed lemurs) were compared and found to be shorter than in humans. Microsatellites of lower complexity were most efficient for the separation of major ethnic groups. The loci of higher complexity showed a leveling of the diversity differences among populations, which could be attributed to higher mutation rates.

The distribution of human genetic diversity: a comparison of mitochondrial, autosomal, and Y-chromosome data

We report a comparison of worldwide genetic variation among 255 individuals by using autosomal, mitochondrial, and Y-chromosome polymorphisms. Variation is assessed by use of 30 autosomal restriction-site polymorphisms (RSPs), 60 autosomal short-tandem-repeat polymorphisms (STRPs), 13 Alu-insertion polymorphisms and one LINE-1 element, 611 bp of mitochondrial control-region sequence, and 10 Y-chromosome polymorphisms. Analysis of these data reveals substantial congruity among this diverse array of genetic systems. With the exception of the autosomal RSPs, in which an ascertainment bias exists, all systems show greater gene diversity in Africans than in either Europeans or Asians. Africans also have the largest total number of alleles, as well as the largest number of unique alleles, for most systems. GST values are 11%-18% for the autosomal systems and are two to three times higher for the mtDNA sequence and Y-chromosome RSPs. This difference is expected because of the lower effective population size of mtDNA and Y chromosomes. A lower value is seen for Y-chromosome STRs, reflecting a relative lack of continental population structure, as a result of rapid mutation and genetic drift. Africa has higher GST values than does either Europe or Asia for all systems except the Y-chromosome STRs and Alus. All systems except the Y-chromosome STRs show less variation between populations within continents than between continents. These results are reassuring in their consistency and offer broad support for an African origin of modern human populations.

Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations

Single-nucleotide polymorphisms (SNPs) may be extremely important for deciphering the impact of genetic variation on complex human diseases. The ultimate value of SNPs for linkage and association mapping studies depends in part on the distribution of SNP allele frequencies and intermarker linkage disequilibrium (LD) across populations. Limited information is available about these distributions on a genomewide scale, particularly for LD. Using 114 SNPs from 33 genes, we compared these distributions in five American populations (727 individuals) of African, European, Chinese, Hispanic, and Japanese descent. The allele frequencies were highly correlated across populations but differed by >20% for at least one pair of populations in 35% of SNPs. The correlation in LD was high for some pairs of populations but not for others (e.g., Chinese American or Japanese American vs. any other population). Regardless of population, average minor-allele frequencies were significantly higher for SNPs in noncoding regions (20%-25%) than for SNPs in coding regions (12%-16%). Interestingly, we found that intermarker LD may be strongest with pairs of SNPs in which both markers are nonconservative substitutions, compared to pairs of SNPs where at least one marker is a conservative substitution. These results suggest that population differences and marker location within the gene may be important factors in the selection of SNPs for use in the study of complex disease with linkage or association mapping methods.

Population bottlenecks and Pleistocene human evolution

We review the anatomical and archaeological evidence for an early population bottleneck in humans and bracket the time when it could have occurred. We outline the subsequent demographic changes that the archaeological evidence of range expansions and contractions address, and we examine how inbreeding effective population size provides an alternative view of past population size change. This addresses the question of other, more recent, population size bottlenecks, and we review nonrecombining and recombining genetic systems that may reflect them. We examine how these genetic data constrain the possibility of significant population size bottlenecks (i.e., of sufficiently small size and/or long duration to minimize genetic variation in autosomal and haploid systems) at several different critical times in human history. Different constraints appear in nonrecombining and recombining systems, and among the autosomal loci most are incompatible with any Pleistocene population size expansions. Microsatellite data seem to show Pleistocene population size expansions, but in aggregate they are difficult to interpret because different microsatellite studies do not show the same expansion. The archaeological data are only compatible with a few of these analyses, most prominently with data from Alu elements, and we use these facts to question whether the view of the past from analysis of inbreeding effective population size is valid. Finally, we examine the issue of whether inbreeding effective population size provides any reasonable measure of the actual past size of the human species. We contend that if the evidence of a population size bottleneck early in the evolution of our lineage is accepted, most genetic data either lack the resolution to address subsequent changes in the human population or do not meet the assumptions required to do so validly. It is our conclusion that, at the moment, genetic data cannot disprove a simple model of exponential population growth following a bottleneck 2 MYA at the origin of our lineage and extending through the Pleistocene. Archaeological and paleontological data indicate that this model is too oversimplified to be an accurate reflection of detailed population history, and therefore we find that genetic data lack the resolution to validly reflect many details of Pleistocene human population change. However, there is one detail that these data are sufficient to address. Both genetic and anthropological data are incompatible with the hypothesis of a recent population size bottleneck. Such an event would be expected to leave a significant mark across numerous genetic loci and observable anatomical traits, but while some subsets of data are compatible with a recent population size bottleneck, there is no consistently expressed effect that can be found across the range where it should appear, and this absence disproves the hypothesis.

Identifying conservation units within captive chimpanzee populations

One of the primary objectives in the captive management of any endangered primate is to preserve as much as possible the genetic diversity that has evolved and still exists in wild gene pools. The rationale for this is based on the theoretical understanding of the relationship between genetic diversity and fitness in response to selection. There remains little consensus, however, as to the type of genetic data that should be used to monitor captive populations. In order to develop a deeper understanding of the degree and nature of genetic diversity among "wild" chimpanzee gene pools, as well as to determine if one type of genetic data is more useful than others, DNA sequence data were generated at three unlinked, nonrepetitive nuclear loci, one polymorphic microsatellite, and the mitochondrial D-loop for 59 unrelated common and pygmy chimpanzees. The results suggest that: 1) data from nuclear loci can be used to differentiate common chimpanzee subspecies; 2) pygmy chimpanzees may have less genetic diversity than common chimpanzees; 3) shared microsatellite alleles do not always indicate identity by descent; and 4) nonrepetitive loci provide unique insights into evolutionary relationships and provide useful information for captive management programs.

Detecting population expansion and decline using microsatellites

This article considers a demographic model where a population varies in size either linearly or exponentially. The genealogical history of microsatellite data sampled from this population can be described using coalescent theory. A method is presented whereby the posterior probability distribution of the genealogical and demographic parameters can be estimated using Markov chain Monte Carlo simulations. The likelihood surface for the demographic parameters is complicated and its general features are described. The method is then applied to published microsatellite data from two populations. Data from the northern hairy-nosed wombat show strong evidence of decline. Data from European humans show weak evidence of expansion.

A microsatellite survey of cattle from a centre of origin: the Near East

Eight humpless cattle breeds from the Near East, three from Europe, one from West Africa and two zebu breeds from India were screened with 20 microsatellite loci. Breeds from the Near East revealed considerable levels of introgression from zebu cattle, which was apparent most in populations from the East and which declined in populations further West. This nonrandom pattern is suggestive of the introduction of zebu cattle from the East. Notwithstanding the overlay of zebu alleles, it was possible to demonstrate that Near Eastern cattle exhibited significantly higher levels of allelic diversity than breeds from other regions, which is consistent with the view that this region represents a primary domestication centre for Bos taurus cattle. The hypothesis that B. taurus and B. indicus cattle have separate domestic origins is also supported by the survey, a large genetic divergence being apparent between the nonhybrid taurine and zebu groups.

Trinucleotide repeats and other microsatellites in yeasts

Microsatellites are direct tandem DNA repeats found in all genomes. A particular class of microsatellites, called trinucleotide repeats, is responsible for a number of neurological disorders in humans. We review here our current state of knowledge on trinucleotide repeat instability, and discuss the molecular mechanisms that may be involved in trinucleotide repeat expansions leading to fatal diseases in humans. We also present original data on microsatellite distribution in several microbial genomes, and on the use of microsatellites as physical markers to accurately and easily genotype yeast strains.

Genomic views of human history

New tools of genomic analysis shed light on historical puzzles. Migrations of ancient peoples, differences in migration patterns of males and females, historical demography of cultures with ancient roots, and patterns of human genetic diversity are increasingly the focus of integrated analysis by historians, anthropologists, and geneticists.

Genetic diversity of two African and sixteen South American populations determined on the basis of six hypervariable loci

A total of 582 individuals (1,164 chromosomes) from two African, eight African-derived South American, five South American Amerindian, and three Brazilian urban populations were studied at four variable number of tandem repeat (VNTR) and two short tandem repeat (STR) hypervariable loci. These two sets of loci did not show distinct allele profiles, which might be expected if different processes promoted their molecular differentiation. The two African groups showed little difference between them, and their intrapopulational variation was similar to those obtained in the African-derived South American communities. The latter showed different degrees of interpopulation variability, despite the fact that they presented almost identical average degrees of non-African admixture. The F(ST) single locus estimates differed in the five sets of populations, probably due to genetic drift, indicating the need to consider population structure in the evaluation of their total variability. A high interpopulational diversity was found among Amerindian populations in relation to Brazilian African-derived isolated communities. This is probably a consequence of the differences in the patterns of gene flow and genetic drift that each of these semi-isolated groups experienced.

Apolipoprotein E and methylenetetrahydrofolate reductase genetic polymorphisms in relation to other risk factors for cardiovascular disease in UK Caucasians and Black South Africans

Genetic polymorphisms for apolipoprotein E (apo E) and methylenetetrahydrofolate reductase (MTHFR) are believed to modulate risk of coronary heart disease (CHD) acting through regulation of lipid and homocysteine metabolism, respectively. The distributions of apo E and MTHFR alleles in Black South Africans, a population with a low CHD incidence, and UK Caucasians from the Cambridge area, with a higher CHD incidence, were therefore compared. Clinically healthy volunteers (207), including 107 UK Caucasians from the Cambridge area and 100 Black South Africans, participated in the study. Apo E and MTHFR genotypes were determined in all of them. Analyses for serum total cholesterol, LDL cholesterol, HDL cholesterol, triglycerides and plasma fibrinogen were carried out in 65 UK Caucasians and 60 Black South Africans. The apo E epsilon4 allele, which is associated with elevated CHD risk, was present in 48% of Black South Africans compared to 20.8% of Caucasians (P < 0.0001); however, both total and LDL cholesterol levels in Black South Africans were 18-32% lower than in Caucasians with similar apo E genotypes. Hyperhomocysteinemia-causing MTHFR 677T variant was detected in only 20% of Black South Africans (no homozygotes) versus 56% of Caucasians with 12% homozygotes (P<0.0001). Our findings suggest that the potentially unfavourable pattern of apo E allele distribution in Black South Africans does not result in increased CHD incidence due to protection by dietary and/or other life style related factors. The exceptionally low frequency of MTHFR mutant homozygotes in this population suggests that this polymorphism should not be regarded as an important CHD risk factor among Black South Africans.

Use of unlinked genetic markers to detect population stratification in association studies

We examine the issue of population stratification in association-mapping studies. In case-control studies of association, population subdivision or recent admixture of populations can lead to spurious associations between a phenotype and unlinked candidate loci. Using a model of sampling from a structured population, we show that if population stratification exists, it can be detected by use of unlinked marker loci. We show that the case-control-study design, using unrelated control individuals, is a valid approach for association mapping, provided that marker loci unlinked to the candidate locus are included in the study, to test for stratification. We suggest guidelines as to the number of unlinked marker loci to use.

Spatial and temporal distribution of the neutral polymorphisms in the last ZFX intron: analysis of the haplotype structure and genealogy

With 10 segregating sites (simple nucleotide polymorphisms) in the last intron (1089 bp) of the ZFX gene we have observed 11 haplotypes in 336 chromosomes representing a worldwide array of 15 human populations. Two haplotypes representing 77% of all chromosomes were distributed almost evenly among four continents. Five of the remaining haplotypes were detected in Africa and 4 others were restricted to Eurasia and the Americas. Using the information about the ancestral state of the segregating positions (inferred from human-great ape comparisons), we applied coalescent analysis to estimate the age of the polymorphisms and the resulting haplotypes. The oldest haplotype, with the ancestral alleles at all the sites, was observed at low frequency only in two groups of African origin. Its estimated age of 740 to 1100 kyr corresponded to the time to the most recent common ancestor. The two most frequent worldwide distributed haplotypes were estimated at 550 to 840 and 260 to 400 kyr, respectively, while the age of the continentally restricted polymorphisms was 120 to 180 kyr and smaller. Comparison of spatial and temporal distribution of the ZFX haplotypes suggests that modern humans diverged from the common ancestral stock in the Middle Paleolithic era. Subsequent range expansion prevented substantial gene flow among continents, separating African groups from populations that colonized Eurasia and the New World.

A View of Modern Human Origins from Y Chromosome Microsatellite Variation

The idea that all modern humans share a recent (within the last 150,000 years) African origin has been proposed and supported on the basis of three observations. Most genetic loci examined to date have (1) shown greater diversity in African populations than in others, (2) placed the first branch between African and all non-African populations in phylogenetic trees, and (3) indicated recent dates for either the molecular coalescence (with the exception of some autosomal and X-chromosomal loci) or for the time of separation between African and non-African populations. We analyze variation at 10 Y chromosome microsatellite loci that were typed in 506 males representing 49 populations and every inhabited continent and find significantly greater Y chromosome diversity in Africa than elsewhere, find the first branch in phylogenetic trees of the continental populations to fall between African and all non-African populations, and date this branching with the (δμ)2 distance measure to 5800–17,400 or 12,800–36,800 years BP depending on the mutation rate used. The magnitude of the excess Y chromosome diversity in African populations appears to result from a greater antiquity of African populations rather than a greater long-term effective population size. These observations are most consistent with a recent African origin for all modern humans.

A view of modern human origins from Y chromosome microsatellite variation

The idea that all modern humans share a recent (within the last 150, 000 years) African origin has been proposed and supported on the basis of three observations. Most genetic loci examined to date have (1) shown greater diversity in African populations than in others, (2) placed the first branch between African and all non-African populations in phylogenetic trees, and (3) indicated recent dates for either the molecular coalescence (with the exception of some autosomal and X-chromosomal loci) or for the time of separation between African and non-African populations. We analyze variation at 10 Y chromosome microsatellite loci that were typed in 506 males representing 49 populations and every inhabited continent and find significantly greater Y chromosome diversity in Africa than elsewhere, find the first branch in phylogenetic trees of the continental populations to fall between African and all non-African populations, and date this branching with the (deltamu)2 distance measure to 5800-17,400 or 12,800-36,800 years BP depending on the mutation rate used. The magnitude of the excess Y chromosome diversity in African populations appears to result from a greater antiquity of African populations rather than a greater long-term effective population size. These observations are most consistent with a recent African origin for all modern humans.

Fine mapping of a polymorphic CA repeat marker on human chromosome 19 and its use in population studies

A highly polymorphic microsatellite (CA)n-marker (CAct685) previously isolated from human chromosome 19 cosmid library was localized near GPI in 19q13.1. For the fine localization of this marker, the hybridization with chromosome 19-specific cosmid libraries assembled in contigs was used. Polymorphism analysis of the marker in 12 populations of Russia and neighboring countries showed 14 alleles containing from 16 to 30 repeat units. Populations belonging to Indo-European, Uralic and Altaic linguistic families demonstrated a great similarity in allele frequency profiles. Differences between these populations were lower for CAct685 than for classical markers. Allele distribution of CAct685 in a Chukchi population belonging to the Chukchi-Kamchatkan linguistic family differs from those in all other populations, that may be typical for Mongoloid population or reflect an ethnic history of Chukchi as a small population. Thus use of the CAct685 marker seems to be effective for analysis of distant peoples.

X chromosome evidence for ancient human histories

Diverse African and non-African samples of the X-linked PDHA1 (pyruvate dehydrogenase E1 alpha subunit) locus revealed a fixed DNA sequence difference between the two sample groups. The age of onset of population subdivision appears to be about 200 thousand years ago. This predates the earliest modern human fossils, suggesting the transformation to modern humans occurred in a subdivided population. The base of the PDHA1 gene tree is relatively ancient, with an estimated age of 1.86 million years, a late Pliocene time associated with early species of Homo. PDHA1 revealed very low variation among non-Africans, but in other respects the data are consistent with reports from other X-linked and autosomal haplotype data sets. Like these other genes, but in conflict with microsatellite and mitochondrial data, PDHA1 does not show evidence of human population expansion.

Genetic evidence for larger African population size during recent human evolution

Genetic evidence suggests that the long-term average effective size of sub-Saharan Africa is larger than other geographic regions. A method is described that allows estimation of relative long-term regional population sizes. This method is applied to 60 microsatellite DNA loci from a sample of 72 sub-Saharan Africans, 63 East Asians, and 120 Europeans. Average heterozygosity is significantly higher in the sub-Saharan African sample. Expected heterozygosity was computed for each region and locus using a population genetic model based on the null hypothesis of equal long-term population sizes. Average residual heterozygosity is significantly higher in the sub-Saharan African sample, indicating that African population size was larger than other regions during recent human evolution. The best fit of the model is with relative population weights of 0.73 for sub-Saharan Africa, 0.09 for East Asia, and 0.18 for Europe. These results are similar to those obtained using craniometric variation for these three geographic regions. These results, combined with inferences from other genetic studies, support a major role of Africa in the origin of modern humans. It is less clear, however, whether complete African replacement is the most appropriate model. An alternative is an African origin with non-African gene flow. While Africa is an important region in recent human evolution, it is not clear whether the gene pool of our species is completely out of Africa or predominately out of Africa.

Genetic differentiation within and between isolated Algerian subpopulations of Barbary macaques (Macaca sylvanus): evidence from microsatellites

This study of wild-living Algerian Barbary macaques (Macaca sylvanus) was designed to examine genetic variability in subpopulations isolated in residual forest patches, in an attempt to obtain data on the effects of habitat fragmentation. The wild population of this species (estimated at a maximum of 15,000) is vulnerable and this study therefore has direct relevance for conservation measures. Data from five microsatellite loci were analysed for 159 individuals from nine different groups living in four isolates in Algeria. Genetic polymorphism was found to be relatively high (4-12 alleles per locus) compared with other genetic markers used in previous studies of this species; mean expected heterozygosity was 65%. The four isolates are all genetically distinct (FST = 0; P < 0.001). Indeed, the results suggest that dispersal is limited even between some social groups within a single isolate. Genetic distances based on models not assuming stepwise mutation (FST and chord distance) gave very similar results and are highly correlated with geographical distances within one isolate but not between isolates. This may indicate that isolation by distance exerts a significant influence within an isolate but that genetic drift prevails between the four isolates. After allowing for variation in sample size, we found no evidence of reduced allelic diversity within small isolates that may have been separated for 250 years or more. The surviving population of Algerian Barbary macaques taken as a whole still shows marked variability in microsatellite alleles, but maintenance of genetic variability over the long term will surely require effective protection of all isolates.

Genealogy reconstruction from short tandem repeat genotypes in an Amazonian population

We have reconstructed partial genealogies in a sample of 44 SW Amazonian Rondonian Surui, in which 45 dinucleotide short tandem repeat polymorphisms had previously been typed. The genotypes of 488 pairs of individuals having an age difference of 13 or greater were compared, and parentage was excluded if a pair failed to share an allele at more than one locus. In order to test the power of this method, we computed the expected distribution of the number of exclusionary loci for such pairs of unrelated individuals, as well as that for individuals with different degrees of relatedness, and compared it to the observed distribution. We estimated that the pairs compared contained approximately 20% of individual pairs with a first-cousin relation or closer. A total of 25 pairs were identified as possible parent-child. In three instances, we could identify two or more children having a common parent; we computed a relatedness coefficient in order to establish whether the children were full or half sibs. The genealogies inferred show that instances of polygyny and polyandry (or, alternatively, serial mating), in addition to apparent monogamy, can be found among the Surui. The Surui sample can be used as a model for paleoanthropological populations, in which the determination of relatedness can provide further insights into the social structure of past populations. We estimate that, depending on the history of the populations and the degree of inbreeding, 10-20 highly informative nuclear loci should be typed in order to infer genealogies with acceptable confidence.

Population substructure and isolation by distance in three continental regions

Isolation by distance and divergence from a shared population history are two sources of population substructure. Isolation by distance erases population history as populations approach migration-drift equilibrium, while diverging populations descended from a single ancestral population will accumulate genetic differences with time. Here I investigate how much of the worldwide genetic diversity from Jorde et al.'s ([1997] Proc. Natl. Acad. Sci. USA 94:3100-3103) 60 tetranucleotide short tandem repeat (STR) data can be explained by isolation by distance. I use Slatkin's measure of population substructure, R(ST), principal components analyses, and Mantel tests to investigate the pattern of genetic diversity at both the intercontinental and intracontinental levels. Geographic distance accounts for almost 60% of world-wide interpopulation genetic relationships. Within continents, the correlations are less, although not significantly so because of wide confidence intervals. These results suggest that populations have not reached migration-drift equilibrium and that there is information in STR data to reconstruct population history. The level of population substructure worldwide is consistent with previous observations, but at the intracontinental level substructure is less. When one examines diversity against distance from the centroid, one sees excess heterozygosity in Africa, a pattern also noted by Stoneking et al. ([1998] Genome Research 7:1061-1071). A larger effective population size in Africa could explain the excess diversity. Greater gene flow in Africa is an unlikely explanation because the African R(ST) value is slightly larger than the Asian and European R(ST)s, pointing to less gene flow and greater substructure among African populations. Furthermore, there are differences in patterns between heterozygosity and allele size variance. Heterozygosity has a higher correlation with distance from the centroid than does allele size variance, and this may reflect demographic history. Kimmel et al. ([1998] Genetics 148:1921-1930) have shown that after a population expansion heterozygosity returns to equilibrium more quickly than does allele size variance. The contrasting patterns between heterozygosity and allele size variance may reflect different times after an expansion. However, simulations and further work need to be done to more thoroughly investigate the possibility that these data reflect population expansions.

Allele frequencies in a worldwide survey of a CA repeat in the first intron of the CFTR gene

A dinucleotide CA repeat within intron 1 of the CFTR gene has recently been identified. We have determined the allele frequencies of this polymorphism in samples from 18 populations covering all major geographical areas, with a total of 1,816 chromosomes. When considering allele distributions, African populations presented a wider range of alleles than other populations and also presented higher expected heterozygosities. Analysis of molecular variance showed that 8.04% of the genetic variance in this locus could be attributed to differences among populations. We concluded that the polymorphism in the CA repeat in intron 1 of the CFTR gene is highly informative in populations from all geographical regions of the world. Thus, it can be applied to family studies of unknown mutations causing cystic fibrosis (CF) and can provide valuable information for genetic counseling. Moreover, its analysis should be included in the haplotypic analysis of known CF mutations.

COMING TO TERMS WITH HUMAN VARIATION

▪ Abstract  Genetics has become the major tool of the life sciences. This is driven partly by technology, and partly by the belief that genes are the ultimate units of biomedical or evolutionary information. The search for variation associated with disease has motivated the Human Genome Project to construct a detailed road map of the entire set of human genetic material, and some additional form of globally representative human genome diversity resource has been proposed for the anthropological purposes of reconstructing human population history. Any such resource raises complex societal and ethical issues as well as scientific ones. However, the amount and complexity of genetic variation has frustrated hopes for simple genetic answers to important biomedical or anthropological questions, and a consequent converging of these differing interests suggests that developing a genetic variation resource will be important in many disciplines.

Mitochondrial and nuclear genetic relationships among Pacific Island and Asian populations

Mitochondrial and autosomal short tandem-repeat (STR) genetic distances among 28 Pacific Island and Asian populations are significantly correlated (r=.25, P<.01) but describe distinct patterns of relationships. Maternally inherited-mtDNA data suggest that Remote Oceanic Islanders originated in island Southeast Asia. In contrast, biparental STR data reveal substantial genetic affinities between Remote Oceanic Islanders and Near Oceanic populations from highland Papua New Guinea and Australia. The low correlation between maternal and biparental genetic markers from the same individuals may reflect differences in genome-effective population sizes or in sex-biased gene flow. To explore these possibilities, we have examined genetic diversity, gene flow, and correlations among genetic, linguistic, and geographic distances within four sets of populations representing potential geographic and cultural spheres of interaction. GST estimates (a measure of genetic differentiation inversely proportional to gene flow) from mtDNA sequences vary between 0.13 and 0.39 and are typically five times greater than GST estimates from STR loci (0.05-0.08). Significant correlations (r>.5, P<.05) between maternal genetic and linguistic distances are coincident with high mtDNA GST estimates (>0.38). Thus, genetic and linguistic distances may coevolve, and their correspondence may be preserved under conditions of genetic isolation. A significant correlation (r=.65, P<.01) between biparental genetic and geographic distances is coincident with a low STR GST estimate (0.05), indicating that isolation by distance is observed under conditions of high nuclear-gene flow. These results are consistent with an initial settlement of Remote Oceania from island Southeast Asia and with extensive postcolonization male-biased gene flow with Near Oceania.

Late Pleistocene human population bottlenecks, volcanic winter, and differentiation of modern humans

The "Weak Garden of Eden" model for the origin and dispersal of modern humans (Harpending et al., 1993) posits that modern humans spread into separate regions from a restricted source, around 100 ka (thousand years ago), then passed through population bottlenecks. Around 50 ka, dramatic growth occurred within dispersed populations that were genetically isolated from each other. Population growth began earliest in Africa and later in Eurasia and is hypothesized to have been caused by the invention and spread of a more efficient Later Stone Age/Upper Paleolithic technology, which developed in equatorial Africa. Climatic and geological evidence suggest an alternative hypothesis for Late Pleistocene population bottlenecks and releases. The last glacial period was preceded by one thousand years of the coldest temperatures of the Later Pleistocene (approximately 71-70 ka), apparently caused by the eruption of Toba, Sumatra. Toba was the largest known explosive eruption of the Quaternary. Toba's volcanic winter could have decimated most modern human populations, especially outside of isolated tropical refugia. Release from the bottleneck could have occurred either at the end of this hypercold phase, or 10,000 years later, at the transition from cold oxygen isotope stage 4 to warmer stage 3. The largest populations surviving through the bottleneck should have been found in the largest tropical refugia, and thus in equatorial Africa. High genetic diversity in modern Africans may thus reflect a less severe bottleneck rather than earlier population growth. Volcanic winter may have reduced populations to levels low enough for founder effects, genetic drift and local adaptations to produce rapid population differentiation. If Toba caused the bottlenecks, then modern human races may have differentiated abruptly, only 70 thousand years ago.

Detection of internal repeats: how common are they?

The exponential growth in the amount of genomic data published in recent years has led to increased efforts in analysing genomes for the presence of repeated sequences, which has in turn fostered the development of novel repeat recognition methods. This has resulted in a deepened understanding of the importance and abundance of protein and nucleotide repeats. In the past year, a shift in focus has taken place--from the significance of repeats to protein structure and function, mostly at the protein domain level, to the implication of generally much shorter repeated fragments in genetic diseases and protein malfunctioning.

Signatures of population expansion in microsatellite repeat data

To examine the signature of population expansion on genetic variability at microsatellite loci, we consider a population that evolves according to the time-continuous Moran model, with growing population size and mutations that follow a general asymmetric stepwise mutation model. We present calculations of expected allele-size variance and homozygosity at a locus in such a model for several variants of growth, including stepwise, exponential, and logistic growth. These calculations in particular prove that population bottleneck followed by growth in size causes an imbalance between allele size variance and heterozygosity, characterized by the variance being transiently higher than expected under equilibrium conditions. This effect is, in a sense, analogous to that demonstrated before for the infinite allele model, where the number of alleles transiently increases after a stepwise growth of population. We analyze a set of data on tetranucleotide repeats that reveals the imbalance expected under the assumption of bottleneck followed by population growth in two out of three major racial groups. The imbalance is strongest in Asians, intermediate in Europeans, and absent in Africans. This finding is consistent with previous findings by others concerning the population expansion of modern humans, with the bottleneck event being most ancient in Africans, most recent in Asians, and intermediate in Europeans. Nevertheless, the imbalance index alone cannot reliably estimate the time of initiation of population expansion.

Heterogeneity of microsatellite mutations within and between loci, and implications for human demographic histories

Microsatellites have been widely used to reconstruct human evolution. However, the efficient use of these markers relies on information regarding the process producing the observed variation. Here, we present a novel approach to the locus-by-locus characterization of this process. By analyzing somatic mutations in cancer patients, we estimated the distributions of mutation size for each of 20 loci. The same loci were then typed in three ethnically diverse population samples. The generalized stepwise mutation model was used to test the predicted relationship between population and mutation parameters under two demographic scenarios: constant population size and rapid expansion. The agreement between the observed and expected relationship between population and mutation parameters, even when the latter are estimated in cancer patients, confirms that somatic mutations may be useful for investigating the process underlying population variation. Estimated distributions of mutation size differ substantially amongst loci, and mutations of more than one repeat unit are common. A new statistic, the normalized population variance, is introduced for multilocus estimation of demographic parameters, and for testing demographic scenarios. The observed population variation is not consistent with a constant population size. Time estimates of the putative population expansion are in agreement with those obtained by other methods.

Using mitochondrial and nuclear DNA markers to reconstruct human evolution

Molecular genetic-data have greatly improved our ability to test hypotheses about human evolution. During the past decade, a large amount of nuclear and mitochondrial data have been collected from diverse human populations. Taken together, these data indicate that modern humans are a relatively young species. African populations show the largest amount of genetic diversity, and they are the most genetically divergent population. Modern human populations expanded in size first on the African continent. These findings support a recent African origin of modern humans, but this conclusion should be tempered by the possible effects of factors such as gene flow, population size differences, and natural selection.

Demographic history and linkage disequilibrium in human populations

In the human genome, linkage disequilibrium (LD)--the non-random association of alleles at chromosomal loci--has been studied mainly in regions surrounding disease genes on affected chromosomes. Consequently, little information is available on the distribution of LD across anonymous genomic regions in the general population. However, demographic history is expected to influence the extent of overall LD across the genome, so a population that has been of constant size will display higher levels of LD than a population that has expanded. In support of this, the extent of LD between anonymous loci on chromosome 4 in chimpanzees (as a model of a population of constant size) has been compared to that in Finns (as a model of an expanded population; refs 8,9) and found to exhibit more LD than in the latter population. In Europe, studies of mitochondrial (mt) DNA sequences have suggested that most populations have experienced expansion, whereas the Saami in northern Fenno-Scandinavia have been of constant size (Table 1). Thus, in northern Europe, populations with radically different demographic histories live in close geographic proximity to each other. We studied the allelic associations between anonymous microsatellite loci on the X chromosome in the Saami and neighbouring populations and found dramatically higher levels of LD in the Saami than in other populations in the region. This indicates that whereas recently expanded populations, such as the Finns, are well suited to map single disease genes affected by recent mutations, populations that have been of constant size, such as the Saami, may be much better suited to map genes for complex traits that are caused by older mutations.

Out of Africa? What do genes tell us?

Genetic diversity patterns in nuclear versus mitochondrial systems and in low versus high mutation rate systems do not support the hypothesis of a recent African origin for all of humanity following a split between Africans and non-Africans 100,000 years ago, nor do genetic distance data. Geographical analyses of nuclear and mitochondrial gene trees do not support the hypothesis of a recent global replacement of humans coming out of Africa, although a local replacement event in Europe is indicated by these analyses and recent studies on Neandertal DNA. The gene tree analyses instead indicate that genetic interchanges have ensured that all of humanity has evolved as a single evolutionary lineage with no major splits.

Evidence for a possible Asian origin of YAP+ Y chromosomes

Images