Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Landscape connectivity and genetic structure in a mainstem and a tributary stonefly (Plecoptera) species using a novel reference genome

Understanding how environmental variation influences population genetic structure can help predict how environmental change influences population connectivity, genetic diversity, and evolutionary potential. We used riverscape genomics modelling to investigate how climatic and habitat variables relate to patterns of genetic variation in two stonefly species, one from mainstem river habitats (Sweltsa coloradensis) and one from tributaries (Sweltsa fidelis) in 40 sites in northwest Montana, USA. We produced a draft genome assembly for S. coloradensis (N50 = 0.251 Mbp, BUSCO > 95% using “insecta_ob9” reference genes). We genotyped 1930 SNPs in 372 individuals for S. coloradensis and 520 SNPs in 153 individuals for S. fidelis. We found higher genetic diversity for S. coloradensis compared to S. fidelis, but nearly identical genetic differentiation among sites within each species (both had global loci median FST = 0.000), despite differences in stream network location. For landscape genomics and testing for selection, we produced a less stringently filtered data set (3454 and 1070 SNPs for S. coloradensis and S. fidelis, respectively). Environmental variables (mean summer precipitation, slope, aspect, mean June stream temperature, land cover type) were correlated with 19 putative adaptive loci for S. coloradensis. but there was only one putative adaptive locus for S. fidelis (correlated with aspect). Interestingly, we also detected potential hybridization between multiple Sweltsa species which has never been previously detected. Studies like ours, that test for adaptive variation in multiple related species are needed to help assess landscape connectivity and the vulnerability of populations and communities to environmental change.

Differential patterns of diversity at neutral and adaptive loci in endangered Rhodeus pseudosericeus populations

Given the fact that threatened species are often composed of isolated small populations, spatial continuity or demography of the populations may be major factors that have shaped the species' genetic diversity. Thus, neutral loci have been the most commonly-used markers in conservation genetics. However, the populations under the influence of different environmental factors may have evolved in response to different selective pressures, which cannot be fully reflected in neutral genetic variation. Rhodeus pseudosericeus, a bitterling species (Acheilognathidae; Cypriniformes) endemic to the Korean Peninsula, are only found in some limited areas of three rivers, Daecheon, Han and Muhan, that flow into the west coast. Here, we genotyped 24 microsatellite loci and two loci (DAB1 and DAB3) of MHC class II peptide-binding β1 domain for 222 individuals collected from seven populations. Our microsatellite analysis revealed distinctive differentiation between the populations of Daecheon and Muhan Rivers and the Han River populations, and populations were structured into two subgroups within the Han River. Apparent positive selection signatures were found in the peptide-binding residues (PBRs) of the MHC loci. The allelic distribution of MHC showed a degree of differentiation between the populations of Daecheon and Muhan Rivers and the Han River populations, partially similar to the results obtained for microsatellites, however showed rather complex patterns among populations in the Han River. Considering the apparent differences in the distribution of supertypes obtained based on the physicochemical differences induced by the polymorphisms of these PBRs, the differentiation in DAB1 between the two regional groups may result in the differences in immune function. No differentiation between these two regions was observed in the supertyping of DAB3, probably indicating that only DAB1 was associated with the response to locally specialized antigenic peptides.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

Temporal dynamics of migration-linked genetic variation are driven by streamflows and riverscape permeability

Landscape permeability is often explored spatially, but may also vary temporally. Landscape permeability, including partial barriers, influences migratory animals that move across the landscape. Partial barriers are common in rivers where barrier passage varies with streamflow. We explore the influence of partial barriers on the spatial and temporal distribution of migration‐linked genotypes of Oncorhynchus mykiss, a salmonid fish with co‐occurring resident and migratory forms, in tributaries to the South Fork Eel River, California, USA, Elder and Fox Creeks. We genotyped >4,000 individuals using RAD‐capture and classified individuals as resident, heterozygous or migratory genotypes using life history‐associated loci. Across four years of study (2014–2017), the permeability of partial barriers varied across dry and wet years. In Elder Creek, the largest waterfall was passable for adults migrating up‐river 4–39 days each year. In this stream, the overall spatial pattern, with fewer migratory genotypes above the waterfall, remained true across dry and wet years (67%–76% of migratory alleles were downstream of the waterfall). We also observed a strong relationship between distance upstream and proportion of migratory alleles. In Fox Creek, the primary barrier is at the mouth, and we found that the migratory allele frequency varied with the annual timing of high flow events. In years when rain events occurred during the peak breeding season, migratory allele frequency was high (60%–68%), but otherwise it was low (30% in two years). We highlight that partial barriers and landscape permeability can be temporally dynamic, and this effect can be observed through changing genotype frequencies in migratory animals.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Telemetry and genetics reveal asymmetric dispersal of a lake-feeding salmonid between inflow and outflow spawning streams at a microgeographic scale

The degree of natal philopatry relative to natal dispersal in animal populations has important demographic and genetic consequences and often varies substantially within species. In salmonid fishes, lakes have been shown to have a strong influence on dispersal and gene flow within catchments; for example, populations spawning in inflow streams are often reproductively isolated and genetically distinct from those spawning in relatively distant outflow streams. Less is known, however, regarding the level of philopatry and genetic differentiation occurring at microgeographic scales, for example, where inflow and outflow streams are separated by very small expanses of lake habitat. Here, we investigated the interplay between genetic differentiation and fine-scale spawning movements of brown trout between their lake-feeding habitat and two spawning streams (one inflow, one outflow, separated by <100 m of lake habitat). Most (69.2%) of the lake-tagged trout subsequently detected during the spawning period were recorded in just one of the two streams, consistent with natal fidelity, while the remainder were detected in both streams, creating an opportunity for these individuals to spawn in both natal and non-natal streams. The latter behavior was supported by genetic sibship analysis, which revealed several half-sibling dyads containing one individual that was sampled as a fry in the outflow and another that was sampled as fry in the inflow. Genetic clustering analyses in conjunction with telemetry data suggested that asymmetrical dispersal patterns were occurring, with natal fidelity being more common among individuals originating from the outflow than the inflow stream. This was corroborated by Bayesian analysis of gene flow, which indicated significantly higher rates of gene flow from the inflow into the outflow than vice versa. Collectively, these results reveal how a combination of telemetry and genetics can identify distinct reproductive behaviors and associated asymmetries in natal dispersal that produce subtle, but nonetheless biologically relevant, population structuring at microgeographic scales.

Differing, multiscale landscape effects on genetic diversity and differentiation in eastern chipmunks

Understanding how habitat loss and fragmentation impact genetic variation is a major goal in landscape genetics, but to date, most studies have focused solely on the correlation between intervening matrix and genetic differentiation at a single spatial scale. Several caveats exist in these study designs, among them is the inability to include measures of genetic diversity in addition to differentiation. Both genetic metrics help predict population persistence, but are expected to function at differing spatial scales, which requires a multiscale investigation. In this study, we sampled two distinct spatial scales in 31 independent landscapes along a gradient of landscape context (i.e., forest amount, configuration, and types of intervening matrix) to investigate how landscape heterogeneity influences genetic diversity and differentiation in the forest-associated eastern chipmunk (Tamias striatus). Overall, quality of intervening matrix was correlated with genetic differentiation at multiple spatial scales, whereas only configuration was associated with regional scale genetic diversity. Habitat amount, in contrast, did not influence genetic differentiation or diversity at either spatial scale. Based on our findings, landscape effects on genetic variation appears to differ based on spatial scale, the type of genetic response variable, and random variation among landscapes, making extrapolation of results from single scale, unreplicated studies difficult. We encourage landscape geneticists to utilize multiscale, replicated landscapes with both genetic diversity, and differentiation to gain a more comprehensive understanding of how habitat loss and fragmentation influence genetic variation.

Standing genetic diversity and selection at functional gene loci are associated with differential invasion success in two non-native fish species

Invasive species are expected to experience a unique combination of high genetic drift due to demographic factors while also experiencing strong selective pressures. The paradigm that reduced genetic diversity should limit the evolutionary potential of invasive species, and thus, their potential for range expansion has received little empirical support, possibly due to the choice of genetic markers. Our goal was to test for effects of genetic drift and selection at functional genetic markers as they relate to the invasion success of two paired invasive goby species, one widespread (successful) and one with limited range expansion (less successful). We genotyped fish using two marker types: single nucleotide polymorphisms (SNPs) in known-function, protein-coding genes and microsatellites to contrast the effects of neutral genetic processes. We identified reduced allelic variation in the invaded range for the less successful tubenose goby. SNPs putatively under selection were responsible for the observed differences in population structure between marker types for round goby (successful) but not tubenose goby (less successful). A higher proportion of functional loci experienced divergent selection for round goby, suggesting increased evolutionary potential in invaded ranges may be associated with round goby's greater invasion success. Genes involved in thermal tolerance were divergent for round goby populations but not tubenose goby, consistent with the hypothesis that invasion success for fish in temperate regions is influenced by capacity for thermal tolerance. Our results highlight the need to incorporate functional genetic markers in studies to better assess evolutionary potential for the improved conservation and management of species.

Functional connectivity in replicated urban landscapes in the land snail (Cornu aspersum)

Urban areas are highly fragmented and thereby exert strong constraints on individual dispersal. Despite this, some species manage to persist in urban areas, such as the garden snail, Cornu aspersum, which is common in cityscapes despite its low mobility. Using landscape genetic approaches, we combined study area replication and multiscale analysis to determine how landscape composition, configuration and connectivity influence snail dispersal across urban areas. At the overall landscape scale, areas with a high percentage of roads decreased genetic differentiation between populations. At the population scale, genetic differentiation was positively linked with building surface, the proportion of borders where wooded patches and roads appeared side by side and the proportion of borders combining wooded patches and other impervious areas. Analyses based on pairwise genetic distances validated the isolation-by-distance and isolation-by-resistance models for this land snail, with an equal fit to least-cost paths and circuit-theory-based models. Each of the 12 landscapes analysed separately yielded specific relations to environmental features, whereas analyses integrating all replicates highlighted general common effects. Our results suggest that urban transport infrastructures facilitate passive snail dispersal. At a local scale, corresponding to active dispersal, unfavourable habitats (wooded and impervious areas) isolate populations. This work upholds the use of replicated landscapes to increase the generalizability of landscape genetics results and shows how multiscale analyses provide insight into scale-dependent processes.

Landscape features along migratory routes influence adaptive genomic variation in anadromous steelhead (Oncorhynchus mykiss)

Organisms typically show evidence of adaptation to features within their local environment. However, many species undergo long-distance dispersal or migration across larger geographic regions that consist of highly heterogeneous habitats. Therefore, selection may influence adaptive genetic variation associated with landscape features at residing sites and along migration routes in migratory species. We tested for genomic adaptation to landscape features at natal spawning sites and along migration paths to the ocean of anadromous steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin. Results from multivariate ordination, gene-environment association and outlier analyses using 24,526 single nucleotide polymorphisms (SNPs) provided evidence that adaptive allele frequencies were more commonly associated with landscape features along migration paths than features at natal sites (91.8% vs. 8.2% of adaptive loci, respectively). Among the 45 landscape variables tested, migration distance to the ocean and mean annual precipitation along migration paths were significantly associated with adaptive genetic variation in three distinct genetic groups. Additionally, variables such as minimum migration water temperature and mean migration slope were significant only in inland stocks of steelhead that migrate up to 1,200 km farther than those near the coast, indicating regional differences in migratory selective pressures. This study provides novel approaches for investigating migratory corridors and some of the first evidence that environment along migration paths can lead to substantial divergent selection. Consequently, our approach to understand genetic adaptation to migration conditions can be applied to other migratory species when migration or dispersal paths are generally known.

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.

Setting spatial conservation priorities despite incomplete data for characterizing metapopulations

Management of spatially structured species poses unique challenges. Despite a strong theoretical foundation, practitioners rarely have sufficient empirical data to evaluate how populations interact. Rather, assumptions about connectivity and source-sink dynamics are often based on incomplete, extrapolated, or modeled data, if such interactions are even considered at all. Therefore, it has been difficult to evaluate whether spatially structured species are meeting conservation goals. We evaluated how estimated metapopulation structure responded to estimates of population sizes and dispersal probabilities and to the set of populations included. We then compared outcomes of alternative management strategies that target conservation of metapopulation processes. We illustrated these concepts for Chinook salmon (Oncorhynchus tshawytscha) in the Snake River, USA. Our description of spatial structure for this metapopulation was consistent with previous characterizations. We found substantial differences in estimated metapopulation structure when we had incomplete information about all populations and when we used different sources of data (three empirical, two modeled) to estimate dispersal, whereas responses to population size estimates were more consistent. Together, these findings suggest that monitoring efforts should target all populations occasionally and populations that play key roles frequently and that multiple types of data should be collected when feasible. When empirical data are incomplete or of uneven quality, analyses using estimates produced from an ensemble of available datasets can help conservation planners and managers weigh near-term options. Doing so, we found trade-offs in connectivity and source dominance in metapopulation-level responses to alternative management strategies that suggest which types of approaches may be inherently less risky.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Identification of landscape features influencing gene flow: How useful are habitat selection models?

Understanding how dispersal patterns are influenced by landscape heterogeneity is critical for modeling species connectivity. Resource selection function (RSF) models are increasingly used in landscape genetics approaches. However, because the ecological factors that drive habitat selection may be different from those influencing dispersal and gene flow, it is important to consider explicit assumptions and spatial scales of measurement. We calculated pairwise genetic distance among 301 Dall's sheep (Ovis dalli dalli) in southcentral Alaska using an intensive noninvasive sampling effort and 15 microsatellite loci. We used multiple regression of distance matrices to assess the correlation of pairwise genetic distance and landscape resistance derived from an RSF, and combinations of landscape features hypothesized to influence dispersal. Dall's sheep gene flow was positively correlated with steep slopes, moderate peak normalized difference vegetation indices (NDVI), and open land cover. Whereas RSF covariates were significant in predicting genetic distance, the RSF model itself was not significantly correlated with Dall's sheep gene flow, suggesting that certain habitat features important during summer (rugged terrain, mid-range elevation) were not influential to effective dispersal. This work underscores that consideration of both habitat selection and landscape genetics models may be useful in developing management strategies to both meet the immediate survival of a species and allow for long-term genetic connectivity.