DIYABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data

Genetic footprint of population fragmentation and contemporary collapse in a freshwater cetacean

Understanding demographic trends and patterns of gene flow in an endangered species is crucial for devising conservation strategies. Here, we examined the extent of population structure and recent evolution of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis). By analysing genetic variation at the mitochondrial and nuclear microsatellite loci for 148 individuals, we identified three populations along the Yangtze River, each one connected to a group of admixed ancestry. Each population displayed extremely low genetic diversity, consistent with extremely small effective size (≤92 individuals). Habitat degradation and distribution gaps correlated with highly asymmetric gene-flow that was inefficient in maintaining connectivity between populations. Genetic inferences of historical demography revealed that the populations in the Yangtze descended from a small number of founders colonizing the river from the sea during the last Ice Age. The colonization was followed by a rapid population split during the last millennium predating the Chinese Modern Economy Development. However, genetic diversity showed a clear footprint of population contraction over the last 50 years leaving only ~2% of the pre-collapsed size, consistent with the population collapses reported from field studies. This genetic perspective provides background information for devising mitigation strategies to prevent this species from extinction.

Inferring the demographic history of an oligophagous grasshopper: Effects of climatic niche stability and host-plant distribution

Understanding the consequences of past environmental changes on the abiotic and biotic components of the landscape and deciphering their impacts on the demographic trajectories of species is a major issue in evolutionary biogeography. In this study, we combine nuclear and mitochondrial genetic data to study the phylogeographical structure and lineage-specific demographic histories of the scrub-legume grasshopper (Chorthippus binotatus binotatus), a montane taxon distributed in the Iberian Peninsula and France that exclusively feeds on certain scrub-legume species. Genetic data and paleo-distribution modelling indicate the presence of four main lineages that seem to have diverged in allopatry and long-term persisted in Iberian and French refugia since the Mid Pleistocene. Comparisons of different demographic hypotheses in an Approximate Bayesian Computation (ABC) framework supported a population bottleneck in the northwestern French clade and paleo-distribution modelling indicate that the populations of this lineage have experienced more severe environmental fluctuations during the last 21 000 years than those from the Iberian Peninsula. Accordingly, we found that nuclear genetic diversity of the populations of scrub-legume grasshopper is positively associated with local stability of suitable habitats defined by both Pleistocene climate changes and historical distributional shifts of host-plant species. Overall, our study highlights the importance of integrating the potential effects of abiotic (i.e. climate and geography) and biotic components (i.e. inter-specific interactions) into the study of the evolutionary and demographic history of specialist taxa with narrow ecological requirements.

Genetic variation in Tertiary relics: The case of eastern-Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.

Genomic signatures of evolution in Nautilus-An endangered living fossil

Living fossils are survivors of previously more diverse lineages that originated millions of years ago and persisted with little morphological change. Therefore, living fossils are model organisms to study both long-term and ongoing adaptation and speciation processes. However, many aspects of living fossil evolution and their persistence in the modern world remain unclear. Here, we investigate three major aspects of the evolutionary history of living fossils: cryptic speciation, population genetics and effective population sizes, using members of the genera Nautilus and Allonautilus as classic examples of true living fossils. For this, we analysed genomewide ddRAD-Seq data for all six currently recognized nautiloid species throughout their distribution range. Our analyses identified three major allopatric Nautilus clades: a South Pacific clade, subdivided into three subclades with no signs of admixture between them; a Coral Sea clade, consisting of two genetically distinct populations with significant admixture; and a widespread Indo-Pacific clade, devoid of significant genetic substructure. Within these major clades, we detected five Nautilus groups, which likely correspond to five distinct species. With the exception of Nautilus macromphalus, all previously described species are at odds with genomewide data, testifying to the prevalence of cryptic species among living fossils. Detailed F_ST analyses further revealed significant genome-wide and locus-specific signatures of selection between species and differentiated populations, which is demonstrated here for the first time in a living fossil. Finally, approximate Bayesian computation (ABC) simulations suggest large effective population sizes, which may explain the low levels of population differentiation commonly observed in living fossils.

Testing Domestication Scenarios of Lima Bean (Phaseolus lunatus L.) in Mesoamerica: Insights from Genome-Wide Genetic Markers

Plant domestication can be seen as a long-term process that involves a complex interplay among demographic processes and evolutionary forces. Previous studies have suggested two domestication scenarios for Lima bean in Mesoamerica: two separate domestication events, one from gene pool MI in central-western Mexico and another one from gene pool MII in the area Guatemala-Costa Rica, or a single domestication from gene pool MI in central-western Mexico followed by post-domestication gene flow with wild populations. In this study we evaluated the genetic structure of the wild gene pool and tested these two competing domestication scenarios of Lima bean in Mesoamerica by applying an ABC approach to a set of genome-wide SNP markers. The results confirm the existence of three gene pools in wild Lima bean, two Mesoamerican gene pools (MI and MII) and the Andean gene pool (AI), and suggest the existence of another gene pool in central Colombia. The results indicate that although both domestication scenarios may be supported by genetic data, higher statistical support was given to the single domestication scenario in central-western Mexico followed by admixture with wild populations. Domestication would have involved strong founder effects reflected in loss of genetic diversity and increased LD levels in landraces. Genomic regions affected by selection were detected and these may harbor candidate genes related to domestication.

Global invasion network of the brown marmorated stink bug, Halyomorpha halys

Human mediated transportation into novel habitats is a prerequisite for the establishment of non-native species that become invasive, so knowledge of common sources may allow prevention. The brown marmorated stink bug (BMSB, Halyomorpha halys) is an East Asian species now established across North America and Europe, that in the Eastern United States of America (US) and Italy is causing significant economic losses to agriculture. After US populations were shown to originate from Northern China, others have tried to source BMSB populations now in Canada, Switzerland, Italy, France, Greece, and Hungary. Due to selection of different molecular markers, however, integrating all the datasets to obtain a broader picture of BMSB’s expansion has been difficult. To address this limitation we focused on a single locus, the barcode region in the cytochrome oxidase I mitochondrial gene, and analyzed representative BMSB samples from across its current global range using an Approximate Bayesian Computation approach. We found that China is the likely source of most non-native populations, with at least four separate introductions in North America and three in Europe. Additionally, we found evidence of one bridgehead event: a likely Eastern US source for the central Italy populations that interestingly share enhanced pest status.

Analysis of spatial genetic variation reveals genetic divergence among populations of Primula veris associated to contrasting habitats

Genetic divergence by environment is a process whereby selection causes the formation of gene flow barriers between populations adapting to contrasting environments and is often considered to be the onset of speciation. Nevertheless, the extent to which genetic differentiation by environment on small spatial scales can be detected by means of neutral markers is still subject to debate. Previous research on the perennial herb Primula veris has shown that plants from grassland and forest habitats showed pronounced differences in phenology and flower morphology, suggesting limited gene flow between habitats. To test this hypothesis, we sampled 33 populations of P. veris consisting of forest and grassland patches and used clustering techniques and network analyses to identify sets of populations that are more connected to each other than to other sets of populations and estimated the timing of divergence. Our results showed that spatial genetic variation had a significantly modular structure and consisted of four well-defined modules that almost perfectly coincided with habitat features. Genetic divergence was estimated to have occurred about 114 generations ago, coinciding with historic major changes in the landscape. Overall, these results illustrate how populations adapting to different environments become structured genetically within landscapes on small spatial scales.

When new human-modified habitats favour the expansion of an amphibian pioneer species: Evolutionary history of the natterjack toad (Bufo calamita) in a coal basin

Human activities affect microevolutionary dynamics by inducing environmental changes. In particular, land cover conversion and loss of native habitats decrease genetic diversity and jeopardize the adaptive ability of populations. Nonetheless, new anthropogenic habitats can also promote the successful establishment of emblematic pioneer species. We investigated this issue by examining the population genetic features and evolutionary history of the natterjack toad (Bufo [Epidalea] calamita) in northern France, where populations can be found in native coastal habitats and coalfield habitats shaped by European industrial history, along with an additional set of European populations located outside this focal area. We predicted contrasting patterns of genetic structure, with newly settled coalfield populations departing from migration-drift equilibrium. As expected, coalfield populations showed a mosaic of genetically divergent populations with short-range patterns of gene flow, and native coastal populations indicated an equilibrium state with an isolation-by-distance pattern suggestive of postglacial range expansion. However, coalfield populations exhibited (i) high levels of genetic diversity, (ii) no evidence of local inbreeding or reduced effective population size and (iii) multiple maternal mitochondrial lineages, a genetic footprint depicting independent colonization events. Furthermore, approximate Bayesian computations suggested several evolutionary trajectories from ancient isolation in glacial refugia during the Pleistocene, with biogeographical signatures of recent expansion probably confounded by human-mediated mixing of different lineages. From an evolutionary and conservation perspective, this study highlights the ecological value of industrial areas, provided that ongoing regional gene flow is ensured within the existing lineage boundaries.

Tracking the return of Aedes aegypti to Brazil, the major vector of the dengue, chikungunya and Zika viruses

BACKGROUND

Aedes aegypti, commonly known as "the yellow fever mosquito", is of great medical concern today primarily as the major vector of dengue, chikungunya and Zika viruses, although yellow fever remains a serious health concern in some regions. The history of Ae. aegypti in Brazil is of particular interest because the country was subjected to a well-documented eradication program during 1940s-1950s. After cessation of the campaign, the mosquito quickly re-established in the early 1970s with several dengue outbreaks reported during the last 30 years. Brazil can be considered the country suffering the most from the yellow fever mosquito, given the high number of dengue, chikungunya and Zika cases reported in the country, after having once been declared "free of Ae. aegypti".

METHODOLOGY/PRINCIPAL FINDINGS

We used 12 microsatellite markers to infer the genetic structure of Brazilian Ae. aegypti populations, genetic variability, genetic affinities with neighboring geographic areas, and the timing of their arrival and spread. This enabled us to reconstruct their recent history and evaluate whether the reappearance in Brazil was the result of re-invasion from neighboring non-eradicated areas or re-emergence from local refugia surviving the eradication program. Our results indicate a genetic break separating the northern and southern Brazilian Ae. aegypti populations, with further genetic differentiation within each cluster, especially in southern Brazil.

CONCLUSIONS/SIGNIFICANCE

Based on our results, re-invasions from non-eradicated regions are the most likely scenario for the reappearance of Ae. aegypti in Brazil. While populations in the northern cluster are likely to have descended from Venezuela populations as early as the 1970s, southern populations seem to have derived more recently from northern Brazilian areas. Possible entry points are also revealed within both southern and northern clusters that could inform strategies to control and monitor this important arbovirus vector.

Deciphering the Routes of invasion of Drosophila suzukii by Means of ABC Random Forest

Deciphering invasion routes from molecular data is crucial to understanding biological invasions, including identifying bottlenecks in population size and admixture among distinct populations. Here, we unravel the invasion routes of the invasive pest Drosophila suzukii using a multi-locus microsatellite dataset (25 loci on 23 worldwide sampling locations). To do this, we use approximate Bayesian computation (ABC), which has improved the reconstruction of invasion routes, but can be computationally expensive. We use our study to illustrate the use of a new, more efficient, ABC method, ABC random forest (ABC-RF) and compare it to a standard ABC method (ABC-LDA). We find that Japan emerges as the most probable source of the earliest recorded invasion into Hawaii. Southeast China and Hawaii together are the most probable sources of populations in western North America, which then in turn served as sources for those in eastern North America. European populations are genetically more homogeneous than North American populations, and their most probable source is northeast China, with evidence of limited gene flow from the eastern US as well. All introduced populations passed through bottlenecks, and analyses reveal five distinct admixture events. These findings can inform hypotheses concerning how this species evolved between different and independent source and invasive populations. Methodological comparisons indicate that ABC-RF and ABC-LDA show concordant results if ABC-LDA is based on a large number of simulated datasets but that ABC-RF out-performs ABC-LDA when using a comparable and more manageable number of simulated datasets, especially when analyzing complex introduction scenarios.

Genetic structure and gene flow of the flea Xenopsylla cheopis in Madagascar and Mayotte

BACKGROUND

The flea Xenopsylla cheopis (Siphonaptera: Pulicidae) is a vector of plague. Despite this insect's medical importance, especially in Madagascar where plague is endemic, little is known about the organization of its natural populations. We undertook population genetic analyses (i) to determine the spatial genetic structure of X. cheopis in Madagascar and (ii) to determine the potential risk of plague introduction in the neighboring island of Mayotte.

RESULTS

We genotyped 205 fleas from 12 sites using nine microsatellite markers. Madagascan populations of X. cheopis differed, with the mean number of alleles per locus per population ranging from 1.78 to 4.44 and with moderate to high levels of genetic differentiation between populations. Three distinct genetic clusters were identified, with different geographical distributions but with some apparent gene flow between both islands and within Malagasy regions. The approximate Bayesian computation (ABC) used to test the predominant direction of flea dispersal implied a recent population introduction from Mayotte to Madagascar, which was estimated to have occurred between 1993 and 2012. The impact of this flea introduction in terms of plague transmission in Madagascar is unclear, but the low level of flea exchange between the two islands seems to keep Mayotte free of plague for now.

CONCLUSION

This study highlights the occurrence of genetic structure among populations of the flea vector of plague, X. cheopis, in Madagascar and suggests that a flea population from Mayotte has been introduced to Madagascar recently. As plague has not been reported in Mayotte, this introduction is unlikely to present a major concern for plague transmission. Nonetheless, evidence of connectivity among flea populations in the two islands indicates a possibility for dispersal by fleas in the opposite direction and thus a risk of plague introduction to Mayotte.

Demographic inference through approximate-Bayesian-computation skyline plots

The skyline plot is a graphical representation of historical effective population sizes as a function of time. Past population sizes for these plots are estimated from genetic data, without a priori assumptions on the mathematical function defining the shape of the demographic trajectory. Because of this flexibility in shape, skyline plots can, in principle, provide realistic descriptions of the complex demographic scenarios that occur in natural populations. Currently, demographic estimates needed for skyline plots are estimated using coalescent samplers or a composite likelihood approach. Here, we provide a way to estimate historical effective population sizes using an Approximate Bayesian Computation (ABC) framework. We assess its performance using simulated and actual microsatellite datasets. Our method correctly retrieves the signal of contracting, constant and expanding populations, although the graphical shape of the plot is not always an accurate representation of the true demographic trajectory, particularly for recent changes in size and contracting populations. Because of the flexibility of ABC, similar approaches can be extended to other types of data, to multiple populations, or to other parameters that can change through time, such as the migration rate.

Genetic variation of the land-locked freshwater shrimp Caridina pseudodenticulata (Decapoda: Atyidae: Caridina) in Taiwan

Caridina pseudodenticulata is a land-locked freshwater shrimp, which is endemic to North and West Taiwan. It is listed as Vulnerable under criteria A1ce in the IUCN Red List. Our study used mitochondrial DNA fragment sequences of the cytochrome oxidase subunit I (COI) gene to examine its genetic structure. In total, 204 sequences were analyzed from 22 populations, and all 39 haplotypes were identified. Our study (1) found that the genetic diversity within population of this species was similar to that of other land-locked freshwater prawn in Taiwan; (2) displayed that during glaciation, the Taiwan Strait was largely above water, an event which created opportunities for colonization Taiwan and dispersal widely; (3) showed that although the phylogenetic analysis showed the lack of a population genetic structure, the AMOVA and haplotypes distribution pattern revealed geographically divided; (4) indicated that C. pseudodenticulata could across to neighbor rivers through the mouths of rivers, which confluent with each other during floods.

Genetic Patterns of Myrceugenia correifolia, a Rare Species of Fog-Dependent Forests of Mediterranean Chile: Is It a Climatic Relict?

Rare species frequently occur in areas with microclimatic conditions that are atypical for their regions, but that were more common in the past, and that probably have operated as climatic refugia for a long time. Myrceugenia correifolia is a rare arboreal species that grows in deep canyons and hilltops of the Coast Range of north-central Chile between 30° and 35°S. In the northern edge of its distribution M. correifolia grows in small patches of fog-dependent forest surrounding by xeric vegetation. These forest formations are thought to be remnants of an ancient and continuous rainforest that according to some authors became fragmented during aridization of the Neogene (Neogene relict) and to others during warm-dry cycles of the Pleistocene (glacial relicts). Here we asked whether the northernmost populations of M. correifolia are Neogene relicts, glacial relicts, or the result of a recent northward colonization. To answer this question we examined genetic diversity and population divergence of M. correifolia using microsatellite markers, tested various competing population history scenarios with an approximate Bayesian computation (ABC) method, and complemented these data with ecological niche modeling (ENM). We detected three genetic clusters with a distinctive latitudinal pattern (north, center, and south) and high levels of differentiation (F_ST = 0.36). Demographic inference supported an admixture event 31 kya between two populations that diverged from an ancient population 139 kya. The admixture time coincides with the beginning of a period of wet conditions in north-central Chile that extended from 33 to 19 kya and was preceded by dry and cold conditions. These results suggest that increased precipitation during glacial periods triggered northward expansion of the range of M. correifolia, with subsequent admixture between populations that remained separated during interglacial periods. Accordingly, ENM models showed that suitable habitats for M. correifolia in north-central Chile were larger and less fragmented during the Last Glacial Maximum than at present, suggesting that northernmost populations of this species are glacial relicts.

The eastern part of the Fertile Crescent concealed an unexpected route of olive (Olea europaea L.) differentiation

BACKGROUND AND AIMS

Olive is considered a native plant of the eastern side of the Mediterranean basin, from where it should have spread westward along the Mediterranean shores, while little is known about its diffusion in the eastern direction.

METHODS

Genetic diversity levels and population genetic structure of a wide set of olive ecotypes and varieties collected from several provinces of Iran, representing a high percentage of the entire olive resources present in the area, was screened with 49 chloroplast and ten nuclear simple sequence repeat markers, and coupled with archaeo-botanical and historical data on Mediterranean olive varieties. Approximate Bayesian Computation was applied to define the demographic history of olives including Iranian germplasm, and species distribution modelling was performed to understand the impact of the Late Quaternary on olive distribution.

KEY RESULTS

The results of the present study demonstrated that: (1) the climatic conditions of the last glacial maximum had an important role on the actual olive distribution, (2) all Iranian olive samples had the same maternal inheritance as Mediterranean cultivars, and (3) the nuclear gene flow from the Mediterranean basin to the Iranian plateau was almost absent, as well as the contribution of subspecies cuspidata to the diversity of Iranian olives.

CONCLUSIONS

Based on this evidence, a new scenario for the origin and distribution of this important fruit crop has been traced. The evaluation of olive trees growing in the eastern part of the Levant highlighted a new perspective on the spread and distribution of olive, suggesting two routes of olive differentiation, one westward, spreading along the Mediterranean basin, and another moving towards the east and reaching the Iranian plateau before its domestication.

A Mesoamerican origin of cherimoya (Annona cherimola Mill.): Implications for the conservation of plant genetic resources

Knowledge on the structure and distribution of genetic diversity is a key aspect to plan and execute an efficient conservation and utilization of the genetic resources of any crop as well as for determining historical demographic inferences. In this work, a large data set of 1,765 accessions of cherimoya (Annona cherimola Mill, Annonaceae), an underutilized fruit tree crop native to the Neotropics and used as a food source by pre-Columbian cultures, was collected from six different countries across the American continent and amplified with nine highly informative microsatellite markers. The structure analyses, fine representation of the genetic diversity and an ABC approach suggest a Mesoamerican origin of the crop, contrary to previous reports, with clear implications for the dispersion of plant germplasm between Central and South America in pre-Columbian times. These results together with the potential distribution of the species in a climatic change context using two different climate models provide new insights for the history and conservation of extant genetic resources of cherimoya that can be applied to other currently underutilized woody perennial crops.

Mitogenomics data reveal effective population size, historical bottlenecks, and the effects of hunting on New Zealand fur seals (Arctocephalus forsteri)

The New Zealand fur seal (Arctocephalus forsteri) passed through a population bottleneck due to commercial sealing during the eighteenth to nineteenth centuries. To facilitate future management options, we reconstructed the demographic history of New Zealand fur seals in a Bayesian framework using maternally inherited, mitochondrial DNA sequences. Mitogenomic data suggested two separate clades (most recent common ancestor 5000 years ago) of New Zealand fur seals that survived large-scale human harvest. Mitochondrial haplotype diversity was high, with 45 singletons identified from 46 individuals although mean nucleotide diversity was low (0.012 ± 0.0061). Variation was not constrained geographically. Analyses of mitogenomes support the hypothesis for a population bottleneck approximately 35 generations ago, which coincides with the peak of commercial sealing. Mitogenomic data are consistent with a pre-human effective population size of approximately 30,000 that first declined to around 10,000 (due to the impact of Polynesian colonization, particularly in the first 100 years of their arrival into New Zealand), and then to 100-200 breeding individuals during peak of commercial sealing.

Extraordinarily rapid speciation in a marine fish

Divergent selection may initiate ecological speciation extremely rapidly. How often and at what pace ecological speciation proceeds to yield strong reproductive isolation is more uncertain. Here, we document a case of extraordinarily rapid speciation associated with ecological selection in the postglacial Baltic Sea. European flounders (Platichthys flesus) in the Baltic exhibit two contrasting reproductive behaviors: pelagic and demersal spawning. Demersal spawning enables flounders to thrive in the low salinity of the Northern Baltic, where eggs cannot achieve neutral buoyancy. We show that demersal and pelagic flounders are a species pair arising from a recent event of speciation. Despite having a parapatric distribution with extensive overlap, the two species are reciprocally monophyletic and show strongly bimodal genotypic clustering and no evidence of contemporary migration, suggesting strong reproductive isolation. Divergence across the genome is weak but shows strong signatures of selection, a pattern suggestive of a recent ecological speciation event. We propose that spawning behavior in Baltic flounders is the trait under ecologically based selection causing reproductive isolation, directly implicating a process of ecological speciation. We evaluated different possible evolutionary scenarios under the approximate Bayesian computation framework and estimate that the speciation process started in allopatry ∼2,400 generations ago, following the colonization of the Baltic by the demersal lineage. This is faster than most known cases of ecological speciation and represents the most rapid event of speciation ever reported for any marine vertebrate.

Ecotypic differentiation under farmers' selection: Molecular insights into the domestication of Pachyrhizus Rich. ex DC. (Fabaceae) in the Peruvian Andes

Understanding the distribution of crop genetic diversity in relation to environmental factors can give insights into the eco‐evolutionary processes involved in plant domestication. Yam beans (Pachyrhizus Rich. ex DC.) are leguminous crops native to South and Central America that are grown for their tuberous roots but are seed‐propagated. Using a landscape genetic approach, we examined correlations between environmental factors and phylogeographic patterns of genetic diversity in Pachyrhizus landrace populations. Molecular analyses based on chloroplast DNA sequencing and a new set of nuclear microsatellite markers revealed two distinct lineages, with strong genetic differentiation between Andean landraces (lineage A) and Amazonian landraces (lineage B). The comparison of different evolutionary scenarios for the diversification history of yam beans in the Andes using approximate Bayesian computation suggests that Pachyrhizus ahipa and Pachyrhizus tuberosus share a progenitor‐derivative relationship, with environmental factors playing an important role in driving selection for divergent ecotypes. The new molecular data call for a revision of the taxonomy of Pachyrhizus but are congruent with paleoclimatic and archeological evidence, and suggest that selection for determinate growth was part of ecophysiological adaptations associated with the diversification of the P. tuberosus–P. ahipa complex during the Mid‐Holocene.

Toll-like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the 'ghost of selection past' and quantify the 'drift debt'

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll-like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent-based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the 'ghost of selection past'. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such 'drift debt' occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.

Combining phylogenetic and demographic inferences to assess the origin of the genetic diversity in an isolated wolf population

The survival of isolated small populations is threatened by both demographic and genetic factors. Large carnivores declined for centuries in most of Europe due to habitat changes, overhunting of their natural prey and direct persecution. However, the current rewilding trends are driving many carnivore populations to expand again, possibly reverting the erosion of their genetic diversity. In this study we reassessed the extent and origin of the genetic variation of the Italian wolf population, which is expanding after centuries of decline and isolation. We genotyped wolves from Italy and other nine populations at four mtDNA regions (control-region, ATP6, COIII and ND4) and 39 autosomal microsatellites. Results of phylogenetic analyses and assignment procedures confirmed in the Italian wolves a second private mtDNA haplotype, which belongs to a haplogroup distributed mostly in southern Europe. Coalescent analyses showed that the unique mtDNA haplotypes in the Italian wolves likely originated during the late Pleistocene. ABC simulations concordantly showed that the extant wolf populations in Italy and in south-western Europe started to be isolated and declined right after the last glacial maximum. Thus, the standing genetic variation in the Italian wolves principally results from the historical isolation south of the Alps.

Shifting Quaternary migration patterns in the Bahamian archipelago: Evidence from the Zamia pumila complex at the northern limits of the Caribbean island biodiversity hotspot

PREMISE OF THE STUDY

The Bahamas archipelago is formed by young, tectonically stable carbonate banks that harbor direct geological evidence of global ice-volume changes. We sought to detect signatures of major changes on gene flow patterns and reconstruct the phylogeographic history of the monophyletic Zamia pumila complex across the Bahamas.

METHODS

Nuclear molecular markers with both high and low mutation rates were used to capture two different time scale signatures and test several gene flow and demographic hypotheses.

KEY RESULTS

Single-copy nuclear genes unveiled apparent ancestral admixture on Andros, suggesting a significant role of this island as main hub of diversity of the archipelago. We detected demographic and spatial expansion of the Zamia pumila complex on both paleo-provinces around the Piacenzian (Pliocene)/Gelasian (Pleistocene). Populations evidenced signatures of different migration models that have occurred at two different times. Populations on Long Island (Z. lucayana) may either represent a secondary colonization of the Bahamas by Zamia or a rapid and early-divergence event of at least one population on the Bahamas.

CONCLUSIONS

Despite changes in migration patterns with global climate, expected heterozygosity with both marker systems remains within the range reported for cycads, but with significant levels of increased inbreeding detected by the microsatellites. This finding is likely associated with reduced gene flow between and within paleo-provinces, accompanied by genetic drift, as rising seas enforced isolation. Our study highlights the importance of the maintenance of the predominant direction of genetic exchange and the role of overseas dispersion among the islands during climate oscillations.

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub-legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west-to-east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.

Genetic evidence challenges the native status of a threatened freshwater fish (Carassius carassius) in England

A fundamental consideration for the conservation of a species is the extent of its native range, that is, regions naturally colonized. However, both natural processes and human-mediated introductions can drive species distribution shifts. Ruling out the human-mediated introduction of a species into a given region is vital for its conservation, but remains a significant challenge in most cases. The crucian carp Carassius carassius (L.) is a threatened freshwater fish thought to be native to much of Europe. However, its native status in England is based only on anecdotal evidence. Here, we devise an approach that can be used to empirically test the native status of English fauna. We use this approach, along with 13 microsatellite loci, population structure analyses, and Approximate Bayesian Computation (ABC), to test hypotheses for the origins of C. carassius in England. Contrary to the current consensus, we find strong support for the human-mediated introduction of C. carassius into England during the 15th century. This result stimulates an interesting and timely debate surrounding motivations for the conservation of species. We discuss this topic, and the potential for continued conservation of C. carassius in England, despite its non-native origins.

The role of the Strait of Gibraltar in shaping the genetic structure of the Mediterranean Grenadier, Coryphaenoides mediterraneus, between the Atlantic and Mediterranean Sea

Population genetic studies of species inhabiting the deepest parts of the oceans are still scarce and only until recently we started to understand how oceanographic processes and topography affect dispersal and gene flow patterns. The aim of this study was to investigate the spatial population genetic structure of the bathyal bony fish Coryphaenoides mediterraneus, with a focus on the Atlantic-Mediterranean transition. We used nine nuclear microsatellites and the mitochondrial cytochrome c oxidase I gene from 6 different sampling areas. No population genetic structure was found within Mediterranean with both marker types (mean ΦST = 0.0960, FST = -0.0003, for both P > 0.05). However, within the Atlantic a contrasting pattern of genetic structure was found for the mtDNA and nuclear markers (mean ΦST = 0.2479, P < 0.001; FST = -0.0001, P > 0.05). When comparing samples from Atlantic and Mediterranean they exhibited high and significant levels of genetic divergence (mean ΦST = 0.7171, FST = 0.0245, for both P < 0.001) regardless the genetic marker used. Furthermore, no shared haplotypes were found between Atlantic and Mediterranean populations. These results suggest very limited genetic exchange between Atlantic and Mediterranean populations of C. mediterraneus, likely due to the shallow bathymetry of the Strait of Gibraltar acting as a barrier to gene flow. This physical barrier not only prevents the direct interactions between the deep-living adults, but also must prevent interchange of pelagic early life stages between the two basins. According to Bayesian simulations it is likely that Atlantic and Mediterranean populations of C. mediterraneus were separated during the late Pleistocene, which is congruent with results for other deep-sea fish from the same region.

Alternative glacial-interglacial refugia demographic hypotheses tested on Cephalocereus columna-trajani (Cactaceae) in the intertropical Mexican drylands

Historic demography changes of plant species adapted to New World arid environments could be consistent with either the Glacial Refugium Hypothesis (GRH), which posits that populations contracted to refuges during the cold-dry glacial and expanded in warm-humid interglacial periods, or with the Interglacial Refugium Hypothesis (IRH), which suggests that populations contracted during interglacials and expanded in glacial times. These contrasting hypotheses are developed in the present study for the giant columnar cactus Cephalocereus columna-trajani in the intertropical Mexican drylands where the effects of Late Quaternary climatic changes on phylogeography of cacti remain largely unknown. In order to determine if the historic demography and phylogeographic structure of the species are consistent with either hypothesis, sequences of the chloroplast regions psbA-trnH and trnT-trnL from 110 individuals from 10 populations comprising the full distribution range of this species were analysed. Standard estimators of genetic diversity and structure were calculated. The historic demography was analysed using a Bayesian approach and the palaeodistribution was derived from ecological niche modelling to determine if, in the arid environments of south-central Mexico, glacial-interglacial cycles drove the genetic divergence and diversification of this species. Results reveal low but statistically significant population differentiation (F_ST = 0.124, P < 0.001), although very clear geographic clusters are not formed. Genetic diversity, haplotype network and Approximate Bayesian Computation (ABC) demographic analyses suggest a population expansion estimated to have taken place in the Last Interglacial (123.04 kya, 95% CI 115.3–130.03). The species palaeodistribution is consistent with the ABC analyses and indicates that the potential area of palaedistribution and climatic suitability were larger during the Last Interglacial and Holocene than in the Last Glacial Maximum. Overall, these results suggest that C. columna-trajani experienced an expansion following the warm conditions of interglacials, in accordance with the GRH.

Tiger on the prowl: Invasion history and spatio-temporal genetic structure of the Asian tiger mosquito Aedes albopictus (Skuse 1894) in the Indo-Pacific

Background

Within the last century, increases in human movement and globalization of trade have facilitated the establishment of several highly invasive mosquito species in new geographic locations with concurrent major environmental, economic and health consequences. The Asian tiger mosquito, Aedes albopictus, is an extremely invasive and aggressive daytime-biting mosquito that is a major public health threat throughout its expanding range.

Methodology/Principal findings

We used 13 nuclear microsatellite loci (on 911 individuals) and mitochondrial COI sequences to gain a better understanding of the historical and contemporary movements of Ae. albopictus in the Indo-Pacific region and to characterize its population structure. Approximate Bayesian computation (ABC) was employed to test competing historical routes of invasion of Ae. albopictus within the Southeast (SE) Asian/Australasian region. Our ABC results show that Ae. albopictus was most likely introduced to New Guinea via mainland Southeast Asia, before colonizing the Solomon Islands via either Papua New Guinea or SE Asia. The analysis also supported that the recent incursion into northern Australia’s Torres Strait Islands was seeded chiefly from Indonesia. For the first time documented in this invasive species, we provide evidence of a recently colonized population (the Torres Strait Islands) that has undergone rapid temporal changes in its genetic makeup, which could be the result of genetic drift or represent a secondary invasion from an unknown source.

Conclusions/Significance

There appears to be high spatial genetic structure and high gene flow between some geographically distant populations. The species' genetic structure in the region tends to favour a dispersal pattern driven mostly by human movements. Importantly, this study provides a more widespread sampling distribution of the species’ native range, revealing more spatial population structure than previously shown. Additionally, we present the most probable invasion history of this species in the Australasian region using ABC analysis.

The Asian tiger mosquito, Aedes albopictus, is an aggressive mosquito that has expanded globally in the last century, chiefly due to the increase of human movements. It is a major public health concern due to its role in transmitting dengue, chikungunya and Zika viruses. Understanding how populations of Ae. albopictus are genetically related and how they have been introduced into new regions is important for controlling them and assessing their disease risk; few studies have explored this in the Indo-Pacific. In our study, we sampled a broader range of populations of Ae. albopictus in the Indo-Pacific to explore genetic patterns and to investigate the likely route of invasion through Australasia. We uncovered clear genetic groups throughout this region, but also found that some geographically distant populations are closely related, likely due to human-associated movements. We also found, that Ae. albopictus likely colonized New Guinea from mainland Southeast (SE) Asia before spreading to the Solomon Islands via either PNG or SE Asia. In contrast, the populations in Australia’s Torres Strait Islands were introduced from Indonesia. Interestingly, we found major genetic changes over time in some Torres Strait populations, less than a decade after its introduction, potentially reflective of a random reduction in genetic diversity (genetic drift) or a secondary invasion.

The population genetic structure of Corythucha ciliata (Say) (Hemiptera: Tingidae) provides insights into its distribution and invasiveness

Corythucha ciliata (Say), an invasive pest originating from North America, causes severe damage on sycamore trees. However, little is known about the population genetics and evolutionary forces underlying the invasiveness of this important pest. In the present study, we use three mitochondrial genes (COI, ND1 and ND5) and nine microsatellite markers to investigate the population genetics of C. ciliata and retrace its spread through China. The results suggest a low level of genetic diversity in Chinese and European populations of C. ciliata. Our results indicate that populations of C. ciliata have obvious genetic structure, and genetic differentiation is not caused by geographic isolation. In median-joining networks, we observed a higher frequency of shared haplotypes in groups 1 and 3. Based on gene flow and approximate Bayesian computation analyses, we discovered that C. ciliata first invaded the east coast of China and subsequently moved inland. Demographic analysis suggested that populations of C. ciliata in China may have undergone a recent bottleneck effect. Finally, our results suggest that population structure, high gene flow and environmental conditions have favored the broad invasiveness of this important pest.

Postglacial recolonization shaped the genetic diversity of the winter moth (Operophtera brumata) in Europe

Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.

Genetic structure and invasion history of the house mouse (Mus musculus domesticus) in Senegal, West Africa: a legacy of colonial and contemporary times

Knowledge of the genetic make-up and demographic history of invasive populations is critical to understand invasion mechanisms. Commensal rodents are ideal models to study whether complex invasion histories are typical of introductions involving human activities. The house mouse Mus musculus domesticus is a major invasive synanthropic rodent originating from South-West Asia. It has been largely studied in Europe and on several remote islands, but the genetic structure and invasion history of this taxon have been little investigated in several continental areas, including West Africa. In this study, we focussed on invasive populations of M. m. domesticus in Senegal. In this focal area for European settlers, the distribution area and invasion spread of the house mouse is documented by decades of data on commensal rodent communities. Genetic variation at one mitochondrial locus and 16 nuclear microsatellite markers was analysed from individuals sampled in 36 sites distributed across the country. A combination of phylogeographic and population genetics methods showed that there was a single introduction event on the northern coast of Senegal, from an exogenous (probably West European) source, followed by a secondary introduction from northern Senegal into a coastal site further south. The geographic locations of these introduction sites were consistent with the colonial history of Senegal. Overall, the marked microsatellite genetic structure observed in Senegal, even between sites located close together, revealed a complex interplay of different demographic processes occurring during house mouse spatial expansion, including sequential founder effects and stratified dispersal due to human transport along major roads.

Rethinking the history of common walnut (Juglans regia L.) in Europe: Its origins and human interactions

Common walnut (Juglans regia L) is an economically important species cultivated worldwide for its high-quality wood and nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of walnut in Europe is a matter of debate, however. In this study, we estimated the genetic diversity and structure of 91 Eurasian walnut populations using 14 neutral microsatellites. By integrating fossil pollen, cultural, and historical data with population genetics, and approximate Bayesian analysis, we reconstructed the demographic history of walnut and its routes of dispersal across Europe. The genetic data confirmed the presence of walnut in glacial refugia in the Balkans and western Europe. We conclude that human-mediated admixture between Anatolian and Balkan walnut germplasm started in the Early Bronze Age, and between western Europe and the Balkans in eastern Europe during the Roman Empire. A population size expansion and subsequent decline in northeastern and western Europe was detected in the last five centuries. The actual distribution of walnut in Europe resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and its human exploitation over the last 5,000 years.

One species or four? Yes!...and, no. Or, arbitrary assignment of lineages to species obscures the diversification processes of Neotropical fishes

Species are fundamental units in many biological disciplines, but there is continuing disagreement as to what species are, how to define them, and even whether the concept is useful. While some of this debate can be attributed to inadequate data and insufficient statistical frameworks in alpha taxonomy, an equal part results from the ambiguity over what species are expected to represent by the many who use them. Here, mtDNA data, microsatellite data, and sequence data from 17 nuclear loci are used in an integrated and quantitative manner to resolve the presence of evolutionary lineages, their contemporary and historical structure, and their correspondence to species, in a species complex of Amazonian peacock "bass" cichlids (Cichla pinima sensu lato). Results suggest that the historical narrative for these populations is more complex than can be portrayed by recognizing them as one, two, or four species: their history and contemporary dynamics cannot be unambiguously rendered as discrete units (taxa) at any level without both choosing the supremacy of one delimitation criterion and obscuring the very information that provides insight into the diversification process. This calls into question the utility of species as a rank, term, or concept, and suggests that while biologists may have a reasonable grasp of the structure of evolution, our methods of conveying these insights need updating. The lack of correspondence between evolutionary phenomena and discrete species should serve as a null hypothesis, and researchers should focus on quantifying the diversity in nature at whatever hierarchical level it occurs.

Population genomic analyses reveal a history of range expansion and trait evolution across the native and invaded range of yellow starthistle (Centaurea solstitialis)

Identifying sources of genetic variation and reconstructing invasion routes for non-native introduced species is central to understanding the circumstances under which they may evolve increased invasiveness. In this study, we used genome-wide single nucleotide polymorphisms to study the colonization history of Centaurea solstitialis in its native range in Eurasia and invasions into the Americas. We leveraged this information to pinpoint key evolutionary shifts in plant size, a focal trait associated with invasiveness in this species. Our analyses revealed clear population genomic structure of potential source populations in Eurasia, including deep differentiation of a lineage found in the southern Apennine and Balkan Peninsulas and divergence among populations in Asia, eastern Europe and western Europe. We found strongest support for an evolutionary scenario in which western European populations were derived from an ancient admixture event between populations from eastern Europe and Asia, and subsequently served as the main genetic 'bridgehead' for introductions to the Americas. Introductions to California appear to be from a single source region, and multiple, independent introductions of divergent genotypes likely occurred into the Pacific Northwest. Plant size has evolved significantly at three points during range expansion, including a large size increase in the lineage responsible for the aggressive invasion of the California interior. These results reveal a long history of colonization, admixture and trait evolution in C. solstitialis, and suggest routes for improving evidence-based management decisions for one of the most ecologically and economically damaging invasive species in the western United States.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

About Ganoderma boninense in oil palm plantations of Sumatra and peninsular Malaysia: Ancient population expansion, extensive gene flow and large scale dispersion ability

Wood rot fungi form one of the main classes of phytopathogenic fungus. The group includes many species, but has remained poorly studied. Many species belonging to the Ganoderma genus are well known for causing decay in a wide range of tree species around the world. Ganoderma boninense, causal agent of oil palm basal stem rot, is responsible for considerable yield losses in Southeast Asian oil palm plantations. In a large-scale sampling operation, 357 sporophores were collected from oil palm plantations spread over peninsular Malaysia and Sumatra and genotyped using 11 SSR markers. The genotyping of these samples made it possible to investigate the population structure and demographic history of G. boninense across the oldest known area of interaction between oil palm and G. boninense. Results show that G. boninense possesses a high degree of genetic diversity and no detectable genetic structure at the scale of Sumatra and peninsular Malaysia. The fact that few duplicate genotypes were found in several studies including this one supports the hypothesis of spore dispersal in the spread of G. boninense. Meanwhile, spatial autocorrelation analysis shows that G. boninense is able to disperse across both short and long distances. These results bring new insight into mechanisms by which G. boninense spreads in oil palm plantations. Finally, the use of approximate Bayesian computation (ABC) modelling indicates that G. boninense has undergone a demographic expansion in the past, probably before the oil palm was introduced into Southeast Asia.

Genetic consequences of Quaternary climatic oscillations in the Himalayas: Primula tibetica as a case study based on restriction site-associated DNA sequencing

The effects of Quaternary climatic oscillations on the demography of organisms vary across regions and continents. In taxa distributed in Europe and North America, several paradigms regarding the distribution of refugia have been identified. By contrast, less is known about the processes that shaped the species' spatial genetic structure in areas such as the Himalayas, which is considered a biodiversity hotspot. Here, we investigated the phylogeographic structure and population dynamics of Primula tibetica by combining genomic phylogeography and species distribution models (SDMs). Genomic data were obtained for 293 samples of P. tibetica using restriction site-associated DNA sequencing (RADseq). Ensemble SDMs were carried out to predict potential present and past distribution ranges. Four distinct lineages were identified. Approximate Bayesian computation analyses showed that each of them have experienced both expansions and bottlenecks since their divergence, which occurred during or across the Quaternary glacial cycles. The two lineages at both edges of the distribution were found to be more vulnerable and responded in different ways to past climatic changes. These results illustrate how past climatic changes affected the demographic history of Himalayan organisms. Our findings highlight the significance of combining genomic approaches with environmental data when evaluating the effects of past climatic changes.

Origin and Dispersal History of Two Colonial Ascidian Clades in the Botryllus schlosseri Species Complex

Human-induced global warming and species introductions are rapidly altering the composition and functioning of Earth’s marine ecosystems. Ascidians (Phylum Chordata, Subphylum Tunicata, Class Ascidiacea) are likely to play an increasingly greater role in marine communities. The colonial ascidian B. schlosseri is a cryptic species complex comprising five genetically divergent clades (A-E). Clade A is a global species, and Clade E has so far been identified in European waters only. Using the largest mitochondrial cytochrome oxidase I datasets yet assembled, we determine the origin and dispersal history of these species. Nucleotide diversity and Approximate Bayesian Computation analyses support a Pacific origin for Clade A, with two likely dispersal scenarios that both show the northwestern Atlantic populations establishing early in the history of the species. Both Discrete Phylogeographic Analysis and Approximate Bayesian Computation support an origin of Clade E on the French side of the English Channel. An unsampled lineage evolved from the French lineage, which reflects the conclusion from the median joining network that not all Clade E lineages have been sampled. This unsampled lineage gave rise to the haplotypes on the English side of the English Channel, which were the ancestors to the Mediterranean and Bay of Biscay populations. Clade E has a wider geographic range than previously thought, and shows evidence of recent range expansion. Both Clade A and Clade E should be considered widespread species: Clade A globally and Clade E within Europe.

Quaternary climate change drives allo-peripatric speciation and refugial divergence in the Dysosma versipellis-pleiantha complex from different forest types in China

Subtropical China harbours the world's most diverse temperate flora, but little is known about the roles of geographical and eco-climatic factors underlying the region's exceptionally high levels of species diversity and endemism. Here we address this key question by investigating the spatio-temporal and ecological processes of divergence within the Dysosma versipellis-pleiantha species complex, endemic to subtropical China. Our cpDNA phylogeny showed that this monophyletic group of understory herbs is derived from a Late Pliocene ancestor of the Qinghai-Tibetan Plateau (QTP)/Southwest China. Genetic and ENM data in conjunction with niche differentiation analyses support that the early divergence of D. versipellis and D. pleiantha proceeded through allo-peripatric speciation, possibly triggered by Early Pleistocene climate change, while subsequent climate-induced cycles of range contractions/expansions enhanced the eco-geographical isolation of both taxa. Furthermore, modelling of population-genetic data indicated that major lineage divergences within D. versipellis likely resulted from long-term allopatric population isolation in multiple localized refugia over the last glacial/interglacial periods, and which in turn fostered endemic species formation (D. difformis, D. majoensis) from within D. versipellis in Southwest China. These findings point to an overriding role of Quaternary climate change in triggering essentially allopatric (incipient) speciation in this group of forest-restricted plant species in subtropical China.