Evolutionary genetics of invasive species

Biological invasions in the Antarctic: extent, impacts and implications

Alien microbes, fungi, plants and animals occur on most of the sub-Antarctic islands and some parts of the Antarctic continent. These have arrived over approximately the last two centuries, coincident with human activity in the region. Introduction routes have varied, but are largely associated with movement of people and cargo in connection with industrial, national scientific program and tourist operations. The large majority of aliens are European in origin. They have both direct and indirect impacts on the functioning of species-poor Antarctic ecosystems, in particular including substantial loss of local biodiversity and changes to ecosystem processes. With rapid climate change occurring in some parts of Antarctica, elevated numbers of introductions and enhanced success of colonization by aliens are likely, with consequent increases in impacts on ecosystems. Mitigation measures that will substantially reduce the risk of introductions to Antarctica and the sub-Antarctic must focus on reducing propagule loads on humans, and their food, cargo, and transport vessels.

Genetic variability and founder effect in the pitcher plant Sarracenia purpurea (Sarraceniaceae) in populations introduced into Switzerland: from inbreeding to invasion

BACKGROUND AND AIMS

The long-lived and mainly outcrossing species Sarracenia purpurea has been introduced into Switzerland and become invasive. This creates the opportunity to study reactions to founder effect and how a species can circumvent deleterious effects of bottlenecks such as reduced genetic diversity, inbreeding and extinction through mutational meltdown, to emerge as a highly invasive plant.

METHODS

A population genetic survey by random amplified polymorphism DNA markers (RAPD) together with historical insights and a field pollination experiment were carried out.

KEY RESULTS

At the regional scale, S. purpurea shows low structure (thetast=0.072) due to a recent founder event and important subsequent growth. Nevertheless, multivariate statistical analyses reveal that, because of a bottleneck that shifted allele frequencies, most of the variability is independent among populations. In one population (Tenasses) the species has become invasive and genetic analysis reveals restricted gene flow and family structure (thetast=0.287). Although inbreeding appears to be high (Fis >0.410 from a Bayesian estimation), a field pollination experiment failed to detect significant inbreeding depression upon F1 seed number and seed weight fitness-traits. Furthermore, crosses between unrelated individuals produced F1 seeds with significantly reduced fitness, thus showing local outbreeding depression.

CONCLUSIONS

The results suggest that, under restricted gene flow among families, the species may not only have rapidly purged deleterious alleles, but also have undergone some form of selection for inbreeding due to co-adaptation between loci.

Genetic analysis of complex demographic scenarios: spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm, to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker, demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Regional population structure of a widely introduced estuarine invertebrate: Nematostella vectensis Stephenson in New England

Nematostella vectensis is an infaunal anemone occurring in salt marshes, lagoons and other estuarine habitats in North America and the United Kingdom. Although it is considered rare and receives protection in England, it is widely distributed and abundant in the United States, particularly along the Atlantic coast. Recent studies suggest that both anthropogenic dispersal and reproductive plasticity may significantly influence the genetic structure of N. vectensis populations. Amplified fragment length polymorphism (AFLP) fingerprinting of individuals from nine populations in the northeastern United States indicates that stable populations are maintained by both asexual and sexual reproduction; in some cases asexually reproducing lineages exist within sexually reproducing populations. F statistics reveal extraordinarily high degrees of genetic differentiation between populations, even those separated by very short distances (less than 100 m). Genetic distances show little to no correlation with geographical distances, consistent with a role for sporadic, geographically discontinuous dispersal coupled with limited gene flow. No single genotype was found at more than one site, despite apparent homogeneity of habitat. In contrast with reported genotypic distributions for Nematostella in the United Kingdom, where a single clonal genotype dominates at multiple sites through southern England, our data thus fail to support the hypothesis of a general-purpose genotype in the northeastern United States. However, they are consistent with important roles for reproductive plasticity, sporadic introductions and complex local population dynamics in determining the global and regional distribution of this species.

Adapting to an invasive species: toxic cane toads induce morphological change in Australian snakes

The arrival of invasive species can devastate natural ecosystems, but the long-term effects of invasion are less clear. If native organisms can adapt to the presence of the invader, the severity of impact will decline with time. In Australia, invasive cane toads (Bufo marinus) are highly toxic to most snakes that attempt to eat them. Because snakes are gape-limited predators with strong negative allometry for head size, maximum relative prey mass (and thus, the probability of eating a toad large enough to be fatal) decreases with an increase in snake body size. Thus, the arrival of toads should exert selection on snake morphology, favoring an increase in mean body size and a decrease in relative head size. We tested these predictions with data from specimens of four species of Australian snakes, collected over >80 years. Geographic information system layers provided data on the duration of toad exposure for each snake population, as well as environmental variables (latitude, precipitation, and temperature). As predicted, two toad-vulnerable species (Pseudechis porphyriacus and Dendrelaphis punctulatus) showed a steady reduction in gape size and a steady increase in body length with time since exposure to toads. In contrast, two species at low risk from toads (Hemiaspis signata and Tropidonophis mairii) showed no consistent change in these morphological traits as a function of the duration of toad exposure. These results provide strong evidence of adaptive changes in native predators as a result of the invasion of toxic prey.

Colonization of America by Drosophila subobscura: association between Odh gene haplotypes, lethal genes and chromosomal arrangements

The colonization of America by Drosophila subobscura has been a unique exper iment in nature that has allowed us to explore the effects of evolution on a continental scale. To analyze this evolutionary event, nucleotide sequences of the Odh (Octanol dehydrogenase) gene were obtained for 43 lethal chromosomal lines from colonizing populations of North America and 5 from South America, in addition to 5 chromosomal lines from Europe with different viabilities and 2 from laboratory marker stocks. Since 10 different Odh haplotypes were found in America, the minimum number of colonizers would be 5 (or 3 mated females). Only one Odh haplotype was found in American O(5) inversions confirming that only one copy of this inversion was included among the sample of colonizers. The same Odh haplotypes were detected in association with the same chromosomal arrangements and with identical lethal genes in both North and South America indicating that exactly the same chromosome types reached both hemispheres. These observations indicate that the two continental colonizations are not independent. They are derived from the same colonization event. The population from which the colonization started should contain the O(5) inversion, a non-negligible frequency of the O(3+4+7) arrangement and all other arrangements found in America. So far the only populations that fulfill all these requirements are those from Greece, indicating that these populations can be considered good candidates as a starting point for an in depth analysis of the origin of the American colonization by D. subobscura.

Genetic variation increases during biological invasion by a Cuban lizard

A genetic paradox exists in invasion biology: how do introduced populations, whose genetic variation has probably been depleted by population bottlenecks, persist and adapt to new conditions? Lessons from conservation genetics show that reduced genetic variation due to genetic drift and founder effects limits the ability of a population to adapt, and small population size increases the risk of extinction. Nonetheless, many introduced species experiencing these same conditions during initial introductions persist, expand their ranges, evolve rapidly and become invasive. To address this issue, we studied the brown anole, a worldwide invasive lizard. Genetic analyses indicate that at least eight introductions have occurred in Florida from across this lizard's native range, blending genetic variation from different geographic source populations and producing populations that contain substantially more, not less, genetic variation than native populations. Moreover, recently introduced brown anole populations around the world originate from Florida, and some have maintained these elevated levels of genetic variation. Here we show that one key to invasion success may be the occurrence of multiple introductions that transform among-population variation in native ranges to within-population variation in introduced areas. Furthermore, these genetically variable populations may be particularly potent sources for introductions elsewhere. The growing problem of invasive species introductions brings considerable economic and biological costs. If these costs are to be mitigated, a greater understanding of the causes, progression and consequences of biological invasions is needed.

Hybridization and contemporary evolution in an introduced cichlid fish from Lake Malawi National Park

Rapidly evolving systems offer the chance to observe genetic and phenotypic change in real time. We exploit a well-characterized introduction of cichlid fish into Lake Malawi National Park to document a short history of habitat colonization and the evolution of genes and colour pattern. In the early 1960s, a fish exporter introduced individuals of Cynotilapia afra to a single site (Mitande Point) of Thumbi West Island and, as late as 1983, the species was confined to this location. In 2001, C. afra had colonized the entire perimeter of Thumbi West. In July of that year, we sampled C. afra individuals from six sites around the island and scored variation in dorsal fin colour as well as allelic diversity at six microsatellite loci. We found that, in two decades, C. afra had diverged into genetically distinct, phenotypically different northern and southern populations. We observed a high proportion of hybrids between the introduced C. afra and the native Metriaclima zebra on the southern coast of Thumbi West, and speculate that hybridization is facilitated by low water clarity at these windward sites. The short history of C. afra at Thumbi West is a microcosm of contemporary evolutionary divergence and may provide the opportunity to study the process from start to finish in genetic detail.

Comparing the freezing susceptibility of third-instar larvae of Gnorimus variabilis (Cetoniidae: Trichiinae) from three distant geographical regions

We compared the freezing susceptibility of three populations of Gnorimus variabilis (L., 1758) (Coleoptera: Cetoniidae) sampled from three distant locations in France. To separate the effects of habitat quality from those of genetics, we reared imagos from all field populations in a common garden experiment until the emergence of the third-instar larvae. The local climate appeared to determine the cold hardiness of the third-instar larvae, which live exclusively in cold seasons. The geographical location had an effect on the frost resistance (supercooling point) of the third-instar larvae of G. variabilis. We found no significant differences between the supercooling points of the populations originating from separate latitudes but from the same longitude, Sare (–10.9 ± 4.1 °C) and Husson (–11.5 ± 3.8 °C). Nonetheless, significant differences occurred between the larvae originating from the same northern latitude but from separate longitudes, Husson and La Robertsau (–17.8 ± 2.9 °C). It is the first time that a highly significant difference in the supercooling point of one stage within a single species has been observed along geographical gradients. Moreover, 19 of the 24 larvae originating from La Robertsau were alive after they were kept for 22 h at –10 °C compared with only 7 of the 24 larvae originating from Sare. When the three populations were reared in the laboratory under the same experimental conditions, the differential cryoresistance was preserved. It is likely that the greater freezing resistance found in the third-instar larvae of G. variabilis from La Robertsau could have a genetic component.

Interspecific hybridization among Hieracium species in New Zealand: evidence from flow cytometry

Hieracium pilosella: (Asteraceae) was accidentally introduced to New Zealand about 100 years ago. Since then it has become an aggressive weed, and an unexpected degree of genetic and genome size variation has been detected; features that might result from interspecies hybridization. We investigated the possibility that H. pilosella has hybridized with related taxa. Of the four other subgenus Pilosella species introduced to New Zealand, H. praealtum is the most abundant and, on morphological and distributional evidence, most likely to be the other parent. Flow cytometry was used to estimate relative genome size for 156 Hieracium plants collected from the wild. Plants assigned to either parental or hybrid morphotypes were found to comprise tetraploid and pentaploid individuals using genome size measurements, and this was confirmed with direct mitotic chromosome counts for a subset of plants. The haploid DNA content of H. praealtum was approximately 22% larger than that of H. pilosella. Putative hybrids that were tetraploid had mean genome sizes equivalent to two H. pilosella and two H. praealtum haploid chromosome sets, implying they were hybrids arising from the fertilization of two reduced gametes. Similar results were obtained from tetraploid hybrids produced by controlled pollination. However, the majority of field hybrids were pentaploid with a genome size equivalent to four H. pilosella and one H. praealtum haploid chromosome sets. We infer that these are not first-generation hybrids but represent successful backcrossing with H. pilosella and/or hybrid-hybrid crossing, and that sexual tetraploid hybrids have been the parents. We note that populations putatively of H. pilosella often comprise apomictic pentaploid hybrids. Significantly, our data indicate the emergence of sexual hybrids that provide further opportunity for gene flow among taxa in this complex.

Ecological determinants of genetic diversity in an expanding population of the shrub Myrica cerifera

The ecological mechanisms that contribute to the acquisition of genetic diversity in an expanding population of the shrub, Myrica cerifera, on an island habitat were investigated. Genealogical reconstruction was used to assess the contribution of early reproductive colonists to subsequent recruitment. In addition, through determination of parentage, the source of recruiting seedlings was identified and the contribution of seed and pollen dispersal into the colonizing sites was inferred. The relative contribution of different sources of gene flow was determined directly and an investigation was made into how variability in breeding patterns may have contributed to observed levels of genetic variability. It was expected that early colonists that could flower would contribute to subsequent recruiting cohorts, and that the limited number of such early reproductive colonists would lead to variance in mating success, inbreeding, or bottlenecks which could reduce genetic diversity and increase genetic differentiation among subsequent recruiting cohorts. Analyses of parentage (with paternity exclusion probability > 95%) for all recruiting plants demonstrated that in fact, there was little contribution by the early reproductive colonists to subsequent cohorts, and that immigration from outside the study sites in the form of seed dispersal accounted for over 94% of the recruitment in the study plots, with pollen dispersal accounting for less than 3% gene flow. No genetic bottleneck or evidence of reproductive skew in the recruiting cohorts were found, suggesting that propagule dispersal was from many source individuals in other established populations.

Dispersal of the Ponto-Caspian amphipod Echinogammarus ischnus: invasion waves from the Pleistocene to the present

The geographical range of the amphipod crustacean Echinogammarus ischnus has expanded over the past century from the Ponto-Caspian region to Western Europe, the Baltic Sea, and the Great Lakes of North America. The present study explores the phylogeographic patterns of this amphipod across its current distribution, based on an examination of nucleotide diversity in the mitochondrial cytochrome c oxidase subunit I (COI) gene. Marked genetic divergence exists among populations of E. ischnus from the Black and Caspian Seas, as well as those from the drainage system of the Black Sea. This divergence suggests the prolonged geographic isolation of these native populations, reflecting the limited dispersal capability of E. ischnus. By contrast, invading populations are characterized by a lack of genetic variation; a single mitochondrial genotype of Black Sea origin has colonized sites from the Rhine River to North America. The dispersal pattern in E. ischnus is very similar to that in the Ponto-Caspian cladoceran Cercopagis pengoi. Despite their contrasting life history strategies, these invading species followed the same route of invasion from the northern Black Sea to the Baltic Sea region, and subsequently to North America.

Chloroplast DNA diversity of Hieracium Pilosella (Asteraceae) introduced to New Zealand: reticulation, hybridization, and invasion

The European hawkweed Hieracium pilosella is a successful invader and a troublesome weed in New Zealand. The systematics of the genus Hieracium is extremely complex and contentious, probably due to recent speciation, hybridization, polyploidy, and diverse reproductive strategies. In the first chloroplast DNA survey of the group, we sequenced 285 plants (including H. pilosella and 12 other species of subgenus Pilosella) from New Zealand and Europe for 900 bp of trnL-trnF. Eleven haplotypes were identified with much sharing among species. Three haplotypes (A, D, G) were found in seven, three, and four species, respectively, but two species (H. lactucella and H. auricula) had single, private haplotypes. Our cpDNA data for subgenus Pilosella are consistent with the group's having incomplete lineage sorting and/or recent reticulate evolution. Six haplotypes were identified in H. pilosella, four of these unique to this taxon in our sample. In New Zealand, haplotype A was common and occurred in plants of different ploidy (i.e., 4×, 5×, 6×), whereas haplotypes C, B, and M were restricted to 4×, 5×, and 6× plants, respectively. The distribution of haplotype variation suggests that some or all of the H. pilosella seeds accidentally introduced into New Zealand probably came from east Europe rather than the United Kingdom and that a minimum of four lineages were introduced. Within New Zealand, hybridization of H. pilosella with a related taxon (probably H. praealtum) has occurred at least three times, involving both obligate sexual tetraploids and facultative apomictic pentaploids of H. pilosella.

The hierarchical spatial distribution of chloroplast DNA polymorphism across the introduced range ofSilene vulgaris

Silene vulgaris was introduced into North America sometime prior to 1800. In order to document the population structure that has developed since that time, collections were made from 56 local populations distributed among 9 geographical regions in eastern North America. Individual plants were characterized for chloroplast DNA (cpDNA) haplotype by restriction fragment size analysis of four noncoding regions of cpDNA amplified by polymerase chain reaction. A total of 19 cpDNA haplotypes were detected using this method. The overall gene diversity of 0.85 is quite similar to the diversity detected in these same regions of cpDNA in a previously published sample of S. vulgaris taken from across much of Europe. The spatial distribution of the North American cpDNA diversity was quantified by hierarchical F-statistics that partitioned the genetic variance into variation among local populations within regions, and variation among regions. The average FST among populations within regions was 0.66 and the FST among regions was 0.09. The among-region variation was due to both differences among regions in the frequency of two most common haplotypes, and to the presence of a number of region-specific haplotypes. In order to test for isolation by distance at the regional level, FST values were calculated for all possible pairs of regions, and regressed against the geographical distance between those regions. There was no evidence for isolation by distance. It is suggested that the local population structure is generated by recent extinction/colonization dynamics, and that the among-region structure reflects demographic events associated with range expansion following introduction to North America.

Microsatellite variation within and among North American lineages of Phragmites australis

Over the past century, the spread of the common reed (Phragmites australis) has had a dramatic impact on wetland communities across North America. Although native populations of Phragmites persist, introduced invasive populations have dominated many sites and it is not clear if the two types can interbreed. This study compares patterns of differentiation in 10 microsatellite loci among North American and European Phragmites individuals with results obtained from sequencing of noncoding chloroplast DNA. Three population lineages (native, introduced and Gulf Coast) were previously identified in North America from chloroplast DNA and similar structuring was found in the nuclear genome. Each lineage was distinguished by unique alleles and allele combinations and the introduced lineage was closely related to its hypothesized source population in Europe. Size homoplasy and diagnostic base substitutions distinguishing lineages were evident at several loci, further emphasizing that native, introduced and Gulf Coast North American Phragmites lineages are genetically distinct. Gene flow between lineages was low and invasive introduced populations do not represent a hybrid population type.

Hybridization and Rorippa austriaca (Brassicaceae) invasion in Germany

Introgressive hybridization between the invasive Rorippa austriaca and the native R. sylvestris in Germany has been studied using chloroplast DNA (trnL intron) and amplified fragment length polymorphism. Three hybrid zones between the invasive and native species were located in the Ruhr Valley (Mülheim) and at the River Main near Würzburg (Randersacker, Winterhausen). In each hybrid zone hybridization was indicated by additivity of region-specific amplified fragment length polymorphism markers proving independent hybridization events. The hybrids were either morphologically intermediate (R. x armoracioides) or were close to R. sylvestris. The trnL intron of R. austriaca is characterized by a species-specific deletion. This diagnostic chloroplast marker of R. austriaca was detected in three individuals of R. sylvestris providing evidence for introgression of the invasive chloroplast into the native species. Bidirectional introgression of R. austriaca markers into R. sylvestris and of R. sylvestris markers into R. austriaca was detected in the amplified fragment length polymorphism analysis. Some of the invasive R. austriaca populations showed high within-population variation. A possible association among introgression, within-population variation and invasion success is discussed. The morphologically intermediate hybrid R. x armoracioides is currently spreading in northern Germany. It forms large populations without its parent species R. austriaca and R. sylvestris. It is concluded that hybridization between invasive R. austriaca and native R. sylvestris may lead to the evolution of a new invasive species R. x armoracioides.

Infection of New- and Old-World Aedes albopictus (Diptera: Culicidae) by the intracellular parasite Wolbachia: implications for host mitochondrial DNA evolution

Wolbachia are cytoplasmically inherited, endosymbiotic bacteria known to infect a wide variety of arthropods. Polymerase chain reaction (PCR) amplification of the Wolbachia surface protein (wsp) gene was used to assay the infection of geographically disparate populations of Aedes albopictus (Skuse) by Wolbachia. Nine North American, four South American, one Hawaiian, and four Old World populations of A. albopictus were all doubly infected with both the wAlbA and wAlbB strains of Wolbachia. A 365-bp region of the wAlbA wsp gene was sequenced from seven geographically disparate host populations, and all sequences were identical. Similarly, a 474-bp region of the wAlbB wsp gene was sequenced from the same populations, and all sequences were identical. These results suggest a role for Wolbachia infection in causing the previously established pattern of low mitochondrial DNA variability, but average nuclear gene diversity, within and among populations of A. albopictus.

Effects of developmental acclimation on adult salinity tolerance in the freshwater-invading copepod Eurytemora affinis

Invasive species are commonly thought to have broad tolerances that enable them to colonize new habitats, but this assumption has rarely been tested. In particular, the relative importance of acclimation (plasticity) and adaptation for invasion success are poorly understood. This study examined effects of short-term and developmental acclimation on adult salinity tolerance in the copepod Eurytemora affinis. This microcrustacean occurs in estuarine and salt marsh habitats but has invaded freshwater habitats within the past century. Effects of short-term acclimation were determined by comparing adult survival in response to acute versus gradual salinity change to low salinity (fresh water). Effects of developmental acclimation on adult tolerance were determined using a split-brood 4 x 2 factorial experimental design for one brackish-water population from Edgartown Great Pond, Massachusetts. Twenty full-sib clutches were split and reared at four salinities (fresh, 5, 10, and 27 practical salinity units [PSU]). On reaching adulthood, clutches from three of the salinity treatments (no survivors at fresh) were split into low- (fresh) and high- (40 PSU) salinity stress treatments, at which survival was measured for 24 h. Short-term acclimation of adults did not appear to have a long-term affect on low-salinity tolerance, given that gradual transfers to fresh water enhanced survival relative to acute transfers in the short term (after 7 h) but not over a longer period of 8 d. Developmental acclimation had contrasting effects on low- versus high-salinity tolerance. Namely, rearing salinity had a significant effect on tolerance of high-salinity (40 PSU) stress but no significant effect on tolerance of low-salinity (freshwater) stress. In addition, there was a significant effect of clutch on survival under freshwater conditions, indicating a genetic component to low-salinity tolerance but no significant clutch effect in response to high salinity. While developmental acclimation might enhance survival at higher salinities, the minimal effect of acclimation and significant effect of clutch on low-salinity tolerance suggest the importance of natural selection during freshwater invasion events.

Genotype-by-environment interaction for salinity tolerance in the freshwater-invading copepod Eurytemora affinis

This study examined the extent of phenotypic plasticity for salinity tolerance and genetic variation in plasticity in the invasive copepod Eurytemora affinis. Euryemora affinis is a species complex inhabiting brackish to hypersaline environments but has invaded freshwater lakes and reservoirs within the past century. Reaction norm experiments were performed on a relatively euryhaline population collected from a brackish lake with fluctuating salinity. Life history traits (hatching rate, survival, and development time) were measured for 20 full-sib clutches that were split and reared at four salinities (fresh, 5, 10, and 27 practical salinity units [PSU]). On average, higher salinities (10 and 27 PSU) were more favorable for larval growth, yielding greater survival and faster development rate. Clutches differed significantly in their response to salinity, with a significant genotype-by-environment interaction for development time. In addition, genetic (clutch) effects were evident in response to low salinity, given that survival in fresh (lake) water was significantly positively correlated with survival at 5 PSU for individual clutches. Clutches raised in fresh water could not survive beyond metamorphosis, suggesting that acclimation to fresh water could not occur in a single generation. Results suggest the importance of natural selection during freshwater invasion events, given the inability of plasticity to generate a freshwater phenotype, and the presence of genetic variation for plasticity upon which natural selection could act.