Reconstructing Native American population history. Nature. 2012;488:370-374. [PMID: 22801491 PMCID: PMC3615710 DOI: 10.1038/nature11258]

9,000 years of genetic continuity in southernmost Africa demonstrated at Oakhurst rockshelter

Southern Africa has one of the longest records of fossil hominins and harbours the largest human genetic diversity in the world. Yet, despite its relevance for human origins and spread around the globe, the formation and processes of its gene pool in the past are still largely unknown. Here, we present a time transect of genome-wide sequences from nine individuals recovered from a single site in South Africa, Oakhurst Rockshelter. Spanning the whole Holocene, the ancient DNA of these individuals allows us to reconstruct the demographic trajectories of the indigenous San population and their ancestors during the last 10,000 years. We show that, in contrast to most regions around the world, the population history of southernmost Africa was not characterized by several waves of migration, replacement and admixture but by long-lasting genetic continuity from the early Holocene to the end of the Later Stone Age. Although the advent of pastoralism and farming substantially transformed the gene pool in most parts of southern Africa after 1,300 BP, we demonstrate using allele-frequency and identity-by-descent segment-based methods that the ‡Khomani San and Karretjiemense from South Africa still show direct signs of relatedness to the Oakhurst hunter-gatherers, a pattern obscured by recent, extensive non-Southern African admixture. Yet, some southern San in South Africa still preserve this ancient, Pleistocene-derived genetic signature, extending the period of genetic continuity until today.

Disparate and parallel craniofacial climatic adaptations in native populations of Asia, North America, and South America

Understanding the impact that climate had in shaping cranial variation is critical for inferring the evolutionary mechanisms that played a role in human diversification. Here, we provide a comprehensive study aiming to analyze the association between climate and cranial variation of high latitude populations living in temperate to cold environments of Asia, North America, and South America. For this, we compiled a large morphometric dataset (N = 2633), which was combined with climatic and genomic data. We tested the influence of climate on the facial skeleton, nasal protrusion, and cranial vault and through multiple statistical tests at two geographical scales: intracontinental and intercontinental. We show that populations living in cold areas share a morphological pattern characterized by an increase in nasal height, facial and orbital heights and widths, a decrease in facial protrusion, and larger, longer, and lower cranial vaults. There are also distinctive features; populations from north Asia present the tallest noses, largest faces, and cranial vaults of the whole sample. Nasal breadth dimensions show small values in Asians, large values in South Americans, and non-significant changes in arctic North America. The morphological pattern in populations living at high latitude may be the result of parallel adaptation, as supported by physiological, morphometric, ecological, and genetic explanations, while the differences in magnitude and phenotypic expression could be due to the diverse population histories, severity of climate, and cultural strategies. Overall, our study shows that climate is a relevant factor shaping modern human morphology and it should be considered when studying modern human evolution and diversification.

Unveiling ancestral threads: Exploring CCR5 ∆32 mutation frequencies in Colombian populations for HIV/AIDS therapeutics

AIDS remains a significant global health challenge since its emergence in 1981, with millions of deaths and new cases every year. The CCR5 ∆32 genetic deletion confers immunity to HIV infection by altering a cell membrane protein crucial for viral entry. Stem cell transplants from homozygous carriers of this mutation to HIV-infected individuals have resulted in viral load reduction and disease remission, suggesting a potential therapeutic avenue. This study aims to investigate the relationship between genetic ancestry and the frequency of the CCR5 ∆32 mutation in Colombian populations, exploring the feasibility of targeted donor searches based on ancestry composition. Utilizing genomic data from the CÓDIGO-Colombia consortium, comprising 532 individuals, the study assessed the presence of the CCR5 ∆32 mutation and examined if the population was on Hardy-Weinberg equilibrium. Individuals were stratified into clusters based on African, American, and European ancestry percentages, with logistic regression analysis performed to evaluate the association between ancestry and mutation frequency. Additionally, global genomic databases were utilized to visualize the worldwide distribution of the mutation. The findings revealed a significant positive association between European ancestry and the CCR5 ∆32 mutation frequency, underscoring its relevance in donor selection. African and American ancestry showed negative but non-significant associations with CCR5 ∆32 frequency, which may be attributed to the study's limitations. These results emphasize the potential importance of considering ancestry in donor selection strategies, reveal the scarcity of potential donors in Colombia, and underscore the need to consider donors from other populations with mainly European ancestry if the CCR5 ∆32 stem cell transplant becomes a routine treatment for HIV/AIDS in Colombia.

Dense SNP-based analyses complement forensic anthropology biogeographical ancestry assessments

Identification of unidentified human remains (UHRs) is crucial yet challenging, especially with traditional forensic techniques. Forensic anthropological examinations can yield ancestry estimations; however, the utility of these estimates is limited by the data points that can be collected from partial remains, complexities of admixture, and variation of phenotypic expression due to environmental effects. While it is generally known that anthropological estimates can be imprecise, the performance of these methods has not been studied at scale. Genome-wide SNP testing is an orthogonal approach for estimating ancestry and offers a unique opportunity to measure the magnitude of anthropological ancestry misattribution. Genomic ancestry inference leverages principal component analysis (PCA) and model-based clustering approaches. This study compares anthropologically determined ancestry with those estimated using genome-wide SNP markers. A dataset of 611 UHR samples with publicly available ancestry assessments from National Missing and Unidentified Persons System (NamUs) was analyzed. The genetic ancestry approach, validated against reference population samples, offers robust ancestry calculations for major population groups. Inconsistency between anthropological and genomic ancestry assignments were observed, particularly for admixed populations. Although forensic anthropological examinations remain valuable, their limitations emphasize the need for refinement and enhancement through the augmentation of SNP-based analyses. Further validation studies are crucial to define the uncertainty associated with both anthropological and genome-based ancestry estimates to resolve cases and aid law enforcement investigations. Additionally, current policy and practices for reporting ancestry for UHRs should be revisited to reduce potential misinformation.

Testing times: disentangling admixture histories in recent and complex demographies using ancient DNA

Our knowledge of human evolutionary history has been greatly advanced by paleogenomics. Since the 2020s, the study of ancient DNA has increasingly focused on reconstructing the recent past. However, the accuracy of paleogenomic methods in resolving questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation remains an open question. We evaluated the performance and behavior of two commonly used methods, qpAdm and the f3-statistic, on admixture inference under a diversity of demographic models and data conditions. We performed two complementary simulation approaches-firstly exploring a wide demographic parameter space under four simple demographic models of varying complexities and configurations using branch-length data from two chromosomes-and secondly, we analyzed a model of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudohaploidization. We observe that population differentiation is the primary factor driving qpAdm performance. Notably, while complex gene flow histories influence which models are classified as plausible, they do not reduce overall performance. Under conditions reflective of the historical period, qpAdm most frequently identifies the true model as plausible among a small candidate set of closely related populations. To increase the utility for resolving fine-scaled hypotheses, we provide a heuristic for further distinguishing between candidate models that incorporates qpAdm model P-values and f3-statistics. Finally, we demonstrate a significant performance increase for qpAdm using whole-genome branch-length f2-statistics, highlighting the potential for improved demographic inference that could be achieved with future advancements in f-statistic estimations.

Performance of qpAdm -based screens for genetic admixture on admixture-graph-shaped histories and stepping-stone landscapes

qpAdm is a statistical tool that is often used in exploratory archaeogenetic studies for finding optimal admixture models of population history. Despite its popularity, qpAdm remains untested on histories in the form of admixture graphs of random topology or stepping-stone landscapes. We analyzed data from such simulations and found that while for admixture-graph-shaped histories there exist simple solutions (temporal stratification) for minimizing false findings of gene flow, in the case of stepping-stone landscapes the method generates results that do not appear suspect but are misleading: feasible qpAdm models are either accurate but simplistic in the context of landscapes, or highly inaccurate in the case of multi-component models. This is largely is due to two reasons: 1) because of complex migration networks that violate the assumptions of the method, there is poor correlation between qpAdm p -values and model optimality in many sections of the parameter space; 2) admixture fraction estimates between 0 and 1 are largely restricted to symmetric source configurations around targets, hence popular [0, 1] model plausibility criteria confound analyses of landscape-type demographies, unless their interpretations are explicitly spatial. For many species/regions/periods archaeogenetic sampling is very sparse and may be random with respect to population density of ancient individuals. In this situation only a specific combination of landscape properties and feasibility criteria allows to efficiently reject highly asymmetric non-optimal models most abundant in random deme sets. This problem may obscure useful signal (rare optimal models) and might be responsible for some claims about rapid long-distance migrations in the literature.

Extended haplotype with rs41524547-G defines the ancestral origin of SCA10

Spinocerebellar ataxia type 10 (SCA10) is a rare autosomal dominant ataxia caused by a large expansion of the (ATTCT)n repeat in ATXN10. SCA10 was described in Native American and Asian individuals which prompted a search for an expanded haplotype to confirm a common ancestral origin for the expansion event. All patients with SCA10 expansions in our cohort share a single haplotype defined at the 5'-end by the minor allele of rs41524547, located ~35 kb upstream of the SCA10 expansion. Intriguingly, rs41524547 is located within the miRNA gene, MIR4762, within its DROSHA cleavage site and just outside the seed sequence for mir4792-5p. The world-wide frequency of rs41524547-G is less than 5% and found almost exclusively in the Americas and East Asia-a geographic distribution that mirrors reported SCA10 cases. We identified rs41524547-G(+) DNA from the 1000 Genomes/International Genome Sample Resource and our own general population samples and identified SCA10 repeat expansions in up to 25% of these samples. The reduced penetrance of these SCA10 expansions may be explained by a young (pre-onset) age at sample collection, a small repeat size, purity of repeat units, or the disruption of miR4762-5p function. We conclude that rs41524547-G is the most robust at-risk SNP allele for SCA10, is useful for screening of SCA10 expansions in population genetics studies and provides the most compelling evidence to date for a single, prehistoric origin of SCA10 expansions sometime prior to or during the migration of individuals across the Bering Land Bridge into the Americas.

Ancient Rapanui genomes reveal resilience and pre-European contact with the Americas

Rapa Nui (also known as Easter Island) is one of the most isolated inhabited places in the world. It has captured the imagination of many owing to its archaeological record, which includes iconic megalithic statues called moai1. Two prominent contentions have arisen from the extensive study of Rapa Nui. First, the history of the Rapanui has been presented as a warning tale of resource overexploitation that would have culminated in a major population collapse-the 'ecocide' theory2-4. Second, the possibility of trans-Pacific voyages to the Americas pre-dating European contact is still debated5-7. Here, to address these questions, we reconstructed the genomic history of the Rapanui on the basis of 15 ancient Rapanui individuals that we radiocarbon dated (1670-1950 CE) and whole-genome sequenced (0.4-25.6×). We find that these individuals are Polynesian in origin and most closely related to present-day Rapanui, a finding that will contribute to repatriation efforts. Through effective population size reconstructions and extensive population genetics simulations, we reject a scenario involving a severe population bottleneck during the 1600s, as proposed by the ecocide theory. Furthermore, the ancient and present-day Rapanui carry similar proportions of Native American admixture (about 10%). Using a Bayesian approach integrating genetic and radiocarbon dates, we estimate that this admixture event occurred about 1250-1430 CE.

Evidence for dynastic succession among early Celtic elites in Central Europe

The early Iron Age (800 to 450 BCE) in France, Germany and Switzerland, known as the 'West-Hallstattkreis', stands out as featuring the earliest evidence for supra-regional organization north of the Alps. Often referred to as 'early Celtic', suggesting tentative connections to later cultural phenomena, its societal and population structure remain enigmatic. Here we present genomic and isotope data from 31 individuals from this context in southern Germany, dating between 616 and 200 BCE. We identify multiple biologically related groups spanning three elite burials as far as 100 km apart, supported by trans-regional individual mobility inferred from isotope data. These include a close biological relationship between two of the richest burial mounds of the Hallstatt culture. Bayesian modelling points to an avuncular relationship between the two individuals, which may suggest a practice of matrilineal dynastic succession in early Celtic elites. We show that their ancestry is shared on a broad geographic scale from Iberia throughout Central-Eastern Europe, undergoing a decline after the late Iron Age (450 BCE to ~50 CE).

Towards cascading genetic risk in Alzheimer's disease

Alzheimer's disease typically progresses in stages, which have been defined by the presence of disease-specific biomarkers: amyloid (A), tau (T) and neurodegeneration (N). This progression of biomarkers has been condensed into the ATN framework, in which each of the biomarkers can be either positive (+) or negative (-). Over the past decades, genome-wide association studies have implicated ∼90 different loci involved with the development of late-onset Alzheimer's disease. Here, we investigate whether genetic risk for Alzheimer's disease contributes equally to the progression in different disease stages or whether it exhibits a stage-dependent effect. Amyloid (A) and tau (T) status was defined using a combination of available PET and CSF biomarkers in the Alzheimer's Disease Neuroimaging Initiative cohort. In 312 participants with biomarker-confirmed A-T- status, we used Cox proportional hazards models to estimate the contribution of APOE and polygenic risk scores (beyond APOE) to convert to A+T- status (65 conversions). Furthermore, we repeated the analysis in 290 participants with A+T- status and investigated the genetic contribution to conversion to A+T+ (45 conversions). Both survival analyses were adjusted for age, sex and years of education. For progression from A-T- to A+T-, APOE-e4 burden showed a significant effect [hazard ratio (HR) = 2.88; 95% confidence interval (CI): 1.70-4.89; P < 0.001], whereas polygenic risk did not (HR = 1.09; 95% CI: 0.84-1.42; P = 0.53). Conversely, for the transition from A+T- to A+T+, the contribution of APOE-e4 burden was reduced (HR = 1.62; 95% CI: 1.05-2.51; P = 0.031), whereas the polygenic risk showed an increased contribution (HR = 1.73; 95% CI: 1.27-2.36; P < 0.001). The marginal APOE effect was driven by e4 homozygotes (HR = 2.58; 95% CI: 1.05-6.35; P = 0.039) as opposed to e4 heterozygotes (HR = 1.74; 95% CI: 0.87-3.49; P = 0.12). The genetic risk for late-onset Alzheimer's disease unfolds in a disease stage-dependent fashion. A better understanding of the interplay between disease stage and genetic risk can lead to a more mechanistic understanding of the transition between ATN stages and a better understanding of the molecular processes leading to Alzheimer's disease, in addition to opening therapeutic windows for targeted interventions.

Coupling and decoupling of ancestral linkages and current cross-border economic activities: Genetics and policy

This paper studies the potential link between the biological evolution of populations and present-day economic interactions by estimating the correlation of shared ancestry among populations with cross-border capital and human flows. To this end, we employ the new concept of genetic distance, based on (dis)similarity of neutral gene alleles, to quantify shared ancestry. We then incorporate the genetic distance measure into an augmented gravity model, traditionally used to analyze the effect of geographical distance on bilateral exchange. Our analysis focuses on bilateral foreign direct investment (FDI) and migration across 135 countries and we use both linear regression techniques as well as the Poisson Pseudo-Maximum Likelihood Estimator to account for any non-linearities in the model. Our results show that a 1% increase in genetic distance reduces FDI flows by 0.08% while controlling for other distance constructs and factors associated with global capital and human movement. Genetic distance also has a negative effect on migration, where a 1% increase in genetic distance reduces migration flows by 0.22%, with all other things remaining constant. Our study, therefore, links shared ancestry with economic behavior, showing how historical connections are associated with current economic exchanges among nations. Additionally, recognizing that ancestral ties are outside human control, we examine policy measures that help nations overcome such distance barriers. Our findings show that strengthening a nation's institutional quality and adherence to the rule of law can effectively mitigate any negative correlation of distance constructs with economic exchanges. These insights suggest that prudent policies to foster a stable business environment are essential for any nation to attract FDI and human capital, even from geographically or genetically distant nations.

GWAS advancements to investigate disease associations and biological mechanisms

Genome-wide association studies (GWAS) have been instrumental in elucidating the genetic architecture of various traits and diseases. Despite the success of GWAS, inherent limitations such as identifying rare and ultra-rare variants, the potential for spurious associations, and in pinpointing causative agents can undermine diagnostic capabilities. This review provides an overview of GWAS and highlights recent advances in genetics that employ a range of methodologies, including Whole Genome Sequencing (WGS), Mendelian Randomization (MR), the Pangenome's high-quality T2T-CHM13 panel, and the Human BioMolecular Atlas Program (HuBMAP), as potential enablers of current and future GWAS research. State of the literature demonstrate the capabilities of these techniques in enhancing the statistical power of GWAS. WGS, with its comprehensive approach, captures the entire genome, surpassing the capabilities of the traditional GWAS technique focused on predefined Single Nucleotide Polymorphism (SNP) sites. The Pangenome's T2T-CHM13 panel, with its holistic approach, aids in the analysis of regions with high sequence identity, such as segmental duplications (SDs). Mendelian Randomization has advanced causative inference, improving clinical diagnostics and facilitating definitive conclusions. Furthermore, spatial biology techniques like HuBMAP, enable 3D molecular mapping of tissues at single-cell resolution, offering insights into pathology of complex traits. This study aims to elucidate and advocate for the increased application of these technologies, highlighting their potential to shape the future of GWAS research.

Ancient genomes reveal insights into ritual life at Chichén Itzá

The ancient city of Chichén Itzá in Yucatán, Mexico, was one of the largest and most influential Maya settlements during the Late and Terminal Classic periods (AD 600-1000) and it remains one of the most intensively studied archaeological sites in Mesoamerica1-4. However, many questions about the social and cultural use of its ceremonial spaces, as well as its population's genetic ties to other Mesoamerican groups, remain unanswered2. Here we present genome-wide data obtained from 64 subadult individuals dating to around AD 500-900 that were found in a subterranean mass burial near the Sacred Cenote (sinkhole) in the ceremonial centre of Chichén Itzá. Genetic analyses showed that all analysed individuals were male and several individuals were closely related, including two pairs of monozygotic twins. Twins feature prominently in Mayan and broader Mesoamerican mythology, where they embody qualities of duality among deities and heroes5, but until now they had not been identified in ancient Mayan mortuary contexts. Genetic comparison to present-day people in the region shows genetic continuity with the ancient inhabitants of Chichén Itzá, except at certain genetic loci related to human immunity, including the human leukocyte antigen complex, suggesting signals of adaptation due to infectious diseases introduced to the region during the colonial period.

Network of large pedigrees reveals social practices of Avar communities

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

Genomic investigation on genes related to mercury metabolism in Amazonian indigenous populations

Studies have identified elevated levels of mercury in Amazonian Indigenous individuals, highlighting them as one of the most exposed to risks. In the unique context of the Brazilian Indigenous population, it is crucial to identify genetic variants with clinical significance to better understand vulnerability to mercury and its adverse effects. Currently, there is a lack of research on the broader genomic profile of Indigenous people, particularly those from the Amazon region, concerning mercury contamination. Therefore, the aim of this study was to assess the genomic profile related to the processes of mercury absorption, distribution, metabolism, and excretion in 64 Indigenous individuals from the Brazilian Amazon. We aimed to determine whether these individuals exhibit a higher susceptibility to mercury exposure. Our study identified three high-impact variants (GSTA1 rs1051775, GSTM1 rs1183423000, and rs1241704212), with the latter two showing a higher frequency in the study population compared to global populations. Additionally, we discovered seven new variants with modifier impact and a genomic profile different from the worldwide populations. These genetic variants may predispose the study population to more harmful mercury exposure compared to global populations. As the first study to analyze broader genomics of mercury metabolism pathways in Brazilian Amazonian Amerindians, we emphasize that our research aims to contribute to public policies by utilizing genomic investigation as a method to identify populations with a heightened susceptibility to mercury exposure.

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Reconstructing the Genetic Relationship between Ancient and Present-Day Siberian Populations

Human populations across a vast area in northern Eurasia, from Fennoscandia to Chukotka, share a distinct genetic component often referred to as the Siberian ancestry. Most enriched in present-day Samoyedic-speaking populations such as Nganasans, its origins and history still remain elusive despite the growing list of ancient and present-day genomes from Siberia. Here, we reanalyze published ancient and present-day Siberian genomes focusing on the Baikal and Yakutia, resolving key questions regarding their genetic history. First, we show a long-term presence of a unique genetic profile in southern Siberia, up to 6,000 yr ago, which distinctly shares a deep ancestral connection with Native Americans. Second, we provide plausible historical models tracing genetic changes in West Baikal and Yakutia in fine resolution. Third, the Middle Neolithic individual from Yakutia, belonging to the Belkachi culture, serves as the best source so far available for the spread of the Siberian ancestry into Fennoscandia and Greenland. These findings shed light on the genetic legacy of the Siberian ancestry and provide insights into the complex interplay between different populations in northern Eurasia throughout history.

Biogeographic Perspectives on Human Genetic Diversification

Modern humans originated in Africa 300,000 yr ago, and before leaving their continent of origin, they underwent a process of intense diversification involving complex demographic dynamics. Upon exiting Africa, different populations emerged on the four other inhabited continents, shaped by the interplay of various evolutionary processes, such as migrations, founder effects, and natural selection. Within each region, continental populations, in turn, diversified and evolved almost independently for millennia. As a backdrop to this diversification, introgressions from archaic species contributed to establishing different patterns of genetic diversity in different geographic regions, reshaping our understanding of our species' variability. With the increasing availability of genomic data, it has become possible to delineate the subcontinental human population structure precisely. However, the bias toward the genomic research focused on populations from the global North has limited our understanding of the real diversity of our species and the processes and events that guided different human groups throughout their evolutionary history. This perspective is part of a series of articles celebrating 40 yr since our journal, Molecular Biology and Evolution, was founded (Russo et al. 2024). The perspective is accompanied by virtual issues, a selection of papers on human diversification published by Genome Biology and Evolution and Molecular Biology and Evolution.

New Y-SNPs in QM3 indigenous populations of Colombia

In evolutionary studies of human populations based on the Y chromosome, the majority of Native Americans belong to the QM3 lineage. Therefore, to study the history of groups inhabiting northern South America, it is necessary to have a higher resolution of the tree. The objective of this work was to identify new SNPs of the QM3 lineage that would allow the evaluation of the phylogenetic relationships between Andean and Amazonian populations of Colombia. Sequences previously obtained from two Y chromosomes of Amazonian populations were used, from which 13 potential SNPs were selected and typed in 171 Amazonian samples from the Vaupés region and in 60 samples from the Pasto, Nasa, Embera, Arhuaco and Kogüi ethnic groups of the Andean region. In addition, the main SNPs/markers (L56, L54, M346, M848, Z780, CTS11780) defining autochthonous Q lineages were typed, along with others defined by different SNPs/markers as reported in the literature (CTS11357, SA05, Z19319, Z5915, and Z19384). It was found that all the new SNPs are present in the Amazonian samples and only 2 of them are shared with the Embera, Nasa and Pasto, but none with the Kogüi and Arhuaco from the northern Andes, in the Colombian Caribbean. Combining the 13 variants of the present study with 14 previously reported and using TMRCA, a new QM3 tree proposal is generated. This method makes it possible to increase the number of sublineages of QM3 with a higher resolution and to detect differences between the different populations of Vaupés in the Amazon, as in the case of the Kubeos and Pisamiras, the latter of which is in grave danger of extinction. These new sublineages are useful for microevolutionary studies of the Amerindian populations of South America.

Peopling of the Americas: A new approach to assessing dental morphological variation in Asian and Native American populations

OBJECTIVES

Through biodistance analyses, anthropologists have used dental morphology to elucidate how people moved into and throughout the Americas. Here, we apply a method that focuses on individuals rather than sample frequencies through the application rASUDAS2, based on a naïve Bayes' algorithm.

MATERIALS AND METHODS

Using the database of C.G. Turner II, we calculated the probability that an individual could be assigned to one of seven biogeographic groups (American Arctic, North & South America, East Asia, Southeast Asia & Polynesia, Australo-Melanesia, Western Eurasia, & Sub-Saharan Africa) through rASUDAS2. The frequency of classifications for each biogeographic group was determined for 1418 individuals from six regions across Asia and the Americas.

RESULTS

Southeast Asians show mixed assignments but rarely to American Arctic or "American Indian." East Asians are assigned to East Asia half the time while 30% are assigned as Native American. People from the American Arctic and North & South America are assigned to Arctic America or non-Arctic America 75%-80% of the time, with 10%-15% classified as East Asian.

DISCUSSION

All Native American groups have a similar degree of morphological affinity to East Asia, as 10%-15% are classified as East Asian. East Asians are classified as Native American in 30% of cases. Individuals in the Western Hemisphere are decreasingly classified as Arctic the farther south they are located. Equivalent levels of classification as East Asian across all Native American groups suggests one divergence between East Asians and the population ancestral to all Native Americans. Non-arctic Native American groups are derived from the Arctic population, which represents the Native American founder group.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f3-statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P-values of qpAdm models, and f3-statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f2-statistics, demonstrating the improved demographic inference that could be achieved with advancements in f-statistic estimations.

Amerindian ancestry proportion as a risk factor for inflammatory bowel diseases: results from a Latin American Andean cohort

Background and aims

Latin American populations remain underrepresented in genetic studies of inflammatory bowel diseases (IBDs). Most genetic association studies of IBD rely on Caucasian, African, and Asian individuals. These associations have yet to be evaluated in detail in the Andean region of South America. We explored the contribution of IBD-reported genetic risk variants to a Chilean cohort and the ancestry contribution to IBD in this cohort.

Methods

A total of 192 Chilean IBD patients were genotyped using Illumina's Global Screening Array. Genotype data were combined with similar information from 3,147 Chilean controls. The proportions of Aymara, African, European, and Mapuche ancestries were estimated using the software ADMIXTURE. We calculated the odds ratios (ORs) and 95% confidence intervals (CIs) for gender, age, and ancestry proportions. We also explored associations with previously reported IBD-risk variants independently and in conjunction with genetic ancestry.

Results

The first and third quartiles of the proportion of Mapuche ancestry in IBD patients were 24.7 and 34.2%, respectively, and the corresponding OR was 2.30 (95%CI 1.52-3.48) for the lowest vs. the highest group. Only one variant (rs7210086) of the 180 reported IBD-risk SNPs was associated with IBD risk in the Chilean cohort (adjusted P = 0.01). This variant is related to myeloid cells.

Conclusion

The type and proportion of Native American ancestry in Chileans seem to be associated with IBD risk. Variants associated with IBD risk in this Andean region were related to myeloid cells and the innate immune response.

Epigenetic aging differences between Wichí and Criollos from Argentina: Insights from genomic history and ecology

Background and objectives

Epigenetic estimators based on DNA methylation levels have emerged as promising biomarkers of human aging. These estimators exhibit natural variations across human groups, but data about indigenous populations remain underrepresented in research. This study aims to investigate differences in epigenetic estimators between two distinct human populations, both residing in the Gran Chaco region of Argentina, the Native-American Wichí, and admixed Criollos who are descendants of intermarriages between Native Americans and the first European colonizers, using a population genetic approach.

Methodology

We analyzed 24 Wichí (mean age: 39.2 ± 12.9 yo) and 24 Criollos (mean age: 41.1 ± 14.0 yo) for DNA methylation levels using the Infinium MethylationEPIC (Illumina) to calculate 16 epigenetic estimators. Additionally, we examined genome-wide genetic variation using the HumanOmniExpress BeadChip (Illumina) to gain insights into the genetic history of these populations.

Results

Our results indicate that Native-American Wichí are epigenetically older compared to Criollos according to five epigenetic estimators. Analyses within the Criollos population reveal that global ancestry does not influence the differences observed, while local (chromosomal) ancestry shows positive associations between specific SNPs located in genomic regions over-represented by Native-American ancestry and measures of epigenetic age acceleration (AgeAccelHannum). Furthermore, we demonstrate that differences in population ecologies also contribute to observed epigenetic differences.

Conclusions and implications

Overall, our study suggests that while the genomic history may partially account for the observed epigenetic differences, non-genetic factors, such as lifestyle and ecological factors, play a substantial role in the variability of epigenetic estimators, thereby contributing to variations in human epigenetic aging.

Demographic modeling of admixed Latin American populations from whole genomes

Demographic models of Latin American populations often fail to fully capture their complex evolutionary history, which has been shaped by both recent admixture and deeper-in-time demographic events. To address this gap, we used high-coverage whole-genome data from Indigenous American ancestries in present-day Mexico and existing genomes from across Latin America to infer multiple demographic models that capture the impact of different timescales on genetic diversity. Our approach, which combines analyses of allele frequencies and ancestry tract length distributions, represents a significant improvement over current models in predicting patterns of genetic variation in admixed Latin American populations. We jointly modeled the contribution of European, African, East Asian, and Indigenous American ancestries into present-day Latin American populations. We infer that the ancestors of Indigenous Americans and East Asians diverged ∼30 thousand years ago, and we characterize genetic contributions of recent migrations from East and Southeast Asia to Peru and Mexico. Our inferred demographic histories are consistent across different genomic regions and annotations, suggesting that our inferences are robust to the potential effects of linked selection. In conjunction with published distributions of fitness effects for new nonsynonymous mutations in humans, we show in large-scale simulations that our models recover important features of both neutral and deleterious variation. By providing a more realistic framework for understanding the evolutionary history of Latin American populations, our models can help address the historical under-representation of admixed groups in genomics research and can be a valuable resource for future studies of populations with complex admixture and demographic histories.

Native American ancestry and breast cancer risk in Colombian and Mexican women: ruling out potential confounding through ancestry-informative markers

BACKGROUND

Latin American and Hispanic women are less likely to develop breast cancer (BC) than women of European descent. Observational studies have found an inverse relationship between the individual proportion of Native American ancestry and BC risk. Here, we use ancestry-informative markers to rule out potential confounding of this relationship, estimating the confounder-free effect of Native American ancestry on BC risk.

METHODS AND STUDY POPULATION

We used the informativeness for assignment measure to select robust instrumental variables for the individual proportion of Native American ancestry. We then conducted separate Mendelian randomization (MR) analyses based on 1401 Colombian women, most of them from the central Andean regions of Cundinamarca and Huila, and 1366 Mexican women from Mexico City, Monterrey and Veracruz, supplemented by sensitivity and stratified analyses.

RESULTS

The proportion of Colombian Native American ancestry showed a putatively causal protective effect on BC risk (inverse variance-weighted odds ratio [OR] = 0.974 per 1% increase in ancestry proportion, 95% confidence interval [CI] 0.970-0.978, p = 3.1 × 10-40). The corresponding OR for Mexican Native American ancestry was 0.988 (95% CI 0.987-0.990, p = 1.4 × 10-44). Stratified analyses revealed a stronger association between Native American ancestry and familial BC (Colombian women: OR = 0.958, 95% CI 0.952-0.964; Mexican women: OR = 0.973, 95% CI 0.969-0.978), and stronger protective effects on oestrogen receptor (ER)-positive BC than on ER-negative and triple-negative BC.

CONCLUSIONS

The present results point to an unconfounded protective effect of Native American ancestry on BC risk in both Colombian and Mexican women which appears to be stronger for familial and ER-positive BC. These findings provide a rationale for personalised prevention programmes that take genetic ancestry into account, as well as for future admixture mapping studies.

Maternal genetic history of ancient Tibetans over the past 4000 years

The settlement of the Tibetan Plateau epitomizes human adaptation to a high-altitude environment that poses great challenges to human activity. Here, we reconstruct a 4000-year maternal genetic history of Tibetans using 128 ancient mitochondrial genome data from 37 sites in Tibet. The phylogeny of haplotypes M9a1a, M9a1b, D4g2, G2a'c, and D4i show that ancient Tibetans share the most recent common ancestor with ancient Middle and Upper Yellow River populations around the Early and Middle Holocene. In addition, the connections between Tibetans and Northeastern Asians vary over the past 4000 years, with a stronger matrilineal connection between the two during 4000 BP-3000 BP, and a weakened connection after 3000 BP, that are coincident with climate change, followed by a reinforced connection after the Tubo period (1400 BP-1100 BP). Besides, an over 4000-year matrilineal continuity is observed in some of the maternal lineages. We also find the maternal genetic structure of ancient Tibetans is correlated to the geography and interactions between ancient Tibetans and ancient Nepal and Pakistan populations. Overall, the maternal genetic history of Tibetans can be characterized as a long-term matrilineal continuity with frequent internal and external population interactions that are dynamically shaped by geography, climate changes, as well as historical events.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

Testing of two SNP array-based genealogy algorithms using extended Han Chinese pedigrees and recommendations for improved performances in forensic practice

Distant genetic relatives can be linked to a crime scene sample by computing identity-by-state (IBS) and identity-by-descent (IBD) shared by individuals. To test the methods of genetic genealogy estimation and optimal the parameters for forensic investigation, a family-based genetic genealogy analysis was performed using a dataset of 262 Han Chinese individuals from 11 families. The dataset covered relative pairs from 1st- to 14th degrees. But the 7th-degree relative is the most distant kinship to be fully investigated, and each individual has ∼200 relatives within the 7th degree. The KING algorithm by calculating IBS and IBD statistics can correctly discriminate the first-degree relationships of monozygotic twin, parent-offspring and full sibling. The inferred relationship was reliable within the fifth-degree, false positive rate <1.8%. The IBD segment algorithm, GERMLINE + ERSA, could provide reliable inference result prolonged to eighth degree. Analysis of IBD segments produced obviously false negative estimations (<27.4%) rather than false positives (0%) within the eighth-degree inferences. We studied different minimum IBD segment threshold settings (changed from >0 to 6 cM); the inferred results did not make much difference. In distant relative analysis, genetically undetectable relationships begin to occur from the sixth degree (second cousin once removed), which means the offspring after seven meiotic divisions may share no ancestor IBD segment at all. Application of KING and GERMLINE + ERSA worked complementarily to ensure accurate inference from first degree to eighth degree. Using simulated low call rate data, the KING algorithm shows better tolerance to marker decrease compared with the GERMLINE + ERSA segment algorithm.

Methods for Assessing Population Relationships and History Using Genomic Data

Genetic data contain a record of our evolutionary history. The availability of large-scale datasets of human populations from various geographic areas and timescales, coupled with advances in the computational methods to analyze these data, has transformed our ability to use genetic data to learn about our evolutionary past. Here, we review some of the widely used statistical methods to explore and characterize population relationships and history using genomic data. We describe the intuition behind commonly used approaches, their interpretation, and important limitations. For illustration, we apply some of these techniques to genome-wide autosomal data from 929 individuals representing 53 worldwide populations that are part of the Human Genome Diversity Project. Finally, we discuss the new frontiers in genomic methods to learn about population history. In sum, this review highlights the power (and limitations) of DNA to infer features of human evolutionary history, complementing the knowledge gleaned from other disciplines, such as archaeology, anthropology, and linguistics.

Gallbladder Cancer Risk and Indigenous South American Mapuche Ancestry: Instrumental Variable Analysis Using Ancestry-Informative Markers

A strong association between the proportion of indigenous South American Mapuche ancestry and the risk of gallbladder cancer (GBC) has been reported in observational studies. Chileans show the highest incidence of GBC worldwide, and the Mapuche are the largest indigenous people in Chile. We set out to assess the confounding-free effect of the individual proportion of Mapuche ancestry on GBC risk and to investigate the mediating effects of gallstone disease and body mass index (BMI) on this association. Genetic markers of Mapuche ancestry were selected based on the informativeness for assignment measure, and then used as instrumental variables in two-sample Mendelian randomization analyses and complementary sensitivity analyses. Results suggested a putatively causal effect of Mapuche ancestry on GBC risk (inverse variance-weighted (IVW) risk increase of 0.8% per 1% increase in Mapuche ancestry proportion, 95% CI 0.4% to 1.2%, p = 6.7 × 10-5) and also on gallstone disease (3.6% IVW risk increase, 95% CI 3.1% to 4.0%), pointing to a mediating effect of gallstones on the association between Mapuche ancestry and GBC. In contrast, the proportion of Mapuche ancestry showed a negative effect on BMI (IVW estimate -0.006 kg/m2, 95% CI -0.009 to -0.003). The results presented here may have significant implications for GBC prevention and are important for future admixture mapping studies. Given that the association between the individual proportion of Mapuche ancestry and GBC risk previously noted in observational studies appears to be free of confounding, primary and secondary prevention strategies that consider genetic ancestry could be particularly efficient.

Admixture mapping of peripheral artery disease in a Dominican population reveals a putative risk locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have a significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors; however, recent work suggests prevalence may differ between sub-groups. Here, we examined a large cohort of diverse adults in the BioMe biobank in New York City. We observed the prevalence of PAD at 1.7% in EAs vs. 8.5% and 9.4% in AAs and HLs, respectively, and among HL sub-groups, the prevalence was found at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs [OR = 3.15 (95% CI 2.33-4.25), p < 6.44 × 10-14]. To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry and PAD in Dominican BioMe participants (N = 1,813) separately from European, African, and Native American (NAT) continental ancestry tracts. The top association with PAD was an NAT ancestry tract at chromosome 2q35 [OR = 1.96 (SE = 0.16), p < 2.75 × 10-05) with 22.6% vs. 12.9% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607, a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. Efforts to reproduce the signal in other Hispanic cohorts were unsuccessful. In summary, we showed how leveraging health system data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a putative risk locus in a Dominican population.

Genomic history of coastal societies from eastern South America

Sambaqui (shellmound) societies are among the most intriguing archaeological phenomena in pre-colonial South America, extending from approximately 8,000 to 1,000 years before present (yr BP) across 3,000 km on the Atlantic coast. However, little is known about their connection to early Holocene hunter-gatherers, how this may have contributed to different historical pathways and the processes through which late Holocene ceramists came to rule the coast shortly before European contact. To contribute to our understanding of the population history of indigenous societies on the eastern coast of South America, we produced genome-wide data from 34 ancient individuals as early as 10,000 yr BP from four different regions in Brazil. Early Holocene hunter-gatherers were found to lack shared genetic drift among themselves and with later populations from eastern South America, suggesting that they derived from a common radiation and did not contribute substantially to later coastal groups. Our analyses show genetic heterogeneity among contemporaneous Sambaqui groups from the southeastern and southern Brazilian coast, contrary to the similarity expressed in the archaeological record. The complex history of intercultural contact between inland horticulturists and coastal populations becomes genetically evident during the final horizon of Sambaqui societies, from around 2,200 yr BP, corroborating evidence of cultural change.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

The genetic history of the Southern Andes from present-day Mapuche ancestry

The southernmost regions of South America harbor some of the earliest evidence of human presence in the Americas. However, connections with the rest of the continent and the contextualization of present-day indigenous ancestries remain poorly resolved. In this study, we analyze the genetic ancestry of one of the largest indigenous groups in South America: the Mapuche. We generate genome-wide data from 64 participants from three Mapuche populations in Southern Chile: Pehuenche, Lafkenche, and Huilliche. Broadly, we describe three main ancestry blocks with a common origin, which characterize the Southern Cone, the Central Andes, and Amazonia. Within the Southern Cone, ancestors of the Mapuche lineages differentiated from those of the Far South during the Middle Holocene and did not experience further migration waves from the north. We find that the deep genetic split between the Central and Southern Andes is followed by instances of gene flow, which may have accompanied the southward spread of cultural traits from the Central Andes, including crops and loanwords from Quechua into Mapudungun (the language of the Mapuche). Finally, we report close genetic relatedness between the three populations analyzed, with the Huilliche characterized additionally by intense recent exchanges with the Far South. Our findings add new perspectives on the genetic (pre)history of South America, from the first settlement through to the present-day indigenous presence. Follow-up fieldwork took these results back to the indigenous communities to contextualize the genetic narrative alongside indigenous knowledge and perspectives.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

A paleogenome from a Holocene individual supports genetic continuity in Southeast Alaska

Many specifics of the population histories of the Indigenous peoples of North America remain contentious owing to a dearth of physical evidence. Only few ancient human genomes have been recovered from the Pacific Northwest Coast, a region increasingly supported as a coastal migration route for the initial peopling of the Americas. Here, we report paleogenomic data from the remains of a ∼3,000-year-old female individual from Southeast Alaska, named Tatóok yík yées sháawat (TYYS). Our results demonstrate at least 3,000 years of matrilineal genetic continuity in Southeast Alaska, and that TYYS is most closely related to ancient and present-day northern Pacific Northwest Coast Indigenous Americans. We find no evidence of Paleo-Inuit (represented by Saqqaq) ancestry in present-day or ancient Pacific Northwest peoples. Instead, our analyses suggest the Saqqaq genome harbors Northern Native American ancestry. This study sheds further light on the human population history of the northern Pacific Northwest Coast.

Comparing Pruning and Thresholding with Continuous Shrinkage Polygenic Score Methods in a Large Sample of Ancestrally Diverse Adolescents from the ABCD Study®

Using individuals' genetic data researchers can generate Polygenic Scores (PS) that are able to predict risk for diseases, variability in different behaviors as well as anthropomorphic measures. This is achieved by leveraging models learned from previously published large Genome-Wide Association Studies (GWASs) associating locations in the genome with a phenotype of interest. Previous GWASs have predominantly been performed in European ancestry individuals. This is of concern as PS generated in samples with a different ancestry to the original training GWAS have been shown to have lower performance and limited portability, and many efforts are now underway to collect genetic databases on individuals of diverse ancestries. In this study, we compare multiple methods of generating PS, including pruning and thresholding and Bayesian continuous shrinkage models, to determine which of them is best able to overcome these limitations. To do this we use the ABCD Study, a longitudinal cohort with deep phenotyping on individuals of diverse ancestry. We generate PS for anthropometric and psychiatric phenotypes using previously published GWAS summary statistics and examine their performance in three subsamples of ABCD: African ancestry individuals (n = 811), European ancestry Individuals (n = 6703), and admixed ancestry individuals (n = 3664). We find that the single ancestry continuous shrinkage method, PRScs (CS), and the multi ancestry meta method, PRScsx Meta (CSx Meta), show the best performance across ancestries and phenotypes.

Admixture Mapping of Peripheral Artery Disease in a Dominican Population Reveals a Novel Risk Locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors, however recent work suggests prevalence may differ between sub-groups. Here we examined a large cohort of diverse adults in the Bio Me biobank in New York City (NYC). We observed the prevalence of PAD at 1.7% in EAs vs 8.5% and 9.4% in AAs and HLs, respectively; and among HL sub-groups, at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs (OR=3.15 (95% CI 2.33-4.25), P <6.44×10 -14 ). To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry (LA) and PAD in Dominican Bio Me participants (N=1,940) separately for European (EUR), African (AFR) and Native American (NAT) continental ancestry tracts. We identified a NAT ancestry tract at chromosome 2q35 that was significantly associated with PAD (OR=2.05 (95% CI 1.51-2.78), P <4.06×10 -6 ) with 22.5% vs 12.5% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607 , a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. In summary, we showed how leveraging health systems data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a novel risk locus in a Dominican population.

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.

The social lives of isolates (and small language families): the case of the Northwest Amazon

The Americas are home to patches of extraordinary linguistic (genealogical) diversity. These high-diversity areas are particularly unexpected given the recent population of the Americas. In this paper, we zoom in on one such area, the Northwest Amazon, and address the question of how the diversity in this area has persisted to the present. We contrast two hypotheses that claim opposite mechanisms for the maintenance of diversity: the isolation hypothesis suggests that isolation facilitates the preservation of diversity, while the integration hypothesis proposes that conscious identity preservation in combination with contact drives diversity maintenance. We test predictions for both hypotheses across four disciplines: biogeography, cultural anthropology, population genetics and linguistics. Our results show signs of both isolation and integration, but they mainly suggest considerable diversity in how groups of speakers have interacted with their surroundings.

Fast, accurate local ancestry inference with FLARE

Local ancestry is the source ancestry at each point in the genome of an admixed individual. Inferred local ancestry is used for admixture mapping and population genetic analyses. We present FLARE (fast local ancestry estimation), a method for local ancestry inference. FLARE achieves high accuracy through the use of an extended Li and Stephens model, and it achieves exceptional computational performance through incorporation of computational techniques developed for genotype imputation. Memory requirements are reduced through on-the-fly compression of reference haplotypes and stored checkpoints. Computation time is reduced through the use of composite reference haplotypes. These techniques allow FLARE to scale to datasets with hundreds of thousands of sequenced individuals and to provide superior accuracy on large-scale data. FLARE is open source and available at https://github.com/browning-lab/flare.

Bayesian inference of admixture graphs on Native American and Arctic populations

Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans.

Social injustice unveiled by genetic analysis: Argentina as a case study

BACKGROUND

The population of the American countries is genetically heterogeneous, whose genesis result from of recent admixture events. In this process, the transoceanic European component displaced the original inhabitants of the continent.

AIM

To investigate whether socially differentiated cohorts exhibit underlying ancestry components within an urban admixed population, two cohorts of individuals inhabiting Argentina were studied. One cohort included genetically unrelated individuals involved in voluntary paternity testing while the other included sexual or blood-crime suspects.

MATERIALS & METHODS

We analyzed over 2500 unrelated individuals: four Native American maternal lineage mtDNA markers in 1024 samples, five Y chromosome haplogroups in 658 male samples, 24 autosomal ancestry informative markers (AIMs) in 205 samples, and 15 autosomal short tandem repeats (STRs) in 1557 samples; countrywide and divided by regions.

RESULTS

While our results confirm a tricontinental ethnic contribution to both cohorts, their proportions showed statistically significant differences, with a higher proportion of Native American ancestry in the cohort linked to violent crimes compared to those in paternity testing. This hallmark was observed with all the marker sets used and at various levels of analysis.

DISCUSSION

Since paternity tests are costly, socio-economic differences might help to interpret our observations. The effect of discrimination against descendants of Native American minorities, and exposure to violent social environments, might link marginal groups to criminality.

CONCLUSION

Our findings underscore the relevance of proper social management since only by improving living conditions, reducing discrimination, promoting education, and providing job opportunities will it be possible to attain equality in a heterogeneous society. Genetic markers proved to be highly informative in unveiling unexpected social differences.

Bidirectional dispersals during the peopling of the North American Arctic

It is thought that Paleo-Inuit were the first people that settled the American Arctic about 5000 BP (before the present) from a migration that crossed Beringia from Northeast Asia. It is theorized that this group initially migrated to the North Slopes of Alaska and subsequently expanded eastward, eventually reaching Greenland. A second circumpolar dispersal of Neo-Inuit from the North Slopes associated with the Thule-Inuk culture has been postulated to have extended eastward around 800 BP, totally replacing the original Paleo-Inuit without admixing. Although generally accepted, this migration scenario is incompatible with previously reported indications of east to west gene flow across the American Arctic. Here we report on the Y-chromosome haplogroup and Y-STR diversity of the four circumpolar populations of the Tuva Republic (N = 24), Northeast Siberia (N = 9), Bethel, Alaska (N = 40), and Barrow, Alaska (N = 31). Four haplogroup lineages (Q-NWT01, Q-M3, Q-M346, and Q-M120) were detected, Q-NWT01 and Q-M3 being the most abundant at 11.11 and 66.67% in Northeast Siberia, 32.50 and 65.00% in Bethel, and 67.74 and 32.26% in Barrow, respectively. The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STYR loci using the AmpFlSTR Yfiler system. Age estimates and diversity values for the Q-NWT01 and Q-M3 mutations suggest extensive movement of male individuals along the entire longitudinal stretch of the American circumpolar region. Throughout the entire region, Q-M3 exhibits a west to east decreasing gradient in age and diversity while Q-NWT01 indicates the opposite with older TMRCA and higher diversity values running from east to west with the most recent estimates in Canada and Alaska. The high age and diversity values in Greenland are congruent with an origin of the Q-NWT01 mutation in the east of the circumpolar range about 2000-3000 ya. This scenario is incompatible with a complete biological replacement starting about 700 BP of Paleo-Inuit like the Dorset by the Thule-Inuit (Neo-Inuit), as is currently thought, and more parsimonious with gene flow carrying the NWT01 mutation from a pre-Thule population to the ancestors of the present-day Inuit.

Modeling of African population history using f -statistics can be highly biased and is not addressed by previously suggested SNP ascertainment schemes

f -statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. These statistics can provide strong evidence for either admixture or cladality, which can be robust to substantial rates of errors or missing data. f -statistics are guaranteed to be unbiased under "SNP ascertainment" (analyzing non-randomly chosen subsets of single nucleotide polymorphisms) only if it relies on a population that is an outgroup for all groups analyzed. However, ascertainment on a true outgroup that is not co-analyzed with other populations is often impractical and uncommon in the literature. In this study focused on practical rather than theoretical aspects of SNP ascertainment, we show that many non-outgroup ascertainment schemes lead to false rejection of true demographic histories, as well as to failure to reject incorrect models. But the bias introduced by common ascertainments such as the 1240K panel is mostly limited to situations when more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans) or non-human outgroups are co-modelled, for example, f 4 -statistics involving one non-African group, two African groups, and one archaic group. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, cannot fix all these problems since for some classes of f -statistics it is not a clean outgroup ascertainment, and in other cases it demonstrates relatively low power to reject incorrect demographic models since it provides a relatively small number of variants common in anatomically modern humans. And due to the paucity of high-coverage archaic genomes, archaic individuals used for ascertainment often act as sole representatives of the respective groups in an analysis, and we show that this approach is highly problematic. By carrying out large numbers of simulations of diverse demographic histories, we find that bias in inferences based on f -statistics introduced by non-outgroup ascertainment can be minimized if the derived allele frequency spectrum in the population used for ascertainment approaches the spectrum that existed at the root of all groups being co-analyzed. Ascertaining on sites with variants common in a diverse group of African individuals provides a good approximation to such a set of SNPs, addressing the great majority of biases and also retaining high statistical power for studying population history. Such a "pan-African" ascertainment, although not completely problem-free, allows unbiased exploration of demographic models for the widest set of archaic and modern human populations, as compared to the other ascertainment schemes we explored.

Global diversity in individualized cortical network topography

Individualized cortical network topography (ICNT) varies between people and exhibits great variability in the association networks in the human brain. However, these findings were mainly discovered in Western populations. It remains unclear whether and how ICNT is shaped by the non-Western populations. Here, we leveraged a multisession hierarchical Bayesian model to define individualized functional networks in White American and Han Chinese populations with data from both US and Chinese Human Connectome Projects. We found that both the size and spatial topography of individualized functional networks differed between White American and Han Chinese groups, especially in the heteromodal association cortex (including the ventral attention, control, language, dorsal attention, and default mode networks). Employing a support vector machine, we then demonstrated that ethnicity-related ICNT diversity can be used to identify an individual's ethnicity with high accuracy (74%, pperm < 0.0001), with heteromodal networks contributing most to the classification. This finding was further validated through mass-univariate analyses with generalized additive models. Moreover, we reveal that the spatial heterogeneity of ethnic diversity in ICNT correlated with fundamental properties of cortical organization, including evolutionary cortical expansion, brain myelination, and cerebral blood flow. Altogether, this case study highlights a need for more globally diverse and publicly available neuroimaging datasets.

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.

Beringia and the peopling of the Western Hemisphere

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Craniovascular variation in four late Holocene human samples from southern South America

Craniovascular traits in the endocranium (traces of middle meningeal vessels and dural venous sinuses, emissary foramina) provide evidence of vascular anatomy in osteological samples. We investigate the craniovascular variation in four South American samples and the effect of artificial cranial modifications (ACM). CT scans of human adult crania from four archeological samples from southern South America (including skulls with ACM) are used for the analyses. The craniovascular features in the four samples are described, skulls with and without ACM are compared, and additionally, South Americans are compared to a previously analyzed sample of Europeans. Of the four South American samples, the Southern Patagonian differs the most, showing the most distinct cranial dimensions, no ACM, and larger diameters of the emissary foramina. Unlike previous studies, we did not find any major differences in craniovascular features between modified and non-modified skulls, except that the skulls with ACM present somewhat smaller foramina. South Americans significantly differed from Europeans, especially in the anteroposterior dominance of the middle meningeal artery, in the pattern of sinus confluence, in the occurrence of enlarged occipito-marginal sinuses, and in foramina frequencies and diameters. Craniovascular morphology is not affected by the cranial size, even in skulls with ACM, indicating a minor or null influence of structural topological factors. Concerning the samples from distinct geographic and climatic environments, it must be evaluated whether the craniovascular morphogenesis might be partially influenced by specific functions possibly associated with thermoregulation, intracranial pressure, and the maintenance of intracranial homeostasis.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing D_ST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.