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

Human mining activity across the ages determines the genetic structure of modern brown trout (Salmo trutta L.) populations

Humans have exploited the earth's metal resources for thousands of years leaving behind a legacy of toxic metal contamination and poor water quality. The southwest of England provides a well-defined example, with a rich history of metal mining dating to the Bronze Age. Mine water washout continues to negatively impact water quality across the region where brown trout (Salmo trutta L.) populations exist in both metal-impacted and relatively clean rivers. We used microsatellites to assess the genetic impact of mining practices on trout populations in this region. Our analyses demonstrated that metal-impacted trout populations have low genetic diversity and have experienced severe population declines. Metal-river trout populations are genetically distinct from clean-river populations, and also from one another, despite being geographically proximate. Using approximate Bayesian computation (ABC), we dated the origins of these genetic patterns to periods of intensive mining activity. The historical split of contemporary metal-impacted populations from clean-river fish dated to the Medieval period. Moreover, we observed two distinct genetic populations of trout within a single catchment and dated their divergence to the Industrial Revolution. Our investigation thus provides an evaluation of contemporary population genetics in showing how human-altered landscapes can change the genetic makeup of a species.

Evolutionary origin and demographic history of an ancient conifer (Juniperus microsperma) in the Qinghai-Tibetan Plateau

All Qinghai-Tibetan Plateau (QTP) endemic species are assumed to have originated recently, although very rare species most likely diverged early. These ancient species provide an excellent model to examine the origin and evolution of QTP endemic plants in response to the QTP uplifts and the climate changes that followed in this high altitude region. In this study, we examined these hypotheses by employing sequence variation from multiple nuclear and chloroplast DNA of 239 individuals of Juniperus microsperma and its five congeners. Both phylogenetic and population genetic analyses revealed that J. microsperma diverged from its sister clade comprising two species with long isolation around the Early Miocene, which corresponds to early QTP uplift. Demographic modeling and coalescent tests suggest that J. microsperma experienced an obvious bottleneck event during the Quaternary when the global climate greatly oscillated. The results presented here support the hypotheses that the QTP uplifts and Quaternary climate changes played important roles in shaping the evolutionary history of this rare juniper.

Widespread primary, but geographically restricted secondary, human introductions of wall lizards, Podarcis muralis

Establishing the introduction pathways of alien species is a fundamental task in invasion biology. The common wall lizard, Podarcis muralis, has been widely introduced outside of its native range in both Europe and North America, primarily through escaped pets or deliberate release of animals from captive or wild populations. Here, we use Bayesian clustering, approximate Bayesian computation (ABC) methods and network analyses to reconstruct the origin and colonization history of 23 non-native populations of wall lizards in England. Our analyses show that established populations in southern England originate from at least nine separate sources of animals from native populations in France and Italy. Secondary introductions from previously established non-native populations were supported for eleven (47%) populations. In contrast to the primary introductions, secondary introductions were highly restricted geographically and appear to have occurred within a limited time frame rather than being increasingly common. Together, these data suggest that extant wall lizard populations in England are the result of isolated accidental and deliberate releases of imported animals since the 1970s, with only local translocation of animals from established non-native populations. Given that populations introduced as recently as 25 years ago show evidence of having adapted to cool climate, discouraging further translocations may be important to prevent more extensive establishment on the south coast of England.

Historical Invasion Records Can Be Misleading: Genetic Evidence for Multiple Introductions of Invasive Raccoons (Procyon lotor) in Germany

Biological invasions provide excellent study systems to understand evolutionary, genetic and ecological processes during range expansions. There is strong evidence for positive effects of high propagule pressure and the associated higher genetic diversity on invasion success, but some species have become invasive despite small founder numbers. The raccoon (Procyon lotor) is often considered as a typical example for such a successful invasion resulting from a small number of founders. The species' largest non-native population in Germany is commonly assumed to stem from a small number of founders and two separate founding events in the 1930s and 1940s. In the present study we analyzed 407 raccoons at 20 microsatellite loci sampled from the invasive range in Western Europe to test if these assumptions are correct. Contrary to the expectations, different genetic clustering methods detected evidence for at least four independent introduction events that gave rise to genetically differentiated subpopulations. Further smaller clusters were either artifacts or resulted from founder events at the range margin and recent release of captive individuals. We also found genetic evidence for on-going introductions of individuals. Furthermore a novel randomization process was used to determine the potential range of founder population size that would suffice to capture all the alleles present in a cluster. Our results falsify the assumption that this species has become widespread and abundant despite being genetically depauperate and show that historical records of species introductions may be misleading.

Genome-enabled analysis of plant-pathogen migration

Trade in plant and plant products has profoundly affected the global distribution and diversity of plant pathogens. Identification of migration pathways can be used to monitor or manage pathogen movement for proactive disease management or quarantine measures. Genomics-based genetic marker discovery is allowing unprecedented collection of population genetic data for plant pathogens. These data can be used for detailed analysis of the ancestry of population samples and therefore for analysis of migration. Reconstruction of migration histories has confirmed previous hypotheses based on observational data and led to unexpected new findings on the origins of pathogens and source populations for past and recent migration. The choice of software for analysis depends on the type of migration being studied and the reproductive mode of the pathogen. Biased sampling and complex population structures are potential challenges to accurate inference of migration pathways.

Tigers of Sundarbans in India: is the population a separate conservation unit?

The Sundarbans tiger inhabits a unique mangrove habitat and are morphologically distinct from the recognized tiger subspecies in terms of skull morphometrics and body size. Thus, there is an urgent need to assess their ecological and genetic distinctiveness and determine if Sundarbans tigers should be defined and managed as separate conservation unit. We utilized nine microsatellites and 3 kb from four mitochondrial DNA (mtDNA) genes to estimate genetic variability, population structure, demographic parameters and visualize historic and contemporary connectivity among tiger populations from Sundarbans and mainland India. We also evaluated the traits that determine exchangeability or adaptive differences among tiger populations. Data from both markers suggest that Sundarbans tiger is not a separate tiger subspecies and should be regarded as Bengal tiger (P. t. tigris) subspecies. Maximum likelihood phylogenetic analyses of the mtDNA data revealed reciprocal monophyly. Genetic differentiation was found stronger for mtDNA than nuclear DNA. Microsatellite markers indicated low genetic variation in Sundarbans tigers (He= 0.58) as compared to other mainland populations, such as northern and Peninsular (Hebetween 0.67- 0.70). Molecular data supports migration between mainland and Sundarbans populations until very recent times. We attribute this reduction in gene flow to accelerated fragmentation and habitat alteration in the landscape over the past few centuries. Demographic analyses suggest that Sundarbans tigers have diverged recently from peninsular tiger population within last 2000 years. Sundarbans tigers are the most divergent group of Bengal tigers, and ecologically non-exchangeable with other tiger populations, and thus should be managed as a separate "evolutionarily significant unit" (ESU) following the adaptive evolutionary conservation (AEC) concept.

Software solutions for the livestock genomics SNP array revolution

Since the beginning of the genomic era, the number of available single nucleotide polymorphism (SNP) arrays has grown considerably. In the bovine species alone, 11 SNP chips not completely covered by intellectual property are currently available, and the number is growing. Genomic/genotype data are not standardized, and this hampers its exchange and integration. In addition, software used for the analyses of these data usually requires not standard (i.e. case specific) input files which, considering the large amount of data to be handled, require at least some programming skills in their production. In this work, we describe a software toolkit for SNP array data management, imputation, genome-wide association studies, population genetics and genomic selection. However, this toolkit does not solve the critical need for standardization of the genotypic data and software input files. It only highlights the chaotic situation each researcher has to face on a daily basis and gives some helpful advice on the currently available tools in order to navigate the SNP array data complexity.

Multiple host shifts by the emerging honeybee parasite, Varroa jacobsoni

Host shifts are a key mechanism of parasite evolution and responsible for the emergence of many economically important pathogens. Varroa destructor has been a major factor in global honeybee (Apis mellifera) declines since shifting hosts from the Asian honeybee (Apis cerana) > 50 years ago. Until recently, only two haplotypes of V. destructor (Korea and Japan) had successfully host shifted to A. mellifera. In 2008, the sister species V. jacobsoni was found for the first time parasitizing A. mellifera in Papua New Guinea (PNG). This recent host shift presents a serious threat to world apiculture but also provides the opportunity to examine host shifting in this system. We used 12 microsatellites to compare genetic variation of V. jacobsoni on A. mellifera in PNG with mites on A. cerana in both PNG and surrounding regions. We identified two distinct lineages of V. jacobsoni reproducing on A. mellifera in PNG. Our analysis indicated independent host shift events have occurred through small numbers of mites shifting from local A. cerana populations. Additional lineages were found in the neighbouring Papua and Solomon Islands that had partially host shifted to A. mellifera, that is producing immature offspring on drone brood only. These mites were likely in transition to full colonization of A. mellifera. Significant population structure between mites on the different hosts suggested host shifted V. jacobsoni populations may not still reproduce on A. cerana, although limited gene flow may exist. Our studies provide further insight into parasite host shift evolution and help characterize this new Varroa mite threat to A. mellifera worldwide.

Intervarietal and intravarietal genetic structure in Douglas-fir: nuclear SSRs bring novel insights into past population demographic processes, phylogeography, and intervarietal hybridization

Douglas-fir (Pseudotsuga menziesii) is one of numerous wide-range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty-eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice-covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.

Olive domestication and diversification in the Mediterranean Basin

Olive (Olea europaea ssp. europaea) is the most important oil fruit crop in temperate areas, but the origin of the cultivated olive remains unclear. The existence of one or several domestication events in the Mediterranean Basin (MB) is still debated. We analyzed a dataset of 387 cultivated and wild accessions that were genotyped at 25 simple-sequence repeat (SSR) loci. The sample represented genetic diversity at the geographic extremes of the MB. We inferred relationships among samples and also applied approximate Bayesian computation to estimate the most probable demographic model of our samples. Cultivated olives clustered into three different gene pools (Q1, Q2 and Q3), corresponding loosely to the west, central and eastern MB, respectively. Q1 consisted primarily of accessions from southern Spain, retained the fingerprint of a genetic bottleneck, and was closely related to accessions from the eastern MB. Q2 showed signs of recent admixture with wild olives and may derive from a local domestication event in the central MB. Overall our results suggest that admixture shaped olive germplasm and perhaps also local domestication events.

Post-fragmentation population structure in a cooperative breeding Afrotropical cloud forest bird: emergence of a source-sink population network

The impact of demographic parameters on the genetic population structure and viability of organisms is a long-standing issue in the study of fragmented populations. Demographic and genetic tools are now readily available to estimate census and effective population sizes and migration and gene flow rates with increasing precision. Here we analysed the demography and genetic population structure over a recent 15-year time span in five remnant populations of Cabanis's greenbul (Phyllastrephus cabanisi), a cooperative breeding bird in a severely fragmented cloud forest habitat. Contrary to our expectation, genetic admixture and effective population sizes slightly increased, rather than decreased between our two sampling periods. In spite of small effective population sizes in tiny forest remnants, none of the populations showed evidence of a recent population bottleneck. Approximate Bayesian modelling, however, suggested that differentiation of the populations coincided at least partially with an episode of habitat fragmentation. The ratio of meta-Ne to meta-Nc was relatively low for birds, which is expected for cooperative breeding species, while Ne /Nc ratios strongly varied among local populations. While the overall trend of increasing population sizes and genetic admixture may suggest that Cabanis's greenbuls increasingly cope with fragmentation, the time period over which these trends were documented is rather short relative to the average longevity of tropical species. Furthermore, the critically low Nc in the small forest remnants keep the species prone to demographic and environmental stochasticity, and it remains open if, and to what extent, its cooperative breeding behaviour helps to buffer such effects.

Frequent introgressions from diploid species contribute to the adaptation of the tetraploid Shepherd's purse (Capsella bursa-pastoris)

Adaptation is the most important ability for organisms to survive in diverse habitats. Animals have the option to escape from stressful environments, but plants do not. In plants, polyploids consist of about 30%-70% angiosperms and 95% ferns, of which some are important crops such as cotton and wheat. How polyploid plants adapt to various habitats has been a fundamental question remained largely unanswered. The tetraploid Shepherd's purse (Capsella bursa-pastoris) is one of the most successful plants on earth and has been distributed across the world, thus being an ideal model system for studying the adaptation of polyploids. We found that there are frequent introgressions from congeneric diploids to Shepherd's purse. Ecological niche modeling suggests that ecological differentiation is evident between the introgressed and non-introgressed C. bursa-pastoris, and the introgressions are a source of adaptation. This result links an evolutionary process to the adaptation of polyploids, and sheds light on the breeding strategy of polyploids as well. We conclude that frequent introgressions from congeneric diploids contributed to the acquisition of adequate genetic variations, thereby allowing C. bursa-pastoris to adapt to various habitats across the world. Our results highlight how a polyploid could have successfully established after it originated.

Assessing the phylogeographic history of the montane caddisfly Thremma gallicum using mitochondrial and restriction-site-associated DNA (RAD) markers

Repeated Quaternary glaciations have significantly shaped the present distribution and diversity of several European species in aquatic and terrestrial habitats. To study the phylogeography of freshwater invertebrates, patterns of intraspecific variation have been examined primarily using mitochondrial DNA markers that may yield results unrepresentative of the true species history. Here, population genetic parameters were inferred for a montane aquatic caddisfly, Thremma gallicum, by sequencing a 658-bp fragment of the mitochondrial CO1 gene, and 12,514 nuclear RAD loci. T. gallicum has a highly disjunct distribution in southern and central Europe, with known populations in the Cantabrian Mountains, Pyrenees, Massif Central, and Black Forest. Both datasets represented rangewide sampling of T. gallicum. For the CO1 dataset, this included 352 specimens from 26 populations, and for the RAD dataset, 17 specimens from eight populations. We tested 20 competing phylogeographic scenarios using approximate Bayesian computation (ABC) and estimated genetic diversity patterns. Support for phylogeographic scenarios and diversity estimates differed between datasets with the RAD data favouring a southern origin of extant populations and indicating the Cantabrian Mountains and Massif Central populations to represent highly diverse populations as compared with the Pyrenees and Black Forest populations. The CO1 data supported a vicariance scenario (north-south) and yielded inconsistent diversity estimates. Permutation tests suggest that a few hundred polymorphic RAD SNPs are necessary for reliable parameter estimates. Our results highlight the potential of RAD and ABC-based hypothesis testing to complement phylogeographic studies on non-model species.

Long-range gene flow and the effects of climatic and ecological factors on genetic structuring in a large, solitary carnivore: the Eurasian lynx

Due to their high mobility, large terrestrial predators are potentially capable of maintaining high connectivity, and therefore low genetic differentiation among populations. However, previous molecular studies have provided contradictory findings in relation to this. To elucidate patterns of genetic structure in large carnivores, we studied the genetic variability of the Eurasian lynx, Lynx lynx throughout north-eastern Europe using microsatellite, mitochondrial DNA control region and Y chromosome-linked markers. Using SAMOVA we found analogous patterns of genetic structure based on both mtDNA and microsatellites, which coincided with a relatively little evidence for male-biased dispersal. No polymorphism for the cytochrome b and ATP6 mtDNA genes and Y chromosome-linked markers were found. Lynx inhabiting a large area encompassing Finland, the Baltic countries and western Russia formed a single genetic unit, while some marginal populations were clearly divergent from others. The existence of a migration corridor was suggested to correspond with distribution of continuous forest cover. The lowest variability (in both markers) was found in lynx from Norway and Białowieża Primeval Forest (BPF), which coincided with a recent demographic bottleneck (Norway) or high habitat fragmentation (BPF). The Carpathian population, being monomorphic for the control region, showed relatively high microsatellite diversity, suggesting the effect of a past bottleneck (e.g. during Last Glacial Maximum) on its present genetic composition. Genetic structuring for the mtDNA control region was best explained by latitude and snow cover depth. Microsatellite structuring correlated with the lynx's main prey, especially the proportion of red deer (Cervus elaphus) in its diet. Eurasian lynx are capable of maintaining panmictic populations across eastern Europe unless they are severely limited by habitat continuity or a reduction in numbers. Different correlations of mtDNA and microsatellite population divergence patterns with climatic and ecological factors may suggest separate selective pressures acting on males and females in this solitary carnivore.

Congruent signals of population history but radically different patterns of genetic diversity between mitochondrial and nuclear markers in a mountain lizard

Historical factors, current population size, population connectivity and selective processes at linked loci contribute to shaping contemporary patterns of neutral genetic diversity. It is now widely acknowledged that nuclear and mitochondrial markers react differently to current demography as well as to past history, so the use of both types of markers is often advocated to gain insight on both historical and contemporary processes. We used 12 microsatellite loci genotyped in 13 populations of a mountain lizard (Iberolacerta bonnali) to test whether the historical scenario favoured by a previous mitochondrial study was also supported by nuclear markers and thereby evaluated the consequences of postglacial range movements on nuclear diversity. Congruent signals of recent history were revealed by nuclear and mitochondrial markers using an Approximate Bayesian computation approach, but contemporary patterns of mtDNA and nuclear DNA diversity were radically different. Although dispersal in this species is probably highly restricted at all spatial scales, colonization abilities have been historically good, suggesting capability for reestablishment of locally extinct populations except in fully disconnected habitats.

Using Ancient DNA to Understand Evolutionary and Ecological Processes

Ancient DNA provides a unique means to record genetic change through time and directly observe evolutionary and ecological processes. Although mostly based on mitochondrial DNA, the increasing availability of genomic sequences is leading to unprecedented levels of resolution. Temporal studies of population genetics have revealed dynamic patterns of change in many large vertebrates, featuring localized extinctions, migrations, and population bottlenecks. The pronounced climate cycles of the Late Pleistocene have played a key role, reducing the taxonomic and genetic diversity of many taxa and shaping modern populations. Importantly, the complex series of events revealed by ancient DNA data is seldom reflected in current biogeographic patterns. DNA preserved in ancient sediments and coprolites has been used to characterize a range of paleoenvironments and reconstruct functional relationships in paleoecological systems. In the near future, genome-level surveys of ancient populations will play an increasingly important role in revealing, calibrating, and testing evolutionary processes.

Contemporary population structure and post-glacial genetic demography in a migratory marine species, the blacknose shark, Carcharhinus acronotus

Patterns of population structure and historical genetic demography of blacknose sharks in the western North Atlantic Ocean were assessed using variation in nuclear-encoded microsatellites and sequences of mitochondrial (mt)DNA. Significant heterogeneity and/or inferred barriers to gene flow, based on microsatellites and/or mtDNA, revealed the occurrence of five genetic populations localized to five geographic regions: the southeastern U.S Atlantic coast, the eastern Gulf of Mexico, the western Gulf of Mexico, Bay of Campeche in the southern Gulf of Mexico and the Bahamas. Pairwise estimates of genetic divergence between sharks in the Bahamas and those in all other localities were more than an order of magnitude higher than between pairwise comparisons involving the other localities. Demographic modelling indicated that sharks in all five regions diverged after the last glacial maximum and, except for the Bahamas, experienced post-glacial, population expansion. The patterns of genetic variation also suggest that the southern Gulf of Mexico may have served as a glacial refuge and source for the expansion. Results of the study demonstrate that barriers to gene flow and historical genetic demography contributed to contemporary patterns of population structure in a coastal migratory species living in an otherwise continuous marine habitat. The results also indicate that for many marine species, failure to properly characterize barriers in terms of levels of contemporary gene flow could in part be due to inferences based solely on equilibrium assumptions. This could lead to erroneous conclusions regarding levels of connectivity in species of conservation concern.

Evolutionary ecotoxicology of perfluoralkyl substances (PFASs) inferred from multigenerational exposure: a case study with Chironomus riparius (Diptera, Chironomidae)

A multigeneration toxicity test on Chironomus riparius was performed with the aim of investigating the evolutionary consequences of exposure to perfluoralkyl substances (perfluorooctane sulfonic acid, PFOS; perfluorooctanoic acid, PFOA; perfluorobutane sulfonate, PFBS). Six-hundred larvae were bred per treatment and per generation until emergence and egg deposition under a nominal concentration of 10μg/L of contaminants. Newborn larvae were used to start the next generation. Evolution of genetic variability was evaluated along a total of 10 consecutive generations based on 5 microsatellite loci. Analysis of life-history traits (survival, sex ratio and reproduction) was also carried out. Rapid genetic variability reduction was observed in all treatments, including controls, across generations due to the test conditions. Nevertheless, an increased mutation rate determined a stronger conservation of genetic variability in PFOS and, at minor extent, in PFBS exposed populations compared to controls. No significant effects were induced by exposure to PFOA. Direct mutagenicity or induced stress conditions may be at the base of increased mutation rate, indicating the potential risk of mutational load caused by exposure to PFOS and PFBS. The test provided the opportunity to evaluate the use of approximate Bayesian computation (ABC) and coalescent approaches in evolutionary ecotoxicology. A weak performance was evidenced for ABC, either in terms of bias or dispersion of effective population sizes and of estimates of mutation rate. On the contrary, coalescent simulations proved the sensitivity of traditional genetic endpoints (i.e. heterozygosity and number of alleles) to the alteration of mutation rate, but not to erosion of genetic effective size.

Genetic signs of multiple colonization events in Baltic ciscoes with radiation into sympatric spring- and autumn-spawners confined to early postglacial arrival

Presence of sympatric populations may reflect local diversification or secondary contact of already distinct forms. The Baltic cisco (Coregonus albula) normally spawns in late autumn, but in a few lakes in Northern Europe sympatric autumn and spring- or winter-spawners have been described. So far, the evolutionary relationships and taxonomic status of these main life history forms have remained largely unclear. With microsatellites and mtDNA sequences, we analyzed extant and extinct spring- and autumn-spawners from a total of 23 Swedish localities, including sympatric populations. Published sequences from Baltic ciscoes in Germany and Finland, and Coregonus sardinella from North America were also included together with novel mtDNA sequences from Siberian C. sardinella. A clear genetic structure within Sweden was found that included two population assemblages markedly differentiated at microsatellites and apparently fixed for mtDNA haplotypes from two distinct clades. All sympatric Swedish populations belonged to the same assemblage, suggesting parallel evolution of spring-spawning rather than secondary contact. The pattern observed further suggests that postglacial immigration to Northern Europe occurred from at least two different refugia. Previous results showing that mtDNA in Baltic cisco is paraphyletic with respect to North American C. sardinella were confirmed. However, the inclusion of Siberian C. sardinella revealed a more complicated pattern, as these novel haplotypes were found within one of the two main C. albula clades and were clearly distinct from those in North American C. sardinella. The evolutionary history of Northern Hemisphere ciscoes thus seems to be more complex than previously recognized.

Genetic evidence for multiple sources of the non-native fish Cichlasoma urophthalmus (Günther; Mayan Cichlids) in southern Florida

The number and diversity of source populations may influence the genetic diversity of newly introduced populations and affect the likelihood of their establishment and spread. We used the cytochrome b mitochondrial gene and nuclear microsatellite loci to identify the sources of a successful invader in southern Florida, USA, Cichlasoma urophthalmus (Mayan cichlid). Our cytochrome b data supported an introduction from Guatemala, while our microsatellite data suggested movement of Mayan Cichlids from the upper Yucatán Peninsula to Guatemala and introductions from Guatemala and Belize to Florida. The mismatch between mitochondrial and nuclear genomes suggests admixture of a female lineage from Guatemala, where all individuals were fixed for the mitochondrial haplotype found in the introduced population, and a more diverse but also relatively small number of individuals from Belize. The Florida cytochrome b haplotype appears to be absent from Belize (0 out of 136 fish screened from Belize had this haplotype). Genetic structure within the Florida population was minimal, indicating a panmictic population, while Mexican and Central American samples displayed more genetic subdivision. Individuals from the Upper Yucatán Peninsula and the Petén region of Guatemala were more genetically similar to each other than to fish from nearby sites and movement of Mayan Cichlids between these regions occurred thousands of generations ago, suggestive of pre-Columbian human transportation of Mayan Cichlids through this region. Mayan Cichlids present a rare example of cytonuclear disequilibrium and reduced genetic diversity in the introduced population that persists more than 30 years (at least 7-8 generations) after introduction. We suggest that hybridization occurred in ornamental fish farms in Florida and may contribute their establishment in the novel habitat. Hybridization prior to release may contribute to other successful invasions.

Museum DNA reveals the demographic history of the endangered Seychelles warbler

The importance of evolutionary conservation – how understanding evolutionary forces can help guide conservation decisions – is widely recognized. However, the historical demography of many endangered species is unknown, despite the fact that this can have important implications for contemporary ecological processes and for extinction risk. Here, we reconstruct the population history of the Seychelles warbler (Acrocephalus sechellensis) – an ecological model species. By the 1960s, this species was on the brink of extinction, but its previous history is unknown. We used DNA samples from contemporary and museum specimens spanning 140 years to reconstruct bottleneck history. We found a 25% reduction in genetic diversity between museum and contemporary populations, and strong genetic structure. Simulations indicate that the Seychelles warbler was bottlenecked from a large population, with an ancestral N_e of several thousands falling to <50 within the last century. Such a rapid decline, due to anthropogenic factors, has important implications for extinction risk in the Seychelles warbler, and our results will inform conservation practices. Reconstructing the population history of this species also allows us to better understand patterns of genetic diversity, inbreeding and promiscuity in the contemporary populations. Our approaches can be applied across species to test ecological hypotheses and inform conservation.

Can novel genetic analyses help to identify low-dispersal marine invasive species?

Genetic methods can be a powerful tool to resolve the native versus introduced status of populations whose taxonomy and biogeography are poorly understood. The genetic study of introduced species is presently dominated by analyses that identify signatures of recent colonization by means of summary statistics. Unfortunately, such approaches cannot be used in low-dispersal species, in which recently established populations originating from elsewhere in the species' native range also experience periods of low population size because they are founded by few individuals. We tested whether coalescent-based molecular analyses that provide detailed information about demographic history supported the hypothesis that a sea squirt whose distribution is centered on Tasmania was recently introduced to mainland Australia and New Zealand through human activities. Methods comparing trends in population size (Bayesian Skyline Plots and Approximate Bayesian Computation) were no more informative than summary statistics, likely because of recent intra-Tasmanian dispersal. However, IMa2 estimates of divergence between putatively native and introduced populations provided information at a temporal scale suitable to differentiate between recent (potentially anthropogenic) introductions and ancient divergence, and indicated that all three non-Tasmanian populations were founded during the period of European settlement. While this approach can be affected by inaccurate molecular dating, it has considerable (albeit largely unexplored) potential to corroborate nongenetic information in species with limited dispersal capabilities.

Back to BaySICS: a user-friendly program for Bayesian Statistical Inference from Coalescent Simulations

Inference of population demographic history has vastly improved in recent years due to a number of technological and theoretical advances including the use of ancient DNA. Approximate Bayesian computation (ABC) stands among the most promising methods due to its simple theoretical fundament and exceptional flexibility. However, limited availability of user-friendly programs that perform ABC analysis renders it difficult to implement, and hence programming skills are frequently required. In addition, there is limited availability of programs able to deal with heterochronous data. Here we present the software BaySICS: Bayesian Statistical Inference of Coalescent Simulations. BaySICS provides an integrated and user-friendly platform that performs ABC analyses by means of coalescent simulations from DNA sequence data. It estimates historical demographic population parameters and performs hypothesis testing by means of Bayes factors obtained from model comparisons. Although providing specific features that improve inference from datasets with heterochronous data, BaySICS also has several capabilities making it a suitable tool for analysing contemporary genetic datasets. Those capabilities include joint analysis of independent tables, a graphical interface and the implementation of Markov-chain Monte Carlo without likelihoods.

Demographic processes shaping genetic variation of the solitarious phase of the desert locust

Between plagues, the solitarious desert locust (Schistocerca gregaria) is generally thought to exist as small populations, which are particularly prone to extinction events in arid regions of Africa and Asia. Given the high genetic structuring observed in one geographical area (the Eritrean coast) by former authors, a metapopulation dynamics model involving repeated extinction and colonization events was favoured. In this study, we assessed the validity of a demographic scenario involving temporary populations of the solitarious phase of the desert locust by analysing large-scale population genetic data. We scored 24 microsatellites in 23 solitarious population samples collected over most of the species range during remission. We found very little genetic structuring and little evidence of declining genetic diversity. A Bayesian clustering method distinguished four genetically differentiated units. Three groups were largely consistent with three population samples which had undergone recent bottleneck events. Nevertheless, the last genetically homogeneous unit included all individuals from the remaining 18 population samples and did not show evidence of demographic disequilibrium. An approximate Bayesian computation treatment indicated a large population size for this main genetic group, moderately reduced between plague and remission but still containing tens of thousands of individuals. Our results diverge from the hypothesis of a classical metapopulation dynamics model. They instead support the scenario in which large populations persist in the solitarious phase of the desert locust.

Phylogeography and post-glacial recolonization in wolverines (Gulo gulo) from across their circumpolar distribution

Interglacial-glacial cycles of the Quaternary are widely recognized in shaping phylogeographic structure. Patterns from cold adapted species can be especially informative - in particular, uncovering additional glacial refugia, identifying likely recolonization patterns, and increasing our understanding of species' responses to climate change. We investigated phylogenetic structure of the wolverine, a wide-ranging cold adapted carnivore, using a 318 bp of the mitochondrial DNA control region for 983 wolverines (n=209 this study, n=774 from GenBank) from across their full Holarctic distribution. Bayesian phylogenetic tree reconstruction and the distribution of observed pairwise haplotype differences (mismatch distribution) provided evidence of a single rapid population expansion across the wolverine's Holarctic range. Even though molecular evidence corroborated a single refugium, significant subdivisions of population genetic structure (0.01< ΦST <0.99, P<0.05) were detected. Pairwise ΦST estimates separated Scandinavia from Russia and Mongolia, and identified five main divisions within North America - the Central Arctic, a western region, an eastern region consisting of Ontario and Quebec/Labrador, Manitoba, and California. These data are in contrast to the nearly panmictic structure observed in northwestern North America using nuclear microsatellites, but largely support the nuclear DNA separation of contemporary Manitoba and Ontario wolverines from northern populations. Historic samples (c. 1900) from the functionally extirpated eastern population of Quebec/Labrador displayed genetic similarities to contemporary Ontario wolverines. To understand these divergence patterns, four hypotheses were tested using Approximate Bayesian Computation (ABC). The most supported hypothesis was a single Beringia incursion during the last glacial maximum that established the northwestern population, followed by a west-to-east colonization during the Holocene. This pattern is suggestive of colonization occurring in accordance with glacial retreat, and supports expansion from a single refugium. These data are significant relative to current discussions on the conservation status of this species across its range.

History of the invasive African olive tree in Australia and Hawaii: evidence for sequential bottlenecks and hybridization with the Mediterranean olive

Humans have introduced plants and animals into new continents and islands with negative effects on local species. This has been the case of the olive that was introduced in Australia, New Zealand and Pacific islands where it became invasive. Two subspecies were introduced in Australia, and each successfully invaded a specific area: the African olive in New South Wales (NSW) and the Mediterranean olive in South Australia. Here, we examine their origins and spread and analyse a large sample of native and invasive accessions with chloroplast and nuclear microsatellites. African olive populations from the invaded range exhibit two South African chlorotypes hence supporting an introduction from South Africa, while populations from South Australia exhibit chlorotypes of Mediterranean cultivars. Congruently, nuclear markers support the occurrence of two lineages in Australia but demonstrate that admixture took place, attesting that they hybridized early after introduction. Furthermore, using an approximate Bayesian computation framework, we found strong support for the serial introduction of the African olive from South Africa to NSW and then from NSW to Hawaii. The taxon experienced successive bottlenecks that did not preclude invasion, meaning that rapid decisions need to be taken to avoid naturalization where it has not established a large population yet.