Eco-immunology of fish invasions: the role of MHC variation

Neutral Forces and Balancing Selection Interplay to Shape the Major Histocompatibility Complex Spatial Patterns in the Striped Hamster in Inner Mongolia: Suggestive of Broad-Scale Local Adaptation

BACKGROUND

The major histocompatibility complex (MHC) plays a key role in the adaptive immune response to pathogens due to its extraordinary polymorphism. However, the spatial patterns of MHC variation in the striped hamster remain unclear, particularly regarding the relative contribution of the balancing selection in shaping MHC spatial variation and diversity compared to neutral forces.

METHODS

In this study, we investigated the immunogenic variation of the striped hamster in four wild populations in Inner Mongolia which experience a heterogeneous parasitic burden. Our goal was to identify local adaptation by comparing the genetic structure at the MHC with that at seven microsatellite loci, taking into account neutral processes.

RESULTS

We observed significant variation in parasite pressure among sites, with parasite burden showing a correlation with temperature and precipitation. Molecular analysis revealed a similar co-structure between MHC and microsatellite loci. We observed lower genetic differentiation at MHC loci compared to microsatellite loci, and no correlation was found between the two.

CONCLUSIONS

Overall, these results suggest a complex interplay between neutral evolutionary forces and balancing selection in shaping the spatial patterns of MHC variation. Local adaptation was not detected on a small scale but may be applicable on a larger scale.

Immunological response of 35 and 42 days old Asian seabass (Lates calcarifer, Bloch 1790) fry following immersion immunization with Streptococcus iniae heat-killed vaccine

Early disease prevention by vaccination requires understanding when fry fish develop specific immunity to a given pathogen. In this research, we explored the immune responses of Asian seabass (Lates calcarifer) at the stages of 35- and 42- days post-hatching (dph) to an immersive heat-killed Streptococcus iniae (Si) vaccine to determine whether fish can produce specific antibodies against the pathogen. The vaccinated fish of each stage (V35 and V42) were immersed with the Si vaccine at 10⁷ CFU/ml for 3 h, whereas the control groups (C35 and C42) were immersed with tryptic soy broth (TSB) in the same manner. Specific antibodies were measured by enzyme-linked immunosorbent assay (ELISA) before and post-immunization (i.e., 0, 7, and 14 days post-immunization, dpi). Expression of innate (TNFα and IL-1β) and adaptive (MHCI, MHCII, CD4, CD8, IgM-like, IgT-like, and IgD-like) immune-related genes were evaluated at the same time points with the addition of 1 dpi. The results showed that a subset of immunized fish from both V35 and V42 fry could elicit specific antibodies (IgM) against Si at 14 dpi. All tested innate and adaptive immune genes upregulated at 7 dpi among fish in V35 group. Interestingly, 42 dph fish appeared to respond to the Si vaccine faster than that of 35 dph, as a significant increase in transcripts was observed in CD4, IL-1β, IgM-like, and IgD-like at 1 dpi; and specific antibody titers of some fish, although not all, were higher than a threshold (p = 0.05) since 7 dpi. In conclusion, this study reveals that 35-42 dph Asian seabass fry can elicit specific immunity to Si immersion vaccine, suggesting that early vaccination of 35 dph fry Asian seabass is feasible.

Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection

Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.

Efficacy testing of an immersion vaccine against Aeromonas salmonicida and immunocompetence in ballan wrasse (Labrus bergylta, Ascanius)

The development of effective vaccines is a critical step towards the domestication of emerging fish species for aquaculture. However, traditional vaccine delivery through intraperitoneal (i.p.) injection requires fish to reach a minimum size and age and therefore cannot provide protection at early developmental stages when infection may occur. This study investigated the effectiveness of immersion vaccination with respect to immunocompetence in a cleaner fish species (ballan wrasse, Labrus bergylta, Ascanius) used in Atlantic salmon farming as an alternative means to control sea lice. The species is susceptible to atypical strains of Aeromonas salmonicida (aAs) at early life stages (<15 g), when i.p. vaccination is not applicable. While immersion vaccination is currently used in commercial hatcheries, the optimal fish size for vaccination, and efficacy of the vaccine delivered by this route has not yet been established. Importantly, efficacy depends on the capability of the species immune system to recognise antigens and process antigens to trigger and produce an adaptive immune response, (process known as immunocompetence). In this study, the efficacy of a polyvalent autogenous vaccine administered by immersion in juvenile ballan wrasse and the subsequent immune response induced was investigated after prime and booster vaccination regimes. In addition, temporal expression (0-150 days post hatch) of adaptive immune genes including major histocompatibility complex (MHC II CD74 molecule) and immunoglobulin M (IgM) was assessed using quantitative PCR (qPCR). Prime and/or boost vaccination by immersion of juvenile ballan wrasse (0.5 g and 1.5 g corresponding to 80 and 170 days post hatch (dph), respectively) did not provide significant protection against aAs vapA V after bath challenge under experimental conditions. Despite no evident protection >80 dph, MHC II and IgM transcripts were first reported at 35 and 75 dph, respectively, suggesting a window of immunocompetence. The results provide important new information on the onset of adaptive immunity in ballan wrasse and highlight that immersion vaccination in the species for protection against aAs should be performed at later developmental stages (>1.5 g) in the hatchery.

SNP analyses and acoustic tagging reveal multiple origins and widespread dispersal of invasive brown trout in the Falkland Islands

Biological invasions are important causes of biodiversity loss, particularly in remote islands. Brown trout (Salmo trutta) have been widely introduced throughout the Southern Hemisphere, impacting endangered native fauna, particularly galaxiid fishes, through predation and competition. However, due to their importance for sport fishing and aquaculture farming, attempts to curtail the impacts of invasive salmonids have generally been met with limited support and the best prospects for protecting native galaxiids is to predict where and how salmonids might disperse. We analysed 266 invasive brown trout from 14 rivers and ponds across the Falkland Islands as well as 32 trout from three potential source populations, using a panel of 592 single nucleotide polymorphisms (SNPs) and acoustic tagging, to ascertain their origins and current patterns of dispersal. We identified four genetically distinct clusters with high levels of genetic diversity and low admixture, likely reflecting the different origins of the invasive brown trout populations. Our analysis suggests that many trout populations in the Falklands may have originated from one of the donor populations analysed (River Wey). The highest genetic diversity was observed in the rivers with the greatest number of introductions and diverse origins, while the lowest diversity corresponded to a location without documented introductions, likely colonized by natural dispersal. High levels of gene flow indicated widespread migration of brown trout across the Falkland Islands, likely aided by anadromous dispersal. This is supported by data from acoustically tagged fish, three of which were detected frequently moving between two rivers ~26 km apart. Our results suggest that, without containment measures, brown trout may invade the last remaining refuges for the native endangered Aplochiton spp. We provide new insights into the origin and dispersal of invasive brown trout in the Falkland Islands that can pave the way for a targeted approach to limit their impact on native fish fauna.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Study on Immune Response of Organs of Epinephelus coioides and Carassius auratus After Immersion Vaccination With Inactivated Vibrio harveyi Vaccine

Immersion vaccination relies on the response of fish mucosa-associated lymphoid tissues, the Crucian carp (Carassius auratus) and Grouper (Epinephelus coioides) were researched in this paper to examine local mucosal immune responses and associated humoral system responses following immersion vaccination. We administered 1.5 × 10⁷ CFU/ml formalin-inactivated Vibrio harveyi cells and measured mucus and serum antibody titers as well as IgM, MHC II mRNA levels in immune organs. The mucosal antibody response preceded the serum response indicating a role for local mucosal immunity in immersion vaccination. IgM and MHC II mRNA levels were relatively greater for the spleen and head kidney indicating the importance and central position of systemic immunity. Expression levels were also high for the gills while skin levels were the lowest. IgM and MHC II mRNA levels were altered over time following vaccination and the hindgut, liver and spleen were similar indicating a close relationship, so the absolute value of r is used to analyze the correlation among different organs immunized. It can be inferred the existence of an internal immune molecular mechanism for Immune synergy hindgut-liver-spleen, from the peak time (14^th day), the relative ratio of genes expression in the same tissues between the immunized grouper and the control group (26 times), and Pearson correlation coefficient (0.8<|r|<1). Injection challenges with live V. harveyi indicated that the relative protection rates for the crucian carp and Grouper was basically the same at 44.4% and 47.4%, respectively. It is believe that crucian carp may be used as a substitute for the valuable grouper in immunity experiment, just from aspect of the relative percent survival (RPS) and how it changes with time. But they were not consistent about the IgM mRNA expression between that of crucian carp and grouper after immersion the Vibrio vaccine.

Comparing raccoon major histocompatibility complex diversity in native and introduced ranges: Evidence for the importance of functional immune diversity for adaptation and survival in novel environments

The adaptive potential of invasive species is related to the genetic diversity of the invader, which is influenced by genetic drift and natural selection. Typically, the genetic diversity of invaders is studied with neutral genetic markers; however, the expectation of reduced diversity has not been consistently supported by empirical studies. Here, we describe and interpret genetic diversity at both neutral microsatellite loci and the immune-related MHC-DRB locus of native and invasive populations of raccoon to better understand of how drift and selection impact patterns of genetic diversity during the invasion process. We found that despite the loss of many MHC (major histocompatibility complex) alleles in comparison with native populations, functional MHC supertypes are preserved in the invasive region. In the native raccoon population, the number of supertypes within individuals was higher than expected under a neutral model. The high level of individual functional divergence may facilitate the adaptation to local conditions in the invasive range. In the invasive populations, we also detected increased population structure at microsatellites compared to the MHC locus, further suggesting that balancing selection is acting on adaptively important regions of the raccoon genome. Finally, we found that alleles known to exhibit resistance to rabies in the native range, Prlo-DRB*4, Prlo-DRB*16 and Prlo-DRB*102, were the most common alleles in the European populations, suggesting directional selection is acting on this locus. Our research shows empirical support for the importance of functional immune diversity for adaptation and survival in novel environments.

Molecular and Structural Characterization of MHC Class II β Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus

This study reports the presence of two distinct MHC class II β genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II β gene.

Contrasting genetic metrics and patterns among naturalized rainbow trout (Oncorhynchus mykiss) in two Patagonian lakes differentially impacted by trout aquaculture

Different pathways of propagation and dispersal of non‐native species into new environments may have contrasting demographic and genetic impacts on established populations. Repeated introductions of rainbow trout (Oncorhynchus mykiss) to Chile in South America, initially through stocking and later through aquaculture escapes, provide a unique setting to contrast these two pathways. Using a panel of single nucleotide polymorphisms, we found contrasting genetic metrics and patterns among naturalized trout in Lake Llanquihue, Chile's largest producer of salmonid smolts for nearly 50 years, and Lake Todos Los Santos (TLS), a reference lake where aquaculture has been prohibited by law. Trout from Lake Llanquihue showed higher genetic diversity, weaker genetic structure, and larger estimates for the effective number of breeders (N_b) than trout from Lake TLS. Trout from Lake TLS were divergent from Lake Llanquihue and showed marked genetic structure and a significant isolation‐by‐distance pattern consistent with secondary contact between documented and undocumented stocking events in opposite shores of the lake. Multiple factors, including differences in propagule pressure, origin of donor populations, lake geomorphology, habitat quality or quantity, and life history, may help explain contrasting genetic metrics and patterns for trout between lakes. We contend that high propagule pressure from aquaculture may not only increase genetic diversity and N_b via demographic effects and admixture, but also may impact the evolution of genetic structure and increase gene flow, consistent with findings from artificially propagated salmonid populations in their native and naturalized ranges.

Effects of invasion history on physiological responses to immune system activation in invasive Australian cane toads

The cane toad (Rhinella marina) has undergone rapid evolution during its invasion of tropical Australia. Toads from invasion front populations (in Western Australia) have been reported to exhibit a stronger baseline phagocytic immune response than do conspecifics from range core populations (in Queensland). To explore this difference, we injected wild-caught toads from both areas with the experimental antigen lipopolysaccharide (LPS, to mimic bacterial infection) and measured whole-blood phagocytosis. Because the hypothalamic-pituitary-adrenal axis is stimulated by infection (and may influence immune responses), we measured glucocorticoid response through urinary corticosterone levels. Relative to injection of a control (phosphate-buffered saline), LPS injection increased both phagocytosis and the proportion of neutrophils in the blood. However, responses were similar in toads from both populations. This null result may reflect the ubiquity of bacterial risks across the toad’s invaded range; utilization of this immune pathway may not have altered during the process of invasion. LPS injection also induced a reduction in urinary corticosterone levels, perhaps as a result of chronic stress.

Genomic patterns of diversity and divergence of two introduced salmonid species in Patagonia, South America

Invasive species have become widespread in aquatic environments throughout the world, yet there are few studies that have examined genomic variation of multiple introduced species in newly colonized environments. In this study, we contrast genomic variation in two salmonid species (anadromous Chinook Salmon, Oncorhynchus tshawytscha, 11,579 SNPs and resident Brook Charr Salvelinus fontinalis, 13,522 SNPs) with differing invasion success after introduction to new environments in South America relative to populations from their native range in North America. Estimates of genetic diversity were not significantly different between introduced and source populations for either species, indicative of propagule pressure that has been shown to maintain diversity in founding populations relative to their native range. Introduced populations also demonstrated higher connectivity and gene flow than those in their native range. Evidence for candidate loci under divergent selection was observed, but was limited to specific introduced populations and was not widely evident. Patterns of genomic variation were consistent with general dispersal potential of each species and therefore also the notion that life history variation may contribute to both invasion success and subsequent genetic structure of these two salmonids in Patagonia.

Temporal and spatial instability in neutral and adaptive (MHC) genetic variation in marginal salmon populations

The role of marginal populations for the long-term maintenance of species’ genetic diversity and evolutionary potential is particularly timely in view of the range shifts caused by climate change. The Centre-Periphery hypothesis predicts that marginal populations should bear reduced genetic diversity and have low evolutionary potential. We analysed temporal stability at neutral microsatellite and adaptive MHC genetic variation over five decades in four marginal Atlantic salmon populations located at the southern limit of the species’ distribution with a complicated demographic history, which includes stocking with foreign and native salmon for at least 2 decades. We found a temporal increase in neutral genetic variation, as well as temporal instability in population structuring, highlighting the importance of temporal analyses in studies that examine the genetic diversity of peripheral populations at the margins of the species’ range, particularly in face of climate change.

Variation in Major Histocompatibility Complex diversity in invasive cane toad populations

Context The cane toad (Rhinella marina), a native species of central and southern America, was introduced to Australia in 1935 as a biocontrol agent after a complex history of prior introductions. The population rapidly expanded and has since spread through much of the Australian landmass, with severe impacts on the endemic wildlife, primarily via toxicity to predators. The invasion process has taken its toll on the cane toad, with changes in the immunological capacity across the Australian invasive population. Aims To investigate the immunogenetic underpinnings of these changes, we studied the diversity of the Major Histocompatiblity Complex (MHC) genes in introduced cane toad populations. Methods We studied the diversity of two MHC genes (the classical class I UA locus and a class II DAB locus) and compared these with neutral microsatellite markers in toads from the Australian site of introduction and the Australian invasion front. We also included toads from Hawai’i, the original source of the Australian toads, to infer founder effect. Key results Diversity across all markers was low across Australian and Hawai’ian samples, consistent with a reduction in genetic diversity through multiple founder effects during the course of the successive translocations. In Australia, allelic diversity at the microsatellite markers and the UA locus was reduced at the invasion front, whereas all three alleles at the DAB locus were maintained in the invasion-front toads. Conclusions Loss of allelic diversity observed at the microsatellite markers and the UA locus could be the result of drift and bottlenecking along the invasion process, however, the persistence of DAB diversity warrants further investigation to disentangle the evolutionary forces influencing this locus. Implications Through the use of different molecular markers, we provide a preliminary description of the adaptive genetic processes occurring in this invasive population. The extremely limited MHC diversity may represent low immunogenetic competence across the Australian population, which could be exploited for invasive species management.

Eco-immunology and bioinvasion: revisiting the evolution of increased competitive ability hypotheses

Immunity is at the core of major theories related to invasion biology. Among them, the evolution of increased competitive ability (EICA) and EICA‐refined hypotheses have been used as a reference work. They postulate that the release from pathogens often experienced during invasion should favour a reallocation of resources from (costly) immune defences to beneficial life‐history traits associated with invasive potential. We review studies documenting immune changes during animal invasions. We describe the designs and approaches that have been applied and discuss some reasons that prevent drawing generalized conclusions regarding EICA hypotheses. We detail why a better assessment of invasion history and immune costs, including immunopathologies and parasite communities, could improve our understanding of the relationships between immunity and invasion success. Finally, we propose new perspectives to revisit the EICA hypotheses. We first emphasize the neutral and adaptive mechanisms involved in immune changes, as well as timing of the later. Such investigation will help decipher whether immune changes are a consequence of pre‐adaptation, or the result of postintroduction adaptations to invasion front conditions. We next bring attention to new avenues of research that remain unexplored, namely age‐dependent immunity and gut microbiota, potential key factors underlying adaptation to invasion front environment and modulating invasion success.

Gyrodactylus spp. diversity in native and introduced minnow (Phoxinus phoxinus) populations: no support for "the enemy release" hypothesis

BACKGROUND

Translocation of native species and introduction of non-native species are potentially harmful to the existing biota by introducing e.g. diseases, parasites and organisms that may negatively affect the native species. The enemy release hypothesis states that parasite species will be lost from host populations when the host is introduced into new environments.

METHODS

We tested the enemy release hypothesis by comparing 14 native and 29 introduced minnow (Phoxinus phoxinus) populations in Norway with regard to the ectoparasitic Gyrodactylus species community and load (on caudal fin). Here, we used a nominal logistic regression on presence/absence of Gyrodactylus spp. and a generalized linear model on the summed number of Gyrodactylus spp. on infected populations, with individual minnow heterozygosity (based on 11 microsatellites) as a covariate. In addition, a sample-based rarefaction analysis was used to test if the Gyrodactylus-species specific load differed between native and introduced minnow populations. An analysis of molecular variance was performed to test for hierarchical population structure between the two groups and to test for signals of population bottlenecks the two-phase model in the Wilcoxon signed-rank test was used. To test for demographic population expansion events in the introduced minnow population, we used the kg-test under a stepwise mutation model.

RESULTS

The native and introduced minnow populations had similar species compositions of Gyrodactylus, lending no support to the enemy release hypothesis. The two minnow groups did not differ in the likelihood of being infected with Gyrodactylus spp. Considering only infected minnow populations it was evident that native populations had a significantly higher mean abundance of Gyrodactylus spp. than introduced populations. The results showed that homozygotic minnows had a higher Gyrodactylus spp. infection than more heterozygotic hosts. Using only infected individuals, the two minnow groups did not differ in their mean number of Gyrodactylus spp. However, a similar negative association between heterozygosity and abundance was observed in the native and introduced group. There was no evidence for demographic bottlenecks in the minnow populations, implying that introduced populations retained a high degree of genetic variation, indicating that the number of introduced minnows may have been large or that introductions have been happening repeatedly. This could partly explain the similar species composition of Gyrodactylus in the native and introduced minnow populations.

CONCLUSIONS

In this study it was observed that native and introduced minnow populations did not differ in their species community of Gyrodactylus spp., lending no support to the enemy release hypothesis. A negative association between individual minnow host heterozygosity and the number of Gyrodactylus spp. was detected. Our results suggest that the enemy release hypothesis does not necessarily limit fish parasite dispersal, further emphasizing the importance of invasive fish species dispersal control.

Temporal Genetic Variance and Propagule-Driven Genetic Structure Characterize Naturalized Rainbow Trout (Oncorhynchus mykiss) from a Patagonian Lake Impacted by Trout Farming

Knowledge about the genetic underpinnings of invasions—a theme addressed by invasion genetics as a discipline—is still scarce amid well documented ecological impacts of non-native species on ecosystems of Patagonia in South America. One of the most invasive species in Patagonia’s freshwater systems and elsewhere is rainbow trout (Oncorhynchus mykiss). This species was introduced to Chile during the early twentieth century for stocking and promoting recreational fishing; during the late twentieth century was reintroduced for farming purposes and is now naturalized. We used population- and individual-based inference from single nucleotide polymorphisms (SNPs) to illuminate three objectives related to the establishment and naturalization of Rainbow Trout in Lake Llanquihue. This lake has been intensively used for trout farming during the last three decades. Our results emanate from samples collected from five inlet streams over two seasons, winter and spring. First, we found that significant intra- population (temporal) genetic variance was greater than inter-population (spatial) genetic variance, downplaying the importance of spatial divergence during the process of naturalization. Allele frequency differences between cohorts, consistent with variation in fish length between spring and winter collections, might explain temporal genetic differences. Second, individual-based Bayesian clustering suggested that genetic structure within Lake Llanquihue was largely driven by putative farm propagules found at one single stream during spring, but not in winter. This suggests that farm broodstock might migrate upstream to breed during spring at that particular stream. It is unclear whether interbreeding has occurred between “pure” naturalized and farm trout in this and other streams. Third, estimates of the annual number of breeders (N_b) were below 73 in half of the collections, suggestive of genetically small and recently founded populations that might experience substantial genetic drift. Our results reinforce the notion that naturalized trout originated recently from a small yet genetically diverse source and that farm propagules might have played a significant role in the invasion of Rainbow Trout within a single lake with intensive trout farming. Our results also argue for proficient mitigation measures that include management of escapes and strategies to minimize unintentional releases from farm facilities.