The peopling of South America and the trans-Andean gene flow of the first settlers.

Genetic and archaeological data indicate that the initial Paleoindian settlers of South America followed two entry routes separated by the Andes and the Amazon rainforest. The interactions between these paths and their impact on the peopling of South America remain unclear. Analysis of genetic variation in the Peruvian Andes and regions located south of the Amazon River might provide clues on this issue. We analyzed mitochondrial DNA variation at different Andean locations and >360,000 autosomal SNPs from 28 Native American ethnic groups to evaluate different trans-Andean demographic scenarios. Our data reveal that the Peruvian Altiplano was an important enclave for early Paleoindian expansions and point to a genetic continuity in the Andes until recent times, which was only marginally affected by gene flow from the Amazonian lowlands. Genomic variation shows a good fit with the archaeological evidence, indicating that the genetic interactions between the descendants of the settlers that followed the Pacific and Atlantic routes were extremely limited.

Cuba: exploring the history of admixture and the genetic basis of pigmentation using autosomal and uniparental markers.

We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.

Charting the Y-chromosome ancestry of present-day Argentinean Mennonites.

Old Order Mennonite communities initially arose in Northern Europe (centered in the Netherlands) and derived from the Anabaptist movement of the 16th century. Mennonites migrated to the New World in the early 18th century, first to North America, and more recently to Mesoamerica and South America. We analyzed Y-chromosome short tandem repeats (STRs) and single nucleotide polymorphisms in males from a community of Mennonites, 'La Nueva Esperanza', which arrived to Argentina in 1985 from colonies in Bolivia and Mexico. Molecular diversity indices coupled with demographic simulations show that Mennonites have a reduced variability when compared with local Argentinean populations and 69 European population samples. Mennonite Y-STR haplotypes were mainly observed in Central Europe. In agreement, multidimensional scaling analyses based on RST genetic distances indicate that Mennonite Y-chromosomes are closely related to Central/Northern Europeans (the Netherlands, Switzerland and Denmark). In addition, statistical inferences made on the most likely geographic origin of Y-chromosome haplotypes point more specifically to the Netherlands as the populations that best represent the majority of the Mennonite Y-chromosomes. Overall, Y-chromosome variation of Mennonites shows the signatures of moderate reduction of variability when compared with source populations, which is in good agreement with their lifestyle in small endogamous demes. These genetic singularities could also help to understand disease conditions that are more prevalent among Mennonites.

Genomic insights on the ethno-history of the Maya and the 'Ladinos' from Guatemala.

BACKGROUND: Guatemala is a multiethnic and multilingual country located in Central America. The main population groups separate 'Ladinos' (mixed Native American-African-Spanish), and Native indigenous people of Maya descent. Among the present-day Guatemalan Maya, there are more than 20 different ethnic groups separated by different languages and cultures. Genetic variation of these communities still remains largely unexplored. The principal aim of this study is to explore the genetic variability of the Maya and 'Ladinos' from Guatemala by means of uniparental and ancestry informative markers (AIMs). RESULTS: Analyses of uniparental genetic markers indicate that Maya have a dominant Native American ancestry (mitochondrial DNA [mtDNA]: 100%; Y-chromosome: 94%). 'Ladino', however, show a clear gender-bias as indicated by the large European ancestry observed in the Y-chromosome (75%) compared to the mtDNA (0%). Autosomal polymorphisms (AIMS) also mirror this marked gender-bias: (i) Native American ancestry: 92% for the Maya vs. 55% for the 'Ladino', and (ii) European ancestry: 8% for the Maya vs. 41% for the 'Ladino'. In addition, the impact of the Trans-Atlantic slave trade on the present-day Guatemalan population is very low (and only occurs in the 'Ladino'; mtDNA: 9%; AIMS: 4%), in part mirroring the fact that Guatemala has a predominant orientation to the Pacific Ocean instead of a Caribbean one. Sequencing of entire Guatemalan mitogenomes has led to improved Native American phylogeny via the addition of new haplogroups that are mainly observed in Mesoamerica and/or the North of South America. CONCLUSIONS: The data reveal the existence of a fluid gene flow in the Mesoamerican area and a predominant unidirectional flow towards South America, most likely occurring during the Pre-Classic (1800 BC-200 AD) and the Classic (200-1000 AD) Eras of the Mesoamerican chronology, coinciding with development of the most distinctive and advanced Mesoamerican civilization, the Maya. Phylogenetic features of mtDNA data also suggest a demographic scenario that is compatible with moderate local endogamy and isolation in the Maya combined with episodes of gene exchange between ethnic groups, suggesting an ethno-genesis in the Guatemalan Maya that is recent and supported on a cultural rather than a biological basis.

Mosaic maternal ancestry in the Great Lakes region of East Africa.

The Great Lakes lie within a region of East Africa with very high human genetic diversity, home of many ethno-linguistic groups usually assumed to be the product of a small number of major dispersals. However, our knowledge of these dispersals relies primarily on the inferences of historical, linguistics and oral traditions, with attempts to match up the archaeological evidence where possible. This is an obvious area to which archaeogenetics can contribute, yet Uganda, at the heart of these developments, has not been studied for mitochondrial DNA (mtDNA) variation. Here, we compare mtDNA lineages at this putative genetic crossroads across 409 representatives of the major language groups: Bantu speakers and Eastern and Western Nilotic speakers. We show that Uganda harbours one of the highest mtDNA diversities within and between linguistic groups, with the various groups significantly differentiated from each other. Despite an inferred linguistic origin in South Sudan, the data from the two Nilotic-speaking groups point to a much more complex history, involving not only possible dispersals from Sudan and the Horn but also large-scale assimilation of autochthonous lineages within East Africa and even Uganda itself. The Eastern Nilotic group also carries signals characteristic of West-Central Africa, primarily due to Bantu influence, whereas a much stronger signal in the Western Nilotic group suggests direct West-Central African ancestry. Bantu speakers share lineages with both Nilotic groups, and also harbour East African lineages not found in Western Nilotic speakers, likely due to assimilating indigenous populations since arriving in the region ~3000 years ago.

No association between typical European mitochondrial variation and prostate cancer risk in a Spanish cohort.

Mitochondrial common variants (mtSNPs) and the haplogroups defined by them have been inconsistently correlated with increased prostate cancer risk. Here we aimed to investigate the influence of the mitochondrial genetic background on prostate cancer. A total of 15 single-nucleotide polymorphisms (SNPs) representing the common European branches of the mtDNA phylogeny were analyzed in a cohort of 620 Spanish prostate cancer patients and 616 matched population-based controls. Association tests were computed on mtSNPs and haplogroups. None of the evaluated mtSNPs or haplogroups were statistically associated with prostate cancer risk in our Spanish cohort. We show that previous association findings do not rest on solid grounds given that all of them (i) were based on underpowered studies, (ii) did not control for population stratification, (iii) lacked replication/confirmation cohorts, and (iv) and did not control for multiple test corrections. Taken together, a critical reassessment of the previous literature and the results obtained in the present study suggest that mtDNA common European variants are not correlated with increases in the risk for prostate cancer.

A melting pot of multicontinental mtDNA lineages in admixed Venezuelans.

The arrival of Europeans in Colonial and post-Colonial times coupled with the forced introduction of sub-Saharan Africans have dramatically changed the genetic background of Venezuela. The main aim of the present study was to evaluate, through the study of mitochondrial DNA (mtDNA) variation, the extent of admixture and the characterization of the most likely continental ancestral sources of present-day urban Venezuelans. We analyzed two admixed populations that have experienced different demographic histories, namely, Caracas (n = 131) and Pueblo Llano (n = 219). The native American component of admixed Venezuelans accounted for 80% (46% haplogroup [hg] A2, 7% hg B2, 21% hg C1, and 6% hg D1) of all mtDNAs; while the sub-Saharan and European contributions made up ∼10% each, indicating that Trans-Atlantic immigrants have only partially erased the native American nature of Venezuelans. A Bayesian-based model allowed the different contributions of European countries to admixed Venezuelans to be disentangled (Spain: ∼38.4%, Portugal: ∼35.5%, Italy: ∼27.0%), in good agreement with the documented history. Seventeen entire mtDNA genomes were sequenced, which allowed five new native American branches to be discovered. B2j and B2k, are supported by two different haplotypes and control region data, and their coalescence ages are 3.9 k.y. (95% C.I. 0-7.8) and 2.6 k.y. (95% C.I. 0.1-5.2), respectively. The other clades were exclusively observed in Pueblo Llano and they show the fingerprint of strong recent genetic drift coupled with severe historical consanguinity episodes that might explain the high prevalence of certain Mendelian and complex multi-factorial diseases in this region.

Reassessing the role of mitochondrial DNA mutations in autism spectrum disorder.

BACKGROUND: There is increasing evidence that impairment of mitochondrial energy metabolism plays an important role in the pathophysiology of autism spectrum disorders (ASD; OMIM number: 209850). A significant proportion of ASD cases display biochemical alterations suggestive of mitochondrial dysfunction and several studies have reported that mutations in the mitochondrial DNA (mtDNA) molecule could be involved in the disease phenotype. METHODS: We analysed a cohort of 148 patients with idiopathic ASD for a number of mutations proposed in the literature as pathogenic in ASD. We also carried out a case control association study for the most common European haplogroups (hgs) and their diagnostic single nucleotide polymorphisms (SNPs) by comparing cases with 753 healthy and ethnically matched controls. RESULTS: We did not find statistical support for an association between mtDNA mutations or polymorphisms and ASD. CONCLUSIONS: Our results are compatible with the idea that mtDNA mutations are not a relevant cause of ASD and the frequent observation of concomitant mitochondrial dysfunction and ASD could be due to nuclear factors influencing mitochondrion functions or to a more complex interplay between the nucleus and the mitochondrion/mtDNA.

The impact of modern migrations on present-day multi-ethnic Argentina as recorded on the mitochondrial DNA genome.

BACKGROUND: The genetic background of Argentineans is a mosaic of different continental ancestries. From colonial to present times, the genetic contribution of Europeans and sub-Saharan Africans has superposed to or replaced the indigenous genetic 'stratum'. A sample of 384 individuals representing different Argentinean provinces was collected and genotyped for the first and the second mitochondrial DNA (mtDNA) hypervariable regions, and selectively genotyped for mtDNA SNPs. This data was analyzed together with additional 440 profiles from rural and urban populations plus 304 from Native American Argentineans, all available from the literature. A worldwide database was used for phylogeographic inferences, inter-population comparisons, and admixture analysis. Samples identified as belonging to hg (hg) H2a5 were sequenced for the entire mtDNA genome. RESULTS: Phylogenetic and admixture analyses indicate that only half of the Native American component in urban Argentineans might be attributed to the legacy of extinct ancestral Argentineans and that the Spanish genetic contribution is slightly higher than the Italian one. Entire H2a5 genomes linked these Argentinean mtDNAs to the Basque Country and improved the phylogeny of this Basque autochthonous clade. The fingerprint of African slaves in urban Argentinean mtDNAs was low and it can be phylogeographically attributed predominantly to western African. The European component is significantly more prevalent in the Buenos Aires province, the main gate of entrance for Atlantic immigration to Argentina, while the Native American component is larger in North and South Argentina. AMOVA, Principal Component Analysis and hgs/haplotype patterns in Argentina revealed an important level of genetic sub-structure in the country. CONCLUSIONS: Studies aimed to compare mtDNA frequency profiles from different Argentinean geographical regions (e.g., forensic and case-control studies) should take into account the important genetic heterogeneity of the country in order to prevent false positive claims of association in disease studies or inadequate evaluation of forensic evidence.

Minisequencing mitochondrial DNA pathogenic mutations.

BACKGROUND: There are a number of well-known mutations responsible of common mitochondrial DNA (mtDNA) diseases. In order to overcome technical problems related to the analysis of complete mtDNA genomes, a variety of different techniques have been proposed that allow the screening of coding region pathogenic mutations. METHODS: We here propose a minisequencing assay for the analysis of mtDNA mutations. In a single reaction, we interrogate a total of 25 pathogenic mutations distributed all around the whole mtDNA genome in a sample of patients suspected for mtDNA disease. RESULTS: We have detected 11 causal homoplasmic mutations in patients suspected for Leber disease, which were further confirmed by standard automatic sequencing. Mutations m.11778G>A and m.14484T>C occur at higher frequency than expected by change in the Galician (northwest Spain) patients carrying haplogroup J lineages (Fisher's Exact test, P-value < 0.01). The assay performs well in mixture experiments of wild:mutant DNAs that emulate heteroplasmic conditions in mtDNA diseases. CONCLUSION: We here developed a minisequencing genotyping method for the screening of the most common pathogenic mtDNA mutations which is simple, fast, and low-cost. The technique is robust and reproducible and can easily be implemented in standard clinical laboratories.

Multiplex SNaPshot for detection of BRCA1/2 common mutations in Spanish and Spanish related breast/ovarian cancer families.

BACKGROUND: It is estimated that 5-10% of all breast cancer are hereditary and attributable to mutations in the highly penetrance susceptibility genes BRCA1 and BRCA2. The genetic analysis of these genes is complex and expensive essentially because their length. Nevertheless, the presence of recurrent and founder mutations allows a pre-screening for the identification of the most frequent mutations found in each geographical region. In Spain, five mutations in BRCA1 and other five in BRCA2 account for approximately 50% of the mutations detected in Spanish families. METHODS: We have developed a novel PCR multiplex SNaPshot reaction that targets all ten recurrent and founder mutations identified in BRCA1 and BRCA2 in Spain to date. RESULTS: The SNaPshot reaction was performed on samples previously analyzed by direct sequencing and all mutations were concordant. This strategy permits the analysis of approximately 50% of all mutations observed to be responsible for breast/ovarian cancer in Spanish families using a single reaction per patient sample. CONCLUSION: The SNaPshot assay developed is sensitive, rapid, with minimum cost per sample and additionally can be automated for high-throughput genotyping. The SNaPshot assay outlined here is not only useful for analysis of Spanish breast/ovarian cancer families, but also e.g. for populations with Spanish ancestry, such as those in Latin America.

Revealing latitudinal patterns of mitochondrial DNA diversity in Chileans.

The territory of Chile is particularly long and narrow, which combined with its mountainous terrain, makes it a unique scenario for human genetic studies. We obtained 995 control region mitochondrial DNA (mtDNA) sequences from Chileans representing populations living at different latitudes of the country from the North to the southernmost region. The majority of the mtDNA profiles are of Native American origin (∼88%). The remaining haplotypes are mostly of recent European origin (∼11%), and only a minor proportion is of recent African ancestry (∼1%). While these proportions are relatively uniform across the country, more structured patterns of diversity emerge when examining the variation from a phylogeographic perspective. For instance, haplogroup A2 reaches ∼9% in the North, and its frequency decreases gradually to ∼1% in the southernmost populations, while the frequency of haplogroup D (sub-haplogroups D1 and D4) follows the opposite pattern: 36% in the southernmost region, gradually decreasing to 21% in the North. Furthermore, there are remarkable signatures of founder effects in specific sub-clades of Native American (e.g. haplogroups D1j and D4p) and European (e.g. haplogroups T2b3 and K1a4a1a+195) ancestry. We conclude that the magnitude of the latitudinal differences observed in the patterns of mtDNA variation might be relevant in forensic casework.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula.

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

SNaPshot typing of mitochondrial DNA coding region variants.

We describe a rapid and robust assay to genotype mitochondrial deoxyribonucleic acid (mtDNA) coding region single nucleotide polymorphism (SNPs) using the SNaPshot (Applied Biosystems, Foster City, CA) minisequencing reaction kit. A protocol for mtDNA SNaPshot typing is described in detail, although we emphasize that this method allows great flexibility in the implementation of whatever set of mtDNA SNPs. We discuss the utility of our selection of mtDNA SNPs for molecular anthropologists and forensic geneticists. Firstly, these SNPs allow allocating common mitochondrial West Eurasian haplotypes into their corresponding branches of the mtDNA skeleton, with special attention to the subdivision of sequences belonging to haplogroup H, the most frequent European haplogroup (40-50%) and the worst phylogenetically characterized in the first and second hypervariable segments (HVS-I/II; by far, the most common segments analyzed by sequencing). Second, the polymorphic positions selected for this multiplex reaction considerably increase the discrimination power of current mitochondrial analysis in the forensic field. The method shows high accuracy and robustness, avoiding both the use of alternative time-consuming classical strategies (i.e., restriction fragment length polymorphism typing) and the requirement of high quantities of DNA template.

The Genomic Legacy of the Transatlantic Slave Trade in the Yungas Valley of Bolivia.

During the period of the Transatlantic Slave Trade (TAST) some enslaved Africans were forced to move to Upper Peru (nowadays Bolivia). At first they were sent to Potosí, but later to the tropical Yungas valley where the Spanish colonizers established a so-called "hacienda system" that was based on slave labor, including African-descendants. Due to their isolation, very little attention has been paid so far to 'Afro-Bolivian' communities either within the research field of TAST or in genetic population studies. In this study, a total of 105 individuals from the Yungas were sequenced for their mitochondrial DNA (mtDNA) control region, and mitogenomes were obtained for a selected subset of these samples. We also genotyped 46 Ancestry Informative Markers (AIM) in order to investigate continental ancestry at the autosomal level. In addition, Y-chromosome STR and SNP data for a subset of the same individuals was also available from the literature. The data indicate that the partitioning of mtDNA ancestry in the Yungas differs significantly from that in the rest of the country: 81% Native American, 18% African, and 1% European. Interestingly, the great majority of 'Afro-descendant' mtDNA haplotypes in the Yungas (84%) concentrates in the locality of Tocaña. This high proportion of African ancestry in the Tocaña is also manifested in the Y-chromosome (44%) and in the autosomes (56%). In sharp contrast with previous studies on the TAST, the ancestry of about 1/3 of the 'Afro-Bolivian' mtDNA haplotypes can be traced back to East and South East Africa, which may be at least partially explained by the Arab slave trade connected to the TAST.

The multiethnic ancestry of Bolivians as revealed by the analysis of Y-chromosome markers.

We have analyzed the specific male genetic component of 226 Bolivians recruited in five different regions ("departments"), La Paz, Cochabamba, Pando, Beni, and Santa Cruz. To evaluate the effect of geography on the distribution of genetic variability, the samples were also grouped into three main eco-geographical regions, namely, Andean, Sub-Andean, and Llanos. All the individuals were genotyped for 17 Y-STR and 32 Y-SNP markers. The average Y-chromosome Native American component in Bolivians is 28%, and it is mainly represented by haplogroup Q1a3a, while the average Y-chromosome European ancestry is 65%, and it is mainly represented by haplogroup R1b1-P25. The data indicate that there exists significant population sub-division in the country in terms of continental ancestry. Thus, the partition of ancestries in Llanos, Sub-Andean, and Andean regions is as follows (respectively): (i) Native American ancestry: 47%, 7%, and 19%, (ii) European ancestry: 46%, 86%, and 75%, and (iii) African ancestry: 7%, 7%, and 6%. The population sub-structure in the country is also well mirrored when inferred from an AMOVA analysis, indicating that among-population variance in the country reaches 9.74-11.15%. This suggests the convenience of using regional datasets for forensic applications in Bolivia, instead of using a global and single country database. By comparing the Y-chromosome patterns with those previously reported on the same individuals on autosomal SNPs and mitochondrial DNA (mtDNA), it becomes clear that Bolivians show a strong gender-bias.

Correction: Uniparental Markers of Contemporary Italian Population Reveals Details on Its Pre-Roman Heritage.

[This corrects the article DOI: 10.1371/journal.pone.0050794.].

Ancestry analysis reveals a predominant Native American component with moderate European admixture in Bolivians.

We have genotyped 46 Ancestry Informative Markers (AIMs) in two of the most populated areas in Bolivia, namely, La Paz (Andean region; n=105), and Chuquisaca (Sub-Andean region; n=73). Using different analytical tools, we inferred admixture proportions of these two American communities by comparing the genetic profiles with those publicly available from the CEPH (Centre d'Etude du Polymorphisme Humain) panel representing three main continental groups (Africa, Europe, and America). By way of simulations, we first evaluated the minimum sample size needed in order to obtain accurate estimates of ancestry proportions. The results indicated that sample sizes above 30 individuals could be large enough to estimate main continental ancestry proportions using the 46 AIMs panel. With the exception of a few individuals, the results also indicated that Bolivians showed a predominantly Native American ancestry with variable levels of European admixture. The proportions of ancestry were statistically different in La Paz and Chuquisaca: the Native American component was 86% and 77% (Mann-Whitney U-test: un-adjusted P-value=2.1×10(-5)), while the European ancestry was 13% and 21% (Mann-Whitney U-test: un-adjusted P-value=3.6×10(-5)), respectively. The African ancestry in Bolivians captured by the AIMs analyzed in the present study was below 2%. The inferred ancestry of Bolivians fits well with previous studies undertaken on haplotype data, indicating a major proportion of Native American lineages. The genetic differences observed in these two groups suggest that forensic genetic analysis should be better performed based on local databases built in the main Bolivian areas.

The genetic legacy of the pre-colonial period in contemporary Bolivians.

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.