Small changes in gene expression of targeted osmoregulatory genes when exposing marine and freshwater threespine stickleback (Gasterosteus aculeatus) to abrupt salinity transfers

Salinity-induced transcriptome profiles in marine and freshwater threespine stickleback after an abrupt 6-hour exposure

Saltwater and freshwater environments have opposing physiological challenges, yet, there are fish species that are able to enter both habitats during short time spans, and as individuals they must therefore adjust quickly to osmoregulatory contrasts. In this study, we conducted an experiment to test for plastic responses to abrupt salinity changes in two populations of threespine stickleback, Gasterosteus aculeatus, representing two ecotypes (freshwater and ancestral saltwater). We exposed both ecotypes to abrupt native (control treatment) and non‐native salinities (0‰ and 30‰) and sampled gill tissue for transcriptomic analyses after 6 h of exposure. To investigate genomic responses to salinity, we analyzed four different comparisons; one for each ecotype (in their control and exposure salinity; (1) and (2), one between ecotypes in their control salinity (3), and the fourth comparison included all transcripts identified in (3) that did not show any expressional changes within ecotype in either the control or the exposed salinity (4)). Abrupt salinity transfer affected the expression of 10 and 1530 transcripts for the saltwater and freshwater ecotype, respectively, and 1314 were differentially expressed between the controls, including 502 that were not affected by salinity within ecotype (fixed expression). In total, these results indicate that factors other than genomic expressional plasticity are important for osmoregulation in stickleback, due to the need for opposite physiological pathways to survive the abrupt change in salinity.

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.

Courtship behavior, nesting microhabitat, and assortative mating in sympatric stickleback species pairs

The maintenance of reproductive isolation in the face of gene flow is a particularly contentious topic, but differences in reproductive behavior may provide the key to explaining this phenomenon. However, we do not yet fully understand how behavior contributes to maintaining species boundaries. How important are behavioral differences during reproduction? To what extent does assortative mating maintain reproductive isolation in recently diverged populations and how important are "magic traits"? Assortative mating can arise as a by-product of accumulated differences between divergent populations as well as an adaptive response to contact between those populations, but this is often overlooked. Here we address these questions using recently described species pairs of three-spined stickleback (Gasterosteus aculeatus), from two separate locations and a phenotypically intermediate allopatric population on the island of North Uist, Scottish Western Isles. We identified stark differences in the preferred nesting substrate and courtship behavior of species pair males. We showed that all males selectively court females of their own ecotype and all females prefer males of the same ecotype, regardless of whether they are from species pairs or allopatric populations. We also showed that mate choice does not appear to be driven by body size differences (a potential "magic trait"). By explicitly comparing the strength of these mating preferences between species pairs and single-ecotype locations, we were able to show that present levels of assortative mating due to direct mate choice are likely a by-product of other adaptations between ecotypes, and not subject to obvious selection in species pairs. Our results suggest that ecological divergence in mating characteristics, particularly nesting microhabitat may be more important than direct mate choice in maintaining reproductive isolation in stickleback species pairs.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Temporal variability in size and growth of Atlantic herring in the Gulf of Maine

Variability in life history traits and structural diversity of commercially exploited fishes in response to stress can impact their population dynamics and sustainability. Using data from a fishery dependent sampling program from 1978 to 2011, we evaluated temporal variability of size and growth of adult Atlantic herring (Clupea harengus) in the Gulf of Maine. We then developed and tested the hypotheses on the links of such temporal changes to population density and environmental factors and found decreases in size and growth potential. Generalized additive models found that density dependence was the main driver of such changes over sea surface temperature and salinity. Our results highlight the importance of density dependent processes in regulating growth and population size structure for Atlantic herring in the Gulf of Maine.

Morphological and physiological traits of Mediterranean sticklebacks living in the Camargue wetland (Rhone river delta)

Three-spined sticklebacks (Gasterosteus aculeatus L.) living at the southern limit of the species distribution range could possess specific morphological and physiological traits that enable these fish to live at the threshold of their physiological capacities. Morphological analysis was carried out on samples of sticklebacks living in different saline habitats of the Camargue area (Rhone delta, northern Mediterranean coast) obtained from 1993 to 2017. Salinity acclimation capacities were also investigated using individuals from freshwater-low salinity drainage canals and from mesohaline-euryhaline lagoons. Fish were maintained in laboratory conditions at salinity values close to those of their respective habitats: low salinity (LS, 5‰) or seawater (SW, 30‰). Fish obtained from a mesohaline brackish water lagoon (BW, 15‰) were acclimated to SW or LS. Oxygen consumption rates and branchial Na⁺ /K⁺ -ATPase (NKA) activity (indicator of fish osmoregulatory capacity) were measured in these LS or SW control fish and in individuals subjected to abrupt SW or LS transfers. At all the studied locations, only the low-plated "leiurus" morphotype showed no spatial or temporal variations in their body morphology. Gill rakers were only longer and denser in fish sampled from the LS-freshwater (FW) drainage canals. All fish presented similar physiological capacities. Oxygen consumption rates were not influenced by salinity challenge except in SW fish transferred to LS immediately and 1 h after transfer. However, and as expected, gill NKA activity was salinity dependent. Sticklebacks of the Camargue area sampled from habitats with contrasted saline conditions are homogenously euryhaline, have low oxygen consumption rates and do not appear to experience significantly greater metabolic costs when challenged with salinity. However, an observed difference in gill raker length and density is most probably related to the nutritional condition of their habitat, indicating that individuals can rapidly acclimatize to different diets.

Draft genome assembly of Tenualosa ilisha, Hilsa shad, provides resource for osmoregulation studies

This study provides the first high-quality draft genome assembly (762.5 Mb) of Tenualosa ilisha that is highly contiguous and nearly complete. We observed a total of 2,864 contigs, with 96.4% completeness with N₅₀ of 2.65 Mbp and the largest contig length of 17.4 Mbp, along with a complete mitochondrial genome of 16,745 bases. A total number of 33,042 protein coding genes were predicted, among these, 512 genes were classified under 61 Gene Ontology (GO) terms, associated with various homeostasis processes. Highest number of genes belongs to cellular calcium ion homeostasis, followed by tissue homeostasis. A total of 97 genes were identified, with 16 GO terms related to water homeostasis. Claudins, Aquaporins, Connexins/Gap junctions, Adenylate cyclase, Solute carriers and Voltage gated potassium channel genes were observed to be higher in number in T. ilisha, as compared to that in other teleost species. Seven novel gene variants, in addition to claudin gene (CLDZ), were found in T. ilisha. The present study also identified two putative novel genes, NKAIN3 and L4AM1, for the first time in fish, for which further studies are required for pinpointing their functions in fish. In addition, 1.6 million simple sequence repeats were mined from draft genome assembly. The study provides a valuable genomic resource for the anadromous Hilsa. It will form a basis for future studies, pertaining to its adaptation mechanisms to different salinity levels during migration, which in turn would facilitate in its domestication.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life‐stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non‐linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life‐stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

Salinity-induced phenotypic plasticity in threespine stickleback sperm activation

Phenotypic expression may be and often is influenced by an organism's developmental environment, referred to as phenotypic plasticity. The sperm cells of teleosts have been found to be inactive in the seminal plasma and are activated by osmotic shock for most fish species, through release in either hypertonic (for marine fish) or hypotonic (for freshwater fish) water. If this is the case, the regulatory system of sperm mobility should be reversed in salt- and freshwater fish. We tested this hypothesis by first activating sperm of salt- and freshwater populations of threespine stickleback in salt- and freshwater. The sperm from saltwater stickleback could be activated in either salinity, which matches the freshwater colonization history of the species, whereas the sperm from the freshwater population acted as predicted by the osmotic shock theory and was activated in freshwater only. As the freshwater population used here was calculated to be thousands of years old, we went on to test whether the trait(s) were plastic and sperm from freshwater males still could be activated in saltwater after individuals were exposed to saltwater. After raising freshwater stickleback in saltwater, we found the mature males to have active sperm in both saltwater and freshwater. Further, we also found the sperm of wild-caught freshwater stickleback to be active in saltwater after exposing those mature males to saltwater for only 2 days. This illustrates that the ability for stickleback sperm to be activated in a range of water qualities is an environmentally induced plastic trait.

RNA sequencing, de novo assembly and differential analysis of the gill transcriptome of freshwater climbing perch Anabas testudineus after 6 days of seawater exposure

To obtain transcriptomic insights into branchial responses to salinity challenge in Anabas testudineus, this study employed RNA sequencing (RNA-Seq) to analyse the gill transcriptome of A. testudineus exposed to seawater (SW) for 6 days compared with the freshwater (FW) control group. A combined FW and SW gill transcriptome was de novo assembled from 169.9 million 101 bp paired-end reads. In silico validation employing 17 A. testudineus Sanger full-length coding sequences showed that 15/17 of them had greater than 80% of their sequences aligned to the de novo assembled contigs where 5/17 had their full-length (100%) aligned and 9/17 had greater than 90% of their sequences aligned. The combined FW and SW gill transcriptome was mapped to 13,780 unique human identifiers at E-value ≤1.0E-20 while 952 and 886 identifiers were determined as up and down-regulated by 1.5 fold, respectively, in the gills of A. testudineus in SW when compared with FW. These genes were found to be associated with at least 23 biological processes. A larger proportion of genes encoding enzymes and transporters associated with molecular transport, energy production, metabolisms were up-regulated, while a larger proportion of genes encoding transmembrane receptors, G-protein coupled receptors, kinases and transcription regulators associated with cell cycle, growth, development, signalling, morphology and gene expression were relatively lower in the gills of A. testudineus in SW when compared with FW. High correlation (R = 0.99) was observed between RNA-Seq data and real-time quantitative PCR validation for 13 selected genes. The transcriptomic sequence information will facilitate development of molecular resources and tools while the findings will provide insights for future studies into branchial iono-osmoregulation and related cellular processes in A. testudineus.

Gene and Blood Analysis Reveal That Transfer from Brackish Water to Freshwater Is More Stressful to the Silverside Odontesthes humensis

Silversides are fish that inhabit marine coastal waters, coastal lagoons, and estuarine regions in southern South America. The freshwater (FW) silversides have the ability to tolerate salinity variations. Odontesthes humensis have similar habitats and biological characteristics of congeneric O. bonariensis, the most studied silverside species and with great economic importance. Studies revealed that O. bonariensis is not fully adapted to FW, despite inhabiting hyposmotic environments in nature. However, there is little information about stressful environments for cultivation of silverside O. humensis. Thus, the aim of this study was to evaluate the stress and osmoregulation responses triggered by the osmotic transfers on silverside O. humensis. Silversides were acclimated to FW (0 ppt) and to brackish water (BW, 10 ppt) and then they were exposed to opposite salinity treatment. Silverside gills and blood were sampled on pre-transfer (D0) and 1, 7, and 15 days (D1, D7, and D15) after changes in environmental salinity, the expression levels of genes atp1a3a, slc12a2b, kcnh1, and hspa1a were determined by quantitative reverse transcription-PCR for evaluation of osmoregulatory and stress responses. Furthermore, glycemia, hematocrit, and osmolality were also evaluated. The expression of atp1a3a was up- and down-regulated at D1 after the FW-BW and BW-FW transfers, respectively. Slc12a2b was up-regulated after FW-BW transfer. Similarly, kcnh1 and hspa1a were up-regulated at D1 after the BW-FW transfer. O. humensis blood osmolality decreased after the exposure to FW. It remained stable after exposure to BW, indicating an efficient hyposmoregulation. The glycemia had a peak at D1 after BW-FW transfer. No changes were observed in hematocrit. The return to the pre-transfer levels at D7 after the significant increases in responses of almost all evaluated molecular and blood parameters indicated that this period is enough for acclimation to the experimental conditions. In conclusion, our results suggest that BW-FW transfer is more stressful to O. humensis than FW-BW transfer and the physiology of O. humensis is only partially adapted to FW.

Ongoing niche differentiation under high gene flow in a polymorphic brackish water threespine stickleback (Gasterosteus aculeatus) population

BACKGROUND

Marine threespine sticklebacks colonized and adapted to brackish and freshwater environments since the last Pleistocene glacial. Throughout the Holarctic, three lateral plate morphs are observed; the low, partial and completely plated morph. We test if the three plate morphs in the brackish water Lake Engervann, Norway, differ in body size, trophic morphology (gill raker number and length), niche (stable isotopes; δ15N, δ13C, and parasites (Theristina gasterostei, Trematoda spp.)), genetic structure (microsatellites) and the lateral-plate encoding Stn382 (Ectodysplasin) gene. We examine differences temporally (autumn 2006/spring 2007) and spatially (upper/lower sections of the lake - reflecting low versus high salinity).

RESULTS

All morphs belonged to one gene pool. The complete morph was larger than the low plated, with the partial morph intermediate. The number of lateral plates ranged 8-71, with means of 64.2 for complete, 40.3 for partial, and 14.9 for low plated morph. Stickleback δ15N was higher in the lower lake section, while δ13C was higher in the upper section. Stickleback isotopic values were greater in autumn. The low plated morph had larger variances in δ15N and δ13C than the other morphs. Sticklebacks in the upper section had more T. gasterostei than in the lower section which had more Trematoda spp. Sticklebacks had less T. gasterostei, but more Trematoda spp. in autumn than spring. Sticklebacks with few and short rakers had more T. gasterostei, while sticklebacks with longer rakers had more Trematoda. spp. Stickleback with higher δ15N values had more T. gasterostei, while sticklebacks with higher δ15N and δ13C values had more Trematoda spp. The low plated morph had fewer Trematoda spp. than other morphs.

CONCLUSIONS

Trait-ecology associations may imply that the three lateral plate morphs in the brackish water lagoon of Lake Engervann are experiencing ongoing divergent selection for niche and migratory life history strategies under high gene flow. As such, the brackish water zone may generally act as a generator of genomic diversity to be selected upon in the different environments where threespine sticklebacks can live.

Selection on the morphology-physiology-performance nexus: Lessons from freshwater stickleback morphs

Conspecifics inhabiting divergent environments frequently differ in morphology, physiology, and performance, but the interrelationships amongst traits and with Darwinian fitness remains poorly understood. We investigated population differentiation in morphology, metabolic rate, and swimming performance in three-spined sticklebacks (Gasterosteus aculeatus L.), contrasting a marine/ancestral population with two distinct freshwater morphotypes derived from it: the "typical" low-plated morph, and a unique "small-plated" morph. We test the hypothesis that similar to plate loss in other freshwater populations, reduction in lateral plate size also evolved in response to selection. Additionally, we test how morphology, physiology, and performance have evolved in concert as a response to differences in selection between marine and freshwater environments. We raised pure-bred second-generation fish originating from three populations and quantified their lateral plate coverage, burst- and critical swimming speeds, as well as standard and active metabolic rates. Using a multivariate QST-FST framework, we detected signals of directional selection on metabolic physiology and lateral plate coverage, notably demonstrating that selection is responsible for the reduction in lateral plate coverage in a small-plated stickleback population. We also uncovered signals of multivariate selection amongst all bivariate trait combinations except the two metrics of swimming performance. Divergence between the freshwater and marine populations exceeded neutral expectation in morphology and in most physiological and performance traits, indicating that adaptation to freshwater habitats has occurred, but through different combinations of traits in different populations. These results highlight both the complex interplay between morphology, physiology and performance in local adaptation, and a framework for their investigation.

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.

Duplicated CFTR isoforms in eels diverged in regulatory structures and osmoregulatory functions

Two cystic fibrosis transmembrane conductance regulator (CFTR) isoforms, CFTRa and CFTRb, were cloned in Japanese eel and their structures and functions were studied in different osmoregulatory tissues in freshwater (FW) and seawater (SW) eels. Molecular phylogenetic results suggested that the CFTR duplication in eels occurred independently of the duplication event in salmonid. CFTRa was expressed in the intestine and kidney and downregulated in both tissues in SW eels, while CFTRb was specifically expressed in the gill and greatly upregulated in SW eels. Structurally, the CFTR isoforms are similar in most functional domains except the regulatory R domain, where the R domain of CFTRa is similar to that of human CFTR but the R domain of CFTRb is unique in having high intrinsic negative charges and fewer phosphorylation sites, suggesting divergence of isoforms in terms of gating properties and hormonal regulation. Immunohistochemical results showed that CFTR was localized on the apical regions of SW ionocytes, suggesting a Cl(-) secretory role as in other teleosts. In intestine and kidney, however, immunoreactive CFTR was mostly found in the cytosolic vesicles in FW eels, indicating that Cl(-) channel activity could be low at basal conditions, but could be rapidly increased by membrane insertion of the stored channels. Guanylin (GN), a known hormone that increases CFTR activity in mammalian intestine, failed to redistribute CFTR and to affect its expression in eel intestine. The results suggested that GN-independent CFTR regulation is present in eel intestine and kidney.

Flexible selection of diversified Na(+)/K(+)-ATPase α-subunit isoforms for osmoregulation in teleosts

BACKGROUND AND METHODS

Multiple Na+/K+-ATPase (NKA) α-subunit isoforms express differentially in response to salinity transfer in teleosts but we observed that the isoform nomenclature is inconsistent with the phylogenetic relationship of NKA α-genes. We cloned the catalytic NKA α-subunit isoforms in eels and medaka, analyzed the time course of their expressions in osmoregulatory tissues after transfer from freshwater (FW) to seawater (SW), and performed phylogenetic analyses to deduce an evolutionary scenario that illustrates how various duplication events have led to the current genomic arrangement of NKA α-genes in teleosts.

RESULTS AND DISCUSSION

Five and six α-subunits were cloned in eels and medaka respectively. In eels, the commonly-reported α1a and α1b isoforms were absent while the α1c isoform was diversified instead (α1c-1, α1c-2, α1c-3, α2, and α3 in eels). Phylogenetic estimation indicated that the α1a and α1b isoforms from salmon, tilapia, and medaka were generated by independent duplication events and thus they are paralogous isoforms. Re-examination of expression changes of known isoforms after salinity challenge revealed that the isoforms selected as predominant SW-types varied among teleost lineages. Diversification of α1 isoforms occurred by various types of gene duplication, or by alternative transcription among tandem genes to form chimeric transcripts, but there is no trend for more α1 copies in euryhaline species. Our data suggest that the isoform switching between FW (α1a predominates) and SW (α1b predominates) that occurs in salmonids is not universal in teleosts. Instead, in eels, α1c-1 was the major α-subunit upregulated gill, intestine, and kidney in SW. Localization of both NKA mRNA and protein showed consistent upregulation in gill and intestine in SW eels, but not in renal distal and collecting tubules, where low transcript expression levels were accompanied by high protein levels, suggesting a tissue-specific translational regulation that determines and fine-tunes the NKA expression. In medaka, α1b was upregulated in SW in anterior intestine while most other α-subunit isoforms were less responsive to salinity changes.

CONCLUSION

By integrating gene expression and phylogenetic results, we propose that the major NKA α-subunits for SW acclimation were not ancestrally selected, but rather were flexibly determined in lineage-specific fashion in teleosts.