Local adaptation and oceanographic connectivity patterns explain genetic differentiation of a marine diatom across the North Sea–Baltic Sea salinity gradient

Drivers of population genetic structure are still poorly understood in marine micro‐organisms. We exploited the North Sea–Baltic Sea transition for investigating the seascape genetics of a marine diatom, Skeletonema marinoi. Eight polymorphic microsatellite loci were analysed in 354 individuals from ten locations to analyse population structure of the species along a 1500‐km‐long salinity gradient ranging from 3 to 30 psu. To test for salinity adaptation, salinity reaction norms were determined for sets of strains originating from three different salinity regimes of the gradient. Modelled oceanographic connectivity was compared to directional relative migration by correlation analyses to examine oceanographic drivers. Population genetic analyses showed distinct genetic divergence of a low‐salinity Baltic Sea population and a high‐salinity North Sea population, coinciding with the most evident physical dispersal barrier in the area, the Danish Straits. Baltic Sea populations displayed reduced genetic diversity compared to North Sea populations. Growth optima of low salinity isolates were significantly lower than those of strains from higher native salinities, indicating local salinity adaptation. Although the North Sea–Baltic Sea transition was identified as a barrier to gene flow, migration between Baltic Sea and North Sea populations occurred. However, the presence of differentiated neutral markers on each side of the transition zone suggests that migrants are maladapted. It is concluded that local salinity adaptation, supported by oceanographic connectivity patterns creating an asymmetric migration pattern between the Baltic Sea and the North Sea, determines genetic differentiation patterns in the transition zone.     

A fundamental difference between macrobiota and microbial eukaryotes: protistan plankton has a species maximum in the freshwater-marine transition zone of the Baltic Sea

Remane's Artenminimum at the horohalinicum is a fundamental concept in ecology to describe and explain the distribution of organisms along salinity gradients. However, a recent metadata analysis challenged this concept for protists, proposing a species maximum in brackish waters. Due to data bias, this literature-based investigation was highly discussed. Reliable data verifying or rejecting the species minimum for protists in brackish waters were critically lacking. Here, we sampled a pronounced salinity gradient along a west–east transect in the Baltic Sea and analysed protistan plankton communities using high-throughput eDNA metabarcoding. A strong salinity barrier at the upper limit of the horohalinicum and 10 psu appeared to select for significant shifts in protistan community structures, with dinoflagellates being dominant at lower salinities, and dictyochophytes and diatoms being keyplayers at higher salinities. Also in vertical water column gradients in deeper basins (Kiel Bight, Arkona and Bornholm Basin) appeared salinity as significant environmental determinant influencing alpha- and beta-diversity patterns. Importantly, alpha-diversity indices revealed species maxima in brackish waters, that is, indeed contrasting Remane's Artenminimum concept. Statistical analyses confirmed salinity as the major driving force for protistan community structuring with high significance. This suggests that macrobiota and microbial eukaryotes follow fundamentally different rules regarding diversity patterns in the transition zone from freshwater to marine waters.     

Scenario simulations of future salinity and ecological consequences in the Baltic Sea and adjacent North Sea areas–implications for environmental monitoring

Substantial ecological changes occurred in the 1970s in the Northern Baltic during a temporary period of low salinity (S). This period was preceded by an episodic increase in the rainfall over the Baltic Sea watershed area. Several climate models, both global and regional, project an increase in the runoff of the Northern latitudes due to proceeding climate change. The aim of this study is to model, firstly, the effects on Baltic Sea salinity of increased runoff due to projected global change and, secondly, the effects of salinity change on the distribution of marine species. The results suggest a critical shift in the S range 5–7, which is a threshold for both freshwater and marine species distributions and diversity. We discuss several topics emphasizing future monitoring, modelling, and fisheries research. Environmental monitoring and modelling are investigated because the developing alternative ecosystems do not necessarily show the same relations to environment quality factors as the retiring ones. An important corollary is that the observed and modelled S changes considered together with species’ ranges indicate what may appear under a future climate. Consequences could include a shift in distribution areas of marine benthic foundation species and some 40–50 other species, affiliated to these. This change would extend over hundreds of kilometres, in the Baltic Sea and the adjacent North Sea areas. Potential cascading effects, in coastal ecology, fish ecology and fisheries would be extensive, and point out the necessity to develop further the “ecosystem approach in the environmental monitoring”.     

A comparative analysis of metabarcoding and morphology‐based identification of benthic communities across different regional seas

In a world of declining biodiversity, monitoring is becoming crucial. Molecular methods, such as metabarcoding, have the potential to rapidly expand our knowledge of biodiversity, supporting assessment, management, and conservation. In the marine environment, where hard substrata are more difficult to access than soft bottoms for quantitative ecological studies, Artificial Substrate Units (ASUs) allow for standardized sampling. We deployed ASUs within five regional seas (Baltic Sea, Northeast Atlantic Ocean, Mediterranean Sea, Black Sea, and Red Sea) for 12–26 months to measure the diversity and community composition of macroinvertebrates. We identified invertebrates using a traditional approach based on morphological characters, and by metabarcoding of the mitochondrial cytochrome oxidase I (COI) gene. We compared community composition and diversity metrics obtained using the two methods. Diversity was significantly correlated between data types. Metabarcoding of ASUs allowed for robust comparisons of community composition and diversity, but not all groups were successfully sequenced. All locations were significantly different in taxonomic composition as measured with both kinds of data. We recovered previously known regional biogeographical patterns in both datasets (e.g., low species diversity in the Black and Baltic Seas, affinity between the Bay of Biscay and the Mediterranean). We conclude that the two approaches provide complementary information and that metabarcoding shows great promise for marine monitoring. However, until its pitfalls are addressed, the use of metabarcoding in monitoring of rocky benthic assemblages should be used in addition to classical approaches rather than instead of them.     

Genetic diversity of Saccharina latissima (Phaeophyceae) along a salinity gradient in the North Sea–Baltic Sea transition zone

The North Sea–Baltic Sea transition zone constitutes a boundary area for the kelp species Saccharina latissima due to a strong salinity gradient operating in the area. Furthermore, the existence of S. latissima there, along Danish waters, is fairly patchy as hard bottom is scarce. In this study, patterns of genetic diversity of S. latissima populations were evaluated along the salinity gradient area of Danish waters (here designated brackish) and were compared to reference sites (here designated marine) outside the gradient area, using microsatellite markers. The results showed that the S. latissima populations were structured into two clusters corresponding to brackish versus marine sites, and that gene flow was reduced both between clusters and between populations within clusters. In addition, results provided empirical evidence that marginal populations of S. latissima in the salinity gradient area exhibited a distinct genetic structure when compared to marine ones. Brackish populations were less diverse, more related, and showed increased differentiation over distance compared to marine populations. The isolation of the brackish S. latissima populations within the salinity gradient area of Danish waters in conjunction with their general low genetic diversity makes these populations vulnerable to ongoing environmental and climate change, predicted to result in declining salinity in the Baltic Sea area that may alter the future distribution and performance of S. latissima in the area.     

Diversity and abundance of “Pelagibacterales” (SAR11) in the Baltic Sea salinity gradient

The candidate order “Pelagibacterales” (SAR11) is one of the most abundant bacterial orders in ocean surface waters and, periodically, in freshwater lakes. The presence of several stable phylogenetic lineages comprising “Pelagibacterales” correlates with the physico-chemical parameters in aquatic environments. A previous amplicon sequencing study covering the bacterial community in the salinity gradient of the Baltic Sea suggested that pelagibacteral subclade SAR11-I was replaced by SAR11-IIIa in the mesohaline region of the Baltic Sea. In this current study, we investigated the cellular abundances of “Pelagibacterales” subclades along the Baltic Sea salinity gradient using catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH). The results obtained with a newly designed probe, which exclusively detected SAR11-IIIa, were compared to CARD-FISH abundances of the marine SAR11-I/II subclade and the freshwater lineage SAR11-IIIb (LD12). The results showed that SAR11-IIIa was abundant in oligohaline–mesohaline conditions (salinities 2.7–13.3), with maximal abundances at a salinity of 7 (up to 35% of total Bacteria, quantified with a universal bacterial probe EUB). As expected, SAR11-I/II was abundant (27% of EUB) in the marine parts of the Baltic Sea, whereas counts of the freshwater lineage SAR11-IIIb were below the detection limit at all stations. The shift from SAR11-IIIa to SAR11-I/II was confirmed in the vertical salinity gradient in the deeper basins of the Baltic Sea. These findings were consistent with an overlapping but defined distribution of SAR11-I/II and SAR11-IIIa in the salinity gradient of the Baltic Sea and suggested the adaptation of SAR11-IIIa for growth and survival in mesohaline conditions.     

Zonation of intertidal macrobenthos in the estuaries of Schelde and Ems

Based on data, collected in 1980–1990, the intertidal benthic macrofauna of the Schelde and Ems estuaries was compared. The spatial occurrence of the benthic macrofauna along the salinity gradient, including the freshwater tidal area was emphasized. Both estuaries appeared to have a very similar species composition, especially at genus level. The higher number of species observed in the Schelde estuary was probably due to a greater habitat diversity. In both estuaries species diversity decreased with distance upstream. The total density did not vary along the estuarine gradient, whereas biomass is highest in the polyhaline zone.In both estuaries distinct intertidal benthic communities were observed along the salinity gradient: a marine community in the polyhaline zone, a brackish community in the mesohaline zone, and a third community in the oligohaline and freshwater tidal zones of the estuary. These three communities were very similar between both estuaries. Their main characteristics were discussed together with the occurrence and distribution of the dominant species.For the Schelde estuary and to a lesser extent also for the Ems estuary, there was evidence that anthropogenic stress had a negative effect on the intertidal macrobenthic communities of the oligohaline/freshwater tidal zone. Only Oligochaeta were dominating, whereas the very euryhaline and/or true limnetic species were missing. In the mesohaline zone, the Schelde estuary was dominated by large numbers of short-living, opportunistic species, whereas in the Ems estuary relatively more stable macrobenthic communities were observed. A comparison with some other European estuaries showed in general similar trends as those observed for the Schelde and Ems estuaries.     

Salinity drives meiofaunal community structure dynamics across the Baltic ecosystem

Coastal benthic biodiversity is under increased pressure from climate change, eutrophication, hypoxia, and changes in salinity due to increase in river runoff. The Baltic Sea is a large brackish system characterized by steep environmental gradients that experiences all of the mentioned stressors. As such it provides an ideal model system for studying the impact of on‐going and future climate change on biodiversity and function of benthic ecosystems. Meiofauna (animals < 1 mm) are abundant in sediment and are still largely unexplored even though they are known to regulate organic matter degradation and nutrient cycling. In this study, benthic meiofaunal community structure was analysed along a salinity gradient in the Baltic Sea proper using high‐throughput sequencing. Our results demonstrate that areas with higher salinity have a higher biodiversity, and salinity is probably the main driver influencing meiofauna diversity and community composition. Furthermore, in the more diverse and saline environments a larger amount of nematode genera classified as predators prevailed, and meiofauna‐macrofauna associations were more prominent. These findings show that in the Baltic Sea, a decrease in salinity resulting from accelerated climate change will probably lead to decreased benthic biodiversity, and cause profound changes in benthic communities, with potential consequences for ecosystem stability, functions and services.     

The benthic macrofauna along the estuarine gradient of the Schelde estuary

The intertidal benthic macrofauna of the Schelde estuary (The Netherlands and Belgium) was sampled in late autumn of 1990 at 50 stations along the whole salinity gradient (between Vlissingen and Dendermonde), including the freshwater tidal part. All stations were situated in sheltered areas with a relatively muddy sediment. Species richness, diversity and total biomass of the benthic macrofauna decreased along the salinity gradient from Vlissingen to Dendermonde, while total density showed no clear trend. Especially the oligohaline and freshwater tidal part of the Schelde estuary was characterized by a very impoverished benthic community, composed only of Oligochaeta. No other species (freshwater, marine or brackish) was observed in this part of the estuary. The marine part had a more diverse macrozoobenthos structure than that of the brackish part. Species found only in the marine zone areCerastoderma edule, Tharyx marioni, Eteone longa, Nephtys hombergii andCapitella capitata. In the brackish part of the estuary,Corophium volutator was a typical, dominant species. However, a lot of the dominant species were common in both the marine and brackish part of the Schelde estuary (e.g. Heteromastus filiformis, Pygospio elegans, Nereis diversicolor, Macoma balthica). The observed gradient in species composition and dominance is compared with some other European estuaries. The marine and brackish part of the Schelde estuary is quite similar to other european estuaries. The freshwater tidal part, however, was more impoverished.     

Cryptic freshwater amoeba species in the bottom sediments of Nivå Bay (Øresund, Baltic Sea)

An aquarium containing a portion of freshly collected brackish-water sediment from Nivå Bay (Øresund, Baltic Sea) was sampled to determine the species of amoebae present and refilled with inorganic salt medium to achieve freshwater salinity. After 2 months incubation, the diversity of amoebae was again studied, and the salinity was restored to the original level. The aquarium was incubated for a further 2 months and the amoeba fauna was studied for the last time. A number of freshwater species appeared in the sample after the first salinity shift, while some marine species disappeared. Most marine species did not re-appear after the salinity in the aquarium was restored, but one marine species not noted previously was recorded. The experiment illustrates the presence of ‘cryptic diversity’ of amoebae in natural habitats and demonstrates that laboratory manipulation of the salinity of a sample prior to inoculation may achieve an increased recovery of species from a brackish-water habitat.     

Nearshore subtidal community structure compared between inner coast and outer coast sites in Southeast Alaska

Processes that structure subarctic marine communities, particularly in glaciated regions, are not well understood. This understanding is needed as a baseline and to manage these communities in the face of future climate-driven changes. This study investigates two coastal regions of Southeast Alaska with the goals to (a) identify and compare patterns of subtidal community structure for macroalgal, fish, macroinvertebrate (>5?cm), and small epibenthic invertebrate (<5?cm) communities between inner coast and outer coast sites and (b) link patterns of community structure to habitat and environmental parameters. Species assemblage and benthic habitat data were used to compare species diversity and community composition at 6?m and 12?m depths at nine inner coast and nine outer coast sites. Multivariate analysis was applied to reduce environmental variables to major gradients, to resolve community structure, and to relate community structure to habitat and environmental variables. Increased salinity and decreased temperature at outer coast sites compared with inner coast sites were associated with community structure, with greater species diversity at outer coast sites at 6?m depth. Invertebrate community composition was associated with benthic habitat, including crust and coralline algae for macroinvertebrates, and algal cover and substrate for small epibenthic invertebrates. This research suggests that marine communities in glaciated regions are strongly influenced by freshwater input and that future climate-driven changes in freshwater input will likely result in marine community composition changes.     

Phylogeographical structure, distribution and genetic variation of the green algae Ulva intestinalis and U. compressa (Chlorophyta) in the Baltic Sea area

The marine algae Ulva intestinalis and U. compressa are morphologically plastic with many overlapping characters and are therefore difficult to distinguish from each other. The present distribution of U. intestinalis and U. compressa is investigated along the salinity gradient in the Baltic Sea area through analyses of internal transcribed spacer (ITS) sequence data. Also, the amount and distribution of intraspecific genetic polymorphism in the ITS region is studied allowing inferences on the phylogeographical pattern and postglacial recolonization of the Baltic Sea area. The data show that of the two species only U. intestinalis occurs in the Baltic Sea. The distribution of U. compressa is more restricted than previously reported, and it was not found in salinities lower than 15 ppt. All of Scandinavia and the Baltic Sea were covered with ice during the last ice age and the organisms in the Baltic Sea must have colonized the area after the ice had started to melt. The genetic diversity of U. intestinalis and U. compressa in the Baltic Sea and the neighbouring area was found to be reduced compared to that in the British Isles. This reduction may be the result of either a historical reduction of diversity or an adaptation of specific clones to the northern environmental conditions.     

COASTAL DIATOM‐ENVIRONMENT RELATIONSHIPS IN THE BRACKISH BALTIC SEA1

High‐quality calibration data sets are required when diatom assemblages are used for monitoring ecological change or reconstructing palaeo‐environments. The quality of such data sets can be validated, in addition to other criteria, by the percentage of significant unimodal species responses as a measure of the length of an environmental gradient. This study presents diatom‐environment relationships analyzed from a robust data set of diatom communities living on submerged stones along a 2,000 km long coastline in the Baltic Sea area, including 524 samples taken at 135 sites and covering a salinity gradient from 0.4 to 11.4. Altogether, 487 diatom taxa belonging to 102 genera were recorded. Detrended canonical correspondence analysis showed that salinity was the overriding environmental factor regulating diatom community composition, while exposure to wave action and nutrient concentrations were of secondary importance. Modeling the abundances of the 58 most common diatom taxa yielded significant relationships with salinity for 57 taxa. Twenty‐three taxa showing monotonic responses were species with optimum distributions in freshwater or marine waters. Thirty‐four taxa showing unimodal responses were brackish‐water species with maximum distributions at different salinities. Separate analyses for small (cell biovolume <1,000 μm³) and large (≥1,000 μm³) taxa yielded similar results. In previous studies along shorter salinity gradients, large and small epilithic diatom taxa responded differently. From our large data, we conclude that counts of large diatom taxa alone seem sufficient for indicating salinity changes in coastal environments with high precision.     

Microphytobenthic Primary Production and Sedimentary Carbohydrates Along Salinity Gradients in the Lagoons of Grado and Marano (Northern Adriatic Sea)

Primary production of the microphytobenthic community and carbohydrates concentrations were studied in the lagoonal system of Grado and Marano, located in the Northern Adriatic coast. Sediment samples were collected along a salinity gradient. Abundance and species composition of the microphytobenthic communities were analysed and the benthic microalgal biomass was estimated as Chlorophyll a (Chl a). Primary production of benthic diatoms was estimated using ¹⁴C-tracer. Extracellular carbohydrates were extracted from the sediment and separated in two operationally defined fractions (colloidal and EDTA-extractable). Salinity was higher in the Grado lagoon, where the benthic microalgal community was mainly composed of marine diatoms. In the Marano lagoon, which has a lower salinity, freshwater species were also found. In both lagoons, photosynthetic efficiency showed an inverse relationship with salinity and a direct relationship with the main biological variables. Photosynthetic activity was directly related to Chl a and abundance of benthic microalgae, suggesting that in the benthic system microalgal community is responsible for primary production. Overall, salinity was also influent on the microphytobenthic primary production, which was greater in the more saline Grado lagoon.     

Effect of salinity on cable bacteria species composition and diversity

Cable bacteria (CB) are Desulfobulbaceae that couple sulphide oxidation to oxygen reduction over centimetre distances by mediating electric currents. Recently, it was suggested that the CB clade is composed of two genera, Ca. Electronema and Ca. Electrothrix, with distinct freshwater and marine habitats respectively. However, only a few studies have reported CB from freshwater sediment, making this distinction uncertain. Here, we report novel data to show that salinity is a controlling factor for the diversity and the species composition within CB populations. CB sampled from a freshwater site (salinity 0.3) grouped into Ca. Electronema and could not grow under brackish conditions (salinity 21), whereas CB from a brackish site (salinity 21) grouped into Ca. Electrothrix and decreased by 93% in activity under freshwater conditions. On a regional scale (Baltic Sea), salinity significantly influenced species richness and composition. However, other environmental factors, such as temperature and quantity and quality of organic matter were also important to explain the observed variation. A global survey of 16S rRNA gene amplicon sequencing revealed that the two genera did not co-occur likely because of competitive exclusion and identified a possible third genus.     

Genomic population structure of freshwater‐resident and anadromous ide (Leuciscus idus) in north‐western Europe

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large‐scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping‐by‐Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north‐western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.     

Rating species sensitivity throughout gradient systems – a consistent approach for the Baltic Sea

Evaluating the state of benthic communities has played an important role in water quality assessments. Indices incorporating species sensitivities, richness and densities are commonly applied. In Europe, the importance of benthic indices has increased in the last years with the implementation of the European Marine Strategy Framework Directive (MSFD) which at the same time demands the applicability of an index across regional scales. To date, environmental variability is rarely considered in benthic indices and most sensitivity rankings have the disadvantages of static values (i.e. the same value in all areas), expert judgement and a limited geographical range.This study presents species sensitivity values calculated along environmental gradients for the Baltic Sea. Sensitivities were calculated according to the procedure of the Benthic Quality Index (BQI). We created a matrix of subregions, classes of salinity, depth and gear to identify comparable subsets for data analysis. Altogether, 19 subsets were defined within the Baltic Sea basins. Sensitivity values were calculated for 329 species out of a total of 678 species that were recorded in this study. Sensitivity values of taxa vary between subsets as it was expected for different environmental conditions. Most sensitivity values can be assigned to species occurring in euhaline and polyhaline waters. Distribution of species with high and low sensitivity values differed along the salinity gradient. In euhaline waters more species with high sensitivity values occurred than species with low sensitivity values, while in mesohaline waters the ratio of high and low sensitivity values among species was almost equal. In oligohaline waters more species with lower sensitivity values were present.For the first time, sensitivity values were calculated for a large number of species using the same method for the entire Baltic Sea. This results in a Baltic-wide comprehensive set of sensitivity values based on a dataset across subregional borders, and divided along environmental gradients and gear type. The same principles can be applied to transient waters from rivers to coastal lagoons as well as to other environments with gradients of, e.g. hydrodynamic characteristics. Publicly available sensitivity values will increase transparency and support the improvement of state assessments under the MSFD.     

Imperfect detection distorts depth‐related trends in marine macrofaunal species richness

Ecology and biodiversity research are underpinned by species richness patterns and their environmental drivers. However, a key topic in this discussion is the accuracy of these patterns which are greatly dependent on species detection probabilities. Due to variations in detection of species, true ecological patterns may be distorted. This is particularly true for subtidal macro‐infaunal communities. We tested three hypothesized relationships between marine benthic macrofaunal diversity and depth using species richness per site estimated with a capture–recapture heterogeneity model that accounts for variable detection probabilities. These metrics were based on data from 42 replicated sites across the continental shelf of the Southern Benguela. Average detection probability decreased with greater depth but species richness increased along the same depth gradient. The conflation of these trends in observed diversity data resulted in a positively near‐linear depth–diversity relationship, while accounting for variable species detection revealed a much stronger relationship. Ignoring species detection in ecosystems with imperfect detection could therefore distort species richness patterns, which has implications for ecological theory, management and conservation.