Introduction history of Drosophila subobscura in the New World: a microsatellite-based survey using ABC methods

Europe as a secondary distribution hub in the worldwide invasion of the potato cyst nematode Globodera rostochiensis

The potato cyst nematode Globodera rostochiensis originates from the Andean Mountain region in South America and has unintentionally been introduced to all inhabited continents. Several studies have examined the population genetic structure of this pest in various countries by using microsatellite markers. However, merging microsatellite data produced from different laboratories is challenging and can introduce uncertainty when interpreting the results. To overcome this challenge and to explore invasion routes of this pest, we have genotyped 22 G. rostochiensis populations from all continents. Within populations, the highest genetic diversity was observed in the South American populations, the European populations showed an intermediate level of genetic diversity and the remaining populations were the less diverse. This confirmed pre-existing knowledge such as a first introduction event from South America to Europe, but the less diverse populations could originate either from South America or from Europe. At the continental scale, STRUCTURE genetic clustering output indicated that North America and Asia have experienced at least two introduction events. Comparing different evolutionary scenarios, the Approximate Bayesian Computation analysis showed that Europe served as a secondary distribution centre for the invasion of G. rostochiensis into all other continents (North America, Africa, Asia and Oceania).

Effective population size of Culex quinquefasciatus under insecticide-based vector management and following Hurricane Harvey in Harris County, Texas

Introduction: Culex quinquefasciatus is a mosquito species of significant public health importance due to its ability to transmit multiple pathogens that can cause mosquito-borne diseases, such as West Nile fever and St. Louis encephalitis. In Harris County, Texas, Cx. quinquefasciatus is a common vector species and is subjected to insecticide-based management by the Harris County Public Health Department. However, insecticide resistance in mosquitoes has increased rapidly worldwide and raises concerns about maintaining the effectiveness of vector control approaches. This concern is highly relevant in Texas, with its humid subtropical climate along the Gulf Coast that provides suitable habitat for Cx. quinquefasciatus and other mosquito species that are known disease vectors. Therefore, there is an urgent and ongoing need to monitor the effectiveness of current vector control programs. Methods: In this study, we evaluated the impact of vector control approaches by estimating the effective population size of Cx. quinquefasciatus in Harris County. We applied Approximate Bayesian Computation to microsatellite data to estimate effective population size. We collected Cx. quinquefasciatus samples from two mosquito control operation areas; 415 and 802, during routine vector monitoring in 2016 and 2017. No county mosquito control operations were applied at area 415 in 2016 and 2017, whereas extensive adulticide spraying operations were in effect at area 802 during the summer of 2016. We collected data for eighteen microsatellite markers for 713 and 723 mosquitoes at eight timepoints from 2016 to 2017 in areas 415 and 802, respectively. We also investigated the impact of Hurricane Harvey's landfall in the Houston area in August of 2017 on Cx. quinquefasciatus population fluctuation. Results: We found that the bottleneck scenario was the most probable historical scenario describing the impact of the winter season at area 415 and area 802, with the highest posterior probability of 0.9167 and 0.4966, respectively. We also detected an expansion event following Hurricane Harvey at area 802, showing a 3.03-fold increase in 2017. Discussion: Although we did not detect significant effects of vector control interventions, we found considerable influences of the winter season and a major hurricane on the effective population size of Cx. quinquefasciatus. The fluctuations in effective population size in both areas showed a significant seasonal pattern. Additionally, the significant population expansion following Hurricane Harvey in 2017 supports the necessity for post-hurricane vector-control interventions.

Experimental evolution on heat tolerance and thermal performance curves under contrasting thermal selection in Drosophila subobscura

Ectotherms can respond to global warming via evolutionary change of their upper thermal limits (CT_max ). Thus, the estimation of CT_max and its evolutionary potential is crucial to determine their vulnerability to global warming. However, CT_max estimations depend on the thermal stress intensity, and it is not completely clear whether its evolutionary capacity can be affected. Here, we performed an artificial selection experiment to increase heat tolerance using fast- and slow-ramping selection protocols in Drosophila subobscura. We found that heat tolerance evolved in both selection protocols, exhibiting similar evolutionary change rates and realized heritabilities. Additionally, we estimated the thermal performance curves (TPC) to evaluate correlated responses to selection on heat tolerance. We detected that thermal optimum increased in fast-ramping selection lines, but with a cost at the thermal performance breadth. Conversely, we did not detect changes in the TPC for the slow-ramping selection lines, indicating that thermal stress intensity has important effects on the evolution of thermal physiology of ectotherms. These findings, together with previous studies in D. subobscura reporting interpopulation variability and significant heritabilities for heat tolerance, suggest that evolutionary change can contribute to insect persistence in thermally changing environments and adaptation to global warming conditions.

Distribution modelling of an introduced species: do adaptive genetic markers affect potential range?

Biological invasions have increased in the last few decades mostly due to anthropogenic causes such as globalization of trade. Because invaders sometimes cause large economic losses and ecological disturbances, estimating their origin and potential geographical ranges is useful. Drosophila subobscura is native to the Old World but was introduced in the New World in the late 1970s and spread widely. We incorporate information on adaptive genetic markers into ecological niche modelling and then estimate the most probable geographical source of colonizers; evaluate whether the genetic bottleneck experienced by founders affects their potential distribution; and finally test whether this species has spread to all its potential suitable habitats worldwide. We find the environmental space occupied by this species in its native and introduced distributions are notably the same, although the introduced niche has shifted slightly towards higher temperature and lower precipitation. The genetic bottleneck of founding individuals was a key factor limiting the spread of this introduced species. We also find that regions in the Mediterranean and north-central Portugal show the highest probability of being the origin of the colonizers. Using genetically informed environmental niche modelling can enhance our understanding of the initial colonization and spread of invasive species, and also elucidate potential areas of future expansions worldwide.

Large-scale genetic admixture suggests high dispersal in an insect pest, the apple fruit moth

Knowledge about population genetic structure and dispersal capabilities is important for the development of targeted management strategies for agricultural pest species. The apple fruit moth, Argyresthia conjugella (Lepidoptera, Yponomeutidae), is a pre-dispersal seed predator. Larvae feed on rowanberries (Sorbus aucuparia), and when rowanberry seed production is low (i.e., inter-masting), the moth switches from laying eggs in rowanberries to apples (Malus domestica), resulting in devastating losses in apple crops. Using genetic methods, we investigated if this small moth expresses any local genetic structure, or alternatively if gene flow may be high within the Scandinavian Peninsula (~850.000 km², 55^o - 69^o N). Genetic diversity was found to be high (n = 669, mean He = 0.71). For three out of ten tetranucleotide STRs, we detected heterozygote deficiency caused by null alleles, but tests showed little impact on the overall results. Genetic differentiation between the 28 sampling locations was very low (average FST = 0.016, P < 0.000). Surprisingly, we found that all individuals could be assigned to one of two non-geographic genetic clusters, and that a third, geographic cluster was found to be associated with 30% of the sampling locations, with weak but significant signals of isolation-by-distance. Conclusively, our findings suggest wind-aided dispersal and spatial synchrony of both sexes of the apple fruit moth over large areas and across very different climatic zones. We speculate that the species may recently have had two separate genetic origins caused by a genetic bottleneck after inter-masting, followed by rapid dispersal and homogenization of the gene pool across the landscape. We suggest further investigations of spatial genetic similarities and differences of the apple fruit moth at larger geographical scales, through life-stages, across inter-masting, and during attacks by the parasitoid wasp (Microgaster politus).

Diurnal and Seasonal Activity Patterns of Drosophilid Species (Diptera: Drosophilidae) Present in Blackberry Agroecosystems With a Focus on Spotted-Wing Drosophila

Drosophilid species with different life histories have been shown to exhibit similar behavioral patterns related to locating and utilizing resources such as hosts, mates, and food sources. Drosophila suzukii (Matsumura) is an invasive species that differs from other frugivorous drosophilids in that females lay eggs in ripe and ripening fruits instead of overripe or rotten fruits. We hypothesized that there may be diurnal and/or seasonal patterns associated with the movement of drosophilid species into and out of crop fields and their attraction to fermentation-odor-based monitoring traps, and that D. suzukii would conform to similar patterns. To test these hypotheses, we deployed passive, 2-headed Malaise traps between crop fields and wooded edges to simultaneously catch flies moving into and out of crop fields. We also deployed monitoring traps with a fermentation-based bait between crop fields and wooded edges and within crop rows. Traps were deployed weekly in June-August in 2014 and 2015 at two commercial blackberry farm in Cleveland County, NC, and were checked hourly for 24 h, except during darkness. Both D. suzukii and other drosophilid species moved between crop fields and wooded edges and were attracted to monitoring traps primarily during the morning and evening hours. Whereas other drosophilids were captured in traps throughout the season, few D. suzukii were caught in traps until early to mid-July in both years and increased as the season progressed. Understanding D. suzukii movement and activity patterns is essential for the development of effective management strategies.

Twin introductions by independent invader mussel lineages are both associated with recent admixture with a native congener in Australia

Introduced species can impose profound impacts on the evolution of receiving communities with which they interact. If native and introduced taxa remain reproductively semi-isolated, human-mediated secondary contact may promote genetic exchange across newly created hybrid zones, potentially impacting native genetic diversity and invasive species spread. Here, we investigate the contributions of recent divergence histories and ongoing (post-introduction) gene flow between the invasive marine mussel, Mytilus galloprovincialis, and a morphologically indistinguishable and taxonomically contentious native Australian taxon, Mytilus planulatus. Using transcriptome-wide markers, we demonstrate that two contemporary M. galloprovincialis introductions into south-eastern Australia originate from genetically divergent lineages from its native range in the Mediterranean Sea and Atlantic Europe, where both introductions have led to repeated instances of admixture between introduced and endemic populations. Through increased genome-wide resolution of species relationships, combined with demographic modelling, we validate that mussels sampled in Tasmania are representative of the endemic Australian taxon (M. planulatus), but share strong genetic affinities to M. galloprovincialis. Demographic inferences indicate late-Pleistocene divergence times and historical gene flow between the Tasmanian endemic lineage and northern M. galloprovincialis, suggesting that native and introduced taxa have experienced a period of historical isolation of at least 100,000 years. Our results demonstrate that many genomic loci and sufficient sampling of closely related lineages in both sympatric (e.g. Australian populations) and allopatric (e.g. northern hemisphere Mytilus taxa) ranges are necessary to accurately (a) interpret patterns of intraspecific differentiation and to (b) distinguish contemporary invasive introgression from signatures left by recent divergence histories in high dispersal marine species. More broadly, our study fills a significant gap in systematic knowledge of native Australian biodiversity and sheds light on the intrinsic challenges for invasive species research when native and introduced species boundaries are not well defined.

Model-Based Species Delimitation: Are Coalescent Species Reproductively Isolated?

A large and growing fraction of systematists define species as independently evolving lineages that may be recognized by analyzing the population genetic history of alleles sampled from individuals belonging to those species. This has motivated the development of increasingly sophisticated statistical models rooted in the multispecies coalescent process. Specifically, these models allow for simultaneous estimation of the number of species present in a sample of individuals and the phylogenetic history of those species using only DNA sequence data from independent loci. These methods hold extraordinary promise for increasing the efficiency of species discovery but require extensive validation to ensure that they are accurate and precise. Whether the species identified by these methods correspond to the species that would be recognized by alternative species recognition criteria (such as measurements of reproductive isolation) is currently an open question and a subject of vigorous debate. Here, we perform an empirical test of these methods by making use of a classic model system in the history of speciation research, flies of the genus Drosophila. Specifically, we use the uniquely comprehensive data on reproductive isolation that is available for this system, along with DNA sequence data, to ask whether Drosophila species inferred under the multispecies coalescent model correspond to those recognized by many decades of speciation research. We found that coalescent based and reproductive isolation-based methods of inferring species boundaries are concordant for 77% of the species pairs. We explore and discuss potential explanations for these discrepancies. We also found that the amount of prezygotic isolation between two species is a strong predictor of the posterior probability of species boundaries based on DNA sequence data, regardless of whether the species pairs are sympatrically or allopatrically distributed. [BPP; Drosophila speciation; genetic distance; multispecies coalescent.]

Coalescent models characterize sources and demographic history of recent round goby colonization of Great Lakes and inland waters

The establishment and spread of aquatic invasive species are ecologically and economically harmful and a source of conservation concern internationally. Processes of species invasion have traditionally been inferred from observational data of species presence/absence and relative abundance. However, genetic-based approaches can provide valuable sources of inference. Restriction site-associated DNA sequencing was used to identify and genotype single nucleotide polymorphism (SNP) loci for Round Gobies (Neogobius melanostomus) (N = 440) from 18 sampling locations in the Great Lakes and in three Michigan, USA, drainages (Flint, Au Sable, and Cheboygan River basins). Sampled rivers differed in size, accessibility, and physical characteristics including man-made dispersal barriers. Population levels of genetic diversity and interpopulation variance in SNP allele frequency were used in coalescence-based approximate Bayesian computation (ABC) to statistically compare models representing competing hypotheses regarding source population, postcolonization dispersal, and demographic history in the Great Lakes and inland waters. Results indicate different patterns of colonization across the three drainages. In the Flint River, models indicate a strong population bottleneck (<3% of contemporary effective population size) and a single founding event from Saginaw Bay led to the colonization of inland river segments. In the Au Sable River, analyses could not distinguish potential source populations, but supported models indicated multiple introductions from one source population. In the Cheboygan River, supported models indicated that colonization likely proceeded from east (Lake Huron source) to west among inland locales sampled in the system. Despite the recent occupancy of Great Lakes and inland habitats, large numbers of loci analyzed in an ABC framework enable statistically supported identification of source populations and reconstruction of the direction of inland spread and demographic history following establishment. Information from analyses can direct management actions to limit the spread of invasive species from identified sources and most probable vectors into additional inland aquatic habitats.

Effects of photoperiod on life-history and thermal stress resistance traits across populations of Drosophila subobscura

INTRODUCTION

Organisms use environmental cues to match their phenotype with the future availability of resources and environmental conditions. Changes in the magnitude and frequency of environmental cues such as photoperiod and temperature along latitudes can be used by organisms to predict seasonal changes. While the role of temperature variation on the induction of plastic and seasonal responses is well established, the importance of photoperiod for predicting seasonal changes is less explored.

MATERIALS AND METHODS

Here we studied changes in life-history and thermal stress resistance traits in Drosophila subobscura in response to variation in photoperiod (6:18, 12:12 and 18:6 light:dark cycles) mimicking seasonal variations in day length. The populations of D. subobscura were collected from five locations along a latitudinal gradient (from North Africa and Europe). These populations were exposed to different photoperiods for two generations, whereafter egg-to-adult viability, productivity, dry body weight, thermal tolerance, and starvation resistance were assessed.

RESULTS

We found strong effects of photoperiod, origin of populations, and their interactions on life-history and stress resistance traits. Thermal resistance varied between the populations and the effect of photoperiod depended on the trait and the method applied for the assessment of thermal resistance.

PERSPECTIVES

Our results show a strong effect of the origin of population and photoperiod on a range of fitness-related traits and provide evidence for local adaptation to environmental cues (photoperiod by population interaction). The findings emphasize an important and often neglected role of photoperiod in studies on thermal resistance and suggest that cues induced by photoperiod may provide some buffer enabling populations to cope with a more variable and unpredictable future climate.

Patterns of geographic variation of thermal adapted candidate genes in Drosophila subobscura sex chromosome arrangements

BACKGROUND

The role of chromosomal arrangements in adaptation is supported by the repeatable clinal variation in inversion frequencies across continents in colonizing species such as Drosophila subobscura. However, there is a lack of knowledge on the genetic variation in genes within inversions, possibly targets of climatic selection, across a geographic latitudinal gradient. In the present study we analysed four candidate loci for thermal adaptation, located close to the breakpoints, in two chromosomal arrangements of the sex (A) chromosome of Drosophila subobscura with different thermal preferences. Individual chromosomes with A₂ (the inverted arrangement considered warm adapted) or A_ST (the standard ancestral arrangement considered cold adapted) were sequenced across four European localities at varying latitudes, up to ~ 2500 Kms apart.

RESULTS

Importantly, we found very low differentiation for each specific arrangement across populations as well as no clinal patterns of genomic variation. This suggests wide gene exchange along the cline. Differentiation between the sex chromosome arrangements was significant in the two more proximal regions relative to the A_ST orientation but not in the distal ones, independently of their location inside or outside the inversion. This can be possibly due to variation in the levels of gene flux and/or selection acting in these regions.

CONCLUSIONS

Gene flow appears to have homogenized the genetic content within-arrangement at a wide geographical scale, despite the expected diverse selective pressures in the specific natural environments of the different populations sampled. It is thus likely that the inversion frequency clines in this species are being maintained by local adaptation in face of gene flow. The differences between arrangements at non-coding regions might be associated with the previously observed differential gene expression in different thermal regimes. Higher resolution genomic scans for individual chromosomal arrangements performed over a large environmental gradient are needed to find the targets of selection and further elucidate the adaptive mechanisms maintaining chromosomal inversion polymorphisms.

Genetic analyses favour an ancient and natural origin of elephants on Borneo

The origin of the elephant on the island of Borneo remains elusive. Research has suggested two alternative hypotheses: the Bornean elephant stems either from a recent introduction in the 17th century or from an ancient colonization several hundreds of thousands years ago. Lack of elephant fossils has been interpreted as evidence for a very recent introduction, whereas mtDNA divergence from other Asian elephants has been argued to favor an ancient colonization. We investigated the demographic history of Bornean elephants using full-likelihood and approximate Bayesian computation analyses. Our results are at odds with both the recent and ancient colonization hypotheses, and favour a third intermediate scenario. We find that genetic data favour a scenario in which Bornean elephants experienced a bottleneck during the last glacial period, possibly as a consequence of the colonization of Borneo, and from which it has slowly recovered since. Altogether the data support a natural colonization of Bornean elephants at a time when large terrestrial mammals could colonise from the Sunda shelf when sea levels were much lower. Our results are important not only in understanding the unique history of the colonization of Borneo by elephants, but also for their long-term conservation.

Biases of STRUCTURE software when exploring introduction routes of invasive species

Population genetic methods are widely used to retrace the introduction routes of invasive species. The unsupervised Bayesian clustering algorithm implemented in STRUCTURE is amongst the most frequently used of these methods, but its ability to provide reliable information about introduction routes has never been assessed. We simulated microsatellite datasets to evaluate the extent to which the results provided by STRUCTURE were misleading for the inference of introduction routes. We focused on an invasion scenario involving one native and two independently introduced populations, because it is the sole scenario that can be rejected when obtaining a particular clustering with a STRUCTURE analysis at K = 2 (two clusters). Results were classified as "misleading" or "non-misleading". We investigated the influence of effective size, bottleneck severity and number of loci on the type and frequency of misleading results. We showed that misleading STRUCTURE results were obtained for 10% of all simulated datasets. Our results highlighted two categories of misleading output. The first occurs when the native population has a low level of diversity. In this case, the two introduced populations may be very similar, despite their independent introduction histories. The second category results from convergence issues in STRUCTURE for K = 2, with strong bottleneck severity and/or large numbers of loci resulting in high levels of differentiation between the three populations. Overall, the risk of being misled by STRUCTURE in the context of introduction routes inferences is moderate, but it is important to remain cautious when low genetic diversity or genuine multimodality between runs are involved.

Loss of dendritic connectivity in southern California's urban riverscape facilitates decline of an endemic freshwater fish

Life history adaptations and spatial configuration of metapopulation networks allow certain species to persist in extreme fluctuating environments, yet long-term stability within these systems relies on the maintenance of linkage habitat. Degradation of such linkages in urban riverscapes can disrupt this dynamic in aquatic species, leading to increased extinction debt in local populations experiencing environment-related demographic flux. We used microsatellites and mtDNA to examine the effects of collapsed network structure in the endemic Santa Ana sucker Catostomus santaanae of southern California, a threatened species affected by natural flood-drought cycles, "boom-and-bust" demography, hybridization and presumed artificial transplantation. Our results show a predominance of drift-mediated processes in shaping population structure and that reverse mechanisms for counterbalancing the genetic effects of these phenomena have dissipated with the collapse of dendritic connectivity. We use approximate Bayesian models to support two cases of artificial transplantation and provide evidence that one of the invaded systems better represents the historic processes that maintained genetic variation within watersheds than any remaining drainages where C. santaanae is considered native. We further show that a stable dry gap in the northern range is preventing genetic dilution of pure C. santaanae persisting upstream of a hybrid assemblage involving a non-native sucker and that local accumulation of genetic variation in the same drainage is influenced by position within the network. This work has important implications for declining species that have historically relied on dendritic metapopulation networks to maintain source-sink dynamics in phasic environments, but no longer possess this capacity in urban-converted landscapes.

Inversion evolutionary rates might limit the experimental identification of inversion breakpoints in non-model species

Chromosomal inversions are structural changes that alter gene order but generally not gene content in the affected region. In Drosophila, extensive cytological studies revealed the widespread character of inversion polymorphism, with evidence for its adaptive character. In Drosophila subobscura, polymorphism affects both its four large autosomal elements and its X (A) chromosome. The characterization of eight of these autosomal inversions breakpoints revealed that most of them originated through the staggered-breaks mechanism. Here, we have performed chromosomal walks to identify the breakpoints of two X-chromosome widely distributed inversions -A₂ and A₁- of D. subobscura. Inversion A₂ is considered a warm-adapted arrangement that exhibits parallel latitudinal clines in the species ancestral distribution area and in both American subcontinents, whereas inversion A₁ is only present in the Palearctic region where it presents an east-west cline. The duplication detected at the A₂ inversion breakpoints is consistent with its origin by the staggered-breaks mechanism. Inversion A₁ breakpoints could not be molecularly identified even though they could be narrowly delimited. This result points to chromosome walking limitations when using as a guide the genome of other species. Limitations stem from the rate of evolution by paracentric inversions, which in Drosophila is highest for the X chromosome.

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 assessment of the effects of streamscape succession on coho salmon Oncorhynchus kisutch colonization in recently deglaciated streams

Measures of genetic diversity within and among populations and historical geomorphological data on stream landscapes were used in model simulations based on approximate Bayesian computation (ABC) to examine hypotheses of the relative importance of stream features (geomorphology and age) associated with colonization events and gene flow for coho salmon Oncorhynchus kisutch breeding in recently deglaciated streams (50-240 years b.p.) in Glacier Bay National Park (GBNP), Alaska. Population estimates of genetic diversity including heterozygosity and allelic richness declined significantly and monotonically from the oldest and largest to youngest and smallest GBNP streams. Interpopulation variance in allele frequency increased with increasing distance between streams (r = 0·435, P < 0·01) and was inversely related to stream age (r = -0·281, P < 0·01). The most supported model of colonization involved ongoing or recent (<10 generations before sampling) colonization originating from large populations outside Glacier Bay proper into all other GBNP streams sampled. Results here show that sustained gene flow from large source populations is important to recently established O. kisutch metapopulations. Studies that document how genetic and demographic characteristics of newly founded populations vary associated with successional changes in stream habitat are of particular importance to and have significant implications for, restoration of declining or repatriation of extirpated populations in other regions of the species' native range.

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.

Gene expression clines reveal local adaptation and associated trade-offs at a continental scale

Local adaptation, where fitness in one environment comes at a cost in another, should lead to spatial variation in trade-offs between life history traits and may be critical for population persistence. Recent studies have sought genomic signals of local adaptation, but often have been limited to laboratory populations representing two environmentally different locations of a species’ distribution. We measured gene expression, as a proxy for fitness, in males of Drosophila subobscura, occupying a 20° latitudinal and 11 °C thermal range. Uniquely, we sampled six populations and studied both common garden and semi-natural responses to identify signals of local adaptation. We found contrasting patterns of investment: transcripts with expression positively correlated to latitude were enriched for metabolic processes, expressed across all tissues whereas negatively correlated transcripts were enriched for reproductive processes, expressed primarily in testes. When using only the end populations, to compare our results to previous studies, we found that locally adaptive patterns were obscured. While phenotypic trade-offs between metabolic and reproductive functions across widespread species are well-known, our results identify underlying genetic and tissue responses at a continental scale that may be responsible for this. This may contribute to understanding population persistence under environmental change.

Range expansion of the Bluetongue vector, Culicoides imicola, in continental France likely due to rare wind-transport events

The role of the northward expansion of Culicoides imicola Kieffer in recent and unprecedented outbreaks of Culicoides-borne arboviruses in southern Europe has been a significant point of contention. We combined entomological surveys, movement simulations of air-borne particles, and population genetics to reconstruct the chain of events that led to a newly colonized French area nestled at the northern foot of the Pyrenees. Simulating the movement of air-borne particles evidenced frequent wind-transport events allowing, within at most 36 hours, the immigration of midges from north-eastern Spain and Balearic Islands, and, as rare events, their immigration from Corsica. Completing the puzzle, population genetic analyses discriminated Corsica as the origin of the new population and identified two successive colonization events within west-Mediterranean basin. Our findings are of considerable importance when trying to understand the invasion of new territories by expanding species.

Recent similarity in distribution ranges does not mean a similar postglacial history: a phylogeographical study of the boreal tree species Alnus incana based on microsatellite and chloroplast DNA variation

We reconstructed the historical pattern of postglacial biogeographic range expansion of the boreal tree species Alnus incana in Europe. To assess population genetic structure and diversity, we performed a combined analysis of nuclear microsatellite loci and chloroplast DNA sequences (65 populations, 1004 individuals). Analysis of haplotype and microsatellite diversity revealed that southeastern refugial populations situated in the Carpathians and the Balkan Peninsula did not spread north and cannot be considered as important source populations for postglacial recolonization of Europe; populations in Eastern Europe did not establish Fennoscandian populations; populations in Fennoscandia and Eastern Europe have no unique genetic cluster, but represent a mix with a predominant cluster typical for Central Europe; and that colonization of Fennoscandia and Eastern Europe took place from Central Europe. Our findings highlight the importance of an effective refugium in Central Europe located outside classical southern refugia confirming the existence of northern refugia for boreal trees in Europe. The postglacial range expansion of A. incana did not follow the model established for Picea abies. Fennoscandian populations are not derived from Eastern European ones, but from Central European ones.

Introduced Drosophila subobscura populations perform better than native populations during an oviposition choice task due to increased fecundity but similar learning ability

The success of invasive species is tightly linked to their fitness in a putatively novel environment. While quantitative components of fitness have been studied extensively in the context of invasive species, fewer studies have looked at qualitative components of fitness, such as behavioral plasticity, and their interaction with quantitative components, despite intuitive benefits over the course of an invasion. In particular, learning is a form of behavioral plasticity that makes it possible to finely tune behavior according to environmental conditions. Learning can be crucial for survival and reproduction of introduced organisms in novel areas, for example, for detecting new predators, or finding mates or oviposition sites. Here we explored how oviposition performance evolved in relation to both fecundity and learning during an invasion, using native and introduced Drosophila subobscura populations performing an ecologically relevant task. Our results indicated that, under comparable conditions, invasive populations performed better during our oviposition task than did native populations. This was because invasive populations had higher fecundity, together with similar cognitive performance when compared to native populations, and that there was no interaction between learning and fecundity. Unexpectedly, our study did not reveal an allocation trade‐off (i.e., a negative relationship) between learning and fecundity. On the contrary, the pattern we observed was more consistent with an acquisition trade‐off, meaning that fecundity could be limited by availability of resources, unlike cognitive ability. This pattern might be the consequence of escaping natural enemies and/or competitors during the introduction. The apparent lack of evolution of learning may indicate that the introduced population did not face novel cognitive challenges in the new environment (i.e., cognitive “pre‐adaptation”). Alternatively, the evolution of learning may have been transient and therefore not detected.

The Wide Potential Trophic Niche of the Asiatic Fruit Fly Drosophila suzukii: The Key of Its Invasion Success in Temperate Europe?

The Asiatic fruit fly Drosophila suzukii has recently invaded Europe and North and South America, causing severe damage to fruit production systems. Although agronomic host plants of that fly are now well documented, little is known about the suitability of wild and ornamental hosts in its exotic area. In order to study the potential trophic niche of D. suzukii with relation to fruit characteristics, fleshy fruits from 67 plant species were sampled in natural and anthropic ecosystems (forests, hedgerows, grasslands, coastal areas, gardens and urban areas) of the north of France and submitted to experimental infestations. A set of fruit traits (structure, colour, shape, skin texture, diameter and weight, phenology) potentially interacting with oviposition choices and development success of D. suzukii was measured. Almost half of the tested plant species belonging to 17 plant families allowed the full development of D. suzukii. This suggests that the extreme polyphagy of the fly and the very large reservoir of hosts producing fruits all year round ensure temporal continuity in resource availability and contribute to the persistence and the exceptional invasion success of D. suzukii in natural habitats and neighbouring cultivated systems. Nevertheless, this very plastic trophic niche is not systematically beneficial to the fly. Some of the tested plants attractive to D. suzukii gravid females stimulate oviposition but do not allow full larval development. Planted near sensitive crops, these "trap plants" may attract and lure D. suzukii, therefore contributing to the control of the invasive fly.

Wing trait-inversion associations in Drosophila subobscura can be generalized within continents, but may change through time

Clinal variation is one of the most emblematic examples of the action of natural selection at a wide geographical range. In Drosophila subobscura, parallel clines in body size and inversions, but not in wing shape, were found in Europe and South and North America. Previous work has shown that a bottleneck effect might be largely responsible for differences in wing trait-inversion association between one European and one South American population. One question still unaddressed is whether the associations found before are present across other populations of the European and South American clines. Another open question is whether evolutionary dynamics in a new environment can lead to relevant changes in wing traits-inversion association. To analyse geographical variation in these associations, we characterized three recently laboratory founded D. subobscura populations from both the European and South American latitudinal clines. To address temporal variation, we also characterized the association at a later generation in the European populations. We found that wing size and shape associations can be generalized across populations of the same continent, but may change through time for wing size. The observed temporal changes are probably due to changes in the genetic content of inversions, derived from adaptation to the new, laboratory environment. Finally, we show that it is not possible to predict clinal variation from intrapopulation associations. All in all this suggests that, at least in the present, wing traits-inversion associations are not responsible for the maintenance of the latitudinal clines in wing shape and size.

The Temporal and Spatial Invasion Genetics of the Western Corn Rootworm (Coleoptera: Chrysomelidae) in Southern Europe

This study describes the genetics of the western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte in southern Europe during the introduction (1996–2001) and establishment/spread (2002–2011) phases of its invasion. The Diabrotica microsatellite core-set was used to perform traditional population genetics analyses. Our results indicated that during the introduction phase genetic diversity and population genetic structure were lower overall as compared to the establishment/spread phase. Unusually high genetic differentiation was found between the Italy and southern Europe comparisons, including high differentiation between Italian populations separated by a short distance during the establishment/spread phase. STRUCTURE analysis revealed two genetic clusters during the introduction phase and two genetic clusters during the establishment/spread phase. However, bottlenecked populations were only detected during the invasion phase. A small but significant isolation by distance effect was noted in both phases. Serbia was the geographic source of WCR to Croatia and Hungary in the introduction phase, while the United States of America was the possible source of WCR to Italy in 2001. These introductory populations were the subsequent source of individuals sampled during the establishment/spread phase. Repeated introductions and admixture events in southern Europe may have resulted in genetically diverse WCR populations that have attained 83% of all known alleles worldwide.

Using ABC and microsatellite data to detect multiple introductions of invasive species from a single source

The introduction of invasive species to new locations (that is, biological invasions) can have major impact on biodiversity, agriculture and public health. As such, determining the routes and modality of introductions with genetic data has become a fundamental goal in molecular ecology. To assist with this goal, new statistical methods and frameworks have been developed, such as approximate Bayesian computation (ABC) for inferring invasion history. Here, we present a model of invasion accounting for multiple introductions from a single source (MISS), a heretofore largely unexplored model. We simulate microsatellite data to evaluate the power of ABC to distinguish between single and multiple introductions from the same source, under a range of demographic parameters. We also apply ABC to microsatellite data from three invasions of bumblebee in New Zealand. In addition, we assess the performance of several methods of summary statistics selection. Our simulated results suggested good ability to distinguish between one- and two-wave models over much but not all of the parameter space tested, independent of summary statistics used. Globally, parameter estimation was good except for bottleneck timing. For one of the bumblebee species, we clearly rejected the MISS model, while for the other two we found inconclusive results. Since a second wave may provide genetic reinforcement to initial colonists, help relieve inbreeding among founders, or increase the hazard of the invasion, its detection may be crucial for managing invasions; we suggest that the MISS model could be considered as a potential model in future theoretical and empirical studies of invasions.

The devil is in the details: genetic variation in introduced populations and its contributions to invasion

The influence of genetic variation on invasion success has captivated researchers since the start of the field of invasion genetics 50 years ago. We review the history of work on this question and conclude that genetic variation-as surveyed with molecular markers-appears to shape invasion rarely. Instead, there is a significant disconnect between marker assays and ecologically relevant genetic variation in introductions. We argue that the potential for adaptation to facilitate invasion will be shaped by the details of genotypes affecting phenotypes, and we highlight three areas in which we see opportunities to make powerful new insights. (i) The genetic architecture of adaptive variation. Traits shaped by large-effect alleles may be strongly impacted by founder events yet more likely to respond to selection when genetic drift is strong. Large-effect loci may be especially relevant for traits involved in biotic interactions. (ii) Cryptic genetic variation exposed during invasion. Introductions have strong potential to uncover masked variation due to alterations in genetic and ecological environments. (iii) Genetic interactions during admixture of multiple source populations. As divergence among sources increases, positive followed by increasingly negative effects of admixture should be expected. Although generally hypothesized to be beneficial during invasion, admixture is most often reported among sources of intermediate divergence, supporting the possibility that incompatibilities among divergent source populations might be limiting their introgression. Finally, we note that these details of invasion genetics can be coupled with comparative demographic analyses to link genetic changes to the evolution of invasiveness itself.

Inferring the demographic history of Drosophila subobscura from nucleotide variation at regions not affected by chromosomal inversions

Drosophila subobscura presents a rich and complex chromosomal inversion polymorphism. It can thus be considered a model system (i) to study the mechanisms originating inversions and how inversions affect the levels and patterns of variation in the inverted regions and (ii) to study adaptation at both the single-gene and chromosomal inversion levels. It is therefore important to infer its demographic history as previous information indicated that its nucleotide variation is not at mutation-drift equilibrium. For that purpose, we sequenced 16 noncoding regions distributed across those parts of the J chromosome not affected by inversions in the studied population and possibly either by other selective events. The pattern of variation detected in these 16 regions is similar to that previously reported within different chromosomal arrangements, suggesting that the latter results would, thus, mainly reflect recent demographic events rather than the partial selective sweep imposed by the origin and frequency increase of inversions. Among the simple demographic models considered in our Approximate Bayesian Computation analysis of variation at the 16 regions, the model best supported by the data implies a population size expansion soon after the penultimate glacial period. This model constitutes a better null model, and it is therefore an important resource for subsequent studies aiming among others to uncover selective events across the species genome. Our results also highlight the importance of introducing the possibility of multiple hits in the coalescent simulations with an outgroup.

Clinal patterns of desiccation and starvation resistance in ancestral and invading populations of Drosophila subobscura

As invading species expand, they eventually encounter physical and biotic stressors that limit their spread. We examine latitudinal and climatic variation in physiological tolerance in one native and two invading populations of Drosophila subobscura. These flies are native to the Palearctic region, but invaded both South and North America around 1980 and spread rapidly across 15° of latitude on each continent. Invading flies rapidly evolved latitudinal clines in chromosome inversion frequencies and in wing size that parallel those of native populations in the Old World. Here we investigate whether flies on all three continents have evolved parallel clines in desiccation and starvation tolerance, such that flies in low-latitude regions (hot, dry) might have increased stress resistance. Starvation tolerance does not vary with latitude or climate on any continent. In contrast, desiccation tolerance varies clinally with latitude on all three continents, although not in parallel. In North American and Europe, desiccation tolerance is inversely related to latitude, as expected. But in South America, desiccation tolerance increases with latitude and is greatest in relatively cool and wet areas. Differences among continents in latitudinal patterns of interspecific-competition potentially influence clinal selection for physiological resistance, but no simple pattern is evident on these continents.

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.

Drosophila suzukii: the genetic footprint of a recent, worldwide invasion

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.