Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gene Flow and Maintenance of Genetic Diversity in Invasive Mosquitofish (Gambusia holbrooki)

Abstract

Genetic analyses contribute to studies of biological invasions by mapping the origin and dispersal patterns of invasive species occupying new territories. Using microsatellite loci, we assessed the genetic diversity and spatial population structure of mosquitofish (Gambusia holbrooki) that had invaded Spanish watersheds, along with the American locations close to the suspected potential source populations. Mosquitofish populations from the Spanish streams that were studied had similar levels of genetic diversity to the American samples; therefore, these populations did not appear to have undergone substantial losses of genetic diversity during the invasion process. Population structure analyses indicated that the Spanish populations fell into four main clusters, which were primarily associated with hydrography. Dispersal patterns indicated that local populations were highly connected upstream and downstream through active dispersal, with an average of 21.5% fish from other locations in each population. After initially introducing fish to one location in a given basin, such dispersal potential might contribute to the spread and colonization of suitable habitats throughout the entire river basin. The two-dimension isolation-by-distance pattern here obtained, indicated that the human-mediated translocation of mosquitofish among the three study basins is a regular occurrence. Overall, both phenomena, high natural dispersal and human translocation, favor gene flow among river basins and the retention of high genetic diversity, which might help retain the invasive potential of mosquitofish populations.

Introduction

Biological invasions are a central component of global change, and a major threat to the biodiversity of freshwater ecosystems [1]. Invaders may change radically the functioning of an ecosystem, to which they are introduced, resulting in the decline or extinction of native species through predation, competition, and habitat alteration [2]. Often, biological invasions begin when humans introduce a few individuals of a species to a new environment. Once established, the new population spreads to neighboring locations by natural dispersal. For example, the Eurasian spread of the topmouth gudgeon (Pseudorasbora parva) began accidentally when humans introduced it during the establishment of new cultured stocks of the common carp (Cyprinus carpio), and was followed by further natural dispersal of short distances [3]. Genetic variation is closely linked to the success of biological invasions [4]. When introductions begin with just a few individuals, reduced genetic diversity is expected during the first stages of the invasion [5], [6]. However, several studies have shown high diversity in populations at later stages of the invasion process (see [7]). The recovery of genetic diversity in invaded territories might result from gene flow between recently established populations within the invaded range that have become increasingly interconnected. Multiple introductions from divergent stocks also contribute towards increasing local diversity in invaded territories [8], particularly if these introductions occur separately in time from the initial founder event (e.g., [9]).

The eastern mosquitofish, Gambusia holbrooki, is one of the most commonly introduced freshwater species [10]. At present, established populations of this species in Europe, Africa, Asia, and Australia are causing local extinction and decline of several native fish and amphibian species [11], [12]. For example, the aggressive behavior of G. holbrooki has caused a decline in the feeding rates and reproductive success of two Iberian endemic fish species, Valencia hispanica and Aphanius iberus [13]. Similarly, Carmona-Catot et al. [14] showed that introduced G. holbrooki were able to competitively displace A. iberus populations. Mosquitofish introductions were originally supported by governmental health agencies to control mosquito populations, which are vectors of various diseases, such as malaria [15]. In Europe, 12 individuals of G. holbrooki were initially introduced into a pond in southern Spain in 1921 [16]. Subsequently, humans spread G. holbrooki throughout the Mediterranean basin [17]. Despite their small size, mosquitofish are extremely successful in new environments [15]. Both the invertivorous diet and wide ecological tolerance of G. holbrooki have probably contributed to its successful integration into Iberian fish communities [18].

Temporal fluctuations in population size reduce the average effective population size (Ne), with reduced effective sizes intensifying the loss of diversity due to genetic drift [19]. Although there is a major decline in the size of mosquitofish populations during winter after the summer flush [20], several studies have indicated high genetic diversity within American populations that often exceeds average values described for freshwater fishes [21]. The high reproductive potential generated by overwintering pregnant females, multiple paternity, and offspring reaching maturity within a few weeks probably contribute towards maintaining large effective population sizes and preventing the loss of population diversity [22][24]. Moreover, gene flow between seasonally isolated demes favors population diversity in large territories [25]. For instance, sporadic individual exchange among close populations prevents divergence among collections within basins in invaded territories [26].

Several models, such as isolated populations with no current migration and metapopulations of ephemeral populations connected by gene flow, may explain the population structure of organisms in linear river systems (e.g., [27]). Native mosquitofish populations usually represent single breeding units [28], while large transects within a river basin are occupied by a single population with ephemeral local subpopulations [29]. Source-sink dynamics are sometimes responsible for the population structure of mosquitofish (e.g., [25], [30]). Along a river system, dominant downstream gene flow increases the genetic diversity of lowland populations [31], [32].

DNA molecular markers contribute to improving our understanding of evolutionary changes that occur during biological invasions [33][36]. Highly polymorphic microsatellite loci provide the discrimination required to address questions about population structure and gene flow [24], [27]. In this study, we used microsatellite loci to evaluate putative losses of genetic variation during introductions of G. holbrooki. We aimed to understand the mechanisms that contribute towards retaining levels of diversity within populations inhabiting invaded rivers compared to populations in native basins.

Materials and Methods

Ethics statement

Animal samples were collected and manipulated under a permit (SF/012/2011) provided by the Agriculture, Fisheries, Food and Environment Department of the Autonomous Community of Catalonia. All work was performed in compliance with and approved by the Ethics Committee of the University of Girona and meets the requirements stated by the Spanish (RD53/2013) and Catalonian (D214/1997) laws of animal care, and experimentation.

Sample collection

A total of 556 G. holbrooki were collected from 15 sites along three watersheds (Muga, Fluvià, and Ter rivers) in northeastern Spain. The largest of these rivers is the Ter, with a basin area of 2955 km2, with its headwaters in the Pyrenees and its upper course being partially snow-fed. The Fluvià (974 km2) and Muga (758 km2) are typical Mediterranean streams with smaller watersheds, and have their headwaters located in mountainous areas. All three rivers are subject to a Mediterranean climate, with severe summer droughts and autumn floods [37]. The Ter and Muga rivers have many small weirs, along with a few large dams that form major barriers, altering connectivity among fish populations, whereas the Fluvià only has weirs. Mosquitofish are currently absent from the upper course of these watersheds; hence, we collected samples from the middle and lower courses of these three watersheds (Fig. 1, Table 1). Sampling sites were shallow areas (<1.5 m depth) along the riverbank, with low water velocity and dense vegetation, usually reed beds (Phragmites australis). We also analyzed 36 individuals from the Potomac River (Washington) and 16 individuals from Brunswick (North Carolina), because it has been suggested that these populations are the closest to the main American source of G. holbrooki individuals that were introduced to Europe [16].

thumbnail
Figure 1. Geographical location of the collection sites.

Sampled sites where G. holbrooki was not found are indicated with empty circles. Grey-scaled pie charts (white, light grey, dark grey, and black) represent mean proportional ancestry of every sampled site attributed to each cluster inferred by STRUCTURE. Watersheds are colored. Dotted lines represent geographical barriers indicated by BARRIER and the letters indicate the order in which the program detected these barriers. Location codes are presented in Table 1. (A detailed map with information about the road network is available at: http://mapsengine.google.com/map/viewer?mid=zgd4mwb-ESLE.kSE9TUfF2uQ4).

https://doi.org/10.1371/journal.pone.0082501.g001

Gambusia holbrooki specimens from Iberian rivers were collected from the riverbank using dip nets. All samples were collected from July to August 2010 and, to minimize any seasonal effects on population demography, only adult individuals born during the spring of the same year were selected by discarding females with a standard length of less than 2.5 cm and more than 3.5 cm, and males with a body length of less than 2.0 cm [38]. Individuals were classified as adult males if a fully formed gonopodium was present, and as females if not. We attempted to collect 20 males and 20 females from each site; however, adult fish availability modified this ratio (Table 1). Whole fish were euthanized by lethal sedation in situ, and then preserved in 96% ethanol until DNA was extracted at the laboratory.

DNA extraction and microsatellite analyses

For each fish that was collected, genomic DNA was isolated from the caudal muscle using the Realpure Genomic DNA extraction toolkit (Durviz SL, Valencia, Spain) following the manufacturer’s instructions. Genomic DNA was stored at –20°C until further use in Polymerase Chain Reactions (PCRs). Variation was analyzed at 11 loci (Pooc-G49, Mf13, and Gafμ3, Gafμ5, Gafμ6, Gafμ7, Gaaf7, Gaaf9, Gaaf10, Gaaf13, and Gaaf15), with two optimized multiplex PCR as described in Diez-del-Molino et al. (submitted). Both multiplex PCR were conducted under the same conditions: 30 µl of reaction volume containing 5–15 ng genomic DNA, 0.34 µM of each primer, 200 µM dNTPs, 1.5 mM MgCl2, and 0.75 units of Taq polymerase. The PCR cycling conditions were as follows: initial denaturation at 94°C for 3 min, followed by 35 cycles of 30 s at 94°C, 90 s at 60°C, 90 s at 72°C, and ending with a final extension of 10 min at 72°C. Forward primers were fluorescently labeled, and genotype peaks were resolved on a 3130 Genetic Analyzer and using GeneMapper 4.0 software (Applied Biosystems, Foster City, CA, USA).

Genetic diversity within locations

Genetic diversity within each study site was estimated from direct counts as the mean observed heterozygosity (HO) and the number of alleles per locus (A). Genetic diversity was also measured using the estimated expected heterozygosity (HE) and allelic richness (r) from allele frequencies using FSTAT 2.9.3 [39]. Using GENEPOP 4.0 [40], we measured the Hardy-Weinberg equilibrium (HWE) at each site, and tested for linkage disequilibrium between all pairs of loci. We corrected for multiple comparisons using the sequential Bonferroni test [41]. The presence of null alleles was detected using MICROCHECKER 2.2.3 [42], and their frequencies were estimated in FREENA [43]. We tested for recent population bottlenecks at the study sites using BOTTLENECK 1.2.02 [44].

Genetic structure within and among rivers

Pairwise population differentiation (FST) and significance values were calculated using FSTAT software. To assess the relevance of stepwise mutations on population differentiation (RST), an allele permutation test was performed with 1000 randomizations that simulated the distribution of allele sizes and RST values using SPAGEDI version 1.1 [45]. Allele richness and gene diversity patterns (HE and FST) within basins were compared among basins using permutation tests in FSTAT (1000 permutations). Non-parametric Wilcoxon signed-rank tests were used to compare allele richness and gene diversity (HE) between upstream and downstream collections within each study basin.

Isolation-by-distance (IBD) within and among watersheds was estimated from the correlation between genetic and geographical distance matrices among sampling sites. We used geographical distances, rather than hydrographical distances; because natural dispersal within linear river basins or human-mediated translocation by road could be involved in the connectivity between locations, given their geographical proximity (less than 50 km on average) and anastomosed road network (Fig. 1). The geographical distances between sample sites were estimated using Google Earth. Pairwise genetic differentiation was linearized as FST/(1-FST) and geographical distance was log-transformed for these analyses [46]. Significance was determined by Mantel tests with 10000 permutations using the IBD Web service 3.15 [47]. Additional information was obtained from the regression analyses of the estimates of the effective number of migrants (Nm) between populations pairs (Nm =  (1 – FST)/4FST) and their geographical distances (both variables log-transformed). Negative relationships indicate IBD, and the slope (b) of the linear regression (Log(Nm)  = a + bLog(d), where d equals geographical distance) is –1 for one-dimensional stepping stones models and –0.5 for the two-dimensional models [48].

The minimum number of homogeneous units (K) over sampled individuals was estimated using the MCMC method in STRUCTURE 2.3.3 [49]. Runs for each possible K (1 to 15) were repeated 10 times. Each run used a burn-in of 40000 iterations, a run length of 100000 iterations, and the model of independent allele frequencies. The most likely value of K was selected following Evanno [50]. The group-level Bayesian analysis in BAPS 5.4 [51] grouped populations that frequently exchanged individuals. BAPS analyses were repeated 10 times, with the maximum number of clusters set to 15. While STRUCTURE results tend to be conservative in the number of clusters detected providing ancestral information related with the history of introductions of the species, BAPS performed better in clustering together populations with recent gene flow (e.g., [52]). In addition, genetic differentiation among populations was depicted by two-dimensional plots from the principal components analysis (PCA) of the allele frequencies matrix in GENALEX 6.4.

Major genetic discontinuities in the study area were assessed using Monmonier’s algorithm in BARRIER 2.2 [53], which detects hidden barriers to gene flow among sites according to their geographical coordinates and relative genetic differentiation (FST). These analyses were conducted using the FST matrices from single-locus information corrected by the presence of null alleles (FREENA software). We identified the main barriers for each locus, and only retained those confirmed by at least six loci.

Analyses of molecular variance (AMOVA) were conducted in ARLEQUIN 3.5 [54]. Two hierarchical models were tested for partitioning the genetic diversity into three levels: within locations, among locations within regions, and among regions. The first AMOVA model assumed a pure hydrographical pattern of population hierarchy (watersheds  =  regions). Another AMOVA grouped locations according to the main clusters identified by STRUCTURE (clusters  =  regions).

Gene flow

Contemporary migration rates among populations were estimated by using the Bayesian inference as implemented in BAYESASS 3.0 software [55], which is a method that does not assume migration-drift or Hardy-Weinberg equilibrium. A total of 5×106 iterations were performed until the MCMC chains reached stationarity (i.e., constant over time). Migration parameters were estimated by sampling every 1000 iterations after a burn-in of 106 iterations. Delta values were adjusted following the BAYESASS manual recommendations. Five runs using different starting points were performed, and the results with the highest likelihood were retained.

Results

Diversity within locations

At the invaded locations, all 11 microsatellite loci were polymorphic, ranging in variability from just two alleles (Mf13, Gaaf15, and Gaaf9 loci) to nine (Gaaf13 locus). Average allelic richness (r) ranged from 2.45 in ES (Fluvià River) to 3.31 in CO (Ter River) (Table 2). The observed heterozygosity (HO) ranged from 0.336 in OF (Fluvià River) to 0.475 in PP (Fluvià River), and the expected heterozygosity (HE) ranged from 0.345 in BA (Fluvià River) to 0.500 in CE (Muga River). At the studied American collections, diversity levels averaged 4.45 for allelic richness, 0.465 for HO, and 0.557 for HE. After adjusting for differences in sample sizes, FSTAT permutation tests demonstrated lower allele richness (P = 0.009) and HE (P = 0.010) in the invaded Spanish locations compared to the potential American sources. Non-significant FST differentiation was detected between males and females (Table 3). Subsequent analyses were then performed pooling both sexes as a single collection for each location.

thumbnail
Table 2. Genetic diversity of Gambusia holbrooki in the study locations.

https://doi.org/10.1371/journal.pone.0082501.t002

thumbnail
Table 3. Genetic (FST corrected by the presence of null alleles, below the diagonal) and geographical distances (km, above the diagonal) between samples.

https://doi.org/10.1371/journal.pone.0082501.t003

Deviations from HWE were detected at seven Spanish locations after Bonferroni correction (Table 2). According to MICROCHECKER, null alleles were responsible for the observed positive FIS values. Significant null allele frequency was estimated at Gafμ5 (q = 0.071 in ES, Fluvià River; q = 0.089 in TO, Ter River), Gafμ6 (q = 0.149 in ES and q = 0.231 in CE, Muga River), Gafμ7 (q = 0.172 in MF, Fluvià River), Gaaf10 (q = 0.081 in CE, Muga River; q = 0.147 in CO, q = 0.126 in TO and q = 0.132 in TV, Ter River), and Gaaf15 (q = 0.231 in VM and q = 0.265 in CE, Muga River). No significant pairwise linkage disequilibria were observed. According to the BOTTLENECK analyses, heterozygosity excess relative to mutation-drift equilibrium was observed at CE (Muga River), BL, and MF (Fluvià River), and VE, and OY (Ter River); however, the allele-shift model test reported additional signals for a bottleneck at VE only.

Population divergence within and among basins

No significant differences between FST and RST estimates were observed (P = 0.345), indicating that local mutations have limited effects on population structure. Significant genetic differentiation was detected between almost all population pairs (Table 3), except for two neighboring sites in the Muga (VM and CE locations) and Fluvià (OF and BA locations) rivers (Table 3). Estimated average genetic differentiation among invaded Spanish locations was FST = 0.1641, with no substantial change in this estimate after correcting for null alleles (FST = 0.1642). The collection OY from the Onyar tributary in the Ter basin had the largest average pairwise FST (0.194).

Despite the different values on estimates of average allele richness, the average expected heterozygosity, and population differentiation observed in the Iberian rivers (Table 4), FSTAT permutation tests only indicated marginally reduced heterozygosity in the Fluvià basin compared to the Muga (P = 0.062) and Ter basins (P = 0.052). Overall, among-basin differences in genetic diversity (allele richness and heterozygosity) were non-significant between upstream and downstream locations. Marginal (P = 0.067) increased divergence among the upstream locations of the three basins was indicated (Table 4). Within-basin comparisons with Wilcoxon signed-rank tests indicated higher allele richness for just the downstream location of the Fluvià River (P = 0.042). Non-significant differences were obtained in the Muga River basin, and an unexpected higher richness was indicated at the upstream CL location compared to the downstream TV in the Ter River basin (P = 0.034). None of these changes in allele richness between upstream and downstream locations resulted in significant differences on the estimated amount of heterozygosity.

thumbnail
Table 4. Genetic diversity patterns within and among the studied locations.

https://doi.org/10.1371/journal.pone.0082501.t004

Genetic and geographical distance matrices were positively correlated across all three watersheds (r = 0.326, P<0.001). Within watersheds, only the Fluvià River displayed a significant correlation (Fluvià River, r = 0.694, P<0.01). Significant negative regression of the log-transformed effective number of migrants and geographical distances were detected for the whole data set of the studied locations in Spain. Furthermore, as the slope of the regression (b = –0.629) was closer to –0.5 than to –1, a two-dimensional stepping stone model better explained the population relationships (Fig. 2). Within watersheds, regression analyses of the number of migrants and distances were only significant for the Fluvià River, where a slope of –1.101 supported a one-dimensional stepping stone model. Marginal significance (P = 0.052) was obtained for the whole Ter River basin (b = –0.455); however, this significance disappeared when the analysis focused on the mainstem of the river (CL, VE, TO, and TV locations).

thumbnail
Figure 2. Linear regression of estimates of the effective number of migrants (Nm) and geographical distances between population pairs (both variables log-transformed, see Methods).

https://doi.org/10.1371/journal.pone.0082501.g002

Although the number of clusters identified varied among analyses, the clusters basically reflected population relationships according to the extent of drainage connection. Evanno’s method indicated four STRUCTURE clusters (Fig. 3a). Cluster 1 grouped individuals from all sites in the Muga river basin, while cluster 2 grouped all samples from the Fluvià River watershed, except for PP. The third cluster grouped all samples from the Ter River populations, except the Onyar tributary, which was assigned to cluster 4. The coastal locations of EP, PP, and TV showed a remarkable degree of cluster admixture. The Monmonier’s algorithm of BARRIER identified four barriers supported by at least six loci (Fig. 1). The first (a) and third (c) barriers reflected the isolation of the three watersheds, with the exception of the PP location in the Fluvià River basin, which was grouped with the Muga River collections. The second (b) and fourth (d) barriers reflected the distinct genetic composition of the OY and BY locations separated from the mainstem of the Ter River by river transects where mosquitofish were not detected during our surveys. Overall, these results mainly agreed with the population relationships depicted by the two principal axes of the PCA analysis (Fig. 4). The first axis explained 52.9% of the allelic variance, and clearly differentiated the Fluvià River collections from the rest, with the exception of PP. The second axis (17.1%) separated the Muga and Ter River basins from the singular population of OY. BAPS identified 10 homogenous units within the study region, basically indicating that each collection represented a single panmictic group (Fig. 3b). Only intra-basin locations that had the largest estimates of current gene flow (Table 5) were grouped together; specifically, VM and CE in the Muga basin, BL and ES and OF and BA in the Fluvià basin, and CL, VE, and TO in the Ter basin.

thumbnail
Figure 3. Bayesian analyses of population structure.

Analyses were carried out with (a) STRUCTURE and (b) BAPS in the Iberian G. holbrooki populations. In (a) each individual is represented as a vertical bar partitioned into segments of different color according to the proportion of the genome belonging to each of the four identified clusters (K = 4). In (b) each location shows a different color according to the cluster to which it belongs. Location codes are presented in Table 1.

https://doi.org/10.1371/journal.pone.0082501.g003

thumbnail
Figure 4. Principal component analysis (PCA) showing the relationships among the studied G. holbrooki populations.

Samples are projected onto the plane formed by the first two principal axes. The first factor explained the 52.9% of total variance, the second 17.1%, and the third 13.5%. Empty circles indicate positive values of the third axis, while grey circles indicate negative values. Location codes are presented in Table 1.

https://doi.org/10.1371/journal.pone.0082501.g004

thumbnail
Table 5. BAYESASS estimated migration rates among locations.

https://doi.org/10.1371/journal.pone.0082501.t005

Hierarchical AMOVAs revealed that the genetic variance was significant at all levels, with most of the variance being attributed to individuals within locations (80.9–82.9%). In the hydrographical model, the variance assigned to divergence among populations within basins (8.8%) was smaller compared to the variance among river basins (10.3%). This pattern reflected the above noted population divergence among and within drainages. In the AMOVA based on the four STRUCTURE clusters, the proportion of genetic variance explained within clusters decreased to 6.8%, while the variance among these clusters increased (12.3%), probably reflecting the distinctiveness of the OY collection from locations in the mainstem of the Ter River.

Overall, contemporary dispersal rates indicated an average of 21.4% of immigrant individuals at each location (range 11–31%, Table 5). As a conservative rule, we only discussed the 5% of highest estimates (11 out of 210 values). Within this framework, the most significant estimates of dispersal rates were basically downstream within rivers. Only at the ES site, and particularly at the OF site, a significant proportion of individuals were upstream immigrants from BA. BAYESASS also indicated the presence of gene flow among river basins, particularly from EP (Muga River) to PP (Fluvià River), and from the Onyar River (OY, Ter basin) to PP (Fluvià River).

Discussion

Genetic diversity and invasive potential retained in Spanish basins

In colonized territories, the level of genetic diversity of the invading species is expected to be reduced compared to original sources as a result of founder events [4], [7], [33], [56]. This phenomenon has been suggested for European mosquitofish populations when compared against the American collections from Florida [57]. For a more accurate evaluation of the effect of founder events on the genetic diversity and evolutionary potential of the Spanish populations, we compared the level of genetic diversity in the Spanish populations against those observed in the American populations considered to be the potential sources of the fish that were introduced to Europe. Previous studies have shown that the haplotype Hol1 is almost fixed in the Spanish collections and the American populations of Brunswick and Potomac River [16]. In addition, the Potomac River collection was the most closely related to European mosquitofish in a survey based on six microsatellites [58]. Genetic diversity declines by a factor of (1-1/2Ne) per generation during a founder effect depending on the effective number (Ne) of introduced individuals (see for instance [4]). Historical records indicate that just 12 individuals of mosquitofish were introduced to Spain [59]. If we consider the best case scenario for diversity retention involving just a single generation founder effect with an effective population size of 12 individuals, the population should have preserved around 95% of the original genetic diversity, or even more if some of the specimens were gravid females, because multiple paternity increases the effective size in mosquitofish populations [23]. With no relevant effects of mutations on population structure (FST = RST), all of the diversity present in the invaded range should be attributed to the population sources of the invasion. Based on the average diversity in the two American collections studied here (HE = 0.522, Table 4), the estimated diversity at these source locations agreed with the estimated total diversity at the Spanish study region (HT = 0.522). This observation indicates minimal, if any, loss of genetic diversity in the introduced populations of the Iberian Peninsula. Nevertheless, a significant reduction in allele richness was detected at the invaded Spanish locations, because this parameter is more sensitive to bottlenecks compared to average heterozygosity [24], [60]. It is therefore likely that introduced Spanish populations have not substantially reduced their evolutionary potential compared to American sources, because the levels of additive variance might still be less sensitive to bottlenecks compared to neutral variation (see reviews in [36], [61]). For example, introduced Australian populations of the guppy (Poecilia reticulata) showed strong genetic bottlenecks in genetic diversity when measured with neutral markers; yet, these populations retained substantial additive variation [62].

Lower neutral genetic variation is often detected in populations at the limit of the distribution range [63]. Reduced diversity and singular mtDNA haplotypes of G. holbrooki in northern American drainages indicated postglacial colonization from refuge populations in Georgia or Florida [64]. A recent work based on microsatellite variation confirmed an important reduction in allele richness (up to 50%) and heterozygosity (up to 30%) of American G. holbrooki populations in North Carolina and northward compared to populations that occurred to the south in South Carolina and Florida [58]. Nevertheless, peripheral populations often display greater stress-adaptation favoring subsistence in unstable environments [65]. Hence, available information about species with broad distributions indicates that less-stable habitats within native ranges serve as frequent sources of invasive populations (see [66]). If this was the case for mosquitofish, marginal populations of the northward range of America might have already acquired the evolutionary changes to be invasive during the postglacial period, as far as substantial additive variation could be retained during related founder effects despite losses of neutral genetic variation. Therefore, the American mosquitofish sources used in the European introduction might represent an “invasive bridgehead”. As defined by Lombaert et al. [67], invasive bridgeheads are particularly successful invasive populations that serve as the source of colonists for remote new territories. Because genetic diversity in the Spanish populations was not significantly reduced during the introduction, enough additive variance to respond to novel selection pressures in these non-native environments was probably conserved, favoring the successful and quick expansion of the species throughout the entire Mediterranean basin documented in historical records [20].

Dispersal patterns and population diversity in invaded locations

Precise historical records are not available about the introduction of mosquitofish to the study basins. While G. holbrooki was first introduced to southwestern Spain in 1921, it was absent from the study basins in 1942 when insecticides (DDT) replaced mosquitofish as the major agent against malaria vectors. Malaria was eradicated in 1964 from Spain; however, mosquito control, including mosquitofish introductions, continued (reviewed in [68]). Mosquitofish were apparently introduced to Lake Banyoles between 1952 and 1964, after they had already become established in other parts of the study watersheds [69]. Given that mosquitoes are abundant in marshlands dominating the lowlands of the three river basins [70], mosquitofish were probably first introduced into these lowland areas. The mosquitofish in the study river basins probably originated from well-established populations in central and southern Catalonia, such as the deltas of Llobregat and Ebro rivers, where mosquitofish were already present by 1942 [68]. In the Ebro River, which is located around 300 km south of our study area, mosquitofish populations exhibit similar levels of total diversity (HT = 0.532), with this diversity mainly being distributed within locations (HE = 0.523)(Diez-del-Molino et al. unpublished).

Significant genetic divergence among study locations indicated the isolation of current mosquitofish populations both within and among the three studied basins. According to Smith et al. [29], American mosquitofish populations along a river basin displayed a pattern of population divergence resulting from genetic drift and gene flow. However, some complex microgeographic patterns were also present as a result of interactions between dramatic demographic fluctuations and breeding structures complicated by multiple insemination and differential sex and cohort dispersal ability [71]. Larger divergence among mosquitofish populations located in the upper reaches of the study basins might be related to founder events during dispersal along the basins, because the contribution of mutations to population structure was not significant. The average level of population diversity at these locations represented 85% of the total genetic diversity in the area. At each location, the stated percentage indicated 4–5 generations of bottlenecks from just 12 individuals in magnitude (1-1/2Ne per generation of diversity losses), or more bottleneck-generations with larger Ne. In fact, signals of recent bottlenecks were detected at 30% of the locations, and affected populations from all river basins.

In the lowlands, differentiation among G. holbrooki populations exhibited similar patterns to that observed for the endemic killifish Aphanius iberus, in which increased gene flow was observed between populations during floods [72], [73]. If flooding also connects mosquitofish populations, this process alone justified why barriers were not detected between the downstream mosquitofish populations in the Muga and Fluvià rivers, because the mouths of both rivers flow out of the same marshland (Aiguamolls de l’Empordà). In addition, substantial reductions in the population size of mosquitofish have been reported to recover within a few months after flooding [74]. Moreover, pregnant females might buffer associated genetic bottlenecks [59], [75]. In the basins studies here, larger population divergence (FST = 0.46) has been detected among remaining native populations of the three-spine stickleback (Gasterosteus aculeatus) [76]. While stickleback remains in unpolluted streams with abundant aquatic vegetation [77], mosquitofish are successful invaders of modified and disturbed habitats, such as ponds, irrigation ditches, and modified stream channels in urban areas [15]; such habitats allow increased gene flow among locations in invaded basins (e.g., Diez-del-Molino et al. unpublished).

BAYESASS indicated current relevant migration rates between some neighboring population pairs in all of the study basins. This phenomenon resulted in BAPS clustering the VM and CE locations in the Muga River, the BL and ES and the BA and OF in the Fluvià River basin, and the CL and VE in the Ter river basin. These location pairs were separated by a distance of 5.3 km on average (range 3.50–8.28), with significant dispersal occurring both downstream and upstream (BA to OF). High positive spatial autocorrelation of allele frequencies at hydrological distances of 6–150 km has been observed in American drainages (reviewed in [29]), indicating gene flow between distant locations within short time scales (few generations). In [78], the authors suggested that mosquitofish may disperse at rates greater than 800 m/day in unimpeded corridors. These observations indicate that the dispersal ability of mosquitofish is sufficient to colonize an entire basin from a single founder effect. Once a single population is established, further active upstream or passive downstream dispersal leads to the founding of new populations, and maintains high enough gene flow to preserve existing genetic diversity throughout all locations along the invaded river basin, and to overcome founder effects.

In the Spanish basins, isolation by distance was indicated among mosquitofish populations through the significant correlation between genetic and geographical distances in the whole territory. The significant negative relationship (b = –0.629) observed between the log-transformed effective number of migrants (Nm) and the geographical distances also supported a two-dimensional stepping stone model of gene flow. Thus, in addition to active linear and flood-mediated dispersal along river basins, human-mediated translocations between road-neighboring populations from separate basins have probably contributed to the spread of mosquitofish in the area (for instance to OY and BY in the Ter river basin). In Australia, unreported and unregulated human-mediated dispersal has led to the introduction of mosquitofish to areas outside of its first sites of introduction [26]. In addition, a higher incidence of aquarists in more densely human populated areas (such as Girona in this study) might contribute to mosquitofish dispersal, through aquaria fish being discarded into urban ponds and river streams. Such practices have contributed to the dispersal of alien poeciliid species in Australia [79] and Spain [18]. Human-mediated translocations of endangered native species, such as A. iberus and G. aculeatus, are forbidden by Spanish laws directed towards protecting biodiversity. However, human-mediated dispersal might represent a major means of promoting gene flow between distant populations of mosquitofish, with such dispersal probably contributing towards retaining the high levels of genetic diversity within the populations of this species throughout the whole territory. It is also likely that human-mediated dispersal plays an important role in the maintenance of the invasive potential of these introduced populations, enabling them to outcompete the native fish.

Acknowledgments

The authors thank Cesca Rivas for assistance with laboratory procedures.

Author Contributions

Conceived and designed the experiments: DDM RA OV NS JGM EGB. Performed the experiments: DDM RA. Analyzed the data: DDM RA JGM. Contributed reagents/materials/analysis tools: DDM GCC RA OV. Wrote the paper: DDM JGM. Collected the samples: GCC EGB. Revised the manuscript: DDM RA OV NS JGM EGB.

References

  1. 1. Dudgeon D, Arthington AH, Gessner MO, Kawabata ZI, Knowler DJ, et al. (2006) Freshwater biodiversity: importance, threats, status and conservation challenges. Biol Rev Camb Philos Soc 81: 163–182.
  2. 2. Levine JM (2008) Biol Invasions. Curr Biol. Volume 18(2): 57–60.
  3. 3. Gozlan RE, Britton JR, Cowx I, Copp GH (2010) Current knowledge on non-native freshwater fish introductions. J Fish Biol 76: 751–786.
  4. 4. Lavergne S, Molofsky J (2007) Increased genetic variation and evolutionary potential drive the success of an invasive grass. PNAS 104: 3883–3888.
  5. 5. Roman J, Darling JA (2007) Paradox lost: genetic diversity and the success of aquatic invasions. . Trends Ecol Evol. 22: 454–464.
  6. 6. Lockwood JL, Cassey P, Blackburn T (2005) The role of propagule pressure in explaining species invasions. Trends Ecol Evol 20: 223–228.
  7. 7. Dlugosch KM, Parker IM (2008) Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol 17: 431–449.
  8. 8. Kolbe JJ, Glor RE, Rodríguez Schettino L, Lara AC, Larson A, et al. (2004) Genetic variation increases during biological invasion by a Cuban lizard. Nature 431: 177–181.
  9. 9. Keller SR, Taylor DR (2010) Genomic admixture increases fitness during a biological invasion. J Evolution Biol 23: 1720–1731.
  10. 10. Lowe S, Browne M, Boudjelas S, De Poorter M. (2000) 100 of the World’s Worst Invasive Alien Species: A selection from the Global Invasive Species Database. ISSG, World Conservation Union (IUCN), 12pp.
  11. 11. Alcaraz C, Bisazza A, García-Berthou E (2008) Salinity mediates the competitive interactions between invasive mosquitofish and an endangered fish. Oecologia 155: 205–213.
  12. 12. Pyke GH (2008) Plague Minnow or Mosquito Fish? A Review of the Biology and Impacts of Introduced Gambusia Species. Annu Rev Ecol Evol S 39: 171–191.
  13. 13. Rincón P, Correas A, Morcillo F, Risueño P, Lobón-Cerviá J (2002) Interaction between the introduced eastern mosquitofish and two autochthonous Spanish toothcarps. J Fish Biol 61: 1560–1585.
  14. 14. Carmona-Catot G, Magellan K, García-Berthou E (2013) Temperature-Specific Competition between Invasive Mosquitofish and an Endangered Cyprinodontid Fish. PLoS ONE 8(1): e54734
  15. 15. Courtenay W, Meffe G (1989) Small fishes in strange places: a review of introduced poeciliids. In: Meffe G, Snelson F, editors. Ecology and Evolution of Live-Bearing Fishes (Poeciliidae). New Jersey: Prentice Hall, Englenwood Cliffs, Vol. 1991. pp. 319–331.
  16. 16. Vidal O, García-Berthou E, Tedesco P, García-Marín JL (2010) Origin and genetic diversity of mosquitofish (Gambusia holbrooki) introduced to Europe. Biol Invasions 12: 841–851.
  17. 17. García-Berthou E, Alcaraz C, Pou-Rovira Q, Zamora L, Coenders G, et al. (2005) Introduction pathways and establishment rates of invasive aquatic species in Europe. Can J Fish Aquat Sci 463: 453–463.
  18. 18. Ribeiro F, Elvira B, Collares-Pereira M (2008) Life-history traits of non-native fishes in Iberian watersheds across several invasion stages: a first approach. Biol Invasions 10: 89–102.
  19. 19. Frankham R (1995) Effective population-size adult-population size ratios in wildlife - a review. Genet Res 66: , 95–107.
  20. 20. Krumholz L (1948) Reproduction in the western mosquitofish, Gambusia affinis affinis (Baird & Girard), and its use in mosquito control. Ecol Monogr 18: 1–43.
  21. 21. Hernandez-Martich J, Smith M (1990) Patterns of genetic variation in eastern mosquitofish (Gambusia holbrooki Girard) from the piedmont and coastal plain of three drainages. Copeia 1990: 619–630.
  22. 22. Echelle A, Wildrick D, Echelle A (1989) Allozyme studies of genetic variation in Poeciliid fishes. In: Meffe G, Snelson F, editors. Ecology and Evolution of Live-Bearing Fishes (Poeciliidae). New Jersey: Prentice Hall, Englenwood Cliffs, Vol. 1991. pp. 217–234.
  23. 23. Zane L, Nelson WS, Jones AG, Avise JC (1999) Microsatellite assessment of multiple paternity in natural populations of a live-bearing fish, Gambusia holbrooki. J Evolution Biol 12: 61–69.
  24. 24. Spencer C, Neigel J, Leberg P (2000) Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks. Mol Ecol 9: 1517–1528.
  25. 25. Smith M, Smith M, Chesser R (1983) Biochemical genetics of mosquitofish. I. Environmental correlates, and temporal and spatial heterogeneity of allele frequencies within a river drainage. Copeia 1983: 182–193.
  26. 26. Ayres RM, Pettigrove VJ, Hoffmann AA (2010) Low diversity and high levels of population genetic structuring in introduced eastern mosquitofish (Gambusia holbrooki) in the greater Melbourne area, Australia. Biol Invasions: 3727–3744.
  27. 27. Barson NJ, Cable J, Van Oosterhout C (2009) Population genetic analysis of microsatellite variation of guppies (Poecilia reticulata) in Trinidad and Tobago: evidence for a dynamic source – sink metapopulation structure, founder events and population bottlenecks. J Evol Biol 22: 485–497
  28. 28. McClenaghan L, Smith M, Smith M (1985) Biochemical genetics of mosquitofish. IV. Changes of allele frequencies through time and space. Evolution 39: 451–460.
  29. 29. Smith M, Scribner K, Hernandez J, Wooten M (1989) Demographic, spatial, and temporal genetic variation in Gambusia. In: Meffe G, Snelson F, editors. Ecology and Evolution of Live-Bearing Fishes (Poeciliidae). New Jersey: Prentice Hall, Englenwood Cliffs, Vol. 220. pp. 235–257.
  30. 30. McElroy T, Kandl K (2011) Temporal Population Genetic Structure of Eastern Mosquitofish in a Dynamic Aquatic Landscape. J Hered 102: 678–687.
  31. 31. Congdon BC (1995) Unidirectional Gene Flow and Maintenance of Genetic Diversity in Mosquitofish Gambusia holbrooki (Teleostei: Poeciliidae). Copeia 1995: 162.
  32. 32. Hernandez-Martich J, Smith M (1997) Downstream gene flow and genetic structure of Gambusia holbrooki (eastern mosquitofish) populations. Heredity 79: 295–301.
  33. 33. Allendorf FW, Lundquist LL (2003) Introduction: population biology, evolution, and control of invasive species. Conserv Biol 17: 24–30.
  34. 34. Stepien CA, Brown JE, Neilson ME, Tumeo MA (2005) Genetic diversity of invasive species in the Great Lakes versus their Eurasian source populations: insights for risk analysis. Risk Anal 25: 1043–1060.
  35. 35. Frankham R (2010) Challenges and opportunities of genetic approaches to biological conservation. Biol Conserv 143: 1919–1927.
  36. 36. Blanchet S (2012) The use of molecular tools in invasion biology: an emphasis on freshwater ecosystems. Fisheries Manag Ecol 19: 120–132.
  37. 37. Trigo RM, Pozo-Vázquez D, Osborn TJ, Castro-Díez Y, Gámiz-Fortis S, et al. (2004) North Atlantic oscillation influence on precipitation, river flow and water resources in the Iberian Peninsula. Int J Climatol 24: 925–944.
  38. 38. Carmona-Catot G, Benito J, García-Berthou E (2011) Comparing latitudinal and upstream-downstream gradients: life history traits of invasive mosquitofish. Divers Distrib 17: 214–224.
  39. 39. Goudet J (2001) FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Available: http://www.unil.ch/izea/softwares/fstat.html.
  40. 40. Raymond M, Rousset F (1995) An exact test for population differentiation. Evolution 49: 1280–1283.
  41. 41. Rice WR (1989) Analyzing tables of statistical tests. Evolution 43: 223–225.
  42. 42. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538.
  43. 43. Chapuis MP, Estoup A (2007) Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol 24: 621–631.
  44. 44. Piry S, Luikart G, Cornuet JM (1999) BOTTLENECK: A Computer Program for Detecting Recent Reductions in the Effective Population Size Using Allele Frequency Data. J Hered 90: 502–503.
  45. 45. Hardy OJ, Vekemans X (2002) SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes 2: 618–620.
  46. 46. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219–1228.
  47. 47. Jensen JL, Bohonak AJ, Kelley ST (2005) Isolation by distance, web service. BMC Genet 6 13: v.3.23 http://ibdws.sdsu.edu/.
  48. 48. Slatkin M (1993) Isolation by distance in equilibrium and non-equilibrium populations. Evolution 47: 264–279.
  49. 49. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
  50. 50. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620.
  51. 51. Corander J, Marttinen P, Sirén J, Tang J (2008) Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinformatics 9: 539.
  52. 52. Vidal O, Sanz N, Araguas RM, Fernández-Cebrián R, Diez-del-Molino D, et al. (2012) SNP diversity in introduced populations of the invasive Gambusia holbrooki. Ecol Freshw Fish 21: 100–108.
  53. 53. Manni F, Guerard E, Heyer E (2004) Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by using Monmonier’s algorithm. Hum Biol 76: 173–190.
  54. 54. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol resources 10: 564–567.
  55. 55. Wilson GA, Rannala B (2003) Bayesian inference of recent migration rates using multilocus genotypes. Genetics 163: 1177–1191.
  56. 56. Suarez AV, Tsutsui ND (2008) The evolutionary consequences of biological invasions. Mol Ecol 17: 351–360.
  57. 57. Grapputo A, Bisazza A, Pilastro A (2006) Invasion success despite reduction of genetic diversity in the European populations of eastern mosquitofish (Gambusia holbrooki). Ital J Zool 73: 67–73.
  58. 58. Sanz N, Araguas RM, Vidal O, Diez-del-Molino D, Fernández-Cebrián R, et al.. (2013) Genetic characterization of the invasive mosquitofish (Gambusia spp.) introduced to Europe: population structure and colonization routes. Biol Invasions. doi: 10.1007/s10530-013-0456-5.
  59. 59. Nájera L (1944) Sobre la identificación de Gambusia holbrooki. Boletín de la Real Sociedad Española de Historia Natural 42: 51–55.
  60. 60. Leberg P (1992) Effects of population bottlenecks on genetic diversity as measured by allozyme electrophoresis. Evolution 46: 477–494.
  61. 61. Hänfling B (2007) Understanding the establishment success of non-indigenous fishes: lessons from population genetics. J Fish Biol71: 115–135.
  62. 62. Lindholm AK, Breden F, Alexander HJ, Chan WK, Thakurta SG, et al. (2005) Invasion success and genetic diversity of introduced populations of guppies Poecilia reticulata in Australia. Mol Ecol 14: 3671–3682.
  63. 63. Eckert CG, Samis KE, Lougheed SC (2008) Genetic variation across species’ geographical ranges: the central-marginal hypothesis and beyond. Mol Ecol 17: 1170–1188.
  64. 64. Scribner KT, Avise JC (1993) Cytonuclear genetic architecture in mosquitofish populations and the possible roles of introgressive hybridization. Mol Ecol 2: 139–149.
  65. 65. Hardie DC, Hutchings JA (2010) Evolutionary ecology at the extremes of species’ ranges. Environ Rev 18: 1–20.
  66. 66. Lee CE, Gelembiuk GW (2008) Evolutionary origins of invasive populations. Evol Appl 1: 427–448.
  67. 67. Lombaert E, Guillemaud T, Cornuet JM, Malausa T, Facon B, et al. (2010) Bridgehead effect in the worldwide invasion of the biocontrol harlequin ladybird. PloS one 5: e9743
  68. 68. Fernández-Astasio B (2002) La erradicación del paludismo en España: aspectos biológicos de la lucha antipalúdica. Tesis doctoral. Universidad Complutense de Madrid, España.
  69. 69. Garcia-Berthou E, Moreno-Amich R (2000) Introduction of exotic fish into a Mediterranean lake over a 90-year period. Fund Appl Limnol 149 (2): 271–284.
  70. 70. Boix D, Gascón S, Sala J, Badosa A, Brucet S, et al. (2008) Patterns of composition and species richness of crustaceans and aquatic insects along environmental gradients in Mediterranean water bodies. Hydrobiologia 597: 53–69.
  71. 71. Kennedy P, Kennedy M, Zimmerman E (1986) Biochemical genetics of mosquitofish. V. Perturbation effects on genetic organization of populations. Copeia 4: 937–945.
  72. 72. Garcia-Marin J, Jorde P, Ryman N, Utter F, Pla C (1991) Management implications of genetic differentiation between native and hatchery populations of brown trout (Salmo trutta) in Spain. Aquaculture 95: 235–249.
  73. 73. Araguas RM, Roldán MI, García-Marín JL, Pla C (2007) Management of gene diversity in the endemic killifish Aphanius iberus: revising Operational Conservation Units. Ecol Freshw Fish 16: 257–266.
  74. 74. Chapman P, Warburton K (2006) Postflood movements and population connectivity in gambusia (Gambusia holbrooki). Ecol Freshw Fish 15: 357–365.
  75. 75. Spencer CC, Neigel JE, Leberg PL (2000) Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks. Mol Ecol 9: 1517–1528.
  76. 76. Araguas RM, Vidal O, Pla C, Sanz N (2011) High genetic diversity of the endangered Iberian three-spined stickleback (Gasterosteus aculeatus) at the Mediterranean edge of its range. Freshwater Biol: 57: , 143–154 doi:10.1111/j.1365-2427.2011.02705.x.
  77. 77. Clavero M, Pou-Rovira Q, Zamora L (2009) Biology and habitat use of three-spined stickleback (Gasterosteus aculeatus) in intermittent Mediterranean streams. Ecol Freshw Fish 18: 550–559.
  78. 78. Alemadi SD, Jenkins DG (2007) Behavioral constraints for the spread of the eastern mosquitofish, Gambusia holbrooki (Poeciliidae). Biol Invasions 10: 59–66.
  79. 79. Lintermans M (2004) Human-assisted dispersal of alien freshwater fish in Australia. New Zeal J Mar Fresh 38 (3): 481–501.