Limited hatchery introgression into wild brook trout (Salvelinus fontinalis) populations despite reoccurring stocking

Human Impacts on Great Lakes Walleye Sander vitreus Structure, Diversity and Local Adaptation

Artificial propagation and wild release may influence the genetic integrity of wild populations. This practice has been prevalent in fisheries for centuries and is often termed 'stocking'. In the Laurentian Great Lakes (Great Lakes here-on), walleye populations faced declines from the 1950s to the 1970s, prompting extensive stocking efforts for restoration. By the mid-2010s, walleye populations showed signs of recovery, but the genetic legacy of stocking on population structure at the genomic level remains unclear. Using a dataset of 45,600 genome-aligned SNP loci genotyped in 1075 walleye individuals, we investigated the genetic impacts of over 50 years of stocking across the Great Lakes. Population structure was associated with both natural geographic barriers and stocking from non-native sources. Admixture between Lake Erie walleye and walleye from the re-populated Tittabawassee River indicate that stocking may have re-distributed putatively adaptive alleles around the Great Lakes. Genome scans identified FST outliers and evidence of selective sweeps, indicating local adaptation of spawning populations is likely. Notably, one genomic region showed strong differentiation between Muskegon River and walleye from the Tittabawassee River, which was re-populated by Muskegon strain walleye, suggesting admixture and selection both impact the observed genetic diversity. Overall, our study underscores how artificial propagation and translocations can significantly alter the evolutionary trajectory of populations. The findings highlight the complex interplay between stocking practices and population genetic diversity, emphasising the need for careful management strategies to preserve the genetic integrity of wild populations amidst conservation efforts.

Evaluation of genetic consequences of stocking on the southern-margin populations of white-spotted charr

Coldwater-adapted freshwater fishes, especially their populations along warm-range margins, are most vulnerable to the climate oscillations associated with global warming. Stocking is a major strategy for avoiding the extinction of these species. However, while stocking can reverse the decline of isolated populations, it may also result in a loss of genetic diversity in the native local population due to the introgressive replacement of hatchery genes. To plan an adequate strategy for conserving locally adapted populations, the genetic impacts of stocking on native lineages should be evaluated from small river branches to wide-ranging drainage areas. We investigated the population genetic structure of white-spotted charr (Salvelinus leucomaenis) within its southern range (Lake Biwa basin, Japan). By applying genome-wide SNP analysis to the population's genetic structure, we assessed the extent of genetic introgression resulting from stocking. White-spotted charr in the Lake Biwa watershed constitutes a distinctive genetic group, within which apparent genetic differentiation was observed. The hatchery-reared fish line commonly used for supplementation stocking in the catchment was discernable from the native population, enabling us to analyze genetic introgression across the entire drainage area. Admixed individuals resulting from hatchery introgression were observed in most of the stocked sites that showed relatively high heterozygosity and nucleotide diversity. However, their genetic differentiation was much lower than that of native populations. The supplementation history as well as the road availability contributed substantially to the introgression of hatchery genes. Populations with the native genetic structure remained in the upstream regions of the tested rivers. However, their heterozygosity and nucleotide diversity were low when compared with that of the populations with hatchery supplementation. Our results shed light on the genetic impacts of stocking on isolated native populations and suggest that conventional supplementation methods cannot preserve a unique biodiversity in the distribution margin.

Competitive interactions affect introgression and population viability amidst maladaptive hybridization

The deliberate release of captive-bred individuals, the accidental escape of domesticated strains, or the invasion of closely related conspecifics into wild populations can all lead to introgressive hybridization, which poses a challenge for conservation and wildlife management. Rates of introgression and the magnitude of associated demographic impacts vary widely across ecological contexts. However, the reasons for this variation remain poorly understood. One rarely considered phenomenon in this context is soft selection, wherein relative trait values determine success in intraspecific competition for a limiting resource. Here we develop an eco-genetic model explicitly focussed on understanding the influence of such competitive interactions on the eco-evolutionary dynamics of wild populations experiencing an influx of foreign/domesticated individuals. The model is applicable to any taxon that experiences natural or human-mediated inputs of locally maladapted genotypes ('intrusion'), in addition to phenotype-dependent competition for a limiting resource (e.g. breeding sites, feeding territories). The effects of both acute and chronic intrusion depended strongly on the relative competitiveness of intruders versus locals. When intruders were competitively inferior, density-dependent regulation limited their reproductive success (ability to compete for limited spawning sites), which prevented strong introgression or population declines from occurring. In contrast, when intruders were competitively superior, this amplified introgression and led to increased maladaptation of the admixed population. This had negative consequences for population size and population viability. The results were sensitive to the intrusion level, the magnitude of reproductive excess, trait heritability and the extent to which intruders were maladapted relative to locals. Our findings draw attention to under-appreciated interactions between phenotype-dependent competitive interactions and maladaptive hybridization, which may be critical to determining the impact captive breeding programmes and domesticated escapees can have on otherwise self-sustaining wild populations.

Transgenerational effects on body size and survival in Brook charr (Salvelinus fontinalis)

Higher temperatures are now observed in several ecosystems and act as new selective agents that shape traits and fitness of individuals. Transgenerational effects may be important in modulating adaptation of future generations and buffering negative impacts of temperature changes. The potential for these effects may be important in freshwater fish species, as temperature is a key abiotic component of their environment. Yet, still, relatively few studies have assessed the presence and importance of transgenerational effects under natural conditions. The purpose of this study was to test how parental thermal conditions influenced offspring growth and survival following stocking in Brook charr (Salvelinus fontinalis). To do so, part of the breeders were exposed to a "cold" treatment while others were exposed to a "warm" treatment during the final steps of gonad maturation (constant 2°C difference between treatments along the seasonal temperature decrease). The impact on offspring of a selection treatment targeting production traits of interest (absence of sexual maturation at 1+, combined with increased growth) in breeders was also evaluated. After 7-8 months of growth in captivity, offspring were stocked in natural lakes. Their growth and survival were assessed about a year later. Offspring from "cold" breeders showed lower survival than those from "warm" breeders and the selection treatment had no effect on survival. However, the selection treatment was linked to lower Fulton's condition index, which, in turn, was positively correlated to survival in lakes. This study highlights the importance of working in ecological/industrial context to fully assess the different impacts of transgenerational effects on traits and survival. Our results also have important implications for stocking practices used to support the sport fishing industry.

Genomics facilitates evaluation and monitoring of McCloud River Redband Trout (Oncorhynchus mykiss stonei)

The McCloud River Redband Trout (MRRT; Oncorhynchus mykiss stonei) is a unique subspecies of rainbow trout that inhabits the isolated Upper McCloud River of Northern California. A major threat to MRRT is introgressive hybridization with non-native rainbow trout from historical stocking and contemporary unauthorized introductions. To help address this concern, we collected RAD-sequencing data on 308 total individuals from MRRT and other California O. mykiss populations and examined population structure using Principal Component and admixture analyses. Our results are consistent with previous studies; we found that populations of MRRT in Sheepheaven, Swamp, Edson, and Moosehead creeks are nonintrogressed. Additionally, we saw no evidence of introgression in Dry Creek, and suggest further investigation to determine if it can be considered a core MRRT conservation population. Sheepheaven Creek was previously thought to be the sole historical lineage of MRRT, but our analysis identified three: Sheepheaven, Edson, and Dry creeks, all of which should be preserved. Finally, we discovered diagnostic and polymorphic SNP markers for monitoring introgression and genetic diversity in MRRT. Collectively, our results provide a valuable resource for the conservation and management of MRRT.

Conservation Genetics of Mediterranean Brown Trout in Central Italy (Latium): A Multi-Marker Approach

Brown trout is considered a complex of incipient species, including several phylogenetic lineages, whose natural distribution in the Mediterranean area has been altered, since the beginning of the 1900s, by massive introductions of domestic strains of Atlantic origin to support fisheries. Introduced trout naturalize in new suitable environments and extensively hybridize with native populations. Here, we characterized putatively neutral and adaptive genetic variability and population structure of Mediterranean brown trout from six river catchments in central peninsular Italy, as revealed by both mitochondrial (Control Region) and nuclear (microsatellites, LDH-C1, major histocompatibility complex) markers. We quantified the admixture of wild populations with hatchery strains and evaluated the effects of domestic trout introductions on shaping population genetics. Our analyses indicated: (1) a composite picture of genetic variability in the area, with the presence of all native Mediterranean trout mitochondrial lineages (“Adriatic”, “Mediterranean”, “marmoratus”), various frequencies of allochthonous genotypes and different rates of introgression among sampling sites; (2) asymmetric mito-nuclear introgression; (3) increasing nuclear marker diversity with increasing levels of admixture across populations; (4) strong population structure coupled with relatively low effective population size. Data allowed the identification of five management units and we propose specific actions to support ongoing and future conservation strategies within the examined area.

Temporal analysis shows relaxed genetic erosion following improved stocking practices in a subarctic transnational brown trout population

Maintaining standing genetic variation is a challenge in human-dominated landscapes. We used genetic (i.e., 16 short tandem repeats) and morphological (i.e., length and weight) measurements of 593 contemporary and historical brown trout (Salmo trutta) samples to study fine-scale and short-term impacts of different management practices. These had changed from traditional breeding practices, using the same broodstock for several years, to modern breeding practices, including annual broodstock replacement, in the transnational subarctic Pasvik River. Using population genetic structure analyses (i.e., Bayesian assignment tests, DAPCs, and PCAs), four historical genetic clusters (E2001A-D), likely representing family lineages resulting from different crosses, were found in zone E. These groups were characterized by consistently lower genetic diversity, higher within-group relatedness, lower effective population size, and significantly smaller body size than contemporary stocked (E2001E) and wild fish (E2001F). However, even current breeding practices are insufficient to prevent genetic diversity loss and morphological changes as demonstrated by on average smaller body sizes and recent genetic bottleneck signatures in the modern breeding stock compared to wild fish. Conservation management must evaluate breeding protocols for stocking programs and assess if these can preserve remaining natural genetic diversity and morphology in brown trout for long-term preservation of freshwater fauna.

Associative Overdominance and Negative Epistasis Shape Genome-Wide Ancestry Landscape in Supplemented Fish Populations

The interplay between recombination rate, genetic drift and selection modulates variation in genome-wide ancestry. Understanding the selective processes at play is of prime importance toward predicting potential beneficial or negative effects of supplementation with domestic strains (i.e., human-introduced strains). In a system of lacustrine populations supplemented with a single domestic strain, we documented how population genetic diversity and stocking intensity produced lake-specific patterns of domestic ancestry by taking the species’ local recombination rate into consideration. We used 552 Brook Charr (Salvelinus fontinalis) from 22 small lacustrine populations, genotyped at ~32,400 mapped SNPs. We observed highly variable patterns of domestic ancestry between each of the 22 populations without any consistency in introgression patterns of the domestic ancestry. Our results suggest that such lake-specific ancestry patterns were mainly due to variable associative overdominance (AOD) effects among populations (i.e., potential positive effects due to the masking of possible deleterious alleles in low recombining regions). Signatures of AOD effects were also emphasized by highly variable patterns of genetic diversity among and within lakes, potentially driven by predominant genetic drift in those small isolated populations. Local negative effects such as negative epistasis (i.e., potential genetic incompatibilities between the native and the introduced population) potentially reflecting precursory signs of outbreeding depression were also observed at a chromosomal scale. Consequently, in order to improve conservation practices and management strategies, it became necessary to assess the consequences of supplementation at the population level by taking into account both genetic diversity and stocking intensity when available.

A novel quantitative framework for riverscape genetics

Riverscape genetics, which applies concepts in landscape genetics to riverine ecosystems, lack appropriate quantitative methods that address the spatial autocorrelation structure of linear stream networks and account for bidirectional geneflow. To address these challenges, we present a general framework for the design and analysis of riverscape genetic studies. Our framework starts with the estimation of pairwise genetic distance at sample sites and the development of a spatially structured ecological network (SSEN) on which riverscape covariates are measured. We then introduce the novel bidirectional geneflow in riverscapes (BGR) model that uses principles of isolation-by-resistance to quantify the effects of environmental covariates on genetic connectivity, with spatial covariance defined using simultaneous autoregressive models on the SSEN and the generalized Wishart distribution to model pairwise distance matrices arising through a random walk model of geneflow. We highlight the utility of this framework in an analysis of riverscape genetics for brook trout (Salvelinus fontinalis) in north central Pennsylvania, USA. Using the fixation index (F_ST ) as the measure of genetic distance, we estimated the effects of 12 riverscape covariates on geneflow by evaluating the relative support of eight competing BGR models. We then compared the performance of the top-ranked BGR model to results obtained from comparable analyses using multiple regression on distance matrices (MRM) and the program STRUCTURE. We found that the BGR model had more power to detect covariate effects, particularly for variables that were only partial barriers to geneflow and/or uncommon in the riverscape, making it more informative for assessing patterns of population connectivity and identifying threats to species conservation. This case study highlights the utility of our modeling framework over other quantitative methods in riverscape genetics, particularly the ability to rigorously test hypotheses about factors that influence geneflow and probabilistically estimate the effect of riverscape covariates, including stream flow direction. This framework is flexible across taxa and riverine networks, is easily executable, and provides intuitive results that can be used to investigate the likely outcomes of current and future management scenarios.

Multiple decades of stocking has resulted in limited hatchery introgression in wild brook trout (Salvelinus fontinalis) populations of Nova Scotia

Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long-term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation-by-distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery-derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long-term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case-by-case basis.

Can genetic assignment tests provide insight on the influence of captive egression on the epizootiology of chronic wasting disease?

Identifying the sources of ongoing and novel disease outbreaks is critical for understanding the diffusion of epizootic diseases. Identifying infection sources is difficult when few physical differences separate individuals with different origins. Genetic assignment procedures show great promise for assessing transmission dynamics in such situations. Here, we use genetic assignment tests to determine the source of chronic wasting disease infections in free-ranging white-tailed deer (Odocoileus virginianus) populations. Natural dispersal is thought to facilitate the geographic diffusion of chronic wasting disease, but egression from captive cervid populations represents an alternative source of infection that is difficult to detect due to physical similarities with wild deer. Simulated reference populations were created based on allele frequencies from 1,912 empirical microsatellite genotypes collected in four sampling subregions and five captive facilities. These reference populations were used to assess the likelihood of ancestry and assignment of 1,861 free-ranging deer (1,834 noninfected and 27 infected) and 51 captive individuals to captive or wild populations. The ancestry (Q) and assignment scores (A) for free-ranging deer to wild populations were high (average Q wild = 0.913 and average A wild = 0.951, respectively), but varied among subregions (Q wild = 0.800-0.947, A wild = 0.857-0.976). These findings suggest that captive egression and admixture are rare, but risk may not be spatially uniform. Ancestry and assignment scores for two free-ranging deer with chronic wasting disease sampled in an area where chronic wasting disease was previously unobserved in free-ranging herds indicated a higher likelihood of assignment and proportion of ancestry attributable to captive populations. While we cannot directly assign these individuals to infected facilities, these findings suggest that rare egression events may influence the epizootiology of chronic wasting disease in free-ranging populations. Continued disease surveillance and genetic analyses may further elucidate the relative disease risk attributable to captive and wild sources.

Genetic diversity, admixture, and hatchery influence in Brook Trout (Salvelinus fontinalis) throughout western New York State

Although Brook Trout are distributed across most of eastern North America, population numbers have declined in many regions due to habitat loss, climate change, and competition with non-native species. In New York State, Brook Trout habitat has been substantially reduced, with many areas showing complete extirpation of Brook Trout populations, predominantly in the western portion of the state. Small, fragmented populations are at risk of genetic diversity loss, inbreeding depression, and reduced fitness, leading to a greater potential for local extirpation. Genetic monitoring is a practical tool that can facilitate further conservation-decision making regarding small populations. In this study, we used 12 microsatellite loci to examine 3,436 sampled Brook Trout, representing 75 sites from the Allegheny, Erie/Niagara, Genesee, Oswego, Lake Ontario, and Susquehanna drainage basins throughout western New York State. Three Brook Trout hatchery strains were also genetically characterized to evaluate the degree of hatchery introgression between wild populations and hatchery strains stocked in the region. Overall, estimates of genetic diversity varied widely: Allelic richness ranged from 2.23 to 7.485, and expected heterozygosity ranged from 0.402 to 0.766. As observed for Brook Trout in other regions, we found a high degree of genetic differentiation among populations, with all comparisons except one showing significant F ST values. Hatchery introgression was found to be minimal, with estimates ranging from 1.96% to 3.10% of wild individuals exhibiting membership proportions to a hatchery strain cluster exceeding 10% (q ≥ 0.10). Results from this investigation can be used to prioritize management efforts for Brook Trout in western New York State and act as a baseline to monitor future population trends.