Strong parallel differential gene expression induced by hatchery rearing weakly associated with methylation signals in adult Coho Salmon (O. kisutch)

Human activities and resource exploitation led to a massive decline of wild salmonid populations, consequently, numerous conservation programs have been developed to supplement wild populations. However, many studies documented reduced fitness of hatchery-born relative to wild fish. Here, by using both RNA sequencing and Whole Genome Bisulfite Sequencing of hatchery and wild-born adult Coho salmon (Oncorhynchus kisutch) originating from two previously studied river systems, we show that early-life hatchery-rearing environment-induced significant and parallel gene expression differentiation is maintained until Coho come back to their natal river for reproduction. A total of 3,643 genes differentially expressed and 859 coexpressed genes were downregulated in parallel in hatchery-born fish from both rivers relative to their wild congeners. Among those genes, 26 displayed a significant relationship between gene expression and the median gene body methylation and 669 single CpGs displayed a significant correlation between methylation level and the associated gene expression. The link between methylation and gene expression was weak suggesting that DNA methylation is not the only player in mediating hatchery-related expression differences. Yet, significant gene expression differentiation was observed despite 18 months spent in a common environment (i.e., the sea). Finally, the differentiation is observed in parallel in two different river systems, highlighting the fact that early-life environment may account for at least some of the reduced fitness of the hatchery salmon in the wild. These results illustrate the relevance and importance of considering both epigenome and transcriptome to evaluate the costs and benefits of large-scale supplementation programs.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

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.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

The role of recombination on genome-wide patterns of local ancestry exemplified by supplemented brook charr populations

Assessing the immediate and long-term evolutionary consequences of human-mediated hybridization is of major concern for conservation biology. Several studies have documented how selection in interaction with recombination modulates introgression at a genome-wide scale, but few have considered the dynamics of this process within and among chromosomes. Here, we used an exploited freshwater fish, the brook charr (Salvelinus fontinalis), for which decades of stocking practices have resulted in admixture between wild populations and an introduced domestic strain, to assess both the temporal dynamics and local chromosomal variation in domestic ancestry. We provide a detailed picture of the domestic ancestry patterns across the genome using about 33,000 mapped single nucleotide polymorphisms genotyped in 611 individuals from 24 supplemented populations. For each lake, we distinguished early- and late-generation hybrids using information regarding admixture tracts. To assess the selective outcomes following admixture we then evaluated the relationship between recombination and admixture proportions at three different scales: the whole genome, chromosomes and within 2-Mb windows. This allowed us to detect a wide range of evolutionary mechanisms varying along the genome, as reflected by the finding of favoured or disfavoured introgression of domestic haplotypes. Among these, the main factor modulating local ancestry was probably the presence of deleterious recessive mutations in the wild populations, which can be efficiently hidden to selection in the presence of long admixture tracts. Overall, our results emphasize the relevance of taking into consideration local ancestry information to assess both the temporal and the chromosomal variation in local admixture ancestry toward better understanding post-hybridization evolutionary outcomes.

Development of a large SNPs resource and a low-density SNP array for brown trout (Salmo trutta) population genetics

Background

The brown trout (Salmo trutta) is an economically and ecologically important species for which population genetic monitoring is frequently performed. The most commonly used genetic markers for this species are microsatellites and mitochondrial markers that lack replicability among laboratories, and a large genome coverage. An alternative that may be particularly efficient and universal is the development of small to large panels of Single Nucleotide Polymorphism markers (SNPs). Here, we used Restriction site Associated DNA sequences (RADs) markers to identify a set of 12,204 informative SNPs positioned on the brown trout linkage map and suitable for population genetics studies. Then, we used this novel resource to develop a cost-effective array of 192 SNPs (96 × 2) evenly spread on this map. This array was tested for genotyping success in five independent rivers occupied by two main brown trout evolutionary lineages (Atlantic -AT- and Mediterranean -ME-) on a total of 1862 individuals. Moreover, inference of admixture rate with domestic strains and population differentiation were assessed for a small river system (the Taurion River, 190 individuals) and results were compared to a panel of 13 microsatellites.

Results

A high genotyping success was observed for all rivers (< 1% of non-genotyped loci per individual), although some initially used SNP failed to be amplified, probably because of mutations in primers, and were replaced. These SNPs permitted to identify patterns of isolation-by-distance for some rivers. Finally, we found that microsatellite and SNP markers yielded very similar patterns for population differentiation and admixture assessments, with SNPs having better ability to detect introgression and differentiation.

Conclusions

The novel resources provided here opens new perspectives for universality and genome-wide studies in brown trout populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5958-9) contains supplementary material, which is available to authorized users.

Genomic consequences of a recent three-way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human-mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make-up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double-digested restriction site-associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple-way admixture with population genomic data.

Genome-wide nucleotide diversity of hatchery-reared Atlantic and Mediterranean strains of brown trout Salmo trutta compared to wild Mediterranean populations

A genome-wide assessment of diversity is provided for wild Mediterranean brown trout Salmo trutta populations from headwater tributaries of the Orb River and from Atlantic and Mediterranean hatchery-reared strains that have been used for stocking. Double-digest restriction-site-associated DNA sequencing (dd-RADseq) was performed and the efficiency of de novo and reference-mapping approaches to obtain individual genotypes was compared. Large numbers of single nucleotide polymorphism (SNP) markers with similar genome-wide distributions were discovered using both approaches (196 639 v. 121 016 SNPs, respectively), with c. 80% of the loci detected de novo being also found with reference mapping, using the Atlantic salmon Salmo salar genome as a reference. Lower mapping density but larger nucleotide diversity (π) was generally observed near extremities of linkage groups, consistent with regions of residual tetrasomic inheritance observed in salmonids. Genome-wide diversity estimates revealed reduced polymorphism in hatchery strains (π = 0·0040 and π = 0·0029 in Atlantic and Mediterranean strains, respectively) compared to wild populations (π = 0·0049), a pattern that was congruent with allelic richness estimated from microsatellite markers. Finally, pronounced heterozygote deficiency was found in hatchery strains (Atlantic F_IS = 0·18; Mediterranean F_IS = 0·42), indicating that stocking practices may affect the genetic diversity in wild populations. These new genomic resources will provide important tools to define better conservation strategies in S. trutta.

Ancestry and adaptive evolution of anadromous, resident, and adfluvial rainbow trout (Oncorhynchus mykiss) in the San Francisco bay area: application of adaptive genomic variation to conservation in a highly impacted landscape

The streams draining of into San Francisco Bay, California, have been impacted by habitat alteration for over 150 years, and roads, dams, water diversions, and other impediments now block the paths of many aquatic migratory species. These changes can affect the genetic structure of fish populations, as well as driving adaptive evolution to novel environmental conditions. Here, we determine the evolutionary relationships of San Francisco Bay Area steelhead/rainbow trout (Oncorhynchus mykiss) populations and show that (i) they are more closely related to native coastal steelhead than to the California Central Valley lineage, with no evidence of introgression by domesticated hatchery rainbow trout, (ii) populations above and below barriers within watersheds are each other's closest relatives, and (iii) adaptive genomic variation associated with migratory life-history traits in O. mykiss shows substantial evolutionary differences between fish above and below dams. These findings support continued habitat restoration and protection of San Francisco Bay Area O. mykiss populations and demonstrate that ecological conditions in novel habitats above barriers to anadromy influence life-history evolution. We highlight the importance of considering the adaptive landscape in conservation and restoration programs for species living in highly modified habitats, particularly with respect to key life-history traits.