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

Diverse Transcriptome Responses to Salinity Change in Atlantic Cod Subpopulations

Adaptation to environmental variation caused by global climate change is a significant aspect of fisheries management and ecology. A reduction in ocean salinity is visible in near-shore areas, especially in the Baltic Sea, where it is affecting the Atlantic cod population. Cod is one of the most significant teleost species, with high ecological and economical value worldwide. The population of cod in the Baltic Sea has been traditionally divided into two subpopulations (western and eastern) existing in higher- and lower-salinity waters, respectively. In recent decades, both Baltic cod subpopulations have declined massively. One of the reasons for the poor condition of cod in the Baltic Sea is environmental factors, including salinity. Thus, in this study, an oligonucleotide microarray was applied to explore differences between Baltic cod subpopulations in response to salinity fluctuations. For this purpose, an exposure experiment was conducted consisting of salinity elevation and reduction, and gene expression was measured in gill tissue. We found 400 differentially expressed genes (DEGs) involved in the immune response, metabolism, programmed cell death, cytoskeleton, and extracellular matrix that showed a subpopulation-dependent pattern. These findings indicate that osmoregulation in Baltic cod is a complex process, and that western and eastern Baltic cod subpopulations respond differently to salinity changes.

Whole Transcriptome Analysis Reveals the Global Molecular Responses of mRNAs, lncRNAs, miRNAs, circRNAs, and Their ceRNA Networks to Salinity Stress in Hong Kong Oysters, Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is an estuarine bivalve with remarkable commercial value in South China, and the increase of salinity in estuaries during the dry season has posed a major threat to the oyster farming. To explore the global transcriptional response to salinity stress, a whole-transcriptome analysis was performed with the gills of oysters in 6‰, 18‰, and 30‰ filtered seawater. Overall, 2243, 194, 371, and 167 differentially expressed mRNAs (DEmRNAs), differentially expressed long non-coding RNAs (DElncRNAs), differentially expressed circular RNAs (DEcircRNAs), and differentially expressed microRNAs (DEmiRNAs) were identified, respectively. Based on GO enrichment and KEGG pathway analysis, these important DEmRNAs, DElncRNAs, DEcircRNAs, and DEmiRNAs were predicted to be mainly involved in amino acids metabolism, microtubule movement, and immune defense. This demonstrated the complexity of dynamic transcriptomic profiles of C. hongkongensis in response to salinity fluctuation. The regulatory relationships of DEmiRNAs-DEmRNAs, DElncRNAs-DEmiRNAs, and DEcircRNAs-DEmiRNAs were also predicted, and finally, a circRNA-associated competing endogenous RNA (ceRNA) network was constructed, consisting of six DEcircRNAs, eight DEmiRNAs, and five DEmRNAs. The key roles of taurine and hypotaurine metabolism and phenylalanine metabolism were highlighted in this ceRNA network, which was consistent with the major contribution of free amino acids to intracellular osmolality and cell volume regulation. Collectively, this study provides comprehensive data, contributing to the exploration of coding and non-coding RNAs in C. hongkongensis salinity response. The results would benefit the understanding of the response mechanism of bivalves against salinity fluctuation, and provide clues for genetic improvement of C. hongkongensis with hyper-salinity tolerance.

Novel lncRNA lncRNA001074 participates in the low salinity-induced response in the sea cucumber Apostichopus japonicus by targeting the let-7/NKAα axis

Salinity fluctuations have severe impacts on sea cucumbers and therefore important consequences in sea cucumber farming. The responses of sea cucumbers to salinity changes are reflected in the expression profiles of multiple genes and non-coding RNAs (ncRNAs). The microRNA (let-7) which is a developmental regulator, the ion transporter gene sodium potassium ATPase gene (NKAα), and the long ncRNA lncRNA001074 were previously shown to be involved in responses to salinity changes in various marine species. To better understand the relationship between ncRNAs and target genes, the let-7/NKAα/lncRNA001074 predicted interaction was investigated in this study using luciferase reporter assays and gene knockdowns in the sea cucumber Apostichopus japonicus. The results showed that NKAα was the target gene of let-7 and NKAα expression levels were inversely correlated with let-7 expression based on the luciferase reporter assays and western blots. The let-7 abundance was negatively regulated by lncRNA001074 and NKAα both in vitro and in vivo. Knockdown of lncRNA001074 led to let-7 overexpression. These results demonstrated that lncRNA001074 binds to the 3'-UTR binding site of let-7 in a regulatory manner. Furthermore, the expression profiles of let-7, NKAα, and lncRNA001074 were analyzed in sea cucumbers after the knockdown of each of these genes. The results found that lncRNA001074 competitively bound let-7 to suppress NKAα expression under low salinity conditions. The downregulation of let-7, in conjunction with the upregulation of lncRNA001074 and NKAα, may be essential for the response to low salinity change in sea cucumbers. Therefore, the dynamic balance of the lncRNA001074, NKAα, and let-7 network might be a potential response mechanism to salinity change in sea cucumbers.

Disparate movement behavior and feeding ecology in sympatric ecotypes of Atlantic cod

Coexistence of ecotypes, genetically divergent population units, is a widespread phenomenon, potentially affecting ecosystem functioning and local food web stability. In coastal Skagerrak, Atlantic cod (Gadus morhua) occur as two such coexisting ecotypes. We applied a combination of acoustic telemetry, genotyping, and stable isotope analysis to 72 individuals to investigate movement ecology and food niche of putative local "Fjord" and putative oceanic "North Sea" ecotypes-thus named based on previous molecular studies. Genotyping and individual origin assignment suggested 41 individuals were Fjord and 31 were North Sea ecotypes. Both ecotypes were found throughout the fjord. Seven percent of Fjord ecotype individuals left the study system during the study while 42% of North Sea individuals left, potentially homing to natal spawning grounds. Home range sizes were similar for the two ecotypes but highly variable among individuals. Fjord ecotype cod had significantly higher δ13C and δ15N stable isotope values than North Sea ecotype cod, suggesting they exploited different food niches. The results suggest coexisting ecotypes may possess innate differences in feeding and movement ecologies and may thus fill different functional roles in marine ecosystems. This highlights the importance of conserving interconnected populations to ensure stable ecosystem functioning and food web structures.

Effects of Sudden Drop in Salinity on Osmotic Pressure Regulation and Antioxidant Defense Mechanism of Scapharca subcrenata

Salinity is an important ecological factor that impacts the growth and survival of aquatic organisms. The salinity of seawater in coastal and estuarine areas is often subject to dynamic changes because of seasonal rainfall and continental runoff. Thus, the current study investigated the effects of sudden changes in salinity on the survival rate and osmotic pressure regulation mechanisms of bottom-sowing seedlings of the economically important ark shell, Scapharca subcrenata. By simulating the sudden changes that occur in seawater salinity after rainstorms, the results showed that the osmotic pressure of the hemolymph and Na⁺, K⁺, Ca²⁺, and Cl^– concentrations first decreased and then increased. When the salinity decreased from 30 to 14‰, hemoglobin, soluble total protein, taurine, and total free amino acid gradually increased; maximum levels of hemoglobin, soluble total protein, and taurine occurred once the salinity increased to 22‰ at 96 h. After 96 h, the total free amino acid content increased until 144 h. The reactive oxygen species (ROS) content and total antioxidant capacity (T-AOC) peaked at 96 h, whereas the expression levels of Mn-superoxide dismutase (MnSOD) and catalase (CAT) increased earlier, indicating that, with continuous ROS generation, antioxidant defense mechanisms were activated to avoid oxidative damage. Expression levels of cathepsin C (CTSC), cathepsin D (CTSD), heat shock protein 20 (HSP20), and heat shock protein 70 (HSP70) were significantly higher than in the control group at 48 h (salinity level 14‰); the expression levels of HSP20, heat shock protein 90 (HSP90), MnSOD, and glutathione peroxidase (GP_x) remained high, indicating that they were still required for osmotic pressure regulation to maintain the dynamic balance between the generation and removal of ROS as the salinity level increased. These results not only add to our basic understanding of the aquatic ecology of S. subcrenata, but also provide a theoretical ground for improving the survival rate of bottom-sowing, propagation, and release of S. subcrenata seedlings.

SNP genotyping reveals substructuring in weakly differentiated populations of Atlantic cod (Gadus morhua) from diverse environments in the Baltic Sea

Atlantic cod (Gadus morhua) is one of the most important fish species in northern Europe for several reasons including its predator status in marine ecosystems, its historical role in fisheries, its potential in aquaculture and its strong public profile. However, due to over-exploitation in the North Atlantic and changes in the ecosystem, many cod populations have been reduced in size and genetic diversity. Cod populations in the Baltic Proper, Kattegat and North Sea have been analyzed using a species specific single nucleotide polymorphism (SNP) array. Using a subset of 8,706 SNPs, moderate genetic differences were found between subdivisions in three traditionally delineated cod management stocks: Kattegat, western and eastern Baltic. However, an FST measure of population differentiation based on allele frequencies from 588 outlier loci for 2 population groups, one including 5 western and the other 4 eastern Baltic populations, indicated high genetic differentiation. In this paper, differentiation has been demonstrated not only between, but also within western and eastern Baltic cod stocks for the first time, with salinity appearing to be the most important environmental factor influencing the maintenance of cod population divergence between the western and eastern Baltic Sea.

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 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.

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.

Transcriptional response of yellow perch to changes in ambient metal concentrations-A reciprocal field transplantation experiment

Recent local adaptation to pollution has been evidenced in several organisms inhabiting environments heavily contaminated by metals. Nevertheless, the molecular mechanisms underlying adaptation to high metal concentrations are poorly understood, especially in fishes. Yellow perch (Perca flavescens) populations from lakes in the mining area of Rouyn-Noranda (QC, Canada) have been faced with metal contamination for about 90 years. Here, we examine gene transcription patterns of fish reciprocally transplanted between a reference and a metal-contaminated lake and also fish caged in their native lake. After four weeks, 111 genes were differentially transcribed in metal-naïve fish transferred to the metal-contaminated lake, revealing a plastic response to metal exposure. Genes involved in the citric cycle and beta-oxidation pathways were under-transcribed, suggesting a potential strategy to mitigate the effects of metal stress by reducing energy turnover. However, metal-contaminated fish transplanted to the reference lake did not show any transcriptomic response, indicating a reduced plastic response capability to sudden reduction in metal concentrations. Moreover, the transcription of other genes, especially ones involved in energy metabolism, was affected by caging. Overall, our results highlight environmental stress response mechanisms in yellow perch at the transcriptomic level and support a rapid adaptive response to metal exposure through genetic assimilation.

Adaptation to Low Salinity Promotes Genomic Divergence in Atlantic Cod (Gadus morhua L.)

How genomic selection enables species to adapt to divergent environments is a fundamental question in ecology and evolution. We investigated the genomic signatures of local adaptation in Atlantic cod (Gadus morhua L.) along a natural salinity gradient, ranging from 35‰ in the North Sea to 7‰ within the Baltic Sea. By utilizing a 12 K SNPchip, we simultaneously assessed neutral and adaptive genetic divergence across the Atlantic cod genome. Combining outlier analyses with a landscape genomic approach, we identified a set of directionally selected loci that are strongly correlated with habitat differences in salinity, oxygen, and temperature. Our results show that discrete regions within the Atlantic cod genome are subject to directional selection and associated with adaptation to the local environmental conditions in the Baltic- and the North Sea, indicating divergence hitchhiking and the presence of genomic islands of divergence. We report a suite of outlier single nucleotide polymorphisms within or closely located to genes associated with osmoregulation, as well as genes known to play important roles in the hydration and development of oocytes. These genes are likely to have key functions within a general osmoregulatory framework and are important for the survival of eggs and larvae, contributing to the buildup of reproductive isolation between the low-salinity adapted Baltic cod and the adjacent cod populations. Hence, our data suggest that adaptive responses to the environmental conditions in the Baltic Sea may contribute to a strong and effective reproductive barrier, and that Baltic cod can be viewed as an example of ongoing speciation.

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

Salinity is one of the key factors that affects metabolism, survival and distribution of fish species, as all fish osmoregulate and euryhaline fish maintain osmotic differences between their extracellular fluid and either freshwater or seawater. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations in both marine and freshwater environments, where the physiological and genomic basis for salinity tolerance adaptation is not fully understood. Therefore, our main objective in this study was to investigate gene expression of three targeted osmoregulatory genes (Na⁺/K⁺-ATPase (ATPA13), cystic fibrosis transmembrane regulator (CFTR) and a voltage gated potassium channel gene (KCNH4) and one stress related heat shock protein gene (HSP70)) in gill tissue from marine and freshwater populations when exposed to non-native salinity for periods ranging from five minutes to three weeks. Overall, the targeted genes showed highly plastic expression profiles, in addition the expression of ATP1A3 was slightly higher in saltwater adapted fish and KCNH4 and HSP70 had slightly higher expression in freshwater. As no pronounced changes were observed in the expression profiles of the targeted genes, this indicates that the osmoregulatory apparatuses of both the marine and landlocked freshwater stickleback population have not been environmentally canalized, but are able to respond plastically to abrupt salinity challenges.

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene-linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high-resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.

Application of microsatellite markers in conservation genetics and fisheries management: recent advances in population structure analysis and conservation strategies

Microsatellites are the most popular and versatile genetic marker with myriads of applications in population genetics, conservation biology, and evolutionary biology. These are the arrays of DNA sequences, consisting of tandemly repeating mono-, di-, tri-, and tetranucleotide units, which are distributed throughout the genomes of most eukaryotic species. Microsatellites are codominant in nature, highly polymorphic, easily typed, and Mendelian inherited, all properties which make them very suitable for the study of population structure and pedigree analysis and capable of detecting differences among closely related species. PCR for microsatellites can be automated for identifying simple sequence repeat polymorphism. Small amount of blood samples or alcohol preserved tissue is adequate for analyzing them. Most of the microsatellites are noncoding, and therefore variations are independent of natural selection. These properties make microsatellites ideal genetic markers for conservation genetics and fisheries management. This review addresses the applications of microsatellite markers in conservation genetics and recent advances in population structure analysis in the context of fisheries management.

Temporal genetic structure in a poecilogonous polychaete: the interplay of developmental mode and environmental stochasticity

Background

Temporal variation in the genetic structure of populations can be caused by multiple factors, including natural selection, stochastic environmental variation, migration, or genetic drift. In benthic marine species, the developmental mode of larvae may indicate a possibility for temporal genetic variation: species with dispersive planktonic larvae are expected to be more likely to show temporal genetic variation than species with benthic or brooded non-dispersive larvae, due to differences in larval mortality and dispersal ability. We examined temporal genetic structure in populations of Pygospio elegans, a poecilogonous polychaete with within-species variation in developmental mode. P. elegans produces either planktonic, benthic, or intermediate larvae, varying both among and within populations, providing a within-species test of the generality of a relationship between temporal genetic variation and larval developmental mode.

Results

In contrast to our expectations, our microsatellite analyses of P. elegans revealed temporal genetic stability in the UK population with planktonic larvae, whereas there was variation indicative of drift in temporal samples of the populations from the Baltic Sea, which have predominantly benthic and intermediate larvae. We also detected temporal variation in relatedness within these populations. A large temporal shift in genetic structure was detected in a population from the Netherlands, having multiple developmental modes. This shift could have been caused by local extiction due to extreme environmental conditions and (re)colonization by planktonic larvae from neighboring populations.

Conclusions

In our study of P. elegans, temporal genetic variation appears to be due to not only larval developmental mode, but also the stochastic environment of adults. Large temporal genetic shifts may be more likely in marine intertidal habitats (e.g. North Sea and Wadden Sea) which are more prone to environmental stochasticity than the sub-tidal Baltic habitats. Sub-tidal and/or brackish (less saline) habitats may support smaller P. elegans populations and these may be more susceptible to the effects of random genetic drift. Moreover, higher frequencies of asexual reproduction and the benthic larval developmental mode in these populations leads to higher relatedness and contributes to drift. Our results indicate that a general relationship between larval developmental mode and temporal genetic variation may not exist.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene-based genome-scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection-as determined with several outlier detection methods-was low (F(ST) = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (F(ST) = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Population genetics provides new insights into biomarker prevalence in dab (Limanda limanda L.): a key marine biomonitoring species

Bioindicators are species for which some quantifiable aspect of its biology, a biomarker, is assumed to be sensitive to ecosystem health. However, there is frequently a lack of information on the underlying genetic and environmental drivers shaping the spatiotemporal variance in prevalence of the biomarkers employed. Here, we explore the relative role of potential variables influencing the spatiotemporal prevalence of biomarkers in dab, Limanda limanda, a species used as a bioindicator of marine contaminants. Firstly, the spatiotemporal genetic structure of dab around UK waters (39 samples across 15 sites for four years: 2005–2008) is evaluated with 16 microsatellites. Two temporally stable groups are identified corresponding to the North and Irish Seas (average between basin = 0.007; = 0.022). Secondly, we examine the association between biomarker prevalence and several variables, including genetic structuring, age and contaminant exposure. Genetic structure had significant interactive effects, together with age and some contaminants, in the prevalence of some of the biomarkers considered, namely hyperpigmentation and liver lesions. The integration of these data sets enhanced our understanding of the relationship between biomarker prevalence, exposure to contaminants and population-specific response, thereby yielding more informative predictive models of response and prospects for environmental remediation.

Heterogeneous Genomic Differentiation in marine threespine sticklebacks: adaptation along an environmental gradient

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on-going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F(ST) ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F(ST) ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers--45% of which were identified as outliers--correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.

Oceanographic connectivity and environmental correlates of genetic structuring in Atlantic herring in the Baltic Sea

Marine fish often show little genetic structuring in neutral marker genes, and Atlantic herring (Clupea harengus) in the Baltic Sea are no exception; historically, very low levels of population differentiation (F ST ≍ 0.002) have been found, despite a high degree of interpopulation environmental heterogeneity in salinity and temperature. Recent exome sequencing and SNP studies have however shown that many loci are under selection in this system. Here, we combined population genetic analyses of a large number of transcriptome-derived microsatellite markers with oceanographic modelling to investigate genetic differentiation and connectivity in Atlantic herring at a relatively fine scale within the Baltic Sea. We found evidence for weak but robust and significant genetic structuring (F ST = 0.008) explainable by oceanographic connectivity. Genetic differentiation was also associated with site differences in temperature and salinity, with the result driven by the locus Her14 which appears to be under directional selection (F ST = 0.08). The results show that Baltic herring are genetically structured within the Baltic Sea, and highlight the role of oceanography and environmental factors in explaining this structuring. The results also have implications for the management of herring fisheries, the most economically important fishery in the Baltic Sea, suggesting that the current fisheries management units may be in need of revision.

Acclimation to different depths by the marine angiosperm Posidonia oceanica: transcriptomic and proteomic profiles

For seagrasses, seasonal and daily variations in light and temperature represent the mains factors driving their distribution along the bathymetric cline. Changes in these environmental factors, due to climatic and anthropogenic effects, can compromise their survival. In a framework of conservation and restoration, it becomes crucial to improve our knowledge about the physiological plasticity of seagrass species along environmental gradients. Here, we aimed to identify differences in transcriptomic and proteomic profiles, involved in the acclimation along the depth gradient in the seagrass Posidonia oceanica, and to improve the available molecular resources in this species, which is an important requisite for the application of eco-genomic approaches. To do that, from plant growing in shallow (-5 m) and deep (-25 m) portions of a single meadow, (i) we generated two reciprocal Expressed Sequences Tags (EST) libraries using a Suppressive Subtractive Hybridization (SSH) approach, to obtain depth/specific transcriptional profiles, and (ii) we identified proteins differentially expressed, using the highly innovative USIS mass spectrometry methodology, coupled with 1D-SDS electrophoresis and labeling free approach. Mass spectra were searched in the open source Global Proteome Machine (GPM) engine against plant databases and with the X!Tandem algorithm against a local database. Transcriptional analysis showed both quantitative and qualitative differences between depths. EST libraries had only the 3% of transcripts in common. A total of 315 peptides belonging to 64 proteins were identified by mass spectrometry. ATP synthase subunits were among the most abundant proteins in both conditions. Both approaches identified genes and proteins in pathways related to energy metabolism, transport and genetic information processing, that appear to be the most involved in depth acclimation in P. oceanica. Their putative rules in acclimation to depth were discussed.