Microsatellite null alleles and estimation of population differentiation

Microsatellite null alleles are commonly encountered in population genetics studies, yet little is known about their impact on the estimation of population differentiation. Computer simulations based on the coalescent were used to investigate the evolutionary dynamics of null alleles, their impact on F(ST) and genetic distances, and the efficiency of estimators of null allele frequency. Further, we explored how the existing method for correcting genotype data for null alleles performed in estimating F(ST) and genetic distances, and we compared this method with a new method proposed here (for F(ST) only). Null alleles were likely to be encountered in populations with a large effective size, with an unusually high mutation rate in the flanking regions, and that have diverged from the population from which the cloned allele state was drawn and the primers designed. When populations were significantly differentiated, F(ST) and genetic distances were overestimated in the presence of null alleles. Frequency of null alleles was estimated precisely with the algorithm presented in Dempster et al. (1977). The conventional method for correcting genotype data for null alleles did not provide an accurate estimate of F(ST) and genetic distances. However, the use of the genetic distance of Cavalli-Sforza and Edwards (1967) corrected by the conventional method gave better estimates than those obtained without correction. F(ST) estimation from corrected genotype frequencies performed well when restricted to visible allele sizes. Both the proposed method and the traditional correction method have been implemented in a program that is available free of charge at http://www.montpellier.inra.fr/URLB/. We used 2 published microsatellite data sets based on original and redesigned pairs of primers to empirically confirm our simulation results.     

msABC: a modification of Hudson's ms to facilitate multi-locus ABC analysis

With the availability of whole-genome sequence data biologists are able to test hypotheses regarding the demography of populations. Furthermore, the advancement of the Approximate Bayesian Computation (ABC) methodology allows the demographic inference to be performed in a simple framework using summary statistics. We present here msABC, a coalescent-based software that facilitates the simulation of multi-locus data, suitable for an ABC analysis. msABC is based on Hudson's ms algorithm, which is used extensively for simulating neutral demographic histories of populations. The flexibility of the original algorithm has been extended so that sample size may vary among loci, missing data can be incorporated in simulations and calculations, and a multitude of summary statistics for single or multiple populations is generated. The source code of msABC is available at http://bio.lmu.de/~pavlidis/msabc or upon request from the authors.     

Positive selection versus demography: evolutionary inferences based on an unusual haplotype structure in Drosophila simulans

Coalescent simulations were used to investigate the possible role of population subdivision and history in shaping nucleotide variation in a recombining 88-kb genomic fragment of Drosophila simulans displaying an unusual large-scale haplotype structure. The multilocus analysis, based on summary statistics using specific demographic null models under recombination, indicates that the observed levels of linkage disequilibrium differed significantly from the values expected under different bottleneck and population admixture scenarios. These results indicate that demography alone may not account for the observed pattern of variation and support the previous claim that the data are better described by a model in which an adaptive mutation has not yet gone to fixation.     

Unexpected persistence on habitat islands: genetic signatures reveal dispersal of a eucalypt-dependent marsupial through a hostile pine matrix

Several factors contribute to the extinction of populations in fragmented habitat but key ones include habitat loss and disruptions to connectivity. Aspects of the ecology of greater gliders (Petauroides volans), along with observations of their response to native forest clearance at a site in southeastern Australia, lead to the prediction in the 1960s that the species would not persist in the replacement exotic pine plantation. However, 35 years later, the species was observed in many remnant native vegetation patches retained within the plantation boundary, albeit at a lower occupancy rate than at matched continuous forest control sites. To determine the role of patch connectivity in persistence of P. volans in remnants, we employed 12 microsatellite markers to genotype individuals from 11 remnants, three contemporary nearby continuous native eucalypt forest sites and a sample collected during native vegetation clearance at the site in the 1960s. Patch samples retained substantially more genetic diversity than expected under an isolation model, suggesting that patches have experienced some immigration. Five putative patch immigrants--two from sampled sites 1- and 7-km distant, and three from unresolved or unsampled localities--were identified via genetic parentage and population assignment analyses. Patch populations displayed varying levels of admixture in Bayesian genetic structure analyses, with the oldest and most geographically isolated ones showing the least admixture, suggesting they have experienced relatively little immigration. Evidence of at least some immigration into patches may explain why P. volans has persisted contrary to expectation in heavily fragmented habitat.     

MARLIN, software to create, run, and analyse spatially realistic simulations

MARLIN is a software to create, run, analyse, and visualize spatially explicit population genetic simulations. It provides an intuitive user interface with which the geographical layout of a metapopulation can be drawn by hand or loaded from a map. Furthermore, the interface allows easy selection of the many different simulation settings. MARLIN then uses the program QuantiNemo to run the simulation in the background. When simulations are finished, MARLIN directly analyses and plots the results, thereby greatly simplifying the simulation workflow. This combination of simulation and analysis makes MARLIN ideal for teaching and for scientists who are interested in doing simulations without having to learn command-line operations. MARLIN is available for computers running Mac OS X and can be downloaded from: http://www.patrickmeirmans.com/software.     

Blind population genetics survey of tropical rainforest trees

Rainforest tree species can be difficult to identify outside of their period of reproduction. Vascular tissues from Carapa spp. individuals were collected during a short field trip in French Guiana and analysed in the laboratory with nuclear and chloroplast markers. Using a Bayesian approach, > 90% of the samples could be assigned to one of two distinct clusters corresponding to previously described species, making it possible to estimate the genetic structure of each species and to identify cases of introgression. We argue that this blind procedure represents a first-choice rather than a fallback option whenever related taxa are investigated.     

Demography or selection on linked cultural traits or genes? Investigating the driver of low mtDNA diversity in the sperm whale using complementary mitochondrial and nuclear genome analyses

Mitochondrial DNA has been heavily utilized in phylogeography studies for several decades. However, underlying patterns of demography and phylogeography may be misrepresented due to coalescence stochasticity, selection, variation in mutation rates and cultural hitchhiking (linkage of genetic variation to culturally‐transmitted traits affecting fitness). Cultural hitchhiking has been suggested as an explanation for low genetic diversity in species with strong social structures, counteracting even high mobility, abundance and limited barriers to dispersal. One such species is the sperm whale, which shows very limited phylogeographic structure and low mtDNA diversity despite a worldwide distribution and large population. Here, we use analyses of 175 globally distributed mitogenomes and three nuclear genomes to evaluate hypotheses of a population bottleneck/expansion vs. a selective sweep due to cultural hitchhiking or selection on mtDNA as the mechanism contributing to low worldwide mitochondrial diversity in sperm whales. In contrast to mtDNA control region (CR) data, mitogenome haplotypes are largely ocean‐specific, with only one of 80 shared between the Atlantic and Pacific. Demographic analyses of nuclear genomes suggest low mtDNA diversity is consistent with a global reduction in population size that ended approximately 125,000 years ago, correlated with the Eemian interglacial. Phylogeographic analysis suggests that extant sperm whales descend from maternal lineages endemic to the Pacific during the period of reduced abundance and have subsequently colonized the Atlantic several times. Results highlight the apparent impact of past climate change, and suggest selection and hitchhiking are not the sole processes responsible for low mtDNA diversity in this highly social species.