The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Genetic Signatures of Positive Selection in Human Populations Adapted to High Altitude in Papua New Guinea

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

Adipose development is consistent across hunter-gatherers and diverges from western references

Despite agreement that humans have evolved to be unusually fat primates, adipose patterning among hunter-gatherers has received little empirical consideration. Here we consider the development of adiposity among four contemporary groups of hunter-gatherers, the Aka, Savanna Pumé, Ju'/Hoansi and Agta using multi-level generalized additive mixed modelling to characterize the growth of tricep skinfolds from early childhood through adolescence. In contrast to references, hunter-gatherers show several consistent patterns: (i) children are lean with little fat accumulation; (ii) no adiposity rebound at 5 years is evident; (iii) girls on average have built 90% of their body size, and reach menarche when adiposity is at its maximum velocity; and (iv) a metabolic trade-off is evident in young, but not older children, such that both boys and girls prioritize skeletal growth during middle childhood, a trade-off that diminishes during adolescence when height velocity increases in pace with fat accumulation. Consistent results across hunter-gatherers living in diverse environments suggest that these patterns reflect a general forager pattern of development. The findings provide a valuable baseline for adipose development not apparent from reference populations. We emphasize both generalized trends among hunter-gatherers, and that inter-populational differences point to the plasticity with which humans organize growth and development.

The need to diversify genomic studies: Insights from Andean highlanders and Amazonians

More globally diverse perspectives are needed in genomic studies and precision medicine practices on non-Europeans. Here, we illustrate this by discussing the distribution of clinically actionable genetic variants involved in drug response in Andean highlanders and Amazonians, considering their environment, history, genetic structure, and historical biases in the perception of biological diversity of Native Americans.

Dietary oily fish intake is inversely associated with moderate-to-severe intracranial artery stenosis in older adults of indigenous Ecuadorian ancestry

BACKGROUND AND AIMS

Information on the association between dietary oily fish intake and intracranial atherosclerosis is limited and contradictory. Inconsistencies might be in part related to heterogeneous designs and differences in race/ethnicity of study populations. We aim to assess whether oily fish intake is inversely associated with intracranial artery stenosis (ICAS) in frequent fish consumers of indigenous ancestry living in coastal Ecuador.

METHODS

The study included 384 participants aged ≥60 years enrolled in the Atahualpa Project Cohort. Dietary oily fish intake was quantified systematically via validated surveys and all participants received a time-of-flight MRA of intracranial vessels. Poisson regression models, adjusted for demographics, level of education, traditional risk factors and severe tooth loss, were fitted to assess the association between amounts of oily fish intake and the number of intracranial arteries with moderate-to-severe (≥50 %) stenosis.

RESULTS

Participants had a mean age of 67.7 ± 7 years, and 56 % were women. The mean oily fish intake was 8.9 ± 5.2 servings/week; 283 (74 %) participants consumed ≥5.2 servings/week (2nd to 4th quartiles of fish intake). Forty-three (11 %) participants had at least one major intracranial artery with moderate-to-severe stenosis. Both univariate and multivariate models showed a significant inverse association between consumption of oily fish in the 2nd to 4th quartiles and ≥50 % stenosis in at least one artery (β: 0.46; 95 % C.I.: 0.27-077, and β: 0.52; 95 % C.I.: 0.30-0.90, respectively).

CONCLUSIONS

Consumption of more than five oily fish servings/week is associated with lower prevalence of moderate-to-severe ICAS in indigenous Ecuadorians.

Genetic portrait of the Amazonian communities of Peru and Bolivia: The legacy of the Takanan-speaking people

During the colonial period in South America, many autochthonous populations were affected by relocation by European missionary reductions and other factors that impacted and reconfigured their genetic makeup. Presently, the descendants of some "reduced" and other isolated groups are distributed in the Amazonian areas of Peru, Bolivia, and Brazil, and among them, speakers of Takanan and Panoan languages. Based on linguistics, these peoples should be closely related, but so far no DNA comparison studies have been conducted to corroborate a genetic relationship. To clarify these questions, we used a set of 15 short tandem repeats of the non-recombining part of the Y-chromosome (Y-STRs) and mitochondrial DNA (mtDNA) control region sequence data. Paternal line comparisons showed the Takanan-speaking peoples from Peru and Bolivia descended from recent common ancestors; one group was related to Arawakan, Jivaroan, and Cocama and the other to Panoan speakers, consistent with linguistics. Also, a genetic affinity for maternal lines was observed between some Takanan speakers and individuals who spoke different Amazonian languages. Our results supported a shared ancestry of Takanan, Panoan, Cocama, and Jivaroan-speaking communities who appeared to be related to each other and came likely from an early Arawak expansion in the western Amazonia of South America.

The multifaceted genomic history of Ashaninka from Amazonian Peru

Despite its crucial location, the western side of Amazonia between the Andes and the source(s) of the Amazon River is still understudied from a genomic and archaeogenomic point of view, albeit possibly harboring essential information to clarify the complex genetic history of local Indigenous groups and their interactions with nearby regions,^{1,2,3,4,5,6,7,8} including central America and the Caribbean.^9,10,11,12 Focusing on this key region, we analyzed the genome-wide profiles of 51 Ashaninka individuals from Amazonian Peru, observing an unexpected extent of genomic variation. We identified at least two Ashaninka subgroups with distinctive genomic makeups, which were differentially shaped by the degree and timing of external admixtures, especially with the Indigenous groups from the Andes and the Pacific coast. On a continental scale, Ashaninka ancestors probably derived from a south-north migration of Indigenous groups moving into the Amazonian rainforest from a southeastern area with contributions from the Southern Cone and the Atlantic coast. These ancestral populations diversified in the variegated geographic regions of interior South America, on the eastern side of the Andes, differentially interacting with surrounding coastal groups. In this complex scenario, we also revealed strict connections between the ancestors of present-day Ashaninka, who belong to the Arawakan language family,¹³ and those Indigenous groups that moved further north into the Caribbean, contributing to the early Ceramic (Saladoid) tradition in the islands.^14,15.

Indigenous people from Amazon show genetic signatures of pathogen-driven selection

Ecological conditions in the Amazon rainforests are historically favorable for the transmission of numerous tropical diseases, especially vector-borne diseases. The high diversity of pathogens likely contributes to the strong selective pressures for human survival and reproduction in this region. However, the genetic basis of human adaptation to this complex ecosystem remains unclear. This study investigates the possible footprints of genetic adaptation to the Amazon rainforest environment by analyzing the genomic data of 19 native populations. The results based on genomic and functional analysis showed an intense signal of natural selection in a set of genes related to Trypanosoma cruzi infection, which is the pathogen responsible for Chagas disease, a neglected tropical parasitic disease native to the Americas that is currently spreading worldwide.

A review of ancestrality and admixture in Latin America and the caribbean focusing on native American and African descendant populations

Genomics can reveal essential features about the demographic evolution of a population that may not be apparent from historical elements. In recent years, there has been a significant increase in the number of studies applying genomic epidemiological approaches to understand the genetic structure and diversity of human populations in the context of demographic history and for implementing precision medicine. These efforts have traditionally been applied predominantly to populations of European origin. More recently, initiatives in the United States and Africa are including more diverse populations, establishing new horizons for research in human populations with African and/or Native ancestries. Still, even in the most recent projects, the under-representation of genomic data from Latin America and the Caribbean (LAC) is remarkable. In addition, because the region presents the most recent global miscegenation, genomics data from LAC may add relevant information to understand population admixture better. Admixture in LAC started during the colonial period, in the 15th century, with intense miscegenation between European settlers, mainly from Portugal and Spain, with local indigenous and sub-Saharan Africans brought through the slave trade. Since, there are descendants of formerly enslaved and Native American populations in the LAC territory; they are considered vulnerable populations because of their history and current living conditions. In this context, studying LAC Native American and African descendant populations is important for several reasons. First, studying human populations from different origins makes it possible to understand the diversity of the human genome better. Second, it also has an immediate application to these populations, such as empowering communities with the knowledge of their ancestral origins. Furthermore, because knowledge of the population genomic structure is an essential requirement for implementing genomic medicine and precision health practices, population genomics studies may ensure that these communities have access to genomic information for risk assessment, prevention, and the delivery of optimized treatment; thus, helping to reduce inequalities in the Western Hemisphere. Hoping to set the stage for future studies, we review different aspects related to genetic and genomic research in vulnerable populations from LAC countries.

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

The loss of biodiversity in Madagascar is contemporaneous with major demographic events

Only 400 km off the coast of East Africa, the island of Madagascar is one of the last large land masses to have been colonized by humans. While many questions surround the human occupation of Madagascar, recent studies raise the question of human impact on endemic biodiversity and landscape transformation. Previous genetic and linguistic analyses have shown that the Malagasy population has emerged from an admixture that happened during the last millennium, between Bantu-speaking African populations and Austronesian-speaking Asian populations. By studying the sharing of chromosome segments between individuals (IBD determination), local ancestry information, and simulated genetic data, we inferred that the Malagasy ancestral Asian population was isolated for more than 1,000 years with an effective size of just a few hundred individuals. This isolation ended around 1,000 years before present (BP) by admixture with a small African population. Around the admixture time, there was a rapid demographic expansion due to intrinsic population growth of the newly admixed population, which coincides with extensive changes in Madagascar's landscape and the extinction of all endemic large-bodied vertebrates. Therefore, our approach can provide new insights into past human demography and associated impacts on ecosystems.

Global and local ancestry modulate APOE association with Alzheimer's neuropathology and cognitive outcomes in an admixed sample

Dementia is more prevalent in Blacks than in Whites, likely due to a combination of environmental and biological factors. Paradoxically, clinical studies suggest an attenuation of APOE ε4 risk of dementia in African ancestry (AFR), but a dearth of neuropathological data preclude the interpretation of the biological factors underlying these findings, including the association between APOE ε4 risk and Alzheimer's disease (AD) pathology, the most frequent cause of dementia. We investigated the interaction between African ancestry, AD-related neuropathology, APOE genotype, and functional cognition in a postmortem sample of 400 individuals with a range of AD pathology severity and lack of comorbid neuropathology from a cohort of community-dwelling, admixed Brazilians. Increasing proportions of African ancestry (AFR) correlated with a lower burden of neuritic plaques (NP). However, for individuals with a severe burden of NP and neurofibrillary tangles (NFT), AFR proportion was associated with worse Clinical Dementia Rating sum of boxes (CDR-SOB). Among APOE ε4 carriers, the association between AFR proportion and CDR-SOB disappeared. APOE local ancestry inference of a subset of 309 individuals revealed that, in APOE ε4 noncarriers, non-European APOE background correlated with lower NP burden and, also, worse cognitive outcomes than European APOE when adjusting by NP burden. Finally, APOE ε4 was associated with worse AD neuropathological burden only in a European APOE background. APOE genotype and its association with AD neuropathology and clinical pattern are highly influenced by ancestry, with AFR associated with lower NP burden and attenuated APOE ε4 risk compared to European ancestry.

Ortholog genes from cactophilic Drosophila provide insight into human adaptation to hallucinogenic cacti

Cultural transformations of lifestyles and dietary practices have been key drivers of human evolution. However, while most of the evidence of genomic adaptations is related to the hunter-gatherer transition to agricultural societies, little is known on the influence of other major cultural manifestations. Shamanism is considered the oldest religion that predominated throughout most of human prehistory and still prevails in many indigenous populations. Several lines of evidence from ethno-archeological studies have demonstrated the continuity and importance of psychoactive plants in South American cultures. However, despite the well-known importance of secondary metabolites in human health, little is known about its role in the evolution of ethnic differences. Herein, we identified candidate genes of adaptation to hallucinogenic cactus in Native Andean populations with a long history of shamanic practices. We used genome-wide expression data from the cactophilic fly Drosophila buzzatii exposed to a hallucinogenic columnar cactus, also consumed by humans, to identify ortholog genes exhibiting adaptive footprints of alkaloid tolerance. Genomic analyses in human populations revealed a suite of ortholog genes evolving under recent positive selection in indigenous populations of the Central Andes. Our results provide evidence of selection in genetic variants related to alkaloids toxicity, xenobiotic metabolism, and neuronal plasticity in Aymara and Quechua populations, suggesting a possible process of gene-culture coevolution driven by religious practices.

Uncovering signals of positive selection in Peruvian populations from three ecological regions

Perú hosts extremely diverse ecosystems which can be broadly classified into three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide SNP-array data from 286 Peruvians living across the three ecoregions and analysed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44 and ICAM1 in the coast; CBLB and PRDM1 in rainforest and the BRD2- HLA-DOA- HLA-DPA1 region in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Potential distributions of pre-Columbian people in Tropical Andean landscapes

Much has yet to be learned of the spatial patterning of pre-Columbian people across the Tropical Andes. Using compiled archaeological data and a suite of environmental variables, we generate an ensemble species distribution model (SDM) that incorporates general additive models, random forest models and Maxent models to reconstruct spatial patterns of pre-Columbian people that inhabited the Tropical Andes east of the continental divide, within the modern countries of Bolivia, Peru and Ecuador. Within this region, here referred to as the eastern Andean flank, elevation, mean annual cloud frequency, distance to rivers and precipitation of the driest quarter are the environmental variables most closely related to human occupancy. Our model indicates that 11.04% of our study area (65 368 km²) was likely occupied by pre-Columbian people. Our model shows that 30 of 351 forest inventory plots, which are used to generate ecological understanding of Andean ecosystems, were likely occupied in the pre-Columbian period. In previously occupied sites, successional trajectories may still be shaping forest dynamics, and those forests may still be recovering from the ecological legacy of pre-Columbian impacts. Our ensemble SDM links palaeo- and neo-ecology and can also be used to guide both future archaeological and ecological studies. This article is part of the theme issue 'Tropical forests in the deep human past'.

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.

Population genetics of PDE4B (Phosphodiesterase-4B) in neglected native americans: implications for cancer pharmacogenetics

PDE4B (phosphodiesterase‐4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US‐Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra‐phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non‐admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B‐intron‐7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra‐phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US‐Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

Peru - Progress in health and sciences in 200 years of independence

Peru celebrates 200 years of independence in 2021. Over this period of independent life, and despite the turbulent socio-political scenarios, from internal armed conflict to economic crisis to political instability over the last 40 years, Peru has experienced major changes on its epidemiological and population health profile. Major advancements in maternal and child health as well as in communicable diseases have been achieved in recent decades, and today Peru faces an increasing burden of non-communicable diseases including mental health conditions. In terms of the configuration of the public health system, Peru has also strived to secure country-wide optimal health care, struggling in particular to improve primary health care and intercultural services. The science and technology infrastructure has also evolved, although the need for substantial investments remains if advancing science is to be a national priority. Climate change will also bring significant challenges to population health given Peru's geographical and microclimates diversity. Looking back over the 200-years of independence, we present a summary of key advances in selected health-related fields, thus serving as the basis for reflections on pending agendas and future challenges, in order to look forward to ensuring the future health and wellbeing of the Peruvian population.

Resumen translated abstract

El Perú cumple 200 años de independencia en 2021. Durante estos dos siglos de vida independiente, junto con periodos sociales y políticos turbulentos, incluyendo un conflicto armado interno, hiperinflación y la inestabilidad política de los últimos 40 años, el Perú ha experimentado importantes cambios en su perfil epidemiológico con repercusiones directas en la salud de la población. En las últimas décadas, los indicadores de salud materno-infantil y de las enfermedades transmisibles muestran mejoría importante, pero el país se enfrenta de manera simultánea a una carga cada vez mayor de enfermedades no transmisibles y de salud mental. En cuanto a los sistemas de salud pública, se han realizado esfuerzos por aumentar la cobertura y calidad de la atención de salud en todo el país, apostándose en particular por mejorar la atención primaria. La ciencia y tecnología relacionadas con la salud también han mejorado, aunque si se quiere que la ciencia sea una prioridad nacional, son necesarias inversiones sustanciales. El cambio climático traerá importantes desafíos para la salud de la población, dada la diversidad geográfica y de microclimas del país. Para conmemorar los 200 años de vida independiente del Perú, presentamos un resumen de avances clave en diversas áreas y temas relacionados con la salud. Este repaso sirve como base para reflexionar sobre agendas y desafíos pendientes y futuros, con el fin de asegurar la salud y el bienestar de la población peruana en las próximas décadas.

Climate and demography drive 7000 years of dietary change in the Central Andes

Explaining the factors that influence past dietary variation is critically important for understanding changes in subsistence, health, and status in past societies; yet systematic studies comparing possible driving factors remain scarce. Here we compile the largest dataset of past diet derived from stable isotope δ¹³C‰ and δ¹⁵N‰ values in the Americas to quantitatively evaluate the impact of 7000 years of climatic and demographic change on dietary variation in the Central Andes. Specifically, we couple paleoclimatic data from a general circulation model with estimates of relative past population inferred from archaeologically derived radiocarbon dates to assess the influence of climate and population on spatiotemporal dietary variation using an ensemble machine learning model capable of accounting for interactions among predictors. Results reveal that climate and population strongly predict diet (80% of δ15N‰ and 66% of δ13C‰) and that Central Andean diets correlate much more strongly with local climatic conditions than regional population size, indicating that the past 7000 years of dietary change was influenced more by climatic than socio-demographic processes. Visually, the temporal pattern suggests decreasing dietary variation across elevation zones during the Late Horizon, raising the possibility that sociopolitical factors overrode the influence of local climatic conditions on diet during that time. The overall findings and approach establish a general framework for understanding the influence of local climate and demography on dietary change across human history.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

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.

Identification of New Helicobacter pylori Subpopulations in Native Americans and Mestizos From Peru

Region-specific Helicobacter pylori subpopulations have been identified. It is proposed that the hspAmerind subpopulation is being displaced from the Americans by an hpEurope population following the conquest. Our study aimed to describe the genomes and methylomes of H. pylori isolates from distinct Peruvian communities: 23 strains collected from three groups of Native Americans (Asháninkas [ASHA, n = 9], Shimaas [SHIM, n = 5] from Amazonas, and Punos from the Andean highlands [PUNO, n = 9]) and 9 modern mestizos from Lima (LIM). Closed genomes and DNA modification calls were obtained using SMRT/PacBio sequencing. We performed evolutionary analyses and evaluated genomic/epigenomic differences among strain groups. We also evaluated human genome-wide data from 74 individuals from the selected Native communities (including the 23 H. pylori strains donors) to compare host and bacterial backgrounds. There were varying degrees of hspAmerind ancestry in all strains, ranging from 7% in LIM to 99% in SHIM. We identified three H. pylori subpopulations corresponding to each of the Native groups and a novel hspEuropePeru which evolved in the modern mestizos. The divergence of the indigenous H. pylori strains recapitulated the genetic structure of Native Americans. Phylogenetic profiling showed that Orthogroups in the indigenous strains seem to have evolved differentially toward epigenomic regulation and chromosome maintenance, whereas OGs in the modern mestizo (LIM) seem to have evolved toward virulence and adherence. The prevalence of cagA⁺/vacA s1i1m1 genotype was similar across populations (p = 0.32): 89% in ASHA, 67% in PUNO, 56% in LIM and 40% in SHIM. Both cagA and vacA sequences showed that LIM strains were genetically differentiated (p < 0.001) as compared to indigenous strains. We identified 642 R-M systems with 39% of the associated genes located in the core genome. We found 692 methylation motifs, including 254 population-specific sequences not previously described. In Peru, hspAmerind is not extinct, with traces found even in a heavily admixed mestizo population. Notably, our study identified three new hspAmerind subpopulations, one per Native group; and a new subpopulation among mestizos that we named hspEuropePeru. This subpopulation seems to have more virulence-related elements than hspAmerind. Purifying selection driven by variable host immune response may have shaped the evolution of Peruvian subpopulations, potentially impacting disease outcomes.