Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Genetic structure and historic demography of endangered unarmoured threespine stickleback at southern latitudes signals a potential new management approach

Habitat loss, flood control infrastructure, and drought have left most of southern California and northern Baja California's native freshwater fish near extinction, including the endangered unarmoured threespine stickleback (Gasterosteus aculeatus williamsoni). This subspecies, an unusual morph lacking the typical lateral bony plates of the G. aculeatus complex, occurs at arid southern latitudes in the eastern Pacific Ocean and survives in only three inland locations. Managers have lacked molecular data to answer basic questions about the ancestry and genetic distinctiveness of unarmoured populations. These data could be used to prioritize conservation efforts. We sampled G. aculeatus from 36 localities and used microsatellites and whole genome data to place unarmoured populations within the broader evolutionary context of G. aculeatus across southern California/northern Baja California. We identified three genetic groups with none consisting solely of unarmoured populations. Unlike G. aculeatus at northern latitudes, where Pleistocene glaciation has produced similar historical demographic profiles across populations, we found markedly different demographics depending on sampling location, with inland unarmoured populations showing steeper population declines and lower heterozygosity compared to low armoured populations in coastal lagoons. One exception involved the only high elevation population in the region, where the demography and alleles of unarmoured fish were similar to low armoured populations near the coast, exposing one of several cases of artificial translocation. Our results suggest that the current "management-by-phenotype" approach, based on lateral plates, is incidentally protecting the most imperilled populations; however, redirecting efforts toward evolutionary units, regardless of phenotype, may more effectively preserve adaptive potential.

Characterizing phenotypic diversity in marine populations of the threespine stickleback

The threespine stickleback (Gasterosteus aculeatus) is an important model for studying the evolution of vertebrate morphology. Sticklebacks inhabit freshwater, brackish, and marine northern hemisphere waters. Anadromous and marine populations (hereafter marine) are assumed to have remained unchanged morphologically from ancestral marine sticklebacks, despite marine environments varying on regional and local scales. Recent studies suggest that genetic and phenotypic structure exists in marine populations, yet the scale of this variation, and its ecological causes remain unclear. Our goal was to assess morphological trait variation in marine stickleback populations around Southern British Columbia (BC) and determine if oceanographic and habitat characteristics were associated with this variation. Between May-July 2019, we sampled 534 sticklebacks from 15 sites around Vancouver Island, a region characterized by a large diversity of oceanographic and habitat features. We characterized trait variation using two-dimensional (2D) geometric morphometric analysis, comparing individuals between oceanographic regions and habitats. We focused on head and body shape. We found that marine sticklebacks varied morphologically among and between regions and habitats, but the variation did not appear to be related to environmental variation. Sexual dimorphism was the largest source of variation, but oceanographic and habitat variables influenced differences between sexes. We concluded that marine sticklebacks offer abundant opportunities for expanding our knowledge of drivers of morphology.

Sexually mediated phenotypic variation within and between sexes as a continuum structured by ecology: The mosaic nature of skeletal variation across body regions in Threespine stickleback (Gasterosteus aculeatus L.)

Ecological character displacement between the sexes, and sexual selection, integrate into a convergent set of factors that produce sexual variation. Ecologically modulated, sexually mediated variation within and between sexes may be a major contributor to the amount of total variation that selection can act on in species. Threespine stickleback (Gasterosteus aculeatus) display rapid adaptive responses and sexual variation in many phenotypic traits. We examined phenotypic variation in the skull, pectoral and pelvic girdles of threespine stickleback from two freshwater and two coastal marine sites on the Sunshine Coast of British Columbia, Canada, using an approach that avoids a priori assumptions about bimodal patterns of variation. We quantified shape and size of the cranial, pectoral and pelvic regions of sticklebacks in marine and freshwater habitats using 3D geometric morphometrics and an index of sexually mediated variation. We show that the expression of phenotypic variation is structured in part by the effects of both habitat marine vs freshwater and the effects of individual sites within each habitat. Relative size exerts variable influence, and patterns of phenotypic variation associated with sex vary among body regions. This fine-grained quantification of sexually mediated variation in the context of habitat difference and different anatomical structures indicates a complex relationship between genetically inferred sex and environmental factors, demonstrating that the interplay between shared genetic background and sexually mediated, ecologically based selective pressures structures the phenotypic expression of complex traits.

Genetic divergence among threespine stickleback that differ in nuptial coloration

Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

Evolution of thermal physiology alters the projected range of threespine stickleback under climate change

Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of its range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographical range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data were either static ("no evolution" models), allowed to evolve at observed evolutionary rates ("evolution" models) or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL; "scaled evolution" models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.

Adaptive divergence and the evolution of hybrid trait mismatch in threespine stickleback

Selection against mismatched traits in hybrids is the phenotypic analogue of intrinsic hybrid incompatibilities. Mismatch occurs when hybrids resemble one parent population for some phenotypic traits and the other parent population for other traits, and is caused by dominance in opposing directions or from segregation of alleles in recombinant hybrids. In this study, we used threespine stickleback fish (Gasterosteus aculeatus L.) to test the theoretical prediction that trait mismatch in hybrids should increase with the magnitude of phenotypic divergence between parent populations. We measured morphological traits in parents and hybrids in crosses between a marine population representing the ancestral form and twelve freshwater populations that have diverged from this ancestral state to varying degrees according to their environments. We found that trait mismatch was greater in more divergent crosses for both F₁ and F₂ hybrids. In the F₁, the divergence–mismatch relationship was caused by traits having dominance in different directions, whereas it was caused by increasing segregating phenotypic variation in the F₂. Our results imply that extrinsic hybrid incompatibilities accumulate as phenotypic divergence proceeds.

Comparing genome scans among species of the stickleback order reveals three different patterns of genetic diversity

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus), and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species' range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (F ST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

Meta-analysis: Congruence of genomic and phenotypic differentiation across diverse natural study systems

Linking genotype to phenotype is a primary goal for understanding the genomic underpinnings of evolution. However, little work has explored whether patterns of linked genomic and phenotypic differentiation are congruent across natural study systems and traits. Here, we investigate such patterns with a meta-analysis of studies examining population-level differentiation at subsets of loci and traits putatively responding to divergent selection. We show that across the 31 studies (88 natural population-level comparisons) we examined, there was a moderate (R 2 = 0.39) relationship between genomic differentiation (F ST ) and phenotypic differentiation (P ST ) for loci and traits putatively under selection. This quantitative relationship between P ST and F ST for loci under selection in diverse taxa provides broad context and cross-system predictions for genomic and phenotypic adaptation by natural selection in natural populations. This context may eventually allow for more precise ideas of what constitutes "strong" differentiation, predictions about the effect size of loci, comparisons of taxa evolving in nonparallel ways, and more. On the other hand, links between P ST and F ST within studies were very weak, suggesting that much work remains in linking genomic differentiation to phenotypic differentiation at specific phenotypes. We suggest that linking genotypes to specific phenotypes can be improved by correlating genomic and phenotypic differentiation across a spectrum of diverging populations within a taxon and including wide coverage of both genomes and phenomes.

Ancient and modern stickleback genomes reveal the demographic constraints on adaptation

Adaptation is typically studied by comparing modern populations with contrasting environments. Individuals persisting in the ancestral habitat are typically used to represent the ancestral founding population; however, it has been questioned whether these individuals are good proxies for the actual ancestors.¹ To address this, we applied a paleogenomics approach² to directly access the ancestral genepool: partially sequencing the genomes of two 11- to 13,000-year-old stickleback recovered from the transitionary layer between marine and freshwater sediments of two Norwegian isolation lakes³ and comparing them with 30 modern stickleback genomes from the same lakes and adjacent marine fjord, in addition to a global dataset of 20 genomes.⁴ The ancient stickleback shared genome-wide ancestry with the modern fjord population, whereas modern lake populations have lost substantial ancestral variation following founder effects, and subsequent drift and selection. Freshwater-adaptive alleles found in one ancient stickleback genome have not risen to high frequency in the present-day population from the same lake. Comparison to the global dataset suggested incomplete adaptation to freshwater in our modern lake populations. Our findings reveal the impact of population bottlenecks in constraining adaptation due to reduced efficacy of selection on standing variation present in founder populations.

Evolution: Predictability and the promise of ancient DNA

Is evolution predictable? Genomes from thousands-of-years-old stickleback suggest that, despite substantial stochasticity in the course of evolution, our understanding of the recent evolutionary past of this species was generally valid.

Mitochondrial volume density and evidence for its role in adaptive divergence in response to thermal tolerance in threespine stickleback

Phenotypic plasticity is predicted to permit persistence in new environments, and may subsequently evolve to enhance fitness. Colonizing environments with lower winter temperatures can lead to the evolution of lower critical thermal minima; the corresponding physiological traits associated with temperature tolerance are predicted to involve mitochondrial function. Threespine stickleback (Gasterosteus aculeatus) have colonized freshwater lakes along the Pacific Northwest. These freshwater populations are known to exhibit cold-induced increases in mitochondrial volume density in pectoral muscle, but whether such plasticity evolved before or after colonization is uncertain. Here, we measure critical thermal minima (CTmin) in one marine and one freshwater population of threespine stickleback, and mitochondrial volume density in pectoral and cardiac tissue of both populations acclimated to different temperature treatments (6.2, 14.5 and 20.6 ℃). Mitochondrial volume density increased with cold acclimation in pectoral muscle; cardiac muscle was non-plastic but had elevated mitochondrial volume densities compared to pectoral muscle across all temperature treatments. There were no differences in the levels of plasticity between marine and freshwater stickleback, but neither were there differences in CTmin. Importantly, marine stickleback exhibited plasticity under low-salinity conditions, suggesting that marine stickleback had at least one necessary phenotype for persistence in freshwater environments before colonization occurred.

A three-pronged approach that leans on Indigenous knowledge for northern fish monitoring and conservation

Investigating whether changes within fish populations may result from harvesting requires a comprehensive approach, especially in more data-sparse northern regions. Our study took a three-pronged approach to investigate walleye population change by combining Indigenous knowledge (IK), phenotypic traits, and genomics. We thank Larson et al. (2020) for their critiques of our study; certainly, there are aspects of their critique that are warranted and merit further investigation. However, we argue that their critique is over-stated and misleading, primarily given that (a) one of three prongs of our research, IK, was dismissed in their assessment of our study's conclusions; (b) our Bayesian size-at-age modeling should help to mitigate sample size issues; (c) their re-analysis of our size-at-age data does not actually refute our results; (d) genomic changes that we observed are nascent; (e) the data file that Larson et al. (2020) used for their genomic re-analysis was not correct; and (f) criteria that Larson et al. (2020) use for their genomic re-analysis were not properly justified.

Threespine Stickleback in Lake Constance: The Ecology and Genomic Substrate of a Recent Invasion

Invasive species can be powerful models for studying contemporary evolution in natural environments. As invading organisms often encounter new habitats during colonization, they will experience novel selection pressures. Threespine stickleback (Gasterosteus aculeatus complex) have recently colonized large parts of Switzerland and are invasive in Lake Constance. Introduced to several watersheds roughly 150 years ago, they spread across the Swiss Plateau (400–800 m a.s.l.), bringing three divergent hitherto allopatric lineages into secondary contact. As stickleback have colonized a variety of different habitat types during this recent range expansion, the Swiss system is a useful model for studying contemporary evolution with and without secondary contact. For example, in the Lake Constance region there has been rapid phenotypic and genetic divergence between a lake population and some stream populations. There is considerable phenotypic variation within the lake population, with individuals foraging in and occupying littoral, offshore pelagic, and profundal waters, the latter of which is a very unusual habitat for stickleback. Furthermore, adults from the lake population can reach up to three times the size of adults from the surrounding stream populations, and are large by comparison to populations globally. Here, we review the historical origins of the threespine stickleback in Switzerland, and the ecomorphological variation and genomic basis of its invasion in Lake Constance. We also outline the potential ecological impacts of this invasion, and highlight the interest for contemporary evolution studies.

Diversity and sexual dimorphism in the head lateral line system in North Sea populations of threespine sticklebacks, Gasterosteus aculeatus (Teleostei: Gasterosteidae)

The mechanosensory lateral line of fishes is a flow sensing system and supports a number of behaviors, e.g. prey detection, schooling or position holding in water currents. Differences in the neuromast pattern of this sensory system reflect adaptation to divergent ecological constraints. The threespine stickleback, Gasterosteus aculeatus, is known for its ecological plasticity resulting in three major ecotypes, a marine type, a migrating anadromous type and a resident freshwater type. We provide the first comparative study of the pattern of the head lateral line system of North Sea populations representing these three ecotypes including a brackish spawning population. We found no distinct difference in the pattern of the head lateral line system between the three ecotypes but significant differences in neuromast numbers. The anadromous and the brackish populations had distinctly less neuromasts than their freshwater and marine conspecifics. This difference in neuromast number between marine and anadromous threespine stickleback points to differences in swimming behavior. We also found sexual dimorphism in neuromast number with males having more neuromasts than females in the anadromous, brackish and the freshwater populations. But no such dimorphism occurred in the marine population. Our results suggest that the head lateral line of the three ecotypes is under divergent hydrodynamic constraints. Additionally, sexual dimorphism points to divergent niche partitioning of males and females in the anadromous and freshwater but not in the marine populations. Our findings imply careful sampling as an important prerequisite to discern especially between anadromous and marine threespine sticklebacks.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.