Primary research | Open | Published:
Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation
Human Genomicsvolume 2, Article number: 81 (2005)
Understanding the distribution of human genetic variation is an important foundation for research into the genetics of common diseases. Some of the alleles that modify common disease risk are themselves likely to be common and, thus, amenable to identification using gene-association methods. A problem with this approach is that the large sample sizes required for sufficient statistical power to detect alleles with moderate effect make gene-association studies susceptible to false-positive findings as the result of population stratification [1, 2]. Such type I errors can be eliminated by using either family-based association tests or methods that sufficiently adjust for population stratification [3–5]. These methods require the availability of genetic markers that can detect and, thus, control for sources of genetic stratification among populations. In an effort to investigate population stratification and identify appropriate marker panels, we have analysed 11,555 single nucleotide polymorphisms in 203 individuals from 12 diverse human populations. Individuals in each population cluster to the exclusion of individuals from other populations using two clustering methods. Higher-order branching and clustering of the populations are consistent with the geographic origins of populations and with previously published genetic analyses. These data provide a valuable resource for the definition of marker panels to detect and control for population stratification in population-based gene identification studies. Using three US resident populations (European-American, African-American and Puerto Rican), we demonstrate how such studies can proceed, quantifying proportional ancestry levels and detecting significant admixture structure in each of these populations.
Substantial progress has been made using genetic markers to elucidate the evolutionary histories of populations, yet this work has primarily been accomplished using large numbers of individuals and small numbers of genetic markers [6, 7]. More recently, studies screening large numbers of markers have demonstrated their effectiveness at clarifying more subtle patterns of population stratification [8, 9]. Such studies can facilitate the exploration of the genetic structure that may exist among and within populations and also provide a valuable source of ancestry informative markers (AIMs) to quantify and adjust for this structure in gene-identification studies.
Here, we analysed 11,555 single nucleotide polymorphism (SNP) markers in 12 population samples using a new microarray-genotyping platform called Whole Genome Sampling Amplification (WGSA; Affymetrix, Santa Clara, CA). Populations were selected to represent a broad spectrum of world variation (Table 1). Four populations stand out as having similarly elevated heterozygosity (Burunge, Spanish, Indian and Altaian). Four groups (Nahua, Quechua, Nasioi and Mbuti) have lower levels of variability, while the two East Asian populations and the Mende are intermediate. These results are largely consistent with expectations for populations known to have experienced restrictions in population size (eg Mbuti, Nasioi, Nahua and Quechua), reducing levels of genetic variability relative to other populations;[9, 10] however, ascertainment bias in terms of the population (s) in which markers were first discovered precludes making strong statements about differences in variability using SNP data . In addition, ascertainment bias -- such as that resulting from both a limited representation of populations and small numbers of individuals in discovey panels-- can lead to deviations in linkage disequilibrium estimates, artefactually elevated FST levels (see Ronald and Akey in this issue of Human Genomics) and higher derived allele frequencies in non-African than in African populations . Despite the general importance of considering ascertainment bias on a number of population genetic parameter estimates, there is no evidence or theory that predicts problems from ascertainment bias on estimates of measures of individual relatedness or deviations from Hardy - Weinberg equilibrium (HWE). Randomly mating populations are expected to show genotype frequencies that are consistent with HWE expectations. The results of tests for HWE are presented as the proportion of loci that have deviations from equilibrium expectations (Table 1). Some populations show slightly higher or lower proportions of significant results. The most notable deviation is seen when all the populations are combined, at which point over half (56 per cent) of the SNPs show significant HWE deviations. These deviations highlight the importance of taking population structure into account in gene-association studies. Extensive admixture structure is created by combining these samples, leading to HWE deviations and, presumably, allelic associations among unlinked markers at loci showing large frequency differences across populations .
The proportion of the total genetic variation due to differences among populations was estimated using FST. Figure 1 shows a histogram of the FST distribution, with the autosomal SNPs plotted separately from the X-linked SNPs. The average level of FST for autosomal SNPs (0.148) is within the range of previously published FST estimates (5 - 15 per cent), confirming the well-known fact that most variability in human populations is observed within populations [4–6, 14, 15]. The average FST observed for the X-chromosomal SNPs (0.224) is substantially higher than that for the autosomal SNPs (p < 0:0001), which is consistent with both the smaller effective population size and the higher levels of natural selection for X-chromosome genes [16–18]. It is also notable that these distributions are not described well by averages, since they are highly skewed and have long tails, highlighting the fact that unlinked loci can have different evolutionary histories .
For calculating heterozygosity and FST, population divisions were assumed to be known and individuals were grouped using ethnic and geographical information. Given the large number of markers in our dataset, population genetic analyses can be performed at the level of the individual, making no presumption of group membership [18, 20]. Two methods were used to investigate clustering among individuals: neighbour-joining trees  and principal coordinates (PCs) analysis, using the allele-sharing distance (ASD) for all pairwise combinations of individuals. Figure 2 shows a neighbour-joining tree of individuals, constructed with the ASD measure matrix, using 11,078 autosomal SNPs. The root of the tree, based on the combined ape out-group, is located between the Mende and Mbuti. This supports an African origin for modern humans. The next group to diverge from the main trunk is the East African Burunge. Most populations have population-specific branches of substantial length, the largest being the Melanesians and the indigenous Americans. In addition to the Burunge, the South Asian Indians and the Altaians have relatively short population-specific branches, consistent with gene flow between these groups and other populations. The largest internal branch separates the three African from the non-African populations, and the next group to diverge is the Spanish, followed by the South Asian Indians. No clear separation of the upper and lower caste populations is seen here (but see Figure 3).
Although trees provide a useful means of illustrating relationships among populations or individuals, they are limited by the assumption of bifurcating topologies. PC analysis is an alternative analytical method, which lacks this assumption. Figure 3a shows the first three PC axes for all populations. As with the tree, individuals from one population cluster tightly, to the exclusion of individuals from other populations. The first PC axis shows a separation of the African and non-African populations, with the Burunge being closer to the non-Africans than either of the other two African populations. The second PC axis shows the indigenous Americans and Melanesians to be on opposite sides of the axis. On the tree, the two indigenous American populations are separated into monophyletic clusters, while the PC analysis shows overlapping clusters. When focusing on the Eurasian populations (Figure 3b), there is a clinal relationship across all three PC axes for these populations, which is consistent with their geographic positions from Spanish in the lower left to Japanese in the upper right. Notable are the near separation of the Indian sample into lower and upper caste, with the upper caste individuals positioned closer to the Spanish [23, 24]. Additionally, the Altaians are intermediate between the East Asians and the Europeans, a finding that is consistent with Y-chromosomal studies showing Central Asian origins for components of the European gene pool .
Another way of exploring the PC analysis results is to examine the pairwise plots of the PC components. As the first four components were significant using the broken stick test, not all can be plotted in three-dimensional space. The six possible pairwise plots are presented in Figures 4a-f.
In addition to these 12 geographically well-defined population samples, we have analysed three cosmopolitan samples collected in the USA (African-Americans, European-Americans and Puerto Ricans). These populations are known to have been subject to both within-continent and among-continent admixture in the recent past. We estimated the individual biogeographical ancestry levels for each person in these three samples. These maximum likelihood estimates of proportional ancestry (Figure 5) show a greater tendency for the European-American subjects to cluster together, by comparison with the other two population samples. The African-Americans and Puerto Ricans both show relatively high levels of variability in individual ancestry levels, with most of the non-African ancestry in the African-Americans being from Europe. The Puerto Ricans show some individuals with more indigenous American ancestry, as well as substantial West African ancestry.
Another test for the presence of admixture structure is based on correlations in individual ancestry indices calculated from independent (unlinked) panels of markers . We tested for significant correlations using two types of individual indices, PCs (Table 2) and biogeographical ancestry (Table 3), calculated separately from the even and odd chromosomal SNPs. To do this, we divided the SNPs into two groups; all of the SNPs on even chromosomes in one group and all of the SNPs on odd chromosomes in the other group. Unless there is structure which is related to the axes of ancestry measured by these indices within a population, no significant relationship between the two estimates is expected . The correlation results on the PC components show that only three (upper caste Indian, Altaian and Nasioi) of the 12 world populations show evidence of population structure. The combined Indian sample (upper caste and lower caste together) also shows significant correlations, while the combined East Asian (Japanese and Chinese) population does not. Alternatively, all three cosmopolitan samples tested (African-American, European-American and Puerto Rican) show significant correlations between the even and odd chromosome PC analyses. Significant correlations are also seen for the estimates of biogeographical ancestry in these three populations (Table 3). It is notable that, not only are there high correlations in the African-American and Puerto Rican samples, but also in the European-American sample, indicating the presence of admixture structure in a population generally assumed to be homogeneous .
The large number of markers used in these analyses has provided an opportunity to assess genetic variation at the level of the individual in a number of populations from around the world. These multilocus genotype data on a large panel of SNPs provide a new level of resolution in the distribution of variation within and among populations. Individuals cluster into groups comprising other individuals from their own or closely-related populations when diverse groups from around the world are analysed [28, 29]. By contrast, when samples from the more cosmopolitan US resident populations are analysed, clustering patterns are less discrete. Substantial levels of variation in ancestry are observed within the African-American and Puerto Rican samples, while smaller, but significant, admixture structure is evident in the European-American sample. Whether the admixture structure in the European-American sample is the result of intraor intercontinental gene flow is an important outstanding question. Thus, although discrete clustering of individuals may be useful in describing some of the variation in diverse, well-defined population samples, continuous measures -- such as biogeographical ancestry or PC indices -- are required to describe the same axes of population structure in populations that have experienced recent admixture.
WGSA technology was used to genotype individuals in this study using the GeneMapping 10 K Array Xba 131 (Affymetrix Inc., Santa Clara, CA). Details of this method have been published elsewhere;[29, 30] in brief, fractions of the genome are obtained by restriction enzyme (Xba I) digestion of genomic DNA, ligated with adaptors and subsequently amplified with a universal primer that is directed to the linker. The amplified target (a smear of polymerase chain reaction products of 400 to 800 base pairs [bps] in length) is fragmented, labelled with terminal transferase and biotin-ddATP and hybridised overnight to synthetic microarrays [31, 32]. Genotypes are called by interpreting signals from allele-specific probes using a model-based algorithm. The accuracy of this method is in excess of 99.5 per cent. SNPs were chosen from The SNP Consortium (TSC) database on the basis of their predicted location on 400 -800 bp fragments generated by in silico digestion of human genome sequences with various restriction enzymes. Predicted SNPs were then assayed against a panel of 108 individuals from diverse populations. If two individuals were observed with each of the three genotypes, and the clustering patterns were acceptable, the SNP was considered to be confirmed and retained as part of the panel.
The population samples used in this study were collected under Internal Review Board approvals from the various institutions involved. The Mbuti population samples were collected in the Ituri Forest, the Mende samples from Sierra Leone. The Cushitic-speaking Burunge samples were collected in Tanzania but are thought to be of Ethiopian descent. The Spanish samples were collected in Valencia in Eastern Spain. The Nasioi were collected in Bougainville, Melanesia. The Altaian samples were collected in Siberia, Russia. The upper and lower caste groups were both sampled from Vishakapatnam, Andhra Pradesh, India. The Chinese (NA17011-NA17020) and Japanese (NA17051-NA17060) samples are from US residents, curated at the Coriell Institute. Quechua were sampled in Lima (n = 9) or Cerro de Pasco, Peru, at 4,338 meters (n = 11). In the former case, the subjects were highland natives, as both parents and grandparents were born on the Altiplano. Quechua subjects were selected to represent a subgroup of subjects with the lowest possible European admixture from a larger total sample of n = 71. Similarly, the Nahua, who were sampled in the city of Tlapa, Guerrero, Mexico, were also selected as a subset of individuals showing low European ancestry, as measured with an independent set of markers. African-Americans (subset of 42 from NA17100-17199) and European-Americans (subset of 42 from NA17200-NA17285) are represented by samples curated at the Coriell Institute. For these analyses, one individual initially classified as 'Caucasian' (Coriell Institute Cat# NA17205) was excluded from the European-American sample, as he/she clusters with South Asians, which, in combination with the lack of monophyletic clustering of South Asians and Spanish in this study, highlights the inappropriateness of the category 'Caucasian' in biomedical research. The Puerto Ricans are women born in Puerto Rico and living in New York City at the time of data collection.
FST was calculated using Weir and Cockerham's unbiased estimator . Pairwise individual genetic distances were estimated using the ASD . The tree of individuals, based on the ASD distance, was constructed using the neighbour-joining method, using the Molecular Evolutionary Genetics Analysis software package (MEGA version 2.1) . The PC analysis was carried out using NTSYS software (Rohlf, F. J. , NTSYS-pc version 1.70). The statistical significance of PC axes was determined using the broken stick model, resulting in four significant axes. These axes together explain 23.6 per cent of the total variation (12.6 per cent by the first axis, 5.5 per cent by the second, 3.5 per cent by the third and 1.9 per cent by the fourth). All pairwise PC axis plots for these four axes are presented in the online supplementary information. The STRUCTURE 2.0  computer program was used to infer the presence of genetic structure in the sample. The analysis was performed both with and without the admixture model for K = 2 to K = 6, the model previously having been determined to show the highest posterior probabilities for these data. A total of 25,000 simulation iterations were run for the burn-in period; 75,000 additional iterations were run to get parameter estimates. Biogeographical ancestry estimates were calculated for the 42 African-American subjects, the 41 European-American subjects and the 20 Puerto Rican women, using the maximum likelihood algorithm previously described, whereby the allele frequencies for the three parental populations were taken to be indigenous American (Nahua and Quechua averaged together), West African (Mende) and European (Spanish). For testing the correlation between subsets of markers, autosomal SNPs were divided by chromosome into those on odd and even chromosomes, respectively. Spearman's correlation coefficient was then calculated for four significant PCs and the three ancestral components between using the odd and even estimates of these statistics. As has been demonstrated, estimates -- even for highly admixed populations -- will be uncorrelated unless there is substantial non-random mating in the population that is related to ancestry .
Risch N, Burchard E, Ziv E, Tang H: Categorization of humans in biomedical research: Genes, race and disease. Genome Biol. 2002, 3: 1-12.
Kittles RA, Chen W, Panguluri RK, et al: CYP3A4-V and prostate cancer in African Americans: Causal or confounding association because of population stratification?. Hum Genet. 2002, 110: 553-560. 10.1007/s00439-002-0731-5.
Hoggart CJ, Parra EJ, Shriver MD, et al: Control of confounding of genetic associations in stratified populations. Am J Hum Genet. 2003, 72: 1492-1504. 10.1086/375613.
Devlin B, Roeder K: Genomic control for association studies. Biometrics. 1999, 55: 997-1004. 10.1111/j.0006-341X.1999.00997.x.
Pritchard JK, Stephens M, Rosenberg NA, Donnelly P: Association mapping in structured populations. Am J Hum Genet. 2000, 67: 170-181. 10.1086/302959.
Cavalli-Sforza LL, Menozzi P, Piazza A: 1994, The History and Geography of Human Genes, Princeton University Press, Princeton, NJ
Ke Y, Su B, Song X, et al: African origin of modern humans in East Asia: A tale of 12,000 Y chromosomes. Science. 2001, 292: 1151-1153. 10.1126/science.1060011.
Rosenberg NA, Pritchard JK, Weber JL, et al: Genetic structure of human populations. Science. 2002, 298: 2381-2385. 10.1126/science.1078311.
Bamshad MJ, Wooding S, Watkins WS, et al: Human population genetic structure and inference of group membership. Am J Hum Genet. 2003, 72: 578-589. 10.1086/368061.
Watkins WS, Rogers AR, Ostler CT, et al: Genetic variation among world populations: Inferences from 100 Alu insertion polymorphisms. Genome Res. 2003, 13: 1167-1618.
Rodgers AR, Jorde LB: Ascertainment bias in estimates of average heterozygosity. Am J Hum Genet. 1996, 58: 1033-1041.
Nielsen R, Signorovitch J: Correcting for ascertainment biases when analyzing SNP data: Applications to the estimation of linkage disequilibrium. Theor Popul Biol. 2003, 63: 245-255. 10.1016/S0040-5809(03)00005-4.
Mountain JL, Cavalli-Sforza LL: Inference of human evolution through cladistic analysis of nuclear DNA restriction polymorphisms. Proc Natl Acad Sci USA. 1994, 91: 6515-6519. 10.1073/pnas.91.14.6515.
Bowcock AM, Ruiz-Linares A, Tomfohrde J, et al: High resolution of human evolutionary trees with polymorphic microsatellites. Nature. 1994, 368: 455-457. 10.1038/368455a0.
Lewontin R: The apportionment of human diversity. Evol Biol. 1972, 6: 381-398.
Charlesworth B, Coyne JA, Barton NH: The relative rates of evolution of sex chromosomes and autosomes. Am Nat. 1987, 130: 113-149. 10.1086/284701.
Payseur BA, Nachman MW: Natural selection at linked sites in humans. Gene. 2002, 300: 31-42. 10.1016/S0378-1119(02)00849-1.
Shriver MD, Kennedy GC, Parra EJ, et al: The genomic distribution of human population substructure in four populations using 8525 SNPs. Hum Genomics. 2004, 1: 274-286.
Cavalli-Sforza LL: Population structure and human evolution. Proc R Soc Lond B Biol Sci. 1966, 164: 362-379. 10.1098/rspb.1966.0038.
Mountain J, Cavalli-Sforza LL: Multilocus genotypes, a tree of individuals and human evolutionary history. Am J Hum Genet. 1997, 61: 705-718. 10.1086/515510.
Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Chakraborty R, Jin L: A unified approach to study hypervariable polymorphisms: Statistical considerations of determining relatedness and population distances. Edited by: Pena SDJ, Jefferys AJ, Epplen J, Chakraborty R. 1993, DNA Fingerprinting: Current State of the Science, EXS, 67: 153-175 Birkhauser, Basel, Switzerland,.
Basu A, Mukherjee N, Roy S, et al: Ethnic India: A genomic view, with special reference to peopling and structure. Genome Res. 2003, 13: 2277-2290. 10.1101/gr.1413403.
Bamshad M, Kivisild T, Watkins WS, et al: Genetic evidence on the origins of Indian caste populations. Genome Res. 2001, 11: 994-1004. 10.1101/gr.GR-1733RR.
Semino O, Passarino G, Oefner PJ, et al: The genetic legacy of Paleolithic Homo sapiens sapiens in extant Europeans: A Y chromosome perspective. Science. 2000, 290: 1155-1159. 10.1126/science.290.5494.1155.
Shriver MD, Parra EJ, Dios S, et al: Skin pigmentation, biogeographical ancestry, and admixture mapping. Hum Genet. 2003, 112: 387-399.
Wacholder S, Rothman N, Caporaso N: Counterpoint: Bias from population stratification is not a major threat to the validity of conclusions from epidemiological studies of common polymorphisms and cancer. Cancer Epidemiol Biomarkers Prev. 2002, 11: 513-520.
Mountain J, Cavalli-Sforza LL: Inference of human evolution through cladistic analysis of nuclear DNA restriction polymorphisms. Proc Natl Acad Sci USA. 1994, 91: 6515-6519. 10.1073/pnas.91.14.6515.
Wilson JF, Weale ME, Smith AC, et al: Population genetic structure of variable drug response. Nat Genet. 2001, 29: 265-269. 10.1038/ng761.
Kennedy GC, Matsuzaki H, Dong S, et al: Large-scale genotyping of complex DNA. Nat Biotechnol. 2003, 21: 1233-1237. 10.1038/nbt869.
Matsuzaki H, Loi H, Dong S, et al: Parallel genotyping of over 10,000 SNPs using a one-primer assay on a high-density oligonucleotide array. Genome Res. 2004, 14: 414-425. 10.1101/gr.2014904.
Chee M, Yang R, Hubbell E, et al: Accessing genetic information with high-density DNA arrays. Science. 1996, 274: 610-614. 10.1126/science.274.5287.610.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population substructure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: Molecular Evolutionary Genetics Analysis software. Bioinformatics. 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244.
Parra EJ, Kittles RA, Shriver MD: Implications of correlations between skin color and genetic ancestry for biomedical research. Nat Genet. 2004, 36: S54-S60. 10.1038/ng1440.
We would like to thank Kateryna Makova, Ken Weiss, Anne Buchanan and S. Malia Fullerton for helpful comments and discussions on this paper. We would also like to acknowledge the trust and generosity of the people and populations who contributed the DNA samples on which this study is based. Additionally, for the collection of the Spanish population sample, we acknowledge Jose Americo Montoro from the blood transfusion centre in Valencia (Servei Valencia de Salut). This work was supported in part by grants: NIH/NHGRI (HG02154) to MDS; NJH/NIA (P30AG/NR15294-01) to JRF; NSF (SBR-9514733 and SBR-9818215) to LBJ;