Complex Demographic History Research Articles

Maternal Genetic Heritage of Southeastern Europe Reveals a New Croatian Isolate and a Novel, Local Sub-Branching in the X2 Haplogroup

High mtDNA variation in Southeastern Europe (SEE) is a reflection of the turbulent and complex demographic history of this area, influenced by gene flow from various parts of Eurasia and a long history of intermixing. Our results of 1035 samples (488 from Croatia, 239 from Bosnia and 130 from Herzegovina, reported earlier, and 97 Slovenians and 81 individuals from Žumberak, reported here for the first time) show that the SEE maternal genetic diversity fits within a broader European maternal genetic landscape. The study also shows that the population of Žumberak, located in the continental part of Croatia, developed some unique mtDNA haplotypes and elevated haplogroup frequencies due to distinctive demographic history and can be considered a moderate genetic isolate. We also report seven samples from the Bosnian population and one Herzegovinian sample designated as X2* individuals that could not be assigned to any of its sublineages (X2a'o) according to the existing X2 phylogeny. In an attempt to clarify the phylogeny of our X2 samples, their mitochondrial DNA has been completely sequenced. We suppose that these lineages are signs of local microdifferentiation processes that occurred in the recent demographic past in this area and could possibly be marked as SEE-specific X2 sublineages.

Genomic Analyses of Three Malaria Vectors Reveals Extensive Shared Polymorphism but Contrasting Population Histories

Anopheles gambiae s.l. are important malaria vectors, but little is known about their genomic variation in the wild. Here, we present inter- and intraspecies analysis of genome-wide RADseq data, in three Anopheles gambiae s.l. species collected from East Africa. The mosquitoes fall into three genotypic clusters representing described species (A. gambiae, A. arabiensis, and A. merus) with no evidence of cryptic breeding units. Anopheles merus is the most divergent of the three species, supporting a recent new phylogeny based on chromosomal inversions. Even though the species clusters are well separated, there is extensive shared polymorphism, particularly between A. gambiae and A. arabiensis. Divergence between A. gambiae and A. arabiensis does not vary across the autosomes but is higher in X-linked inversions than elsewhere on X or on the autosomes, consistent with the suggestion that this inversion (or a gene within it) is important in reproductive isolation between the species. The 2La/2L+a inversion shows no more evidence of introgression between A. gambiae and A. arabiensis than the rest of the autosomes. Population differentiation within A. gambiae and A. arabiensis is weak over approximately 190–270 km, implying no strong barriers to dispersal. Analysis of Tajima’s D and the allele frequency spectrum is consistent with modest population increases in A. arabiensis and A. merus, but a more complex demographic history of expansion followed by contraction in A. gambiae. Although they are less than 200 km apart, the two A. gambiae populations show evidence of different demographic histories.

Population genomic analysis reveals a rich speciation and demographic history of orang-utans (Pongo pygmaeus and Pongo abelii).

To gain insights into evolutionary forces that have shaped the history of Bornean and Sumatran populations of orang-utans, we compare patterns of variation across more than 11 million single nucleotide polymorphisms found by previous mitochondrial and autosomal genome sequencing of 10 wild-caught orang-utans. Our analysis of the mitochondrial data yields a far more ancient split time between the two populations (∼3.4 million years ago) than estimates based on autosomal data (0.4 million years ago), suggesting a complex speciation process with moderate levels of primarily male migration. We find that the distribution of selection coefficients consistent with the observed frequency spectrum of autosomal non-synonymous polymorphisms in orang-utans is similar to the distribution in humans. Our analysis indicates that 35% of genes have evolved under detectable negative selection. Overall, our findings suggest that purifying natural selection, genetic drift, and a complex demographic history are the dominant drivers of genome evolution for the two orang-utan populations.

High Genetic Diversity and Distinctiveness of Rear-Edge Climate Relicts Maintained by Ancient Tetraploidisation for Alnus glutinosa

Populations located at the rear-edge of a species’ distribution may have disproportionate ecological and evolutionary importance for biodiversity conservation in a changing global environment. Yet genetic studies of such populations remain rare. This study investigates the evolutionary history of North-African low latitude marginal populations of Alnus glutinosa Gaertn., a European tree species that plays a significant ecological role as a keystone of riparian ecosystems. We genotyped 551 adults from 19 populations located across North Africa at 12 microsatellite loci and applied a coalescent-based simulation approach to reconstruct the demographic and evolutionary history of these populations. Surprisingly, Moroccan trees were tetraploids demonstrating a strong distinctiveness of these populations within a species otherwise known as diploid. Best-fitting models of demographic reconstruction revealed the relict nature of Moroccan populations that were found to have withstood past climate change events and to be much older than Algerian and Tunisian populations. This study highlights the complex demographic history that can be encountered in rear-edge distribution margins that here consist of both old stable climate relict and more recent populations, distinctively diverse genetically both quantitatively and qualitatively. We emphasize the high evolutionary and conservation value of marginal rear-edge populations of a keystone riparian species in the context of on-going climate change in the Mediterranean region.

Insights from population genetics for range limits of a widely distributed native plant

Species' range limits are determined by a combination of history, ecology, and genetics. While much of our understanding of range limits comes from ecological studies, molecular population genetic methods can provide insight into the contribution of historical range shifts and gene flow to current range limits. We sequenced nine nuclear loci in 68 individuals from four populations (two within the range interior, one at the western range edge, and one at northern range edge) of the native annual legume Chamaecrista fasciculata in the Upper Midwest of North America. Molecular diversity was greatest in an interior population and significantly reduced at both the western and northern range edges. We found no molecular genetic evidence for historical demographic expansion or contraction in any of the populations. Coalescent simulations indicate that the absence of a genetic signal for expansion or contraction could be due to inadequate statistical power or to a bottleneck that was weak and old. Populations were significantly differentiated and clustered into three demes, with the fourth population containing individuals that were assigned to the two nearest populations. We found evidence that recent migration between the northern edge and the most proximate interior population has been slight, within the distribution of values that theoretical models predict would alleviate genetic load without impeding adaptation. When coupled with results from transplant experiments, these results suggest that ecological-genetic factors other than gene flow currently limit range expansion of C. fasciculata. Further, they highlight challenges in using molecular data to make inferences about species' distributions that have complex demographic histories.

Cryptic Genetic Diversity within the Anopheles nili group of Malaria Vectors in the Equatorial Forest Area of Cameroon (Central Africa)

BackgroundThe Anopheles nili group of mosquitoes includes important vectors of human malaria in equatorial forest and humid savannah regions of sub-Saharan Africa. However, it remains largely understudied, and data on its populations’ bionomics and genetic structure are crucially lacking. Here, we used a combination of nuclear (i.e. microsatellite and ribosomal DNA) and mitochondrial DNA markers to explore and compare the level of genetic polymorphism and divergence among populations and species of the group in the savannah and forested areas of Cameroon, Central Africa.Principal FindingsAll the markers provided support for the current classification within the An. nili group. However, they revealed high genetic heterogeneity within An. nili s.s. in deep equatorial forest environment. Nuclear markers showed the species to be composed of five highly divergent genetic lineages that differed by 1.8 to 12.9% of their Internal Transcribed Spacer 2 (ITS2) sequences, implying approximate divergence time of 0.82 to 5.86 million years. However, mitochondrial data only detected three major subdivisions, suggesting different evolutionary histories of the markers.Conclusions/SignificanceThis study enlightened additional cryptic genetic diversity within An. nili s.s. in the deep equatorial forest environment of South Cameroon, reflecting a complex demographic history for this major vector of malaria in this environment. These preliminary results should be complemented by further studies which will shed light on the distribution, epidemiological importance and evolutionary history of this species group in the African rainforest, providing opportunities for in-depth comparative studies of local adaptation and speciation in major African malaria vectors.

Mitochondrial DNA diversity in two ethnic groups in Southeastern Kenya: Perspectives from the northeastern periphery of the bantu expansion

The Bantu languages are widely distributed throughout sub-Saharan Africa. Genetic research supports linguists and historians who argue that migration played an important role in the spread of this language family, but the genetic data also indicates a more complex process involving substantial gene flow with resident populations. In order to understand the Bantu expansion process in east Africa, mtDNA hypervariable region I variation in 352 individuals from the Taita and Mijikenda ethnic groups was analyzed, and we evaluated the interactions that took place between the Bantu- and non-Bantu-speaking populations in east Africa. The Taita and Mijikenda are Bantu-speaking agropastoralists from southeastern Kenya, at least some of whose ancestors probably migrated into the area as part of Bantu migrations that began around 3,000 BCE. Our analyses indicate that they show some distinctive differences that reflect their unique cultural histories. The Taita are genetically more diverse than the Mijikenda with larger estimates of genetic diversity. The Taita cluster with other east African groups, having high frequencies of haplogroups from that region, while the Mijikenda have high frequencies of central African haplogroups and cluster more closely with central African Bantu-speaking groups. The non-Bantu speakers who lived in southeastern Kenya before Bantu speaking groups arrived were at least partially incorporated into what are now Bantu-speaking Taita groups. In contrast, gene flow from non-Bantu speakers into the Mijikenda was more limited. These results suggest a more complex demographic history where the nature of Bantu and non-Bantu interactions varied throughout the area.

Polymorphism at genes involved in salt tolerance in Arabidopsis thaliana (Brassicaceae)

Genes involved in relevant functions for environmental adaptation can be considered primary candidates for their variation having been shaped by natural selection. Detecting recent selective events through their footprint on nucleotide variation constitutes a challenging task in species with a complex demographic history such as Arabidopsis thaliana. We have surveyed nucleotide variation in this species at nine genes involved in salt tolerance. The available genomewide information for this species has allowed us to contrast the levels and patterns of variation detected at the candidate genes with empirical distributions obtained from noncandidate regions. We sequenced nine genes involved in salt tolerance (~32 kb) in 20 ecotypes of A. thaliana and analyzed polymorphism and divergence at the individual gene and multilocus levels. Variation at the nine genes studied was characterized by a generalized skew toward polymorphisms with low-frequency variants. Except for genes RCD1 and NHX8, this pattern was similar to that generally detected in the A. thaliana genome and could thus be primarily explained by the species demographic history. The more extreme deviation at the NHX8 gene and its excess of polymorphism relative to divergence points to the recent action of selection on this gene. The analysis of nucleotide polymorphism and divergence at nine genes involved in salt tolerance provided little evidence for the recent action of positive selection. Only the signals detected at NHX8 from both polymorphism and divergence were suggestive of the putative contribution of this gene to local adaptation.

Genetic characterization of northeastern Italian population isolates in the context of broader European genetic diversity

Population genetic studies on European populations have highlighted Italy as one of genetically most diverse regions. This is possibly due to the country's complex demographic history and large variability in terrain throughout the territory. This is the reason why Italy is enriched for population isolates, Sardinia being the best-known example. As the population isolates have a great potential in disease-causing genetic variants identification, we aimed to genetically characterize a region from northeastern Italy, which is known for isolated communities. Total of 1310 samples, collected from six geographically isolated villages, were genotyped at >145 000 single-nucleotide polymorphism positions. Newly genotyped data were analyzed jointly with the available genome-wide data sets of individuals of European descent, including several population isolates. Despite the linguistic differences and geographical isolation the village populations still show the greatest genetic similarity to other Italian samples. The genetic isolation and small effective population size of the village populations is manifested by higher levels of genomic homozygosity and elevated linkage disequilibrium. These estimates become even more striking when the detected substructure is taken into account. The observed level of genetic isolation in Friuli-Venezia Giulia region is more extreme according to several measures of isolation compared with Sardinians, French Basques and northern Finns, thus proving the status of an isolate.

Inferring the History of Population Size Change from Genome-Wide SNP Data

Dense, genome-wide single-nucleotide polymorphism (SNP) data can be used to reconstruct the demographic history of human populations. However, demographic inferences from such data are complicated by recombination and ascertainment bias. We introduce two new statistics, allele frequency-identity by descent (AF-IBD) and allele frequency-identity by state (AF-IBS), that make use of linkage disequilibrium information and show defined relationships to the time of coalescence. These statistics, when conditioned on the derived allele frequency, are able to infer complex population size changes. Moreover, the AF-IBS statistic, which is based on genome-wide SNP data, is robust to varying ascertainment conditions. We constructed an efficient approximate Bayesian computation (ABC) pipeline based on AF-IBD and AF-IBS that can accurately estimate demographic parameters, even for fairly complex models. Finally, we applied this ABC approach to genome-wide SNP data and inferred the demographic histories of two human populations, Yoruba and French. Our results suggest a rather stable ancestral population size with a mild recent expansion for Yoruba, whereas the French seemingly experienced a long-lasting severe bottleneck followed by a drastic population growth. This approach should prove useful for new insights into populations, especially those with complex demographic histories.

Comparative Nucleotide Diversity Across North American and European Populus Species

Nucleotide polymorphisms in two North American balsam poplars (Populus trichocarpa Torr. & Gray and P. balsamifera L.; section Tacamahaca), and one Eurasian aspen (P. tremula L.; section Populus) were compared using nine loci involved in defense, stress response, photoperiodism, freezing tolerance, and housekeeping. Nucleotide diversity varied among species and was highest for P. tremula (θ(w) = 0.005, π(T) = 0.007) as compared to P. balsamifera (θ(w) = 0.004, π(T) = 0.005) or P. trichocarpa (θ(w) = 0.002, π(T) = 0.003). Across species, the defense and the stress response loci accounted for the majority of the observed level of nucleotide diversity. In general, the studied loci did not deviate from neutral expectation either at the individual locus (non-significant normalized Fay and Wu's H) or at the multi-locus level (non-significant HKA test). Using molecular clock analysis, section Tacamahaca probably shared a common ancestor with section Populus approximately 4.5 million year ago. Divergence between the two closely related balsam poplars was about 0.8 million years ago, a pattern consistent with an isolation-with-migration (IM) model. As expected, P. tremula showed a five-fold higher substitution rate (2 × 10(-8) substitution/site/year) compared to the North American species (0.4 × 10(-8) substitution/site/year), probably reflecting its complex demographic history. Linkage disequilibrium (LD) varied among species with a more rapid decay in the North American species (<400 bp) in comparison to P. tremula (≫400 bp). The similarities in nucleotide diversity pattern and LD decay of the two balsam poplar species likely reflects the recent time of their divergence.

The complete mitochondrial genome sequence of the filarial nematode Wuchereria bancrofti from three geographic isolates provides evidence of complex demographic history

Mitochondrial (mt) genome sequences have enabled comparison of population genetics and evolution for numerous free-living and parasitic nematodes. Here we define the complete mt genome of Wuchereria bancrofti through analysis of isolates from Papua New Guinea, India and West Africa. Sequences were assembled for each isolate and annotated with reference to the mt genome sequence for Brugia malayi. The length of the W. bancrofti mt genome is approximately 13,637 nucleotides, contains 2 ribosomal RNAs (rrns), 22 transfer RNAs (trns), 12 protein-coding genes, and is characterized by a 74.6% AT content. The W. bancrofti mt gene order is identical to that reported for Onchocerca volvulus, Dirofilaria immitis, Setaria digitata and B. malayi. In addition to using translational start codons identified previously in the mt protein-coding genes of other filarial nematodes, W. bancrofti appears to be unique in using TGT as a translational start codon. Similarly, use of incomplete stop codons in mt protein-coding genes appears to be more common in W. bancrofti than in other human filarial parasites. The complete mt genome sequence reported here provides new genetic markers for investigating phylogenetic and geographic relationships between isolates, and assessing population diversity within endemic regions. The sequence polymorphism enables new strategies to monitor the progress of public health interventions to control and eliminate this important human parasite. We illustrate the utility of this sequence and single nucleotide polymorphisms by inferring the divergence times between the three W. bancrofti isolates, suggesting predictions into their origin and migration.

Complex Population Dynamics and the Coalescent Under Neutrality

Estimates of the coalescent effective population size Ne can be poorly correlated with the true population size. The relationship between Ne and the population size is sensitive to the way in which birth and death rates vary over time. The problem of inference is exacerbated when the mechanisms underlying population dynamics are complex and depend on many parameters. In instances where nonparametric estimators of Ne such as the skyline struggle to reproduce the correct demographic history, model-based estimators that can draw on prior information about population size and growth rates may be more efficient. A coalescent model is developed for a large class of populations such that the demographic history is described by a deterministic nonlinear dynamical system of arbitrary dimension. This class of demographic model differs from those typically used in population genetics. Birth and death rates are not fixed, and no assumptions are made regarding the fraction of the population sampled. Furthermore, the population may be structured in such a way that gene copies reproduce both within and across demes. For this large class of models, it is shown how to derive the rate of coalescence, as well as the likelihood of a gene genealogy with heterochronous sampling and labeled taxa, and how to simulate a coalescent tree conditional on a complex demographic history. This theoretical framework encapsulates many of the models used by ecologists and epidemiologists and should facilitate the integration of population genetics with the study of mathematical population dynamics.

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.

Multi-Locus Analysis Reveals A Different Pattern of Genetic Diversity for Mitochondrial and Nuclear DNA between Wild and Domestic Pigs in East Asia

BackgroundA major reduction of genetic diversity in mtDNA occurred during the domestication of East Asian pigs. However, the extent to which genetic diversity has been lost in the nuclear genome is uncertain. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationships and demographic history of domestic pigs and their ancestors, wild boars, we investigated 14 nuclear markers (including 8 functional genes, 2 pseudogenes and 4 intergenic regions) from 11 different chromosomes in East Asia-wide samples and pooled them with previously obtained mtDNA data for a combined analysis.Principal FindingsThe results indicated that domestic pigs and wild boars possess comparable levels of nucleotide diversity across the nuclear genome, which is inconsistent with patterns that have been found in mitochondrial genome.ConclusionsThis incongruence between the mtDNA and nuclear genomes is suggestive of a large-scale backcross between male wild boars and female domestic pigs in East Asia. Our data reveal the impacts of founder effects and backcross on the pig genome and help us better understand the complex demographic histories of East Asian pigs, which will be useful for future work on artificial selection.

Historical and contemporary population genetic connectivity of the European short‐snouted seahorse Hippocampus hippocampus and implications for management

This first genetic study of Hippocampus hippocampus covers the species' entire geographic range and employs two mtDNA markers (control region and cytochrome b) to establish patterns of population structuring. A total of 255 specimens from 21 locations were used to obtain 89 concatenated haplotypes. The common haplotype was present in all but one population, however, most haplotypes were unique. The haplotype network had a star-like construction, suggesting expansion from a bottleneck event. F(ST) and AMOVA revealed population subdivision into three geographic regions (English Channel + Bay of Biscay, Mediterranean Sea + Atlantic Ocean Iberian coast + Macaronesian Islands, and West Africa) with barriers to gene flow indentified at Cape Finisterre and the Cape Verde frontal zone. Neutrality tests and nested clade analysis suggest a complex demographic history, with both historic events and contemporary processes shaping patterns of genetic differentiation. The genetic population subdivision detected in this study indicates that H. hippocampus should be managed as three separate units. This is especially pertinent as H. hippocampus populations within the West African region are the only ones known to be specifically targeted for exploitation.

Comment on "Minimum Effective Volume of Local Anesthetic for Shoulder Analgesia by Ultrasound-Guided Block at Root C7 With Assessment of Pulmonary Function"

The observation that even a tiny sample of genome sequences from a natural population contains a plethora of information about the history of the population has enticed researchers to use these data to fit complex demographic histories and make detailed inference about the changes a population has experienced through time. Unfortunately, the standard assumptions required to make these inferences are often violated by natural populations in such ways as to produce specious results. This paper examines two phenomena of particular concern: when a sample is drawn from a single sub-population of a larger meta-population these models infer a spurious recent population decline, and when a genome contains loci under weak or recessive purifying selection these models infer a spurious recent population expansion.

Patterns of Genetic Variation Within and Between Gibbon Species

Gibbons are small, arboreal, highly endangered apes that are understudied compared with other hominoids. At present, there are four recognized genera and approximately 17 species, all likely to have diverged from each other within the last 5–6 My. Although the gibbon phylogeny has been investigated using various approaches (i.e., vocalization, morphology, mitochondrial DNA, karyotype, etc.), the precise taxonomic relationships are still highly debated. Here, we present the first survey of nuclear sequence variation within and between gibbon species with the goal of estimating basic population genetic parameters. We gathered ∼60 kb of sequence data from a panel of 19 gibbons representing nine species and all four genera. We observe high levels of nucleotide diversity within species, indicative of large historical population sizes. In addition, we find low levels of genetic differentiation between species within a genus comparable to what has been estimated for human populations. This is likely due to ongoing or episodic gene flow between species, and we estimate a migration rate between Nomascus leucogenys and N. gabriellae of roughly one migrant every two generations. Together, our findings suggest that gibbons have had a complex demographic history involving hybridization or mixing between diverged populations.

The Genetic Structure of Domestic Rabbits

Understanding the genetic structure of domestic species provides a window into the process of domestication and motivates the design of studies aimed at making links between genotype and phenotype. Rabbits exhibit exceptional phenotypic diversity, are of great commercial value, and serve as important animal models in biomedical research. Here, we provide the first comprehensive survey of nucleotide polymorphism and linkage disequilibrium (LD) within and among rabbit breeds. We resequenced 16 genomic regions in population samples of both wild and domestic rabbits and additional 35 fragments in 150 rabbits representing six commonly used breeds. Patterns of genetic variation suggest a single origin of domestication in wild populations from France, supporting historical records that place rabbit domestication in French monasteries. Levels of nucleotide diversity both within and among breeds were ~0.2%, but only 60% of the diversity present in wild populations from France was captured by domestic rabbits. Despite the recent origin of most breeds, levels of population differentiation were high (F(ST) = 17.9%), but the majority of polymorphisms were shared and thus transferable among breeds. Coalescent simulations suggest that domestication began with a small founding population of less than 1,200 individuals. Taking into account the complex demographic history of domestication with two successive bottlenecks, two loci showed deviations that were consistent with artificial selection, including GPC4, which is known to be associated with growth rates in humans. Levels of diversity were not significantly different between autosomal and X-linked loci, providing no evidence for differential contributions of males and females to the domesticated gene pool. The structure of LD differed substantially within and among breeds. Within breeds, LD extends over large genomic distances. Markers separated by 400 kb typically showed r(2) higher than 0.2, and some LD extended up to 3,200 kb. Much less LD was found among breeds. This advantageous LD structure holds great promise for reducing the interval of association in future mapping studies.

Genetic diversity in India and the inference of Eurasian population expansion

BackgroundGenetic studies of populations from the Indian subcontinent are of great interest because of India's large population size, complex demographic history, and unique social structure. Despite recent large-scale efforts in discovering human genetic variation, India's vast reservoir of genetic diversity remains largely unexplored.ResultsTo analyze an unbiased sample of genetic diversity in India and to investigate human migration history in Eurasia, we resequenced one 100-kb ENCODE region in 92 samples collected from three castes and one tribal group from the state of Andhra Pradesh in south India. Analyses of the four Indian populations, along with eight HapMap populations (692 samples), showed that 30% of all SNPs in the south Indian populations are not seen in HapMap populations. Several Indian populations, such as the Yadava, Mala/Madiga, and Irula, have nucleotide diversity levels as high as those of HapMap African populations. Using unbiased allele-frequency spectra, we investigated the expansion of human populations into Eurasia. The divergence time estimates among the major population groups suggest that Eurasian populations in this study diverged from Africans during the same time frame (approximately 90 to 110 thousand years ago). The divergence among different Eurasian populations occurred more than 40,000 years after their divergence with Africans.ConclusionsOur results show that Indian populations harbor large amounts of genetic variation that have not been surveyed adequately by public SNP discovery efforts. Our data also support a delayed expansion hypothesis in which an ancestral Eurasian founding population remained isolated long after the out-of-Africa diaspora, before expanding throughout Eurasia.