Local adaptation to salinity in the three-spined stickleback?

Diving into broad-scale and high-resolution population genomics to decipher drivers of structure and climatic vulnerability in a marine invertebrate

Species with widespread distributions play a crucial role in our understanding of climate change impacts on population structure. In marine species, population structure is often governed by both high connectivity potential and selection across strong environmental gradients. Despite the complexity of factors influencing marine populations, studying species with broad distribution can provide valuable insights into the relative importance of these factors and the consequences of climate-induced alterations across environmental gradients. We used the northern shrimp Pandalus borealis and its wide latitudinal distribution to identify current drivers of population structure and predict the species' vulnerability to climate change. A total of 1514 individuals sampled across 24° latitude were genotyped at high geographic (54 stations) and genetic (14,331 SNPs) resolutions to assess genetic variation and environmental correlations. Four populations were identified in addition to finer substructure associated with local adaptation. Geographic patterns of neutral population structure reflected predominant oceanographic currents, while a significant proportion of the genetic variation was associated with gradients in salinity and temperature. Adaptive landscapes generated using climate projections suggest a larger genomic offset in the southern extent of the P. borealis range, where shrimp had the largest adaptive standing genetic variation. Our genomic results combined with recent observations point to further deterioration in southern regions and an impending vulnerable status in the regions at higher latitudes for P. borealis. They also provide rare insights into the drivers of population structure and climatic vulnerability of a widespread meroplanktonic species, which is crucial to understanding future challenges associated with invertebrates essential to ecosystem functioning.

A major 6 Mb superlocus is involved in pyrethroid resistance in the common bed bug Cimex lectularius

In the last few years, the bed bug Cimex lectularius has been an increasing problem worldwide, mainly due to the development of insecticide resistance to pyrethroids. The characterization of resistance alleles is a prerequisite to improve surveillance and resistance management. To identify genomic variants associated with pyrethroid resistance in Cimex lectularius, we compared the genetic composition of two recent and resistant populations with that of two ancient-susceptible strains using a genome-wide pool-seq design. We identified a large 6 Mb "superlocus" showing particularly high genetic differentiation and association with the resistance phenotype. This superlocus contained several clustered resistance genes and was also characterized by a high density of structural variants (inversions, duplications). The possibility that this superlocus constitutes a resistance "supergene" that evolved after the clustering of alleles adapted to insecticide and after reduction in recombination is discussed.

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.

Effects of salinity stress on methylation of the liver genome and complement gene in large yellow croaker (Larimichthys crocea)

Salinity is an important environmental factor that affects the yield and quality of large yellow croaker (Larimichthys crocea) during aquaculture. Here, whole-genome bisulfite sequencing (WGBS), RNA-seq, bisulfite sequencing PCR (BSP), quantitative real-time PCR (qPCR), and dual luciferase reporter gene detection technologies were used to analyze the DNA methylation characteristics and patterns of the liver genome, the expression and methylation levels of important immune genes in large yellow croaker in response to salinity stress. The results of WGBS showed that the cytosine methylation of CG type was dominant, CpGIsland and repeat regions were important regions where DNA methylation occurred, and the DNA methylation in upstream 2k (2000bp upstream of the promoter) and repeat regions had different changes in the liver tissue of large yellow croaker in the response to the 12‰, 24‰, 36‰ salinity stress of 4 w (weeks). In the combined analysis of WGBS and transcriptome, the complement and coagulation cascade pathways were significantly enriched, in which the complement-related genes C7, C3, C5, C4, C1R, MASP1, and CD59 were mainly changed in response to salinity stress. In the studied area of MASP1 gene promoter, the methylation levels of many CpG sites as well as total cytosine were strongly negatively correlated with mRNA expression level. Methylation function analysis of MASP1 promoter further proved that DNA methylation could inhibit the activity of MASP1 promoter, indicating that salinity may affect the expressions of complement-related genes by DNA methylation of gene promoter region.

Landscape genomics of the American lobster (Homarus americanus)

In marine species experiencing intense fishing pressures, knowledge of genetic structure and local adaptation represent a critical information to assist sustainable management. In this study, we performed a landscape genomics analysis in the American lobster to investigate the issues pertaining to the consequences of making use of putative adaptive loci to reliably infer population structure and thus more rigorously delineating biological management units in marine exploited species. Toward this end, we genotyped 14,893 single nucleotide polymorphism (SNPs) in 4190 lobsters sampled across 96 sampling sites distributed along 1000 km in the northwest Atlantic in both Canada and the USA. As typical for most marine species, we observed a weak, albeit highly significant genetic structure. We also found that adaptive genetic variation allows detecting fine-scale population structure not resolved by neutral genetic variation alone. Using the recent genome assembly of the American lobster, we were able to map and annotate several SNPs located in functional genes potentially implicated in adaptive processes such as thermal stress response, salinity tolerance and growth metabolism pathways. Taken together, our study indicates that weak population structure in high gene flow systems can be resolved at various spatial scales, and that putatively adaptive genetic variation can substantially enhance the delineation of biological management units of marine exploited species.

Freshwater adaptation in prickly sculpin (Pisces: Cottidae): intraspecific comparisons reveal evidence for water pH and Na+ concentration driving diversity in gill H+-ATPase and ionoregulation

Phenotypic divergence is a hallmark of adaptive radiation. One example involves differentiation in physiological traits involved in ion regulation among species with contrasting life-styles and living in distinct environments. Differentiation in ion regulation and its ecological implications among populations within species are, however, less well understood. To address this knowledge gap, we collected prickly sculpin (Cottus asper) from distinct habitat types including coastal rivers connected to estuaries, coastal lakes, and interior lakes, all from British Columbia, Canada. We tested for differences in plasma Na+ and Cl-, gill Na+/K+-ATPase and H+-ATPase activities and protein abundance as well as changes in body mass, and arterial blood pH in fish sampled from the field and acclimated to two different freshwater conditions in the laboratory including artificial lake water (ALW) and ion-poor water (IPW). We also tested for associations between environmental water chemistry and the physiological characteristics associated with ion regulation. Transfer to IPW resulted upregulation in gill Na+/K+-ATPase and H+-ATPase activities as well as increases in gill H+-ATPase protein expression level in each habitat compared with the common ALW treatment. Despite the presence of population-within-habitat type differences, significant habitat-type effects were revealed in most of the ion regulation characteristics examined under different acclimation conditions. Significantly lower plasma Cl- was detected in fish from coastal rivers compared to fish from the other two habitat types during the IPW treatment, which was also significantly lower compared with ALW. Similarly, gill Na+/K+-ATPase activity was lower in the coastal river populations in IPW than in fish from coastal and interior lakes, which was not in accordance with the protein expressed in the gill. For gill H+-ATPase, fish from interior lake populations had the highest level of activity across all habitat types under all conditions, which was related to the protein levels in the gill. The activity of gill H+-ATPase was positively correlated with the combined effect of water Na+ and pH under the ALW treatment. Our results suggest that variation in habitat may be an important factor driving differences in gill Na+/K+-ATPase and H+-ATPase activities across populations of C. asper. Further, the combined effect of water Na+ and pH may have played a key role in physiological adaptation in C. asper during post-glacial freshwater colonization and dispersal.

Metabolic costs associated with seawater acclimation in a euryhaline teleost, the fourspine stickleback (Apeltes quadracus)

The cost of osmoregulation in teleosts has been debated for decades, with estimates ranging from one to 30 % of routine metabolic rate. The variation in the energy budget appears to be greater for euryhaline fish due to their ability to withstand dynamic salinity levels. In this study, a time course of metabolic and physiological responses of the euryhaline fourspine stickleback (Apeltes quadracus) acclimated to freshwater (FW) and then exposed to seawater (SW) was examined. There was 18% mortality in the first 3 days following exposure to SW, with no mortalities in the FW control group. Gill Na⁺/K⁺-ATPase (NKA) activity, an index of osmoregulatory capacity, increased 2.6-fold in SW fish peaking on days 7 and 14. Gill citrate synthase activity, an index of aerobic capacity, was 50-62% greater in SW than FW fish and peaked on day 7. Tissue water content was significantly lower in the SW fish on day 1 only, returning to FW levels by day 3. Routine metabolic rate was decreased within 24 h of SW exposure and was maintained slightly (8-22%) but significantly lower in SW compared to FW water controls throughout the 2-week experiment. These results indicate that elevated salinity resulted in increased SW osmoregulatory and aerobic capacity in the gill, but with a reduced whole animal metabolic rate to this euryhaline species.

Repeated Genetic Targets of Natural Selection Underlying Adaptation of Fishes to Changing Salinity

Ecological transitions across salinity boundaries have led to some of the most important diversification events in the animal kingdom, especially among fishes. Adaptations accompanying such transitions include changes in morphology, diet, whole-organism performance, and osmoregulatory function, which may be particularly prominent since divergent salinity regimes make opposing demands on systems that maintain ion and water balance. Research in the last decade has focused on the genetic targets underlying such adaptations, most notably by comparing populations of species that are distributed across salinity boundaries. Here, we synthesize research on the targets of natural selection using whole-genome approaches, with a particular emphasis on the osmoregulatory system. Given the complex, integrated and polygenic nature of this system, we expected that signatures of natural selection would span numerous genes across functional levels of osmoregulation, especially salinity sensing, hormonal control, and cellular ion exchange mechanisms. We find support for this prediction: genes coding for V-type, Ca2+, and Na+/K+-ATPases, which are key cellular ion exchange enzymes, are especially common targets of selection in species from six orders of fishes. This indicates that while polygenic selection contributes to adaptation across salinity boundaries, changes in ATPase enzymes may be of particular importance in supporting such transitions.

Epigenetic variation associated with responses to different habitats in the context of genetic divergence in Phragmites australis

The mechanisms underlying heritable phenotypic divergence associated with adaptation in response to environmental stresses may involve both genetic and epigenetic variations. Several prior studies have revealed even higher levels of epigenetic variation than genetic variation. However, few population-level studies have explored the effects of epigenetic variation on species with high levels of genetic diversity distributed across different habitats. Using AFLP and methylation-sensitive AFLP markers, we tested the hypothesis that epigenetic variation may contribute to differences in plants occupying different habitats when genetic variation alone cannot fully explain adaptation. As a cosmopolitan invasive species, Phragmites australis (common reed) together with high genetic diversity and remarkable adaptability has been suggested as a model for responses to global change and indicators of environmental fluctuations. We found high levels of genetic and epigenetic diversity and significant genetic/epigenetic structure within each of 12 studied populations sampled from four natural habitats of P. australis. Possible adaptive epigenetic variation was suggested by significant correlations between DNA methylation-based epigenetic differentiation and adaptive genetic divergence in populations across the habitats. Meanwhile, various AMOVAs indicated that some epigenetic differences may respond to various local habitats. A partial Mantel test was used to tease out the correlations between genetic/epigenetic variation and habitat after controlling for the correlation between genetic and epigenetic variations. We found that epigenetic diversity was affected mostly by soil nutrient availability, suggesting that at least some epigenetic differentiation occurred independently of genetic variation. We also found stronger correlations between epigenetic variation and phenotypic traits than between genetic variation and such traits. Overall, our findings indicate that genetically based differentiation correlates with heterogeneous habitats, while epigenetic variation plays an important role in ecological differentiation in natural populations of P. australis. In addition, our results suggest that when assessing global change responses of plant species, intraspecific variation needs to be considered.

Sperm adaptation in relation to salinity in three goby species

In externally fertilizing species, the gametes of both males and females are exposed to the influences of the environment into which they are released. Sperm are sensitive to abiotic factors such as salinity, but they are also affected by biotic factors such as sperm competition. In this study, the authors compared the performance of sperm of three goby species, the painted goby, Pomatoschistus pictus, the two-spotted goby, Pomatoschistus flavescens, and the sand goby, Pomatoschistus minutus. These species differ in their distributions, with painted goby having the narrowest salinity range and sand goby the widest. Moreover, data from paternity show that the two-spotted goby experiences the least sperm competition, whereas in the sand goby sperm competition is ubiquitous. The authors took sperm samples from dissected males and exposed them to high salinity water (31 PSU) representing the North Sea and low salinity water (6 PSU) representing the brackish Baltic Sea Proper. They then used computer-assisted sperm analysis to measure the proportion of motile sperm and sperm swimming speed 10 min and 20 h after sperm activation. The authors found that sperm performance depended on salinity, but there seemed to be no relationship to the species' geographical distribution in relation to salinity range. The species differed in the proportion of motile sperm, but there was no significant decrease in sperm motility during 20 h. The sand goby was the only species with motile sperm after 72 h.

Phenotypic stability in scalar calcium of freshwater fish across a wide range of aqueous calcium availability in nature

Spatial environmental gradients can promote adaptive differences among conspecific populations as a result of local adaptation or phenotypic plasticity. Such divergence can be opposed by various constraints, including gene flow, limited genetic variation, temporal fluctuations, or developmental constraints. We focus on the constraint that can be imposed when some populations are found in locations characterized by low levels of an essential nutrient. We use scales of wild fish to investigate phenotypic effects of spatial variation in a potentially limiting nutrient-calcium. If scale calcium (we use "scalar" calcium for consistency with the physiology literature) simply reflects environmental calcium availability, we expect higher levels of scalar calcium in fish from calcium-rich water, compared to fish from calcium-poor water. To consider this "passive response" scenario, we analyzed scalar calcium concentrations from three native fish species (Lepomis gibbosus, Percina caprodes, and Perca flavescens) collected at multiple sites across a dissolved calcium gradient in the Upper St. Lawrence River. Contradicting the "passive response" scenario, we did not detect strong or consistent relationships between scalar calcium and water calcium. Instead, for a given proportional increase in water calcium across the wide environmental gradient, the corresponding proportional change in scalar calcium was much smaller. We thus favor the alternative "active homeostasis" scenario, wherein fish from calcium-poor water are better able to uptake, mobilize, and deposit calcium than are fish from calcium-rich water. We further highlight the importance of studying functional traits, such as scales, in their natural setting as opposed to only laboratory studies.

Evolutionary mismatch along salinity gradients in a Neotropical water strider

The evolution of local adaptation is crucial for the in situ persistence of populations in changing environments. However, selection along broad environmental gradients could render local adaptation difficult, and might even result in maladaptation. We address this issue by quantifying fitness trade-offs (via common garden experiments) along a salinity gradient in two populations of the Neotropical water strider Telmatometra withei-a species found in both fresh (FW) and brackish (BW) water environments across Panama. We found evidence for local adaptation in the FW population in its home FW environment. However, the BW population showed only partial adaptation to the BW environment, with a high magnitude of maladaptation along naturally occurring salinity gradients. Indeed, its overall fitness was ~60% lower than that of the ancestral FW population in its home environment, highlighting the role of phenotypic plasticity, rather than local adaptation, in high salinity environments. This suggests that populations seemingly persisting in high salinity environments might in fact be maladapted, following drastic changes in salinity. Thus, variable selection imposed by salinization could result in evolutionary mismatch, where the fitness of a population is displaced from its optimal environment. Understanding the fitness consequences of persisting in fluctuating salinity environments is crucial to predict the persistence of populations facing increasing salinization. It will also help develop evolutionarily informed management strategies in the context of global change.

Pipefish Locally Adapted to Low Salinity in the Baltic Sea Retain Phenotypic Plasticity to Cope With Ancestral Salinity Levels

Genetic adaptation and phenotypic plasticity facilitate the migration into new habitats and enable organisms to cope with a rapidly changing environment. In contrast to genetic adaptation that spans multiple generations as an evolutionary process, phenotypic plasticity allows acclimation within the life-time of an organism. Genetic adaptation and phenotypic plasticity are usually studied in isolation, however, only by including their interactive impact, we can understand acclimation and adaptation in nature. We aimed to explore the contribution of adaptation and plasticity in coping with an abiotic (salinity) and a biotic (Vibriobacteria) stressor using six different populations of the broad-nosed pipefishSyngnathus typhlethat originated from either high [14–17 Practical Salinity Unit (PSU)] or low (7–11 PSU) saline environments along the German coastline of the Baltic Sea. We exposed wild caught animals, to either high (15 PSU) or low (7 PSU) salinity, representing native and novel salinity conditions and allowed animals to mate. After male pregnancy, offspring was split and each half was exposed to one of the two salinities and infected withVibrio alginolyticusbacteria that were evolved at either of the two salinities in a fully reciprocal design. We investigated life-history traits of fathers and expression of 47 target genes in mothers and offspring. Pregnant males originating from high salinity exposed to low salinity were highly susceptible to opportunistic fungi infections resulting in decreased offspring size and number. In contrast, no signs of fungal infection were identified in fathers originating from low saline conditions suggesting that genetic adaptation has the potential to overcome the challenges encountered at low salinity. Offspring from parents with low saline origin survived better at low salinity suggesting genetic adaptation to low salinity. In addition, gene expression analyses of juveniles indicated patterns of local adaptation,trans-generational plasticity and developmental plasticity. In conclusion, our study suggests that pipefish are locally adapted to the low salinity in their environment, however, they are retaining phenotypic plasticity, which allows them to also cope with ancestral salinity levels and prevailing pathogens.

Genetic response to human-induced habitat changes in the marine environment: A century of evolution of European sprat in Landvikvannet, Norway

Habitat changes represent one of the five most pervasive threats to biodiversity. However, anthropogenic activities also have the capacity to create novel niche spaces to which species respond differently. In 1880, one such habitat alterations occurred in Landvikvannet, a freshwater lake on the Norwegian coast of Skagerrak, which became brackish after being artificially connected to the sea. This lake is now home to the European sprat, a pelagic marine fish that managed to develop a self-recruiting population in barely few decades. Landvikvannet sprat proved to be genetically isolated from the three main populations described for this species; that is, Norwegian fjords, Baltic Sea, and the combination of North Sea, Kattegat, and Skagerrak. This distinctness was depicted by an accuracy self-assignment of 89% and a highly significant F ST between the lake sprat and each of the remaining samples (average of ≈0.105). The correlation between genetic and environmental variation indicated that salinity could be an important environmental driver of selection (3.3% of the 91 SNPs showed strong associations). Likewise, Isolation by Environment was detected for salinity, although not for temperature, in samples not adhering to an Isolation by Distance pattern. Neighbor-joining tree analysis suggested that the source of the lake sprat is in the Norwegian fjords, rather than in the Baltic Sea despite a similar salinity profile. Strongly drifted allele frequencies and lower genetic diversity in Landvikvannet compared with the Norwegian fjords concur with a founder effect potentially associated with local adaptation to low salinity. Genetic differentiation (F ST) between marine and brackish sprat is larger in the comparison Norway-Landvikvannet than in Norway-Baltic, which suggests that the observed divergence was achieved in Landvikvannet in some 65 generations, that is, 132 years, rather than gradually over thousands of years (the age of the Baltic Sea), thus highlighting the pace at which human-driven evolution can happen.

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.

Invertebrate turnover along gradients of anthropogenic salinisation in rivers of two German regions

Rising salinity in freshwater ecosystems can affect community composition. Previous studies mainly focused on changes in freshwater communities along gradients of absolute levels of electrical conductivity (EC). However, both geogenic and anthropogenic drivers contribute to the EC level and taxa may regionally be adapted to geogenic EC levels. Therefore, we examined the turnover in freshwater invertebrates along gradients of anthropogenic EC change in two regions of Germany. The anthropogenic change of EC was estimated as the difference between the measured EC and the modeled background EC driven by geochemical and climate variables. Turnover in freshwater invertebrates (β-diversity) was estimated using the Jaccard index (JI). We found that invertebrate turnover between EC gradient categories is generally greater than 47%, with a maximum of approximately 70% in sites with a more than 0.4 mS cm^-1 change compared to the baseline (i.e. no difference between predicted and measured EC). The invertebrates Amphinemura sp., Anomalopterygella chauviniana and Leuctra sp. were reliable indicators of low EC change, whereas Potamopyrgus antipodarum indicated sites with the highest EC change. Variability within categories of EC change was slightly lower than within categories of absolute EC. Elevated nutrient concentrations that are often linked to land use may have contributed to the observed change of the invertebrate richness and can exacerbate effects of EC on communities in water. Overall, our study suggests that the change in EC, quantified as the difference between measured EC and modeled background EC, can be used to examine the response of invertebrate communities to increasing anthropogenic salinity concentrations in rivers. However, due to the strong correlation between EC change and observed EC in our study regions, the response to these two variables was very similar. Further studies in areas where EC change and observed EC are less correlated are required. In addition, such studies should consider the change in specific ions.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation

Epigenetic inheritance has been proposed to contribute to adaptation and acclimation via two information channels: (i) inducible epigenetic marks that enable transgenerational plasticity and (ii) noninducible epigenetic marks resulting from random epimutations shaped by selection. We studied both postulated channels by sequencing methylomes and genomes of Baltic three-spined sticklebacks (Gasterosteus aculeatus) along a salinity cline. Wild populations differing in salinity tolerance revealed differential methylation (pop-DMS) at genes enriched for osmoregulatory processes. A two-generation experiment demonstrated that 62% of these pop-DMS were noninducible by salinity manipulation, suggesting that they are the result of either direct selection or associated genomic divergence at cis- or trans-regulatory sites. Two-thirds of the remaining inducible pop-DMS increased in similarity to patterns detected in wild populations from corresponding salinities. The level of similarity accentuated over consecutive generations, indicating a mechanism of transgenerational plasticity. While we can attribute natural DNA methylation patterns to the two information channels, their interplay with genomic variation in salinity adaptation is still unresolved.

Population Genomics in Rhamdia quelen (Heptapteridae, Siluriformes) Reveals Deep Divergence and Adaptation in the Neotropical Region

Rhamdia quelen, a Neotropical fish with hybridization between highly divergent mitochondrial DNA (mtDNA) lineages, represents an interesting evolutionary model. Previous studies suggested that there might be demographic differences between coastal lagoons and riverine environments, as well as divergent populations that could be reproductively isolated. Here, we investigated the genetic diversity pattern of this taxon in the Southern Neotropical Basin system that includes the La Plata Basin, Patos-Merin lagoon basin and the coastal lagoons draining to the SW Atlantic Ocean, through a population genomics approach using 2b-RAD-sequencing-derived single nucleotide polymorphisms (SNPs). The genomic scan identified selection footprints associated with divergence and suggested local adaptation environmental drivers. Two major genomic clusters latitudinally distributed in the Northern and Southern basins were identified, along with consistent signatures of divergent selection between them. Population structure based on the whole set of loci and on the presumptive neutral vs. adaptive loci showed deep genomic divergence between the two major clusters. Annotation of the most consistent SNPs under divergent selection revealed some interesting candidate genes for further functional studies. Moreover, signals of adaptation to a coastal lagoon environment mediated by purifying selection were found. These new insights provide a better understanding of the complex evolutionary history of R. quelen in the southernmost basin of the Neotropical region.

Intraspecific variation in the response of the estuarine European isopod Cyathura carinata (Krøyer, 1847) to ocean acidification

In the present study the model isopod, Cyathura carinata were exposed to four pH_NIST treatments (control: 7.9; 7.5, 7, 6.5) in order to determine the tolerance and pH threshold value this estuarine species withstand under future acidification scenarios. Seawater acidification significantly affected the lifespan of C. carinata, where population density was remarkably reduced at the lowest pH treatment. The longevity, survivorship and swimming activity (related to the acquisition of energy) of these isopods decreased with decreasing pH. Furthermore, to determine the possible metabolic plasticity of this species, the swimming activity, the Na+/K + -ATPase activity (relevant for osmoregulation process), and the RNA:DNA ratio (an indicator of fitness) were measure from two populations of C. carinata, one inhabiting a stable environment (pH_NIST 7.5-8.0) and one inhabiting a fluctuating pCO₂ regimes (pH 3.3-8.5) subjected to three pH treatments (7.9, 7.0 and 6.5). The population from high fluctuating pCO₂ conditions showed capacity to withstand to pH 6.5, as well as higher longevity and metabolic plasticity, when compared with the population from the habitat with slight pCO₂ variation. These results indicate that Cyathura population from stable environments could be vulnerable to ocean acidification because it could trigger detrimental effects on its survival energy budget, and growth. However, ocean acidification has limited effect on the energy budget and survival of C. carinata population from highly variable habitats, suggesting that they are able to cope with the elevated energy demand. The difference showed between populations is likely an indication of genetic differentiation in tolerance to ocean acidification, possibly attributable to local adaptations, which could provide the raw material necessary for adaptation to future conditions. In addition, our results suggest that when assessing marine crustacean responses to changing environments on a global scale, variability in population and metabolic responses need to be considered.

Ecological factors and morphological traits are associated with repeated genomic differentiation between lake and stream stickleback

The repeated evolution of similar phenotypes in independent populations (i.e. parallel or convergent evolution) provides an opportunity to identify genetic and ecological factors that influence the process of adaptation. Threespine stickleback fish ( Gasterosteus aculeatus) are an excellent model for such studies, as they have repeatedly adapted to divergent habitats across the Northern hemisphere. Here, we use genomic, ecological and morphological data from 16 independent pairs of stickleback populations adapted to divergent lake and stream habitats. We combine a population genomic approach to identify regions of the genome that are likely under selection in these divergent habitats with an association mapping approach to identify regions of the genome that underlie variation in ecological factors and morphological traits. Over 37% of genomic windows are repeatedly differentiated across lake-stream pairs. Similarly, many genomic windows are associated with variation in abiotic factors, diet items and morphological phenotypes. Both the highly differentiated windows and candidate trait windows are non-randomly distributed across the genome and show some overlap. However, the overlap is not significant on a genome-wide scale. Together, our data suggest that adaptation to divergent food resources and predation regimes are drivers of differentiation in lake-stream stickleback, but that additional ecological factors are also important. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species

Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).

Divergent transcriptomic signatures in response to salinity exposure in two populations of an estuarine fish

In estuary and coastal systems, human demand for freshwater, climate change-driven precipitation variability, and extreme weather impact salinity levels, reducing connectivity between mesohaline coastal fish populations and potentially contributing to genomic divergence. We examined gill transcriptome responses to salinity in wild-caught juveniles from two populations of Sacramento splittail (Pogonichthys macrolepidotus), a species of conservation concern that is endemic to the San Francisco Estuary, USA, and the lower reaches of its tributaries. Recent extreme droughts have led to salinities above the tolerance limits for this species, creating a migration barrier between these populations, which potentially contributed to population divergence. We identified transcripts involved in a conserved response to salinity; however, the more salinity-tolerant San Pablo population had greater transcriptome plasticity (3.6-fold more transcripts responded than the Central Valley population) and a response consistent with gill remodeling after 168 hr of exposure to elevated salinity. The reorganization of the gill in response to changing osmotic gradients is a process critical for acclimation and would facilitate enhanced salinity tolerance. We detected an upregulation of receptors that control the Wnt (wingless-type) cell signaling pathway that may be required for an adaptive response to increases in salinity, patterns not observed in the relatively salinity-sensitive Central Valley population. We detected 62 single nucleotide polymorphisms (SNPs) in coding regions of 26 transcripts that differed between the populations. Eight transcripts that contained SNPs were associated with immune responses, highlighting the importance of diversity in immune gene sequences as a defining characteristic of genomic divergence between these populations. Our data demonstrate that these populations have divergent transcriptomic responses to salinity, which is consistent with observed physiological differences in salinity tolerance.

Transgenerational plasticity and selection shape the adaptive potential of sticklebacks to salinity change

In marine climate change research, salinity shifts have been widely overlooked. While widespread desalination effects are expected in higher latitudes, salinity is predicted to increase closer to the equator. We took advantage of the steep salinity gradient of the Baltic Sea as a space-for-time design to address effects of salinity change on populations. Additionally, genetic diversity, a prerequisite for adaptive responses, is reduced in Baltic compared to Atlantic populations. On the one hand, adaptive transgenerational plasticity (TGP) might buffer the effects of environmental change, which may be of particular importance under reduced genetic variation. On the other hand, physiological trade-offs due to environmental stress may hamper parental provisioning to offspring thereby intensifying the impact of climate change across generations (nonadaptive TGP). Here, we studied both hypothesis of adaptive and nonadaptive TGP in the three-spined stickleback (Gasterosteus aculeatus) fish model along the strong salinity gradient of the Baltic Sea in a space-for-time experiment. Each population tolerated desalination well, which was not altered by parental exposure to low salinity. Despite a common marine ancestor, populations locally adapted to low salinity lost their ability to cope with fully marine conditions, resulting in lower survival and reduced relative fitness. Negative transgenerational effects were evident in early life stages, but disappeared after selection via mortality occurred during the first 12-30 days posthatch. Modeling various strengths of selection, we showed that nonadaptive transgenerational plasticity accelerated evolution by increasing directional selection within the offspring generation. Qualitatively, when genetic diversity is large, we predict that such effects will facilitate rapid adaptation and population persistence, while below a certain threshold populations suffer a higher risk of local extinction. Overall, our results suggest that transgenerational plasticity and selection are not independent of each other and thereby highlight a current gap in TGP studies.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Ecological variation along the salinity gradient in the Baltic Sea Area and its consequences for reproduction in the common goby

Although it has become clear that sexual selection may shape mating systems and drive speciation, the potential constraints of environmental factors on processes and outcomes of sexual selection are largely unexplored. Here, we investigate the geographic variation of such environmental factors, more precisely the quality and quantity of nest resources (bivalve shells) along a salinity gradient in the Baltic Sea Area (Baltic Sea, Sounds and Belts, and Kattegat). We further test whether we find any salinity-associated morphological differences in body size between populations of common gobies Pomatoschistus microps, a small marine fish with a resource-based mating system. In a geographically expansive field study, we sampled 5 populations of P. microps occurring along the salinity gradient (decreasing from West to East) in the Baltic Sea Area over 3 consecutive years. Nest resource quantity and quality decreased from West to East, and a correlation between mussel size and male body size was detected. Population density, sex ratios, mating- and reproductive success as well as brood characteristics also differed between populations but with a less clear relation to salinity. With this field study we shed light on geographic variation of distinct environmental parameters possibly acting on population differentiation. We provide insights on relevant ecological variation, and draw attention to its importance in the framework of context-dependent plasticity of sexual selection.

Variation in salinity tolerance between and within anadromous subpopulations of pike (Esox 1ucius)

Environmental heterogeneity is a key determinant of genetic and phenotypic diversity. Stable and homogenous environments tends to result in evolution of specialism and local adaptations, while temporally unpredictable environments may maintain a diversity of specialists, promote generalist strategies, or favour diversified bet hedging strategies. We compared salinity tolerance between two anadromous subpopulations of pike (Esox Lucius) that utilize freshwater spawning sites with different salinity regimes. Eggs from each population were artificially fertilized and incubated in a salinity gradient (0, 3, 5, 7, and 9 psu) using a split-brood design. Effects on embryonic development, hatching success, survival of larvae, and fry body length were compared between populations and families. The population naturally spawning in the stable freshwater habitat showed signs of specialization for freshwater spawning. The population exposed to fluctuating selective pressure in a spawning area with occasional brackish water intrusions tolerated higher salinities and displayed considerable variation in reaction norms. Genetic differences and plasticity of salinity tolerance may enable populations to cope with changes in salinity regimes associated with future climate change. That geographically adjacent subpopulations can constitute separate units with different genetic characteristics must be considered in management and conservation efforts to avoid potentially negative effects of genetic admixture on population fitness and persistence.

Kidney morphology and candidate gene expression shows plasticity in sticklebacks adapted to divergent osmotic environments

Novel physiological challenges in different environments can promote the evolution of divergent phenotypes, either through plastic or genetic changes. Environmental salinity serves as a key barrier to the distribution of nearly all aquatic organisms, and species diversification is likely to be enabled by adaptation to alternative osmotic environments. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations found both in marine and freshwater environments. It has evolved both highly plastic and locally adapted phenotypes due to salinity-derived selection, but the physiological and genetic basis of adaptation to salinity is not fully understood. We integrated comparative cellular morphology of the kidney, a key organ for osmoregulation, and candidate gene expression to explore the underpinnings of evolved variation in osmotic plasticity within two populations of sticklebacks from distinct salinity zones in the Baltic Sea: the high salinity Kattegat, representative of the ancestral marine habitat, and the low salinity Bay of Bothnia. A common-garden experiment revealed that kidney morphology in the ancestral high salinity population had a highly plastic response to salinity conditions, whereas this plastic response was reduced in the low salinity population. Candidate gene expression in kidney tissue revealed a similar pattern of population-specific differences, with a higher degree of plasticity in the native high salinity population. Together these results suggest that renal cellular morphology has become canalized to low salinity, and that these structural differences may have functional implications for osmoregulation.

Salinized rivers: degraded systems or new habitats for salt-tolerant faunas?

Anthropogenic salinization of rivers is an emerging issue of global concern, with significant adverse effects on biodiversity and ecosystem functioning. Impacts of freshwater salinization on biota are strongly mediated by evolutionary history, as this is a major factor determining species physiological salinity tolerance. Freshwater insects dominate most flowing waters, and the common lotic insect orders Ephemeroptera (mayflies), Plecoptera (stoneflies) and Trichoptera (caddisflies) are particularly salt-sensitive. Tolerances of existing taxa, rapid adaption, colonization by novel taxa (from naturally saline environments) and interactions between species will be key drivers of assemblages in saline lotic systems. Here we outline a conceptual framework predicting how communities may change in salinizing rivers. We envision that a relatively small number of taxa will be saline-tolerant and able to colonize salinized rivers (e.g. most naturally saline habitats are lentic; thus potential colonizers would need to adapt to lotic environments), leading to depauperate communities in these environments.

Competitive advantage and higher fitness in native populations of genetically structured planktonic diatoms

It has been shown that the planktonic diatom Skeletonema from neighbouring areas are genetically differentiated despite absence of physical dispersal barriers. We revisited two sites, Mariager Fjord and Kattegat, NE Atlantic, and isolated new strains. Microsatellite genotyping and F-statistics revealed that the populations were genetically differentiated. An experiment was designed to investigate if populations are locally adapted and have a native competitive advantage. Ten strains from each location were grown individually in native and foreign water to investigate differences in produced biomass. Additionally, we mixed six pairs, one strain from each site, and let them grow together in native and foreign water. Strains from Mariager Fjord and Kattegat produced higher biomass in native water. In the competition experiment, strains from both sites displayed higher relative abundance and demonstrated competitive advantage in their native water. The cause of the differentiated growth is unknown, but could possibly be attributed to differences in silica concentration or viruses in the two water types. Our data show that dispersal potential does not influence the genetic structure of the populations. We conclude that genetic adaptation has not been overruled by gene flow, but instead the responses to different selection conditions are enforcing the observed genetic structure.

Morphological differences between habitats are associated with physiological and behavioural trade-offs in stickleback (Gasterosteus aculeatus)

Local specialization can be advantageous for individuals and may increase the resilience of the species to environmental change. However, there may be trade-offs between morphological responses and physiological performance and behaviour. Our aim was to test whether habitat-specific morphology of stickleback (Gasterosteus aculeatus) interacts with physiological performance and behaviour at different salinities. We rejected the hypothesis that deeper body shape of fish from habitats with high predation pressure led to decreases in locomotor performance. However, there was a trade-off between deeper body shape and muscle quality. Muscle of deeper-bodied fish produced less force than that of shallow-bodied saltmarsh fish. Nonetheless, saltmarsh fish had lower swimming performance, presumably because of lower muscle mass overall coupled with smaller caudal peduncles and larger heads. Saltmarsh fish performed better in saline water (20 ppt) relative to freshwater and relative to fish from freshwater habitats. However, exposure to salinity affected shoaling behaviour of fish from all habitats and shoals moved faster and closer together compared with freshwater. We show that habitat modification can alter phenotypes of native species, but local morphological specialization is associated with trade-offs that may reduce its benefits.

Genetic basis for body size variation between an anadromous and two derived lacustrine populations of threespine stickleback Gasterosteus aculeatus in southwest Alaska

Body size is a highly variable trait among geographically separated populations. Size-assortative reproductive isolation has been linked to recent adaptive radiations of threespine stickleback (Gasterosteus aculeatus) into freshwater, but the genetic basis of the commonly found size difference between anadromous and derived lacustrine sticklebacks has not been tested. We studied the genetic basis of size differences between recently diverging stickleback lineages in southwest Alaska using a common environment experiment. We crossed stickleback within one anadromous (Naknek River) and one lake (Pringle Lake) population and between the anadromous and two lake populations (Pringle and JoJo Lakes), and raised them in a salinity of 4–6 ppt. The F1 anadromous and freshwater forms differed significantly in size, whereas hybrids were intermediate or exhibited dominance toward the anadromous form. Additionally, the size of freshwater F1s differed from their wild counterparts, with within-population F1s from Pringle Lake growing larger than their wild counterparts, while there was no size difference between lab-raised and wild anadromous fish. Sexual dimorphism was always present in anadromous fish, but not in freshwater, and not always in the hybrid crosses. These results, along with parallel changes among anadromous and freshwater forms in other regions, suggest that this heritable trait is both plastic and may be under divergent and/or sexual selection.

Small-Scale Habitat-Specific Variation and Adaptive Divergence of Photosynthetic Pigments in Different Alkali Soils in Reed Identified by Common Garden and Genetic Tests

Flexibility of photosynthetic pigment traits is an important adaptive mechanism through which plants can increase mean fitness in a variable environment. Unlike morphological traits in plants, photosythesis has been shown to exhibit phenotypic plasticity, responding rapidly to environmental conditions. Meanwhile, local adaptation at small scales is considered to be rare. Thus, detecting the small-scale adaptive divergence of photosynthetic pigments presents a challenge. Leaf concentrations of photosynthetic pigments under stressful conditions may be reduced or maintained. Concentrations of some pigments and/or ratio of Chlorophyll a (Chla) to Chlorophyll b (Chlb) do not change markedly in some species, such as the common reed, Phragmites australis, a cosmopolitan grass and common invader. Little is known about photosynthetic responses of this plant to varying levels of alkali salt. Few studies have attempted to account for the relationship between pigment accumulation and leaf position in wild plant populations in grasslands. In this study, photosynthetic pigment concentrations and the total Chl(a+b)/Car ratio decreased as the growing season progressed and were shown to be significantly lower in the habitat with a higher soil pH value and less moisture when compared between habitats. The Chla/Chlb ratio did not differ significantly between habitats, although it increased significantly over time. Leaves in the middle position may be functionally important in the response to soil conditions because only pigment concentrations and the Chl(a+b)/Car ratio of those leaves varied between habitats significantly. The outlier loci, used to evaluate molecular signatures of selection, were detected by Arlequin, Bayescan, and Bayenv analyses. In the simulated habitats of common garden, the local genotypes had higher values of Chla, Chlb, Chl(a+b), Car in their home habitat than did genotypes originating from the other habitat. Q_ST-F_ST comparisons provided evidence of divergent selection. It appears likely that soil moisture, pH and electric conductivity drove local adaptation. Combined approaches that utilize information on phenotypes from field and common garden experiments, genome-wide markers, and environmental data will be the most informative for understanding the adaptive nature of the intraspecific divergence.

Evolution of stickleback in 50 years on earthquake-uplifted islands

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Common garden experiments in the genomic era: new perspectives and opportunities

The study of local adaptation is rendered difficult by many evolutionary confounding phenomena (for example, genetic drift and demographic history). When complex traits are involved in local adaptation, phenomena such as phenotypic plasticity further hamper evolutionary biologists to study the complex relationships between phenotype, genotype and environment. In this perspective paper, we suggest that the common garden experiment, specifically designed to deal with phenotypic plasticity, has a clear role to play in the study of local adaptation, even (if not specifically) in the genomic era. After a quick review of some high-throughput genotyping protocols relevant in the context of a common garden, we explore how to improve common garden analyses with dense marker panel data and recent statistical methods. We then show how combining approaches from population genomics and genome-wide association studies with the settings of a common garden can yield to a very efficient, thorough and integrative study of local adaptation. Especially, evidence from genomic (for example, genome scan) and phenotypic origins constitute independent insights into the possibility of local adaptation scenarios, and genome-wide association studies in the context of a common garden experiment allow to decipher the genetic bases of adaptive traits.

Population genomic evidence for adaptive differentiation in Baltic Sea three-spined sticklebacks

Background

The degree of genetic differentiation among populations experiencing high levels of gene flow is expected to be low for neutral genomic sites, but substantial divergence can occur in sites subject to directional selection. Studies of highly mobile marine fish populations provide an opportunity to investigate this kind of heterogeneous genomic differentiation, but most studies to this effect have focused on a relatively low number of genetic markers and/or few populations. Hence, the patterns and extent of genomic divergence in high-gene-flow marine fish populations remain poorly understood.

Results

We here investigated genome-wide patterns of genetic variability and differentiation in ten marine populations of three-spined stickleback (Gasterosteus aculeatus) distributed across a steep salinity and temperature gradient in the Baltic Sea, by utilizing >30,000 single nucleotide polymorphisms obtained with a pooled RAD-seq approach. We found that genetic diversity and differentiation varied widely across the genome, and identified numerous fairly narrow genomic regions exhibiting signatures of both divergent and balancing selection. Evidence was uncovered for substantial genetic differentiation associated with both salinity and temperature gradients, and many candidate genes associated with local adaptation in the Baltic Sea were identified.

Conclusions

The patterns of genetic diversity and differentiation, as well as candidate genes associated with adaptation, in Baltic Sea sticklebacks were similar to those observed in earlier comparisons between marine and freshwater populations, suggesting that similar processes may be driving adaptation to brackish and freshwater environments. Taken together, our results provide strong evidence for heterogenic genomic divergence driven by local adaptation in the face of gene flow along an environmental gradient in the post-glacially formed Baltic Sea.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0130-8) contains supplementary material, which is available to authorized users.

Reciprocal osmotic challenges reveal mechanisms of divergence in phenotypic plasticity in the killifish Fundulus heteroclitus

The killifish Fundulus heteroclitus is an estuarine species with broad physiological plasticity enabling acclimation to diverse stressors. Previous work suggests freshwater populations expanded their physiology to accommodate low salinity environments, however, it is unknown if this compromises their tolerance to high salinity. We employed a comparative approach to investigate the mechanisms of a derived freshwater phenotype and the fate of an ancestral euryhaline phenotype after invasion of a freshwater environment. We compared physiological and transcriptomic responses to high and low salinity stress in fresh and brackish water populations and found an enhanced plasticity to low salinity in the freshwater population coupled with a reduced ability to acclimate to high salinity. Transcriptomic data identified genes with a conserved common response, a conserved salinity dependent response, and responses associated with population divergence. Conserved common acclimation responses revealed stress responses and alterations in cell-cycle regulation as important mechanisms in the general osmotic response. Salinity-specific responses included the regulation of genes involved in ion transport, intracellular calcium, energetic processes, and cellular remodeling. Genes diverged between populations were primarily those showing salinity-specific expression and included those regulating polyamine homeostasis and cell cycle. Additionally, when populations were matched with their native salinity, expression patterns were consistent with the concept of “transcriptomic resilience,” suggesting local adaptation. These findings provide insight into the fate of a plastic phenotype after a shift in environmental salinity and help to reveal mechanisms allowing for euryhalinity.

Relaxed selection causes microevolution of seawater osmoregulation and gene expression in landlocked Alewives

Ecological transitions from marine to freshwater environments have been important in the creation of diversity among fishes. Evolutionary changes associated with these transitions likely involve modifications of osmoregulatory function. In particular, relaxed selection on hypo-osmoregulation should strongly affect animals that transition into novel freshwater environments. We used populations of the Alewife (Alosa pseudoharengus) to study evolutionary shifts in hypo-osmoregulatory capacity and ion regulation associated with freshwater transitions. Alewives are ancestrally anadromous, but multiple populations in Connecticut have been independently restricted to freshwater lakes; these landlocked populations complete their entire life cycle in freshwater. Juvenile landlocked and anadromous Alewives were exposed to three salinities (1, 20 and 30 ppt) in small enclosures within the lake. We detected strong differentiation between life history forms: landlocked Alewives exhibited reduced seawater tolerance and hypo-osmoregulatory performance compared to anadromous Alewives. Furthermore, gill Na(+)/K(+)-ATPase activity and transcription of genes for seawater osmoregulation (NKCC-Na(+)/K(+)/2Cl(-) cotransporter and CFTR-cystic fibrosis transmembrane conductance regulator) exhibited reduced responsiveness to seawater challenge. Our study demonstrates that adaptations of marine-derived species to completely freshwater life cycles involve partial loss of seawater osmoregulatory performance mediated through changes to ion regulation in the gill.