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.

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations.

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Continental Origin for Q Haplogroup Patrilineages in Argentina and Paraguay.

Haplogroup Q originated in Eurasia around 30,000 years ago. It is present in Y-chromosomes from Asia and Europe at rather low frequencies. Since America is undoubtedly one of the continents where this haplogroup is highly represented, it has been defined as one of the founding haplogroups. Its M3 clade has been early described as the most frequent, with pan-American representation. However, it was also possible to find several other haplogroup Q clades at low frequencies. Numerous mutations have been described for haplogroup Q, allowing analysis of its variability and assignment of its geographic origin. We have analyzed 442 samples of unrelated men from Argentina and Paraguay belonging to haplogroup Q; here we report specifically on 27 Q (xM3) lineages. We tested 3 single-nucleotide polymorphisms (SNPs) by amplified product-length polymorphism (APLP) analysis, 3 SNPs for restriction fragment length polymorphism (RFLP) analysis, 15 SNPs by Sanger sequencing, and 17 short tandem repeats (STRs). Our approach allowed us to identify five subhaplogroups. Q-M3 and Q-CTS2730/Z780 are undoubtedly autochthonous lineages and represent the most frequent subhaplogroups, with significant representation in self-defined aboriginal populations, and their autochthonous status has been previously described. The aim of present work was to identify the continental origin of the remaining Q lineages. Thus, we analyzed the STR haplotypes for the samples and compared them with haplotypes described by other authors for the rest of the world. Even when haplogroup Q lineages have been extensively studied in America, some of them could have their origin in post-Columbian human migration from Europe and Middle East.

Ancient human genomes suggest three ancestral populations for present-day Europeans.

We sequenced the genomes of a ∼7,000-year-old farmer from Germany and eight ∼8,000-year-old hunter-gatherers from Luxembourg and Sweden. We analysed these and other ancient genomes with 2,345 contemporary humans to show that most present-day Europeans derive from at least three highly differentiated populations: west European hunter-gatherers, who contributed ancestry to all Europeans but not to Near Easterners; ancient north Eurasians related to Upper Palaeolithic Siberians, who contributed to both Europeans and Near Easterners; and early European farmers, who were mainly of Near Eastern origin but also harboured west European hunter-gatherer related ancestry. We model these populations' deep relationships and show that early European farmers had ∼44% ancestry from a 'basal Eurasian' population that split before the diversification of other non-African lineages.

GBStools: A Statistical Method for Estimating Allelic Dropout in Reduced Representation Sequencing Data.

Reduced representation sequencing methods such as genotyping-by-sequencing (GBS) enable low-cost measurement of genetic variation without the need for a reference genome assembly. These methods are widely used in genetic mapping and population genetics studies, especially with non-model organisms. Variant calling error rates, however, are higher in GBS than in standard sequencing, in particular due to restriction site polymorphisms, and few computational tools exist that specifically model and correct these errors. We developed a statistical method to remove errors caused by restriction site polymorphisms, implemented in the software package GBStools. We evaluated it in several simulated data sets, varying in number of samples, mean coverage and population mutation rate, and in two empirical human data sets (N = 8 and N = 63 samples). In our simulations, GBStools improved genotype accuracy more than commonly used filters such as Hardy-Weinberg equilibrium p-values. GBStools is most effective at removing genotype errors in data sets over 100 samples when coverage is 40X or higher, and the improvement is most pronounced in species with high genomic diversity. We also demonstrate the utility of GBS and GBStools for human population genetic inference in Argentine populations and reveal widely varying individual ancestry proportions and an excess of singletons, consistent with recent population growth.

Genetic variation in populations from central Argentina based on mitochondrial and Y chromosome DNA evidence.

We present new data and analysis on the genetic variation of contemporary inhabitants of central Argentina, including a total of 812 unrelated individuals from 20 populations. Our goal was to bring new elements for understanding micro-evolutionary and historical processes that generated the genetic diversity of the region, using molecular markers of uniparental inheritance (mitochondrial DNA and Y chromosome). Almost 76% of the individuals show mitochondrial lineages of American origin. The Native American haplogroups predominate in all surveyed localities, except in one. The larger presence of Eurasian maternal lineages were observed in the plains (Pampas) of the southeast, whereas the African lineages are more frequent in northern Córdoba. On the other hand, the analysis of 258 male samples reveals that 92% of them present Eurasian paternal lineages, 7% carry Native American haplogroups, and only 1% of the males show African lineages. The maternal lineages have high genetic diversity homogeneously distributed throughout central Argentina, probably as result of a recent common origin and sustained gene flow. Migratory events that occurred in colonial and recent times should have contributed to hiding any traces of differentiation that might have existed in the past. The analysis of paternal lineages showed also homogeneous distribution of the variation together with a drastic reduction of the native male population.

Reconstructing Native American population history.

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved. One contentious issue is whether the settlement occurred by means of a single migration or multiple streams of migration from Siberia. The pattern of dispersals within the Americas is also poorly understood. To address these questions at a higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. Here we show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call 'First American'. However, speakers of Eskimo-Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan speakers on both sides of the Panama isthmus, who have ancestry from both North and South America.

The mitochondrial DNA history of a former native American village in northern Uruguay.

OBJECTIVES: In 1828, between 8,000 and 15,000 Indians from the Jesuit Missions were brought to Uruguay. There, they were settled in a village, presently named Bella Unión, in the northwest corner of the country. According to historic sources, the Indians abandoned the settlement shortly thereafter, with the village subsequently repopulated by "criollos" and immigrants from abroad. As a first approach to reconstruct the genetic history of the population, data about the living population genetic structure will be used. Based on the analysis of the maternal lineages of the inhabitants of Bella Unión, and of those from two nearby villages, we expect to partially answer what happened with the first and subsequent inhabitants. METHODS: We analyzed the maternal lineages of the present inhabitants of Bella Unión and neighboring localities through the sequencing of the mitochondrial DNA control region. RESULTS: A total of 64.3%, 5.7%, and 30% of the mtDNAs were of Native, African, and West Eurasian origin, respectively. These figures are quite similar to that of the population of Tacuarembó, which is located in northeastern Uruguay. The four main Native American founding haplogroups were detected, with B2 being the most frequent, while some rare subhaplogroups (B2h, C1b2, D1f1) were also found. When compared with other Native American sequences, near- matches most consistently pointed to an Amazonian Indian origin which, when considered with historical evidence, suggested a probable Guaraní-Missionary-related origin. CONCLUSIONS: The data support the existence of a relationship between the historic and present inhabitants of the extreme northwest Uruguay, with a strong contribution of Native Americans to the mitochondrial DNA diversity observed there.

Maternal admixture and population structure in Mexican-Mestizos based on mtDNA haplogroups.

The maternal ancestry (mtDNA) has important applications in different research fields, such as evolution, epidemiology, identification, and human population history. This is particularly interesting in Mestizos, which constitute the main population in Mexico (∼93%) resulting from post-Columbian admixture between Spaniards, Amerindians, and African slaves, principally. Consequently, we conducted minisequencing analysis (SNaPshot) of 11 mitochondrial single-nucleotide polymorphisms in 742 Mestizos of 10 populations from different regions in Mexico. The predominant maternal ancestry was Native American (92.9%), including Haplogroups A, B, C, and D (47, 23.7, 15.9, and 6.2%, respectively). Conversely, European and African ancestries were less frequent (5.3 and 1.9%, respectively). The main characteristics of the maternal lineages observed in Mexican-Mestizos comprised the following: 1) contrasting geographic gradient of Haplogroups A and C; 2) increase of European lineages toward the Northwest; 3) low or absent, but homogeneous, African ancestry throughout the Mexican territory; 4) maternal lineages in Mestizos roughly represent the genetic makeup of the surrounding Amerindian groups, particularly toward the Southeast, but not in the North and West; 5) continuity over time of the geographic distribution of Amerindian lineages in Mayas; and 6) low but significant maternal population structure (FST = 2.8%; P = 0.0000). The average ancestry obtained from uniparental systems (mtDNA and Y-chromosome) in Mexican-Mestizos was correlated with previous ancestry estimates based on autosomal systems (genome-wide single-nucleotide polymorphisms and short tandem repeats). Finally, the comparison of paternal and maternal lineages provided additional information concerning the gender bias admixture, mating patterns, and population structure in Mestizos throughout the Mexican territory.

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.

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.

Phylogeography of mitochondrial haplogroup D1: an early spread of subhaplogroup D1j from Central Argentina.

We analyzed the patterns of variation of haplogroup D1 in central Argentina, including new data and published information from other populations of South America. Almost 28% (107/388) of the individuals sampled in the region belong to haplogroup D1, whereas more than 52% of them correspond to the recently described subhaplogroup D1j (Bodner et al.: Genome Res 22 (2012) 811-820), defined by the presence of additional transitions at np T152C-C16242T-T16311C to the nodal D1 motif. This lineage was found at high frequencies across a wide territory with marked geographical-ecological differences. Additionally, 12 individuals present the mutation C16187T that defines the recently named subhaplogroup D1g (Bodner et al.: Genome Res 22 (2012) 811-820), previously described in populations of Patagonia and Tierra del Fuego. Based on our results and additional data already published, we postulate that the most likely origin of subhaplogroup D1j is the region of Sierras Pampeanas, which occupies the center and part of the northwestern portion of Argentina. The extensive yet restricted geographical distribution, the relatively large internal diversity, and the absence or low incidence of D1j in other regions of South America suggest the existence of an ancient metapopulation covering the Sierras Pampeanas, being this lineage its genetic signature. Further support for a scenario of local origin for D1j in the Sierras Pampeanas stems from the fact that early derivatives from a putative ancestral lineage carrying the transitions T16311C-T152C have only been found in this region, supporting the hypothesis that it might represent an ancestral motif previous to the appearance of D1j-specific change C16242T.

Human Y chromosome sequences from Q Haplogroup reveal a South American settlement pre-18,000 years ago and a profound genomic impact during the Younger Dryas.

The settlement of the Americas has been the focus of incessant debate for more than 100 years, and open questions regarding the timing and spatial patterns of colonization still remain today. Phylogenetic studies with complete human Y chromosome sequences are used as a highly informative tool to investigate the history of human populations in a given time frame. To study the phylogenetic relationships of Native American lineages and infer the settlement history of the Americas, we analyzed Y chromosome Q Haplogroup, which is a Pan-American haplogroup and represents practically all Native American lineages in Mesoamerica and South America. We built a phylogenetic tree for Q Haplogroup based on 102 whole Y chromosome sequences, of which 13 new Argentine sequences were provided by our group. Moreover, 1,072 new single nucleotide polymorphisms (SNPs) that contribute to its resolution and diversity were identified. Q-M848 is known to be the most frequent autochthonous sub-haplogroup of the Americas. The present is the first genomic study of Q Haplogroup in which current knowledge on Q-M848 sub-lineages is contrasted with the historical, archaeological and linguistic data available. The divergence times, spatial structure and the SNPs found here as novel for Q-Z780, a less frequent sub-haplogroup autochthonous of the Americas, provide genetic support for a South American settlement before 18,000 years ago. We analyzed how environmental events that occurred during the Younger Dryas period may have affected Native American lineages, and found that this event may have caused a substantial loss of lineages. This could explain the current low frequency of Q-Z780 (also perhaps of Q-F4674, a third possible sub-haplogroup autochthonous of the Americas). These environmental events could have acted as a driving force for expansion and diversification of the Q-M848 sub-lineages, which show a spatial structure that developed during the Younger Dryas period.

mtDNA and Y-chromosome diversity in Aymaras and Quechuas from Bolivia: different stories and special genetic traits of the Andean Altiplano populations.

Two Bolivian samples belonging to the two main Andean linguistic groups (Aymaras and Quechuas) were studied for mtDNA and Y-chromosome uniparental markers to evaluate sex-specific differences and give new insights into the demographic processes of the Andean region. mtDNA-coding polymorphisms, HVI-HVII control regions, 17 Y-STRs, and three SNPs were typed in two well-defined populations with adequate size samples. The two Bolivian samples showed more genetic differences for the mtDNA than for the Y-chromosome. For the mtDNA, 81% of Aymaras and 61% of Quechuas presented haplogroup B2. Native American Y-chromosomes were found in 97% of Aymaras (89% hg Q1a3a and 11% hg Q1a3*) and 78% of Quechuas (100% hg Q1a3a). Our data revealed high diversity values in the two populations, in agreement with other Andean studies. The comparisons with the available literature for both sets of markers indicated that the central Andean area is relatively homogeneous. For mtDNA, the Aymaras seemed to have been more isolated throughout time, maintaining their genetic characteristics, while the Quechuas have been more permeable to the incorporation of female foreigners and Peruvian influences. On the other hand, male mobility would have been widespread across the Andean region according to the homogeneity found in the area. Particular genetic characteristics presented by both samples support a past common origin of the Altiplano populations in the ancient Aymara territory, with independent, although related histories, with Peruvian (Quechuas) populations.

Paleogenetical study of pre-Columbian samples from Pampa Grande (Salta, Argentina).

Ancient DNA recovered from 21 individuals excavated from burial sites in the Pampa Grande (PG) region (Salta province) of North-Western Argentina (NWA) was analyzed using various genetic markers (mitochondrial DNA, autosomal STRs, and Y chromosomal STRs). The results were compared to ancient and modern DNA from various populations in the Andean and North Argentinean regions, with the aim of establishing their relationships with PG. The mitochondrial haplogroup frequencies described (11% A, 47% B, and 42% D) presented values comparable to those found for the ancient Andean populations from Peru and San Pedro de Atacama. On the other hand, mitochondrial and Y chromosomal haplotypes were specific to PG, as they did not match any other of the South American populations studied. The described genetic diversity indicates homogeneity in the genetic structure of the ancient Andean populations, which was probably facilitated by the intense exchange network in the Andean zone, in particular among Tiwanaku, San Pedro de Atacama, and NWA. The discovery of haplotypes unique to PG could be due to a loss of genetic diversity caused by recent events affecting the autochthonous populations (establishment of the Inca Empire in the region, colonization by the Europeans).

Fine-scale genomic analyses of admixed individuals reveal unrecognized genetic ancestry components in Argentina.

Similarly to other populations across the Americas, Argentinean populations trace back their genetic ancestry into African, European and Native American ancestors, reflecting a complex demographic history with multiple migration and admixture events in pre- and post-colonial times. However, little is known about the sub-continental origins of these three main ancestries. We present new high-throughput genotyping data for 87 admixed individuals across Argentina. This data was combined to previously published data for admixed individuals in the region and then compared to different reference panels specifically built to perform population structure analyses at a sub-continental level. Concerning the Native American ancestry, we could identify four Native American components segregating in modern Argentinean populations. Three of them are also found in modern South American populations and are specifically represented in Central Andes, Central Chile/Patagonia, and Subtropical and Tropical Forests geographic areas. The fourth component might be specific to the Central Western region of Argentina, and it is not well represented in any genomic data from the literature. As for the European and African ancestries, we confirmed previous results about origins from Southern Europe, Western and Central Western Africa, and we provide evidences for the presence of Northern European and Eastern African ancestries.

Median network analysis of defectively sequenced entire mitochondrial genomes from early and contemporary disease studies.

Sequence analysis of the mitochondrial genome has become a routine method in the study of mitochondrial diseases. Quite often, the sequencing efforts in the search of pathogenic or disease-associated mutations are affected by technical and interpretive problems, caused by sample mix-up, contamination, biochemical problems, incomplete sequencing, misdocumentation and insufficient reference to previously published data. To assess data quality in case studies of mitochondrial diseases, it is recommended to compare any mtDNA sequence under consideration to their phylogenetically closest lineages available in the Web. The median network method has proven useful for visualizing potential problems with the data. We contrast some early reports of complete mtDNA sequences to more recent total mtDNA sequencing efforts in studies of various mitochondrial diseases. We conclude that the quality of complete mtDNA sequences generated in the medical field in the past few years is somewhat unsatisfactory and may even fall behind that of pioneer manual sequencing in the early nineties. Our study provides a paradigm for an a posteriori evaluation of sequence quality and for detection of potential problems with inferring a pathogenic status of a particular mutation.

Brief communication: Restricted geographic distribution for Y-Q* paragroup in South America.

We analyzed 21 paragroup Q* Y chromosomes from South American aboriginal and urban populations. Our aims were to evaluate the phylogenetic status, geographic distribution, and genetic diversity in these groups of chromosomes and compare the degree of genetic variation in relation to Q1a3a haplotypes. All Q* chromosomes from our series and five samples from North American Q* presented the derivate state for M346, that is present upstream to M3, and determined Q1a3* paragroup. We found a restrictive geographic distribution and low frequency of Q1a3* in South America. We assumed that this low frequency could be reflecting extreme drift effects. However, several estimates of gene diversity do not support the existence of a severe bottleneck. The mean haplotype diversity expected was similar to that for South American Q1a3* and Q1a3a (0.478 and 0.501, respectively). The analysis of previous reports from other research groups and this study shows the highest frequencies of Q* for the West Corner and the Grand Chaco regions of South America. At present, there is no information on whether the phylogenetic status of Q* paragoup described in previous reports is similar to that of Q1a3* paragroup though our results support this possibility.

Population structure in Argentina.

We analyzed 391 samples from 12 Argentinian populations from the Center-West, East and North-West regions with the Illumina Human Exome Beadchip v1.0 (HumanExome-12v1-A). We did Principal Components analysis to infer patterns of populational divergence and migrations. We identified proportions and patterns of European, African and Native American ancestry and found a correlation between distance to Buenos Aires and proportion of Native American ancestry, where the highest proportion corresponds to the Northernmost populations, which is also the furthest from the Argentinian capital. Most of the European sources are from a South European origin, matching historical records, and we see two different Native American components, one that spreads all over Argentina and another specifically Andean. The highest percentages of African ancestry were in the Center West of Argentina, where the old trade routes took the slaves from Buenos Aires to Chile and Peru. Subcontinentaly, sources of this African component are represented by both West Africa and groups influenced by the Bantu expansion, the second slightly higher than the first, unlike North America and the Caribbean, where the main source is West Africa. This is reasonable, considering that a large proportion of the ships arriving at the Southern Hemisphere came from Mozambique, Loango and Angola.