Implications of Large-Effect Loci for Conservation: A Review and Case Study with Pacific Salmon

Temporal allele frequency changes in large-effect loci reveal potential fishing impacts on salmon life-history diversity

Fishing has the potential to influence the life-history traits of exploited populations. However, our understanding of how fisheries can induce evolutionary genetic changes remains incomplete. The discovery of large-effect loci linked with ecologically important life-history traits, such as age at maturity in Atlantic salmon (Salmo salar), provides an opportunity to study the impacts of temporally varying fishing pressures on these traits. A 93-year archive of fish scales from wild Atlantic salmon catches from the northern Baltic Sea region allowed us to monitor variation in adaptive genetic diversity linked with age at maturity of wild Atlantic salmon populations. The dataset consisted of samples from both commercial and recreational fisheries that target salmon on their spawning migration. Using a genotyping-by-sequencing approach (GT-seq), we discovered strong within-season allele frequency changes at the vgll3 locus linked with Atlantic salmon age at maturity: fishing in the early season preferentially targeted the vgll3 variant linked with older maturation. We also found within-season temporal variation in catch proportions of different wild Atlantic salmon subpopulations. Therefore, selective pressures of harvesting may vary depending on the seasonal timing of fishing, which has the potential to cause evolutionary changes in key life-history traits and their diversity. This knowledge can be used to guide fisheries management to reduce the effects of fishing practices on salmon life-history diversity. Thus, this study provides a tangible example of using genomic approaches to infer, monitor and help mitigate human impacts on adaptively important genetic variation in nature.

Genetic parentage reveals the (un)natural history of Central Valley hatchery steelhead

Populations composed of individuals descended from multiple distinct genetic lineages often feature significant differences in phenotypic frequencies. We considered hatchery production of steelhead, the migratory anadromous form of the salmonid species Oncorhynchus mykiss, and investigated how differences among genetic lineages and environmental variation impacted life history traits. We genotyped 23,670 steelhead returning to the four California Central Valley hatcheries over 9 years from 2011 to 2019, confidently assigning parentage to 13,576 individuals to determine age and date of spawning and rates of iteroparity and repeat spawning within each year. We found steelhead from different genetic lineages showed significant differences in adult life history traits despite inhabiting similar environments. Differences between coastal and Central Valley steelhead lineages contributed to significant differences in age at return, timing of spawning, and rates of iteroparity among programs. In addition, adaptive genomic variation associated with life history development in this species varied among hatchery programs and was associated with the age of steelhead spawners only in the coastal lineage population. Environmental variation likely contributed to variations in phenotypic patterns observed over time, as our study period spanned both a marine heatwave and a serious drought in California. Our results highlight evidence of a strong genetic component underlying known phenotypic differences in life history traits between two steelhead lineages.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

A landmark environmental law looks ahead

In late December 1973, the United States enacted what some would come to call "the pitbull of environmental laws." In the 50 years since, the formidable regulatory teeth of the Endangered Species Act (ESA) have been credited with considerable successes, obliging agencies to draw upon the best available science to protect species and habitats. Yet human pressures continue to push the planet toward extinctions on a massive scale. With that prospect looming, and with scientific understanding ever changing, Science invited experts to discuss how the ESA has evolved and what its future might hold. -Brad Wible.

Where to now with the evolutionarily significant unit?

The designation of units for conservation has been a necessary but challenging objective since conservation efforts began. Most species are divided, typically by environment, into populations with independent evolutionary trajectories. There are practical conservation objectives for defining these boundaries. Separate genetic clusters provide future evolutionary potential as environments change, and individuals in isolated populations may lose fitness when population size is reduced. The history of the effort to define units is briefly reviewed here, but I focus on finding a process that may facilitate uniform and effective application at a time when conservation urgency is great. I propose a refinement of the designated unit concept, distinguishing between conservation units (CUs) and evolutionarily sustaining conservation units (ESCUs).

Multifaceted framework for defining conservation units: An example from Atlantic salmon (Salmo salar) in Canada

Conservation units represent important components of intraspecific diversity that can aid in prioritizing and protecting at-risk populations, while also safeguarding unique diversity that can contribute to species resilience. In Canada, identification and assessments of conservation units is done by the Committee on the Status of Endangered Wildlife in Canada (COSEWIC). COSEWIC can recognize conservation units below the species level (termed "designatable units"; DUs) if the unit has attributes that make it both discrete and evolutionarily significant. There are various ways in which a DU can meet criteria of discreteness and significance, and increasing access to "big data" is providing unprecedented information that can directly inform both criteria. Specifically, the incorporation of genomic data for an increasing number of non-model species is informing more COSEWIC assessments; thus, a repeatable, robust framework is needed for integrating these data into DU characterization. Here, we develop a framework that uses a multifaceted, weight of evidence approach to incorporate multiple data types, including genetic and genomic data, to inform COSEWIC DUs. We apply this framework to delineate DUs of Atlantic salmon (Salmo salar, L.), an economically, culturally, and ecologically significant species, that is also characterized by complex hierarchical population structure. Specifically, we focus on an in-depth example of how our approach was applied to a previously data limited region of northern Canada that was defined by a single large DU. Application of our framework with newly available genetic and genomic data led to subdividing this DU into three new DUs. Although our approach was developed to meet criteria of COSEWIC, it is widely applicable given similarities in the definitions of a conservation unit.

Detecting parallel polygenic adaptation to novel evolutionary pressure in wild populations: a case study in Atlantic cod (Gadus morhua)

Populations can adapt to novel selection pressures through dramatic frequency changes in a few genes of large effect or subtle shifts in many genes of small effect. The latter (polygenic adaptation) is expected to be the primary mode of evolution for many life-history traits but tends to be more difficult to detect than changes in genes of large effect. Atlantic cod (Gadus morhua) were subjected to intense fishing pressure over the twentieth century, leading to abundance crashes and a phenotypic shift toward earlier maturation across many populations. Here, we use spatially replicated temporal genomic data to test for a shared polygenic adaptive response to fishing using methods previously applied to evolve-and-resequence experiments. Cod populations on either side of the Atlantic show covariance in allele frequency change across the genome that are characteristic of recent polygenic adaptation. Using simulations, we demonstrate that the degree of covariance in allele frequency change observed in cod is unlikely to be explained by neutral processes or background selection. As human pressures on wild populations continue to increase, understanding and attributing modes of adaptation using methods similar to those demonstrated here will be important in identifying the capacity for adaptive responses and evolutionary rescue. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Whole genome resequencing identifies local adaptation associated with environmental variation for redband trout

Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low-coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L, ROCK1, SIX6). Genotype-environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high-elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high-elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

Assortative mating for reproductive timing affects population recruitment and resilience in a quantitative genetic model

Quantitative models that simulate the inheritance and evolution of fitness-linked traits offer a method for predicting how environmental or anthropogenic perturbations can affect the dynamics of wild populations. Random mating between individuals within populations is a key assumption of many such models used in conservation and management to predict the impacts of proposed management or conservation actions. However, recent evidence suggests that non-random mating may be underestimated in wild populations and play an important role in diversity-stability relationships. Here we introduce a novel individual-based quantitative genetic model that incorporates assortative mating for reproductive timing, a defining attribute of many aggregate breeding species. We demonstrate the utility of this framework by simulating a generalized salmonid lifecycle, varying input parameters, and comparing model outputs to theoretical expectations for several eco-evolutionary, population dynamic scenarios. Simulations with assortative mating systems resulted in more resilient and productive populations than those that were randomly mating. In accordance with established ecological and evolutionary theory, we also found that decreasing the magnitude of trait correlations, environmental variability, and strength of selection each had a positive effect on population growth. Our model is constructed in a modular framework so that future components can be easily added to address pressing issues such as the effects of supportive breeding, variable age structure, differential selection by sex or age, and fishery interactions on population growth and resilience. With code published in a public Github repository, model outputs may easily be tailored to specific study systems by parameterizing with empirically generated values from long-term ecological monitoring programs.

Whole-genome resequencing reveals genetic differences and the genetic basis of parapodium number in Russian and Chinese Apostichopus japonicus

BACKGROUND

Apostichopus japonicus is an economically important species in the global aquaculture industry. Russian A. japonicus, mainly harvested in the Vladivostok region, exhibits significant phenotypic differentiation, including in many economically important traits, compared with Chinese A. japonicus owing to differences in their habitat. However, both the genetic basis for the phenotypic divergence and the population genetic structure of Russian and Chinese A. japonicus are unknown.

RESULT

In this study, 210 individuals from seven Russian and Chinese A. japonicus populations were sampled for whole-genome resequencing. The genetic structure analysis differentiated the Russian and Chinese A. japonicus into two groups. Population genetic analyses indicated that the Russian population showed a high degree of allelic linkage and had undergone stronger positive selection compared with the Chinese populations. Gene ontology terms enriched among candidate genes with group selection analysis were mainly involved in immunity, such as inflammatory response, antimicrobial peptides, humoral immunity, and apoptosis. Genome-wide association analysis yielded eight single-nucleotide polymorphism loci significantly associated with parapodium number, and these loci are located in regions with a high degree of genomic differentiation between the Chinese and Russia populations. These SNPs were associated with five genes. Gene expression validation revealed that three of these genes were significantly differentially expressed in individuals differing in parapodium number. AJAP08772 and AJAP08773 may directly affect parapodium production by promoting endothelial cell proliferation and metabolism, whereas AJAP07248 indirectly affects parapodium production by participating in immune responses.

CONCLUSIONS

This study, we performed population genetic structure and GWAS analysis on Chinese and Russian A. japonicus, and found three candidate genes related to the number of parapodium. The results provide an in-depth understanding of the differences in the genetic structure of A. japonicus populations in China and Russia, and provide important information for subsequent genetic analysis and breeding of this species.

Population structure and adaptive differentiation in the sea cucumber Apostichopus californicus and implications for spatial resource management

A growing body of evidence suggests that spatial population structure can develop in marine species despite large population sizes and high gene flow. Characterizing population structure is important for the effective management of exploited species, as it can be used to identify appropriate scales of management in fishery and aquaculture contexts. The California sea cucumber, Apostichopus californicus, is one such exploited species whose management could benefit from further characterization of population structure. Using restriction site-associated DNA (RAD) sequencing, we developed 2075 single nucleotide polymorphisms (SNPs) to quantify genetic structure over a broad section of the species' range along the North American west coast and within the Salish Sea, a region supporting the Washington State A. californicus fishery and developing aquaculture production of the species. We found evidence for population structure (global fixation index (FST) = 0.0068) with limited dispersal driving two patterns of differentiation: isolation-by-distance and a latitudinal gradient of differentiation. Notably, we found detectable population differences among collection sites within the Salish Sea (pairwise FST = 0.001-0.006). Using FST outlier detection and gene-environment association, we identified 10.2% of total SNPs as putatively adaptive. Environmental variables (e.g., temperature, salinity) from the sea surface were more correlated with genetic variation than those same variables measured near the benthos, suggesting that selection on pelagic larvae may drive adaptive differentiation to a greater degree than selection on adults. Our results were consistent with previous estimates of and patterns in population structure for this species in other extents of the range. Additionally, we found that patterns of neutral and adaptive differentiation co-varied, suggesting that adaptive barriers may limit dispersal. Our study provides guidance to decision-makers regarding the designation of management units for A. californicus and adds to the growing body of literature identifying genetic population differentiation in marine species despite large, nominally connected populations.

Understanding Local Adaptation to Prepare Populations for Climate Change

Adaptation within species to local environments is widespread in nature. Better understanding this local adaptation is critical to conserving biodiversity. However, conservation practices can rely on species’ trait averages or can broadly assume homogeneity across the range to inform management. Recent methodological advances for studying local adaptation provide the opportunity to fine-tune efforts for managing and conserving species. The implementation of these advances will allow us to better identify populations at greatest risk of decline because of climate change, as well as highlighting possible strategies for improving the likelihood of population persistence amid climate change. In the present article, we review recent advances in the study of local adaptation and highlight ways these tools can be applied in conservation efforts. Cutting-edge tools are available to help better identify and characterize local adaptation. Indeed, increased incorporation of local adaptation in management decisions may help meet the imminent demands of managing species amid a rapidly changing world.