Biases in demographic modelling affect our understanding of recent divergence

Abstract

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (Ne) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Schopthalmus maximus) by testing 16 Isolation with Migration (IM) and 16 Secondary Contact (SC) scenarios, modelling changes in Ne as well as the effects of linked selection and barrier loci. Failure to account for changes in Ne resulted in selecting SC models with long periods of isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for Ne changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbots data islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by high-lighting the potential effects of unmodeleld changes in Ne on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modelling.

Genomic data reveal the biogeographical and demographic history of Ammospiza sparrows in northeast tidal marshes

“…Several factors (e.g. ancestral size changes and undetected bottlenecks, Momigliano et al, 2021;unac-counted gene flow from 'ghost' populations, Lynch & Sethuraman, 2019; low coverage or sequencing error, Gravel et al, 2011) could be contributing to discordance between our parameter estimates. We also note here that methodological approaches that fail to account for gene flow, such as SNAPP, can result in topological uncertainty (Thom et al, 2018), which could also be contributing to the observed discordance between our parameter estimates.…”

Deep connections: divergence histories with gene flow in mesophotic Agaricia corals

“…Although results from dadi gave the greatest support to the two-epoch models of divergence in isolation followed by secondary contacts (Tables 1 and 2), if population expansions or bottlenecks have occurred, divergence time can be over-estimated and results can favour the secondary contact model (Momigliano, Florin, & Merila, 2021). As RAD datasets with de novo assemblies are prone to high error rates in low frequency SNP calling (Andrews, Good, Miller, Luikart, & Hohenlohe, 2016), we did not include singletons and doubletons in the analyses and thus were unable to accurately model population size changes and differentiate between the different migration scenarios (i.e ., ancient migration vs. continuous migration vs. secondary contact).…”

Estimating Daily Dew Point Temperature Based on Local and Cross-Station Meteorological Data Using CatBoost Algorithm

“…From another point of view, the correct choice of parameters is very important. Moreover, the experience of the scholar determines the amount of work involved in the research [48][49][50]. In Fig.…”