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

Urbanization drives adaptive evolution in a Neotropical bird

Urbanization has dramatic impacts on natural habitats and such changes may potentially drive local adaptation of urban populations. Behavioral change has been specifically shown to facilitate the fast adaptation of birds to changing environments, but few studies have investigated the genetic mechanisms of this process. Such investigations could provide insights into questions about both evolutionary theory and management of urban populations. In this study, we investigated whether local adaptation has occurred in urban populations of a Neotropical bird species, Coereba flaveola, specifically addressing whether observed behavioral adaptations are correlated to genetic signatures of natural selection. To answer this question, we sampled 24 individuals in urban and rural environments, and searched for selected loci through a genome-scan approach based on RADseq genomic data, generated and assembled using a reference genome for the species. We recovered 46 loci as putative selection outliers, and 30 of them were identified as associated with biological processes possibly related to urban adaptation, such as the regulation of energetic metabolism, regulation of genetic expression, and changes in the immunological system. Moreover, genes involved in the development of the nervous system showed signatures of selection, suggesting a link between behavioral and genetic adaptations. Our findings, in conjunction with similar results in previous studies, support the idea that cities provide a similar selective pressure on urban populations and that behavioral plasticity may be enhanced through genetic changes in urban populations.

Optimization of an amplicon sequencing-based microsatellite panel and protocol for stock identification and kinship inference of lake trout (Salvelinus namaycush)

Kinship-based methods of population assessment such as close-kin mark-recapture require accurate and efficient genotyping methods capable of resolving complex relationships among kin. Inference of such relationships can be difficult using biallelic loci due to the large number of markers required to obtain the necessary power. Sequencing-based microsatellite panels offer an efficient alternative, combining high polymorphism with efficient next-generation methods. Here we construct, optimize, and test one such panel for lake trout (Salvelinus namaycush) using a combination of previously-published loci adapted for sequencing and de novo loci mined from a genome assembly. We performed three rounds of primer optimization, yielding a final panel of 131 loci, followed by testing with two different levels of PCR multiplexing (all primers in one or two groups) and two different reaction volumes (5 and 10 μL). Our results showed that the use of the largest multiplex and smallest reaction volume did not substantially change results, allowing significant cost and time savings. To test panel accuracy, we used both a set of 153 known-origin samples from origins of management interest and a series of hatchery crosses representing nine families with parent-offspring, half-sibling, and largely-unrelated pairs. Our results indicate that sequencing-based microsatellite panels can efficiently and accurately provide the information required for a population genetics analyses including population assignment, calculation of between-population F _ST, and kinship-based population estimation techniques. Such techniques are seeing increasing applications for a wide range of taxa; our findings should provide insight and guidance for the development of the necessary molecular resources.

Further evidence from common garden rearing experiments of heritable traits separating lean and siscowet lake charr (Salvelinus namaycush) ecotypes

Genetic evidence of selection for complex and polygenically regulated phenotypes can easily become masked by neutral population genetic structure and phenotypic plasticity. Without direct evidence of genotype‐phenotype associations it can be difficult to conclude to what degree a phenotype is heritable or a product of environment. Common garden laboratory studies control for environmental stochasticity and help to determine the mechanism that regulate traits. Here we assess lipid content, growth, weight, and length variation in full and hybrid F₁ crosses of deep and shallow water sympatric lake charr ecotypes reared for nine years in a common garden experiment. Redundancy analysis (RDA) and quantitative‐trait‐loci (QTL) genomic scans are used to identify associations between genotypes at 19,714 single nucleotide polymorphisms (SNPs) aligned to the lake charr genome and individual phenotypes to determine the role that genetic inheritance plays in ecotype phenotypic diversity. Lipid content, growth, length, and weight differed significantly among lake charr crosses throughout the experiment suggesting that pedigree plays a large roll in lake charr development. Polygenic scores of 15 SNPs putatively associated with lipid content and/or condition factor indicated that ecotype distinguishing traits are polygenically regulated and additive. A QTL identified on chromosome 38 contained >200 genes, some of which were associated with lipid metabolism and growth, demonstrating the complex nature of ecotype diversity. The results of our common garden study further indicate that lake charr ecotypes observed in nature are predetermined at birth and that ecotypes differ fundamentally in lipid metabolism and growth.

Quantifying nonlinear temporal effects of ethanol preservation on round goby (Neogobius melanostomus) anatomical traits

Geometric morphometrics provides a powerful means of evaluating differences in phenotypic traits among specimens. However, inferences of trait variability can be confounded when measurements are based on preserved samples. We evaluated effects of ethanol preservation on morphology over a 22-week time period for a Laurentian Great Lakes invasive fish, round goby (Neogobius melanostomus, Pallas 1814), using sets of 17 lateral and six dorsal landmarks. We tested whether ethanol preservation affected the magnitude of inter-population variation between individuals collected from lake and river habitats. Generalized least square regression determined that length did not significantly vary through the preservation time series for fish from either population, while mass decreased significantly. Body shape variation was summarized using principal component analysis, which revealed that most preservation-associated changes occurred in the first 14 days. The lateral shape experienced a large magnitude change during the first 24 h in ethanol then only minor changes for the remainder of the study. The dorsal shape began to revert to pre-preservation measurements about 14 days following preservation. Additionally, differences in shape were apparent between the two populations throughout the experiment; however, the magnitude of differences between populations varied depending on whether dorsal or lateral landmarks were considered. Our study demonstrates that tissue responses to ethanol preservation can be more complex than a simple loss of mass, resulting in difficult to predict consequences for geometric morphometric analyses, including variable responses depending on the anatomical region being analyzed.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Detecting genetic signals of selection in heavily bottlenecked reindeer populations by comparing parallel founder events

Founder populations are of special interest to both evolutionary and conservation biologists, but the detection of genetic signals of selection in these populations is challenging due to their demographic history. Geographically separated founder populations likely to have been subjected to similar selection pressures provide an ideal but rare opportunity to overcome these challenges. Here we take advantage of such a situation generated when small, isolated founder populations of reindeer were established on the island of South Georgia, and using this system we look for empirical evidence of selection overcoming strong genetic drift. We generated a 70 k ddRADseq single nucleotide polymorphism database for the two parallel reindeer founder populations and screened for signatures of soft sweeps. We find evidence for a genomic region under selection shared among the two populations, and support our findings with Wright-Fisher model simulations to assess the power and specificity of interpopulation selection scans-namely Bayescan, OutFLANK, PCadapt and a newly developed scan called Genome Wide Differentiation Scan (GWDS)-in the context of pairwise source-founder comparisons. Our simulations indicate that loci under selection in small founder populations are most probably detected by GWDS, and strengthen the hypothesis that the outlier region represents a true locus under selection. We explore possible, relevant functional roles for genes in linkage with the detected outlier loci.

Among-individual diet variation within a lake trout ecotype: Lack of stability of niche use

In a polyphenic species, differences in resource use are expected among ecotypes, and homogeneity in resource use is expected within an ecotype. Yet, using a broad resource spectrum has been identified as a strategy for fishes living in unproductive northern environments, where food is patchily distributed and ephemeral. We investigated whether specialization of trophic resources by individuals occurred within the generalist piscivore ecotype of lake trout from Great Bear Lake, Canada, reflective of a form of diversity. Four distinct dietary patterns of resource use within this lake trout ecotype were detected from fatty acid composition, with some variation linked to spatial patterns within Great Bear Lake. Feeding habits of different groups within the ecotype were not associated with detectable morphological or genetic differentiation, suggesting that behavioral plasticity caused the trophic differences. A low level of genetic differentiation was detected between exceptionally large-sized individuals and other piscivore individuals. We demonstrated that individual trophic specialization can occur within an ecotype inhabiting a geologically young system (8,000-10,000 yr BP), a lake that sustains high levels of phenotypic diversity of lake trout overall. The characterization of niche use among individuals, as done in this study, is necessary to understand the role that individual variation can play at the beginning of differentiation processes.

Genetic Structure and Population Demography of White-Spotted Charr in the Upstream Watershed of a Large Dam

White-spotted charr (Salvelinus leucomaenis leucomaenis) is an anadromous fish that has been severely harmed by human land-use development, particularly through habitat fragmentation. However, the anthropogenic impacts on populations of this species have not been evaluated, except those on small dammed-off populations. Using multiplexed ISSR genotyping by sequencing, we investigated the genetic structure of white-spotted charr in four tributaries in the upper section of the Kanayama Dam in the Sorachi River, Hokkaido Island, Japan. There were no distinct genetic structures (FST = 0.014), probably because some active individuals migrate frequently among tributaries. By model-flexible demographic simulation, historical changes in the effective population size were inferred. The result indicates that the population size has decreased since the end of the last glacial period, with three major population decline events, including recent declines that were probably associated with recent human activities. Nevertheless, populations in the watershed upstream of the Kanayama Dam are still expected to be at low risk of immediate extinction, owing to the large watershed size and the limited number of small check dams. An effective conservation measure for sustaining the white-spotted charr population is to maintain high connectivity between tributaries, such as by providing fishways in check dams during construction.

Seascape genomics reveals population isolation in the reef-building honeycomb worm, Sabellaria alveolata (L.)

BACKGROUND

Under the threat of climate change populations can disperse, acclimatise or evolve in order to avoid fitness loss. In light of this, it is important to understand neutral gene flow patterns as a measure of dispersal potential, but also adaptive genetic variation as a measure of evolutionary potential. In order to assess genetic variation and how this relates to environment in the honeycomb worm (Sabellaria alveolata (L.)), a reef-building polychaete that supports high biodiversity, we carried out RAD sequencing using individuals from along its complete latitudinal range. Patterns of neutral population genetic structure were compared to larval dispersal as predicted by ocean circulation modelling, and outlier analyses and genotype-environment association tests were used to attempt to identify loci under selection in relation to local temperature data.

RESULTS

We genotyped 482 filtered SNPs, from 68 individuals across nine sites, 27 of which were identified as outliers using BAYESCAN and ARLEQUIN. All outlier loci were potentially under balancing selection, despite previous evidence of local adaptation in the system. Limited gene flow was observed among reef-sites (FST = 0.28 ± 0.10), in line with the low dispersal potential identified by the larval dispersal models. The North Atlantic reef emerged as a distinct population and this was linked to high local larval retention and the effect of the North Atlantic Current on dispersal.

CONCLUSIONS

As an isolated population, with limited potential for natural genetic or demographic augmentation from other reefs, the North Atlantic site warrants conservation attention in order to preserve not only this species, but above all the crucial functional ecological roles that are associated with their bioconstructions. Our study highlights the utility of using seascape genomics to identify populations of conservation concern.

Contrasting patterns of local adaptation along climatic gradients between a sympatric parasitic and autotrophic tree species

Sympatric tree species are subject to similar climatic drivers, posing a question as to whether they display comparable adaptive responses. However, no study has explicitly examined local adaptation of co-occurring parasitic and autotrophic plant species to the abiotic environment. Here we test the hypotheses that a generalist parasitic tree would display a weaker signal of selection and that genomic variation would associate with fewer climatic variables (particularly precipitation) but have similar spatial patterns to a sympatric autotrophic tree species. To test these hypotheses, we collected samples from 17 sites across the range of two tree species, the hemiparasite Nuytsia floribunda (n = 264) and sympatric autotroph Melaleuca rhaphiophylla (n = 272). We obtained 5,531 high-quality genome-wide single nucleotide polymorphisms (SNPs) for M. rhaphiophylla and 6,727 SNPs for N. floribunda using DArTseq genome scan technology. Population differentiation and environmental association approaches were used to identify signals of selection. Generalized dissimilarly modelling was used to detect climatic and spatial patterns of local adaptation across climatic gradients. Overall, 322 SNPs were identified as putatively adaptive for the autotroph, while only 57 SNPs were identified for the parasitic species. We found genomic variation to associate with different sets of bioclimatic variables for each species, with precipitation relatively less important for the parasite. Spatial patterns of predicted adaptive variability were different and indicate that co-occurring species with disparate life history traits may not respond equally to selective pressures (i.e., temperature and precipitation). Together, these findings provide insight into local adaptation of sympatric parasitic and autotrophic tree species to abiotic environments.

Post-glacial establishment of locally adapted fish populations over a steep salinity gradient

Studies of colonization of new habitats that appear from rapidly changing environments are interesting and highly relevant to our understanding of divergence and speciation. Here, we analyse phenotypic and genetic variation involved in the successful establishment of a marine fish (sand goby, Pomatoschistus minutus) over a steep salinity drop from 35 PSU in the North Sea (NE Atlantic) to two PSU in the inner parts of the post-glacial Baltic Sea. We first show that populations are adapted to local salinity in a key reproductive trait, the proportion of motile sperm. Thereafter, we show that genome variation at 22,190 single nucleotide polymorphisms (SNPs) shows strong differentiation among populations along the gradient. Sequences containing outlier SNPs and transcriptome sequences, mapped to a draft genome, reveal associations with genes with relevant functions for adaptation in this environment but without overall evidence of functional enrichment. The many contigs involved suggest polygenic differentiation. We trace the origin of this differentiation using demographic modelling and find the most likely scenario is that at least part of the genetic differentiation is older than the Baltic Sea and is a result of isolation of two lineages prior to the current contact over the North Sea-Baltic Sea transition zone.

Mapping of Adaptive Traits Enabled by a High-Density Linkage Map for Lake Trout

Understanding the genomic basis of adaptative intraspecific phenotypic variation is a central goal in conservation genetics and evolutionary biology. Lake trout (Salvelinus namaycush) are an excellent species for addressing the genetic basis for adaptive variation because they express a striking degree of ecophenotypic variation across their range; however, necessary genomic resources are lacking. Here we utilize recently-developed analytical methods and sequencing technologies to (1) construct a high-density linkage and centromere map for lake trout, (2) identify loci underlying variation in traits that differentiate lake trout ecophenotypes and populations, (3) determine the location of the lake trout sex determination locus, and (4) identify chromosomal homologies between lake trout and other salmonids of varying divergence. The resulting linkage map contains 15,740 single nucleotide polymorphisms (SNPs) mapped to 42 linkage groups, likely representing the 42 lake trout chromosomes. Female and male linkage group lengths ranged from 43.07 to 134.64 centimorgans, and 1.97 to 92.87 centimorgans, respectively. We improved the map by determining coordinates for 41 of 42 centromeres, resulting in a map with 8 metacentric chromosomes and 34 acrocentric or telocentric chromosomes. We use the map to localize the sex determination locus and multiple quantitative trait loci (QTL) associated with intraspecific phenotypic divergence including traits related to growth and body condition, patterns of skin pigmentation, and two composite geomorphometric variables quantifying body shape. Two QTL for the presence of vermiculations and spots mapped with high certainty to an arm of linkage group Sna3, growth related traits mapped to two QTL on linkage groups Sna1 and Sna12, and putative body shape QTL were detected on six separate linkage groups. The sex determination locus was mapped to Sna4 with high confidence. Synteny analysis revealed that lake trout and congener Arctic char (Salvelinus alpinus) are likely differentiated by three or four chromosomal fissions, possibly one chromosomal fusion, and 6 or more large inversions. Combining centromere mapping information with putative inversion coordinates revealed that the majority of detected inversions differentiating lake trout from other salmonids are pericentric and located on acrocentric and telocentric linkage groups. Our results suggest that speciation and adaptive divergence within the genus Salvelinus may have been associated with multiple pericentric inversions occurring primarily on acrocentric and telocentric chromosomes. The linkage map presented here will be a critical resource for advancing conservation oriented genomic research on lake trout and exploring chromosomal evolution within and between salmonid species.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Gene expression vs. sequence divergence: comparative transcriptome sequencing among natural Rhinolophus ferrumequinum populations with different acoustic phenotypes

BACKGROUND

Although the sensory drive hypothesis can explain the geographic variation in echolocation frequencies of some bat species, the molecular mechanisms underlying this phenomenon are still unclear. The three lineages of greater horseshoe bat (Rhinolophus ferrumequinum) in China (northeast, central-east, and southwest) have significant geographic variation in resting frequencies (RF) of echolocation calls. Because their cochleae have an acoustic fovea that is highly sensitive to a narrow range of frequencies, we reported the transcriptomes of cochleae collected from three genetic lineages of R. ferrumequinum, which is an ideal organism for studying geographic variation in echolocation signals, and tried to understand the mechanisms behind this bat phenomenon by analyzing gene expression and sequence variation.

RESULTS

A total of 8190 differentially expressed genes (DEGs) were identified. We identified five modules from all DEGs that were significantly related to RF or forearm length (FL). DEGs in the RF-related modules were significantly enriched in the gene categories involved in neural activity, learning, and response to sound. DEGs in the FL-related modules were significantly enriched in the pathways related to muscle and actin functions. Using 21,945 single nucleotide polymorphisms, we identified 18 candidate unigenes associated with hearing, five of which were differentially expressed among the three populations. Additionally, the gene ERBB4, which regulates diverse cellular processes in the inner ear such as cell proliferation and differentiation, was in the largest module. We also found 49 unigenes that were under positive selection from 4105 one-to-one orthologous gene pairs between the three R. ferrumequinum lineages and three other Chiroptera species.

CONCLUSIONS

The variability of gene expression and sequence divergence at the molecular level might provide evidence that can help elucidate the genetic basis of geographic variation in echolocation signals of greater horseshoe bats.

Evolution and diversity of two cisco forms in an outlet of glacial Lake Algonquin

The diversity of Laurentian Great Lakes ciscoes (Coregonus artedi, sensu lato) arose via repeated local adaptive divergence including deepwater ciscoes that are now extirpated or threatened. The nigripinnis form, or Blackfin Cisco, is extirpated from the Great Lakes and remains only in Lake Nipigon. Putative nigripinnis populations were recently discovered in sympatry with artedi in a historical drainage system of glacial Lake Algonquin, the precursor of lakes Michigan and Huron. Given the apparent convergence on Great Lakes form, we labeled this form blackfin. Here, we test the hypothesis that nigripinnis may have colonized this area from the Great Lakes as a distinct lineage. It would then represent a relict occurrence of the historical diversity of Great Lakes ciscoes. Alternatively, blackfin could have evolved in situ in several lakes. We captured more than 600 individuals in the benthic or pelagic habitat in 14 lakes in or near Algonquin Provincial Park (Ontario, Canada). Fish were compared based on habitat, morphology, and genetic variation at 6,676 SNPs. Contrary to our expectations, both cisco and blackfin belonged to an Atlantic lineage that colonized the area from the east, not from the Great Lakes. Sympatric cisco and blackfin were closely related while fish from different lakes were genetically differentiated, strongly suggesting the repeated in situ origin of each form. Across lakes, there was a continuum of ecological, morphological, and genetic differentiation that could be associated with alternative resources and lake characteristics. This study uncovers a new component of cisco diversity in inland lakes of Canada that evolved independently from ciscoes of the Laurentian Great lakes. The diversity of cisco revealed in this study and across their Canadian range presents a challenge for designating conservation units at the intraspecific level within the framework of the Committee on the Status of Endangered Wildlife in Canada (COSEWIC).

Impacts of acidification on brown trout Salmo trutta populations and the contribution of stocking to population recovery and genetic diversity

Anthropogenic acidification in SW-Scotland, from the early 19th Century onwards, led to the extinction of several loch (lake) brown trout (Salmo trutta) populations and substantial reductions in numbers in many others. Higher altitude populations with no stocking influence, which are isolated above natural and artificial barriers and subjected to the greatest effect of acidification, exhibited the least intrapopulation genetic diversity (34% of the allelic richness of the populations accessible to anadromous S. trutta). These, however, were characterised by the greatest interpopulation divergence (highest pairwise D_EST 0.61 and F_ST 0.53 in contemporary samples) based on 16 microsatellite loci and are among the most differentiated S. trutta populations in NW-Europe. Five lochs above impassable waterfalls, where S. trutta were thought to be extinct, are documented as having been stocked in the late 1980s or 1990s. All five lochs now support self-sustaining S. trutta populations; three as a direct result of restoration stocking and two adjoining lochs largely arising from a small remnant wild population in one, but with some stocking input. The genetically unique Loch Grannoch S. trutta, which has been shown to have a heritable increased tolerance to acid conditions, was successfully used as a donor stock to restore populations in two acidic lochs. Loch Fleet S. trutta, which were re-established from four separate donor sources in the late 1980s, showed differential contribution from these ancestors and a higher genetic diversity than all 17 natural loch populations examined in the area. Genetically distinct inlet and outlet spawning S. trutta populations were found in this loch. Three genetically distinct sympatric populations of S. trutta were identified in Loch Grannoch, most likely representing recruitment from the three main spawning rivers. A distinct genetic signature of Loch Leven S. trutta, the progenitor of many Scottish farm strains, facilitated detection of stocking with these strains. One artificially created loch was shown to have a population genetically very similar to Loch Leven S. trutta. In spite of recorded historical supplemental stocking with Loch Leven derived farm strains, much of the indigenous S. trutta genetic diversity in the area remains intact, aside from the effects of acidification induced bottlenecks. Overall genetic diversity and extant populations have been increased by allochthonous stocking.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

The role of ecotype-environment interactions in intraspecific trophic niche partitioning subsequent to stocking

Worldwide, stocking of fish represents a valuable tool for conservation and maintenance of species exploited by recreational fishing. Releases of hatchery-reared fish are more and more recognized to have numerous demographic, ecological, and genetic impacts on wild populations. However, consequences on intraspecific trophic relationships have rarely been investigated. In this study, we assessed the impacts of supplementation stocking and resulting introgressive hybridization on the trophic niches occupied by stocked, local, and hybrid lake trout (Salvelinus namaycush) within populations of piscivorous and planktivorous ecotypes stocked from a wild piscivorous source population. We compared trophic niches using stable isotope analysis (δ¹³ C and δ¹⁵ N) and trophic position among the three genetic origins. Putative genetic effects were tested with phenotype-genotype association of "life history" ecological traits (body size, growth rate, condition index, and trophic niche) and genotypes (RADseq SNP markers) using redundant discriminant analysis (RDA). Results showed that sympatry resulting from the stocking of contrasting ecotypes is a risk factor for niche partitioning. Planktivorous populations are more susceptible to niche partitioning, by competitive exclusion of the local fish from a littoral niche to an alternative pelagic/profundal niche. Observed niche partitioning is probably a manifestation of competitive interactions between ecotypes. Our results emphasize that ecotypic variation should be considered for more efficient management and conservation practices and in order to mitigate negative impact of supplementation stocking.

The ecology of sexual dimorphism in size and shape of the freshwater blenny Salaria fluviatilis

Sexual selection is considered the major cause of sexual dimorphism, but recent observations suggest that natural selection may play a more important role in the evolution of sex differentiation than previously recognized. Therefore, studying the trade-offs between natural selection and sexual selection is crucial to a better understanding of the ecology underlying the evolution of sexual dimorphism. The freshwater blenny Salaria fluviatilis, a fish inhabiting lakes and rivers around the Mediterranean Sea, displays strong sexual dimorphism in size, shape, and behavior (i.e., larger body and head size for males and higher swimming requirements for females during the reproductive period). We tested for differences in sexual dimorphism in size and shape between the populations from lake and river habitats with the goal of identifying the trade-offs between natural and sexual selection that underlie variations in sexual dimorphism in this species. Our results show i) differences in sexual size dimorphism (SSizeD) in accordance to Rensch's rule (i.e., larger individuals in rivers associated with higher SSizeD), and ii) a decrease in shape differentiation between males and females in lake populations. Together, this suggests that the different environmental conditions between lake and river habitats (e.g., resource limitations, predation pressure, water velocity) affect the relative importance of sexual selection in the display of sexual dimorphism within the species. This study highlights the importance of considering the environmental conditions to which populations are exposed to better understand the ecology underlying the evolution of sexual dimorphism.

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.

Supplementation stocking of Lake Trout (Salvelinus namaycush) in small boreal lakes: Ecotypes influence on growth and condition

Supplementation stocking is a commonly used management tool to sustain exploited fish populations. Possible negative consequences of supplementation on local stocks are a concern for the conservation of wild fish populations. However, the direct impacts of supplementation on life history traits of local populations have rarely been investigated. In addition, intraspecific hybridization between contrasting ecotypes (planktivorous and piscivorous) has been seldom considered in supplementation plans. Here, we combined genetic (genotype-by-sequencing analysis) and life history traits to document the effects of supplementation on maximum length, growth rates, body condition and genetic admixture in stocked populations of two Lake Trout ecotypes from small boreal lakes in Quebec and Ontario, Canada. In both ecotypes, the length of stocked individuals was greater than local individuals and, in planktivorous-stocked populations, most stocked fish exhibited a planktivorous-like growth while 20% of fish exhibited piscivorous-like growth. The body condition index was positively related to the proportion of local genetic background, but this pattern was only observed in stocked planktivorous populations. We conclude that interactions and hybridization between contrasting ecotypes is a risk that could result in deleterious impacts and possible outbreeding depression. We discuss the implications of these findings for supplementation stocking.

From top to bottom: Do Lake Trout diversify along a depth gradient in Great Bear Lake, NT, Canada?

Depth is usually considered the main driver of Lake Trout intraspecific diversity across lakes in North America. Given that Great Bear Lake is one of the largest and deepest freshwater systems in North America, we predicted that Lake Trout intraspecific diversity to be organized along a depth axis within this system. Thus, we investigated whether a deep-water morph of Lake Trout co-existed with four shallow-water morphs previously described in Great Bear Lake. Morphology, neutral genetic variation, isotopic niches, and life-history traits of Lake Trout across depths (0-150 m) were compared among morphs. Due to the propensity of Lake Trout with high levels of morphological diversity to occupy multiple habitat niches, a novel multivariate grouping method using a suite of composite variables was applied in addition to two other commonly used grouping methods to classify individuals. Depth alone did not explain Lake Trout diversity in Great Bear Lake; a distinct fifth deep-water morph was not found. Rather, Lake Trout diversity followed an ecological continuum, with some evidence for adaptation to local conditions in deep-water habitat. Overall, trout caught from deep-water showed low levels of genetic and phenotypic differentiation from shallow-water trout, and displayed higher lipid content (C:N ratio) and occupied a higher trophic level that suggested an potential increase of piscivory (including cannibalism) than the previously described four morphs. Why phenotypic divergence between shallow- and deep-water Lake Trout was low is unknown, especially when the potential for phenotypic variation should be high in deep and large Great Bear Lake. Given that variation in complexity of freshwater environments has dramatic consequences for divergence, variation in the complexity in Great Bear Lake (i.e., shallow being more complex than deep), may explain the observed dichotomy in the expression of intraspecific phenotypic diversity between shallow- vs. deep-water habitats. The ambiguity surrounding mechanisms driving divergence of Lake Trout in Great Bear Lake should be seen as reflective of the highly variable nature of ecological opportunity and divergent natural selection itself.

Life histories and ecotype conservation in an adaptive vertebrate: Genetic constitution of piscivorous brown trout covaries with habitat stability

Ecotype variation in species exhibiting different life history strategies may reflect heritable adaptations to optimize reproductive success, and potential for speciation. Traditionally, ecotypes have, however, been defined by morphometrics and life history characteristics, which may be confounded with individual plasticity. Here, we use the widely distributed and polytypic freshwater fish species brown trout (Salmo trutta) as a model to study piscivorous life history and its genetic characteristics in environmentally contrasting habitats; a large lake ecosystem with one major large and stable tributary, and several small tributaries. Data from 550 fish and 13 polymorphic microsatellites (He = 0.67) indicated ecotype-specific genetic differentiation (θ = 0.0170, p < .0001) among Bayesian assigned small riverine resident and large, lake migrating brown trout (>35 cm), but only in the large tributary. In contrast, large trout did not constitute a distinct genetic group in small tributaries, or across riverine sites. Whereas life history data suggest a small, river resident and a large migratory piscivorous ecotype in all studied tributaries, genetic data indicated that a genetically distinct piscivorous ecotype is more likely to evolve in the large and relatively more stable river habitat. In the smaller tributaries, ecotypes apparently resulted from individual plasticity. Whether different life histories and ecotypes result from individual plasticity or define different genetic types, have important consequence for conservation strategies.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.