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

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.

"What We Know and What We Do Not Know about Evolutionary Genetic Adaptation to High Altitude Hypoxia in Andean Aymaras"

Three well-studied populations living at high altitudes are Tibetans, Andeans (Aymaras and Quechuas), and Ethiopians. Unlike Tibetans and Ethiopians who have similar hemoglobin (Hb) levels as individuals living at sea level, Aymara Hb levels increase when living at higher altitudes. Our previous whole genome study of Aymara people revealed several selected genes that are involved in cardiovascular functions, but their relationship with Hb levels was not elucidated. Here, we studied the frequencies of known evolutionary-selected variants in Tibetan and Aymara populations and their correlation with high Hb levels in Aymara. We genotyped 177 Aymaras at three different altitudes: 400 m (Santa Cruz), 4000 m (La Paz), and 5000 m (Chorolque), and correlated the results with the elevation of residence. Some of the Tibetan-selected variants also exist in Aymaras, but at a lower prevalence. Two of 10 Tibetan selected variants of EPAS1 were found (rs13005507 and rs142764723) and these variants did not correlate with Hb levels. Allele frequencies of 5 Aymara selected SNPs (heterozygous and homozygous) at 4000 m (rs11578671_BRINP3, rs34913965_NOS2, rs12448902_SH2B1, rs10744822_TBX5, and rs487105_PYGM) were higher compared to Europeans. The allelic frequencies of rs11578671_BRINP3, rs34913965_NOS2, and rs10744822_SH2B1 were significantly higher for Aymaras living at 5000 m than those at 400 m elevation. Variant rs11578671, close to the BRINP3 coding region, correlated with Hb levels in females. Variant rs34913965 (NOS2) correlated with leukocyte counts. Variants rs12448902 (SH2B1) and rs34913965 (NOS2) associated with higher platelet levels. The correlation of these SNPs with blood cell counts demonstrates that the selected genetic variants in Aymara influence hematopoiesis and cardiovascular effects.

An approach on the migratory processes in the north of Chile based on Y chromosome analysis

OBJECTIVES

Northern Chile is an area characterized by a complex cultural and demographic trajectory. During the last few centuries, this complex trajectory has become the destination of intra- and intercontinental migratory waves. In this study, we analyzed the Y chromosome to evaluate how migratory and admixture patterns have affected the genetic composition of the populations in northern Chile compared with other populations of the country.

METHODS

A total of 311 people from urban (Antofagasta and Calama), rural (Azapa and Camarones), and Native (Aymara and Atacameño) populations from northern Chile were characterized by 26 SNPs and the STR DYS393 of the Y chromosome, along with 69 individuals from Native populations (Mapuche, Pehuenche, and Huilliche) from southern Chile. In addition to characterizing the paternal lineages, multivariate analyses were performed to compare with published data from other Chilean populations.

RESULTS

Both the Antofagasta and Calama populations show differences compared with the rest of the Chilean population. On one side, Antofagasta shows a high diversity of non-Amerindian lineages, including the highest value for haplogroup I (12%) for all Chileans populations. Otherwise, Calama has the highest value of any Chilean urban population (31.9%) for Amerindian lineages, including the only Q-M3 sub-lineage detected in the entire sample. Regarding the Native population, Aymara presents the highest percentage of Q-M3 (94.4%).

CONCLUSIONS

The Y chromosome haplogroup distribution allowed us to identify recent migratory processes typical of the northern populations studied. These have shaped the demographic and cultural dynamics of local and migrant groups in the territory.

Uniparental Lineages from the Oldest Indigenous Population of Ecuador: The Tsachilas

Together with Cayapas, the Tsachilas constitute the oldest population in the country of Ecuador and, according to some historians, they are the last descendants of the ancient Yumbos. Several anthropological issues underlie the interest towards this peculiar population: the uncertainty of their origin, their belonging to the Barbacoan linguistic family, which is still at the center of an intense linguistic debate, and the relations of their Yumbo ancestors with the Inca invaders who occupied their ancient territory. Our contribution to the knowledge of their complex past was the reconstruction of their genetic maternal and paternal inheritance through the sequencing of 70 entire mitochondrial genomes and the characterization of the non-recombinant region of the Y chromosome in 26 males. For both markers, we built comprehensive datasets of various populations from the surrounding geographical area, northwestern South America, NW, with a known linguistic affiliation, and we could then compare our sample against the overall variability to infer relationships with other Barbacoan people and with other NW natives. We found contrasting patterns of genetic diversity for the two markers, but generally, our results indicated a possible common origin between the Tsachilas, the Chachi, and other Ecuadorian and Colombian Barbacoans and are suggestive of an interesting ancient linkage to the Inca invaders in Yumbo country.

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.

Genetic assessment reveals no population substructure and divergent regional and sex-specific histories in the Chachapoyas from northeast Peru

Many native populations in South America have been severely impacted by two relatively recent historical events, the Inca and the Spanish conquest. However decisive these disruptive events may have been, the populations and their gene pools have been shaped markedly also by the history prior to the conquests. This study focuses mainly on the Chachapoya peoples that inhabit the montane forests on the eastern slopes of the northern Peruvian Andes, but also includes three distinct neighboring populations (the Jívaro, the Huancas and the Cajamarca). By assessing mitochondrial, Y-chromosomal and autosomal diversity in the region, we explore questions that have emerged from archaeological and historical studies of the regional culture (s). These studies have shown, among others, that Chachapoyas was a crossroads for Coast-Andes-Amazon interactions since very early times. In this study, we examine the following questions: 1) was there pre-Hispanic genetic population substructure in the Chachapoyas sample? 2) did the Spanish conquest cause a more severe population decline on Chachapoyan males than on females? 3) can we detect different patterns of European gene flow in the Chachapoyas region? and, 4) did the demographic history in the Chachapoyas resemble the one from the Andean area? Despite cultural differences within the Chachapoyas region as shown by archaeological and ethnohistorical research, genetic markers show no significant evidence for past or current population substructure, although an Amazonian gene flow dynamic in the northern part of this territory is suggested. The data also indicates a bottleneck c. 25 generations ago that was more severe among males than females, as well as divergent population histories for populations in the Andean and Amazonian regions. In line with previous studies, we observe high genetic diversity in the Chachapoyas, despite the documented dramatic population declines. The diverse topography and great biodiversity of the northeastern Peruvian montane forests are potential contributing agents in shaping and maintaining the high genetic diversity in the Chachapoyas region.

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.

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.

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.

Mapping the genomic mosaic of two 'Afro-Bolivians' from the isolated Yungas valleys

Background

Unraveling the ancestry of ‘Afro-American’ communities is hampered by the complex demographic processes that took place during the Transatlantic Slave Trade (TAST) and the (post-)colonization periods. ‘Afro-Bolivians’ from the subtropical Yungas valleys constitute small and isolated communities that live surrounded by the predominant Native American community of Bolivia. By genotyping >580,000 SNPs in two ‘Afro-Bolivians’, and comparing these genomic profiles with data compiled from more than 57 African groups and other reference ancestral populations (n = 1,161 in total), we aimed to disentangle the complex admixture processes undergone by ‘Afro-Bolivians’.

Results

The data indicate that these two genomes constitute a complex mosaic of ancestries that is approximately 80 % of recent African origin; the remaining ~20 % being European and Native American. West-Central Africa contributed most of the African ancestry to ‘Afro-Bolivians’, and this component is related to populations living along the Atlantic coast (i.e. Senegal, Ghana, Nigeria). Using tract length distribution of genomic segments attributable to distinct ancestries, we could date the time of admixture in about 400 years ago. This time coincides with the maximum importation of slaves to Bolivia to compensate the diminishing indigenous labor force needed for the development of the National Mint of Potosí.

Conclusions

Overall, the data indicate that the genome of ‘Afro-Bolivians’ was shaped by a complex process of admixture occurring in America among individuals originating in different West-Central African populations; their genomic mosaics received additional contributions of Europeans and local Native Americans (e.g. Aymaras).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2520-x) contains supplementary material, which is available to authorized users.

Exploring the Y Chromosomal Ancestry of Modern Panamanians

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama’s population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

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

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

A novel candidate region for genetic adaptation to high altitude in Andean populations

Humans living at high altitude (≥2,500 meters above sea level) have acquired unique abilities to survive the associated extreme environmental conditions, including hypoxia, cold temperature, limited food availability and high levels of free radicals and oxidants. Long-term inhabitants of the most elevated regions of the world have undergone extensive physiological and/or genetic changes, particularly in the regulation of respiration and circulation, when compared to lowland populations. Genome scans have identified candidate genes involved in altitude adaption in the Tibetan Plateau and the Ethiopian highlands, in contrast to populations from the Andes, which have not been as intensively investigated. In the present study, we focused on three indigenous populations from Bolivia: two groups of Andean natives, Aymara and Quechua, and the low-altitude control group of Guarani from the Gran Chaco lowlands. Using pooled samples, we identified a number of SNPs exhibiting large allele frequency differences over 900,000 genotyped SNPs. A region in chromosome 10 (within the cytogenetic bands q22.3 and q23.1) was significantly differentiated between highland and lowland groups. We resequenced ~1.5 Mb surrounding the candidate region and identified strong signals of positive selection in the highland populations. A composite of multiple signals like test localized the signal to FAM213A and a related enhancer; the product of this gene acts as an antioxidant to lower oxidative stress and may help to maintain bone mass. The results suggest that positive selection on the enhancer might increase the expression of this antioxidant, and thereby prevent oxidative damage. In addition, the most significant signal in a relative extended haplotype homozygosity analysis was localized to the SFTPD gene, which encodes a surfactant pulmonary-associated protein involved in normal respiration and innate host defense. Our study thus identifies two novel candidate genes and associated pathways that may be involved in high-altitude adaptation in Andean populations.

Prevalence, clinical profile, iron status, and subject-specific traits for excessive erythrocytosis in andean adults living permanently at 3,825 meters above sea level

BACKGROUND

Excessive erythrocytosis (EE) is a prevalent condition in populations living at high altitudes (> 2,500 m above sea level). Few large population-based studies have explored the association between EE and multiple subject-specific traits including oxygen saturation, iron status indicators, and pulmonary function.

METHODS

We enrolled a sex-stratified and age-stratified sample of 1,065 high-altitude residents aged ≥ 35 years from Puno, Peru (3,825 m above sea level) and conducted a standardized questionnaire and physical examination that included spirometry, pulse oximetry, and a blood sample for multiple clinical markers. Our primary objectives were to estimate the prevalence of EE, characterize the clinical profile and iron status indicators of subjects with EE, and describe subject-specific traits associated with EE.

RESULTS

Overall prevalence of EE was 4.5% (95% CI, 3.3%-6.0%). Oxygen saturation was significantly lower among EE than non-EE group subjects (85.3% vs 90.1%, P < .001) but no difference was found in iron status indicators between both groups (P > .09 for all values). In multivariable logistic regression, we found that age ≥ 65 years (OR = 2.45, 95% CI, 1.16-5.09), male sex (3.86, 1.78-9.08), having metabolic syndrome (2.66, 1.27-5.75) or being overweight (5.20, 1.95-16.77), pulse oximetry < 85% (14.90, 6.43-34.90), and % predicted FVC < 80% (13.62, 4.40-41.80) were strongly associated with EE. Attributable fractions for EE were greatest for being overweight (26.7%), followed by male sex (21.5%), pulse oximetry < 85% (16.4%), having metabolic syndrome (14.4%), and % predicted FVC < 80% (9.3%).

CONCLUSIONS

We found a lower prevalence of EE than in previous reports in the Peruvian Andes. Although the presence of hypoxemia and decreased vital capacity were strongly associated with excessive erythrocytosis, being overweight or having metabolic syndrome were associated with an important fraction of cases in our study population.

Population data of 15 autosomal STR markers from Afro-Bolivians of Nor Yungas Province (Bolivia)

Allele frequencies and forensic parameters for 15 autosomal loci included in the AmpFlSTR® Identifiler kit were estimated in a sample of 57 unrelated Afro-descendants from Nor Yungas (Bolivia). Buccal swabs samples were obtained from voluntary donors, after consent was given. All loci were in Hardy-Weinberg equilibrium after Bonferroni correction. D21S11 was the most informative locus, while the least discriminating locus was D3S1358. The combined power of discrimination and the combined probability of exclusion were >0.99999999 and >0.99997, respectively. The multidimensional scaling (MDS) plot generated by Rst matrix supported that Afro-Bolivians of Nor Yungas preserved a stronger African descent compared to other admixed Latin American populations. These results amplified the Bolivian databases of autosomal STR loci and may provide a useful tool for human identification tests and population genetic studies.

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.

Association between Y haplogroups and autosomal AIMs reveals intra-population substructure in Bolivian populations

For the correct evaluation of the weight of genetic evidence in a forensic context, databases must reflect the structure of the population, with all possible groups being represented. Countries with a recent history of admixture between strongly differentiated populations are usually highly heterogeneous and sub-structured. Bolivia is one of these countries, with a high diversity of ethnic groups and different levels of admixture (among Native Americans, Europeans and Africans) across the territory. For a better characterization of the male lineages in Bolivia, 17 Y-STR and 42 Y-SNP loci were genotyped in samples from La Paz and Chuquisaca. Only European and Native American Y-haplogroups were detected, and no sub-Saharan African chromosomes were found. Significant differences were observed between the two samples, with a higher frequency of European lineages in Chuquisaca than in La Paz. A sample belonging to haplogroup Q1a3a1a1-M19 was detected in La Paz, in a haplotype background different from those previously found in Argentina. This result supports an old M19 North-south dispersion in South America, possibly via two routes. When comparing the ancestry of each individual assessed through his Y chromosome with the one estimated using autosomal AIMs, (a) increased European ancestry in individuals with European Y chromosomes and (b) higher Native American ancestry in the carriers of Native American Y-haplogroups were observed, revealing an association between autosomal and Y-chromosomal markers. The results of this study demonstrate that a sub-structure does exist in Bolivia at both inter- and intrapopulation levels, a fact which must be taken into account in the evaluation of forensic genetic evidence.

The Andean adaptive toolkit to counteract high altitude maladaptation: genome-wide and phenotypic analysis of the Collas

During their migrations out of Africa, humans successfully colonised and adapted to a wide range of habitats, including extreme high altitude environments, where reduced atmospheric oxygen (hypoxia) imposes a number of physiological challenges. This study evaluates genetic and phenotypic variation in the Colla population living in the Argentinean Andes above 3500 m and compares it to the nearby lowland Wichí group in an attempt to pinpoint evolutionary mechanisms underlying adaptation to high altitude hypoxia. We genotyped 730,525 SNPs in 25 individuals from each population. In genome-wide scans of extended haplotype homozygosity Collas showed the strongest signal around VEGFB, which plays an essential role in the ischemic heart, and ELTD1, another gene crucial for heart development and prevention of cardiac hypertrophy. Moreover, pathway enrichment analysis showed an overrepresentation of pathways associated with cardiac morphology. Taken together, these findings suggest that Colla highlanders may have evolved a toolkit of adaptative mechanisms resulting in cardiac reinforcement, most likely to counteract the adverse effects of the permanently increased haematocrit and associated shear forces that characterise the Andean response to hypoxia. Regulation of cerebral vascular flow also appears to be part of the adaptive response in Collas. These findings are not only relevant to understand the evolution of hypoxia protection in high altitude populations but may also suggest new avenues for medical research into conditions where hypoxia constitutes a detrimental factor.

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

The genetic history of indigenous populations of the Peruvian and Bolivian Altiplano: the legacy of the Uros

The Altiplano region of the South American Andes is marked by an inhospitable climate to which the autochthonous human populations adapted and then developed great ancient civilizations, such as the Tiwanaku culture and the Inca Empire. Since pre-Columbian times, different rulers established themselves around the Titicaca and Poopo Lakes. By the time of the arrival of Spaniards, Aymara and Quechua languages were predominant on the Altiplano under the rule of the Incas, although the occurrence of other spoken languages, such as Puquina and Uruquilla, suggests the existence of different ethnic groups in this region. In this study, we focused on the pre-Columbian history of the autochthonous Altiplano populations, particularly the Uros ethnic group, which claims to directly descend from the first settlers of the Andes, and some linguists suggest they might otherwise be related to Arawak speaking groups from the Amazon. Using phylogeographic, population structure and spatial genetic analyses of Y-chromosome and mtDNA data, we inferred the genetic relationships among Uros populations (Los Uros from Peru, Uru-Chipaya and Uru-Poopo from Bolivia), and compared their haplotype profiles with eight Aymara, nine Quechua and two Arawak (Machiguenga and Yanesha) speaking populations from Peru and Bolivia. Our results indicated that Uros populations stand out among the Altiplano populations, while appearing more closely related to the Aymara and Quechua from Lake Titicaca and surrounding regions than to the Amazon Arawaks. Moreover, the Uros populations from Peru and Bolivia are genetically differentiated from each other, indicating a high heterogeneity in this ethnic group. Finally, our results support the distinctive ancestry for the Uros populations of Peru and Bolivia, which are likely derived from ancient Andean lineages that were partially replaced during more recent farming expansion events and the establishment of complex civilizations in the Andes.

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.

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

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

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

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

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

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

On the origins, rapid expansion and genetic diversity of Native Americans from hunting-gatherers to agriculturalists

Given the importance of Y-chromosome haplogroup Q to better understand the source populations of contemporary Native Americans, we studied 8 biallelic and 17 microsatellite polymorphisms on the background of 128 Q Y-chromosomes from geographically targeted populations. The populations examined in this study include three from the Tuva Republic in Central Asia (Bai-Tai, Kungurtug, and Toora-Hem, n = 146), two from the northeastern tip of Siberia (New Chaplino and Chukchi, n = 32), and two from Mesoamerica (Mayans from Yucatan, Mexico n = 72, and Mayans from the Guatemalan Highlands, n = 43). We also see evidence of a dramatic Mesoamerican post-migration population growth in the ubiquitous and diverse Y-STR profiles of the Mayan and other Mesoamerican populations. In the case of the Mayans, this demographic growth was most likely fueled by the agricultural- and trade-based subsistence adopted during the Pre-Classic, Classic and Post-Classic periods of their empire. The limited diversity levels observed in the Altaian and Tuvinian regions of Central Asia, the lowest of all populations examined, may be the consequence of bottleneck events fostered by the spatial isolation and low effective population size characteristic of a nomadic lifestyle. Furthermore, our data illustrate how a sociocultural characteristic such as mode of subsistence may be of impact on the genetic structure of populations. We analyzed our genetic data using Multidimensional Scaling Analysis of populations, Principal Component Analysis of individuals, Median-joining networks of M242, M346, L54, and M3 individuals, age estimations based on microsatellite variation utilizing genealogical and evolutionary mutation rates/generation times and estimation of Y- STR average gene diversity indices.

Y-STR diversity and sex-biased gene flow among Caribbean populations

In the present study, we report, for the first time, the allele and haplotype frequencies of 17 Y-STR (Y-filer) loci in the populations of Haiti, Jamaica and the Bahamas (Abaco, Eleuthera, Exuma, Grand Bahama, Long Island and New Providence). This investigation was undertaken to assess the paternal genetic structure of the abovementioned Caribbean islands. A total of 607 different haplotypes were identified among the 691 males examined, of which 537 (88.5%) were unique. Haplotype diversities (HD) ranged from 0.989 in Long Island to 1.000 in Grand Bahama, with limited haplotype sharing observed among these Caribbean collections. Discriminatory capacity (DC) values were also high, ranging from 79.1% to 100% in Long Island and Grand Bahama, respectively, illustrating the capacity of this set of markers to differentiate between patrilineal related individuals within each population. Phylogenetic comparison of the Bahamian, Haitian and Jamaican groups with available African, European, East Asian and Native American populations reveals strong genetic ties with the continental African collections, a finding that corroborates our earlier work using autosomal STR and Y-chromosome binary markers. In addition, various degrees of sex-biased gene flow exhibiting disproportionately higher European paternal (as compared to autosomal) influences were detected in all Caribbean islands genotyped except for Abaco and Eleuthera. We attribute the presence or absence of asymmetric gene flow to unique, island specific demographic events and family structures.

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.

Apolipoprotein E/C1/C4/C2 gene cluster diversity in two native Andean populations: Aymaras and Quechuas

The APOE/C1/C4/C2 gene cluster presents high relevance in lipid metabolism and, therefore, has important epidemiological implications. Here, we study for the first time the variation patterns of 25 polymorphisms (10 short tandem repeats, STRs, and 15 single nucleotide polymorphismas, SNPs) in two native Andean samples from Bolivia (45 Aymaras and 45 Quechuas) as well as one European sample (n = 41) as external reference. We estimated diversity parameters, linkage disequilibrium patterns, population structure, and possible selective effects. In general, diversity was low and could be partly attributed to selection (probably due to its physiological importance), since the APOE/C1/C4/C2 region was highly conserved compared to the flanking genes in both Bolivians and Europeans. Moreover, the lower gene diversity in Bolivians compared to Europeans for some markers might indicate different demographic histories. Regarding the APOE isoforms, in addition to ɛ3 (94%) and ɛ4 (5%), isoform ɛ2 (1%) was also detected in Bolivians. In relation to previous hypotheses, our results support that genetic drift or founder effects rather than selection for increased cholesterol absorption are the main factors that have shaped the distribution of APOE isoforms observed in South America.

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.

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

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