Understanding admixture patterns in supplemented populations: a case study combining molecular analyses and temporally explicit simulations in Atlantic salmon

Assessing the impact of pedigree attributes on the validity of quantitative genetic parameter estimates

Investigating the evolution of complex traits in nature requires accurate assessment of their genetic basis. Quantitative genetic (QG) modeling is frequently applied to estimate the additive genetic variance (VA) in traits, combining phenotypic and pedigree data from a sample of individuals. Whether reconstructed from social links or molecular markers, empirical pedigrees differ in completeness, genealogical error rates, and other attributes that can impact QG estimation. Here we investigate this impact using human genealogical data for 6 French-Canadian (FC) populations originating from the same genetic founding source but differing in their pedigrees' attributes. First, we simulated phenotypic values along pedigrees and under different trait architectures and "true" parameter values (e.g., VA). Then we fitted mixed effects "animal" models to these simulated data, to assess how QG estimation was impacted by pedigree attributes. Our results show that pedigree size and depth were important determinants of the precision, but not accuracy, of genetic parameter estimates. In contrast, pedigree completeness and entropy, 2 attributes related to the density of genealogical links, were not clearly associated with the performance of parameter estimation. Noticeably, a slight increase in the genealogical error rate was sufficient to cause a detectable underestimation of VA. Including maternal genetic effects in the simulations led to a slight underestimation of VA with pedigrees of smaller size and depth. Despite originating from the same genetic source, the 6 pedigrees yielded wide variations in QG estimates under identical conditions. These findings highlight the importance of sensitivity analyses in pedigree-based genetic studies on natural populations.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Captive-bred Atlantic salmon released into the wild have fewer offspring than wild-bred fish and decrease population productivity

The release of captive-bred animals into the wild is commonly practised to restore or supplement wild populations but comes with a suite of ecological and genetic consequences. Vast numbers of hatchery-reared fish are released annually, ostensibly to restore/enhance wild populations or provide greater angling returns. While previous studies have shown that captive-bred fish perform poorly in the wild relative to wild-bred conspecifics, few have measured individual lifetime reproductive success (LRS) and how this affects population productivity. Here, we analyse data on Atlantic salmon from an intensely studied catchment into which varying numbers of captive-bred fish have escaped/been released and potentially bred over several decades. Using a molecular pedigree, we demonstrate that, on average, the LRS of captive-bred individuals was only 36% that of wild-bred individuals. A significant LRS difference remained after excluding individuals that left no surviving offspring, some of which might have simply failed to spawn, consistent with transgenerational effects on offspring survival. The annual productivity of the mixed population (wild-bred plus captive-bred) was lower in years where captive-bred fish comprised a greater fraction of potential spawners. These results bolster previous empirical and theoretical findings that intentional stocking, or non-intentional escapees, threaten, rather than enhance, recipient natural populations.

Long-term stocking practices threaten the original genetic diversity of the southernmost European populations of Atlantic salmon Salmo salar

Many Atlantic salmon Salmo salar populations in Europe are threatened by previous stocking with foreign hatchery strains. Temporal patterns of genetic characteristics of salmon from northern Spain, the southernmost European populations, were compared before and after species decline and heavy stocking with specimens from northern Europe. Eleven microsatellite loci were analysed in archival (scales from 1958-1960) and contemporary (2007-2008) samples from the River Sella. Temporal analyses revealed a similar heterozygosity between archival and contemporary samples, despite a drastic decrease in population abundance, while the contemporary sample showed a higher allelic richness due to the occurrence of foreign alleles. Considering only the alleles with at least 4% frequency in the archival sample, 2 alleles exclusive to the River Sella were absent in the contemporary sample, and 14 alleles showed a decrease of at least 4% frequency. Four alleles common in Scotland showed a high occurrence in the contemporary sample, so they are good candidates as markers of introgression of foreign genes. The heavy stocking with non-native Scottish broodstocks between 1970 and 1990 caused the introgression found in the contemporary sample when compared with the pristine population. An abrupt decrease was evident when the estimates of effective number of breeders were adjusted to take into account overlapping generations (NbAdj), effective population size (NeAdj) estimated from NbAdj, and number of breeders estimated using the sibship assignment method (NbSIB). The very low effective size values found in the contemporary sample, together with the detrimental synergy between genetic drift and high rates of introgression, represent a severe risk for the conservation of native salmon.

Combining population genomics and forward simulations to investigate stocking impacts: A case study of Muskellunge (Esox masquinongy) from the St. Lawrence River basin

Understanding the genetic and evolutionary impacts of stocking on wild fish populations has long been of interest as negative consequences such as reduced fitness and loss of genetic diversity are commonly reported outcomes. In an attempt to sustain a fishery, managers implemented nearly five decades of extensive stocking of over a million Muskellunge (Esox masquinongy), a native species in the Lower St. Lawrence River (Québec, Canada). We investigated the effect of this stocking on population genetic structure and allelic diversity in the St. Lawrence River in addition to tributaries and several stocked inland lakes. Using genotype by sequencing, we genotyped 643 individuals representing 22 locations and combined this information with forward simulations to investigate the genetic consequences of long-term stocking. Individuals native to the St. Lawrence watershed were genetically differentiated from stocking sources and tributaries, and inland lakes were naturally differentiated from the main river. Empirical data and simulations within the St. Lawrence River revealed weak stocking effects on admixture patterns. Our data suggest that the genetic structure associated with stocked fish was diluted into its relatively large effective population size. This interpretation is also consistent with a hypothesis that selection against introgression was in operation and relatively efficient within the large St. Lawrence River system. In contrast, smaller populations from adjacent tributaries and lakes displayed greater stocking-related admixture that resulted in comparatively higher heterozygosity than the St. Lawrence. Finally, individuals from inland lakes that were established by stocking maintained a close affinity with their source populations. This study illustrated a benefit of combining extensive genomic data with forward simulations for improved inference regarding population-level genetic effects of long-term stocking, and its relevance for fishery management decision making.

Chromosomal fusion and life history-associated genomic variation contribute to within-river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST  = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST  = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250-Kb region of chromosome 9, including the gene six6 that has previously been linked to age-at-maturity and run-timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.

The demographic history of Atlantic salmon (Salmo salar) across its distribution range reconstructed from approximate Bayesian computations

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo-genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.

Predicting the genetic impact of stocking in Brook Charr (Salvelinus fontinalis) by combining RAD sequencing and modeling of explanatory variables

In fisheries management, intensive stocking programs are commonly used to enhance population abundance and maintain stock productivity. However, such practices are increasingly raising concerns as multiple studies documented adverse genetic and evolutionary impacts of stocking on wild populations. Improvement of stocking management relies on a better understanding of the dynamic of introgressive hybridization between wild and domestic population and on assessment of the genetic state of wild populations after stocking cessation. In Québec, Canada, over five million captive-reared Brook Charr (Salvelinus fontinalis) are stocked every year to support recreational fishing activities. Here, we investigated how variation in stocking history and environmental variables, including water temperature, pH, and dissolved oxygen, may influence the impact of stocking practices on the genetic integrity of wild Brook Charr populations. We collected DNA samples (n = 862, average of 30 individuals per lake) from 29 lakes that underwent different stocking intensity through time and also collected environmental parameters for each sampled lake. An average of 4,580 high-quality filtered SNPs was obtained for each population using genotyping by sequencing (GBS), which were then used to quantify the mean domestic membership of each sampled population. An exhaustive process of model selection was conducted to obtain a best-fitted model that explained 56% of the variance observed in mean domestic genetic membership. The number of years since the mean year of stocking was the best explanatory variable to predict variation in mean domestic genetic membership whereas environmental characteristics had little influence on observed patterns of admixture. Our model predictions also revealed that each sampled wild population could potentially return to a wild genetic state (absence of domestic genetic background) after stocking cessation. Overall, our study provides new insights on factors determining level of introgressive hybridization and suggests that stocking impacts could be reversible with time.

Atlantic salmon Salmo salar in the chalk streams of England are genetically unique

Recent research has identified genetic groups of Atlantic salmon Salmo salar that show association with geological and environmental boundaries. This study focuses on one particular subgroup of the species inhabiting the chalk streams of southern England, U.K. These fish are genetically distinct from other British and European S. salar populations and have previously demonstrated markedly low admixture with populations in neighbouring regions. The genetic population structure of S. salar occupying five chalk streams was explored using 16 microsatellite loci. The analysis provides evidence of the genetic distinctiveness of chalk-stream S. salar in southern England, in comparison with populations from non-chalk regions elsewhere in western Europe. Little genetic differentiation exists between the chalk-stream populations and a pattern of isolation by distance was evident. Furthermore, evidence of temporal stability of S. salar populations across the five chalk streams was found. This work provides new insights into the temporal stability and lack of genetic population sub-structuring within a unique component of the species' range of S. salar.

Loss of genetic integrity and biological invasions result from stocking and introductions of Barbus barbus: insights from rivers in England

Anthropogenic activities, including the intentional releases of fish for enhancing populations (stocking), are recognized as adversely impacting the adaptive potential of wild populations. Here, the genetic characteristics of European barbel Barbus barbus were investigated using 18 populations in England, where it is indigenous to eastern‐flowing rivers and where stocking has been used to enhance these populations. Invasive populations are also present in western‐flowing rivers following introductions of translocated fish. Two genetic clusters were evident in the indigenous range, centered on catchments in northeast and southeast England. However, stocking activities, including the release of hatchery‐reared fish, have significantly reduced the genetic differentiation across the majority of this range. In addition, in smaller indigenous rivers, populations appeared to mainly comprise fish of hatchery origin. In the nonindigenous range, genetic data largely aligned to historical stocking records, corroborating information that one particular river (Kennet) in southeast England was the original source of most invasive B. barbus in England. It is recommended that these genetic outputs inform management measures to either restore or maintain the original genetic diversity of the indigenous rivers, as this should help ensure populations can maintain their ability to adapt to changing environmental conditions. Where stocking is considered necessary, it is recommended that only broodstock from within the catchment is used.

Genetic and phenotypic changes in an Atlantic salmon population supplemented with non-local individuals: a longitudinal study over 21 years

While introductions and supplementations using non-native and potentially domesticated individuals may have dramatic evolutionary effects on wild populations, few studies documented the evolution of genetic diversity and life-history traits in supplemented populations. Here, we investigated year-to-year changes from 1989 to 2009 in genetic admixture at 15 microsatellite loci and in phenotypic traits in an Atlantic salmon (Salmo salar) population stocked during the first decade of this period with two genetically and phenotypically distinct source populations. We detected a pattern of temporally increasing introgressive hybridization between the stocked population and both source populations. The proportion of fish returning to the river after a single winter at sea (versus several ones) was higher in fish assigned to the main source population than in local individuals. Moreover, during the first decade of the study, both single-sea-winter and multi-sea-winter (MSW) fish assigned to the main source population were smaller than local fish. During the second decade of the study, MSW fish defined as hybrids were lighter and smaller than fish from parental populations, suggesting outbreeding depression. Overall, this study suggests that supplementation with non-local individuals may alter not only the genetic diversity of wild populations but also life-history traits of adaptive significance.

Loss of genetic integrity in wild lake trout populations following stocking: insights from an exhaustive study of 72 lakes from Québec, Canada

Stocking represents the most important management tool worldwide to increase and sustain commercial and recreational fisheries in a context of overexploitation. Genetic impacts of this practice have been investigated in many studies, which examined population and individual admixture, but few have investigated determinants of these processes. Here, we addressed these questions from the genotyping at 19 microsatellite loci of 3341 adult lake trout (Salvelinus namaycush) from 72 unstocked and stocked lakes. Results showed an increase in genetic diversity and a twofold decrease in the extent of genetic differentiation among stocked populations when compared to unstocked. Stocked populations were characterized by significant admixture at both population and individual levels. Moreover, levels of admixture in stocked populations were strongly correlated with stocking intensity and a threshold value of total homogenization between source and stocked populations was identified. Our results also suggest that under certain scenarios, the genetic impacts of stocking could be of short duration. Overall, our study emphasizes the important alteration of the genetic integrity of stocked populations and the need to better understand determinants of admixture to optimize stocking strategies and to conserve the genetic integrity of wild populations.

Source-sink estimates of genetic introgression show influence of hatchery strays on wild chum salmon populations in Prince William Sound, Alaska

The extent to which stray, hatchery-reared salmon affect wild populations is much debated. Although experiments show that artificial breeding and culture influence the genetics of hatchery salmon, little is known about the interaction between hatchery and wild salmon in a natural setting. Here, we estimated historical and contemporary genetic population structures of chum salmon (Oncorhynchus keta) in Prince William Sound (PWS), Alaska, with 135 single nucleotide polymorphism (SNP) markers. Historical population structure was inferred from the analysis of DNA from fish scales, which had been archived since the late 1960’s for several populations in PWS. Parallel analyses with microsatellites and a test based on Hardy-Weinberg proportions showed that about 50% of the fish-scale DNA was cross-contaminated with DNA from other fish. These samples were removed from the analysis. We used a novel application of the classical source-sink model to compare SNP allele frequencies in these archived fish-scales (1964–1982) with frequencies in contemporary samples (2008–2010) and found a temporal shift toward hatchery allele frequencies in some wild populations. Other populations showed markedly less introgression, despite moderate amounts of hatchery straying. The extent of introgression may reflect similarities in spawning time and life-history traits between hatchery and wild fish, or the degree that hybrids return to a natal spawning area. The source-sink model is a powerful means of detecting low levels of introgression over several generations.

Atlantic salmon populations invaded by farmed escapees: quantifying genetic introgression with a Bayesian approach and SNPs

Background

Many native Atlantic salmon populations have been invaded by domesticated escapees for three decades or longer. However, thus far, the cumulative level of gene-flow that has occurred from farmed to wild salmon has not been reported for any native Atlantic salmon population. The aim of the present study was to investigate temporal genetic stability in native populations, and, quantify gene-flow from farmed salmon that caused genetic changes where they were observed. This was achieved by genotyping historical and contemporary samples from 20 populations covering all of Norway with recently identified single nucleotide polymorphism markers that are collectively diagnostic for farmed and wild salmon. These analyses were combined with analysis of farmed salmon and implementation of Approximate Bayesian computation based simulations.

Results

Five of the populations displayed statistically significant temporal genetic changes. All five of these populations became more similar to a pool of farmed fish with time, strongly suggesting introgression of farmed fish as the primary cause. The remaining 15 populations displayed weak or non-significant temporal genetic changes. Estimated introgression of farmed fish ranged from 2-47% per population using approximate Bayesian computation. Thus, some populations exhibited high degrees of farmed salmon introgression while others were more or less unaffected. The observed frequency of escapees in each population was moderately correlated with estimated introgression per population R² = 0.47 P < 0.001. Genetic isolation by distance existed within the historical and contemporary data sets, however, the among-population level of divergence decreased with time.

Conclusions

This is the first study to quantify cumulative introgression of farmed salmon in any native Atlantic salmon population. The estimations demonstrate that the level of introgression has been population-specific, and that the level of introgression is not solely predicted by the frequency of escapees observed in the population. However, some populations have been strongly admixed with farmed salmon, and these data provide policy makers with unique information to address this situation.

Changes in the genetic structure of Atlantic salmon populations over four decades reveal substantial impacts of stocking and potential resiliency

While the stocking of captive-bred fish has been occurring for decades and has had substantial immediate genetic and evolutionary impacts on wild populations, its long-term consequences have only been weakly investigated. Here, we conducted a spatiotemporal analysis of 1428 Atlantic salmon sampled from 1965 to 2006 in 25 populations throughout France to investigate the influence of stocking on the neutral genetic structure in wild Atlantic salmon (Salmo salar) populations. On the basis of the analysis of 11 microsatellite loci, we found that the overall genetic structure among populations dramatically decreased over the period studied. Admixture rates among populations were highly variable, ranging from a nearly undetectable contribution from donor stocks to total replacement of the native gene pool, suggesting extremely variable impacts of stocking. Depending on population, admixture rates either increased, remained stable, or decreased in samples collected between 1998 and 2006 compared to samples from 1965 to 1987, suggesting either rising, long-lasting or short-term impacts of stocking. We discuss the potential mechanisms contributing to this variability, including the reduced fitness of stocked fish and persistence of wild locally adapted individuals.