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

Using uniparental genetic profiles to unravel the complexity of Argentine admixed populations

Latin American countries are distinguished by their highly admixed populations, characterized by a significant preservation of Native American matrilineal ancestry. This contrasts with the paternal lineages, which exhibit different patterns due to pronounced sex-biased mating practices during the colonial period. Uniparental genetic markers have been instrumental in population genetics, facilitating the reconstruction of human settlement histories and serving forensic identification purposes. The primary objective of this study was to investigate the diversity and structure of lineage markers in Argentina and compare them with other admixed populations in South America. For this study, we analyzed Y-STR and mtDNA haplotypes from 5202 unrelated individuals, providing a detailed description of the observed variability in both markers. Additionally, we conducted a genetic distance analysis, incorporating data from bibliographic sources across Argentina and South America. In pairwise comparisons among provinces, higher FST values were found in mtDNA haplotypes than in Y-STR haplotypes. This allows for more provinces to be grouped by similarity when using Y-STR data. These differences were also evident in the multidimensional scaling (MDS) analysis between South American countries. Y-STR haplotypes showed greater similarity to European haplotypes, whereas mtDNA haplotypes exhibited greater dispersion. Thus, the comprehensive compilation of haplotypes in this study, including those integrated from our research and those cited in existing literature, provides an in-depth understanding of the inherent genetic complexities within Argentina.

A genomic history of the North Pontic Region from the Neolithic to the Bronze Age

The North Pontic Region was the meeting point of the farmers of Old Europe and the foragers and pastoralists of the Eurasian steppe1,2, and the source of migrations deep into Europe3-5. Here we report genome-wide data from 81 prehistoric North Pontic individuals to understand the genetic makeup of its people. North Pontic foragers had ancestry from Balkan and Eastern hunter-gatherers6 as well as European farmers and, occasionally, Caucasus hunter-gatherers. During the Eneolithic period, a wave of migrants from the Caucasus-Lower Volga area7 bypassed local foragers to mix in equal parts with Trypillian farmers, forming the people of the Usatove culture around 4500 BCE. A temporally overlapping wave of migrants from the Caucasus-Lower Volga blended with foragers instead of farmers to form Serednii Stih people7. The third wave was the Yamna-descendants of the Serednii Stih who formed by mixture around 4000 BCE and expanded during the Early Bronze Age (3300 BCE). The temporal gap between Serednii Stih and the Yamna is bridged by a genetically Yamna individual from Mykhailivka, Ukraine (3635-3383 BCE), a site of archaeological continuity across the Eneolithic-Bronze Age transition and a likely epicentre of Yamna formation. Each of these three waves of migration propagated distinctive ancestries while also incorporating outsiders, a flexible strategy that may explain the success of the peoples of the North Pontic in spreading their genes and culture across Eurasia3-5,8-10.

The genetic origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300 BC across the steppes north of the Black and Caspian Seas, and by 3000 BC it reached its maximal extent, ranging from Hungary in the west to Kazakhstan in the east. To localize Yamnaya origins among the preceding Eneolithic people, we assembled ancient DNA from 435 individuals, demonstrating three genetic clines. A Caucasus-lower Volga (CLV) cline suffused with Caucasus hunter-gatherer1 ancestry extended between a Caucasus Neolithic southern end and a northern end at Berezhnovka along the lower Volga river. Bidirectional gene flow created intermediate populations, such as the north Caucasus Maikop people, and those at Remontnoye on the steppe. The Volga cline was formed as CLV people mixed with upriver populations of Eastern hunter-gatherer2 ancestry, creating hypervariable groups, including one at Khvalynsk. The Dnipro cline was formed when CLV people moved west, mixing with people with Ukraine Neolithic hunter-gatherer ancestry3 along the Dnipro and Don rivers to establish Serednii Stih groups, from whom Yamnaya ancestors formed around 4000 BC and grew rapidly after 3750-3350 BC. The CLV people contributed around four-fifths of the ancestry of the Yamnaya and, entering Anatolia, probably from the east, at least one-tenth of the ancestry of Bronze Age central Anatolians, who spoke Hittite4,5. We therefore propose that the final unity of the speakers of 'proto-Indo-Anatolian', the language ancestral to both Anatolian and Indo-European people, occurred in CLV people some time between 4400 BC and 4000 BC.

Genetic Ancestry and 21-Gene Oncotype DX Breast Cancer Recurrence Scores

BACKGROUND

Racial and ethnic breast cancer disparities persist. This may be reflected by differences in Oncotype DX recurrence scores (RS), which are higher for Black women. This study assesses the association between ancestry and RS.

METHODS

Stage I-III hormone receptor (HR)-positive, human epidermal growth factor receptor 2 (HER2)-negative breast cancer patients with ancestry and RS data were prospectively identified from 2017 to 2021. RS were grouped into low, intermediate, and high categories. Multinomial regression determined the association between ancestry and RS controlling for ancestry and estrogen receptor (ER), progesterone receptor (PR), and HER2 expression.

RESULTS

Of 174 patients, 28 (16.1%) self-identified as non-Hispanic White, 107 (61.5%) self-identified as Hispanic White, 19 (10.9%) self-identified as non-Hispanic Black, and 9 (5.2%) self-identified as Hispanic Black. Ninety-four (54.0%) patients had low RS, 51 (29.3%) had intermediate RS, and 29 (16.7%) had high RS. On multivariable analyses, West African ancestry was associated with increased odds of intermediate (odds ratio [OR] 1.02, 95% confidence interval [CI] 1.00-1.04, p = 0.039) and high (OR 1.03, 95% CI 1.00-1.06, p = 0.022) RS. East Asian ancestry was associated with decreased odds of intermediate RS (OR 0.78, 95% CI 0.60-1.00, p = 0.048). Increasing ER (OR 0.43, 95% CI 0.23-0.82, p = 0.011), PR (OR 0.20, 95% CI 0.11-0.34, p < 0.001), and HER2 (OR 0.24, 95% CI 0.09-0.63, p = 0.004) expression were associated with lower odds of high RS.

CONCLUSION

Increasing West African ancestry is associated with increased odds of high and intermediate RS, while increasing East Asian ancestry is associated with lower odds of intermediate RS. These findings require validation but suggest ancestry may represent a biological biomarker and that assays guiding adjuvant therapy may require ancestry-based calibration in HR+/HER2- tumors.

The impact of Indigenous American-like ancestry on risk of acute lymphoblastic leukemia in Hispanic/Latino children

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer, with Hispanic/Latino children having a higher incidence of ALL than other racial/ethnic groups. Genetic variants, particularly ones found enriched in Indigenous American (IA)-like ancestry and inherited by Hispanics/Latinos, may contribute to this disparity. In this study, we characterized the impact of IA-like ancestry on overall ALL risk and the frequency and effect size of known risk alleles in a large cohort of self-reported Hispanic/Latino individuals. We also performed genome-wide admixture mapping analysis to identify potentially novel ALL risk loci. We found that global IA ancestry was positively associated with ALL risk, but the association was not significant after adjusting for socio-economic indicators. In a series of local ancestry analyses, we uncovered that at known ALL risk loci, increasing copies of the IA-like haplotype were positively and significantly associated with ALL case-control status. Further, the IA-like haplotype had ~1.33 times the odds of harboring the risk allele compared to non-IA-like haplotypes. We found no evidence of interaction between genotype and ancestry (local or global) in relation to ALL risk. Admixture mapping identified association signals on chromosomes 2 (2q21.2), 7 (7p12.2), 10 (10q21.2), and 15 (15q22.31); however, only the variants at 7p12.2 and 10q21.2 replicated in additional cohorts. Taken together, our results suggest that increased risk of ALL in Hispanic/Latino children may be conferred by higher frequency of risk alleles within IA-like ancestry, which can be leveraged as targets of new precision health strategies and therapeutics.

The genomic origin of early maize in eastern North America

Indigenous maize varieties from eastern North America have played an outsized role in breeding programs, yet their early origins are not fully understood. We generated paleogenomic data to reconstruct how maize first reached this region and how it was selected during the process. Genomic ancestry analyses reveal recurrent movements northward from different parts of Mexico, likely culminating in at least two dispersals from the US Southwest across the Great Plains to the Ozarks and beyond. We find that 1,000-year-old Ozark specimens carry a highly differentiated wx1 gene, which is involved in the synthesis of amylose, highlighting repeated selective pressures on the starch metabolic pathway throughout maize's domestication. This population shows a close affinity with the lineage that ultimately became the Northern Flints, a major contributor to modern commercial maize.

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.

High-resolution genomic history of early medieval Europe

Many known and unknown historical events have remained below detection thresholds of genetic studies because subtle ancestry changes are challenging to reconstruct. Methods based on shared haplotypes1,2 and rare variants3,4 improve power but are not explicitly temporal and have not been possible to adopt in unbiased ancestry models. Here we develop Twigstats, an approach of time-stratified ancestry analysis that can improve statistical power by an order of magnitude by focusing on coalescences in recent times, while remaining unbiased by population-specific drift. We apply this framework to 1,556 available ancient whole genomes from Europe in the historical period. We are able to model individual-level ancestry using preceding genomes to provide high resolution. During the first half of the first millennium CE, we observe at least two different streams of Scandinavian-related ancestry expanding across western, central and eastern Europe. By contrast, during the second half of the first millennium CE, ancestry patterns suggest the regional disappearance or substantial admixture of these ancestries. In Scandinavia, we document a major ancestry influx by approximately 800 CE, when a large proportion of Viking Age individuals carried ancestry from groups related to central Europe not seen in individuals from the early Iron Age. Our findings suggest that time-stratified ancestry analysis can provide a higher-resolution lens for genetic history.

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.

Inbreeding and Gallbladder Cancer Risk: Homozygosity Associations Adjusted for Indigenous American Ancestry, BMI, and Genetic Risk of Gallstone Disease

Latin Americans have a rich genetic make-up that translates into heterogeneous fractions of the autosomal genome in runs of homozygosity (FROH) and heterogeneous types and proportions of indigenous American ancestry. While autozygosity has been linked to several human diseases, very little is known about the relationship between inbreeding, genetic ancestry, and cancer risk in Latin Americans. Chile has one of the highest incidences of gallbladder cancer (GBC) in the world, and we investigated the association between inbreeding, GBC, gallstone disease (GSD), and body mass index (BMI) in 4029 genetically admixed Chileans. We calculated individual FROH above 1.5 Mb and weighted polygenic risk scores for GSD, and applied multiple logistic regression to assess the association between homozygosity and GBC risk. We found that homozygosity was due to a heterogeneous mixture of genetic drift and consanguinity in the study population. Although we found no association between homozygosity and overall GBC risk, we detected interactions of FROH with sex, age, and genetic risk of GSD that affected GBC risk. Specifically, the increase in GBC risk per 1% FROH was 19% in men (p-value = 0.002), 30% in those under 60 years of age (p-value = 0.001), and 12% in those with a genetic risk of GSD above the median (p-value = 0.01). The present study highlighted the complex interplay between inbreeding, genetic ancestry, and genetic risk of GSD in the development of GBC. The applied methodology and our findings underscored the importance of considering the population-specific genetic architecture, along with sex- and age-specific effects, when investigating the genetic basis of complex traits in Latin Americans.

Hemochromatosis neural archetype reveals iron disruption in motor circuits

Our understanding of brain iron regulation and its disruption in disease is limited. Excess iron affects motor circuitry, contributing to Parkinson's disease (PD) risk. The molecular mechanisms regulating central iron levels, beyond a few well-known genes controlling peripheral iron, remain unclear. We generated scores based on the archetypal brain iron accumulation observed in magnetic resonance imaging scans of individuals with excessive dietary iron absorption and hemochromatosis risk. Genome-wide analysis revealed that this score is highly heritable, identifying loci associated with iron homeostasis, and driven by peripheral iron levels. Our score predicted gait abnormalities and showed a U-shaped relationship with PD risk, identifying individuals with threefold increased risk. These results establish a hormetic relationship between brain iron and PD risk, where central iron levels are strongly determined by genetics via peripheral iron. This framework combining forward and reverse genetics is a powerful study design to understand genomic drivers underlying high dimensional phenotypes.

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.

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.

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.

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.

Bioarchaeological Perspectives on Late Antiquity in Dalmatia: Paleogenetic, Dietary, and Population Studies of the Hvar - Radošević burial site

Late Antiquity Dalmatia was a time and place of political unrest in the Roman Empire that influenced the lives of those in that region. The Late Antique burial site of Hvar - Radošević, spanning the 3rd to 5th centuries CE, is located on the Croatian Dalmatian island of Hvar. Given the time frame and location on a busy marine trade route, the study of this burial site offers us a glimpse into the lives of the Late Antique population living on this island. It comprises 33 individuals, with 17 buried within a confined grave tomb and the remaining individuals buried in separate locations in the tomb's proximity. The study aims to provide a new perspective on the lives of people on the island during those times by studying ancestry, population structure, possible differences within the buried population, dietary habits, and general health. The genetic analysis of the ancestral origins of the individuals buried at Hvar - Radošević revealed a diverse population reflective of the era's genetic variability. The identification of genetic outliers suggests a range of ancestries from distinct regions of the Roman Empire, possibly linked to trade routes associated with the Late Antique port in ancient Hvar. Stable isotope ratio analysis (δ13C and δ15N) indicated a diet mainly consisting of C3 plants, with minimal consumption of marine foods. High childhood mortality rates, physiological stress markers, and dental diseases suggest a low quality of life in the population. Assessment of kinship and dietary patterns revealed no discernible distinctions between individuals buried within the tomb and those buried outside, indicative of an absence of differential burial practices based on social status and familial ties among this specific buried population.

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.

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.