Stickleback Populations Research Articles

An experimental test of how parasites of predators can influence trophic cascades and ecosystem functioning.

Parasites can shape the structure and function of ecosystems by influencing both the density and traits of their hosts. Such changes in ecosystems are particularly likely when the host is a predator that mediates the dynamics of trophic cascades. Here, we experimentally tested how parasite load of a small predatory fish, the threespine stickleback, can affect the occurrence and strength of trophic cascades and ecosystem functioning. In a factorial mesocosm experiment, we manipulated the density of stickleback (low vs. high), and the level of parasite load (natural vs. reduced). In addition, we used two stickleback populations from different lineages: an eastern European lineage with a more pelagic phenotype (Lake Constance) and a western European lineage with a more benthic phenotype (Lake Geneva). We found that stickleback caused trophic cascades in the pelagic but not the benthic food chain. Evidence for pelagic trophic cascades was stronger in treatments where parasite load of stickleback was reduced with an antihelmintic medication, and where fish originated from Lake Constance (i.e., the more pelagic lineage). A structural equationmodel revealed that differences in stickleback lineage and parasite load were most likely to impact trophic cascades via changes in the composition, rather than overall biomass, of zooplankton communities. Overall, our results provide experimental evidence that parasites of predators can influence the cascading effects of fish on lower trophic levels with consequences on ecosystem functioning.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish.

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination “hotspots” with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

Genome Characterization, Prevalence, and Transmission Mode of a Novel Picornavirus Associated with the Threespine Stickleback Fish (Gasterosteus aculeatus).

The complete genome sequence of an RNA virus was assembled from RNA sequencing of virus particles purified from threespine stickleback intestine tissue samples. This new virus is most closely related to the Eel picornavirus and can be assigned to the genus Potamipivirus in the family Picornaviridae Its unique genetic properties are enough to establish a new species, dubbed the Threespine Stickleback picornavirus (TSPV). Due to their broad geographic distribution throughout the Northern Hemisphere and parallel adaptation to freshwater, threespine sticklebacks have become a model in evolutionary ecology. Further analysis using diagnostic PCRs revealed that TSPV is highly prevalent in both anadromous and freshwater populations of threespine sticklebacks, infects almost all fish tissues, and is transmitted vertically to offspring obtained from in vitro fertilization in laboratory settings. Finally, TSPV was found in Sequence Reads Archives of transcriptome of Gasterosteus aculeatus, further demonstrating its wide distribution and unsought prevalence in samples. It is thus necessary to test the impact of TSPV on the biology of threespine sticklebacks, as this widespread virus could interfere with the behavioral, physiological, or immunological studies that employ this fish as a model system.IMPORTANCE The threespine stickleback species complex is an important model system in ecological and evolutionary studies because of the large number of isolated divergent populations that are experimentally tractable. For similar reasons, its coevolution with the cestode parasite Schistocephalus solidus, its interaction with gut microbes, and the evolution of its immune system are of growing interest. Herein we describe the discovery of an RNA virus that infects both freshwater and anadromous populations of sticklebacks. We show that the virus is transmitted vertically in laboratory settings and found it in Sequence Reads Archives, suggesting that experiments using sticklebacks were conducted in the presence of the virus. This discovery can serve as a reminder that the presence of viruses in wild-caught animals is possible, even when animals appear healthy. Regarding threespine sticklebacks, the impact of Threespine Stickleback picornavirus (TSPV) on the fish biology should be investigated further to ensure that it does not interfere with experimental results.

The role of variation and plasticity in parental care during the adaptive radiation of three-spine sticklebacks.

Phenotypic plasticity might influence evolutionary processes such as adaptive radiations. Plasticity in parental care might be especially effective in facilitating adaptive radiations if it allows populations to persist in novel environments. Here, we test the hypothesis that behavioral plasticity by parents in response to predation risk facilitated the adaptive radiation of three-spine sticklebacks. We compared the behavior of fathers across multiple ancestral (marine) and derived (freshwater) stickleback populations that differ in time since establishment. We measured behavioral plasticity in fathers in response to a predator found only in freshwater environments, simulating conditions marine males experience when colonizing freshwater. The antipredator behavior of males from newly established freshwater populations was intermediate between marine populations and well-established freshwater populations. In contrast to our predictions, on average, there was greater behavioral plasticity in derived freshwater populations than in ancestral marine populations. However, we found greater individual variation in behavioral reaction norms in marine populations compared to well-established freshwater populations, with newly established freshwater populations intermediate. This suggests that standing variation in behavioral reaction norms within ancestral populations might provide different evolutionary trajectories, and illustrates how plasticity can contribute to adaptive radiations.

Specificity of resistance and geographic patterns of virulence in a vertebrate host-parasite system

BackgroundHost genotype - parasite genotype co-evolutionary dynamics are influenced by local biotic and abiotic environmental conditions. This results in spatially heterogeneous selection among host populations. How such heterogeneous selection influences host resistance, parasite infectivity and virulence remains largely unknown. We hypothesized that different co-evolutionary trajectories of a vertebrate host-parasite association result in specific virulence patterns when assessed on a large geographic scale. We used two reference host populations of three-spined sticklebacks and nine strains of their specific cestode parasite Schistocephalus solidus from across the Northern Hemisphere for controlled infection experiments. Host and parasite effects on infection phenotypes including host immune gene expression were determined.ResultsS. solidus strains grew generally larger in hosts coming from a population with high parasite diversity and low S. solidus prevalence (DE hosts). Hosts from a population with low parasite diversity and high S. solidus prevalence (NO hosts) were better able to control the parasite’s growth, regardless of the origin of the parasite. Host condition and immunological parameters converged upon infection and parasite growth showed the same geographic pattern in both host types.ConclusionOur results suggest that NO sticklebacks evolved resistance against a variety of S. solidus strains, whereas DE sticklebacks are less resistant against S. solidus. Our data provide evidence that differences in parasite prevalence can cause immunological heterogeneity and that parasite size, a proxy for virulence and resistance, is, on a geographic scale, determined by main effects of the host and the parasite and less by an interaction of both genotypes.

Modelling population dynamics in mesocosms using an individual-based model coupled to a bioenergetics model

Developing population models is of great interest as these models enable to extrapolate toxicity observed at the molecular or individual levels to the population level, and thus improve environmental risk assessment of chemical substances. For this purpose, accounting for natural variations of environmental conditions and prey dynamics along with the chemical stress is needed. In this study, an individual-based model (IBM) coupled to a Dynamic Energy Budget (DEB) model was developed in order to predict three-spined stickleback population dynamics in semi-controlled stream experiments (in mesocosms). Datasets obtained in mesocosms offer the opportunity to develop and evaluate the model predictions. The most sensitive parameters of the DEB-IBM were identified by sensitivity analyses and were calibrated based on data from two independent mesocosm experiments. The predictive capacities of our model were subsequently evaluated using three independent mesocosm datasets under different environmental scenarios. Finally, our model was applied to a theoretical case of toxic effects to show an example of application of the model in a regulatory context. While the most uncertain population processes (in particular, competition for food in mesocosms) in our three-spined stickleback DEB-IBM could be modelled more accurately, our model can already serve to assess the impacts of toxicants at the population level by improving the analyses of mesocosm experiments, in decreasing the uncertainty of the experimental results. Therefore, in a second step, it could be used to predict the consequences on viability of a population exposed to a contaminant under various environmental and exposure scenarios.

A Single Interacting Species Leads to Widespread Parallel Evolution of the Stickleback Genome

Biotic interactions are potent, widespread causes of natural selection and divergent phenotypic evolution and can lead to genetic differentiation with gene flow among wild populations ("isolation by ecology") [1-4]. Biotic selection has been predicted to act onmore genes than abiotic selection thereby drivinggreater adaptation [5]. However, difficulties in isolating the genome-wide effect of single biotic agents of selection have limited our ability to identify and quantify the number and type of genetic regions responding to biotic selection [6-9]. We identified geographically interspersed lakes in which threespine stickleback fish (Gasterosteus aculeatus) have repeatedly adapted to the presence or absence of a single member of the ecological community, prickly sculpin (Cottus asper), a fish that is both a competitor and a predator of the stickleback [10]. Whole-genome sequencing revealed that sculpin presence or absence accounted for the majority of genetic divergence among stickleback populations, more so than geography. The major axis of genomic variation within and between the two lake types was correlated with multiple traits, indicating parallel natural selection across a gradient of biotic environments. A large proportion of the genome-about 1.8%, encompassing more than 600 genes-differentiated stickleback from the two biotic environments. Divergence occurred in 141 discrete genomic clumps located mainly in regions of low recombination, suggesting that genes brought to lakes by the colonizing ancestral population often evolved together in linked blocks. Strong selection and a wealth of standing genetic variation explain how a single member of the biotic community can have such a rapid and profound evolutionary impact.

The evolution of sex determination associated with a chromosomal inversion

Sex determination is a fundamentally important and highly diversified biological process, yet the mechanisms behind the origin of this diversity are mostly unknown. Here we suggest that the evolution of sex determination systems can be driven by a chromosomal inversion. We show that an XY system evolved recently in particular nine-spined stickleback (Pungitius pungitius) populations, which arose from ancient hybridization between two divergent lineages. Our phylogenetic and genetic mapping analyses indicate that the XY system is formed in a large inversion that is associated with hybrid sterility between the divergent lineages. We suggest that a new male-determining gene evolved in the inversion in response to selection against impaired male fertility in a hybridized population. Given that inversions are often associated with hybrid incompatibility in animals and plants, they might frequently contribute to the diversification of sex determination systems.

Final Preferred Temperature of Juvenile Three-Spined Stickleback Gasterosteus aculeatus (Linnaeus, 1758)

The preferred temperatures of juvenile three-spined stickleback from bodies of freshwater of the White Sea coastal region are studied under experimental conditions. The final preferred temperature in the young-of-the-year is 15.2°C, regardless of the acclimation temperature (11, 16, and 20°C). The increased duration of experiments (10 days) made it possible to track the process of the temperature selection as a whole, rather than observing its particular phases. A comparison of the original results with the literature data reveal possible differences in the temperature characteristics, both in geographically remote populations of the three-spined stickleback and between its freshwater and marine forms.

Parasite-induced colour alteration of intermediate hosts increases ingestion by suitable final host species

Abstract Parasites with complex life cycles often alter the phenotypic appearance of their intermediate hosts in order to facilitate ingestion by the final host. However, such manipulation can be costly as it might increase ingestion by less suitable or dead-end hosts as well. Species-specific parasitic manipulation is a way to enhance the transmission to suitable final hosts. Here, we experimentally show that the altered body colouration of the intermediate host Gammarus pulex caused by its acanthocephalan parasite Pomphorhynchus laevis differently affects predation by different fish species (barbel, perch, ruffe, brown trout and two populations of three-spined stickleback) depending on their suitability to act as final host. Species that were responsive to colour manipulation in a predation experiment were more susceptible to infection with P. laevis than unresponsive species. Furthermore, three-spined stickleback from different populations responded to parasite manipulation in opposite directions. Such increased ingestion of the intermediate host by preferred and suitable hosts suggests fine-tuned adaptive parasitic manipulation and sheds light on the ongoing evolutionary arms race between hosts and manipulative parasites.

Carry-over effects of early thermal conditions on somatic and germline oxidative damages are mediated by compensatory growth in sticklebacks.

Most studies of climate change impacts focus on the effects of summer temperatures, which can immediately impact fitness of breeders, but winter temperatures are expected to have a greater impact on development and growth of animals with long-lasting consequences. Exposure to warmer temperatures can increase cellular oxidative damage in ectotherms. Yet, it is unknown whether thermal stress during early life has prolonged effects on oxidative status during adulthood. In an experiment using F1 fish originated from a wild three-spined stickleback population at the southern edge of its European distribution, we examined whether experimental thermal conditions experienced in winter had carry-over effects on oxidative status and telomere length, a marker of accumulated stress, in the soma and germline during adulthood. For this, oxidative DNA damage, enzymatic antioxidant activities and telomere length were measured three months after the termination of the temperature manipulation. In addition, we tested whether such delayed effects, if any, were due to individuals' compensatory growth after experiencing unfavourable growth conditions in winter. Warm acclimation during winter induced increased levels of oxidative DNA damage in muscle and sperm and increased enzymatic antioxidant defences in muscle during the breeding season. Telomere length of adult fish was not influenced by thermal conditions experienced during early life. Winter temperature manipulation influenced fish to alter the temporal pattern of growth trajectories across the juvenile and adult stages. Fish reared in warm winter conditions grew at a slower rate than the controls during the period of temperature manipulation then accelerated body mass gain to catch up during the breeding season. Faster somatic growth during the breeding season incurred a higher cost in terms of oxidative damage in the warm-treated individuals. For the first time, we experimentally show the long-lasting detrimental effects of thermal stress on and the positive link between catch-up growth and oxidative DNA damage in the soma and germline. Winter temperature increases due to climate change can reduce fertility and survival of fish by inducing catch-up growth. The detrimental effects of winter climate change may accumulate across generations through the pre-mutagenic DNA damage in the germline.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1. The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB, that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB, suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature.

Evolution of steroid hormones in reproductive females of the threespine stickleback fish

Hormones play a prominent role in animal development, mediating the expression of traits and coordinating phenotypic responses to the environment. Their role as physiological integrators has implications for how populations respond to natural selection and can impact the speed and direction of evolutionary change. However, many emerging and established fish models with the potential to be ecologically or evolutionarily informative are small-bodied, making hormone sampling through traditional methods (whole-body or plasma) lethal or highly disruptive. Sampling methodology has thus restricted study design, often limiting sample sizes, and has prevented the study of at-risk/endangered populations. We utilize water-borne hormone sampling, a minimally invasive method of measuring the rate of steroid hormone release across the gills and further validate this method in a novel, evolutionary context. First, we compare water-borne hormone measures of cortisol with those quantified from plasma and whole-body samples collected from the same individuals to establish the relationship between concentrations quantified via the three methods. We then compare the release of steroid hormones in three populations of threespine stickleback to establish the sensitivity of this tool in measuring within-individual and between-individual variation in biologically relevant contexts (reproductive stages), and in assessing differences among populations with distinct evolutionary histories. We demonstrate a strong positive relationship between cortisol concentrations measured with water-borne, plasma, and whole-body collection techniques. Tracking estradiol and testosterone throughout clutch production in females produced anticipated patterns associated with growing and maturing eggs, with divergence in estradiol production in one population. Additionally, differences among populations in cortisol levels at ovulation paralleled the relative presence of a social stressor, and thus expected energetic needs within each population. We confirm that water-borne hormone sampling is sufficiently sensitive to capture biologically relevant fluctuations in steroid hormones between environmental contexts and demonstrate that among-population differences are detectable. This technique can be applied broadly to small fish to answer important ecological and evolutionary questions. By linking population variation in hormones and the multivariate phenotype, this technique will help elucidate both proximate mechanisms underlying phenotypic development and variation, and the way hormone networks alter evolutionary responses to selection.

Otolith development in wild populations of stickleback: Jones & Hynes method does not apply to most populations.

This paper critiques Jones & Hynes (1950) findings by analysing sequential samples of otoliths from three wild populations of Gasterosteus aculeatus from North Uist, Scotland and Nottingham, England. Contrary to Jones & Hynes (1950), but coincident with the finding of later researchers, our results showed that no true translucent ring formed in the otolith of G. aculeatus during their first summer. The first translucent ring was probably starting to be formed by the end of summer and was completed by the end of their first winter. There was no second opaque ring in the otoliths of G. aculeatus before they passed their first winter. The second opaque ring was just starting to appear by early April of the second year in the southern population i.e. Nottingham, but later, by May, in the northern populations i.e. North Uist. Formation of the opaque ring in G. aculeatus mostly occurs in spring and summer, with younger fish starting earlier. In contrast, the formation of translucent rings is mostly during autumn and winter, but can be more widespread through the year, possibly as a result of reproductive investment.

Convergent evolution of gene expression in two high-toothed stickleback populations.

Changes in developmental gene regulatory networks enable evolved changes in morphology. These changes can be in cis regulatory elements that act in an allele-specific manner, or changes to the overall trans regulatory environment that interacts with cis regulatory sequences. Here we address several questions about the evolution of gene expression accompanying a convergently evolved constructive morphological trait, increases in tooth number in two independently derived freshwater populations of threespine stickleback fish (Gasterosteus aculeatus). Are convergently evolved cis and/or trans changes in gene expression associated with convergently evolved morphological evolution? Do cis or trans regulatory changes contribute more to gene expression changes accompanying an evolved morphological gain trait? Transcriptome data from dental tissue of ancestral low-toothed and two independently derived high-toothed stickleback populations revealed significantly shared gene expression changes that have convergently evolved in the two high-toothed populations. Comparing cis and trans regulatory changes using phased gene expression data from F1 hybrids, we found that trans regulatory changes were predominant and more likely to be shared among both high-toothed populations. In contrast, while cis regulatory changes have evolved in both high-toothed populations, overall these changes were distinct and not shared among high-toothed populations. Together these data suggest that a convergently evolved trait can occur through genetically distinct regulatory changes that converge on similar trans regulatory environments.

Opsin expression predicts male nuptial color in threespine stickleback.

Theoretical models of sexual selection suggest that male courtship signals can evolve through the build‐up of genetic correlations between the male signal and female preference. When preference is mediated via increased sensitivity of the signal characteristics, correlations between male signal and perception/sensitivity are expected. When signal expression is limited to males, we would expect to find signal‐sensitivity correlations in males. Here, we document such a correlation within a breeding population of threespine stickleback mediated by differences in opsin expression. Males with redder nuptial coloration express more long‐wavelength‐sensitive (LWS) opsin, making them more sensitive to orange and red. This correlation is not an artifact of shared tuning to the optical microhabitat. Such correlations are an essential feature of many models of sexual selection, and our results highlight the potential importance of opsin expression variation as a substrate for signal‐preference evolution. Finally, these results suggest a potential sensory mechanism that could drive negative frequency‐dependent selection via male–male competition and thus maintain variation in male nuptial color.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback.

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Adaptation in temporally variable environments: stickleback armor in periodically breaching bar-built estuaries.

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar-built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re-form within a few weeks or months, thereby re-isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within-population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among-population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within-population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within-population variation while reducing, but not eliminating, among-population variation driven by spatial environmental variation.

Liver Protein Comparisons of Warm‐adapted Versus Cold‐adapted Populations of Threespine Sticklebacks (Gasterosteus aculeatus)

This study investigated the molecular phenotypes of a warm‐adapted population of threespine sticklebacks (Gasterosteus aculeatus). We created a liver protein assay consisting of 1,016 proteins with Skyline (MacCoss Lab) from four populations (N=24) using data dependent acquisition (DDA) shotgun proteomics. This assay will be used to perform targeted proteomics on future comparisons, which will allow for more precise identification and quantification of the targeted proteins. A population of sticklebacks from Baja California, Mexico is surviving and breeding at temperatures near 30 degrees Celsius, which is close to the upper thermal limit for this species. Understanding what this population is doing at the molecular level to function at these high temperatures may give insight into mechanisms that have allowed for their survival, as well as the potential for other populations to adapt to higher temperatures. Livers from sticklebacks of the Baja California population were then compared to three cold‐adapted populations using gel‐free quantitative proteomics by liquid chromatography‐mass spectrometry (LC‐MS/MS). A resident marine population (Bodega Harbor, CA), a freshwater population (Lake Solano, CA), and an anadromous population (Anchorage, Alaska) were used in these comparisons. Six Baja California threespine sticklebacks and six fish from each of the other three groups were compared in four separate sets (N=96). Approximately 20 proteins were significantly more abundant in the Baja California population than the other three populations across all four sets. Our results suggest that certain molecular chaperones and redox/antioxidant proteins may contribute to the survival of the Baja California population. Our study lays the foundation for highly accurate quantification of protein differences in future comparisons using targeted proteomics. This study also highlights proteins that are consistently found at higher abundance in the warm‐adapted population and thus may be important for liver adaptation to warmer temperatures.Support or Funding InformationFunded by NSF Grant IOS‐1355098.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Interactive effects of salinity and temperature acclimation on gill morphology and gene expression in threespine stickleback

Colonization of freshwater habitats from marine environments exposes organisms to novel combinations of temperature and salinity, but little is known about physiological responses to the interactive effects of these stressors. Here, we examined the effects of temperature (14 versus 4 °C) and salinity (11 versus 0.3 ppt) on gill gene expression in marine, anadromous, and freshwater populations of threespine stickleback (Gasterosteus aculeatus). Expression of the epithelial calcium channel was not affected by temperature or salinity, but had significantly higher expression in the freshwater ecotype. The combination of low temperature and low salinity had non-additive effects on the expression of the Na+/H+ exchanger. Fish exposed to the combination of low temperature and low salinity had expression levels similar to fish exposed to either factor in isolation. Expression of Na+,K+-ATPase α-subunit was greater in fish exposed to low temperature and low salinity than in fish exposed to the factors separately, and this effect was the most pronounced in the marine ecotype. We also examined the interactive effects of salinity and temperature on gill morphology in the marine ecotype, and observed non-additive effects. Low temperature increased the size of the interlamellar cell mass in fish held at 11 ppt, but not at 0.3 ppt, and the effect of low salinity was in the opposite direction in fish at high and low temperatures. These data demonstrate interactive effects of temperature and salinity and highlight that overwintering in cold freshwater was likely a physiological challenge for marine stickleback as they colonized freshwater following the last glaciation.