Interpopulation differences in expression of candidate genes for salinity tolerance in winter migrating anadromous brown trout (Salmo trutta L.)

The influence of HSP inducers on salinity stress in sterlet sturgeon (Acipenser ruthenus): In vitro study on HSP expression, immune responses, and antioxidant capacity

Heat shock proteins (HSPs) play a crucial role in antioxidant systems, immune responses, and enzyme activation during stress conditions. Salinity changes can cause stress and energy expenditure in fish, resulting in mortality, especially in fingerlings. The purpose of this study was to examine the relationship between salinity and HSPs in stressed fish by assessing the effects of various HSP inducers (HSPis), including Pro-Tex® (800 mM), amygdalin (80 mM), and a novel synthetic compound derived from pirano piranazole (80 µM), on isolated cells from Sterlet Sturgeon (Acipenser ruthenus) exposed to 13 ‰ salinity (S13). After liver, kidney, and gill cells were cultured, the HSPi compounds were treated in vitro in the presence and absence of salinity. The expression patterns of HSP27, HSP70, and HSP90 were assessed by Western blotting. Biochemical enzymes (aspartate aminotransferase, alanine aminotransferase, alkaline phosphatase, and lactate dehydrogenase), cortisol levels, and immune parameters (component 3, immunoglobulin M, and lysozyme) were measured before and after treatment with HSPis and HSPi + S13. According to these findings, HSPis positively modulate HSP expression, immune responses, and antioxidant levels. Furthermore, they increased in vitro cell survival by maintaining cortisol levels and biochemical enzyme activities in A. ruthenus under saline conditions (P < 0.0001). In conclusion, HSPis can increase A. ruthenus resistance to salinity stress. However, the results also indicated that these compounds can reverse the adverse effects of salinity. The effectiveness of this approach depends on further research into the effects of these ecological factors on the health status of the species, especially in vivo and in combination with other stresses.

The ecology and fishery of the vendace (Coregonus albula) in the Baltic Sea

Brackish water ecosystems often have high primary production, intermediate salinities, and fluctuating physical conditions and therefore provide challenging environments for many of their inhabitants. This is especially true of the Baltic Sea, which is a large body of brackish water under strong anthropogenic influence. One freshwater species that is able to cope under these conditions in the northern Baltic Sea is the vendace (Coregonus albula), a small salmonid fish. Here, we review the current knowledge of its ecology and fishery in this brackish water environment. The literature shows that, by competing for resources with other planktivores and being an important prey for a range of larger species, C. albula plays a notable role in the northern Baltic Sea ecosystem. It also sustains significant fisheries in the coastal waters of Sweden and Finland. We identify the need to better understand these C. albula populations in terms of the predator-prey interactions, distributions of anadromous and sea spawning populations and other putative (eco)morphs, extent of gene exchange between the populations, and effects of climate change on their future. In this regard, we recommend strengthening C. albula-related research and management efforts by improved collaboration and coordination between research institutions, other governmental agencies, and fishers, as well as by harmonization of fishery policies across national borders.

Biological responses of stellate sturgeon fingerlings (Acipenser stellatus) immersed in HSP inducer to salinity changes

Changes in salinity is a stressful and energy-consuming process in fish which give rise to mortalities, especially in fish fingerlings that are more sensitive during the early stages of their life. In the present study, the effects of three salinities, 3‰ (downstream of river), 8‰ (estuarine), and 13‰ (the maximum salinity in the Caspian Sea), on HSP70 gene expression, cortisol level, immune response (lysozyme, complement C3, IgM), and antioxidant enzyme activities (SOD, CAT, T-AOC) of the stellate sturgeon fingerlings in the presence of HSP inducer compound (TEX-OE®) were evaluated. Our results showed that levels of plasma cortisol and heat shock protein (HSP70) in Acipenser stellatus fingerlings increased due to salinity changes. In the presence of the HSP inducer, HSP70 expression in both gill and liver was significantly increased, whereas cortisol level was notably decreased. Exposure to salinity changes resulted in an increase in antioxidant defense activities (SOD, CAT, and T-AOC) and immune response (lysozyme, IgM, and C3) in the presence of an HSP inducer. In conclusion, an HSP-inducing compounds can have a positive effect in strengthening the immunity and antioxidant system of sturgeon fingerlings by increasing the expression of the HSP70 gene against salinity fluctuations and generally increase the body's physiological tolerance.

Sea trout (Salmo trutta) straying rate decreases as distance from river mouth increases

For both conservation and management purposes, it is important to identify the natal origin of migratory individuals entering a river, particularly in genetically spatially structured species like brown trout (Salmo trutta) where the migrant ecotype (called sea trout) can originate from different populations. Nonetheless, little attention has been paid to the spatial distribution of non-local fish at the level of an entire river catchment. The objective of the study was to quantify the proportion of non-local sea trout entering a river catchment (i.e., straying rate) and estimate the spatial extent of their upstream migration. Here, the authors considered dispersal in three distinct rivers, taking advantage of 10 sampling sites. Sea trout, either trapped or rod-caught, were genotyped and genetically assigned to their source populations using appropriate baselines. Based on 1437 sea trout fish classified as local or non-local, the authors empirically demonstrate that straying rate declines in each river as distance from the coast increases in a non-linear fashion. Straying rate exceeds 50% near the mouth, and then decreases gradually to reach <10% 40-50 km inland. A similar spatial pattern is found in the three rivers investigated suggesting an underlying common behaviour of non-local sea trout. The data and results presented here suggest that straying in is far more constrained spatially than first expected. The majority of non-local sea trout were found within the first 25 km of the estuary in the three rivers investigated.

Salinity significantly affects intestinal microbiota and gene expression in striped catfish juveniles

In the present study, juvenile striped catfish (Pangasianodon hypophthalmus), a freshwater fish species, have been chronically exposed to a salinity gradient from freshwater to 20 psu (practical salinity unit) and were sampled at the beginning (D20) and the end (D34) of exposure. The results revealed that the intestinal microbial profile of striped catfish reared in freshwater conditions were dominated by the phyla Bacteroidetes, Firmicutes, Proteobacteria, and Verrucomicrobia. Alpha diversity measures (observed OTUs (operational taxonomic units), Shannon and Faith's PD (phylogenetic diversity)) showed a decreasing pattern as the salinities increased, except for the phylogenetic diversity at D34, which was showing an opposite trend. Furthermore, the beta diversity between groups was significantly different. Vibrio and Akkermansia genera were affected differentially with increasing salinity, the former being increased while the latter was decreased. The genus Sulfurospirillium was found predominantly in fish submitted to salinity treatments. Regarding the host response, the fish intestine likely contributed to osmoregulation by modifying the expression of osmoregulatory genes such as nka1a, nka1b, slc12a1, slc12a2, cftr, and aqp1, especially in fish exposed to 15 and 20 psu. The expression of heat shock proteins (hsp) hsp60, hsp70, and hsp90 was significantly increased in fish reared in 15 and 20 psu. On the other hand, the expression of pattern recognition receptors (PRRs) were inhibited in fish exposed to 20 psu at D20. In conclusion, the fish intestinal microbiota was significantly disrupted in salinities higher than 10 psu and these effects were proportional to the exposure time. In addition, the modifications of intestinal gene expression related to ion exchange and stressful responses may help the fish to adapt hyperosmotic environment. KEY POINTS: • It is the first study to provide detailed information on the gut microbiota of fish using the amplicon sequencing method. • Salinity environment significantly modified the intestinal microbiota of striped catfish. • Intestinal responses may help the fish adapt to hyperosmotic environment.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

mRNA profile provides novel insights into stress adaptation in mud crab megalopa, Scylla paramamosain after salinity stress

Background

Mud crab, Scylla paramamosain, a euryhaline crustacean species, mainly inhabits the Indo-Western Pacific region. Wild mud crab spawn in high-salt condition and the salinity reduced with the growth of the hatching larvae. When the larvae grow up to megalopa, they migrate back to estuaries and coasts in virtue of the flood tide, settle and recruit adult habitats and metamorphose into the crablet stage. Adult crab can even survive in a wide salinity of 0–35 ppt. To investigate the mRNA profile after salinity stress, S. paramamosain megalopa were exposed to different salinity seawater (low, 14 ppt; control, 25 ppt; high, 39 ppt).

Results

Firstly, from the expression profiles of Na+/K+/2Cl- cotransporter, chloride channel protein 2, and ABC transporter, it turned out that the 24 h might be the most influenced duration in the short-term stress. We collected megalopa under different salinity for 24 h and then submitted to mRNA profiling. Totally, 57.87 Gb Clean Data were obtained. The comparative genomic analysis detected 342 differentially expressed genes (DEGs). The most significantly DEGs include gamma-butyrobetaine dioxygenase-like, facilitated trehalose transporter Tret1, sodium/potassium-transporting ATPase subunit alpha, rhodanese 1-like protein, etc. And the significantly enriched pathways were lysine degradation, choline metabolism in cancer, phospholipase D signaling pathway, Fc gamma R-mediated phagocytosis, and sphingolipid signaling pathway. The results indicate that in the short-term salinity stress, the megalopa might regulate some mechanism such as metabolism, immunity responses, osmoregulation to adapt to the alteration of the environment.

Conclusions

This study represents the first genome-wide transcriptome analysis of S. paramamosain megalopa for studying its stress adaption mechanisms under different salinity. The results reveal numbers of genes modified by salinity stress and some important pathways, which will provide valuable resources for discovering the molecular basis of salinity stress adaptation of S. paramamosain larvae and further boost the understanding of the potential molecular mechanisms of salinity stress adaptation for crustacean species.

Intra-Specific Difference in the Effect of Salinity on Physiological Performance in European Perch (Perca fluviatilis) and Its Ecological Importance for Fish in Estuaries

Changes in environmental salinity challenge fish homeostasis and may affect physiological performance, such as swimming capacity and metabolism, which are important for foraging, migration, and escaping predators in the wild. The effects of salinity stress on physiological performance are largely species specific, but may also depend on intra-specific differences in physiological capabilities of sub-populations. We measured critical swimming speed (U_crit) and metabolic rates during swimming and at rest at salinities of 0 and 10 in European perch (Perca fluviatilis) from a low salinity tolerance population (LSTP) and a high salinity tolerance population (HSTP). U_crit of LSTP was significantly reduced at a salinity of 10 yet was unaffected by salinity change in HSTP. We did not detect a significant cost of osmoregulation, which should theoretically be apparent from the metabolic rates during swimming and at rest at a salinity of 0 compared to at a salinity of 10 (iso-osmotic). Maximum metabolic rates were also not affected by salinity, indicating a modest tradeoff between respiration and osmoregulation (osmo-respiratory compromise). Intra-specific differences in effects of salinity on physiological performance are important for fish species to maintain ecological compatibility in estuarine environments, yet render these sub-populations vulnerable to fisheries. The findings of the present study are therefore valuable knowledge in conservation and management of estuarine fish populations.

Association Mapping Based on a Common-Garden Migration Experiment Reveals Candidate Genes for Migration Tendency in Brown Trout

A better understanding of the environmental and genetic contribution to migratory behavior and the evolution of traits linked to migration is crucial for fish conservation and fisheries management. Up to date, a few genes with unequivocal influence on the adoption of alternative migration strategies have been identified in salmonids. Here, we used a common garden set-up to measure individual migration distances of generally highly polymorphic brown trout Salmo trutta from two populations. Fish from the assumedly resident population showed clearly shorter migration distances than the fish from the assumed migratory population at the ages of 2 and 3 years. By using two alternative analytical pipelines with 22186 and 18264 SNPs obtained through RAD-sequencing, we searched for associations between individual migration distance, and both called genotypes and genotype probabilities. None of the SNPs showed statistically significant individual effects on migration after correction for multiple testing. By choosing a less stringent threshold, defined as an overlap of the top 0.1% SNPs identified by the analytical pipelines, GAPIT and Angsd, we identified eight candidate genes that are potentially linked to individual migration distance. While our results demonstrate large individual and population level differences in migration distances, the detected genetic associations were weak suggesting that migration traits likely have multigenic control.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Proteomic Analysis of the Chlorophyta Dunaliella New Strain AL-1 Revealed Global Changes of Metabolism during High Carotenoid Production

The green microalgae Dunaliella genus is known for the production of high added value molecules. In this study, strain AL-1 was isolated from the Sebkha of Sidi El Hani (Sousse, Tunisia). This isolate was identified both morphologically and genetically via 18S rRNA gene sequence as a member of the genus Dunaliella. Strain AL-1 was found to be closely related to Dunaliella salina, Dunaliella quartolecta and Dunaliella polymorpha with more than 97% similarity. Response surface methodology was used to maximize carotenoid production by strain AL-1 by optimizing its growth conditions. The highest carotenoid content was obtained at salinity: 51, light intensity: 189.89 μmol photons·m^-2·s^-1, and nitrogen: 60 mg·L^-1. Proteomic profiling, using two-dimensional gel electrophoresis, was performed from standard and optimized cultures. We detected 127 protein spots which were significantly differentially expressed between standard and optimized cultures. Among them 16 protein spots were identified with mass spectrometry and grouped into different functional categories using KEGG (Kyoto Encyclopedia of Genes and Genomes) such as photosynthetic Calvin cycle, regulation/defense, energy metabolism, glycolysis, and cellular processes. The current study could be of great interest in providing information on the effect of stressful conditions in microalgae carotenoid production.

Arctic grayling (Thymallus arcticus) in saltwater: a response to Blair et al. (2016)

Although not well known, Arctic grayling can move through saline waters and are captured regularly in nearshore coastal waters in Arctic Canada and Alaska with salinities up to 18 ppt. We highlight the implications this has for Blair et al. (2016), a paper recently published in Conservation Physiology.

Fungal and cyanobacterial gene expression in a lichen symbiosis: Effect of temperature and location

Organisms have evolved different cellular mechanisms to deal with environmental stress, primarily through complex molecular mechanisms including protein refolding and DNA repair. As mutualistic symbioses, lichens offer the possibility of analyzing molecular stress responses in a particularly tight interspecific relationship. We study the widespread cyanolichen Peltigera membranacea, a key player in carbon and nitrogen cycling in terrestrial ecosystems at northern latitudes. We ask whether increased temperature is reflected in mRNA levels of selected damage control genes, and do the response patterns show geographical associations? Using real-time PCR quantification of 38 transcripts, differential expression was demonstrated for nine cyanobacterial and nine fungal stress response genes (plus the fungal symbiosis-related lec2 gene) when the temperature was increased from 5 °C to 15 °C and 25 °C. Principle component analysis (PCA) revealed two gene groups with different response patterns. Whereas a set of cyanobacterial DNA repair genes and the fungal lec2 (PC1 group) showed an expression drop at 15 °C vs. 5 °C, most fungal candidates (PC2 group) showed increased expression at 25 °C vs. 5 °C. PC1 responses also correlated with elevation. The correlated downregulation of lec2 and cyanobacterial DNA repair genes suggests a possible interplay between the symbionts warranting further studies.

Reduced salinity tolerance in the Arctic grayling (Thymallus arcticus) is associated with rapid development of a gill interlamellar cell mass: implications of high-saline spills on native freshwater salmonids

Arctic grayling (Thymallus arcticus) are salmonids that have a strict freshwater existence in post-glacial North America. Oil and gas development is associated with production of high volumes of hypersaline water. With planned industrial expansion into northern areas of Canada and the USA that directly overlap grayling habitat, the threat of accidental saline water release poses a significant risk. Despite this, we understand little about the responses of grayling to hypersaline waters. We compared the physiological responses and survivability of Arctic grayling and rainbow trout (Oncorhynchus mykiss) to tolerate an acute transfer to higher saline waters. Arctic grayling and rainbow trout were placed directly into 17 ppt salinity and sampled at 24 and 96 h along with control animals in freshwater at 24 h. Serum sodium, chloride and osmolality levels increased significantly in grayling at both 24 and 96 h time points, whereas trout were able to compensate for the osmoregulatory disturbance by 96 h. Sodium-potassium ATPase mRNA expression responses to salinity were also compared, demonstrating the inability of the grayling to up-regulate the seawater isoform nkaα1b. Our results demonstrated a substantially lower salinity tolerance in grayling. We also found a significant salinity-induced morphological gill remodelling by Arctic grayling, as demonstrated by the rapid growth of an interlamellar cell mass by 24 h that persisted at 96 h. We visualized and quantified the appearance of the interlamellar cell mass as a response to high salinity, although the functional significance remains to be understood fully. Compared with rainbow trout, which are used as an environmental regulatory species, Arctic grayling are unable to compensate for the osmotic stressors that would result from a highly saline produced water spill. Given these new data, collaboration between fisheries and the oil and gas industry will be vital in the long-term conservation strategies with regard to the Arctic grayling in their native habitat.

Identification of potential markers and sensitive tissues for low or high salinity stress in an intertidal mud crab (Macrophthalmus japonicus)

Macrophthalmus japonicus is an intertidal mud crab is an ecologically important species in Korea, can tolerate a wide range of natural and anthropogenic stressors. Environmental changes especially salinity cause physiological stress to the marine habitats. Differential gene transcription of M. japonicus tissues provided information about tissue specific responses against salinity. Five potential genes were identified and their transcription levels were determined quantitatively comparison to seawater (SW: 31 ± 1 psu) in M. japonicus gills and hepatopancreas after exposed them to different salinities. Ecdysteroid receptor (Mj-EcR), trypsin (Mj-Tryp), arginine kinase (Mj-AK), lipopolysaccharide and β-1,3-glucan binding protein (Mj-LGBP) and peroxinectin (Mj-Prx) in hepatopancreas up-regulated against different salinities. In contrast, the gills, Mj-EcR, Mj-Tryp and Mj-AK showed late up-regulated responses to 40 psu compared to SW. All genes except Mj-LGBP showed up regulation in the gills as time dependent manner. These genes can be considered as potential markers to assess responses in salinity changes. This study suggests hepatopancreas is a suitable tissue for transcriptional, biochemical and physiological responses analysis on M. japonicus in low and high salinity stress.

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

An integrated transcriptomic and comparative genomic analysis of differential gene expression in Arctic charr (Salvelinus alpinus) following seawater exposure

High-throughput RNA sequencing was employed to compare expression profiles in two Arctic charr (Salvelinus alpinus) families post seawater exposure to identify genes and biological processes involved in hypo-osmoregulation and regulation of salinity tolerance. To further understand the genetic architecture of hypo-osmoregulation, the genomic organization of differentially expressed (DE) genes was also analysed. Using a de novo gill transcriptome assembly we found over 2300 contigs to be DE. Major transporters from the seawater mitochondrion-rich cell (MRC) complex were up-regulated in seawater. Expression ratios for 257 differentially expressed contigs were highly correlated between families, suggesting they are strictly regulated. Based on expression profiles and known molecular pathways we inferred that seawater exposure induced changes in methylation states and elevated peroxynitrite formation in gill. We hypothesized that concomitance between DE immune genes and the transition to a hypo-osmoregulatory state could be related to Cl- sequestration by antimicrobial defence mechanisms. Gene Ontology analysis revealed that cell division genes were up-regulated, which could reflect the proliferation of ATP1α1b-type seawater MRCs. Comparative genomics analyses suggest that hypo-osmoregulation is influenced by the relative proximities among a contingent of genes on Arctic charr linkage groups AC-4 and AC-12 that exhibit homologous affinities with a region on stickleback chromosome Ga-I. This supports the hypothesis that relative gene location along a chromosome is a property of the genetic architecture of hypo-osmoregulation. Evidence of non-random structure between hypo-osmoregulation candidate genes was found on AC-1/11 and AC-28, suggesting that interchromosomal rearrangements played a role in the evolution of hypo-osmoregulation in Arctic charr.

Transcriptomics of salinity tolerance capacity in Arctic charr (Salvelinus alpinus): a comparison of gene expression profiles between divergent QTL genotypes

Osmoregulatory capabilities have played an important role in the evolution, dispersal, and diversification of vertebrates. To better understand the genetic architecture of hypo-osmoregulation in fishes and to determine which genes and biological processes affect intraspecific variation in salinity tolerance, we used mRNA sequence libraries from Arctic charr gill tissue to compare gene expression profiles in fish exhibiting divergent salinity tolerance quantitative trait locus (QTL) genotypes. We compared differentially expressed genes with QTL positions to gain insight about the nature of the underlying polymorphisms and examined gene expression within the context of genome organization to gain insight about the evolution of hypo-osmoregulation in fishes. mRNA sequencing of 18 gill tissue libraries produced 417 million reads, and the final reduced de novo transcriptome assembly consisted of 92,543 contigs. Families contained a similar number of differentially expressed contigs between high and low salinity tolerance capacity groups, and log2 expression ratios ranged from 10.4 to -8.6. We found that intraspecific variation in salinity tolerance capacity correlated with differential expression of immune response genes. Some differentially expressed genes formed clusters along linkage groups. Most clusters comprised gene pairs, though clusters of three, four, and eight genes were also observed. We postulated that conserved synteny of gene clusters on multiple ancestral and teleost chromosomes may have been preserved via purifying selection. Colocalization of QTL with differentially expressed genes suggests that polymorphisms in cis-regulatory elements are part of a majority of QTL.

Whole genome scanning and association mapping identified a significant association between growth and a SNP in the IFABP-a gene of the Asian seabass

BACKGROUND

Aquaculture is the quickest growing sector in agriculture. However, QTL for important traits have been only identified in a few aquaculture species. We conducted QTL mapping for growth traits in an Asian seabass F(2) family with 359 individuals using 123 microsatellites and 22 SNPs, and performed association mapping in four populations with 881 individuals.

RESULTS

Twelve and nine significant QTL, as well as 14 and 10 suggestive QTL were detected for growth traits at six and nine months post hatch, respectively. These QTL explained 0.9-12.0% of the phenotypic variance. For body weight, two QTL intervals at two stages were overlapped while the others were mapped onto different positions. The IFABP-a gene located in a significant QTL interval for growth on LG5 was cloned and characterized. A SNP in exon 3 of the gene was significantly associated with growth traits in different populations.

CONCLUSIONS

The results of QTL mapping for growth traits suggest that growth at different stages was controlled by some common QTL and some different QTL. Positional candidate genes and association mapping suggest that the IFABP-a is a strong candidate gene for growth. Our data supply a basis for fine mapping QTL, marker-assisted selection and further detailed analysis of the functions of the IFABP-a gene in fish growth.

Identification and characterization of differentially expressed transcripts in the gills of freshwater prawn (Macrobrachium rosenbergii) under salt stress

The giant freshwater prawn, Macrobrachium rosenbergii, is an economically important species. It is a euryhaline shrimp, surviving in wide-range salinity conditions. A change in gene expression has been suggested as an important component for stress management. To better understand the osmoregulatory mechanisms mediated by the gill, a subtractive and suppressive hybridization (SSH) tool was used to identify expressed transcripts linked to adaptations in saline water. A total of 117 transcripts represented potentially expressed under salinity conditions. BLAST analysis identified 22% as known genes, 9% as uncharacterized showing homologous to unannotated ESTs, and 69% as unknown sequences. All the identified known genes representing broad spectrum of biological pathways were particularly linked to stress tolerance including salinity tolerance. Expression analysis of 10 known genes and 7 unknown/uncharacterized genes suggested their upregulation in the gills of prawn exposed to saline water as compared to control indicating that these are likely to be associated with salinity acclimation. Rapid amplification of cDNA ends (RACE) was used for obtaining full-length cDNA of MRSW-40 clone that was highly upregulated during salt exposure. The sequenced ESTs presented here will have potential implications for future understanding about salinity acclimation and/or tolerance of the prawn.

Differences in salinity tolerance and gene expression between two populations of Atlantic cod (Gadus morhua) in response to salinity stress

Populations of marine fish, even from contrasting habitats, generally show low genetic differentiation at neutral genetic markers. Nevertheless, there is increasing evidence for differences in gene expression among populations that may be ascribed to adaptive divergence. Studying variation in salinity tolerance and gene expression among Atlantic cod (Gadus morhua) from two populations distributed across a steep salinity gradient, we observed high mortality (45% North Sea cod and 80% Baltic Sea cod) in a reciprocal common garden setup. Quantitative RT-PCR assays for expression of hsp70 and Na/K-ATPase α genes demonstrated significant differences in gene regulation within and between populations and treatment groups despite low sample sizes. Most interesting are the significant differences observed in expression of the Na/K-ATPase α gene in gill tissue between North Sea and Baltic cod. The findings strongly suggest that Atlantic cod are adapted to local saline conditions, despite relatively low levels of neutral genetic divergence between populations.

The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus)

BACKGROUND

The capacity to maintain internal ion homeostasis amidst changing conditions is particularly important for teleost fishes whose reproductive cycle is dependent upon movement from freshwater to seawater. Although the physiology of seawater osmoregulation in mitochondria-rich cells of fish gill epithelium is well understood, less is known about the underlying causes of inter- and intraspecific variation in salinity tolerance. We used a genome-scan approach in Arctic charr (Salvelinus alpinus) to map quantitative trait loci (QTL) correlated with variation in four salinity tolerance performance traits and six body size traits. Comparative genomics approaches allowed us to infer whether allelic variation at candidate gene loci (e.g., ATP1α1b, NKCC1, CFTR, and cldn10e) could have underlain observed variation.

RESULTS

Combined parental analyses yielded genome-wide significant QTL on linkage groups 8, 14 and 20 for salinity tolerance performance traits, and on 1, 19, 20 and 28 for body size traits. Several QTL exhibited chromosome-wide significance. Among the salinity tolerance performance QTL, trait co-localizations occurred on chromosomes 1, 4, 7, 18 and 20, while the greatest experimental variation was explained by QTL on chromosomes 20 (19.9%), 19 (14.2%), 4 (14.1%) and 12 (13.1%). Several QTL localized to linkage groups exhibiting homeologous affinities, and multiple QTL mapped to regions homologous with the positions of candidate gene loci in other teleosts. There was no gene × environment interaction among body size QTL and ambient salinity.

CONCLUSIONS

Variation in salinity tolerance capacity can be mapped to a subset of Arctic charr genomic regions that significantly influence performance in a seawater environment. The detection of QTL on linkage group 12 was consistent with the hypothesis that variation in salinity tolerance may be affected by allelic variation at the ATP1α1b locus. IGF2 may also affect salinity tolerance capacity as suggested by a genome-wide QTL on linkage group 19. The detection of salinity tolerance QTL in homeologous regions suggests that candidate loci duplicated from the salmonid-specific whole-genome duplication may have retained their function on both sets of homeologous chromosomes. Homologous affinities suggest that loci affecting salinity tolerance in Arctic charr may coincide with QTL for smoltification and salinity tolerance traits in rainbow trout. The effects of body size QTL appear to be independent of changes in ambient salinity.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.

Differential expression of the heat shock protein Hsp70 in natural populations of the tilapia, Sarotherodon melanotheron, acclimatised to a range of environmental salinities

Background

The relationship between environmental variation and induction of heat shock proteins (Hsps) has been much documented under experimental conditions. However, very little is known about such induction in natural populations acclimatised to prevailing environmental conditions. Furthermore, while induction of stress proteins has been well documented in response to environmental contaminants and thermal stressors, little is known about whether factors, such as extreme salinity, are also potential inductors. The black-chinned tilapia Sarotherodon melanotheron is unusual for its ability to colonise estuarine environments in West Africa that are characterised by extremely high salinities. The relationships between mRNA levels of the 70 kDa heat shock protein (Hsp70) and Na⁺, K⁺-ATPase1α (Naka) in the gills, environmental salinity, and a life-history trait (condition factor) were investigated in wild populations of this species sampled from three locations in the Saloum estuary, at salinities ranging from 40 to 100 psu.

Results

The highest Hsp70 and Naka mRNA levels, and the poorest condition factors were recorded in the most saline sampling site (100 psu). The Hsp70 and Naka mRNA were correlated amongst themselves and showed a direct positive correlation with environmental salinity, and a negative correlation with fish condition factor. Thus, the Hsp70 is constitutively overexpressed by S. melanotheron acclimatised to extreme hypersalinity.

Conclusions

These results indicate that, although S. melanotheron can colonise extremely saline environments, the overexpression of Hsp70 combined with the higher Naka mRNA expression reveals that this represents a chronic stress. The induction of Hsp70 was, therefore, a biomarker of chronic hyper-osmotic stress which presumably can be linked to the impaired growth performance and precocious reproduction that have been demonstrated in the populations at the extremely saline sites.

Screening for coping style increases the power of gene expression studies

Background

Individuals of many vertebrate species show different stress coping styles and these have a striking influence on how gene expression shifts in response to a variety of challenges.

Principal Findings

This is clearly illustrated by a study in which common carp displaying behavioural predictors of different coping styles (characterised by a proactive, adrenaline-based or a reactive, cortisol-based response) were subjected to inflammatory challenge and specific gene transcripts measured in individual brains. Proactive and reactive fish differed in baseline gene expression and also showed diametrically opposite responses to the challenge for 80% of the genes investigated.

Significance

Incorporating coping style as an explanatory variable can account for some the unexplained variation that is common in gene expression studies, can uncover important effects that would otherwise have passed unnoticed and greatly enhances the interpretive value of gene expression data.