Benchmarking the Mantel test and derived methods for testing association between distance matrices.

Testing the association between objects is central in ecology, evolution, and quantitative sciences in general. Two types of variables can describe the relationships between objects: point variables (measured on individual objects), and distance variables (measured between pairs of objects). The Mantel test and derived methods have been extensively used for distance variables. Yet, these methods have been criticized due to low statistical power and inflated type I error when spatial autocorrelation is present. Here, we assessed the statistical power between different types of tested variables and the type I error rate over a wider range of autocorrelation intensities than previously assessed, both on univariate and multivariate data. We also illustrated the performance of distance matrix statistics through computational simulations of genetic diversity. We show that the Mantel test and derived methods are not affected by inflated type I error when spatial autocorrelation affects only one variable when investigating correlations, or when either the response or the explanatory variable(s) is affected by spatial autocorrelation while investigating causal relationships. As previously noted, with autocorrelation affecting more variables, inflated type I error could be reduced by modifying the significance threshold. Additionally, the Mantel test has no problem of statistical power when the hypothesis is formulated in terms of distance variables. We highlight that transformation of variable types should be avoided because of the potential information loss and modification of the tested hypothesis. We propose a set of guidelines to help choose the appropriate method according to the type of variables and defined hypothesis.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses.

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

Past human expansions shaped the spatial pattern of Neanderthal ancestry.

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

Ancient mitochondrial diversity reveals population homogeneity in Neolithic Greece and identifies population dynamics along the Danubian expansion axis.

The aim of the study is to investigate mitochondrial diversity in Neolithic Greece and its relation to hunter-gatherers and farmers who populated the Danubian Neolithic expansion axis. We sequenced 42 mitochondrial palaeogenomes from Greece and analysed them together with European set of 328 mtDNA sequences dating from the Early to the Final Neolithic and 319 modern sequences. To test for population continuity through time in Greece, we use an original structured population continuity test that simulates DNA from different periods by explicitly considering the spatial and temporal dynamics of populations. We explore specific scenarios of the mode and tempo of the European Neolithic expansion along the Danubian axis applying spatially explicit simulations coupled with Approximate Bayesian Computation. We observe a striking genetic homogeneity for the maternal line throughout the Neolithic in Greece whereas population continuity is rejected between the Neolithic and present-day Greeks. Along the Danubian expansion axis, our best-fitting scenario supports a substantial decrease in mobility and an increasing local hunter-gatherer contribution to the gene-pool of farmers following the initial rapid Neolithic expansion. Οur original simulation approach models key demographic parameters rather than inferring them from fragmentary data leading to a better understanding of this important process in European prehistory.

Challenging Ancient DNA Results About Putative HLA Protection or Susceptibility to Yersinia pestis.

In a recent article, Immel et al. (Immel A, Key FM, Szolek A, Barquera R, Robinson MK, Harrison GF, Palmer WH, Spyrou MA, Susat J, Krause-Kyora B, et al. 2021. Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes. Mol Biol Evol. 38:4059-4076) extracted DNA from 36 individuals dead from plague in Ellwangen, Southern Germany, during the 16th century. By comparing their human leukocyte antigen (HLA) genotypes with those of 50 present-day Ellwangen inhabitants, the authors reported a significant decrease of HLA-B*51:01 and HLA-C*06:02 and a significant increase of HLA-DRB1*13:01/13:02 frequencies from ancient to modern populations. After comparing these frequencies with a larger sample of 8,862 modern Germans and performing simulations of natural selection, they concluded that these changes had been driven by natural selection. In an attempt to provide more evidence on such stimulating results, we explored the HLA frequency patterns over all of Europe, we predicted binding affinities of HLA-B/C/DRB1 alleles to 106,515 Yersinia pestis-derived peptides, and we performed forward simulations of HLA genetic profiles under neutrality. Our analyses do not sustain the conclusions of HLA protection or susceptibility to plague based on ancient DNA.

Simulated patterns of mitochondrial diversity are consistent with partial population turnover in Bronze Age Central Europe.

OBJECTIVES: The analysis of ancient mitochondrial DNA from osteological remains has challenged previous conclusions drawn from the analysis of mitochondrial DNA from present populations, notably by revealing an absence of genetic continuity between the Neolithic and modern populations in Central Europe. Our study investigates how to reconcile these contradictions at the mitochondrial level using a modeling approach. MATERIALS AND METHODS: We used a spatially explicit computational framework to simulate ancient and modern DNA sequences under various evolutionary scenarios of post Neolithic demographic events and compared the genetic diversity of the simulated and observed mitochondrial sequences. We investigated which-if any-scenarios were able to reproduce statistics of genetic diversity similar to those observed, with a focus on the haplogroup N1a, associated with the spread of early Neolithic farmers. RESULTS: Demographic fluctuations during the Neolithic transition or subsequent demographic collapses after this period, that is, due to epidemics such as plague, are not sufficient to explain the signal of population discontinuity detected on the mitochondrial DNA in Central Europe. Only a scenario involving a substantial genetic input due to the arrival of migrants after the Neolithic transition, possibly during the Bronze Age, is compatible with observed patterns of genetic diversity. DISCUSSION: Our results corroborate paleogenomic studies, since out of the alternative hypotheses tested, the best one that was able to recover observed patterns of mitochondrial diversity in modern and ancient Central European populations was one were immigration of populations from the Pontic steppes during the Bronze Age was explicitly simulated.

Spatially explicit paleogenomic simulations support cohabitation with limited admixture between Bronze Age Central European populations.

The Bronze Age is a complex period of social, cultural and economic changes. Recent paleogenomic studies have documented a large and rapid genetic change in early Bronze Age populations from Central Europe. However, the detailed demographic and genetic processes involved in this change are still debated. Here we have used spatially explicit simulations of genomic components to better characterize the demographic and migratory conditions that may have led to this change. We investigated various scenarios representing the expansion of pastoralists from the Pontic steppe, potentially linked to the Yamnaya cultural complex, and their interactions with local populations in Central Europe, considering various eco-evolutionary factors, such as population admixture, competition and long-distance dispersal. Our results do not support direct competition but rather the cohabitation of pastoralists and farmers in Central Europe, with limited gene flow between populations. They also suggest occasional long-distance migrations accompanying the expansion of pastoralists and a demographic decline in both populations following their initial contact. These results link recent archaeological and paleogenomic observations and move further the debate of genomic changes during the early Bronze Age.

Air temperature influences early Covid-19 outbreak as indicated by worldwide mortality.

The Covid-19 outbreak has triggered a global crisis that is challenging governments, health systems and the scientific community worldwide. A central question in the Covid-19 pandemic is whether climatic factors have influenced its progression. To address this question, we used mortality rates during the first three weeks of recorded mortality in 144 countries, during the first wave of the pandemic. We examined the effect of climatic variables, along with the proportion of the population older than 64 years old, the number of beds in hospitals, and the timing and strength of the governmental travel measures to control the spread of the disease. Our first model focuses on air temperature as the central climatic factor and explains 67% of the variation in mortality rate, with 37% explained by the fixed variables considered and 31% explained by country-specific variations. We show that mortality rate is negatively influenced by warmer air temperature. Each additional Celsius degree decreases mortality rate by ~5%. Our second model is centred on the UV Index and follows the same trend as air temperature, explaining 69% of the variation in mortality rate. These results are robust to the exclusion of countries with low incomes, as well as to the exclusion of low- and medium-income countries. We also show that the proportion of vulnerable age classes and access to healthcare are critical factors impacting the mortality rate of this disease. The effects of air temperature at an early stage of the Covid-19 outbreak is a key factor to understand the primary spread of this pandemic, and should be considered in projecting subsequent waves.

Projecting introgression from domestic cats into European wildcats in the Swiss Jura.

Hybridization between wild and domesticated organisms is a worldwide conservation issue. In the Jura Mountains, threatened European wildcats (Felis silvestris) have been demographically spreading for approximately the last 50 years, but this recovery is coupled with hybridization with domestic cats (Felis catus). Here, we project the pattern of future introgression using different spatially explicit scenarios to model the interactions between the two species, including competition and different population sizes. We project the fast introgression of domestic cat genes into the wildcat population under all scenarios if hybridization is not severely restricted. If the current hybridization rate and population sizes remain unchanged, we expect the loss of genetic distinctiveness between wild and domestic cats at neutral nuclear, mitochondrial and Y chromosome markers in one hundred years. However, scenarios involving a competitive advantage for wildcats and a future increase in the wildcat population size project a slower increase in introgression. We recommend that future studies assess the fitness of these hybrids and better characterize their ecological niche and their ecological interactions with parental species to elucidate effective conservation measures.

Harmonizing hybridization dissonance in conservation.

A dramatic increase in the hybridization between historically allopatric species has been induced by human activities. However, the notion of hybridization seems to lack consistency in two respects. On the one hand, it is inconsistent with the biological species concept, which does not allow for interbreeding between species, and on the other hand, it is considered either as an evolutionary process leading to the emergence of new biodiversity or as a cause of biodiversity loss, with conservation implications. In the first case, we argue that conservation biology should avoid the discussion around the species concept and delimit priorities of conservation units based on the impact on biodiversity if taxa are lost. In the second case, we show that this is not a paradox but an intrinsic property of hybridization, which should be considered in conservation programmes. We propose a novel view of conservation guidelines, in which human-induced hybridization may also be a tool to enhance the likelihood of adaptation to changing environmental conditions or to increase the genetic diversity of taxa affected by inbreeding depression. The conservation guidelines presented here represent a guide for the development of programmes aimed at protecting biodiversity as a dynamic evolutionary system.

Influence of Paleolithic range contraction, admixture and long-distance dispersal on genetic gradients of modern humans in Asia.

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long-distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

SPLATCHE3: simulation of serial genetic data under spatially explicit evolutionary scenarios including long-distance dispersal.

SUMMARY: SPLATCHE3 simulates genetic data under a variety of spatially explicit evolutionary scenarios, extending previous versions of the framework. The new capabilities include long-distance migration, spatially and temporally heterogeneous short-scale migrations, alternative hybridization models, simulation of serial samples of genetic data and a large variety of DNA mutation models. These implementations have been applied independently to various studies, but grouped together in the current version. AVAILABILITY AND IMPLEMENTATION: SPLATCHE3 is written in C++ and is freely available for non-commercial use from the website http://www.splatche.com/splatche3. It includes console versions for Linux, MacOs and Windows and a user-friendly GUI for Windows, as well as detailed documentation and ready-to-use examples.

Hybridization and introgression during density-dependent range expansion: European wildcats as a case study.

Interbreeding between historically allopatric species with incomplete reproductive barriers may result when species expand their range. The genetic consequences of such hybridization depend critically on the dynamics of the range expansion. Hybridization models during range expansion have been developed but assume dispersal to be independent from neighboring population densities. However, organisms may disperse because they are attracted by conspecifics or because they prefer depopulated areas. Here, through spatially explicit simulations, we assess the effect of various density-dependent dispersal modes on the introgression between two species. We find huge introgression from the local species into the invasive one with all dispersal modes investigated, even when the hybridization rate is relatively low. This represents a general expectation for neutral genes even if the dispersal modes differ in colonization times and amount of introgression. Invasive individuals attracted by conspecifics need more time to colonize the whole area and are more introgressed by local genes, whereas the opposite is found for solitary individuals. We applied our approach to a recent expansion of European wildcats in the Jura Mountains and the hybridization with domestic cats. We show that the simulations explained better the observed level of introgression at nuclear, mtDNA, and Y chromosome markers, when using solitary dispersal for wildcats instead of random or gregarious dispersal, in accordance with ecological knowledge. Thus, use of density-dependent dispersal models increases the predictive power of the approach.

Bayesian estimation of partial population continuity using ancient DNA and spatially explicit simulations.

The retrieval of ancient DNA from osteological material provides direct evidence of human genetic diversity in the past. Ancient DNA samples are often used to investigate whether there was population continuity in the settlement history of an area. Methods based on the serial coalescent algorithm have been developed to test whether the population continuity hypothesis can be statistically rejected by analysing DNA samples from the same region but of different ages. Rejection of this hypothesis is indicative of a large genetic shift, possibly due to immigration occurring between two sampling times. However, this approach is only able to reject a model of full continuity model (a total absence of genetic input from outside), but admixture between local and immigrant populations may lead to partial continuity. We have recently developed a method to test for population continuity that explicitly considers the spatial and temporal dynamics of populations. Here, we extended this approach to estimate the proportion of genetic continuity between two populations, using ancient genetic samples. We applied our original approach to the question of the Neolithic transition in Central Europe. Our results confirmed the rejection of full continuity, but our approach represents an important step forward by estimating the relative contribution of immigrant farmers and of local hunter-gatherers to the final Central European Neolithic genetic pool. Furthermore, we show that a substantial proportion of genes brought by the farmers in this region were assimilated from other hunter-gatherer populations along the way from Anatolia, which was not detectable by previous continuity tests. Our approach is also able to jointly estimate demographic parameters, as we show here by finding both low density and low migration rate for pre-Neolithic hunter-gatherers. It provides a useful tool for the analysis of the numerous ancient DNA data sets that are currently being produced for many different species.

Consequences of diverse evolutionary processes on american genetic gradients of modern humans.

European genetic gradients of modern humans were initially interpreted as a consequence of the demic diffusion of expanding Neolithic farmers. However, recent studies showed that these gradients may also be influenced by other evolutionary processes such as population admixture or range contractions. Genetic gradients were observed in the Americas, although their specific evolutionary causes were not investigated. Here we extended the approach used to study genetic gradients in Europe to analyze the influence of diverse evolutionary scenarios on American genetic gradients. Using extensive computer simulations, we evaluated the impact of (i) admixture between expansion waves of modern humans, (ii) the presence of ice-sheets during the last glacial maximum (LGM) and (iii) long-distance dispersal (LDD) events, on the genetic gradients (detected by principal component analysis) of the entire continent, North America and South America. The specific simulation of North and South America showed that genetic gradients are usually orthogonal to the direction of range expansions-either expansions from Bering or posterior re-expansions to recolonize northern regions after ice sheets melting-and we suggest that they result from allele surfing processes. Conversely, our results on the entire continent show a northwest-southeast gradient obtained with any scenario, which we interpreted as a consequence of isolation by distance along the long length of the continent. These findings suggest that distinct genetic gradients can be detected at different regions of the Americas and that subcontinent regions present gradients more sensible to evolutionary and environmental factors (such as LDD and the LGM) than the whole continent.

Effect of hybridization with genome exclusion on extinction risk.

Human-induced habitat changes may lead to the breakdown of reproductive barriers between distantly related species. This phenomenon may result in fertile first-generation hybrids (F1 ) that exclude the genome of one parental species during gametogenesis, thus disabling introgression. The species extinction risk associated with hybridization with genome exclusion is largely underappreciated because the phenomenon produces only F1 hybrid phenotype, leading to the misconception that hybrids are sterile and potentially of minor conservation concern. We used a simulation model that integrates the main genetic, demographic, and ecological processes to examine the dynamics of hybridization with genome exclusion. We showed that this mode of hybridization may lead to extremely rapid extinction when the process of genome exclusion is unbalanced between the interbreeding species and when the hybridization rate is not negligible. The coexistence of parental species was possible in some cases of asymmetrical genome exclusion, but show this equilibrium was highly vulnerable to environmental variation. Expanding the exclusive habitat of the species at risk allowed its persistence. Our results highlight the extent of possible extinction risk due to hybridization with genome exclusion and suggest habitat management as a promising conservation strategy. In anticipation of serious threats to biodiversity due to hybridization with genome exclusion, we recommend a detailed assessment of the reproductive status of hybrids in conservation programs. We suggest such assessments include the inspection of genetic content in hybrid gametes.

Cryptic Biological Invasions: a General Model of Hybridization.

The dispersal of non-native genes due to hybridization is a form of cryptic invasion with growing concern in evolution and conservation. This includes the spread of transgenic genes and antibiotic resistance. To investigate how genes and phenotypes are transmitted, we developed a general model that, for the first time, considers concurrently: multiple loci, quantitative and qualitative gene expression, assortative mating, dominance/recessivity inheritance and density-dependent demographic effects. Selection acting on alleles or genotypes can also be incorporated. Our results reveal that the conclusions about how hybridization threatens a species can be biased if they are based on single-gene models, while considering two or more genes can correct this bias. We also show that demography can amplify or balance the genetic effects, evidencing the need of jointly incorporating both processes. By implementing our model in a real case, we show that mallard ducks introduced in New Zealand benefit from hybridization to replace native grey-ducks. Total displacement can take a few generations and occurs by interspecific competition and by competition between hybrids and natives, demonstrating how hybridization may facilitate biological invasions. We argue that our general model represents a powerful tool for the study of a wide range of biological and societal questions.

Investigating population continuity with ancient DNA under a spatially explicit simulation framework.

BACKGROUND: Recent advances in sequencing technologies have allowed for the retrieval of ancient DNA data (aDNA) from skeletal remains, providing direct genetic snapshots from diverse periods of human prehistory. Comparing samples taken in the same region but at different times, hereafter called "serial samples", may indicate whether there is continuity in the peopling history of that area or whether an immigration of a genetically different population has occurred between the two sampling times. However, the exploration of genetic relationships between serial samples generally ignores their geographical locations and the spatiotemporal dynamics of populations. Here, we present a new coalescent-based, spatially explicit modelling approach to investigate population continuity using aDNA, which includes two fundamental elements neglected in previous methods: population structure and migration. The approach also considers the extensive temporal and geographical variance that is commonly found in aDNA population samples. RESULTS: We first showed that our spatially explicit approach is more conservative than the previous (panmictic) approach and should be preferred to test for population continuity, especially when small and isolated populations are considered. We then applied our method to two mitochondrial datasets from Germany and France, both including modern and ancient lineages dating from the early Neolithic. The results clearly reject population continuity for the maternal line over the last 7500 years for the German dataset but not for the French dataset, suggesting regional heterogeneity in post-Neolithic migratory processes. CONCLUSIONS: Here, we demonstrate the benefits of using a spatially explicit method when investigating population continuity with aDNA. It constitutes an improvement over panmictic methods by considering the spatiotemporal dynamics of genetic lineages and the precise location of ancient samples. The method can be used to investigate population continuity between any pair of serial samples (ancient-ancient or ancient-modern) and to investigate more complex evolutionary scenarios. Although we based our study on mitochondrial DNA sequences, diploid molecular markers of different types (DNA, SNP, STR) can also be simulated with our approach. It thus constitutes a promising tool for the analysis of the numerous aDNA datasets being produced, including genome wide data, in humans but also in many other species.

Ancestry and demography and descendants of Iron Age nomads of the Eurasian Steppe.

During the 1st millennium before the Common Era (BCE), nomadic tribes associated with the Iron Age Scythian culture spread over the Eurasian Steppe, covering a territory of more than 3,500 km in breadth. To understand the demographic processes behind the spread of the Scythian culture, we analysed genomic data from eight individuals and a mitochondrial dataset of 96 individuals originating in eastern and western parts of the Eurasian Steppe. Genomic inference reveals that Scythians in the east and the west of the steppe zone can best be described as a mixture of Yamnaya-related ancestry and an East Asian component. Demographic modelling suggests independent origins for eastern and western groups with ongoing gene-flow between them, plausibly explaining the striking uniformity of their material culture. We also find evidence that significant gene-flow from east to west Eurasia must have occurred early during the Iron Age.

Early Neolithic genomes from the eastern Fertile Crescent.

We sequenced Early Neolithic genomes from the Zagros region of Iran (eastern Fertile Crescent), where some of the earliest evidence for farming is found, and identify a previously uncharacterized population that is neither ancestral to the first European farmers nor has contributed substantially to the ancestry of modern Europeans. These people are estimated to have separated from Early Neolithic farmers in Anatolia some 46,000 to 77,000 years ago and show affinities to modern-day Pakistani and Afghan populations, but particularly to Iranian Zoroastrians. We conclude that multiple, genetically differentiated hunter-gatherer populations adopted farming in southwestern Asia, that components of pre-Neolithic population structure were preserved as farming spread into neighboring regions, and that the Zagros region was the cradle of eastward expansion.