Bayesian inferences suggest that Amazon Yunga Natives diverged from Andeans less than 5000 ybp: implications for South American prehistory

A genomic perspective on South American human history

It has generally been accepted that the current indigenous peoples of the Americas are derived from ancestors from northeastern Asia. The latter were believed to have spread into the American continent by the end of the Last Glacial Maximum. In this sense, a joint and in-depth study of the earliest settlement of East Asia and the Americas is required to elucidate these events accurately. The first Americans underwent an adaptation process to the Americas' vast environmental diversity, mediated by biological and cultural evolution and niche construction, resulting in enormous cultural diversity, a wealth of domesticated species, and extensive landscape modifications. Afterward, in the Late Holocene, the advent of intensive agricultural food production systems, sedentism, and climate change significantly reshaped genetic and cultural diversity across the continent, particularly in the Andes and Amazonia. Furthermore, starting around the end of the 15th century, European colonization resulted in massive extermination of indigenous peoples and extensive admixture. Thus, the present review aims to create a comprehensive picture of the main events involved in the formation of contemporary South American indigenous populations and the dynamics responsible for shaping their genetic diversity by integrating current genetic data with evidence from archeology, linguistics and other disciplines.

NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also release a genealogies simulator software used for different tests performed in this study.

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.

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 genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

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.

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

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

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

Population genetics of immune-related multilocus copy number variation in Native Americans

While multiallelic copy number variation (mCNV) loci are a major component of genomic variation, quantifying the individual copy number of a locus and defining genotypes is challenging. Few methods exist to study how mCNV genetic diversity is apportioned within and between populations (i.e. to define the population genetic structure of mCNV). These inferences are critical in populations with a small effective size, such as Amerindians, that may not fit the Hardy-Weinberg model due to inbreeding, assortative mating, population subdivision, natural selection or a combination of these evolutionary factors. We propose a likelihood-based method that simultaneously infers mCNV allele frequencies and the population structure parameter f, which quantifies the departure of homozygosity from the Hardy-Weinberg expectation. This method is implemented in the freely available software CNVice, which also infers individual genotypes using information from both the population and from trios, if available. We studied the population genetics of five immune-related mCNV loci associated with complex diseases (beta-defensins, CCL3L1/CCL4L1, FCGR3A, FCGR3B and FCGR2C) in 12 traditional Native American populations and found that the population structure parameters inferred for these mCNVs are comparable to but lower than those for single nucleotide polymorphisms studied in the same populations.

East of the Andes: The genetic profile of the Peruvian Amazon populations

OBJECTIVES

Assuming that the differences between the Andes and the Amazon rainforest at environmental and historical levels have influenced the distribution patterns of genes, languages, and cultures, the maternal and paternal genetic reconstruction of the Peruvian Amazon populations was used to test the relationships within and between these two extreme environments.

MATERIALS AND METHODS

We analyzed four Peruvian Amazon communities (Ashaninka, Huambisa, Cashibo, and Shipibo) for both Y chromosome (17 STRs and 8 SNPs) and mtDNA data (control region sequences, two diagnostic sites of the coding region, and one INDEL), and we studied their variability against the rest of South America.

RESULTS

We detected a high degree of genetic diversity in the Peruvian Amazon people, both for mtDNA than for Y chromosome, excepting for Cashibo people, who seem to have had no exchanges with their neighbors, in contrast with the others communities. The genetic structure follows the divide between the Andes and the Amazon, but we found a certain degree of gene flow between these two environments, as particularly emerged with the Y chromosome descent cluster's (DCs) analysis.

DISCUSSION

The Peruvian Amazon is home to an array of populations with differential rates of genetic exchanges with their neighbors and with the Andean people, depending on their peculiar demographic histories. We highlighted some successful Y chromosome lineages expansions originated in Peru during the pre-Columbian history which involved both Andeans and Amazon Arawak people, showing that at least a part of the Amazon rainforest did not remain isolated from those exchanges.

MtDNA and Y-chromosomal diversity in the Chachapoya, a population from the northeast Peruvian Andes-Amazon divide

OBJECTIVES

The ancient Chachapoya were an aggregate of several ethnic groups that shared a common language, religion, and material culture. They inhabited a territory at the juncture of the Andes and the Amazon basin. Their position between those ecozones most likely influenced their genetic composition. We attempted to better understand their population history by assessing the contemporary genetic diversity in the Chachapoya and three of their immediate neighbors (Huancas, Jivaro, and Cajamarca). We inferred signatures of demographic history and genetic affinities, and contrasted the findings with data from other populations on local and continental scales.

METHODS

We studied mitochondrial DNA (mtDNA; hypervariable segment [HVSI and HVSII]) and Y chromosome (23 short tandem repeats (STRs)) marker data in 382 modern individuals. We used Sanger sequencing for mtDNA and a commercially available kit for Y-chromosomal STR typing.

RESULTS

The Chachapoya had affinities with various populations of Andean and Amazonian origin. When examining the Native American component, the Chachapoya displayed high levels of genetic diversity. Together with other parameters, for example, large Tajima's D and Fu's Fs, the data indicated no drastic reduction of the population size in the past.

CONCLUSION

The high level of diversity in the Chachapoya, the lack of evidence of drift in the past, and genetic affinities with a broad range of populations in the Americas reflects an intricate population history in the region. The new genetic data from the Chachapoya indeed seems to point to a genetic complexity that is not yet resolved but beginning to be elucidated. Am. J. Hum. Biol. 28:857-867, 2016. © 2016Wiley Periodicals, Inc.

New native South American Y chromosome lineages

Many single-nucleotide polymorphisms (SNPs) in the non-recombining region of the human Y chromosome have been described in the last decade. High-coverage sequencing has helped to characterize new SNPs, which has in turn increased the level of detail in paternal phylogenies. However, these paternal lineages still provide insufficient information on population history and demography, especially for Native Americans. The present study aimed to identify informative paternal sublineages derived from the main founder lineage of the Americas-haplogroup Q-L54-in a sample of 1841 native South Americans. For this purpose, we used a Y-chromosomal genotyping multiplex platform and conventional genotyping methods to validate 34 new SNPs that were identified in the present study by sequencing, together with many Y-SNPs previously described in the literature. We updated the haplogroup Q phylogeny and identified two new Q-M3 and three new Q-L54*(xM3) sublineages defined by five informative SNPs, designated SA04, SA05, SA02, SA03 and SA29. Within the Q-M3, sublineage Q-SA04 was mostly found in individuals from ethnic groups belonging to the Tukanoan linguistic family in the northwest Amazon, whereas sublineage Q-SA05 was found in Peruvian and Bolivian Amazon ethnic groups. Within Q-L54*, the derived sublineages Q-SA03 and Q-SA02 were exclusively found among Coyaima individuals (Cariban linguistic family) from Colombia, while Q-SA29 was found only in Maxacali individuals (Jean linguistic family) from southeast Brazil. Furthermore, we validated the usefulness of several published SNPs among indigenous South Americans. This new Y chromosome haplogroup Q phylogeny offers an informative paternal genealogy to investigate the pre-Columbian history of South America.Journal of Human Genetics advance online publication, 31 March 2016; doi:10.1038/jhg.2016.26.

Pharmacogenomic assessment of Mexican and Peruvian populations

BACKGROUND

Clinically relevant polymorphisms often demonstrate population-specific allele frequencies. Central and South America remain largely uncategorized in the context of pharmacogenomics.

MATERIALS & METHODS

We assessed 15 polymorphisms from 12 genes (ABCB1 3435C>T, ABCG2 Q141K, CYP1B1*3, CYP2C19*2, CYP3A4*1B, CYP3A5*3C, ERCC1 N118N, ERCC2 K751Q, GSTP1 I105V, TPMT 238G>C, TPMT 460G>A, TPMT 719A>G, TYMS TSER, UGT1A1*28 and UGT1A1 -3156G>A) in 81 Peruvian and 95 Mexican individuals.

RESULTS

Six polymorphism frequencies differed significantly between the two populations: ABCB1 3435C>T, CYP1B1*3, GSTP1 I105V, TPMT 460G>A, UGT1A1*28 and UGT1A1 -3156G>A. The pattern of observed allele frequencies for all polymorphisms could not be accurately estimated from any single previously studied population.

CONCLUSION

This highlights the need to expand the scope of geographic data for use in pharmacogenomics studies.

Revisiting the Diego Blood Group System in Amerindians: Evidence for Gene-Culture Comigration

Six decades ago the DI*A allele of the Diego blood group system was instrumental in proving Native American populations originated from Siberia. Since then, it has received scant attention. The present study was undertaken to reappraise distribution of the DI*A allele in 144 Native American populations based on current knowledge. Using analysis of variance tests, frequency distribution was studied according to geographical, environmental, and cultural parameters. Frequencies were highest in Amazonian populations. In contrast, DI*A was undetectable in subarctic, Fuegian, Panamanian, Chaco and Yanomama populations. Closer study revealed a correlation that this unequal distribution was correlated with language, suggesting that linguistic divergence was a driving force in the expansion of DI*A among Native Americans. The absence of DI*A in circumpolar Eskimo-Aleut and Na-Dene speakers was consistent with a late migratory event confined to North America. Distribution of DI*A in subtropical areas indicated that gene and culture exchanges were more intense within than between ecozones. Bolstering the utility of classical genetic markers in biological anthropology, the present study of the expansion of Diego blood group genetic polymorphism in Native Americans shows strong evidence of gene-culture comigration.