Land, language, and loci: mtDNA in Native Americans and the genetic history of Peru

Insights into the genetic histories and lifeways of Machu Picchu's occupants

Machu Picchu originally functioned as a palace within the estate of the Inca emperor Pachacuti between ~1420 and 1532 CE. Before this study, little was known about the people who lived and died there, where they came from or how they were related to the inhabitants of the Inca capital of Cusco. We generated genome-wide data for 34 individuals buried at Machu Picchu who are believed to have been retainers or attendants assigned to serve the Inca royal family, as well as 34 individuals from Cusco for comparative purposes. When the ancient DNA results are contextualized using historical and archaeological data, we conclude that the retainer population at Machu Picchu was highly heterogeneous with individuals exhibiting genetic ancestries associated with groups from throughout the Inca Empire and Amazonia. The results suggest a diverse retainer community at Machu Picchu in which people of different genetic backgrounds lived, reproduced, and were interred together.

Sequencing errors in Native American mitogenomes: impact on clade definitions, haplogroup assignation, and beyond

Available evidence allows the interpretation that some cases of absence of otherwise expected variation, based on phylogenetic expectations in mitogenomes of Native American origin, are due to artificial recombination rather than to homoplasy, while other more complex scenarios involving combination of original Cambridge Reference Sequence mistakes plus incomplete or incorrect scoring of variation are also showed. Several instances of mismatched control and coding regions as well as partially duplicated HV2 are observed in Peruvians, while intra-haplogroup chimaeras of different D1 subhaplogroups are referred to in Mexican Native Americans. A revised definition for haplogroup B2h is proposed, and preventive quality control measures are suggested.

Uncovering signals of positive selection in Peruvian populations from three ecological regions

Perú hosts extremely diverse ecosystems which can be broadly classified into three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide SNP-array data from 286 Peruvians living across the three ecoregions and analysed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44 and ICAM1 in the coast; CBLB and PRDM1 in rainforest and the BRD2- HLA-DOA- HLA-DPA1 region in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

HLA-G 3'UTR haplotype frequencies in highland and lowland South Native American populations

Human Leukocyte Antigen (HLA)-G participates in several biological processes, including reproduction, vascular remodeling, immune tolerance, and hypoxia response. HLA-G is a potential candidate gene for high altitude adaptation since its expression is modulated in both micro and macro environment under hypoxia and constant cellular stress. Besides the promoter region, the HLA-G 3'untranslated region (UTR) influences HLA-G expression patterns through several post-transcriptional mechanisms. Currently, the 3'UTR genetic diversity in terms of altitude adaptation of Native American populations is still unexplored, particularly at high altitude ecoregions. Here, we evaluated 288 Native Americans from 9 communities located in the Andes [highland (HL); ≥2,500 m (range = 2,838-4,433 m)] and 8 populations located in lowland (LL) regions [<2,500 m (range = 80-431 m); Amazonian tropical forest, Brazilian central plateau, and Chaco] of South America. In total, nine polymorphic sites and ten haplotypes were observed. The most frequent haplotypes (UTR-1, UTR-2, and UTR-3) accounted for ∼ 77% of haplotypes found in LL, while in the HL, the same haplotypes reach ∼ 93%. Also, a remarkable high frequency of putative ancestral UTR-5 haplotype was observed in LL (21.5%), while in HL UTR-2 reaches up to 47%. Further, UTR-2 frequency positively correlates with altitude-related variables, while a negative correlation for UTR-5 was observed. From an evolutionary perspective, we observed a tendency towards balancing selection in HL and LL populations thus suggesting that haplotypes of ancient and more derived alleles may have been co-opted for relatively recent adaptations such as those experienced by modern humans in the highland and lowland of South America. We also discuss how long-term balancing selection can be a reservoir of genetic variants that can be positively selected. Finally, our study provides some pieces of evidence that HLA-G 3'UTR haplotypes may have contributed to high altitude adaptation in the Andes.

Demographic history and selection at HLA loci in Native Americans

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

Dissecting the Pre-Columbian Genomic Ancestry of Native Americans along the Andes-Amazonia Divide

Extensive European and African admixture coupled with loss of Amerindian lineages makes the reconstruction of pre-Columbian history of Native Americans based on present-day genomes extremely challenging. Still open questions remain about the dispersals that occurred throughout the continent after the initial peopling from the Beringia, especially concerning the number and dynamics of diffusions into South America. Indeed, if environmental and historical factors contributed to shape distinct gene pools in the Andes and Amazonia, the origins of this East-West genetic structure and the extension of further interactions between populations residing along this divide are still not well understood.

To this end, we generated new high-resolution genome-wide data for 229 individuals representative of one Central and ten South Amerindian ethnic groups from Mexico, Peru, Bolivia, and Argentina. Low levels of European and African admixture in the sampled individuals allowed the application of fine-scale haplotype-based methods and demographic modeling approaches. These analyses revealed highly specific Native American genetic ancestries and great intragroup homogeneity, along with limited traces of gene flow mainly from the Andes into Peruvian Amazonians. Substantial amount of genetic drift differentially experienced by the considered populations underlined distinct patterns of recent inbreeding or prolonged isolation. Overall, our results support the hypothesis that all non-Andean South Americans are compatible with descending from a common lineage, while we found low support for common Mesoamerican ancestors of both Andeans and other South American groups. These findings suggest extensive back-migrations into Central America from non-Andean sources or conceal distinct peopling events into the Southern Continent.

Análisis de ADN mitocondrial en una muestra de restos óseos precolombinos de Norte de Santander, Colombia (Área Cultural Chitarera)

Los análisis de ADN antiguo (ADNa) han incrementado en los últimos años permitiendo conocer la diversidad genética de las poblaciones precolombinas. En Colombia, existen pocos registros arqueológicos de la población prehispánica del Norte de Santander habitada en el siglo XVI por el grupo Chitarero. Por este motivo, nos propusimos analizar la diversidad genética a partir de secuencias de la región HVRI del ADNmt y determinar sus posibles relaciones con otras comunidades tanto antiguas como contemporáneas. Se analizaron siete individuos precolombinos asociados a este grupo prehispánico, recuperados en los municipios de Cácota y Silos en el departamento de Norte de Santander de los Andes Orientales colombianos, siguiendo criterios estrictos de autenticidad para el ADNa. En todos los individuos se logró identificar el haplogrupo B caracterizado por el polimorfismo en la posición 16217C, siendo éste uno de los más frecuentes en comunidades precolombinas y contemporáneas de los Andes Suramericanos. Este hallazgo indica que este grupo poblacional se encuentra estrechamente emparentado por línea materna, con posibles índices de endogamia, con una probable densidad demográfica baja y una baja diversidad genética, similares a lo observado en comunidades pertenecientes a periodos anteriores como el Formativo. Este grupo precolombino exhibe una de las diversidades genéticas más bajas reportadas en las poblaciones pertenecientes a la familia lingüistica Chibcha. Estos resultados genéticos coinciden con los planteamientos sobre el grupo Chitarero de pertenecer a comunidades pequeñas independientes, con asentamientos dispersos, apartados unos de otros. 

Evolutionary genomic dynamics of Peruvians before, during, and after the Inca Empire

Through the Peruvian Genome Project we generate and analyze the genomes of 280 individuals where the majority have >90% Native American ancestry and explore questions at the interface of evolutionary genetics, history, anthropology, and medicine. This is the most extensive sampling of high-coverage Native American and mestizo whole genomes to date. We estimate an initial peopling of Peru was rapid and began by 12,000 y ago. In addition, the mestizo populations exhibit admixture between Native American groups prior to their Spanish admixture and was likely influenced by the Inca Empire and Spanish conquest. Our results address important Native American population history questions and establish a dataset beneficial to address the underrepresentation of Native American ancestry in sequencing studies.

Native Americans from the Amazon, Andes, and coastal geographic regions of South America have a rich cultural heritage but are genetically understudied, therefore leading to gaps in our knowledge of their genomic architecture and demographic history. In this study, we sequence 150 genomes to high coverage combined with an additional 130 genotype array samples from Native American and mestizo populations in Peru. The majority of our samples possess greater than 90% Native American ancestry, which makes this the most extensive Native American sequencing project to date. Demographic modeling reveals that the peopling of Peru began ∼12,000 y ago, consistent with the hypothesis of the rapid peopling of the Americas and Peruvian archeological data. We find that the Native American populations possess distinct ancestral divisions, whereas the mestizo groups were admixtures of multiple Native American communities that occurred before and during the Inca Empire and Spanish rule. In addition, the mestizo communities also show Spanish introgression largely following Peruvian Independence, nearly 300 y after Spain conquered Peru. Further, we estimate migration events between Peruvian populations from all three geographic regions with the majority of between-region migration moving from the high Andes to the low-altitude Amazon and coast. As such, we present a detailed model of the evolutionary dynamics which impacted the genomes of modern-day Peruvians and a Native American ancestry dataset that will serve as a beneficial resource to addressing the underrepresentation of Native American ancestry in sequencing studies.

Population genetic structure of traditional populations in the Peruvian Central Andes and implications for South American population history

Molecular-based characterizations of Andean peoples are traditionally conducted in the service of elucidating continent-level evolutionary processes in South America. Consequently, genetic variation among "western" Andean populations is often represented in relation to variation among "eastern" Amazon and Orinoco River Basin populations. This west-east contrast in patterns of population genetic variation is typically attributed to large-scale phenomena, such as dual founder colonization events or differing long-term microevolutionary histories. However, alternative explanations that consider the nature and causes of population genetic diversity within the Andean region remain underexplored. Here we examine population genetic diversity in the Peruvian Central Andes using data from the mtDNA first hypervariable region and Y-chromosome short tandem repeats among 17 newly sampled populations and 15 published samples. Using this geographically comprehensive data set, we first reassessed the currently accepted pattern of western versus eastern population genetic structure, which our results ultimately reject: mtDNA population diversities were lower, rather than higher, within Andean versus eastern populations, and only highland Y-chromosomes exhibited significantly higher within-population diversities compared with eastern groups. Multiple populations, including several highland samples, exhibited low genetic diversities for both genetic systems. Second, we explored whether the implementation of Inca state and Spanish colonial policies starting at about ad 1400 could have substantially restructured population genetic variation and consequently constitute a primary explanation for the extant pattern of population diversity in the Peruvian Central Andes. Our results suggest that Peruvian Central Andean population structure cannot be parsimoniously explained as the sole outcome of combined Inca and Spanish policies on the region's population demography: highland populations differed from coastal and lowland populations in mtDNA genetic structure only; highland groups also showed strong evidence of female-biased gene flow and/or effective sizes relative to other Peruvian ecozones. Taken together, these findings indicate that population genetic structure in the Peruvian Central Andes is considerably more complex than previously reported and that characterizations of and explanations for genetic variation may be best pursued within more localized regions and defined time periods.

Mitochondrial DNA variability among six South American Amerindian villages from the Pano linguistic group

Although scattered throughout a large geographic area, the members of the Pano linguistic group present strong ethnic, linguistic, and cultural homogeneity, a feature that causes them to be considered components of a same "Pano" tribe. Nevertheless, the genetic homogeneity between Pano villages has not yet been examined. To study the genetic structure of the Pano linguistic group, four major Native American mitochondrial DNA (mtDNA) founder haplogroups were analyzed in 77 Amerindians from six villages of four Pano tribes (Katukina, Kaxináwa, Marúbo, and Yaminawa) located in the Brazilian Amazon. The central position of these tribes in the continent makes them relevant for attempts to reconstruct population movements in South America. Except for a single individual that presented an African haplogroup L, all remaining individuals presented one of the four Native American haplogroups. Significant heterogeneity was observed across the six Pano villages. Although Amerindian populations are usually characterized by considerable interpopulational diversity, the high heterogeneity level observed is unexpected if the strong ethnic, linguistic, and cultural homogeneity of the Pano linguistic group is taken into account. The present findings indicate that the ethnic, linguistic, and cultural homogeneity does not imply genetic homogeneity. Even though the genetic heterogeneity uncovered may be a female-specific process, the most probable explanation for that is the joint action of isolation and genetic drift as major factors influencing the genetic structure of the Pano linguistic group.

Genetic diversity of a late prehispanic group of the Quebrada de Humahuaca, northwestern Argentina

This palaeogenetic study focused on the analysis of a late prehispanic Argentinean group from the Humahuaca valley, with the main aim of reconstructing its (micro)evolutionary history. The Humahuaca valley, a natural passageway from the eastern plains to the highlands, was the living environment of Andean societies whose cultural but especially biological diversity is still poorly understood. We analyzed the DNA extracted from 39 individuals who populated this upper valley during the Regional Development period (RDP) (between the 11th and 15th centuries CE), to determine their maternal and paternal genetic ancestry. Some mitochondrial and Y-chromosomal haplotypes specific to the Andean region are consistent with an origin in the highlands of Central Andes. On the other hand, a significant genetic affinity with contemporary admixed communities of the Chaco area was detected. Expectedly, recent demographic events, such as the expansion of the Inca Empire or the European colonization, have changed the original mitochondrial gene pool of the ancient Humahuaca Valley community. Finally, we identified a particular geographical organization of the prehispanic populations of Northwestern Argentina. Our results suggest that the communities of the region were divided between two different spheres of interaction, which is consistent with assumptions made by means of craniometric traits.

Mitochondrial variation among the Aymara and the signatures of population expansion in the central Andes

OBJECTIVES

The exploitation of marine resources and intensive agriculture led to a marked population increase early in central Andean prehistory. Constant historic and prehistoric population movements also characterize this region. These features undoubtedly affected regional genetic variation, but the exact nature of these effects remains uncertain.

METHODS

Mitochondrial DNA (mtDNA) hypervariable region I sequence variation in 61 Aymara individuals from La Paz, Bolivia, was analyzed and compared to sequences from 47 other South American populations to test hypotheses of whether increased female effective population size and gene flow influenced the mtDNA variation among central Andean populations.

RESULTS

The Aymara and Quechua were genetically diverse showing evidence of population expansion and large effective population size, and a demographic expansion model fits the mtDNA variation found among central Andean populations well. Estimated migration rates and the results of AMOVA and multidimensional scaling analysis suggest that female gene flow was also an important factor, influencing genetic variation among the central Andeans as well as lowland populations from western South America. mtDNA variation in south central Andes correlated better with geographic proximity than with language, and fit a population continuity model.

CONCLUSION

The mtDNA data suggests that the central Andeans experienced population expansion, most likely because of rapid demographic expansion after introduction of intensive agriculture, but roles of female gene flow need to be further explored.

Identification of Polynesian mtDNA haplogroups in remains of Botocudo Amerindians from Brazil

There is a consensus that modern humans arrived in the Americas 15,000-20,000 y ago during the Late Pleistocene, most probably from northeast Asia through Beringia. However, there is still debate about the time of entry and number of migratory waves, including apparent inconsistencies between genetic and morphological data on Paleoamericans. Here we report the identification of mitochondrial sequences belonging to haplogroups characteristic of Polynesians in DNA extracted from ancient skulls of the now extinct Botocudo Indians from Brazil. The identification of these two Polynesian haplogroups was confirmed in independent replications in Brazil and Denmark, ensuring reliability of the data. Parallel analysis of 12 other Botocudo individuals yielded only the well-known Amerindian mtDNA haplogroup C1. Potential scenarios to try to help understand these results are presented and discussed. The findings of this study may be relevant for the understanding of the pre-Columbian and/or post-Columbian peopling of the Americas.

Genetic differences between Chibcha and Non-Chibcha speaking tribes based on mitochondrial DNA (mtDNA) haplogroups from 21 Amerindian tribes from Colombia

We analyzed the frequency of four mitochondrial DNA haplogroups in 424 individuals from 21 Colombian Amerindian tribes. Our results showed a high degree of mtDNA diversity and genetic heterogeneity. Frequencies of mtDNA haplogroups A and C were high in the majority of populations studied. The distribution of these four mtDNA haplogroups from Amerindian populations was different in the northern region of the country compared to those in the south. Haplogroup A was more frequently found among Amerindian tribes in northern Colombia, while haplogroup D was more frequent among tribes in the south. Haplogroups A, C and D have clinal tendencies in Colombia and South America in general. Populations belonging to the Chibcha linguistic family of Colombia and other countries nearby showed a strong genetic differentiation from the other populations tested, thus corroborating previous findings. Genetically, the Ingano, Paez and Guambiano populations are more closely related to other groups of south eastern Colombia, as also inferred from other genetic markers and from archeological data. Strong evidence for a correspondence between geographical and linguistic classification was found, and this is consistent with evidence that gene flow and the exchange of customs and knowledge and language elements between groups is facilitated by close proximity.

Mitochondrial DNA control region data reveal high prevalence of Native American lineages in Jujuy province, NW Argentina

Mitochondrial control region (16024-576) sequences were generated from 180 individuals of four population nuclei from the province of Jujuy (NW Argentina), located at different altitudes above sea level. The frequency at which a randomly selected mtDNA profile would be expected to occur in the general population (random match probability) was estimated at 0.011, indicating a relatively high diversity. Analysis of the haplogroup distribution revealed that Native American lineages A2 (13.9%), B (56.7%), C1 (17.8%), D1 (8.9%) and D4h3a (1.1%) accounted for more than 98% of the total mtDNA haplogroup diversity in the sample examined. We detected a certain degree of genetic heterogeneity between two subpopulations located at different points along the altitudinal gradient (Valles and Puna), suggesting that altitude above sea level cannot be ruled out as a factor promoting divergences in mtDNA haplogroup frequencies, since altitude is closely associated with human living conditions, and consequently, with low demographic sizes and the occurrence of genetic drift processes in human communities. In all, mitochondrial DNA database obtained for Jujuy province strongly points to the need for creating local mtDNA databases, to avoid bias in forensic estimations caused by genetic substructuring of the populations.

An alternative model for the early peopling of southern South America revealed by analyses of three mitochondrial DNA haplogroups

After several years of research, there is now a consensus that America was populated from Asia through Beringia, probably at the end of the Pleistocene. But many details such as the timing, route(s), and origin of the first settlers remain uncertain. In the last decade genetic evidence has taken on a major role in elucidating the peopling of the Americas. To study the early peopling of South America, we sequenced the control region of mitochondrial DNA from 300 individuals belonging to indigenous populations of Chile and Argentina, and also obtained seven complete mitochondrial DNA sequences. We identified two novel mtDNA monophyletic clades, preliminarily designated B2l and C1b13, which together with the recently described D1g sub-haplogroup have locally high frequencies and are basically restricted to populations from the extreme south of South America. The estimated ages of D1g and B2l, about ~15,000 years BP, together with their similar population dynamics and the high haplotype diversity shown by the networks, suggests that they probably appeared soon after the arrival of the first settlers and agrees with the dating of the earliest archaeological sites in South America (Monte Verde, Chile, 14,500 BP). One further sub-haplogroup, D4h3a5, appears to be restricted to Fuegian-Patagonian populations and reinforces our hypothesis of the continuity of the current Patagonian populations with the initial founders. Our results indicate that the extant native populations inhabiting South Chile and Argentina are a group which had a common origin, and suggest a population break between the extreme south of South America and the more northern part of the continent. Thus the early colonization process was not just an expansion from north to south, but also included movements across the Andes.

Uniparental genetic markers in South Amerindians

A comprehensive review of uniparental systems in South Amerindians was undertaken. Variability in the Y-chromosome haplogroups were assessed in 68 populations and 1,814 individuals whereas that of Y-STR markers was assessed in 29 populations and 590 subjects. Variability in the mitochondrial DNA (mtDNA) haplogroup was examined in 108 populations and 6,697 persons, and sequencing studies used either the complete mtDNA genome or the highly variable segments 1 and 2. The diversity of the markers made it difficult to establish a general picture of Y-chromosome variability in the populations studied. However, haplogroup Q1a3a* was almost always the most prevalent whereas Q1a3* occurred equally in all regions, which suggested its prevalence among the early colonizers. The STR allele frequencies were used to derive a possible ancient Native American Q-clade chromosome haplotype and five of six STR loci showed significant geographic variation. Geographic and linguistic factors moderately influenced the mtDNA distributions (6% and 7%, respectively) and mtDNA haplogroups A and D correlated positively and negatively, respectively, with latitude. The data analyzed here provide rich material for understanding the biological history of South Amerindians and can serve as a basis for comparative studies involving other types of data, such as cultural data.

Ancient DNA reveals kinship burial patterns of a pre-Columbian Andean community

BACKGROUND

A detailed genetic study of the pre-Columbian population inhabiting the Tompullo 2 archaeological site (department Arequipa, Peru) was undertaken to resolve the kin relationships between individuals buried in six different chullpas. Kin relationships were an important factor shaping the social organization in the pre-Columbian Andean communities, centering on the ayllu, a group of relatives that shared a common land and responsibilities. The aim of this study was to evaluate whether this Andean model of a social organization had an influence on mortuary practices, in particular to determine whether chullpas served as family graves.

RESULTS

The remains of forty-one individuals were analyzed with both uniparental (mtDNA, Y-chromosome) and biparental (autosomal microsatellites) markers. Reproducible HVRI sequences, autosomal and Y chromosomal STR profiles were obtained for 24, 16 and 11 individuals, respectively. Mitochondrial DNA diversity was comparable to that of ancient and contemporary Andean populations. The Tompullo 2 population exhibited the closest relationship with the modern population from the same region. A kinship analysis revealed complex pattern of relations within and between the graves. However mean relatedness coefficients regarding the pairs of individuals buried in the same grave were significantly higher than those regarding pairs buried in different graves. The Y chromosome profiles of 11 males suggest that only members of one male line were buried in the same grave.

CONCLUSIONS

Genetic investigation of the population that inhabited Tompullo 2 site shows continuity between pre-Columbian and modern Native Amerindian populations inhabiting the Arequipa region. This suggests that no major demographic processes have influenced the mitochondrial DNA diversity of these populations during the past five hundred years. The kinship analysis involving uni- and biparental markers suggests that the community that inhabited the Tompullo 2 site was organized into extended family groups that were buried in different graves. This finding is in congruence with known models of social organization of Andean communities.

Rapid coastal spread of First Americans: novel insights from South America's Southern Cone mitochondrial genomes

It is now widely agreed that the Native American founders originated from a Beringian source population ∼15–18 thousand years ago (kya) and rapidly populated all of the New World, probably mainly following the Pacific coastal route. However, details about the migration into the Americas and the routes pursued on the continent still remain unresolved, despite numerous genetic, archaeological, and linguistic investigations. To examine the pioneering peopling phase of the South American continent, we screened literature and mtDNA databases and identified two novel mitochondrial DNA (mtDNA) clades, here named D1g and D1j, within the pan-American haplogroup D1. They both show overall rare occurrences but local high frequencies, and are essentially restricted to populations from the Southern Cone of South America (Chile and Argentina). We selected and completely sequenced 43 D1g and D1j mtDNA genomes applying highest quality standards. Molecular and phylogeographic analyses revealed extensive variation within each of the two clades and possibly distinct dispersal patterns. Their age estimates agree with the dating of the earliest archaeological sites in South America and indicate that the Paleo-Indian spread along the entire longitude of the American double continent might have taken even <2000 yr. This study confirms that major sampling and sequencing efforts are mandatory for uncovering all of the most basal variation in the Native American mtDNA haplogroups and for clarification of Paleo-Indian migrations, by targeting, if possible, both the general mixed population of national states and autochthonous Native American groups, especially in South America.

Candida theae sp. nov., a new anamorphic beverage-associated member of the Lodderomyces clade

Four strains representing a novel yeast species belonging to the genus Candida were independently isolated in Taiwan and Ecuador. Two strains (G17(T) and G31) were isolated in Taiwan, by pellet precipitation from plastic-bottled tea drinks produced in Indonesia, while two additional strains (CLQCA 10-049 and CLQCA 10-062) were recovered from ancient chicha fermentation vessels found in tombs in Quito, Ecuador. These four strains were morphologically, and phylogenetically identical to each other. No sexual reproduction was observed on common sporulation media. Large-subunit (LSU) rRNA gene sequence analysis revealed the four strains to belong to the Lodderomyces clade, closely related to members of the Candida parapsilosis species complex. The four strains, which have identical LSU D1/D2 sequences, differ from their closest phylogenetic neighbors, Candida orthopsilosis and Candida parapsilosis, by 6-9 nt substitutions, respectively. Physiologically, the four strains are similar to Candida parapsilosis, although they can be distinguished from their closest relative by the assimilation of arbutin, nitrite, and creatine. The Indonesian and Ecuadorian strain sets can also be distinguished from one another based on ITS sequencing, differing by 4 substitutions in ITS1 and 1 single nucleotide indel in ITS2. Collectively, the results indicate that the four strains represent a previously unrecognized species of Candida. The name Candida theae sp. nov. is proposed to accommodate these strains, with G-17(T) (BCRC 23242(T)=CBS 12239(T)=ATCC MYA-4746(T)) designated as the type strain.

Diachronic investigations of mitochondrial and Y-chromosomal genetic markers in pre-Columbian Andean highlanders from South Peru

This study examines the reciprocal effects of cultural evolution, and population dynamics in pre-Columbian southern Peru by the analysis of DNA from pre-Columbian populations that lived in the fringe area between the Andean highlands and the Pacific coast. The main objective is to reveal whether the transition from the Middle Horizon (MH: 650-1000 AD) to the Late Intermediate Period (LIP: 1000-1400 AD) was accompanied or influenced by population dynamic processes. Tooth samples from 90 individuals from several archaeological sites, dating to the MH and LIP, in the research area were collected to analyse mitochodrial, and Y-chromosomal genetic markers. Coding region polymorphisms were successfully analysed and replicated for 72 individuals, as were control region sequences for 65 individuals and Y-chromosomal single nucleotide polymorphisms (SNPs) for 19 individuals, and these were compared to a large set of ancient and modern indigenous South American populations. The diachronic comparison of the upper valley samples from both time periods reveals no genetic discontinuities accompanying the cultural dynamic processes. A high genetic affinity for other ancient and modern highland populations can be observed, suggesting genetic continuity in the Andean highlands at the latest from the MH. A significant matrilineal differentiation to ancient Peruvian coastal populations can be observed suggesting a differential population history.

Hierarchical modeling of genome-wide Short Tandem Repeat (STR) markers infers native American prehistory

This study examines a genome-wide dataset of 678 Short Tandem Repeat loci characterized in 444 individuals representing 29 Native American populations as well as the Tundra Netsi and Yakut populations from Siberia. Using these data, the study tests four current hypotheses regarding the hierarchical distribution of neutral genetic variation in native South American populations: (1) the western region of South America harbors more variation than the eastern region of South America, (2) Central American and western South American populations cluster exclusively, (3) populations speaking the Chibchan-Paezan and Equatorial-Tucanoan language stock emerge as a group within an otherwise South American clade, (4) Chibchan-Paezan populations in Central America emerge together at the tips of the Chibchan-Paezan cluster. This study finds that hierarchical models with the best fit place Central American populations, and populations speaking the Chibchan-Paezan language stock, at a basal position or separated from the South American group, which is more consistent with a serial founder effect into South America than that previously described. Western (Andean) South America is found to harbor similar levels of variation as eastern (Equatorial-Tucanoan and Ge-Pano-Carib) South America, which is inconsistent with an initial west coast migration into South America. Moreover, in all relevant models, the estimates of genetic diversity within geographic regions suggest a major bottleneck or founder effect occurring within the North American subcontinent, before the peopling of Central and South America.

Pre-Columbian population dynamics in coastal southern Peru: A diachronic investigation of mtDNA patterns in the Palpa region by ancient DNA analysis

Alternative models have been proposed to explain the formation and decline of the south Peruvian Nasca culture, ranging from migration or invasion to autochthonous development and ecological crisis. To reveal to what extent population dynamic processes accounted for cultural development in the Nasca mainland, or were influenced by them, we analyzed ancient mitochondrial DNA of 218 individuals, originating from chronologically successive archaeological sites in the Palpa region, the Paracas Peninsula, and the Andean highlands in southern Peru. The sampling strategy allowed a diachronic analysis in a time frame from approximately 800 BC to 800 AD. Mitochondrial coding region polymorphisms were successfully analyzed and replicated for 130 individuals and control region sequences (np 16021-16408) for 104 individuals to determine Native American mitochondrial DNA haplogroups and haplotypes. The results were compared with ancient and contemporary Peruvian populations to reveal genetic relations of the archaeological samples. Frequency data and statistics show clear proximity of the Nasca populations to the populations of the preceding Paracas culture from Palpa and the Peninsula, and suggest, along with archaeological data, that the Nasca culture developed autochthonously in the Rio Grande drainage. Furthermore, the influence of changes in socioeconomic complexity in the Palpa area on the genetic diversity of the local population could be observed. In all, a strong genetic affinity between pre-Columbian coastal populations from southern Peru could be determined, together with a significant differentiation from ancient highland and all present-day Peruvian reference populations, best shown in the differential distribution of mitochondrial haplogroups.

The late Pleistocene colonization of South America: an interdisciplinary perspective

In this article, we briefly review scenarios for the first peopling of the New World, including the Clovis First/Single Origin, Tripartite, and Dual Migration models. We then discuss bioanthropological, molecular genetic, archaeological and biogeographic evidence for the peopling of the continent. To conclude, we draw on different lines of evidence to make inferences concerning the timing, direction, and type of early human dispersion in South America.

Genetic and linguistic coevolution in Northern Island Melanesia

Recent studies have detailed a remarkable degree of genetic and linguistic diversity in Northern Island Melanesia. Here we utilize that diversity to examine two models of genetic and linguistic coevolution. The first model predicts that genetic and linguistic correspondences formed following population splits and isolation at the time of early range expansions into the region. The second is analogous to the genetic model of isolation by distance, and it predicts that genetic and linguistic correspondences formed through continuing genetic and linguistic exchange between neighboring populations. We tested the predictions of the two models by comparing observed and simulated patterns of genetic variation, genetic and linguistic trees, and matrices of genetic, linguistic, and geographic distances. The data consist of 751 autosomal microsatellites and 108 structural linguistic features collected from 33 Northern Island Melanesian populations. The results of the tests indicate that linguistic and genetic exchange have erased any evidence of a splitting and isolation process that might have occurred early in the settlement history of the region. The correlation patterns are also inconsistent with the predictions of the isolation by distance coevolutionary process in the larger Northern Island Melanesian region, but there is strong evidence for the process in the rugged interior of the largest island in the region (New Britain). There we found some of the strongest recorded correlations between genetic, linguistic, and geographic distances. We also found that, throughout the region, linguistic features have generally been less likely to diffuse across population boundaries than genes. The results from our study, based on exceptionally fine-grained data, show that local genetic and linguistic exchange are likely to obscure evidence of the early history of a region, and that language barriers do not particularly hinder genetic exchange. In contrast, global patterns may emphasize more ancient demographic events, including population splits associated with the early colonization of major world regions.

Mixed integer linear programming for maximum-parsimony phylogeny inference

Reconstruction of phylogenetic trees is a fundamental problem in computational biology. While excellent heuristic methods are available for many variants of this problem, new advances in phylogeny inference will be required if we are to be able to continue to make effective use of the rapidly growing stores of variation data now being gathered. In this paper, we present two integer linear programming (ILP) formulations to find the most parsimonious phylogenetic tree from a set of binary variation data. One method uses a flow-based formulation that can produce exponential numbers of variables and constraints in the worst case. The method has, however, proven extremely efficient in practice on datasets that are well beyond the reach of the available provably efficient methods, solving several large mtDNA and Y-chromosome instances within a few seconds and giving provably optimal results in times competitive with fast heuristics than cannot guarantee optimality. An alternative formulation establishes that the problem can be solved with a polynomial-sized ILP. We further present a web server developed based on the exponential-sized ILP that performs fast maximum parsimony inferences and serves as a front end to a database of precomputed phylogenies spanning the human genome.

Genetic variation and population structure in native Americans

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians--signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.

Uniparental (mtDNA, Y-chromosome) polymorphisms in French Guiana and two related populations--implications for the region's colonization

Blood samples collected in four Amerindian French Guiana populations (Palikur, Emerillon, Wayampi and Kali'na) in the early 1980s were screened for selected mtDNA and Y-chromosome length polymorphisms, and sequenced for the mtDNA hypervariable segment I (HVS-I). In addition, two other Amerindian populations (Apalaí and Matsiguenga) were examined for the same markers to establish the genetic relationships in the area. Strong dissimilarities were observed in the distribution of the founding Amerindian haplogroups, and significant p-values were obtained from F(ST) genetic distances. Interpopulation similarities occurred mainly due to geography. The Palikur did not show obvious genetic similarity to the Matsiguenga, who speak the same language and live in a region from where they could have migrated to French Guiana. The African-origin admixture observed in the Kali'na probably derives from historical contacts they had with the Bushinengue (Noir Marron), a group of escaped slaves who now lead independent lives in a nearby region. This analysis has identified significant clues about the Amerindian peopling of the North-East Amazonian region.

Mitochondrial DNA diversity of the Amerindian populations living in the Andean Piedmont of Bolivia: Chimane, Moseten, Aymara and Quechua

BACKGROUND

Chimane, Moseten Aymara and Quechua are Amerindian populations living in the Bolivian Piedmont, a characteristic ecoregion between the eastern slope of the Andean mountains and the Amazonian Llanos de Moxos. In both neighbouring areas, dense and complex societies have developed over the centuries. The Piedmont area is especially interesting from a human peopling perspective since there is no clear evidence regarding the genetic influence and peculiarities of these populations. This land has been used extensively as a territory of economic and cultural exchange between the Andes and Amazonia, however Chimane and Moseten populations have been sufficiently isolated from their neighbour groups to be recognized as distinct populations. Genetic information suggests that evolutionary processes, such as genetic drift, natural selection and genetic admixture have formed the history of the Piedmont populations.

AIM

The objective of this study is to characterize the genetic diversity of the Piedmont populations, analysing the sequence variability of the HVR-I control region in the mitochondrial DNA (mtDNA). Haplogroup mtDNA data available from the whole of Central and South America were utilized to determine the relationship of the Piedmont populations with other Amerindian populations.

SUBJECT AND METHODS

Hair pulls were obtained in situ, and DNA from non-related individuals was extracted using a standard Chelex 100 method. A 401 bp DNA fragment of HVR-I region was amplified using standard procedures. Two independent 401 and 328 bp DNA fragments were sequenced separately for each sample. The sequence analyses included mismatch distribution and mean pairwise differences, median network analyses, AMOVA and principal component analyses. The genetic diversity of DNA sequences was measured and compared with other South Amerindian populations.

RESULTS

The genetic diversity of 401 nucleotide mtDNA sequences, in the hypervariable Control Region, from positions 16 000-16 400, was characterized in a sample of 46 Amerindians living in the Piedmont area in the Beni Department of Bolivia. The results obtained indicate that the genetic diversity in the area is higher than that observed in other American groups living in much larger areas and despite the reduced size of the studied area the human groups analysed show high levels of inter-group variability. In addition, results show that Amerindian populations living in the Piedmont are genetically more related to those in the Andean than in the Amazonian populations.

Biological relationship between Central and South American Chibchan speaking populations: evidence from mtDNA

We examined mitochondrial DNA (mtDNA) haplogroup and haplotype diversity in 188 individuals from three Chibchan (Kogi, Arsario, and Ijka) populations and one Arawak (Wayuú) group from northeast Colombia to determine the biological relationship between lower Central American and northern South American Chibchan speakers. mtDNA haplogroups were obtained for all individuals and mtDNA HVS-I sequence data were obtained for 110 samples. Resulting sequence data were compared to 16 other Caribbean, South, and Central American populations using diversity measures, neutrality test statistics, sudden and spatial mismatch models, intermatch distributions, phylogenetic networks, and a multidimensional scaling plot. Our results demonstrate the existence of a shared maternal genetic structure between Central American Chibchan, Mayan populations and northern South American Chibchan-speakers. Additionally, these results suggest an expansion of Chibchan-speakers into South America associated with a shift in subsistence strategies because of changing ecological conditions that occurred in the region between 10,000-14,000 years before present.

Melanesian mtDNA complexity

Melanesian populations are known for their diversity, but it has been hard to grasp the pattern of the variation or its underlying dynamic. Using 1,223 mitochondrial DNA (mtDNA) sequences from hypervariable regions 1 and 2 (HVR1 and HVR2) from 32 populations, we found the among-group variation is structured by island, island size, and also by language affiliation. The more isolated inland Papuan-speaking groups on the largest islands have the greatest distinctions, while shore dwelling populations are considerably less diverse (at the same time, within-group haplotype diversity is less in the most isolated groups). Persistent differences between shore and inland groups in effective population sizes and marital migration rates probably cause these differences. We also add 16 whole sequences to the Melanesian mtDNA phylogenies. We identify the likely origins of a number of the haplogroups and ancient branches in specific islands, point to some ancient mtDNA connections between Near Oceania and Australia, and show additional Holocene connections between Island Southeast Asia/Taiwan and Island Melanesia with branches of haplogroup E. Coalescence estimates based on synonymous transitions in the coding region suggest an initial settlement and expansion in the region at approximately 30-50,000 years before present (YBP), and a second important expansion from Island Southeast Asia/Taiwan during the interval approximately 3,500-8,000 YBP. However, there are some important variance components in molecular dating that have been overlooked, and the specific nature of ancestral (maternal) Austronesian influence in this region remains unresolved.

The β-globin gene cluster distribution revisited—Patterns in Native American populations

New frequencies for the beta-globin gene cluster haplotypes are presented for the Aché (N = 82 individuals), Guarani (N = 76), and Kaingang (N = 54), three Native South American populations that live in an area between parallels 20 degrees S and 30 degrees S not covered by previous studies at this locus. The haplotype frequencies obtained for the three populations are within the interval observed for 28 other Native American populations. The Aché show much less haplotypes (five) than the other two populations (9-10), the haplotype prevalences being more similar to those of the Guarani than to the Kaingang. The Native American total heterozygosity was about half (0.41) that obtained for the African populations (0.71), but was not much different from those obtained for other continents. A geographical pattern was disclosed in South America by mapping the frequencies of the most common haplotype (haplotype 2), and by means of spatial correlation analysis. The analysis of molecular variance (AMOVA) and pairwise F(ST) data suggest three distinct sectors for the genetic landscape of Native South America: the Andes, the Center/Southeast region, and the Amazon.

Mitochondrial DNA analysis of ancient Peruvian highlanders

Ancient DNA recovered from 57 individuals excavated by Hiram Bingham at the rural communities of Paucarcancha, Patallacta, and Huata near the famed Inca royal estate and ritual site of Machu Picchu was analyzed by polymerase chain reaction, and the results were compared with ancient and modern DNA from various Central Andean areas to test their hypothesized indigenous highland origins. The control and coding regions of the mitochondrial DNA (mtDNA) of 35 individuals in this group were sequenced, and the haplogroups of each individual were determined. The frequency data for the haplogroups of these samples show clear proximity to those of modern Quechua and Aymara populations in the Peruvian and Bolivian highlands, and contrast with those of pre-Hispanic individuals of the north coast of Peru that we defined previously. Our study suggests a strong genetic affinity between sampled late pre-Hispanic individuals and modern Andean highlanders. A previous analysis of the Machu Picchu osteological collection suggests that the residents there were a mixed group of natives from various coastal and highland regions relocated by the Inca state for varied purposes. Overall, our study indicates that the sampled individuals from Paucarcancha and Patallacta were indigenous highlanders who provided supportive roles for nearby Machu Picchu.

Is the genetic structure of Gran Chaco populations unique? Interregional perspectives on native South American mitochondrial DNA variation

This study reevaluates the hypothesis in Demarchi et al. (2001 Am. J. Phys. Anthropol. 115:199-203) that Gran Chaco peoples demonstrate a unique pattern of genetic diversity due to a distinct regional population history. Specifically, they found populations in the central part of the Gran Chaco, or Central Chaco, to have higher within- and lower between-population mitochondrial DNA (mtDNA) haplogroup frequency variation compared to populations in other South American regions. To test this hypothesis of regional uniqueness, we applied analytical and simulation methods to mtDNA first hypervariable (HVI) region sequence data from a broad set of comparative South and Central American population samples. Contrary to the results of Demarchi et al. (2001 Am. J. Phys. Anthropol. 115:199-203), we found that the Gran Chaco's regional within-population diversity is about average among regions, and populations are highly differentiated from each other. When we limited the scale of analysis to the Central Chaco, a more localized subregion of the Gran Chaco, our results fell more in line with the original findings of Demarchi et al. (2001 Am. J. Phys. Anthropol. 115:199-203). Still, we conclude that neither the Gran Chaco regional pattern, nor the Central Chaco subregional pattern, is unique within South America. Nonetheless, the Central Chaco pattern accords well with the area's history, including pre-European contact lifeways and the documented historical use of the area as an interregional crossroads. However, we cannot exclude post-European contact disruption of traditional mating networks as an equally plausible explanation for the observed diversity pattern. Finally, these results additionally inform broader models of South American genetic diversity. While other researchers proposed an east-west continental division in patterns of genetic variation (e.g., Fuselli et al. 2003 Mol. Biol. Evol. 20:1682-1691), we found that in the geographically intermediate Central Chaco, a strict east-west divide in genetic variation breaks down. We suggest that future genetic characterizations of the continent, and subsequent interpretations of evolutionary history, involve a broad regional sampling of South American populations.

Phantom mutation hotspots in human mitochondrial DNA

Phantom mutations are systematic artifacts generated in the course of the sequencing process. Contra common belief these artificial mutations are nearly ubiquitous in sequencing results, albeit at frequencies that may vary dramatically. The amount of artifacts depends not only on the sort of automated sequencer and sequencing chemistry employed, but also on other lab-specific factors. An experimental study executed on four samples under various combinations of sequencing conditions revealed a number of phantom mutations occurring at the same sites of mitochondrial DNA (mtDNA) repeatedly. To confirm these and identify further hotspots for artifacts, > 5000 mtDNA electropherograms were screened for artificial patterns. Further, > 30 000 published hypervariable segment I sequences were compared at potential hotspots for phantom mutations, especially for variation at positions 16085 and 16197. Resequencing of several samples confirmed the artificial nature of these and other polymorphisms in the original publications. Single-strand sequencing, as typically executed in medical and anthropological studies, is thus highly vulnerable to this kind of artifacts. In particular, phantom mutation hotspots could easily lead to misidentification of somatic mutations and to misinterpretations in all kinds of clinical mtDNA studies.