Arrival of Paleo-Indians to the southern cone of South America: new clues from mitogenomes

Novel B2 mitogenomes from Continental southern Patagonia's Late Holocene: New insights into the peopling of the Southern Cone

OBJECTIVES

The main aim of this study is to discuss the migratory processes and peopling dynamics that shaped the genetic variability of populations during the settlement of the Southern Cone, through the analysis of complete mitogenomes of individuals from southern Patagonia.

MATERIALS AND METHODS

Complete mitogenomes were sequenced through massively parallel sequencing from two late Holocene individuals (SAC 1-1-3 and SAC 1-1-4) buried in the same chenque at Salitroso Lake Basin (Santa Cruz province, Argentina). To evaluate matrilineal phylogenetic affinities with other haplotypes, maximum likelihood and Bayesian phylogenetic reconstructions were performed, as well as a haplotype median-joining network.

RESULTS

The mitogenomes were assigned to haplogroups B2 and B2b, exhibiting an average depth of 54X and 89X (≥1X coverage of 98.6% and 100%), and a high number of nucleotide differences among them. The phylogenetic analyses showed a relatively close relationship between the haplotype found in SAC 1-1-4 and those retrieved from a Middle Holocene individual from Laguna Chica (Buenos Aires province), and from a group of individuals from the Peruvian coast. For the SAC 1-1-3, no clear affiliations to any other haplotype were established.

DISCUSSION

The large divergence between the haplotypes presented in this study suggests either a highly variable founder gene pool, or a later enrichment by frequent biological contact with other populations. Our results underline the persistence of genetic signals related to the first waves of peopling in South America, suggesting that the regional settlement of the southern end of the continent has been much more complex than initially thought.

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 (14C) 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 (14C) 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.

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.

Ancient mitogenomes from the Southern Pampas of Argentina reflect local differentiation and limited extra-regional linkages after rapid initial colonization

OBJECTIVE

This study aims to contribute to the recovery of Indigenous evolutionary history in the Southern Pampas region of Argentina through an analysis of ancient complete mitochondrial genomes.

MATERIALS AND METHODS

We generated DNA data for nine complete mitogenomes from the Southern Pampas, dated to between 2531 and 723 cal BP. In combination with previously published ancient mitogenomes from the region and from throughout South America, we documented instances of extra-regional lineage-sharing, and estimated coalescent ages for local lineages using a Bayesian method with tip calibrations in a phylogenetic analysis.

RESULTS

We identified a novel mitochondrial haplogroup, B2b16, and two recently defined haplogroups, A2ay and B2ak1, as well as three local haplotypes within founder haplogroups C1b and C1d. We detected lineage-sharing with ancient and contemporary individuals from Central Argentina, but not with ancient or contemporary samples from North Patagonian or Littoral regions of Argentina, despite archeological evidence of cultural interactions with the latter regions. The estimated coalescent age of these shared lineages is ~10,000 years BP.

DISCUSSION

The history of the human populations in the Southern Pampas is temporally deep, exhibiting long-term continuity of mitogenome lineages. Additionally, the identification of highly localized mtDNA clades accords with a model of relatively rapid initial colonization of South America by Indigenous communities, followed by more local patterns of limited gene flow and genetic drift in various South American regions, including the Pampas.

Evaluating population histories in Patagonia and Tierra del Fuego, Chile, using ancient mitochondrial and Y-chromosomal DNA

OBJECTIVES

This study aims to characterize the genetic histories of ancient hunter-gatherer groups in Fuego-Patagonia (Chile) with distinct Marine, Terrestrial, and Mixed Economy subsistence strategies. Mitochondrial (mtDNA) and Y-chromosome data were generated to test three hypotheses. H₀: All individuals were drawn from the same panmictic population; H₁: Terrestrial groups first populated the region and gave rise to highly specialized Marine groups by ~7,500 cal BP; or H₂: Marine and Terrestrial groups represent distinct ancestral lineages who migrated independently into the region.

METHODS

Ancient DNA was extracted from the teeth of 50 Fuegian-Patagonian individuals dating from 6,895 cal BP to after European arrival, and analyzed alongside other individuals from previous studies. Individuals were assigned to Marine, Terrestrial, and Mixed Economy groups based on archeological context and stable isotope diet inferences, and mtDNA (HVR1/2) and Y-chromosome variation was analyzed.

RESULTS

Endogenous aDNA was obtained from 49/50 (98%) individuals. Haplotype diversities, F_ST comparisons, and exact tests of population differentiation showed that Marine groups were significantly different from Terrestrial groups based on mtDNA (p < 0.05). No statistically significant differences were found between Terrestrial and Mixed Economy groups. Demographic simulations support models in which Marine groups diverged from the others by ~14,000 cal BP. Y-chromosome results showed similar patterns but were not statistically significant due to small sample sizes and allelic dropout.

DISCUSSION

These results support the hypothesis that Marine and Terrestrial economic groups represent distinct ancestral lineages who diverged during the time populations were expanding in the Americas, and may represent independent migrations into Fuego-Patagonia.

Native American mitochondrial lineages in admixed populations from Chile: Detecting recent migrations during post-Columbian times using geographically restricted lineages

OBJECTIVES

To analyze the mitochondrial diversity in three admixed populations and evaluate the historical migration effect of native southern population movement to Santiago (capital of Chile). The intensity of migration was quantified using three mitochondrial lineages restricted to South-Central native groups.

METHODS

D-loop sequences were genotyped in 550 unrelated individuals from San Felipe-Los Andes (n = 108), Santiago (n = 217), and Concepción (n = 225). Sequence processing, alignment, and haplogroup inference were carried out, and different genetic structure analyses were performed for haplogroup frequencies and D-loop sequences.

RESULTS

The Native lineages B2i2, C1b13, and D1g were the most frequent haplogroups, especially in Santiago (71.8%). Despite the distance, this city showed a high-genetic affinity with southern populations, including Concepción (~500 km distant) and native groups, rather than with those from San Felipe-Los Andes (<100 km distant). In fact, there was a negative correlation between geographical and genetic distance among these cities (r corr = -0.5593, p value = 0.8387). Network analysis revealed shared haplotypes between Santiago, Concepción, and other southern populations. Finally, we found lineages from Concepción acting as ancestral nodes in the northern clade.

CONCLUSIONS

Considering the geographic distances from these cities, the results were not consistent with a model of genetic isolation by geographic distance, revealing the effects of a historical migration process from the south to the capital. We also show evidence of possible north-to-south migration during admixture onset in Concepción and in addition, we were able to identify previously unreported mitochondrial diversity in urban populations that became lost in Native groups post-European contact.

Uniparental origins of the admixed Argentine Patagonia

OBJECTIVES

We aimed to contribute to the understanding of the ancient geographic origins of the uniparentally inherited markers in modern admixed Argentinian populations from central Patagonia with new information provided for the city of Trelew. We attempted to highlight the importance of combining different genetic markers when studying population history.

METHODS

The mtDNA control region sequence was typified in 89 individuals and 12 Y-STR and 15 Y-SNP loci were analyzed in 66 males. With these data, analysis of molecular variance and Network analyses were carried out. We exhaustively compared the modern data with ancient mtDNA information. Finally, we tested the differences in continental origins estimated by uniparental and previously published biparental markers.

RESULTS

Native American mtDNAs (53.9%) increased when maternal ancestors were born in the northern (81.8%) and southern (58.5%) regions of Argentina or in Chile (77.8%). Population substructure was only observed for Y-chromosome haplotypes. Some mtDNA haplogroups have been present in the area for at least ca. 2762-2430 and ca. 500 (D1g and D1g4 haplogroups) and ca. 6736 and ca. 6620 (C1b and C1c haplogroups) years, respectively. In contrast, haplogroups B2i2 and C1b13, frequent in modern Patagonia populations, had not been found in previous ancient DNA studies of the region.

CONCLUSIONS

The results suggest that Native American ancestry is well preserved in the region. Trelew samples had characteristic native mtDNA haplogroups previously described in Chilean and Argentine Patagonian populations, but not observed in ancient samples until now. These findings support the idea that these lineages have a recent regional origin. Finally, the estimated proportions of continental ancestry depend on the genetic marker analyzed.

Ancient mitochondrial genomes from the Argentinian Pampas inform the early peopling of the Southern Cone of South America

The Southern Cone of South America (SCSA) is a key region for investigations about the peopling of the Americas. However, little is known about the eastern sector, the Argentinian Pampas. We analyzed 18 mitochondrial genomes-7 of which are novel-from human skeletal remains from 3 Early to Late Holocene archaeological sites. The Pampas present a distinctive genetic makeup compared to other Middle to Late Holocene pre-Columbian SCSA populations. We also report the earliest individuals carrying SCSA-specific mitochondrial haplogroups D1j and D1g from Early and Middle Holocene, respectively. Using these deep calibration time points in Bayesian phylogenetic reconstructions, we suggest that the first settlers of the Pampas were part of a single and rapid dispersal ∼15,600 years ago. Finally, we propose that present-day genetic differences between the Pampas and the rest of the SCSA are due to founder effects, genetic drift, and a partial population replacement ∼9,000 years ago.

Ancient and modern mitogenomes from Central Argentina: new insights into population continuity, temporal depth and migration in South America

The inverted triangle shape of South America places Argentina territory as a geographical crossroads between the two principal peopling streams that followed either the Pacific or the Atlantic coasts, which could have then merged in Central Argentina (CA). Although the genetic diversity from this region is therefore crucial to decipher past population movements in South America, its characterization has been overlooked so far. We report 92 modern and 22 ancient mitogenomes spanning a temporal range of 5000 years, which were compared with a large set of previously reported data. Leveraging this dataset representative of the mitochondrial diversity of the subcontinent, we investigate the maternal history of CA populations within a wider geographical context. We describe a large number of novel clades within the mitochondrial DNA tree, thus providing new phylogenetic interpretations for South America. We also identify several local clades of great temporal depth with continuity until the present time, which stem directly from the founder haplotypes, suggesting that they originated in the region and expanded from there. Moreover, the presence of lineages characteristic of other South American regions reveals the existence of gene flow to CA. Finally, we report some lineages with discontinuous distribution across the Americas, which suggest the persistence of relic lineages likely linked to the first population arrivals. The present study represents to date the most exhaustive attempt to elaborate a Native American genetic map from modern and ancient complete mitochondrial genomes in Argentina and provides relevant information about the general process of settlement in South America.

Distribution of maternal lineages in hunter-gatherer societies of the southern coast of Tierra del Fuego, Argentina

OBJECTIVES

The aim of this work is to explore the maternal genetic diversity of hunter-gatherers of the southern Tierra del Fuego, specifically the north coast of Beagle Channel, the Península Mitre, and Isla de los Estados through ancient mitochondrial DNA analysis.

MATERIALS AND METHODS

The hypervariable regions 1 and 2 of the mitochondrial genome of five individuals from the north coast of Beagle Channel, six individuals from Península Mitre, and one individual from Isla de los Estados were analyzed. Through diversity statistics, Analysis of Molecular Variance (AMOVA), and Median Joining networks analyses, maternal relationships in the region were evaluated and phylogenetic similarities between ancient and contemporary populations of Tierra del Fuego were determined.

RESULTS

The mitochondrial DNA lineages from the ancient individuals analyzed reveals the presence of subclades C1b and D1g. Pattern of decreasing genetic diversity toward the South is observed. The AMOVAs performed found no statistically significant differences between individuals of the north coast of Beagle Channel and Península Mitre-Isla de los Estados, and modern Yámana populations. Median joining network of haplotypes of clades C1 and D1g, show the same results.

DISCUSSION

Ethnohistoric and ethnographic records of Península Mitre show that this region was occupied during the 19th century by Haush or Manekenk populations, although their biological, cultural, and subsistence characterization is unclear. We explore their maternal lineages and encounter low levels of genetic diversity and the absence of population differentiation with modern Yámana groups. We suggest that Península Mitre-Isla de los Estado was part of the same hunting and gathering populations as those of the Beagle Channel.

First analysis of mitochondrial lineages from the eastern Pampa-Patagonia transition during the final late Holocene

OBJECTIVE

Studies on population genetics have become highly relevant for understanding the evolutionary history of human settlement in southern South America. The eastern Pampa-Patagonia transition is an area that stands out due to its complex population dynamics, especially during the last about 1,000 years BP. The aim of this work is to characterize the maternal lineages of individuals buried in the Paso Alsina 1 archaeological site (ca. 500 years BP) through the analysis of mitochondrial genetic variability, in order to discuss the population models previously proposed for the southern cone of South America.

METHODS

Mitochondrial HyperVariable Region I sequences were analyzed on teeth belonging to 20 adult individuals. Statistical analyses were carried out to compare the interpopulation and intrapopulation molecular variability between the results obtained in this work and those previously published data from pre-Hispanic human groups. D1 haplotype network was constructed drawing from data on ancient and extant population group samples.

RESULTS

Thirteen sequences (65%) were obtained from the 20 analyzed samples. The maternal lineages or subhaplogroups identified were D1g (69.24%), C1 (15.38%), D1 (7.69%), and D1j (7.69%). There was low haplotype variability within the site; some individuals could be matrilineally related.

DISCUSSION

The subhaplogroups registered in Paso Alsina 1 site are in accordance with those reported for ancient and contemporary Patagonian populations. The results suggest that an initial nucleus of individuals carrying mostly subhaplogroup D1g settled in northern Patagonia, from which local diversity of this matrilineage could have arisen. The existence of gene flow in the final late Holocene with groups from Northern Andean Patagonia, as well as from Central Argentina, is proposed. The D1j variant probably developed in the latter region.

Whole Genome Sequence, Variant Discovery and Annotation in Mapuche-Huilliche Native South Americans

Whole human genome sequencing initiatives help us understand population history and the basis of genetic diseases. Current data mostly focuses on Old World populations, and the information of the genomic structure of Native Americans, especially those from the Southern Cone is scant. Here we present annotation and variant discovery from high-quality complete genome sequences of a cohort of 11 Mapuche-Huilliche individuals (HUI) from Southern Chile. We found approximately 3.1 × 106 single nucleotide variants (SNVs) per individual and identified 403,383 (6.9%) of novel SNVs events. Analyses of large-scale genomic events detected 680 copy number variants (CNVs) and 4,514 structural variants (SVs), including 398 and 1,910 novel events, respectively. Global ancestry composition of HUI genomes revealed that the cohort represents a sample from a marginally admixed population from the Southern Cone, whose main genetic component derives from Native American ancestors. Additionally, we found that HUI genomes contain variants in genes associated with 5 of the 6 leading causes of noncommunicable diseases in Chile, which may have an impact on the risk of prevalent diseases in Chilean and Amerindian populations. Our data represents a useful resource that can contribute to population-based studies and for the design of early diagnostics or prevention tools for Native and admixed Latin American populations.

Analysis of the human Y-chromosome haplogroup Q characterizes ancient population movements in Eurasia and the Americas

BACKGROUND

Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America.

RESULTS

Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a "coastal and inland routes scenario" for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe.

CONCLUSIONS

We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Cultural Innovations Influence Patterns of Genetic Diversity in Northwestern Amazonia

Human populations often exhibit contrasting patterns of genetic diversity in the mtDNA and the nonrecombining portion of the Y-chromosome (NRY), which reflect sex-specific cultural behaviors and population histories. Here, we sequenced 2.3 Mb of the NRY from 284 individuals representing more than 30 Native American groups from Northwestern Amazonia (NWA) and compared these data to previously generated mtDNA genomes from the same groups, to investigate the impact of cultural practices on genetic diversity and gain new insights about NWA population history. Relevant cultural practices in NWA include postmarital residential rules and linguistic exogamy, a marital practice in which men are required to marry women speaking a different language. We identified 2,969 SNPs in the NRY sequences, only 925 of which were previously described. The NRY and mtDNA data showed different sex-specific demographic histories: female effective population size has been larger than that of males through time, which might reflect larger variance in male reproductive success. Both markers show an increase in lineage diversification beginning ∼5,000 years ago, which may reflect the intensification of agriculture, technological innovations, and the expansion of regional trade networks documented in the archaeological evidence. Furthermore, we find similar excesses of NRY versus mtDNA between-population divergence at both the local and continental scale, suggesting long-term stability of female versus male migration. We also find evidence of the impact of sociocultural practices on diversity patterns. Finally, our study highlights the importance of analyzing high-resolution mtDNA and NRY sequences to reconstruct demographic history, since this can differ considerably between sexes.

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. 

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.

Genomic insights into the origin and diversification of late maritime hunter-gatherers from the Chilean Patagonia

Patagonia was the last region of the Americas reached by humans who entered the continent from Siberia ∼15,000-20,000 y ago. Despite recent genomic approaches to reconstruct the continental evolutionary history, regional characterization of ancient and modern genomes remains understudied. Exploring the genomic diversity within Patagonia is not just a valuable strategy to gain a better understanding of the history and diversification of human populations in the southernmost tip of the Americas, but it would also improve the representation of Native American diversity in global databases of human variation. Here, we present genome data from four modern populations from Central Southern Chile and Patagonia (n = 61) and four ancient maritime individuals from Patagonia (∼1,000 y old). Both the modern and ancient individuals studied in this work have a greater genetic affinity with other modern Native Americans than to any non-American population, showing within South America a clear structure between major geographical regions. Native Patagonian Kawéskar and Yámana showed the highest genetic affinity with the ancient individuals, indicating genetic continuity in the region during the past 1,000 y before present, together with an important agreement between the ethnic affiliation and historical distribution of both groups. Lastly, the ancient maritime individuals were genetically equidistant to a ∼200-y-old terrestrial hunter-gatherer from Tierra del Fuego, which supports a model with an initial separation of a common ancestral group to both maritime populations from a terrestrial population, with a later diversification of the maritime groups.

The Paleo-Indian Entry into South America According to Mitogenomes

Recent and compelling archaeological evidence attests to human presence ∼14.5 ka at multiple sites in South America and a very early exploitation of extreme high-altitude Andean environments. Considering that, according to genetic evidence, human entry into North America from Beringia most likely occurred ∼16 ka, these archeological findings would imply an extremely rapid spread along the double continent. To shed light on this issue from a genetic perspective, we first completely sequenced 217 novel modern mitogenomes of Native American ancestry from the northwestern area of South America (Ecuador and Peru); we then evaluated them phylogenetically together with other available mitogenomes (430 samples, both modern and ancient) from the same geographic area and, finally, with all closely related mitogenomes from the entire double continent. We detected a large number (N = 48) of novel subhaplogroups, often branching into further subclades, belonging to two classes: those that arose in South America early after its peopling and those that instead originated in North or Central America and reached South America with the first settlers. Coalescence age estimates for these subhaplogroups provide time boundaries indicating that early Paleo-Indians probably moved from North America to the area corresponding to modern Ecuador and Peru over the short time frame of ∼1.5 ka comprised between 16.0 and 14.6 ka.

Enclaves of genetic diversity resisted Inca impacts on population history

The Inca Empire is claimed to have driven massive population movements in western South America, and to have spread Quechua, the most widely-spoken language family of the indigenous Americas. A test-case is the Chachapoyas region of northern Peru, reported as a focal point of Inca population displacements. Chachapoyas also spans the environmental, cultural and demographic divides between Amazonia and the Andes, and stands along the lowest-altitude corridor from the rainforest to the Pacific coast. Following a sampling strategy informed by linguistic data, we collected 119 samples, analysed for full mtDNA genomes and Y-chromosome STRs. We report a high indigenous component, which stands apart from the network of intense genetic exchange in the core central zone of Andean civilization, and is also distinct from neighbouring populations. This unique genetic profile challenges the routine assumption of large-scale population relocations by the Incas. Furthermore, speakers of Chachapoyas Quechua are found to share no particular genetic similarity or gene-flow with Quechua speakers elsewhere, suggesting that here the language spread primarily by cultural diffusion, not migration. Our results demonstrate how population genetics, when fully guided by the archaeological, historical and linguistic records, can inform multiple disciplines within anthropology.

[Mitochondrial DNA diversity in prehispanic bone remains on the eastern Colombian Andes]

INTRODUCTION

DNA extracted from ancient human bones allows to analyze the genetic makeup of pre-Columbian populations and to determine the dynamics that gave rise to the diversity of contemporary populations.

OBJECTIVE

To determine the genetic diversity of skeletal remains associated with the Templo del Sol (Sun Temple) and their relationship with other contemporary and ancient communities of America.

MATERIALS AND METHODS

We analyzed 13 individuals belonging to the pre-Columbian Muisca Period (IX-XVI centuries AD) from the vicinities of the Templo del Sol (Sun Temple) (Sogamoso, Boyacá) in the eastern Colombian Andes. Mitochondrial DNA was amplified and RFLPs were performed in order to type the four traditional Amerindian haplogroups (A, B, C and D). In addition, autosomal markers including amelogenin and Y-chromosome STRs were amplified.

RESULTS

Among the observed mitochondrial lineages, haplogroup A was the most frequent, followed by haplogroups B and C; no evidence of haplogroup D was found. The genetic variation analysis indicated a similar diversity of pre-ColumbianMuiscas to that of contemporary populations belonging to the Chibcha linguistic family from Colombia and Central America. Molecular sexing was accomplished and it was compared to osteological data. With only one exception, anthropological and molecular data were consistent.

CONCLUSIONS

Our results contribute new genetic elements supporting the hypothesis of Central American origin of the Chibcha groups of the Cundiboyacense plateau, and allowed sex typing and kinship evaluations.

Native American Genomics and Population Histories

Studies of Native American genetic diversity and population history have been transformed over the last decade by important developments in anthropological genetics. During this time, researchers have adopted new DNA technologies and computational approaches for analyzing genomic data, and they have become increasingly sensitive to the views of research participants and communities. As new methods are applied to long-standing questions, and as more research is conducted in collaboration with indigenous communities, we are gaining new insights into the history and diversity of indigenous populations. This review discusses the recent methodological advances and genetic studies that have improved our understanding of Native American genomics and population histories. We synthesize current knowledge about Native American genomic variation and build a model of population history in the Americas. We also discuss the broader implications of this research for anthropology and related disciplines, and we highlight challenges and other considerations for future research.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

A South American Prehistoric Mitogenome: Context, Continuity, and the Origin of Haplogroup C1d

Based on mitochondrial DNA (mtDNA), it has been estimated that at least 15 founder haplogroups peopled the Americas. Subhaplogroup C1d3 was defined based on the mitogenome of a living individual from Uruguay that carried a lineage previously identified in hypervariable region I sequences from ancient and modern Uruguayan individuals. When complete mitogenomes were studied, additional substitutions were found in the coding region of the mitochondrial genome. Using a complete ancient mitogenome and three modern mitogenomes, we aim to clarify the ancestral state of subhaplogroup C1d3 and to better understand the peopling of the region of the Río de la Plata basin, as well as of the builders of the mounds from which the ancient individuals were recovered. The ancient mitogenome, belonging to a female dated to 1,610±46 years before present, was identical to the mitogenome of one of the modern individuals. All individuals share the mutations defining subhaplogroup C1d3. We estimated an age of 8,974 (5,748–12,261) years for the most recent common ancestor of C1d3, in agreement with the initial peopling of the geographic region. No individuals belonging to the defined lineage were found outside of Uruguay, which raises questions regarding the mobility of the prehistoric inhabitants of the country. Moreover, the present study shows the continuity of Native lineages over at least 6,000 years.

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.

Demographic History of Indigenous Populations in Mesoamerica Based on mtDNA Sequence Data

The genetic characterization of Native American groups provides insights into their history and demographic events. We sequenced the mitochondrial D-loop region (control region) of 520 samples from eight Mexican indigenous groups. In addition to an analysis of the genetic diversity, structure and genetic relationship between 28 Native American populations, we applied Bayesian skyline methodology for a deeper insight into the history of Mesoamerica. AMOVA tests applying cultural, linguistic and geographic criteria were performed. MDS plots showed a central cluster of Oaxaca and Maya populations, whereas those from the North and West were located on the periphery. Demographic reconstruction indicates higher values of the effective number of breeding females (Nef) in Central Mesoamerica during the Preclassic period, whereas this pattern moves toward the Classic period for groups in the North and West. Conversely, Nef minimum values are distributed either in the Lithic period (i.e. founder effects) or in recent periods (i.e. population declines). The Mesomerican regions showed differences in population fluctuation as indicated by the maximum Inter-Generational Rate (IGRmax): i) Center-South from the lithic period until the Preclassic; ii) West from the beginning of the Preclassic period until early Classic; iii) North characterized by a wide range of temporal variation from the Lithic to the Preclassic. Our findings are consistent with the genetic variations observed between central, South and Southeast Mesoamerica and the North-West region that are related to differences in genetic drift, structure, and temporal survival strategies (agriculture versus hunter-gathering, respectively). Interestingly, although the European contact had a major negative demographic impact, we detect a previous decline in Mesoamerica that had begun a few hundred years before.

Epidemiological, clinical and biochemical characterization of the p.(Ala359Asp) SMPD1 variant causing Niemann-Pick disease type B

Niemann-Pick disease type B (NPDB) is a rare, inherited lysosomal storage disorder that occurs due to variants in the sphingomyelin phosphodiesterase 1 (SMPD1) gene and the resultant deficiency of acid sphingomyelinase (ASM) activity. While numerous variants causing NPDB have been described, only a small number have been studied in any detail. Herein, we describe the frequency of the p.(Ala359Asp) variant in the healthy Chilean population, and determine the haplotype background of homozygous patients to establish if this variant originated from a common founder. Genomic DNA samples from 1691 healthy individuals were analyzed for the p.(Ala359Asp) variant. The frequency of p.(Ala359Asp) was found to be 1/105.7, predicting a disease incidence of 1/44 960 in Chile, higher than the incidence estimated by the number of confirmed NPDB cases. We also describe the clinical characteristics of 13 patients homozygous for p.(Ala359Asp) and all of them had moderate to severe NPDB disease. In addition, a conserved haplotype and shared 280 Kb region around the SMPD1 gene was observed in the patients analyzed, indicating that the variant originated from a common ancestor. The haplotype frequency and mitochondrial DNA analysis suggest an Amerindian origin for the variant. To assess the effect of the p.(Ala359Asp) variant, we transfected cells with the ASM-p.(Ala359Asp) cDNA and the activity was only 4.2% compared with the wild-type cDNA, definitively demonstrating the causative effect of the variant on ASM function. Information on common variants such as p.(Ala359Asp) is essential to guide the successful implementation for future therapies and benefit to patients.

Human settlement history between Sunda and Sahul: a focus on East Timor (Timor-Leste) and the Pleistocenic mtDNA diversity

BACKGROUND

Distinct, partly competing, "waves" have been proposed to explain human migration in(to) today's Island Southeast Asia and Australia based on genetic (and other) evidence. The paucity of high quality and high resolution data has impeded insights so far. In this study, one of the first in a forensic environment, we used the Ion Torrent Personal Genome Machine (PGM) for generating complete mitogenome sequences via stand-alone massively parallel sequencing and describe a standard data validation practice.

RESULTS

In this first representative investigation on the mitochondrial DNA (mtDNA) variation of East Timor (Timor-Leste) population including >300 individuals, we put special emphasis on the reconstruction of the initial settlement, in particular on the previously poorly resolved haplogroup P1, an indigenous lineage of the Southwest Pacific region. Our results suggest a colonization of southern Sahul (Australia) >37 kya, limited subsequent exchange, and a parallel incubation of initial settlers in northern Sahul (New Guinea) followed by westward migrations <28 kya.

CONCLUSIONS

The temporal proximity and possible coincidence of these latter dispersals, which encompassed autochthonous haplogroups, with the postulated "later" events of (South) East Asian origin pinpoints a highly dynamic migratory phase.

The impact of modelling rate heterogeneity among sites on phylogenetic estimates of intraspecific evolutionary rates and timescales

Phylogenetic analyses of DNA sequence data can provide estimates of evolutionary rates and timescales. Nearly all phylogenetic methods rely on accurate models of nucleotide substitution. A key feature of molecular evolution is the heterogeneity of substitution rates among sites, which is often modelled using a discrete gamma distribution. A widely used derivative of this is the gamma-invariable mixture model, which assumes that a proportion of sites in the sequence are completely resistant to change, while substitution rates at the remaining sites are gamma-distributed. For data sampled at the intraspecific level, however, biological assumptions involved in the invariable-sites model are commonly violated. We examined the use of these models in analyses of five intraspecific data sets. We show that using 6-10 rate categories for the discrete gamma distribution of rates among sites is sufficient to provide a good approximation of the marginal likelihood. Increasing the number of gamma rate categories did not have a substantial effect on estimates of the substitution rate or coalescence time, unless rates varied strongly among sites in a non-gamma-distributed manner. The assumption of a proportion of invariable sites provided a better approximation of the asymptotic marginal likelihood when the number of gamma categories was small, but had minimal impact on estimates of rates and coalescence times. However, the estimated proportion of invariable sites was highly susceptible to changes in the number of gamma rate categories. The concurrent use of gamma and invariable-site models for intraspecific data is not biologically meaningful and has been challenged on statistical grounds; here we have found that the assumption of a proportion of invariable sites has no obvious impact on Bayesian estimates of rates and timescales from intraspecific data.

The first peopling of South America: new evidence from Y-chromosome haplogroup Q

Recent progress in the phylogenetic resolution of the Y-chromosome phylogeny permits the male demographic dynamics and migratory events that occurred in Central and Southern America after the initial human spread into the Americas to be investigated at the regional level. To delve further into this issue, we examined more than 400 Native American Y chromosomes (collected in the region ranging from Mexico to South America) belonging to haplogroup Q – virtually the only branch of the Y phylogeny observed in modern-day Amerindians of Central and South America – together with 27 from Mongolia and Kamchatka. Two main founding lineages, Q1a3a1a-M3 and Q1a3a1-L54(xM3), were detected along with novel sub-clades of younger age and more restricted geographic distributions. The first was also observed in Far East Asia while no Q1a3a1-L54(xM3) Y chromosome was found in Asia except the southern Siberian-specific sub-clade Q1a3a1c-L330. Our data not only confirm a southern Siberian origin of ancestral populations that gave rise to Paleo-Indians and the differentiation of both Native American Q founding lineages in Beringia, but support their concomitant arrival in Mesoamerica, where Mexico acted as recipient for the first wave of migration, followed by a rapid southward migration, along the Pacific coast, into the Andean region. Although Q1a3a1a-M3 and Q1a3a1-L54(xM3) display overlapping general distributions, they show different patterns of evolution in the Mexican plateau and the Andean area, which can be explained by local differentiations due to demographic events triggered by the introduction of agriculture and associated with the flourishing of the Great Empires.

Reconciling migration models to the Americas with the variation of North American native mitogenomes

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

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.