Uveal melanoma (UM) is a rare but frequently metastasizing cancer. Genome-wide association studies have identified three common genome-wide significant germline risk loci. Here, we perform a genome-wide association study on 401 new cases and conduct a meta-analysis with three independent previously published cohorts for a total sample size of 2,426 cases. We confirm the three previously identified risk loci and identify four additional genome-wide significant loci. We find that eye pigmentation-decreasing variants are systematically associated with increased UM risk and that selection for lighter pigmentation in the past 5,000 years explains about 73% of the difference in UM incidence between Northern and Southern Europe. We find evidence of effect size heterogeneity at significant loci across cohorts, in particular, a weaker association between eye pigmentation and UM in a Finnish cohort. Finally, we confirm differential effect sizes between uveal melanoma cases with and without loss of chromosome 3, the major determinant of metastatic risk. Our study identifies novel germline risk factors for UM and highlights genetic and environmental heterogeneity in its etiology.

• Keywords
• eye pigmentation, genome-wide association study, uveal melanoma,
• MeSH
• Genome-Wide Association Study * MeSH
• Genetic Heterogeneity * MeSH
• Genetic Predisposition to Disease * MeSH
• Genetic Loci * MeSH
• Polymorphism, Single Nucleotide MeSH
• Humans MeSH
• Melanoma * genetics   epidemiology   pathology   MeSH
• Uveal Neoplasms * genetics   epidemiology   pathology   MeSH
• Risk Factors MeSH
• Uveal Melanoma MeSH
• Check Tag
• Humans MeSH
• Female MeSH
• Publication type
• Journal Article MeSH
• Meta-Analysis MeSH

BACKGROUND: Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes. RESULTS: We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes. CONCLUSIONS: This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

• Keywords
• Flanders, Low countries, Medieval, Migration, Palaeo-genomics, Plague, Urbanization,
• MeSH
• History, Medieval MeSH
• Genetic Variation MeSH
• Genome, Human MeSH
• Genomics MeSH
• Humans MeSH
• Plague epidemiology   history   genetics   MeSH
• Genetics, Population MeSH
• Urbanization * history   MeSH
• Check Tag
• History, Medieval MeSH
• Humans MeSH
• Publication type
• Journal Article MeSH
• Historical Article MeSH
• Geographicals
• Belgium MeSH

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on 2D stepping stone landscapes and in situations with low prestudy odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture graph-shaped and stepping stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: (1) prestudy odds are low and fall rapidly with increasing model complexity; (2) complex migration networks violate the assumptions of the method; hence, there is poor correlation between qpAdm P-values and model optimality, contributing to low but nonzero false-positive rate and low power; and (3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly nonoptimal models have estimates in the same interval, contributing to the false-positive rate. We also reinterpret large sets of qpAdm models from 2 studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: (1) temporal stratification of targets and proxy sources in the case of admixture graph-shaped histories, (2) focused exploration of few models for increasing prestudy odds; and (3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false-positive rate.

• Keywords
• qpAdm, admixture graphs, archaeogenetics, genetic admixture, simulation, stepping stone models,
• MeSH
• Humans MeSH
• Models, Genetic * MeSH
• Genetics, Population * methods   MeSH
• Software * MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH

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.

• MeSH
• White People genetics   history   MeSH
• History, Ancient MeSH
• European People MeSH
• Phylogeny MeSH
• Haplotypes MeSH
• Humans MeSH
• Human Migration * history   MeSH
• DNA, Mitochondrial analysis   MeSH
• Genetics, Population * MeSH
• Rivers MeSH
• DNA, Ancient analysis   MeSH
• Gene Flow MeSH
• Check Tag
• History, Ancient MeSH
• Humans MeSH
• Publication type
• Journal Article MeSH
• Historical Article MeSH
• Names of Substances
• DNA, Mitochondrial MeSH
• DNA, Ancient MeSH

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on two-dimensional stepping-stone landscapes and in situations with low pre-study odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture-graph-shaped and stepping-stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: 1) pre-study odds are low and fall rapidly with increasing model complexity; 2) complex migration networks violate the assumptions of the method, hence there is poor correlation between qpAdm p-values and model optimality, contributing to low but non-zero false positive rate and low power; 3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly non-optimal models have estimates in the same interval, contributing to the false positive rate. We also re-interpret large sets of qpAdm models from two studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: 1) temporal stratification of targets and proxy sources in the case of admixture-graph-shaped histories; 2) focused exploration of few models for increasing pre-study odds; 3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false positive rate.

• Keywords
• admixture graphs, archaeogenetics, genetic admixture, qpAdm, simulation, stepping-stone models,
• Publication type
• Journal Article MeSH
• Preprint MeSH

The Linearbandkeramik (LBK) Neolithic communities were the first to spread farming across large parts of Europe. We report genome-wide data for 250 individuals: 178 individuals from whole-cemetery surveys of the Alföld Linearbankeramik Culture eastern LBK site of Polgár-Ferenci-hát, the western LBK site of Nitra Horné Krškany and the western LBK settlement and massacre site of Asparn-Schletz, as well as 48 LBK individuals from 16 other sites and 24 earlier Körös and Starčevo individuals from 17 more sites. Here we show a systematically higher percentage of western hunter-gatherer ancestry in eastern than in western LBK sites, showing that these two distinct LBK groups had different genetic trajectories. We find evidence for patrilocality, with more structure across sites in the male than in the female lines and a higher rate of within-site relatives for males. At Asparn-Schletz we find almost no relatives, showing that the massacred individuals were from a large population, not a small community.

• MeSH
• History, Ancient MeSH
• Genetic Variation * MeSH
• Humans MeSH
• Farmers * history   MeSH
• Agriculture * history   MeSH
• Check Tag
• History, Ancient MeSH
• Humans MeSH
• Male MeSH
• Female MeSH
• Publication type
• Journal Article MeSH
• Historical Article MeSH
• Geographicals
• Europe MeSH

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.

• Keywords
• f-statistics, aDNA, admixture, ancient DNA, archaeogenetics, paleogenomics, qpAdm,
• MeSH
• Demography MeSH
• Genomics * methods   MeSH
• Models, Genetic MeSH
• Paleontology * methods   MeSH
• Software MeSH
• Data Accuracy MeSH
• Publication type
• Journal Article MeSH

During the Hungarian Conquest in the 10th century CE, the early medieval Magyars, a group of mounted warriors from Eastern Europe, settled in the Carpathian Basin. They likely introduced the Hungarian language to this new settlement area, during an event documented by both written sources and archaeological evidence. Previous archaeogenetic research identified the newcomers as migrants from the Eurasian steppe. However, genome-wide ancient DNA from putative source populations has not been available to test alternative theories of their precise source. We generated genome-wide ancient DNA data for 131 individuals from candidate archaeological contexts in the Circum-Uralic region in present-day Russia. Our results tightly link the Magyars to people of the Early Medieval Karayakupovo archaeological horizon on both the European and Asian sides of the southern Urals. Our analyes show that ancestors of the people of the Karayakupovo archaeological horizon were established in the Southern Urals by the Iron Age and that their descendants persisted locally in the Volga-Kama region until at least the 14th century.

• Publication type
• Journal Article MeSH
• Preprint MeSH

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.

• Publication type
• Journal Article MeSH
• Preprint MeSH

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.

• Keywords
• aDNA, admixture, archaeogenetics, f-statistics, paleogenomics, qpAdm,
• Publication type
• Journal Article MeSH
• Preprint MeSH