Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela

Bioanthropological analysis of human remains from the archaic and classic period discovered in Puyil cave, Mexico

OBJECTIVES

Determine the geographic place of origin and maternal lineage of prehistoric human skeletal remains discovered in Puyil Cave, Tabasco State, Mexico, located in a region currently populated by Olmec, Zoque and Maya populations.

MATERIALS AND METHODS

All specimens were radiocarbon (14 C) dated (beta analytic), had dental modifications classified, and had an analysis of 13 homologous reference points conducted to evaluate artificial cranial deformation (ACD). Following DNA purification, hypervariable region I (HVR-1) of the mitogenome was amplified and Sanger sequenced. Finally, Next Generation Sequencing (NGS) was performed for total DNA. Mitochondrial DNA (mtDNA) variants and haplogroups were determined using BioEdit 7.2 and IGV software and confirmed with MITOMASTER and WebHome softwares.

RESULTS

Radiocarbon dating (14 C) demonstrated that the inhabitants of Puyil Cave lived during the Archaic and Classic Periods and displayed tabular oblique and tabular mimetic ACD. These pre-Hispanic remains exhibited five mtDNA lineages: A, A2, C1, C1c and D4. Network analysis revealed a close genetic affinity between pre-Hispanic Puyil Cave inhabitants and contemporary Maya subpopulations from Mexico and Guatemala, as well as individuals from Bolivia, Brazil, the Dominican Republic, and China.

CONCLUSIONS

Our results elucidate the dispersal of pre-Hispanic Olmec and Maya ancestors and suggest that ACD practices are closely related to Olmec and Maya practices. Additionally, we conclude that ACD has likely been practiced in the region since the Middle-Archaic Period.

Genetic Diversity in the Lesser Antilles and Its Implications for the Settlement of the Caribbean Basin

Historical discourses about the Caribbean often chronicle West African and European influence to the general neglect of indigenous people's contributions to the contemporary region. Consequently, demographic histories of Caribbean people prior to and after European contact are not well understood. Although archeological evidence suggests that the Lesser Antilles were populated in a series of northward and eastern migratory waves, many questions remain regarding the relationship of the Caribbean migrants to other indigenous people of South and Central America and changes to the demography of indigenous communities post-European contact. To explore these issues, we analyzed mitochondrial DNA and Y-chromosome diversity in 12 unrelated individuals from the First Peoples Community in Arima, Trinidad, and 43 unrelated Garifuna individuals residing in St. Vincent. In this community-sanctioned research, we detected maternal indigenous ancestry in 42% of the participants, with the remainder having haplotypes indicative of African and South Asian maternal ancestry. Analysis of Y-chromosome variation revealed paternal indigenous American ancestry indicated by the presence of haplogroup Q-M3 in 28% of the male participants from both communities, with the remainder possessing either African or European haplogroups. This finding is the first report of indigenous American paternal ancestry among indigenous populations in this region of the Caribbean. Overall, this study illustrates the role of the region's first peoples in shaping the genetic diversity seen in contemporary Caribbean populations.

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.

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.

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.

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.

Recovering mitochondrial DNA lineages of extinct Amerindian nations in extant homopatric Brazilian populations

Background

Brazilian Amerindians have experienced a drastic population decrease in the past 500 years. Indeed, many native groups from eastern Brazil have vanished. However, their mitochondrial mtDNA haplotypes, still persist in Brazilians, at least 50 million of whom carry Amerindian mitochondrial lineages. Our objective was to test whether, by analyzing extant rural populations from regions anciently occupied by specific Amerindian groups, we could identify potentially authentic mitochondrial lineages, a strategy we have named 'homopatric targeting'.

Results

We studied 173 individuals from Queixadinha, a small village located in a territory previously occupied by the now extinct Botocudo Amerindian nation. Pedigree analysis revealed 74 unrelated matrilineages, which were screened for Amerindian mtDNA lineages by restriction fragment length polymorphism. A cosmopolitan control group was composed of 100 individuals from surrounding cities. All Amerindian lineages identified had their hypervariable segment HVSI sequenced, yielding 13 Amerindian haplotypes in Queixadinha, nine of which were not present in available databanks or in the literature. Among these haplotypes, there was a significant excess of haplogroup C (70%) and absence of haplogroup A lineages, which were the most common in the control group. The novelty of the haplotypes and the excess of the C haplogroup suggested that we might indeed have identified Botocudo lineages. To validate our strategy, we studied teeth extracted from 14 ancient skulls of Botocudo Amerindians from the collection of the National Museum of Rio de Janeiro. We recovered mtDNA sequences from all the teeth, identifying only six different haplotypes (a low haplotypic diversity of 0.8352 ± 0.0617), one of which was present among the lineages observed in the extant individuals studied.

Conclusions

These findings validate the technique of homopatric targeting as a useful new strategy to study the peopling and colonization of the New World, especially when direct analysis of genetic material is not possible.

Mitochondrial DNA diversity in a population from Santa Catarina (Brazil): predominance of the European input

The state of Santa Catarina (Brazil) is known to have represented a cultural crossroads in South America due to several historic migrations mainly from Europe and Africa. We set out to scrutinize whether the genetic imprint of these migrations could be traced through analysis of the matrilineal gene pool of the Catarinenses. The entire control region of the mitochondrial DNA was studied in 80 healthy and maternally unrelated individuals. The analysis of haplogroup distribution revealed that this population is extremely heterogeneous, showing the coexistence of matrilineal lineages with three different phylogeographic origins. European lineages are the most frequent due mainly to the impact of relatively recent migratory waves from Europe. In spite of this, Native American lineages and African lineages incorporated with the slave trade are also present in noticeable proportions. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9930) and power of discrimination (0.9806). Thus, analysis of the entire mitochondrial DNA control region emerges as a valuable tool for forensic genetic purposes in this highly admixed population, an attribute common to several present-day Latin American populations.

Mitochondrial echoes of first settlement and genetic continuity in El Salvador

Background

From Paleo-Indian times to recent historical episodes, the Mesoamerican isthmus played an important role in the distribution and patterns of variability all around the double American continent. However, the amount of genetic information currently available on Central American continental populations is very scarce. In order to shed light on the role of Mesoamerica in the peopling of the New World, the present study focuses on the analysis of the mtDNA variation in a population sample from El Salvador.

Methodology/Principal Findings

We have carried out DNA sequencing of the entire control region of the mitochondrial DNA (mtDNA) genome in 90 individuals from El Salvador. We have also compiled more than 3,985 control region profiles from the public domain and the literature in order to carry out inter-population comparisons. The results reveal a predominant Native American component in this region: by far, the most prevalent mtDNA haplogroup in this country (at ∼90%) is A2, in contrast with other North, Meso- and South American populations. Haplogroup A2 shows a star-like phylogeny and is very diverse with a substantial proportion of mtDNAs (45%; sequence range 16090–16365) still unobserved in other American populations. Two different Bayesian approaches used to estimate admixture proportions in El Salvador shows that the majority of the mtDNAs observed come from North America. A preliminary founder analysis indicates that the settlement of El Salvador occurred about 13,400±5,200 Y.B.P.. The founder age of A2 in El Salvador is close to the overall age of A2 in America, which suggests that the colonization of this region occurred within a few thousand years of the initial expansion into the Americas.

Conclusions/Significance

As a whole, the results are compatible with the hypothesis that today's A2 variability in El Salvador represents to a large extent the indigenous component of the region. Concordant with this hypothesis is also the observation of a very limited contribution from European and African women (∼5%). This implies that the Atlantic slave trade had a very small demographic impact in El Salvador in contrast to its transformation of the gene pool in neighbouring populations from the Caribbean facade.

The mtDNA ancestry of admixed Colombian populations

A total of 185 individuals from Colombia were sequenced for the first hypervariable region (HVS-I) of the mitochondrial DNA (mtDNA) genome, and a subset of these individuals were additionally genotyped for the second hypervariable segment (HVS-II). These individuals were collected according to their "self-reported ethnicity" in Colombia, comprising "Mestizos," "Mulatos," and "Afro-Colombians." We used databases containing more than 4,300 Native American lineages, 6,800 Africans, and 15,600 Europeans for population comparisons and phylogeographic inferences. We observe that Mulatos and Afro-Colombians have a dominant African mtDNA component, whereas Mestizos carry predominantly Native American haplotypes. All the populations analyzed have high diversity indices and there are no signatures of dramatic genetic drift episodes. Central and South America are the main candidate source populations of the Colombian Native American lineages, whereas west-central, southwest, and southeast Africa are the main original mtDNA sources for the African Colombian mtDNAs. We found that our results differ from those obtained in other studies for the same "population groups" in terms of haplogroup frequencies. This observation leads us to conclude that (i) self-reported ancestry is not a reliable proxy to indicate an individual's "ethnicity" in Colombia, (ii) our results do not support the use of outmoded race descriptions (Mestizos, Mulatos, etc.) mainly because these labels do not correspond to any genetically homogeneous population group, and (iii) studies relying on these terms to describe the population group of the individual, which then treat them as genetically homogeneous, carry a high risk of type I error (false positives) in medical studies in this country and of misinterpretation of the frequency of observed variation in forensic casework.

Indentured migration and differential gender gene flow: the origin and evolution of the East-Indian community of Limón, Costa Rica

After the emancipation of African slaves in the Caribbean, the labor void left by out-migrating former slaves was filled by in-migrating indentured servants from prepartition India and China. In some areas of the Caribbean such as Trinidad, Suriname, and Guyana, the East-Indian migrants formed large communities. In this article, we report a study based on mtDNA and Y-chromosomal markers of a small East-Indian community from Limón, Costa Rica. The purpose of the project is to determine the place of origin in the Indian subcontinent of the ancestors of our group and the contributions to its gene pool through gene flow by members of other ethnic groups. Both Y-chromosome and mtDNA suggest that the Indo-Costa Ricans descend from migrants primarily from Central India. While both paternal and maternal markers indicate that this group is overwhelmingly of Indian origin, they also indicate that males and females of African, European, and Amerindian origin contributed to it differently. We discuss our results in the historical context of the virtual extinction of Amerindian Caribbean groups, the forced migration of African slaves to the Caribbean, and the gene flow between Amerindians, Europeans, East-Indians, and Africans that eventually produced the Caribbean's currently diverse gene pool.

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.

Mosaics of ancient mitochondrial DNA: positive indicators of nonauthenticity

Research into ancient mitochondrial DNA is plagued by contamination, post mortem damage, and other artefacts. The stringent set of controls suggested by Cooper and Poinar a few years ago are, however, rarely followed in practice, and even when applied carefully, these criteria need not be sufficient to guarantee authenticity. The fairly relaxed prerequisites now common for ancient population studies have opened the door for all kinds of contamination and sequencing errors to enter ancient mtDNA data. To reject or question authenticity of particular sequencing results a posteriori, one can follow similar strategies of focused database comparisons that have proven to be effective and successful in the case of flawed modern mtDNA data.