Inference of Gene Flow in the Process of Speciation: An Efficient Maximum-Likelihood Method for the Isolation-with-Initial-Migration Model

Abstract

The isolation-with-migration (IM) model is commonly used to make inferences about gene flow during speciation, using polymorphism data. However, it has been reported that the parameter estimates obtained by fitting the IM model are very sensitive to the model’s assumptions—including the assumption of constant gene flow until the present. This article is concerned with the isolation-with-initial-migration (IIM) model, which drops precisely this assumption. In the IIM model, one ancestral population divides into two descendant subpopulations, between which there is an initial period of gene flow and a subsequent period of isolation. We derive a very fast method of fitting an extended version of the IIM model, which also allows for asymmetric gene flow and unequal population sizes. This is a maximum-likelihood method, applicable to data on the number of segregating sites between pairs of DNA sequences from a large number of independent loci. In addition to obtaining parameter estimates, our method can also be used, by means of likelihood-ratio tests, to distinguish between alternative models representing the following divergence scenarios: (a) divergence with potentially asymmetric gene flow until the present, (b) divergence with potentially asymmetric gene flow until some point in the past and in isolation since then, and (c) divergence in complete isolation. We illustrate the procedure on pairs of Drosophila sequences from ∼30,000 loci. The computing time needed to fit the most complex version of the model to this data set is only a couple of minutes. The R code to fit the IIM model can be found in the supplementary files of this article.

Whole genome data reveal the complex history of a diverse ecological community

“…[32]) there is no requirement for an outgroup. Analogous composite likelihood strategies based on subsamples have been used in phylogenetic network analysis [M12] and to infer divergence and continuous gene flow [M13, M14]. We maximized InCL for all parameters ( N E , N e , N W , T 1 , T 2 , T adm , f ) included in the best supported model using a Nelder-Mead simplex algorithm implemented in the Mathematica function NMaximize.…”

Demographic and genetic approaches to study dispersal in wild animal populations: A methodological review

“…Recent model development has relaxed some of the previous assumptions using different variants of the IM model. For instance, it is now possible to include both asymmetric migration and variable population size (Costa & Wilkinson‐Herbots, ). It is also now possible to infer complex histories using joint information from the blockwise site frequency spectrum and linkage disequilibrium (Beeravolu Reddy, Hickerson, Frantz, & Lohse, ) or to perform exact calculation of the joint allele frequency spectrum under an IM model using Markov chain representation of the coalescence (Kern & Hey, ).…”

A legacy of geo‐climatic complexity and genetic divergence along the lower Colorado River: Insights from the geological record and 33 desert‐adapted animals

“…The genetic data in this study cannot adjudicate between these scenarios, but future studies using whole-genome sequencing have promise. Methods such as δaδi (Gutenkunst, Hernandez, Williamson, & Bustamante, 2009) and iiM (Costa & Wilkinson-Herbots, 2017) can be used to assess the fit of speciation models, such as speciation in isolation, with initial migration or with continuous migration. Methods to differentiate between incomplete lineage sorting and post-vicariance admixture can be achieved by statistically comparing consistency of tree topologies (e.g.…”