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

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Genome-Wide RADseq Reveals Genetic Differentiation of Wild and Cultured Populations of Large Yellow Croaker

Larimichthys crocea (also known as the large yellow croaker) is one of the most economically important marine fishes in China, and research on the ecology and genetics of this species is of immense significance. In this study, we performed restriction site-associated DNA sequencing (RAD-seq) of 54 individuals collected from four sites in China to analyze the genetic structure and diversity of large yellow croaker at the genome level. It revealed that the large yellow croaker populations in the Ningde and Zhoushan coastal waters can be clearly distinguished. Different genetic diversity indices were used to analyze the genetic diversity of the large yellow croaker, which showed that there was a differentiation trend between the wild and farmed populations in Ningde. Moreover, we identified genetically differentiated genomic regions between the populations. GO gene enrichment analysis identified genes that are related to fatty acid metabolism and growth. These findings enhance our understanding of genetic differentiation and adaptation to different living environments, providing a theoretical basis for the preservation and restoration of the genetic resources of the large yellow croaker.

Is now the time? Review of genetic rescue as a conservation tool for brook trout

Brook trout populations have been declining throughout their native range in the east coast of the United States. Many populations are now distributed in small, isolated habitat patches where low genetic diversity and high rates of inbreeding reduce contemporary viability and long-term adaptive potential. Although human-assisted gene flow could theoretically improve conservation outcomes through genetic rescue, there is widespread hesitancy to use this tool to support brook trout conservation. Here, we review the major uncertainties that have limited genetic rescue from being considered as a viable conservation tool for isolated brook trout populations and compare the risks of genetic rescue with other management alternatives. Drawing on theoretical and empirical studies, we discuss methods for implementing genetic rescue in brook trout that could yield long-term evolutionary benefits while avoiding negative fitness effects associated with outbreeding depression and the spread of maladapted alleles. We also highlight the potential for future collaborative efforts to accelerate our understanding of genetic rescue as a viable tool for conservation. Ultimately, while we acknowledge that genetic rescue is not without risk, we emphasize the merits that this tool offers for protecting and propagating adaptive potential and improving species' resilience to rapid environmental change.

Environmental selection, rather than neutral processes, best explain regional patterns of diversity in a tropical rainforest fish

To conserve the high functional and genetic variation in hotspots such as tropical rainforests, it is essential to understand the forces driving and maintaining biodiversity. We asked to what extent environmental gradients and terrain structure affect morphological and genomic variation across the wet tropical distribution of an Australian rainbowfish, Melanotaenia splendida splendida. We used an integrative riverscape genomics and morphometrics framework to assess the influence of these factors on both putative adaptive and non-adaptive spatial divergence. We found that neutral genetic population structure was largely explainable by restricted gene flow among drainages. However, environmental associations revealed that ecological variables had a similar power to explain overall genetic variation, and greater power to explain body shape variation, than the included neutral covariables. Hydrological and thermal variables were the strongest environmental predictors and were correlated with traits previously linked to heritable habitat-associated dimorphism in rainbowfishes. In addition, climate-associated genetic variation was significantly associated with morphology, supporting heritability of shape variation. These results support the inference of evolved functional differences among localities, and the importance of hydroclimate in early stages of diversification. We expect that substantial evolutionary responses will be required in tropical rainforest endemics to mitigate local fitness losses due to changing climates.

Species delimitation using genomic data to resolve taxonomic uncertainties in a speciation continuum of pelagic seabirds

Speciation is a continuous and complex process shaped by the interaction of numerous evolutionary forces. Despite the continuous nature of the speciation process, the implementation of conservation policies relies on the delimitation of species and evolutionary significant units (ESUs). Puffinus shearwaters are globally distributed and threatened pelagic seabirds. Due to remarkable morphological status the group has been under intense taxonomic debate for the past three decades. Here, we use double digest Restriction-Site Associated DNA sequencing (ddRAD-Seq) to genotype species and subspecies of North Atlantic and Mediterranean Puffinus shearwaters across their entire geographical range. We assess the phylogenetic relationships and population structure among and within the group, evaluate species boundaries, and characterise the genomic landscape of divergence. We find that current taxonomies are not supported by genomic data and propose a more accurate taxonomy by integrating genomic information with other sources of evidence. Our results show that several taxon pairs are at different stages of a speciation continuum. Our study emphasises the potential of genomic data to resolve taxonomic uncertainties, which can help to focus management actions on relevant taxa, even if they do not necessarily coincide with the taxonomic rank of species.

Give and Take: Effects of Genetic Admixture on Mutation Load in Endangered Florida Panthers

Genetic admixture is a biological event inherent to genetic rescue programs aimed at the long-term conservation of endangered wildlife. Although the success of such programs can be measured by the increase in genetic diversity and fitness of subsequent admixed individuals, predictions supporting admixture costs to fitness due to the introduction of novel deleterious alleles are necessary. Here, we analyzed nonsynonymous variation from conserved genes to quantify and compare levels of mutation load (i.e., proportion of deleterious alleles and genotypes carrying these alleles) among endangered Florida panthers and non-endangered Texas pumas. Specifically, we used canonical (i.e., non-admixed) Florida panthers, Texas pumas, and F1 (canonical Florida x Texas) panthers dating from a genetic rescue program and Everglades National Park panthers with Central American ancestry resulting from an earlier admixture event. We found neither genetic drift nor selection significantly reduced overall proportions of deleterious alleles in the severely bottlenecked canonical Florida panthers. Nevertheless, the deleterious alleles identified were distributed into a disproportionately high number of homozygous genotypes due to close inbreeding in this group. Conversely, admixed Florida panthers (either with Texas or Central American ancestry) presented reduced levels of homozygous genotypes carrying deleterious alleles but increased levels of heterozygous genotypes carrying these variants relative to canonical Florida panthers. Although admixture is likely to alleviate the load of standing deleterious variation present in homozygous genotypes, our results suggest introduced novel deleterious alleles (temporarily present in heterozygous state) in genetically rescued populations could potentially be expressed in subsequent generations if their effective sizes remain small.

Genetic divergence and range expansion in a western North Pacific coral

Coral poleward range expansions have recently been observed in response to warming oceans. Range expansion can lead to reduced genetic diversity and increased frequency of deleterious mutations that were rare in core populations, potentially limiting the ability for adaptation and persistence in novel environments. Successful expansions that overcome these founder effects and colonize new habitat have been attributed to multiple introductions from different sources, hybridization with native populations, or rapid adaptive evolution. Here, we investigate population genomic patterns of the reef-building coral Acropora hyacinthus along a latitudinal cline that includes a well-established range expansion front in Japan using 2b-RAD sequencing. A total of 184 coral samples were collected across seven sites spanning from ~24°N to near its northern range front at ~33°N. We uncover the presence of three cryptic lineages of A. hyacinthus, which occupy discrete reefs within this region. Only one lineage is present along the expansion front and we find evidence for its historical occupation of marginal habitats. Within this lineage we also find evidence of bottleneck pressures associated with expansion events including higher clonality, increased linkage disequilibrium, and lower genetic diversity in range edge populations compared to core populations. Asymmetric migration between populations was also detected with lower migration from edge sites. Lastly, we describe genomic signatures of local adaptation potentially attributed to lower winter temperatures experienced at the more recently expanded northern populations. Together these data illuminate the genomic consequences of range expansion in a coral and highlight how adaptation to discrete environments along expansion fronts may facilitate further range expansion in this temperate coral lineage.

Fitness Effects of Mutations: An Assessment of PROVEAN Predictions Using Mutation Accumulation Data

Predicting fitness in natural populations is a major challenge in biology. It may be possible to leverage fast-accumulating genomic data sets to infer the fitness effects of mutant alleles, allowing evolutionary questions to be addressed in any organism. In this paper, we investigate the utility of one such tool, called PROVEAN. This program compares a query sequence with existing data to provide an alignment-based score for any protein variant, with scores categorized as neutral or deleterious based on a pre-set threshold. PROVEAN has been used widely in evolutionary studies, for example, to estimate mutation load in natural populations, but has not been formally tested as a predictor of aggregate mutational effects on fitness. Using three large published data sets on the genome sequences of laboratory mutation accumulation lines, we assessed how well PROVEAN predicted the actual fitness patterns observed, relative to other metrics. In most cases, we find that a simple count of the total number of mutant proteins is a better predictor of fitness than the number of proteins with variants scored as deleterious by PROVEAN. We also find that the sum of all mutant protein scores explains variation in fitness better than the number of mutant proteins in one of the data sets. We discuss the implications of these results for studies of populations in the wild.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Fish out of water: Genomic insights into persistence of rainbowfish populations in the desert

How populations of aquatic fauna persist in extreme desert environments is an enigma. Individuals often breed and disperse during favorable conditions. Theory predicts that adaptive capacity should be low in small populations, such as in desert fishes. We integrated satellite-derived surface water data and population genomic diversity from 20,294 single-nucleotide polymorphisms across 344 individuals to understand metapopulation persistence of the desert rainbowfish (Melanotaenia splendida tatei) in central Australia. Desert rainbowfish showed very small effective population sizes, especially at peripheral populations, and low connectivity between river catchments. Yet, there was no evidence of population-level inbreeding and a signal of possible adaptive divergence associated with aridity was detected. Candidate genes for local adaptation included functions related to environmental cues and stressful conditions. Eco-evolutionary modeling showed that positive selection in refugial subpopulations combined with connectivity during flood periods can enable retention of adaptive diversity. Our study suggests that adaptive variation can be maintained in small populations and integrate with neutral metapopulation processes to allow persistence in the desert.

One Step Away From Extinction: A Population Genomic Analysis of A Narrow Endemic, Tropical Plant Species

Intraspecific genetic variation plays a fundamental role in maintaining the evolutionary potential of wild populations. Hence, the assessment of genetic diversity patterns becomes essential to guide biodiversity conservation policies, particularly for threatened species. To inform management strategies for conservation of Mimosa catharinensis - a narrow endemic, critically endangered plant species - we identified 1,497 unlinked SNP markers derived from a reduced representation sequencing method (i.e., double digest restriction site associated DNA sequencing, or ddRADseq). This set of molecular markers was employed to assess intrapopulation genetic parameters and the demographic history of one extremely small population of M. catharinensis (N=33) located in the Brazilian Atlantic Forest. Contrary to what is expected for narrow endemic and threatened species with small population sizes, we observed a moderate level of genetic diversity for M. catharinensis [uH E(0%missing data)=0.205, 95% CI (0.160, 0.250); uH E(30%missing data)=0.233, 95% CI (0.174, 0.292)]. Interestingly, M. catharinensis, which is a lianescent shrub with no indication of seed production for at least two decades, presented high levels of outcrossing [t (0%missing data)=0.883, SE±0.0483; t (30%missing data)=0.909, SE±0.011] and an apparent absence of inbreeding [F (0%missing data)=-0.145, 95% CI (-0.189, -0.101); F (30%missing data)=-0.105, 95% CI (-0.199, -0.011)]. However, the reconstruction of demographic history of M. catharinensis indicated that the population should be suffered a recent bottleneck. Our population genomic study tackles a central issue in evolution and conservation biology and we expect that it will be useful to help safeguard the remaining genetic diversity reported for this unique genetic resource.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Genetic diversity in a long-lived mammal is explained by the past's demographic shadow and current connectivity

Within-species genetic diversity is crucial for the persistence and integrity of populations and ecosystems. Conservation actions require an understanding of factors influencing genetic diversity, especially in the context of global change. Both population size and connectivity are factors greatly influencing genetic diversity; the relative importance of these factors can, however, change through time. Hence, quantifying the degree to which population size or genetic connectivity are shaping genetic diversity, and at which ecological time scale (past or present), is challenging, yet essential for the development of efficient conservation strategies. In this study, we estimated the genetic diversity of 42 colonies of Rhinolophus hipposideros, a long-lived mammal vulnerable to global change, sampling locations spanning its continental northern range. Here, we present an integrative approach that disentangles and quantifies the contribution of different connectivity measures in addition to contemporary colony size and historic bottlenecks in shaping genetic diversity. In our study, the best model explained 64% of the variation in genetic diversity. It included historic bottlenecks, contemporary colony size, connectivity and a negative interaction between the latter two. Contemporary connectivity explained most genetic diversity when considering a 65 km radius around the focal colonies, emphasizing the large geographic scale at which the positive impact of connectivity on genetic diversity is most profound and hence, the minimum scale at which conservation should be planned. Our results highlight that the relative importance of the two main factors shaping genetic diversity varies through time, emphasizing the relevance of disentangling them to ensure appropriate conservation strategies.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (P_ST-F_ST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. P_ST-F_ST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

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.

Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog

Temperature is a critical driver of ectotherm life‐history strategies, whereby a warmer environment is associated with increased growth, reduced longevity and accelerated senescence. Increasing evidence indicates that thermal adaptation may underlie such life‐history shifts in wild populations. Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) can help uncover the molecular mechanisms of temperature‐driven variation in growth, longevity and senescence. However, our understanding of these mechanisms is still limited, which reduces our ability to predict the response of non‐model ectotherms to global temperature change.

In this study, we examined the potential role of thermal adaptation in clinal shifts of life‐history traits (i.e. life span, senescence rate and recruitment) in the Columbia spotted frog Rana luteiventris along a broad temperature gradient in the western United States.

We took advantage of extensive capture–recapture datasets of 20,033 marked individuals from eight populations surveyed annually for 14–18 years to examine how mean annual temperature and precipitation influenced demographic parameters (i.e. adult survival, life span, senescence rate, recruitment and population growth). After showing that temperature was the main climatic predictor influencing demography, we used RAD‐seq data (50,829 SNPs and 6,599 putative CNVs) generated for 352 individuals from 31 breeding sites to identify the genomic signatures of thermal adaptation.

Our results showed that temperature was negatively associated with annual adult survival and reproductive life span and positively associated with senescence rate. By contrast, recruitment increased with temperature, promoting the long‐term viability of most populations. These temperature‐dependent demographic changes were associated with strong genomic signatures of thermal adaptation. We identified 148 SNP candidates associated with temperature including three SNPs located within protein‐coding genes regulating resistance to cold and hypoxia, immunity and reproduction in ranids. We also identified 39 CNV candidates (including within 38 transposable elements) for which normalized read depth was associated with temperature.

Our study indicates that both SNPs and structural variants are associated with temperature and could eventually be found to play a functional role in clinal shifts in senescence rate and life‐history strategies in R. luteiventris. These results highlight the potential role of different sources of molecular variation in the response of ectotherms to environmental temperature variation in the context of global warming.

Limited transgenerational effects of environmental temperatures on thermal performance of a cold-adapted salmonid

The capacity of ectotherms to cope with rising temperatures associated with climate change is a significant conservation concern as the rate of warming is likely too rapid to allow for adaptative responses in many populations. Transgenerational plasticity (TGP), if present, could potentially buffer some of the negative impacts of warming on future generations. We examined TGP in lake trout to assess their inter-generational potential to cope with anticipated warming. We acclimated adult lake trout to cold (10°C) or warm (17°C) temperatures for several months, then bred them to produce offspring from parents within a temperature treatment (cold-acclimated and warm-acclimated parents) and between temperature treatments (i.e. reciprocal crosses). At the fry stage, offspring were also acclimated to cold (11°C) or warm (15°C) temperatures. Thermal performance was assessed by measuring their critical thermal maximum (CTM) and the change in metabolic rate during an acute temperature challenge. From this dataset, we also determined their resting and peak (highest achieved, thermally induced) metabolic rates. There was little variation in offspring CTM or peak metabolic rate, although cold-acclimated offspring from warm-acclimated parents exhibited elevated resting metabolic rates without a corresponding increase in mass or condition factor, suggesting that transgenerational effects can be detrimental when parent and offspring environments mismatch. These results suggest that the limited TGP in thermal performance of lake trout is unlikely to significantly influence population responses to projected increases in environmental temperatures.

The ghosts of propagation past: haplotype information clarifies the relative influence of stocking history and phylogeographic processes on contemporary population structure of walleye (Sander vitreus)

Stocking of fish is an important tool for maintaining fisheries but can also significantly alter population genetic structure and erode the portfolio of within-species diversity that is important for promoting resilience and adaptability. Walleye (Sander vitreus) are a highly valued sportfish in the midwestern United States, a region characterized by postglacial recolonization from multiple lineages and an extensive history of stocking. We leveraged genomic data and recently developed analytical approaches to explore the population structure of walleye from two midwestern states, Minnesota and Wisconsin. We genotyped 954 walleye from 23 populations at ~20,000 loci using genotyping by sequencing and tested for patterns of population structure with single-SNP and microhaplotype data. Populations from Minnesota and Wisconsin were highly differentiated from each other, with additional substructure found in each state. Population structure did not consistently adhere to drainage boundaries, as cases of high intra-drainage and low inter-drainage differentiation were observed. Low genetic structure was observed between populations from the upper Wisconsin and upper Chippewa river watersheds, which are found as few as 50 km apart and were likely homogenized through historical stocking. Nevertheless, we were able to differentiate these populations using microhaplotype-based co-ancestry analysis, providing increased resolution over previous microsatellite studies and our other single SNP-based analyses. Although our results illustrate that walleye population structure has been influenced by past stocking practices, native ancestry still exists in most populations and walleye populations may be able to purge non-native alleles and haplotypes in the absence of stocking. Our study is one of the first to use genomic tools to investigate the influence of stocking on population structure in a nonsalmonid fish and outlines a workflow leveraging recently developed analytical methods to improve resolution of complex population structure that will be highly applicable in many species and systems.

Genomic signatures of drift and selection driven by predation and human pressure in an insular lizard

Genomic divergence was studied in 10 small insular populations of the endangered Balearic Islands lizard (Podarcis lilfordi) using double digest restriction-site associated DNA sequencing. The objectives were to establish levels of divergence among populations, investigate the impact of population size on genetic variability and to evaluate the role of different environmental factors on local adaptation. Analyses of 72,846 SNPs supported a highly differentiated genetic structure, being the populations with the lowest population size (Porros, Foradada and Esclatasang islets) the most divergent, indicative of greater genetic drift. Outlier tests identified ~ 2% of loci as candidates for selection. Genomic divergence-Enviroment Association analyses were performed using redundancy analyses based on SNPs putatively under selection, detecting predation and human pressure as the environmental variables with the greatest explanatory power. Geographical distributions of populations and environmental factors appear to be fundamental drivers of divergence. These results support the combined role of genetic drift and divergent selection in shaping the genetic structure of these endemic island lizard populations.

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

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

Demographic history and genomics of local adaptation in blue tit populations

Understanding the genomic processes underlying local adaptation is a central aim of modern evolutionary biology. This task requires identifying footprints of local selection but also estimating spatio‐temporal variations in population demography and variations in recombination rate and in diversity along the genome. Here, we investigated these parameters in blue tit populations inhabiting deciduous versus evergreen forests, and insular versus mainland areas, in the context of a previously described strong phenotypic differentiation. Neighboring population pairs of deciduous and evergreen habitats were weakly genetically differentiated (F_ST = 0.003 on average), nevertheless with a statistically significant effect of habitat type on the overall genetic structure. This low differentiation was consistent with the strong and long‐lasting gene flow between populations inferred by demographic modeling. In turn, insular and mainland populations were moderately differentiated (F_ST = 0.08 on average), in line with the inference of moderate ancestral migration, followed by isolation since the end of the last glaciation. Effective population sizes were large, yet smaller on the island than on the mainland. Weak and nonparallel footprints of divergent selection between deciduous and evergreen populations were consistent with their high connectivity and the probable polygenic nature of local adaptation in these habitats. In turn, stronger footprints of divergent selection were identified between long isolated insular versus mainland birds and were more often found in regions of low recombination, as expected from theory. Lastly, we identified a genomic inversion on the mainland, spanning 2.8 Mb. These results provide insights into the demographic history and genetic architecture of local adaptation in blue tit populations at multiple geographic scales.

Small coastal streams-Critical reservoirs of genetic diversity for trout (Salmo trutta L.) in the face of increasing anthropogenic stressors

We used microsatellite markers to investigate levels and structuring of genetic diversity in trout (Salmo trutta L.) sampled from 16 rivers along the south coast of Cornwall in southwest England. This region is characterized by many small coastal streams with a few larger catchments. At a regional level, genetic structuring of contemporary populations has been influenced by a combination of events, including the last Ice Age and also more recent human activities over the last millennium. All populations are shown to have gone through strong genetic bottlenecks, coinciding with increased exploitation of mineral resources within catchments, beginning during the Medieval period. At more local levels, contemporary human-induced habitat fragmentation, such as weir and culvert construction, has disproportionally affected trout populations in the smaller catchments within the study area. However, where small catchments are relatively unaffected by such activities, they can host trout populations with diversity levels comparable to those found in larger rivers in the region. We also predict significant future loses of diversity and heterozygosity in the trout populations inhabiting small, isolated catchments. Our study highlights how multiple factors, especially the activity of humans, have and continue to affect the levels and structuring of genetic diversity in trout over long timescales.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Conservation and Management of Salmon in the Age of Genomics

Salmon were among the first nonmodel species for which systematic population genetic studies of natural populations were conducted, often to support management and conservation. The genomics revolution has improved our understanding of the evolutionary ecology of salmon in two major ways: (a) Large increases in the numbers of genetic markers (from dozens to 10⁴-10⁶) provide greater power for traditional analyses, such as the delineation of population structure, hybridization, and population assignment, and (b) qualitatively new insights that were not possible with traditional genetic methods can be achieved by leveraging detailed information about the structure and function of the genome. Studies of the first type have been more common to date, largely because it has taken time for the necessary tools to be developed to fully understand the complex salmon genome. We expect that the next decade will witness many new studies that take full advantage of salmonid genomic resources.

No evidence of inbreeding depression in fast declining herds of migratory caribou

Identifying inbreeding depression early in small and declining populations is essential for management and conservation decisions. Correlations between heterozygosity and fitness (HFCs) provide a way to identify inbreeding depression without prior knowledge of kinship among individuals. In Northern Quebec and Labrador, the size of two herds of migratory caribou (Rivière-George, RG and Rivière-aux-Feuilles, RAF) has declined by one to two orders of magnitude in the last three decades. This raises the question of a possible increase in inbreeding depression originating from, and possibly contributing to, the demographic decline in those populations. Here, we tested for the association of genomic inbreeding indices (estimated with 22,073 SNPs) with body mass and survival in 400 caribou sampled in RG and RAF herds between 1996 and 2016. We found no association of individual heterozygosity or inbreeding coefficient with body mass or annual survival. Furthermore, those genomic inbreeding indices remained stable over the period monitored. These results suggest that the rapid and intense demographic decline of the herds did not cause inbreeding depression in those populations. Although we found no evidence for HFCs, if demographic decline continues, it is possible that such inbreeding depression would be triggered.

Conservation through the lens of (mal)adaptation: Concepts and meta-analysis

Evolutionary approaches are gaining popularity in conservation science, with diverse strategies applied in efforts to support adaptive population outcomes. Yet conservation strategies differ in the type of adaptive outcomes they promote as conservation goals. For instance, strategies based on genetic or demographic rescue implicitly target adaptive population states whereas strategies utilizing transgenerational plasticity or evolutionary rescue implicitly target adaptive processes. These two goals are somewhat polar: adaptive state strategies optimize current population fitness, which should reduce phenotypic and/or genetic variance, reducing adaptability in changing or uncertain environments; adaptive process strategies increase genetic variance, causing maladaptation in the short term, but increase adaptability over the long term. Maladaptation refers to suboptimal population fitness, adaptation refers to optimal population fitness, and (mal)adaptation refers to the continuum of fitness variation from maladaptation to adaptation. Here, we present a conceptual classification for conservation that implicitly considers (mal)adaptation in the short-term and long-term outcomes of conservation strategies. We describe cases of how (mal)adaptation is implicated in traditional conservation strategies, as well as strategies that have potential as a conservation tool but are relatively underutilized. We use a meta-analysis of a small number of available studies to evaluate whether the different conservation strategies employed are better suited toward increasing population fitness across multiple generations. We found weakly increasing adaptation over time for transgenerational plasticity, genetic rescue, and evolutionary rescue. Demographic rescue was generally maladaptive, both immediately after conservation intervention and after several generations. Interspecific hybridization was adaptive only in the F1 generation, but then rapidly leads to maladaptation. Management decisions that are made to support the process of adaptation must adequately account for (mal)adaptation as a potential outcome and even as a tool to bolster adaptive capacity to changing conditions.

Genetic signatures of small effective population sizes and demographic declines in an endangered rattlesnake, Sistrurus catenatus

Endangered species that exist in small isolated populations are at elevated risk of losing adaptive variation due to genetic drift. Analyses that estimate short-term effective population sizes, characterize historical demographic processes, and project the trajectory of genetic variation into the future are useful for predicting how levels of genetic diversity may change. Here, we use data from two independent types of genetic markers (single nucleotide polymorphisms [SNPs] and microsatellites) to evaluate genetic diversity in 17 populations spanning the geographic range of the endangered eastern massasauga rattlesnake (Sistrurus catenatus). First, we use SNP data to confirm previous reports that these populations exhibit high levels of genetic structure (overall Fst = 0.25). Second, we show that most populations have contemporary Ne estimates <50. Heterozygosity-fitness correlations in these populations provided no evidence for a genetic cost to living in small populations, though these tests may lack power. Third, model-based demographic analyses of individual populations indicate that all have experienced declines, with the onset of many of these declines occurring over timescales consistent with anthropogenic impacts (<200 years). Finally, forward simulations of the expected loss of variation in relatively large (Ne = 50) and small (Ne = 10) populations indicate they will lose a substantial amount of their current standing neutral variation (63% and 99%, respectively) over the next 100 years. Our results argue that drift has a significant and increasing impact on levels of genetic variation in isolated populations of this snake, and efforts to assess and mitigate associated impacts on adaptive variation should be components of the management of this endangered reptile.

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.

Next-generation conservation genetics and biodiversity monitoring

This special issue of Evolutionary Applications consists of 10 publications investigating the use of next-generation tools and techniques in population genetic analyses and biodiversity assessment. The special issue stems from a 2016 Next Generation Genetic Monitoring Workshop, hosted by the National Institute for Mathematical and Biological Synthesis (NIMBioS) in Tennessee, USA. The improved accessibility of next-generation sequencing platforms has allowed molecular ecologists to rapidly produce large amounts of data. However, with the increased availability of new genomic markers and mathematical techniques, care is needed in selecting appropriate study designs, interpreting results in light of conservation concerns, and determining appropriate management actions. This special issue identifies key attributes of successful genetic data analyses in biodiversity evaluation and suggests ways to improve analyses and their application in current population and conservation genetics research.