Martin Laporte

RAD sequencing reveals within‐generation polygenic selection in response to anthropogenic organic and metal contamination in North Atlantic Eels

Measuring the effects of selection on the genome imposed by human‐altered environment is currently a major goal in ecological genomics. Given the polygenic basis of most phenotypic traits, quantitative genetic theory predicts that selection is expected to cause subtle allelic changes among covarying loci rather than pronounced changes at few loci of large effects. The goal of this study was to test for the occurrence of polygenic selection in both North Atlantic eels (European Eel, Anguilla anguilla and American Eel, A. rostrata), using a method that searches for covariation among loci that would discriminate eels from ‘control’ vs. ‘polluted’ environments and be associated with specific contaminants acting as putative selective agents. RAD‐seq libraries resulted in 23 659 and 14 755 filtered loci for the European and American Eels, respectively. A total of 142 and 141 covarying markers discriminating European and American Eels from ‘control’ vs. ‘polluted’ sampling localities were obtained using the Random Forest algorithm. Distance‐based redundancy analyses (db‐RDAs) were used to assess the relationships between these covarying markers and concentration of 34 contaminants measured for each individual eel. PCB153, 4′4′DDE and selenium were associated with covarying markers for both species, thus pointing to these contaminants as major selective agents in contaminated sites. Gene enrichment analyses suggested that sterol regulation plays an important role in the differential survival of eels in ‘polluted’ environment. This study illustrates the power of combining methods for detecting signals of polygenic selection and for associating variation of markers with putative selective agents in studies aiming at documenting the dynamics of selection at the genomic level and particularly so in human‐altered environments.

Seascape genomics provides evidence for thermal adaptation and current-mediated population structure in American lobster (Homarus americanus).

Investigating how environmental features shape the genetic structure of populations is crucial for understanding how they are potentially adapted to their habitats, as well as for sound management. In this study, we assessed the relative importance of spatial distribution, ocean currents and sea surface temperature (SST) on patterns of putatively neutral and adaptive genetic variation among American lobster from 19 locations using population differentiation (PD) approaches combined with environmental association (EA) analyses. First, PD approaches (using bayescan, arlequin and outflank) found 28 outlier SNPs putatively under divergent selection and 9770 neutral SNPs in common. Redundancy analysis revealed that spatial distribution, ocean current‐mediated larval connectivity and SST explained 31.7% of the neutral genetic differentiation, with ocean currents driving the majority of this relationship (21.0%). After removing the influence of spatial distribution, no SST were significant for putatively neutral genetic variation whereas minimum annual SST still had a significant impact and explained 8.1% of the putatively adaptive genetic variation. Second, EA analyses (using Pearson correlation tests, bayescenv and lfmm) jointly identified seven SNPs as candidates for thermal adaptation. Covariation at these SNPs was assessed with a spatial multivariate analysis that highlighted a significant temperature association, after accounting for the influence of spatial distribution. Among the 505 candidate SNPs detected by at least one of the three approaches, we discovered three polymorphisms located in genes previously shown to play a role in thermal adaptation. Our results have implications for the management of the American lobster and provide a foundation on which to predict how this species will cope with climate change.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs.

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size‐matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming “dwarf” ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I–IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.

Investigating genomic and phenotypic parallelism between piscivorous and planktivorous lake trout (Salvelinus namaycush) ecotypes by means of RADseq and morphometrics analyses

Repeated adaptive ecological diversification has commonly been reported in fish and has often been associated with trophic niche diversity. The main goal of this study was to investigate the extent of parallelism in the genomic and phenotypic divergence between piscivorous and planktivorous lake trout ecotypes from Laurentian Shield lakes, Canada. This was achieved by documenting the extent of morphological differentiation using geometric morphometrics and linear measurements as well as the pattern of genomic divergence by means of RADseq genotyping (3925 filtered SNPs) in 12 lakes. Our results indicate that the two ecotypes evolved distinct body shape and several linear measurements in parallel. Neutral genetic differentiation was pronounced between all isolated populations (Mean FST = 0.433), indicating no or very limited migration and pronounced genetic drift. Significant genetic differentiation also suggested partial reproductive isolation between ecotypes in the two lakes where they are found in sympatry. Combining different outlier detection methods, we identified 48 SNPs putatively under divergent selection between ecotypes, among which 10 could be annotated and related to functions such as developmental processes and ionic regulation. Finally, our results indicate that parallel morphological divergence is accompanied by both parallel and nonparallel genomic divergence, which is associated with the use of different trophic niches between ecotypes. The results are also discussed in the context of management and conservation of this highly exploited species throughout northern North America.

Parallel epigenetic modifications induced by hatchery rearing in a Pacific Salmon

Significance Captive rearing is known to impact the fitness of individuals released in the wild, but the relative role of genetic vs. nongenetic underlying processes is still debated. We measured genome-wide methylation profiles to document epigenetic differences between Pacific salmon originating from a hatchery and their natural-born congeners in two geographically distant rivers. Our results provide evidence that the epigenetic modifications induced by hatchery rearing provide a potential explanatory mechanism for reduced fitness of hatchery-reared salmon once released in the wild. Wild stocks of Pacific salmonids have experienced sharp declines in abundance over the past century. Consequently, billions of fish are released each year for enhancing abundance and sustaining fisheries. However, the beneficial role of this widely used management practice is highly debated since fitness decrease of hatchery-origin fish in the wild has been documented. Artificial selection in hatcheries has often been invoked as the most likely explanation for reduced fitness, and most studies to date have focused on finding signatures of hatchery-induced selection at the DNA level. We tested an alternative hypothesis, that captive rearing induces epigenetic reprogramming, by comparing genome-wide patterns of methylation and variation at the DNA level in hatchery-reared coho salmon (Oncorhynchus kisutch) with those of their wild counterparts in two geographically distant rivers. We found a highly significant proportion of epigenetic variation explained by the rearing environment that was as high as the one explained by the river of origin. The differentially methylated regions show enrichment for biological functions that may affect the capacity of hatchery-born smolts to migrate successfully in the ocean. Shared epigenetic variation between hatchery-reared salmon provides evidence for parallel epigenetic modifications induced by hatchery rearing in the absence of genetic differentiation between hatchery and natural-origin fish for each river. This study highlights epigenetic modifications induced by captive rearing as a potential explanatory mechanism for reduced fitness in hatchery-reared salmon.

Draft genome of the American eel (ANGUILLA rostrata)

Freshwater eels (Anguilla sp.) have large economic, cultural, ecological and aesthetic importance worldwide, but they suffered more than 90% decline in global stocks over the past few decades. Proper genetic resources, such as sequenced, assembled and annotated genomes, are essential to help plan sustainable recoveries by identifying physiological, biochemical and genetic mechanisms that caused the declines or that may lead to recoveries. Here, we present the first sequenced genome of the American eel. This genome contained 305 043 contigs (N50 = 7397) and 79 209 scaffolds (N50 = 86 641) for a total size of 1.41 Gb, which is in the middle of the range of previous estimations for this species. In addition, protein‐coding regions, including introns and flanking regions, are very well represented in the genome, as 95.2% of the 458 core eukaryotic genes and 98.8% of the 248 ultra‐conserved subset were represented in the assembly and a total of 26 564 genes were annotated for future functional genomics studies. We performed a candidate gene analysis to compare three genes among all three freshwater eel species and, congruent with the phylogenetic relationships, Japanese eel (A. japanica) exhibited the most divergence. Overall, the sequenced genome presented in this study is a crucial addition to the presently available genetic tools to help guide future conservation efforts of freshwater eels.

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.

Convergence in organ size but not energy metabolism enzyme activities among wild Lake Whitefish (Coregonus clupeaformis) species-pairs

The repeated evolution of similar phenotypes by similar mechanisms can be indicative of local adaptation, constraints or biases in the evolutionary process. Little is known about the incidence of physiological convergence in natural populations, so here we test whether energy metabolism in ‘dwarf’ and ‘normal’ Lake Whitefish evolves by similar mechanisms. Prior genomic and transcriptomic studies have found that divergence in energy metabolism is key to local adaptation in whitefish species pairs, but that distinct genetic and transcriptomic changes often underlie phenotypic evolution among lakes. Here, we predicted that traits at higher levels of biological organization, including the activities of energy metabolism enzymes (the product of enzyme concentration and turnover rate) and the relative proportions of metabolically active tissues (heart, liver, skeletal muscle), would show greater convergence than genetic and transcriptomic variation. We compared four whitefish species pairs and found convergence in organ size whereby all dwarf whitefish populations have a higher proportion of red skeletal muscle, three have relatively larger livers and two have relatively larger ventricles than normal fish. On the other hand, hepatic and muscle enzyme activities showed little convergence and were largely dependent on lake of origin. Only the most genetically divergent species pair (Cliff Lake) displayed white muscle enzyme activities matching results from laboratory‐reared normal and dwarf whitefish. Overall, these data show convergence in the evolution of organ size, but not in the activities of candidate enzymes of energy metabolism, which may have evolved mainly as a consequence of demographic or ecological differences among lakes.

Impact of supplementation on deleterious mutation distribution in an exploited salmonid

Deleterious mutations have important implications for the evolutionary trajectories of populations. While several studies recently investigated the dynamics of deleterious mutations in wild populations, no study has yet explored the fate of deleterious mutations in a context of populations managed by supplementation. Here, based on a dataset of nine wild and 15 supplemented Lake Trout populations genotyped at 4,982 single nucleotide polymorphisms (SNP)s by means of genotype by sequencing (GBS), we explored the effect of supplementation on the frequency of putatively deleterious variants. Three main findings are consequential for the management of fish populations. First, an increase in neutral genetic diversity in stocked populations compared with unstocked ones was observed. Second, putatively deleterious mutations were more likely to be found in unstocked than in stocked populations, suggesting a lower efficiency to purge deleterious mutations in unstocked lakes. Third, a population currently used as a major source for supplementation is characterized by several fixed putatively deleterious alleles. Therefore, other source populations with lower abundance of putatively deleterious mutations should be favored as sources of supplementation. We discuss management implications of our results, especially pertaining to the joint identification of neutral and deleterious mutations that could help refining the choice of source and sink populations for supplementation in order to maximize their evolutionary potential and to limit their mutation load.