Learning about modes of speciation by computational approaches

Biases in Demographic Modeling Affect Our Understanding of Recent Divergence

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 (N_e) 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 (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in N_e as well as the effects of linked selection and barrier loci. Failure to account for changes in N_e resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for N_e 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 turbot 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 highlighting the potential effects of unmodeled changes in N_e 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 modeling.

The effect of undetected recombination on genealogy sampling and inference under an isolation-with-migration model

Many methods for fitting demographic models to data sets of aligned sequences rely upon an assumption that the data have a branching coalescent history without recombination within regions or loci. To mitigate the effects of the failure of this assumption, a common approach is to filter data and sample regions that pass the four-gamete criterion for recombination, an approach that allows data to run, but that is expected to detect only a minority of recombination events. A series of empirical tests of this approach were conducted using computer simulations with and without recombination for a variety of isolation-with-migration (IM) model for two and three populations. Only the IMa3 program was used, but the general results should apply to related genealogy-sampling-based methods for IM models or subsets of IM models. It was found that the details of sampling intervals that pass a four-gamete filter have a moderate effect, and that schemes that use the longest intervals, or that use overlapping intervals, gave poorer results. A simple approach of using a random nonoverlapping interval returned the smallest difference between results with and without recombination, with the mean difference between parameter estimates usually less than 20% of the true value (usually much less). However, the posterior probability distributions for migration rates were flatter with recombination, suggesting that filtering based on the four-gamete criterion, while necessary for methods like these, leads to reduced resolution on migration. A distinct, alternative approach, of using a finite sites mutation model and not filtering the data, performed quite poorly.

Can secondary contact following range expansion be distinguished from barriers to gene flow?

Secondary contact is the reestablishment of gene flow between sister populations that have diverged. For instance, at the end of the Quaternary glaciations in Europe, secondary contact occurred during the northward expansion of the populations which had found refugia in the southern peninsulas. With the advent of multi-locus markers, secondary contact can be investigated using various molecular signatures including gradients of allele frequency, admixture clines, and local increase of genetic differentiation. We use coalescent simulations to investigate if molecular data provide enough information to distinguish between secondary contact following range expansion and an alternative evolutionary scenario consisting of a barrier to gene flow in an isolation-by-distance model. We find that an excess of linkage disequilibrium and of genetic diversity at the suture zone is a unique signature of secondary contact. We also find that the directionality index ψ, which was proposed to study range expansion, is informative to distinguish between the two hypotheses. However, although evidence for secondary contact is usually conveyed by statistics related to admixture coefficients, we find that they can be confounded by isolation-by-distance. We recommend to account for the spatial repartition of individuals when investigating secondary contact in order to better reflect the complex spatio-temporal evolution of populations and species.

The search for loci under selection: trends, biases and progress

Detecting genetic variants under selection using F_ST outlier analysis (OA) and environmental association analyses (EAAs) are popular approaches that provide insight into the genetic basis of local adaptation. Despite the frequent use of OA and EAA approaches and their increasing attractiveness for detecting signatures of selection, their application to field-based empirical data have not been synthesized. Here, we review 66 empirical studies that use Single Nucleotide Polymorphisms (SNPs) in OA and EAA. We report trends and biases across biological systems, sequencing methods, approaches, parameters, environmental variables and their influence on detecting signatures of selection. We found striking variability in both the use and reporting of environmental data and statistical parameters. For example, linkage disequilibrium among SNPs and numbers of unique SNP associations identified with EAA were rarely reported. The proportion of putatively adaptive SNPs detected varied widely among studies, and decreased with the number of SNPs analysed. We found that genomic sampling effort had a greater impact than biological sampling effort on the proportion of identified SNPs under selection. OA identified a higher proportion of outliers when more individuals were sampled, but this was not the case for EAA. To facilitate repeatability, interpretation and synthesis of studies detecting selection, we recommend that future studies consistently report geographical coordinates, environmental data, model parameters, linkage disequilibrium, and measures of genetic structure. Identifying standards for how OA and EAA studies are designed and reported will aid future transparency and comparability of SNP-based selection studies and help to progress landscape and evolutionary genomics.

A neutral view of the evolving genomic architecture of speciation

Analyses of genomewide polymorphism data have begun to shed light on speciation and adaptation. Genome scans to identify regions of the genome that are unusually different between populations or species, possibly due to divergent natural or sexual selection, are widespread in speciation genomics. Theoretical and empirical work suggests that such outlier regions may grow faster than linearly during speciation with gene flow due to a rapid transition between low and high reproductive isolation. We investigate whether this pattern could be attributed to neutral processes by simulating genomes under neutral evolution with varying amounts and timing of gene flow. Under both neutral evolution and divergent selection, simulations with little or no gene flow, or with a long allopatric period after its cessation, resulted in faster than linear growth of the proportion of the genome lying in outlier regions. Without selection, higher recent gene flow erased differentiation; with divergent selection, these same scenarios produced nonlinear growth to a plateau. Our results suggest that, given a history of gene flow, the growth of the divergent genome is informative about selection during divergence, but that in many scenarios, this pattern does not easily distinguish neutral and non-neutral processes during speciation with gene flow.

Limited genomic consequences of hybridization between two African clawed frogs, Xenopus gilli and X. laevis (Anura: Pipidae)

The Cape platanna, Xenopus gilli, an endangered frog, hybridizes with the African clawed frog, X. laevis, in South Africa. Estimates of the extent of gene flow between these species range from pervasive to rare. Efforts have been made in the last 30 years to minimize hybridization between these two species in the west population of X. gilli, but not the east populations. To further explore the impact of hybridization and the efforts to minimize it, we examined molecular variation in one mitochondrial and 13 nuclear genes in genetic samples collected recently (2013) and also over two decades ago (1994). Despite the presence of F₁ hybrids, none of the genomic regions we surveyed had evidence of gene flow between these species, indicating a lack of extensive introgression. Additionally we found no significant effect of sampling time on genetic diversity of populations of each species. Thus, we speculate that F₁ hybrids have low fitness and are not backcrossing with the parental species to an appreciable degree. Within X. gilli, evidence for gene flow was recovered between eastern and western populations, a finding that has implications for conservation management of this species and its threatened habitat.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

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

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.

Wolbachia in the Drosophila yakuba Complex: Pervasive Frequency Variation and Weak Cytoplasmic Incompatibility, but No Apparent Effect on Reproductive Isolation

Three hybridizing species-the clade [(Drosophila yakuba, D. santomea), D. teissieri]-comprise the yakuba complex in the D. melanogaster subgroup. Their ranges overlap on Bioko and São Tomé, islands off west Africa. All three species are infected with Wolbachia-maternally inherited, endosymbiotic bacteria, best known for manipulating host reproduction to favor infected females. Previous analyses reported no cytoplasmic incompatibility (CI) in these species. However, we discovered that Wolbachia from each species cause intraspecific and interspecific CI. In D teissieri, analyses of F1 and backcross genotypes show that both host genotype and Wolbachia variation modulate CI intensity. Wolbachia-infected females seem largely protected from intraspecific and interspecific CI, irrespective of Wolbachia and host genotypes. Wolbachia do not affect host mating behavior or female fecundity, within or between species. The latter suggests little apparent effect of Wolbachia on premating or gametic reproductive isolation (RI) between host species. In nature, Wolbachia frequencies varied spatially for D. yakuba in 2009, with 76% (N = 155) infected on São Tomé, and only 3% (N = 36) infected on Bioko; frequencies also varied temporally in D. yakuba and D. santomea on São Tomé between 2009 and 2015. These temporal frequency fluctuations could generate asymmetries in interspecific mating success, and contribute to postzygotic RI. However, the fluctuations in Wolbachia frequencies that we observe also suggest that asymmetries are unlikely to persist. Finally, we address theoretical questions that our empirical findings raise about Wolbachia persistence when conditions fluctuate, and about the stable coexistence of Wolbachia and host variants that modulate Wolbachia effects.

Identifying Genetic Signatures of Natural Selection Using Pooled Population Sequencing in Picea abies

The joint inference of selection and past demography remain a costly and demanding task. We used next generation sequencing of two pools of 48 Norway spruce mother trees, one corresponding to the Fennoscandian domain, and the other to the Alpine domain, to assess nucleotide polymorphism at 88 nuclear genes. These genes are candidate genes for phenological traits, and most belong to the photoperiod pathway. Estimates of population genetic summary statistics from the pooled data are similar to previous estimates, suggesting that pooled sequencing is reliable. The nonsynonymous SNPs tended to have both lower frequency differences and lower F_ST values between the two domains than silent ones. These results suggest the presence of purifying selection. The divergence between the two domains based on synonymous changes was around 5 million yr, a time similar to a recent phylogenetic estimate of 6 million yr, but much larger than earlier estimates based on isozymes. Two approaches, one of them novel and that considers both F_ST and difference in allele frequencies between the two domains, were used to identify SNPs potentially under diversifying selection. SNPs from around 20 genes were detected, including genes previously identified as main target for selection, such as PaPRR3 and PaGI.

The genomics of ecological vicariance in threespine stickleback fish

Populations occurring in similar habitats and displaying similar phenotypes are increasingly used to explore parallel evolution at the molecular level. This generally ignores the possibility that parallel evolution can be mimicked by the fragmentation of an ancestral population followed by genetic exchange with ecologically different populations. Here we demonstrate such an ecological vicariance scenario in multiple stream populations of threespine stickleback fish divergent from a single adjacent lake population. On the basis of demographic and population genomic analyses, we infer the initial spread of a stream-adapted ancestor followed by the emergence of a lake-adapted population, that selective sweeps have occurred mainly in the lake population, that adaptive lake–stream divergence is maintained in the face of gene flow from the lake into the streams, and that this divergence involves major inversion polymorphisms also important to marine-freshwater stickleback divergence. Overall, our study highlights the need for a robust understanding of the demographic and selective history in evolutionary investigations.

Threespine stickleback fish are adapted to lake and stream habitats in Central Europe. Here, the authors show colonization of a lake basin by a stream-adapted ancestor, followed by the emergence of a lake-adapted population in the face of gene flow across lake–stream boundaries.

The Genomic Signature of Population Reconnection Following Isolation: From Theory to HIV

Ease of worldwide travel provides increased opportunities for organisms not only to colonize new environments but also to encounter related but diverged populations. Such events of reconnection and secondary contact of previously isolated populations are widely observed at different time scales. For example, during the quaternary glaciation, sea water level fluctuations caused temporal isolation of populations, often to be followed by secondary contact. At shorter time scales, population isolation and reconnection of viruses are commonly observed, and such events are often associated with epidemics and pandemics. Here, using coalescent theory and simulations, we describe the temporal impact of population reconnection after isolation on nucleotide differences and the site frequency spectrum, as well as common summary statistics of DNA variation. We identify robust genomic signatures of population reconnection after isolation. We utilize our development to infer the recent evolutionary history of human immunodeficiency virus 1 (HIV-1) in Asia and South America, successfully retrieving the successive HIV subtype colonization events in these regions. Our analysis reveals that divergent HIV-1 subtype populations are currently admixing in these regions, suggesting that HIV-1 may be undergoing a process of homogenization, contrary to popular belief.

On the occurrence of false positives in tests of migration under an isolation-with-migration model

The population genetic study of divergence is often carried out using a Bayesian genealogy sampler, like those implemented in ima2 and related programs, and these analyses frequently include a likelihood ratio test of the null hypothesis of no migration between populations. Cruickshank and Hahn (2014, Molecular Ecology, 23, 3133-3157) recently reported a high rate of false-positive test results with ima2 for data simulated with small numbers of loci under models with no migration and recent splitting times. We confirm these findings and discover that they are caused by a failure of the assumptions underlying likelihood ratio tests that arises when using marginal likelihoods for a subset of model parameters. We also show that for small data sets, with little divergence between samples from two populations, an excellent fit can often be found by a model with a low migration rate and recent splitting time and a model with a high migration rate and a deep splitting time.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Comparative species divergence across eight triplets of spiny lizards (Sceloporus) using genomic sequence data

Species divergence is typically thought to occur in the absence of gene flow, but many empirical studies are discovering that gene flow may be more pervasive during species formation. Although many examples of divergence with gene flow have been identified, few clades have been investigated in a comparative manner, and fewer have been studied using genome-wide sequence data. We contrast species divergence genetic histories across eight triplets of North American Sceloporus lizards using a maximum likelihood implementation of the isolation–migration (IM) model. Gene flow at the time of species divergence is modeled indirectly as variation in species divergence time across the genome or explicitly using a migration rate parameter. Likelihood ratio tests (LRTs) are used to test the null model of no gene flow at speciation against these two alternative gene flow models. We also use the Akaike information criterion to rank the models. Hundreds of loci are needed for the LRTs to have statistical power, and we use genome sequencing of reduced representation libraries to obtain DNA sequence alignments at many loci (between 340 and 3,478; mean = 1,678) for each triplet. We find that current species distributions are a poor predictor of whether a species pair diverged with gene flow. Interrogating the genome using the triplet method expedites the comparative study of species divergence history and the estimation of genetic parameters associated with speciation.

Speciation of two desert poplar species triggered by Pleistocene climatic oscillations

Despite the evidence that the Pleistocene climatic fluctuations have seriously affected the distribution of intraspecific diversity, less is known on its impact on interspecific divergence. In this study, we aimed to test the hypothesis that the divergence of two desert poplar species Populus euphratica Oliv. and P. pruinosa Schrenk. occurred during the Pleistocene. We sequenced 11 nuclear loci in 60 individuals from the two species to estimate the divergence time between them and to test whether gene flow occurred after species separation. Divergence time between the two species was estimated to be 0.66-1.37 million years ago (Ma), a time at which glaciation was at its maximum in China and deserts developed widely in central Asia. Isolation-with-Migration model also indicated that the two species had diverged in the presence of gene flow. We also detected evidence of selection at GO in P. euphratica and to a lesser extent at PhyB2. Together, these results underscore the importance of Pleistocene climate oscillations in triggering plant speciation as a result of habitats divergence.

Parallel evolution of local adaptation and reproductive isolation in the face of gene flow

Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat-associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky-shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

Origin and demographic history of the endemic Taiwan spruce (Picea morrisonicola)

Taiwan spruce (Picea morrisonicola) is a vulnerable conifer species endemic to the island of Taiwan. A warming climate and competition from subtropical tree species has limited the range of Taiwan spruce to the higher altitudes of the island. Using seeds sampled from an area in the central mountain range of Taiwan, 15 nuclear loci were sequenced in order to measure genetic variation and to assess the long-term genetic stability of the species. Genetic diversity is low and comparable to other spruce species with limited ranges such as Picea breweriana, Picea chihuahuana, and Picea schrenkiana. Importantly, analysis using approximate Bayesian computation (ABC) provides evidence for a drastic decline in the effective population size approximately 0.3–0.5 million years ago (mya). We used simulations to show that this is unlikely to be a false-positive result due to the limited sample used here. To investigate the phylogenetic origin of Taiwan spruce, additional sequencing was performed in the Chinese spruce Picea wilsonii and combined with previously published data for three other mainland China species, Picea purpurea, Picea likiangensis, and P. schrenkiana. Analysis of population structure revealed that P. morrisonicola clusters most closely with P. wilsonii, and coalescent analyses using the program MIMAR dated the split to 4–8 mya, coincidental to the formation of Taiwan. Considering the population decrease that occurred after the split, however, led to a much more recent origin.

A genome-wide survey of genetic variation in gorillas using reduced representation sequencing

All non-human great apes are endangered in the wild, and it is therefore important to gain an understanding of their demography and genetic diversity. Whole genome assembly projects have provided an invaluable foundation for understanding genetics in all four genera, but to date genetic studies of multiple individuals within great ape species have largely been confined to mitochondrial DNA and a small number of other loci. Here, we present a genome-wide survey of genetic variation in gorillas using a reduced representation sequencing approach, focusing on the two lowland subspecies. We identify 3,006,670 polymorphic sites in 14 individuals: 12 western lowland gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla beringei graueri). We find that the two species are genetically distinct, based on levels of heterozygosity and patterns of allele sharing. Focusing on the western lowland population, we observe evidence for population substructure, and a deficit of rare genetic variants suggesting a recent episode of population contraction. In western lowland gorillas, there is an elevation of variation towards telomeres and centromeres on the chromosomal scale. On a finer scale, we find substantial variation in genetic diversity, including a marked reduction close to the major histocompatibility locus, perhaps indicative of recent strong selection there. These findings suggest that despite their maintaining an overall level of genetic diversity equal to or greater than that of humans, population decline, perhaps associated with disease, has been a significant factor in recent and long-term pressures on wild gorilla populations.

Molecular hyperdiversity and evolution in very large populations

The genomic density of sequence polymorphisms critically affects the sensitivity of inferences about ongoing sequence evolution, function and demographic history. Most animal and plant genomes have relatively low densities of polymorphisms, but some species are hyperdiverse with neutral nucleotide heterozygosity exceeding 5%. Eukaryotes with extremely large populations, mimicking bacterial and viral populations, present novel opportunities for studying molecular evolution in sexually reproducing taxa with complex development. In particular, hyperdiverse species can help answer controversial questions about the evolution of genome complexity, the limits of natural selection, modes of adaptation and subtleties of the mutation process. However, such systems have some inherent complications and here we identify topics in need of theoretical developments. Close relatives of the model organisms Caenorhabditis elegans and Drosophila melanogaster provide known examples of hyperdiverse eukaryotes, encouraging functional dissection of resulting molecular evolutionary patterns. We recommend how best to exploit hyperdiverse populations for analysis, for example, in quantifying the impact of noncrossover recombination in genomes and for determining the identity and micro-evolutionary selective pressures on noncoding regulatory elements.

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π(sil)=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F(st)=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

Demographic histories of adaptively diverged riparian and non-riparian species of Ainsliaea (Asteraceae) inferred from coalescent analyses using multiple nuclear loci

BACKGROUND

Understanding demographic histories, such as divergence time, patterns of gene flow, and population size changes, in ecologically diverging lineages provide implications for the process and maintenance of population differentiation by ecological adaptation. This study addressed the demographic histories in two independently derived lineages of flood-resistant riparian plants and their non-riparian relatives [Ainsliaea linearis (riparian) and A. apiculata (non-riparian); A. oblonga (riparian) and A. macroclinidioides (non-riparian); Asteraceae] using an isolation-with-migration (IM) model based on variation at 10 nuclear DNA loci.

RESULTS

The highest posterior probabilities of the divergence time parameters were estimated to be ca. 25,000 years ago for A. linearis and A. apiculata and ca. 9000 years ago for A. oblonga and A. macroclinidioides, although the confidence intervals of the parameters had broad ranges. The likelihood ratio tests detected evidence of historical gene flow between both riparian/non-riparian species pairs. The riparian populations showed lower levels of genetic diversity and a significant reduction in effective population sizes compared to the non-riparian populations and their ancestral populations.

CONCLUSIONS

This study showed the recent origins of flood-resistant riparian plants, which are remarkable examples of plant ecological adaptation. The recent divergence and genetic signatures of historical gene flow among riparian/non-riparian species implied that they underwent morphological and ecological differentiation within short evolutionary timescales and have maintained their species boundaries in the face of gene flow. Comparative analyses of adaptive divergence in two sets of riparian/non-riparian lineages suggested that strong natural selection by flooding had frequently reduced the genetic diversity and size of riparian populations through genetic drift, possibly leading to fixation of adaptive traits in riparian populations. The two sets of riparian/non-riparian lineages showed contrasting patterns of gene flow and genetic differentiation, implying that each lineage showed different degrees of reproductive isolation and that they had experienced unique evolutionary and demographic histories in the process of adaptive divergence.

Introgression of wing pattern alleles and speciation via homoploid hybridization in Heliconius butterflies: a review of evidence from the genome

The diverse Müllerian mimetic wing patterns of neotropical Heliconius (Nymphalidae) have been proposed to be not only aposematic signals to potential predators, but also intra- and interspecific recognition signals that allow the butterflies to maintain their specific identities, and which perhaps drive the process of speciation, as well. Adaptive features under differential selection that also serve as cues for assortative mating have been referred to as 'magic traits', which can drive ecological speciation. Such traits are expected to exhibit allelic differentiation between closely related species with ongoing gene flow, whereas unlinked neutral traits are expected to be homogenized to a greater degree by introgression. However, recent evidence suggests that interspecific hybridization among Heliconius butterflies may have resulted in adaptive introgression of these very same traits across species boundaries, and in the evolution of new species by homoploid hybrid speciation. The theory and data supporting various aspects of the apparent paradox of 'magic trait' introgression are reviewed, with emphasis on population genomic comparisons of Heliconius melpomene and its close relatives.

Revising the human mutation rate: implications for understanding human evolution

It is now possible to make direct measurements of the mutation rate in modern humans using next-generation sequencing. These measurements reveal a value that is approximately half of that previously derived from fossil calibration, and this has implications for our understanding of demographic events in human evolution and other aspects of population genetics. Here, we discuss the implications of a lower-than-expected mutation rate in relation to the timescale of human evolution.

A European phylogeography of Rhinanthus minor compared to Rhinanthus angustifolius: unexpected splits and signs of hybridization

Rhinanthus minor and Rhinanthus angustifolius (Orobanchaceae) are annual hemiparasites, which occur sympatrically in Europe and are known to hybridize. We studied chloroplast and nuclear (amplified fragment length polymorphism [AFLP]) diversity in R. minor and compared genetic structuring in this species with R. angustifolius by analyzing the AFLP data for both species simultaneously. The AFLP data revealed that populations in Italy, Greece, and southeast Russia initially identified as R. minor were so distant from the other R. minor populations that they probably belong to another, yet unidentified taxon, and we refer to them as Rhinanthus sp. R. minor s.s. showed a clear geographic genetic structure in both the chloroplast DNA (cpDNA) and nuclear genome. The simultaneous analysis of both species shed new light on the previously published findings for R. angustifolius, because some populations now turned out to belong to R. minor. The admixture analysis revealed very few individuals of mixed R. minor–R.angustifolius ancestry in the natural populations in the west of Europe, while admixture levels were higher in the east. The combined haplotype network showed that haplotype H1 was shared among all species and is likely to be ancestral. H2 was more abundant in R. angustifolius and H3 in R. minor, and the latter probably arose from H1 in this species in the east of Europe. The occurrence of H3 in R. angustifolius may be explained by introgression from R. minor, but without interspecific admixture, these are likely to have been old hybridization events. Our study underlines the importance of including related species in phylogeographic studies.

Speciation history of three closely related oak gall wasps, Andricus mukaigawae, A. kashiwaphilus, and A. pseudoflos (Hymenoptera: Cynipidae) inferred from nuclear and mitochondrial DNA sequences

The Andricus mukaigawae complex of oak gall wasps is composed of cyclically parthenogenetic species: A. mukaigawae and Andricus kashiwaphilus, and a parthenogenetic species, Andricus pseudoflos. The component species differ in life history, host plant, karyotype, and asexual gall shape, although little difference is found in the external morphology of asexual adults. To understand the speciation history of this species complex, DNA sequences of one mitochondrial region and nine nuclear gene regions were investigated. The genetic relationship among the species suggested that a loss of sex occurred after host shift. Unexpectedly, two or three distinct groups in the parthenogenetic species, A. pseudoflos, were revealed by both mitochondrial and nuclear DNA data. Gene flow in nuclear genes from the species not infected by Wolbachia (A. kashiwaphilus) to the species infected by it (A. mukaigawae) was suggested by a method based on coalescent simulations. On the other hand, gene flow in mitochondrial genes was suggested to be in the opposite direction. These findings indicate possible involvement of Wolbachia infection in the speciation process of the A. mukaigawae complex.

Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae)

PREMISE OF THE STUDY

Plant populations that face new environments adapt and diverge simultaneously, and both processes leave footprints in their genetic diversity. Arabidopsis lyrata is an excellent species for studying these processes. Pairs of populations and subspecies of A. lyrata represent different stages of divergence. These populations are also known to be locally adapted and display various stages of emerging reproductive isolation.

METHODS

We used nucleotide diversity data from 19 loci to estimate divergence times and levels of diversity among nine A. lyrata populations. Traditional distance-based methods and model-based clustering analysis were used to supplement pairwise coalescence-based analysis of divergence.

KEY RESULTS

Estimated divergence times varied from 130,000 generations between North American and European subspecies to 39,000 generations between central European and Scandinavian populations. In concordance with previous studies, the highest level of diversity was found in Central Europe and the lowest in North America and a diverged Russian Karhumäki population. Local adaptation among Northern and central European populations has emerged during the last 39,000 generations. Populations that are reproductively isolated by prezygotic mechanisms have been separated for a longer time period of ∼70,000 generations but still have shared nucleotide polymorphism.

CONCLUSIONS

In A. lyrata, reproductively isolated populations started to diverge ∼70,000 generations ago and more closely related, locally adapted populations have been separate lineages for ∼39,000 generations. However, based on the posterior distribution of divergence times, the processes leading to reproductive isolation and local adaptation are likely to temporally coincide.

Demography and speciation history of the homoploid hybrid pine Pinus densata on the Tibetan Plateau

Pinus densata is an ecologically successful homoploid hybrid that inhabits vast areas of heterogeneous terrain on the south-eastern Tibetan Plateau as a result of multiple waves of colonization. Its region of origin, route of colonization onto the plateau and the directions of introgression with its parental species have previously been defined, but little is known about the isolation and divergence history of its populations. In this study, we surveyed nucleotide polymorphism over eight nuclear loci in 19 representative populations of P. densata and its parental species. Using this information and coalescence simulations, we assessed the historical changes in its population size, gene flow and divergence in time and space. The results indicate a late Miocene origin for P. densata associated with the recent uplift of south-eastern Tibet. The subsequent differentiation between geographical regions of this species began in the late Pliocene and was induced by regional topographical changes and Pleistocene glaciations. The ancestral P. densata population had a large effective population size but the central and western populations were established by limited founders, suggesting that there were severe bottlenecks during the westward migration out of the ancestral hybrid zone. After separating from their ancestral populations, population expansion occurred in all geographical regions especially in the western range. Gene flow in P. densata was restricted to geographically neighbouring populations, resulting in significant differentiation between regional groups. The new information on the divergence and demographic history of P. densata reported herein enhances our understanding of its speciation process on the Tibetan Plateau.

The distribution of the coalescence time and the number of pairwise nucleotide differences in a model of population divergence or speciation with an initial period of gene flow

This paper is concerned with a model of "isolation with an initial period of migration", where a panmictic ancestral population split into n descendant populations which exchanged migrants symmetrically at a constant rate for a period of time and subsequently became completely isolated. In the limit as the population split occurred an infinitely long time ago, the model becomes an "isolation after migration" model, describing completely isolated descendant populations which arose from a subdivided ancestral population. The probability density function of the coalescence time of a pair of genes and the probability distribution of the number of pairwise nucleotide differences are derived for both models. Whilst these are theoretical results of interest in their own right, they also give an exact analytical expression for the likelihood, for data consisting of the numbers of nucleotide differences between pairs of DNA sequences where each pair is at a different, independent locus. The behaviour of the distribution of the number of pairwise nucleotide differences under these models is illustrated and compared to the corresponding distributions under the "isolation with migration" and "complete isolation" models. It is shown that the distribution of the number of nucleotide differences between a pair of DNA sequences from different descendant populations in the model of "isolation with an initial period of migration" can be quite different from that under the "isolation with migration model", even if the average migration rate over time (and hence the total number of migrants) is the same in both scenarios. It is also illustrated how the results can be extended to other demographic scenarios that can be described by a combination of isolated panmictic populations and "symmetric island" models.

Evaluating the role of contracting and expanding rainforest in initiating cycles of speciation across the Isthmus of Panama

Climatic and geological changes across time are presumed to have shaped the rich biodiversity of tropical regions. However, the impact climatic drying and subsequent tropical rainforest contraction had on speciation has been controversial because of inconsistent palaeoecological and genetic data. Despite the strong interest in examining the role of climatic change on speciation in the Neotropics there has been few comparative studies, particularly, those that include non-rainforest taxa. We used bird species that inhabit humid or dry habitats that dispersed across the Panamanian Isthmus to characterize temporal and spatial patterns of speciation across this barrier. Here, we show that these two assemblages of birds exhibit temporally different speciation time patterns that supports multiple cycles of speciation. Evidence for these cycles is further corroborated by the finding that both assemblages consist of 'young' and 'old' species, despite dry habitat species pairs being geographically more distant than pairs of humid habitat species. The matrix of humid and dry habitats in the tropics not only allows for the maintenance of high species richness, but additionally this study suggests that these environments may have promoted speciation. We conclude that differentially expanding and contracting distributions of dry and humid habitats was probably an important contributor to speciation in the tropics.

Pleistocene speciation with and without gene flow in Euphaea damselflies of subtropical and tropical East Asian islands

Climatic oscillations during the Pleistocene period could have had a profound impact on the origin of tropical species by the alternation of allopatric isolation and interpopulation gene flow cycles. However, whether tropical speciation involves strictly allopatric isolation, or proceeds in the face of homogenizing gene flow, is relatively unclear. Here, we investigated geographical modes of speciation in four closely related Euphaea damselfly species endemic to the subtropical and tropical East Asian islands using coalescent analyses of a multilocus data set. The reconstructed phylogenies demonstrated distinct species status for each of the four species and the existence of two sister species pairs, Euphaea formosa/E. yayeyamana and E. decorata/E. ornata. The species divergence time of the sibling Euphaea damselflies dates back to within the last one Mya of the Middle to Lower Pleistocene. The speciation between the populous E. formosa of Taiwan and the less numerous E. yayeyamana of the Yaeyama islands occurred despite significant bidirectional, asymmetric gene flow, which is strongly inconsistent with a strictly allopatric model. In contrast, speciation of the approximately equal-sized populations of E. decorata of the southeast Asian mainland and E. ornata of Hainan is inferred to have involved allopatric divergence without gene flow. Our findings suggest that differential selection of natural or sexual environments is a prominent driver of species divergence in subtropical E. formosa and E. yayeyamana; whereas for tropical E. decorata and E. ornata at lower latitudes, allopatric isolation may well be a pivotal promoter of species formation.

Population genetics and objectivity in species diagnosis

Species as evolutionary lineages are expected to show greater evolutionary independence from one another than are populations within species. Two measures of evolutionary independence that stem from the study of isolation-with-migration models, one reflecting the amount of gene exchange and one reflecting the time of separation, were drawn from the literature for a large number of pairs of closely related species and pairs of populations within species. Both measures, for gene flow and time, showed broadly overlapping distributions for pairs of species and for pairs of populations within species. Species on average show more time and less gene flow than populations, but the similarity of the distributions argues against there being a qualitative difference associated with species status, as compared to populations. The two measures of evolutionary independence were similarly correlated with F(ST) estimates, which in turn also showed similar distributions for species comparisons relative to population comparisons. The measures of gene flow and separation time were examined for the capacity to discriminate intraspecific differences from interspecific differences. If used together, the two measures could be used to develop an objective (in the sense of being repeatable) measure for species diagnosis.

Pleistocene climatic oscillations and the speciation history of an alpine endemic and a widespread arctic-alpine plant

Even in cases in which geographic isolation appears to have driven the speciation of regional endemics, range shifts during the Pleistocene climatic oscillations may also have influenced their evolutionary history. Elucidating speciation history can provide novel insights into evolutionary dynamics following climatic oscillations. We demonstrated a sister relationship between the Japanese alpine endemic Cardamine nipponica and the currently allopatric, widespread arctic-alpine Cardamine bellidifolia (Brassicaceae) based on internal transcribed spacer (ITS) sequences and 10 other nuclear genes. Speciation history was inferred using demographic parameters under the isolation with migration model. The estimated demographic parameters showed that the population size of C. nipponica was similar to that of C. bellidifolia and that gene flow occurred exclusively from C. nipponica to C. bellidifolia after speciation. The inferred speciation history, which included gene flow, suggests that geographic barriers between the peripheral C. nipponica and the widespread C. bellidifolia were reduced during the Pleistocene. The asymmetric introgression implies that genetic isolation may have been involved in the speciation of C. nipponica. Our results suggest that even currently allopatric species may not have diverged solely under geographic isolation, and that their evolutionary history may have been influenced by Pleistocene range dynamics.

Molecular phylogenetics: principles and practice

Phylogenies are important for addressing various biological questions such as relationships among species or genes, the origin and spread of viral infection and the demographic changes and migration patterns of species. The advancement of sequencing technologies has taken phylogenetic analysis to a new height. Phylogenies have permeated nearly every branch of biology, and the plethora of phylogenetic methods and software packages that are now available may seem daunting to an experimental biologist. Here, we review the major methods of phylogenetic analysis, including parsimony, distance, likelihood and Bayesian methods. We discuss their strengths and weaknesses and provide guidance for their use.

Fish migrate underground: the example of Delminichthys adspersus (Cyprinidae)

Complex aquatic systems of karst harbour a rich but little-investigated biodiversity. In Croatia and Bosnia-Herzegovina karst, temporal springs are inhabited by a group of minnow-like fishes that retreat to the associated ground water during dry seasons and spend several months underground. The most abundant species in this group is Delminichthys adspersus (Heckel 1843), which also has the most fragmented distribution range. To determine the population composition and dispersal patterns, and to detect potential underground migration, a large genetic data set comprising 544 specimens of D. adspersus covering most of its distribution area was analysed. Analysis of mitochondrial cytochrome b sequences (∼1000 bp) and eight microsatellite loci showed that D. adspersus comprises at least three subpopulations with gene flow occurring among them. Coalescent-based analysis revealed a complex migration pattern, with several unidirectional dispersal paths, including between temporal springs that share no surface connection. The results of this study suggest the existence of recurrent underground migration of fish in a karst environment and demonstrate the complexity of its hydrological network. The findings are relevant to conservation strategies for endemic karst organisms and karst ecosystems as a whole.

Heterogeneity in genetic diversity among non-coding loci fails to fit neutral coalescent models of population history

Inferring aspects of the population histories of species using coalescent analyses of non-coding nuclear DNA has grown in popularity. These inferences, such as divergence, gene flow, and changes in population size, assume that genetic data reflect simple population histories and neutral evolutionary processes. However, violating model assumptions can result in a poor fit between empirical data and the models. We sampled 22 nuclear intron sequences from at least 19 different chromosomes (a genomic transect) to test for deviations from selective neutrality in the gadwall (Anas strepera), a Holarctic duck. Nucleotide diversity among these loci varied by nearly two orders of magnitude (from 0.0004 to 0.029), and this heterogeneity could not be explained by differences in substitution rates alone. Using two different coalescent methods to infer models of population history and then simulating neutral genetic diversity under these models, we found that the observed among-locus heterogeneity in nucleotide diversity was significantly higher than expected for these simple models. Defining more complex models of population history demonstrated that a pre-divergence bottleneck was also unlikely to explain this heterogeneity. However, both selection and interspecific hybridization could account for the heterogeneity observed among loci. Regardless of the cause of the deviation, our results illustrate that violating key assumptions of coalescent models can mislead inferences of population history.

Pleistocene speciation in the genus Populus (salicaceae)

The macroevolutionary consequences of recent climate change remain controversial, and there is little paleobotanical or morphological evidence that Pleistocene (1.8-0.12 Ma) glacial cycles acted as drivers of speciation, especially among lineages with long generation times, such as trees. We combined genetic and ecogeographic data from 2 closely related North American tree species, Populus balsamifera and P. trichocarpa (Salicacaeae), to determine if their divergence coincided with and was possibly caused by Pleistocene climatic events. We analyzed 32 nuclear loci from individuals of P. balsamifera and P. trichocarpa to produce coalescent-based estimates of the divergence time between the 2 species. We coupled the coalescent analyses with paleodistribution models to assess the influence of climate change on species' range. Furthermore, measures of niche overlap were used to investigate patterns of ecological differentiation between species. We estimated the divergence date of P. balsamifera and P. trichocarpa at approximately 75 Ka, which corresponds closely with the onset of Marine Isotope Stage 4 (∼76 Ka) and a rapid increase in global ice volume. Significance tests of niche overlap, in conjunction with genetic estimates of migration, suggested that speciation occurred in allopatry, possibly resulting from the environmental effects of Pleistocene glacial cycles. Our results indicate that the divergence of keystone tree species, which have shaped community diversity in northern North American ecosystems, was recent and may have been a consequence of Pleistocene-era glaciation and climate change.