Genetic differentiation of brackish water populations of cod Gadus morhua in the southern Baltic, inferred from genotyping using SNP-arrays

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

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.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s)

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Population Genetic Study on the European Flounder (Platichthys flesus) from the Southern Baltic Sea Using SNPs and Microsatellite Markers

The European flounder (Platichthys flesus), which is closely related to the recently discovered Baltic flounder (Platichthys solemdali), is currently the third most commercially fished species in the Baltic Sea. According to the available data from the Polish Fisheries Monitoring Center and fishermen's observations, the body condition indices of the species in the Baltic Sea have declined in recent years. The aim of the present study was to obtain information on the current patterns of genetic variability and the population structure of the European flounder and to verify whether the Baltic flounder is present in the southern Baltic Sea. Moreover, we aimed to verify whether the observed decline in the body condition indices of the species in the Baltic Sea might be associated with adaptive alterations in its gene pool due to increased fishing pressure. For this purpose, 190 fish were collected from four locations along the central coastline of Poland, i.e., Mechelinki, Władysławowo, the Vistula Lagoon in 2018, and the Słupsk Bank in 2020. The fish were morphologically analyzed and then genetically screened by the application of nineteen microsatellite DNA and two diagnostic SNP markers. The examined European flounder specimens displayed a high level of genetic diversity (PIC = 0.832-0.903, I = 2.579-2.768). A lack of significant genetic differentiation (Fst = 0.004, p > 0.05) was observed in all the examined fish, indicating that the European flounder in the sampled area constitutes a single genetic cluster. A significant deficiency in heterozygotes (Fis = 0.093, p < 0.05) and overall deviations from Hardy-Weinberg expectations (H-WE) were only detected in fish sampled from the Słupsk Bank. The estimated effective population size (Ne) among the sampled fish groups varied from 712 (Słupsk Bank) to 10,115 (Władysławowo and Mechelinki). However, the recorded values of the Garza-Williamson indicator (M = 0.574-0.600) and the lack of significant (p > 0.05) differences in Heq > He under the SMM model did not support the species' population size changes in the past. The applied SNP markers did not detect the presence of the Baltic flounder among the fish sampled from the studied area. The analysis of an association between biological traits and patterns of genetic diversity did not detect any signs of directional selection or density-dependent adaptive changes in the gene pool of the examined fish that might be caused by increased fishing pressure.

Single Nucleotide Polymorphism Markers with Applications in Conservation and Exploitation of Aquatic Natural Populations

Simple Summary

In recent decades, societies, states and local authorities have become increasingly aware that for effective long-term management and protection of aquatic ecosystems and populations, it is necessary to take into account the genetic changes occurring in these populations. One type of high-resolution molecular marker suitable for studying the neutral and adaptive genetic diversity of populations is single nucleotide polymorphism (SNP). This review is an attempt to show the benefits of using SNPs to recognize natural populations of aquatic animals and detect the threats to them from accidentally or intentionally released farm animals, fishery and global climate changes. It is postulated that conservation actions should protect not only pristine natural populations that are endangered or overfished, but also protect populations of non-threatened species from unnecessarily released semi-domesticated animals. The enhancement of natural populations with farmed material usually reduces their genetic diversity. Experimental size-selective catches of artificially created populations have caused evolutionary changes in the life cycles of fishes. However, fishery-induced evolution in natural populations is difficult to observe. The negative measurable effects on populations can be expected when the number of breeding individuals is reduced below 100, which occurs very rarely in the sea and more often in fragmented freshwater streams, ponds and seasonal rivers.

Abstract

An increasing number of aquatic species have been studied for genetic polymorphism, which extends the knowledge on their natural populations. One type of high-resolution molecular marker suitable for studying the genetic diversity of large numbers of individuals is single nucleotide polymorphism (SNP). This review is an attempt to show the range of applications of SNPs in studies of natural populations of aquatic animals. In recent years, SNPs have been used in the genetic analysis of wild and enhanced fish and invertebrate populations in natural habitats, exploited migratory species in the oceans, migratory anadromous and freshwater fish and demersal species. SNPs have been used for the identification of species and their hybrids in natural environments, to study the genetic consequences of restocking for conservation purposes and the negative effects on natural populations of fish accidentally escaping from culture. SNPs are very useful for identifying genomic regions correlated with phenotypic variants relevant for wildlife protection, management and aquaculture. Experimental size-selective catches of populations created in tanks have caused evolutionary changes in life cycles of fishes. The research results have been discussed to clarify whether the fish populations in natural conditions can undergo changes due to selective harvesting targeting the fastest-growing fishes.

Identification of Fish Species and Targeted Genetic Modifications Based on DNA Analysis: State of the Art

Food adulteration is one of the most serious problems regarding food safety and quality worldwide. Besides misleading consumers, it poses a considerable health risk associated with the potential non-labeled allergen content. Fish and fish products are one of the most expensive and widely traded commodities, which predisposes them to being adulterated. Among all fraud types, replacing high-quality or rare fish with a less valuable species predominates. Because fish differ in their allergen content, specifically the main one, parvalbumin, their replacement can endanger consumers. This underlines the need for reliable, robust control systems for fish species identification. Various methods may be used for the aforementioned purpose. DNA-based methods are favored due to the characteristics of the target molecule, DNA, which is heat resistant, and the fact that through its sequencing, several other traits, including the recognition of genetic modifications, can be determined. Thus, they are considered to be powerful tools for identifying cases of food fraud. In this review, the major DNA-based methods applicable for fish meat and product authentication and their commercial applications are discussed, the possibilities of detecting genetic modifications in fish are evaluated, and future trends are highlighted, emphasizing the need for comprehensive and regularly updated online database resources.

The use of various statistical methods for authenticity and detection of adulteration in fish and seafood

Various methodologies including genetic analyses, morphometrics, proteomics, lipidomics, metabolomics, etc. are now used or being developed to authenticate fish and seafood. Such techniques usually lead to the generation of enormous amounts of data. The analysis and interpretation of this information can be particularly challenging. Statistical techniques are therefore commonly used to assist in analyzing these data, visualizing trends and differences and extracting conclusions. This review article aims at presenting and discussing statistical methods used in studies on fish and seafood authenticity and adulteration, allowing researchers to consider their options based on previous successes/failures but also offering some recommendations about the future of such techniques. Techniques such as PCA, AMOVA and FST statistics, that allow the differentiation of genetic groups, or techniques such as MANOVA that allow large data sets of morphometric characteristics or elemental differences to be analyzed are discussed. Furthermore, methods such as cluster analysis, DFA, CVA, CDA and heatmaps/Circos plots that allow samples to be differentiated based on their geographical origin are also reviewed and their advantages and disadvantages as found in past studies are given. Finally, mathematical simulations and modeling are presented in a detailed review of studies using them, together with their advantages and limitations.

The use of molecular markers in the verification of fish and seafood authenticity and the detection of adulteration

The verification of authenticity and detection of food mislabeling are elements that have been of high importance for centuries. During the last few decades there has been an increasing consumer demand for the verification of food identity and the implementation of stricter controls around these matters. Fish and seafood are among the most easily adulterated foodstuffs mainly due to the significant alterations of the species' morphological characteristics that occur during the different types of processing, which render the visual identification of the animals impossible. Even simple processes, such as filleting remove very important morphological elements and suffice to prevent the visual identification of species in marketed products. Novel techniques have therefore been developed that allow identification of species, the differentiation between species and also the differentiation of individuals that belong to the same species but grow in different populations and regions. Molecular markers have been used during the last few decades to fulfill this purpose and several improvements have been implemented rendering their use applicable to a commercial scale. The reliability, accuracy, reproducibility, and time-and cost-effectiveness of these techniques allowed them to be established as routine methods in the industry and research institutes. This review article aims at presenting the most important molecular markers used for the authentication of fish and seafood. The most important techniques are described, and the results of numerous studies are outlined and discussed, allowing interested parties to easily access and compare information about several techniques and fish/seafood species.

Climate change driven hyposalinity as a selective agent in the littoral mesoherbivore Idotea balthica

Climate change will include a decrease in seawater salinity in the Baltic Sea. We quantified the effects of the projected future desalination on survival of the early life stage of the littoral herbivore Idotea balthica. We collected egg-bearing Idotea from three range-margin Baltic Sea populations, we exposed half of each brood to either current (6‰) or future salinity (3.5‰). We genotyped a subsample of each brood to analyse patterns of allelic change and to identify genomic regions targeted by selection. The survival was overall reduced by hyposalinity and broods varied in response to hyposalinity implying genetic variation in tolerance, with a stronger decrease in genetic diversity in future salinity. Finally, we identified proteins with crucial roles in basic cellular functions. This study indicates that projected future northern Baltic Sea hyposalinity will not just hamper I. balthica survival, but its selective pressure may also affect genetic diversity and cell physiology.

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.

Assessing SNP-markers to study population mixing and ecological adaptation in Baltic cod

Atlantic cod (Gadus morhua) is a species of great ecological and economical importance in the Baltic Sea. Here, two genetically differentiated stocks, the western and the eastern Baltic cod, display substantial mechanical mixing, hampering our understanding of cod ecology and impeding stock assessments and management. Based on whole-genome re-sequencing data from reference samples obtained from the study area, we designed two different panels of Single Nucleotide Polymorphisms markers (SNPs), which take into account the exceptional genome architecture of cod. A minimum panel of 20 diagnostic SNPs and an extended panel (20 diagnostic and 18 biologically informative SNPs, 38 in total) were developed and validated to distinguish unambiguously between the western and the eastern Baltic cod stocks and to enable studies of local adaptation to the specific environment in the Baltic Sea, respectively. We tested both panels on cod sampled from the southern Baltic Sea (n = 603) caught in 2015 and 2016. Genotyping results showed that catches from the mixing zone in the Arkona Sea, were composed of similar proportions of individuals of the western and the eastern stock. Catches from adjacent areas to the east, the Bornholm Basin and Gdańsk Deep, were exclusively composed of eastern Baltic cod, whereas catches from adjacent western areas (Belt Sea and Öresund) were composed of western Baltic cod. Interestingly, the two Baltic cod stocks showed strong genetic differences at loci associated with life-history trait candidate genes, highlighting the species’ potential for ecological adaptation even at small geographical scales. The minimum and the extended panel of SNP markers presented in this study provide powerful tools for future applications in research and fisheries management to further illuminate the mixing dynamics of cod in the Baltic Sea and to better understand Baltic cod ecology.

Genetic analyses reveal complex dynamics within a marine fish management area

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Alternatively spliced variants in Atlantic cod (Gadus morhua) support response to variable salinity environment

Analysis of gill transcriptome of the Atlantic cod from the Baltic Sea demonstrated that alternatively spliced (AS) variants may be actively involved in the process of adaptation to altered salinity. Some AS variants of different genes, like phospholipase A2 group IVC (PLA2G4C), appeared only in fish exposed to altered salinity, while other isoforms of the same genes were present in all experimental groups. Novel sequence arrangements represent 89% of all AS in the Baltic cod compared to the Atlantic population. Profiles of modified pathways suggest that regulation by AS can afford specific changes of genes expressed in response to the environment. The AS variants appear to be involved in the response to stress by modifications of signalling in apoptosis pathways, an innate immunological response and pro-inflammatory process. Present results support the hypothesis that developing new AS variants could support genome complexity and reinforce the ability to fast adapt to local environments.

Phylogeographic mitogenomics of Atlantic cod Gadus morhua: Variation in and among trans-Atlantic, trans-Laurentian, Northern cod, and landlocked fjord populations

The historical phylogeography, biogeography, and ecology of Atlantic cod (Gadus morhua) have been impacted by cyclic Pleistocene glaciations, where drops in sea temperatures led to sequestering of water in ice sheets, emergence of continental shelves, and changes to ocean currents. High-resolution, whole-genome mitogenomic phylogeography can help to elucidate this history. We identified eight major haplogroups among 153 fish from 14 populations by Bayesian, parsimony, and distance methods, including one that extends the species coalescent back to ca. 330 kya. Fish from the Barents and Baltic Seas tend to occur in basal haplogroups versus more recent distribution of fish in the Northwest Atlantic. There was significant differentiation in the majority of trans-Atlantic comparisons (ΦST = .029-.180), but little or none in pairwise comparisons within the Northwest Atlantic of individual populations (ΦST = .000-.060) or defined management stocks (ΦST = .000-.023). Monte Carlo randomization tests of population phylogeography showed significantly nonrandom trans-Atlantic phylogeography versus absence of such structure within various partitions of trans-Laurentian, Northern cod (NAFO 2J3KL) and other management stocks, and Flemish Cap populations. A landlocked meromictic fjord on Baffin Island comprised multiple identical or near-identical mitogenomes in two major polyphyletic clades, and was significantly differentiated from all other populations (ΦST = .153-.340). The phylogeography supports a hypothesis of an eastern origin of genetic diversity ca. 200-250 kya, rapid expansion of a western superhaplogroup comprising four haplogroups ca. 150 kya, and recent postglacial founder populations.

Historic changes in length distributions of three Baltic cod (Gadus morhua) stocks: Evidence of growth retardation

Understanding how combinations of fishing effort and selectivity affect productivity is central to fisheries research. We investigate the roles of fishing regulation in comparison with ecosystem status for Baltic Sea cod stock productivity, growth performance, and population stability. This case study is interesting because three cod populations with different exploitation patterns and stock status are located in three adjacent but partially, ecologically different areas. In assessing stock status, growth, and productivity, we use survey information and rather basic stock parameters without relying on age readings. Because there is an urgent interest of better understanding of the current development of the Eastern Baltic cod stock, we argue that our approach represents partly a novel way of interpreting monitoring information together with catch data in a simplified yet more informative way. Our study reports how the Eastern and Western Baltic cod have gone toward more truncated size structures between 1991 and 2016, in particular for the Eastern Baltic cod, whereas the Öresund cod show no trend. We suggest that selective fishing may disrupt fish population dynamic stability and that lower natural productivity might amplify the effects of selective fishing. In support of earlier findings on a density‐dependent growth of Eastern Baltic cod, management is advised to acknowledge that sustainable exploitation levels for Eastern Baltic cod are much more limited than perceived in regular assessments. Of more general importance, our results emphasize the need to embrace a more realistic view on what ecosystems can produce regarding tractable fish biomass to facilitate a more ecosystem‐based fisheries management.

Selection by higher-order effects of salinity and bacteria on early life-stages of Western Baltic spring-spawning herring

Habitat stratification by abiotic and biotic factors initiates divergence of populations and leads to ecological speciation. In contrast to fully marine waters, the Baltic Sea is stratified by a salinity gradient that strongly affects fish physiology, distribution, diversity and virulence of important marine pathogens. Animals thus face the challenge to simultaneously adapt to the concurrent salinity and cope with the selection imposed by the changing pathogenic virulence. Western Baltic spring‐spawning herring (Clupea harengus) migrate to spawning grounds characterized by different salinities to which herring are supposedly adapted. We hypothesized that herring populations do not only have to cope with different salinity levels but that they are simultaneously exposed to higher‐order effects that accompany the shifts in salinity, that is induced pathogenicity of Vibrio bacteria in lower saline waters. To experimentally evaluate this, adults of two populations were caught in their spawning grounds and fully reciprocally crossed within and between populations. Larvae were reared at three salinity levels, representing the spawning ground salinity of each of the two populations, or Atlantic salinity conditions resembling the phylogenetic origin of Clupea harengus. In addition, larvae were exposed to a Vibrio spp. infection. Life‐history traits and gene expression analysis served as response variables. Herring seem adapted to Baltic Sea conditions and cope better with low saline waters. However, upon a bacterial infection, herring larvae suffer more when kept at lower salinities implying reduced resistance against Vibrio or higher Vibrio virulence. In the context of recent climate change with less saline marine waters in the Baltic Sea, such interactions may constitute key future stressors.

Construction of a high-density genetic map and QTL mapping for pearl quality-related traits in Hyriopsis cumingii

A high-density genetic map is essential for quantitative trait locus (QTL) fine mapping. In this study, 4,508 effective single nucleotide polymorphism markers (detected using specific-locus amplified fragment sequencing) and 475 microsatellites were mapped to 19 linkage groups (LGs) using a family with 157 individuals. The map spanned 2,713 cM, with an average of 259 markers and 79 loci per LG and an average inter-marker distance of 1.81 cM. To identify QTLs for pearl quality traits, 26 putatively significant QTLs were detected for 10 traits, including, three for shell width, seven for body weight, two for shell weight, two for margin mantle weight, five for inner mantle weight, and seven for shell nacre colour. Among them, five QTLs associated with shell nacre colour were mapped to LG17 and explained 19.7% to 22.8% of the trait variation; this suggests that some important genes or loci determine shell nacre colour in LG17. The linkage map and mapped QTLs for shell nacre colour would be useful for improving the quality of Hyriopsis cumingii via marker-assisted selection.

Transcriptome analysis of gill tissue of Atlantic cod Gadus morhua L. from the Baltic Sea

The Atlantic cod (Gadus morhua L.) is one of the most ecologically and economically important marine fish species in the North Atlantic Ocean. Using Roche GS-FLX 454 pyrosequencing technique 962,516 reads, representing 379Mbp of the Baltic cod transcriptome, were obtained. Data was assembled into 14,029 contigs of which 100% displayed homology to the Atlantic cod transcriptome. Despite a high similarity between transcripts, evidence for significant differences between Baltic and Atlantic cod was found.