Incorporating non-equilibrium dynamics into demographic history inferences of a migratory marine species

ELC was supported while writing this paper by a EU Horizon 2020 Marie Slodowska Curie Fellowship, project BEHAVIOUR-CONNECT, by a Newton Fellowship from the Royal Society of London and Bayesian statistical training was supported by National Science Foundation (award DEB- 1145200). Laboratory analyses conducted by ELC were funded by a small grant from the British Ecological Society 5076 / 6118 and Bayesian analysis was supported by training from the National Science Foundation under Grant No. DEB-1145200. OEG was supported by the Marine Alliance for Science and Technology for Scotland (MASTS) funded by the Scottish Founding Council (grant reference HR09011). Genetic data from the South African right whale samples were generated by MB and PJP with the support of UC Berkeley, University of Stockholm and University of Groningen. Computational Biology analyses were supported by the University of St Andrews Bioinformatics Unit which is funded by a Wellcome Trust ISSF award.Understanding how dispersal and gene flow link geographically separated populations over evolutionary history is challenging, particularly in migratory marine species. In southern right whales (SRWs, Eubalaena australis), patterns of genetic diversity are likely influenced by the glacial climate cycle and recent history of whaling. Here we use a dataset of mitochondrial DNA (mtDNA) sequences (n=1,327) and nuclear markers (17 microsatellite loci, n=222) from major wintering grounds to investigate circumpolar population structure, historical demography, and effective population size. Analyses of nuclear genetic variation identify two population clusters that correspond to the South Atlantic and Indo-Pacific ocean basins that have similar effective breeder estimates. In contrast, all wintering grounds show significant differentiation for mtDNA, but no sex-biased dispersal was detected using the microsatellite genotypes. An approximate Bayesian computation (ABC) approach with microsatellite markers compared scenarios with gene flow through time, or isolation and secondary contact between ocean basins, while modeling declines in abundance linked to whaling. Secondary-contact scenarios yield the highest posterior probabilities, implying that populations in different ocean basins were largely isolated and came into secondary contact within the last 25,000 years, but the role of whaling in changes in genetic diversity and gene flow over recent generations could not be resolved. We hypothesis that these findings are driven by factors that promote isolation, such as female philopatry, and factors that could promote dispersal, such oceanographic changes. These findings highlight the application of ABC approaches to infer connectivity in mobile species with complex population histories and currently low levels of differentiation.PostprintPeer reviewe

The trials and triumphs of genius

A digital copy of a pamphlet. 1 PDF. 30 pages

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the 'species trait' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics