Genetic diversity and admixture patterns in Indian populations

India is a diverse land whose population holds the history of waves of human dispersal. Recent studies suggest two major ancestral contributions to most of the Indian sub-populations. However, present day Indians are thought to contain huge genetic diversity derived consequent to multiple cultural, linguistic and geographical variations. Genome-wide survey of individuals from current North (N-In) and South (S-In) India along with populations from HapMap Phase III and Indian sub-populations from HUGO Pan-Asian SNP Consortium is performed. Multivariate analysis (MDS and PCA) was carried out after merging data from the current study and other consortia. Indian sub-populations clustered separately from populations of major global geographical regions in MDS and PCA in a loose agglomeration except for two Indian subpopulations clustering near far eastern populations. F_st values indicated diversity among Indian sub-populations which was substantiated by STRUCTURE analysis suggesting the possibility of additional admixture events.     

The Austroasiatic Munda population from India and Its enigmatic origin: a HLA diversity study

The Austroasiatic linguistic family disputes its origin between two geographically distant regions of Asia, India, and Southeast Asia, respectively. As genetic studies based on classical and gender-specific genetic markers provided contradictory results to this debate thus far, we investigated the HLA diversity (HLA-A, -B, and -DRB1 loci) of an Austroasiatic Munda population from Northeast India and its relationships with other populations from India and Southeast Asia. Because molecular methods currently used to test HLA markers often provide ambiguous results due to the high complexity of this polymorphism, we applied two different techniques (reverse PCR-SSO typing on microbeads arrays based on Luminex technology, and PCR-SSP typing) to type the samples. After validating the resulting frequency distributions through the original statistical method described in our companion article ( Nunes et al. 2011 ), we compared the HLA genetic profile of the sampled Munda to those of other Asiatic populations, among which Dravidian and Indo-European-speakers from India and populations from East and Southeast Asia speaking languages belonging to different linguistic families. We showed that the Munda from Northeast India exhibit a peculiar genetic profile with a reduced level of HLA diversity compared to surrounding Indian populations. They also exhibit less diversity than Southeast Asian populations except at locus DRB1. Several analyses using genetic distances indicate that the Munda are much more closely related to populations from the Indian subcontinent than to Southeast Asian populations speaking languages of the same Austroasiatic linguistic family. On the other hand, they do not share a closer relationship with Dravidians compared with Indo-Europeans, thus arguing against the idea that the Munda share a common and ancient Indian origin with Dravidians. Our results do not favor either a scenario where the Munda would be representative of an ancestral Austroasiatic population giving rise to an eastward Austroasiatic expansion to Southeast Asia. Rather, their peculiar genetic profile is better explained by a decrease in genetic diversity through genetic drift from an ancestral population having a genetic profile similar to present-day Austroasiatic populations from Southeast Asia (thus suggesting a possible southeastern origin), followed by intensive gene flow with neighboring Indian populations. This conclusion is in agreement with archaeological and linguistic information. The history of the Austroasiatic family represents a fascinating example where complex interactions among culturally distinct human populations occurred in the past.     

Human population genetic structure and diversity inferred from polymorphic L1(LINE-1) and Alu insertions

BACKGROUND/AIMS: The L1 retrotransposable element family is the most successful self-replicating genomic parasite of the human genome. L1 elements drive replication of Alu elements, and both have had far-reaching impacts on the human genome. We use L1 and Alu insertion polymorphisms to analyze human population structure. METHODS: We genotyped 75 recent, polymorphic L1 insertions in 317 individuals from 21 populations in sub-Saharan Africa, East Asia, Europe and the Indian subcontinent. This is the first sample of L1 loci large enough to support detailed population genetic inference. We analyzed these data in parallel with a set of 100 polymorphic Alu insertion loci previously genotyped in the same individuals. RESULTS AND CONCLUSION: The data sets yield congruent results that support the recent African origin model of human ancestry. A genetic clustering algorithm detects clusters of individuals corresponding to continental regions. The number of loci sampled is critical: with fewer than 50 typical loci, structure cannot be reliably discerned in these populations. The inclusion of geographically intermediate populations (from India) reduces the distinctness of clustering. Our results indicate that human genetic variation is neither perfectly correlated with geographic distance (purely clinal) nor independent of distance (purely clustered), but a combination of both: stepped clinal.     

Very low genetic variability in the Indian queenless ant Diacamma indicum

Abstract We developed microsatellite markers and combined them with mitochondrial markers to analyse the population genetic structure of the queenless ant Diacamma indicum. This species, lacking winged queens, is likely to have a restricted female dispersal but exhibits various life history traits suggesting higher dispersal abilities than the other Diacamma species. Only 4 of 11 microsatellites were polymorphic and only 1 had more than 4 alleles over 166 individuals originating from 7 populations from the south of India. Only one mitochondrial DNA (mtDNA) haplotype was detected throughout India (including one population in the north) and Sri Lanka. Such a level of polymorphism is particularly low compared with other Diacamma species having much smaller ranges in the south of India. A strong genetic differentiation was observed between populations separated by more than a few kilometres. We also analysed the genetic differentiation between the Indian populations and two populations from the Japanese island of Okinawa, which are morphologically similar and might belong to the same species. The genetic differentiation was high for both markers, suggesting an absence of ongoing gene flow between these populations.     

Normative genetic profiles of RAAS pathway gene polymorphisms in North Indian and South Indian populations

Population-based genetic association studies, popularly known as case-control studies, have continued to be the most preferred method for deciphering the genetic basis of various complex diseases, even in the post-human genome sequencing era. However, interpopulation differences in allele, genotype, and haplotype frequencies and linkage disequilibrium patterns lead to inconsistent results in candidate gene association studies. Therefore, for any meaningful disease association study, knowledge of the normative genetic background of the baseline population is a prerequisite. In addition, such genetic variation data also provide a ready-made menu of allele frequencies and linkage disequilibrium patterns of various polymorphisms in specific candidate genes in a particular population, which is a useful reference for further genetic association studies. Such genetic variation data are lacking for the Indian population, which represents about one-sixth of the world's population. In the present study we have reported the allele, genotype, and haplotype frequencies, Hardy-Weinberg equilibrium status, and linkage disequilibrium patterns of 12 polymorphisms in six candidate genes from the renin-angiotensin-aldosterone system among Indians. Because of their different history of origin, the Indian population is broadly divided into two subpopulations: North Indians (Caucasian Europeans) and South Indians (Dravidians). Considering this well-documented difference in gene pools, we have presented a comparative account of the normative genetic data of North Indian and South Indian populations with at least four individuals of urban and suburban origin from each of the representative states of northern and southern India.     

The Genetics of Bene Israel from India Reveals Both Substantial Jewish and Indian Ancestry

The Bene Israel Jewish community from West India is a unique population whose history before the 18^th century remains largely unknown. Bene Israel members consider themselves as descendants of Jews, yet the identity of Jewish ancestors and their arrival time to India are unknown, with speculations on arrival time varying between the 8th century BCE and the 6th century CE. Here, we characterize the genetic history of Bene Israel by collecting and genotyping 18 Bene Israel individuals. Combining with 486 individuals from 41 other Jewish, Indian and Pakistani populations, and additional individuals from worldwide populations, we conducted comprehensive genome-wide analyses based on F_ST, principal component analysis, ADMIXTURE, identity-by-descent sharing, admixture linkage disequilibrium decay, haplotype sharing and allele sharing autocorrelation decay, as well as contrasted patterns between the X chromosome and the autosomes. The genetics of Bene Israel individuals resemble local Indian populations, while at the same time constituting a clearly separated and unique population in India. They are unique among Indian and Pakistani populations we analyzed in sharing considerable genetic ancestry with other Jewish populations. Putting together the results from all analyses point to Bene Israel being an admixed population with both Jewish and Indian ancestry, with the genetic contribution of each of these ancestral populations being substantial. The admixture took place in the last millennium, about 19–33 generations ago. It involved Middle-Eastern Jews and was sex-biased, with more male Jewish and local female contribution. It was followed by a population bottleneck and high endogamy, which can lead to increased prevalence of recessive diseases in this population. This study provides an example of how genetic analysis advances our knowledge of human history in cases where other disciplines lack the relevant data to do so.     

Why sampling scheme matters: the effect of sampling scheme on landscape genetic results

There has been a recent trend in genetic studies of wild populations where researchers have changed their sampling schemes from sampling pre-defined populations to sampling individuals uniformly across landscapes. This reflects the fact that many species under study are continuously distributed rather than clumped into obvious “populations”. Once individual samples are collected, many landscape genetic studies use clustering algorithms and multilocus genetic data to group samples into subpopulations. After clusters are derived, landscape features that may be acting as barriers are examined and described. In theory, if populations were evenly sampled, this course of action should reliably identify population structure. However, genetic gradients and irregularly collected samples may impact the composition and location of clusters. We built genetic models where individual genotypes were either randomly distributed across a landscape or contained gradients created by neighbor mating for multiple generations. We investigated the influence of six different sampling protocols on population clustering using program STRUCTURE, the most commonly used model-based clustering method for multilocus genotype data. For models where individuals (and their alleles) were randomly distributed across a landscape, STRUCTURE correctly predicted that only one population was being sampled. However, when gradients created by neighbor mating existed, STRUCTURE detected multiple, but different numbers of clusters, depending on sampling protocols. We recommend testing for fine scale autocorrelation patterns prior to sample clustering, as the scale of the autocorrelation appears to influence the results. Further, we recommend that researchers pay attention to the impacts that sampling may have on subsequent population and landscape genetic results. The U.S. Government's right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

The origin of Indian Star tortoises (Geochelone elegans) based on nuclear and mitochondrial DNA analysis: A story of rescue and repatriation

The Indian Star tortoise (Geochelone elegans) belongs to the family Testunidae and is distributed in southwest India and Sri Lanka. In addition to facing loss of its natural habitat, the species is also illegally traded as food and as an exotic pet internationally. Here we report DNA-based analyses for identification and repatriation of these tortoises into their natural habitat. We have attempted to establish the geographical origin of these tortoises rescued from smugglers, by comparing the microsatellite and mitochondrial markers of rescued animals with animals of known provenance. Star tortoises exhibited strong genetic structure in India. The populations from western India were genetically distinct at microsatellite and mitochondrial loci from southern populations. The rescued individuals had similar multilocus genotypes and mitochondrial DNA haplotypes as the reference individuals from south India. However, the precise geographic origin of many of the rescued samples remains unresolved, because we could not assign them to southern populations and the Neighbor-Joining cluster analysis indicated that some of rescued tortoises formed distinct clusters. These data strongly suggest that the rescued group of tortoises is composed of a mix of individuals from differentiated source populations that are probably located in southern India and possibly Sri Lanka. Our study provides valuable information based on molecular markers for the assessment of genetic diversity in Indian Star tortoises.     

A population genetic study of six VNTR loci in three ethnically defined populations.

To investigate the population genetic characteristics of VNTR polymorphisms in human populations, we have studied the allele frequency distribution of six VNTR loci (D1S57, RB1, D1S77, D1S61, alpha-globin 5'HVR, D1S76) in three well-defined populations (Kachari of Northeast India; Dogrib Indian of Canada; and New Guinea Highlander of Papua New Guinea). Even though the number of alleles sampled is limited, 48 to 92 alleles per locus per population, significant variation is noticed in the number of alleles per locus for all the populations. Using alternate summary measures, we have observed that genotype distributions at the six VNTR loci apparently conform to their respective Hardy-Weinberg predictions. Multilocus genotype profiles of the individuals in each of the three populations suggest that the VNTR alleles are independently segregating with the exception of the two linked loci D1S76 and D1S77. Lack of fit of all VNTR loci to one particular model of mutational change, either the Infinite Allele Model or the Stepwise Mutation Model, suggests more than one mechanism for production of new VNTR alleles. This study also indicates that increased heterozygosity at VNTR loci in comparison to protein and blood group loci may lead to more accurate estimates of genetic distance.     

Multiple origins of the mtDNA 9-bp deletion in populations of South India

The origins and genetic affinities of the more than 500 tribal populations living in South Asia are widely disputed. This may reflect differential contributions that continental populations have made to tribal groups in South Asia. We assayed for the presence of the intergenic COII/tRNALys 9-bp deletion in human mtDNA in 646 individuals from 12 caste and 14 tribal populations of South India and compared them to individuals from Africa, Europe, and Asia. The 9-bp deletion is observed in four South Indian tribal populations, the Irula, Yanadi, Siddi, and Maria Gond, and in the Nicobarese. Length polymorphisms of the 9-bp motif are present in the Santal, Khonda Dora, and Jalari, all of whom live in a circumscribed region on the eastern Indian coast. Phylogenetic analyses of mtDNA control region sequence from individuals with the 9-bp deletion indicate that it has arisen independently in some Indian tribal populations. Other 9-bp deletion haplotypes are likely to be of Asian and African origin, implying multiple origins of the 9-bp deletion in South India. These results demonstrate varying genetic affinities of different South Indian tribes to continental populations and underscore the complex histories of the tribal populations living in South Asia. Am J Phys Anthropol 109:147–158, 1999. © 1999 Wiley-Liss, Inc.     

Integrating Linguistics,Social Structure,and Geography to Model Genetic Diversity within India

India represents an intricate tapestry of population substructure shaped by geography, language, culture, and social stratification. Although geography closely correlates with genetic structure in other parts of the world, the strict endogamy imposed by the Indian caste system and the large number of spoken languages add further levels of complexity to understand Indian population structure. To date, no study has attempted to model and evaluate how these factors have interacted to shape the patterns of genetic diversity within India. We merged all publicly available data from the Indian subcontinent into a data set of 891 individuals from 90 well-defined groups. Bringing together geography, genetics, and demographic factors, we developed Correlation Optimization of Genetics and Geodemographics to build a model that explains the observed population genetic substructure. We show that shared language along with social structure have been the most powerful forces in creating paths of gene flow in the subcontinent. Furthermore, we discover the ethnic groups that best capture the diverse genetic substructure using a ridge leverage score statistic. Integrating data from India with a data set of additional 1,323 individuals from 50 Eurasian populations, we find that Indo-European and Dravidian speakers of India show shared genetic drift with Europeans, whereas the Tibeto-Burman speaking tribal groups have maximum shared genetic drift with East Asians.     

Why the Indian Subcontinent Holds the Key to Global Tiger Recovery

With only ~3,000 wild individuals surviving restricted to just 7% of their historical range, tigers are now a globally threatened species. Therefore, conservation efforts must prioritize regions that harbor more tigers, as well try to capture most of the remaining genetic variation and habitat diversity. Only such prioritization based on demographic, genetic, and ecological considerations can ensure species recovery and retention of evolutionary flexibility in the face of ongoing global changes. Although scientific understanding of ecological and demographic aspects of extant wild tiger populations has improved recently, little is known about their genetic composition and variability. We sampled 73 individual tigers from 28 reserves spread across a diversity of habitats in the Indian subcontinent to obtain 1,263 bp of mitochondrial DNA and 10 microsatellite loci. Our analyses reveals that Indian tigers retain more than half of the extant genetic diversity in the species. Coalescent simulations attribute this high genetic diversity to a historically large population size of about 58,200 tigers for peninsular India south of the Gangetic plains. Furthermore, our analyses indicate a precipitous, possibly human-induced population crash ~200 years ago in India, which is in concordance with historical records. Our results suggest that only 1.7% (with an upper limit of 13% and a lower limit of 0.2%) of tiger numbers in historical times remain now. In the global conservation context our results suggest that, based on genetic, demographic, and ecological considerations, the Indian subcontinent holds the key to global survival and recovery of wild tigers.     

The Impact of Divergence Time on the Nature of Population Structure: An Example from Iceland

The Icelandic population has been sampled in many disease association studies, providing a strong motivation to understand the structure of this population and its ramifications for disease gene mapping. Previous work using 40 microsatellites showed that the Icelandic population is relatively homogeneous, but exhibits subtle population structure that can bias disease association statistics. Here, we show that regional geographic ancestries of individuals from Iceland can be distinguished using 292,289 autosomal single-nucleotide polymorphisms (SNPs). We further show that subpopulation differences are due to genetic drift since the settlement of Iceland 1100 years ago, and not to varying contributions from different ancestral populations. A consequence of the recent origin of Icelandic population structure is that allele frequency differences follow a null distribution devoid of outliers, so that the risk of false positive associations due to stratification is minimal. Our results highlight an important distinction between population differences attributable to recent drift and those arising from more ancient divergence, which has implications both for association studies and for efforts to detect natural selection using population differentiation.     

The 9-bp deletion between the mitochondrial lysine tRNA and COII genes in tribal populations of India

A 9-base-pair (bp) deletion located between the lysine tRNA (MTTK) and COII (MTCOX*2) genes in the human mitochondrial genome is a valuable marker for tracing population relationships. Previous research has shown that the 9-bp deletion is associated with two major clusters of control region sequences; one occurs in sub-Saharan Africa, while the other is associated with Asian populations and populations of Asian origin. We surveyed 898 individuals from 16 tribal populations in India and found 6 individuals with the 9-bp deletion. Sequences of the first hypervariable segment (HV1) of the mtDNA control region from these 9-bp deletion-bearing mtDNAs were compared to those previously reported from Asian and African populations. Phylogenetic analysis indicates three distinct clusters of tribal Indian 9-bp deletion mtDNA types. One cluster, found in northeast India, includes southeast Asian and Indonesian mtDNA types. The remaining two clusters appear to have unique origins in southern India. These data provide further evidence of past migrations from Asia into the northeast corner of the Indian subcontinent.     

Population and genomic lessons from genetic analysis of two Indian populations

Indian demographic history includes special features such as founder effects, interpopulation segregation, complex social structure with a caste system and elevated frequency of consanguineous marriages. It also presents a higher frequency for some rare mendelian disorders and in the last two decades increased prevalence of some complex disorders. Despite the fact that India represents about one-sixth of the human population, deep genetic studies from this terrain have been scarce. In this study, we analyzed high-density genotyping and whole-exome sequencing data of a North and a South Indian population. Indian populations show higher differentiation levels than those reported between populations of other continents. In this work, we have analyzed its consequences, by specifically assessing the transferability of genetic markers from or to Indian populations. We show that there is limited genetic marker portability from available genetic resources such as HapMap or the 1,000 Genomes Project to Indian populations, which also present an excess of private rare variants. Conversely, tagSNPs show a high level of portability between the two Indian populations, in contrast to the common belief that North and South Indian populations are genetically very different. By estimating kinship from mates and consanguinity in our data from trios, we also describe different patterns of assortative mating and inbreeding in the two populations, in agreement with distinct mating preferences and social structures. In addition, this analysis has allowed us to describe genomic regions under recent adaptive selection, indicating differential adaptive histories for North and South Indian populations. Our findings highlight the importance of considering demography for design and analysis of genetic studies, as well as the need for extending human genetic variation catalogs to new populations and particularly to those with particular demographic histories.     

Shared and unique components of human population structure and genome-wide signals of positive selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.     

Identifying the source of species invasions: sampling intensity vs. genetic diversity

Population geneticists and community ecologists have long recognized the importance of sampling design for uncovering patterns of diversity within and among populations and in communities. Invasion ecologists increasingly have utilized phylogeographical patterns of mitochondrial or chloroplast DNA sequence variation to link introduced populations with putative source populations. However, many studies have ignored lessons from population genetics and community ecology and are vulnerable to sampling errors owing to insufficient field collections. A review of published invasion studies that utilized mitochondrial or chloroplast DNA markers reveals that insufficient sampling could strongly influence results and interpretations. Sixty per cent of studies sampled an average of less than six individuals per source population, vs. only 45% for introduced populations. Typically, far fewer introduced than source populations were surveyed, although they were sampled more intensively. Simulations based on published data forming a comprehensive mtDNA haplotype data set highlight and quantify the impact of the number of individuals surveyed per source population and number of putative source populations surveyed for accurate assignment of introduced individuals. Errors associated with sampling a low number of individuals are most acute when rare source haplotypes are dominant or fixed in the introduced population. Accuracy of assignment of introduced individuals is also directly related to the number of source populations surveyed and to the degree of genetic differentiation among them ( F _ST). Incorrect interpretations resulting from sampling errors can be avoided if sampling design is considered before field collections are made.     

Genetic landscape of the people of India: a canvas for disease gene exploration

Analyses of frequency profiles of markers on disease or drug-response related genes in diverse populations are important for the dissection of common diseases. We report the results of analyses of data on 405 SNPs from 75 such genes and a 5.2 Mb chromosome, 22 genomic region in 1871 individuals from diverse 55 endogamous Indian populations. These include 32 large (> 10 million individuals) and 23 isolated populations, representing a large fraction of the people of India. We observe high levels of genetic divergence between groups of populations that cluster largely on the basis of ethnicity and language. Indian populations not only overlap with the diversity of HapMap populations, but also contain population groups that are genetically distinct. These data and results are useful for addressing stratification and study design issues in complex traits especially for heterogeneous populations. For complete authors’ list of the Indian Genome Variation Database Consortium please see Appendix. For correspondence. E-mail: skb@igib.res.in; Functional Genomics Unit, Institute of Genomics and Integrative Biology (CSIR), Mall Road, New Delhi 110 007, India.     

Origins and Divergence of the Roma (Gypsies)

The identification of a growing number of novel Mendelian disorders and private mutations in the Roma (Gypsies) points to their unique genetic heritage. Linguistic evidence suggests that they are of diverse Indian origins. Their social structure within Europe resembles that of the jatis of India, where the endogamous group, often defined by profession, is the primary unit. Genetic studies have reported dramatic differences in the frequencies of mutations and neutral polymorphisms in different Romani populations. However, these studies have not resolved ambiguities regarding the origins and relatedness of Romani populations. In this study, we examine the genetic structure of 14 well-defined Romani populations. Y-chromosome and mtDNA markers of different mutability were analyzed in a total of 275 individuals. Asian Y-chromosome haplogroup VI-68, defined by a mutation at the M82 locus, was present in all 14 populations and accounted for 44.8% of Romani Y chromosomes. Asian mtDNA-haplogroup M was also identified in all Romani populations and accounted for 26.5% of female lineages in the sample. Limited diversity within these two haplogroups, measured by the variation at eight short-tandem-repeat loci for the Y chromosome, and sequencing of the HVS1 for the mtDNA are consistent with a small group of founders splitting from a single ethnic population in the Indian subcontinent. Principal-components analysis and analysis of molecular variance indicate that genetic structure in extant endogamous Romani populations has been shaped by genetic drift and differential admixture and correlates with the migrational history of the Roma in Europe. By contrast, social organization and professional group divisions appear to be the product of a more recent restitution of the caste system of India.     

Genetic Population Structure Analysis in New Hampshire Reveals Eastern European Ancestry

Genetic structure due to ancestry has been well documented among many divergent human populations. However, the ability to associate ancestry with genetic substructure without using supervised clustering has not been explored in more presumably homogeneous and admixed US populations. The goal of this study was to determine if genetic structure could be detected in a United States population from a single state where the individuals have mixed European ancestry. Using Bayesian clustering with a set of 960 single nucleotide polymorphisms (SNPs) we found evidence of population stratification in 864 individuals from New Hampshire that can be used to differentiate the population into six distinct genetic subgroups. We then correlated self-reported ancestry of the individuals with the Bayesian clustering results. Finnish and Russian/Polish/Lithuanian ancestries were most notably found to be associated with genetic substructure. The ancestral results were further explained and substantiated using New Hampshire census data from 1870 to 1930 when the largest waves of European immigrants came to the area. We also discerned distinct patterns of linkage disequilibrium (LD) between the genetic groups in the growth hormone receptor gene (GHR). To our knowledge, this is the first time such an investigation has uncovered a strong link between genetic structure and ancestry in what would otherwise be considered a homogenous US population.