Maximum-likelihood inference of population size contractions from microsatellite data.

Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method (MSVAR) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package.

Genetic structure of populations of whale sharks among ocean basins and evidence for their historic rise and recent decline.

This study presents genetic evidence that whale sharks, Rhincodon typus, are comprised of at least two populations that rarely mix and is the first to document a population expansion. Relatively high genetic structure is found when comparing sharks from the Gulf of Mexico with sharks from the Indo-Pacific. If mixing occurs between the Indian and Atlantic Oceans, it is not sufficient to counter genetic drift. This suggests whale sharks are not all part of a single global metapopulation. The significant population expansion we found was indicated by both microsatellite and mitochondrial DNA. The expansion may have happened during the Holocene, when tropical species could expand their range due to sea-level rise, eliminating dispersal barriers and increasing plankton productivity. However, the historic trend of population increase may have reversed recently. Declines in genetic diversity are found for 6 consecutive years at Ningaloo Reef in Australia. The declines in genetic diversity being seen now in Australia may be due to commercial-scale harvesting of whale sharks and collision with boats in past decades in other countries in the Indo-Pacific. The study findings have implications for models of population connectivity for whale sharks and advocate for continued focus on effective protection of the world's largest fish at multiple spatial scales.

ABLE: blockwise site frequency spectra for inferring complex population histories and recombination.

We introduce ABLE (Approximate Blockwise Likelihood Estimation), a novel simulation-based composite likelihood method that uses the blockwise site frequency spectrum to jointly infer past demography and recombination. ABLE is explicitly designed for a wide variety of data from unphased diploid genomes to genome-wide multi-locus data (for example, RADSeq) and can also accommodate arbitrarily large samples. We use simulations to demonstrate the accuracy of this method to infer complex histories of divergence and gene flow and reanalyze whole genome data from two species of orangutan. ABLE is available for download at https://github.com/champost/ABLE .

Do spatially-implicit estimates of neutral migration comply with seed dispersal data in tropical forests?

Neutral community models have shown that limited migration can have a pervasive influence on the taxonomic composition of local communities even when all individuals are assumed of equivalent ecological fitness. Notably, the spatially implicit neutral theory yields a single parameter I for the immigration-drift equilibrium in a local community. In the case of plants, seed dispersal is considered as a defining moment of the immigration process and has attracted empirical and theoretical work. In this paper, we consider a version of the immigration parameter I depending on dispersal limitation from the neighbourhood of a community. Seed dispersal distance is alternatively modelled using a distribution that decreases quickly in the tails (thin-tailed Gaussian kernel) and another that enhances the chance of dispersal events over very long distances (heavily fat-tailed Cauchy kernel). Our analysis highlights two contrasting situations, where I is either mainly sensitive to community size (related to ecological drift) under the heavily fat-tailed kernel or mainly sensitive to dispersal distance under the thin-tailed kernel. We review dispersal distances of rainforest trees from field studies and assess the consistency between published estimates of I based on spatially-implicit models and the predictions of the kernel-based model in tropical forest plots. Most estimates of I were derived from large plots (10-50 ha) and were too large to be accounted for by a Cauchy kernel. Conversely, a fraction of the estimates based on multiple smaller plots (1 ha) appeared too small to be consistent with reported ranges of dispersal distances in tropical forests. Very large estimates may reflect within-plot habitat heterogeneity or estimation problems, while the smallest estimates likely imply other factors inhibiting migration beyond dispersal limitation. Our study underscores the need for interpreting I as an integrative index of migration limitation which, besides the limited seed dispersal, possibly includes habitat filtering or fragmentation.