Skip to main content

Importance of plasticity and local adaptation for coping with changing salinity in coastal areas: a test case with barnacles in the Baltic Sea

Abstract

Background

Salinity plays an important role in shaping coastal marine communities. Near-future climate predictions indicate that salinity will decrease in many shallow coastal areas due to increased precipitation; however, few studies have addressed this issue. The ability of ecosystems to cope with future changes will depend on species’ capacities to acclimatise or adapt to new environmental conditions. Here, we investigated the effects of a strong salinity gradient (the Baltic Sea system – Baltic, Kattegat, Skagerrak) on plasticity and adaptations in the euryhaline barnacle Balanus improvisus. We used a common-garden approach, where multiple batches of newly settled barnacles from each of three different geographical areas along the Skagerrak-Baltic salinity gradient were exposed to corresponding native salinities (6, 15 and 30 PSU), and phenotypic traits including mortality, growth, shell strength, condition index and reproductive maturity were recorded.

Results

We found that B. improvisus was highly euryhaline, but had highest growth and reproductive maturity at intermediate salinities. We also found that low salinity had negative effects on other fitness-related traits including initial growth and shell strength, although mortality was also lowest in low salinity. Overall, differences between populations in most measured traits were weak, indicating little local adaptation to salinity. Nonetheless, we observed some population-specific responses – notably that populations from high salinity grew stronger shells in their native salinity compared to the other populations, possibly indicating adaptation to differences in local predation pressure.

Conclusions

Our study shows that B. improvisus is an example of a true brackish-water species, and that plastic responses are more likely than evolutionary tracking in coping with future changes in coastal salinity.

Background

Future climate-driven changes in the marine environment are projected to include decreased salinity in many coastal areas due to increased precipitation and enhanced freshwater run-off [13]. Most climate change research to date has focused on effects of increasing temperatures and ocean acidification, rather than salinity [46]. Salinity plays an important role in shaping the distribution of marine species and future alterations in salinity may pose major ecological challenges to organisms inhabiting coastal areas [7]. Changes in salinity are known to have direct or indirect effects on survival, metabolism, growth, reproduction, and/or osmotic balance in aquatic organisms [812]. Future salinity shifts in coastal areas may therefore impose strong selection on species inhabiting these areas, leading to changes in species composition and the evolution of new adaptations. Understanding how marine organisms respond to environmental changes and how rapidly new adaptations can evolve is key for predicting how ecosystems will respond to global environmental change in the future.

In the Baltic Sea, salinity forms a strong environmental gradient that determines the composition and distribution of species [13, 14]. It is one of the world’s largest semi-enclosed brackish seas and reached its present brackish state about 8000 years ago [15, 16]. Since then, the Baltic Sea has been colonised by organisms from both freshwater and marine environments [16, 17]. The sea surface salinity along this gradient ranges from < 3 PSU (practical salinity units) in the northern Bothnian Bay (Figure 1) to approximately 30 PSU at the border to the North Sea [13]. Limited water exchange with surrounding seas and almost no tidal flow make salinity conditions in the Baltic Sea relatively stable compared to many other coastal areas [18, 19]. Species diversity and within-species genetic diversity are markedly lower inside the brackish Baltic Sea compared to adjacent areas, which to a large extent is attributed to these marginal environmental conditions creating specific selection pressures and promoting the evolution of local adaptations [2022]. Recent modelling of future climate scenarios in the Baltic Sea indicates that increased precipitation will lead to reduced salinity [3, 23], which could potentially result in dramatic shifts in species’ distributions. Such shifts would also likely influence ecological interactions and ecosystem functions.

Figure 1
figure 1

Map of the study area. The Skagerrak-Baltic Sea region is characterized by a strong salinity gradient. The sites where the barnacles were collected (shown as stars) include: Stockholm (59.286 E; 18.708 N), Kiel fjord (54.327 E; 10.185 N) and Tjärnö (58.881E; 11.134 N).

Organisms can respond to changes in salinity through plastic responses involving molecular (e.g. gene expression), physiological, morphological, behavioural or life history trait changes, without genotypic changes occurring through selection [24]. Here, we mainly consider phenotype plasticity in the context of “phenotypic buffering”, i.e. a plasticity that maintains a functional phenotype despite external disturbances [25]. From an evolutionary perspective, it has been suggested that phenotypic plasticity can act as a buffer against selection that may delay or even preclude the evolution of local adaptations. However, it can also be considered a source of novel opportunities upon which selection can act, e.g. by facilitating the transition to new fitness levels while maintaining a large population size [2628]. In either case, phenotypic plasticity is likely to play an important role in the ability of species to tolerate future environmental changes [29]. Being highly plastic is often assumed to involve a cost, which would explain why plasticity is not universal [30, 31]. Consequently, it has been suggested that in more stable environments, where (costly) high plasticity is not required, adaptations to a narrower range of environmental conditions may be selected for, resulting in a loss of phenotypic plasticity [3234]. Adaptation, in contrast to phenotypic plasticity, always involves a genetic change, resulting from selection acting on fitness-related phenotypic traits. Although there is mounting evidence that local adaptation may be common in the marine environment, there are few studies on the relative role of evolutionary responses versus plasticity to changing salinity in marine invertebrates [35], but see [36].

Barnacles (Crustacea; Cirripeda) generally display high tolerance to fluctuating environments. They make up an important component of intertidal and estuarine communities, where they are frequently exposed to variation in salinity, as well as fluctuating temperature and desiccation [37]. The bay barnacle Balanus (Amphibalanus) improvisus[38] is the only barnacle species recorded from the Baltic Sea, but it is also found in other coastal areas and estuaries worldwide [39]. B. improvisus is believed to originate from the eastern North America and was first recorded in the Baltic Sea in the mid 1800’s [40]. Compared to most other barnacle species, B. improvisus can tolerate a very wide range of salinities, from <1 to 35 PSU, however, salinities below 3 PSU have been reported to have long-term negative impacts on fitness [4143]. The molecular and cellular mechanisms behind this broad salinity tolerance are poorly known, although it has been proposed that B. improvisus can alternate between being a passive osmoconformer in high salinities (>15 PSU) and actively osmoregulate in salinities below 15 PSU [41]. It was recently shown that alternatively spliced variants of Na+/K+ ATPases were differentially expressed in response to reduced salinity giving further support for plasticity through active osmoregulation [44].

If there is indeed a cost to the euryhalinity (maintaining high physiological plasticity in relation to salinity) of B. improvisus, then this capacity may be reduced, or even lost, in populations living in a stable and low salinity environment such as the Baltic Sea. Previous studies comparing the effects of salinity in different barnacle species found that B. improvisus had a broad tolerance and fed actively in most salinity conditions, and that the response depended on previous salinity history [45]. More specifically, B. improvisus from the Baltic Sea showed highest activity in low salinities (6 PSU), compared to populations of B. improvisus from England (higher salinities), potentially indicating local adaptation to the Baltic environment [45]. This may indicate that rapid evolution had occurred in Baltic Sea populations of B. improvisus over a timescale of a few hundred years after colonisation. However, the study by Davenport [45] only tested short-term exposures to different salinities without appropriate acclimation of barnacles collected from the field, which limits the ability to draw conclusions about the degree of local adaptation. By performing long-term controlled experiments, starting from early post-settlement and following barnacles until sexual maturity, we can better understand how salinity affects a wide range of phenotypic traits and reveal which traits may be more exposed to selection in different salinity conditions, potentially resulting in local adaptations.

Here, we investigated whether the broad distribution of B. improvisus along the strong salinity gradient of the Baltic Sea can be explained by physiological plasticity or if local adaptations to various salinity regimes have evolved in populations along the gradient. We used a common-garden experiment in which multiple barnacle populations sampled along the salinity gradient were grown from immediate post-settlement to reproductive maturity under different salinity treatments. We measured several fitness-related traits including growth, mortality, shell strength and reproductive maturity.

Methods

Sampling of brood stock barnacles and larval culture

Several hundred adult barnacles (Balanus (Amphibalanus) improvisus) were collected on settling panels from three different populations in: i) the Stockholm archipelago, Baltic Sea (Sweden, 4-6 PSU [19]; 59.286 E; 18.708 N), ii) the Kiel Fjord, Baltic Sea (Germany, 14-17 PSU [46]; 54.327 E; 10.185 N) and iii) Tjärnö, Swedish west coast (Sweden, 22-30 PSU [47]; 58.881E; 11.134 N) during July to August 2011 (Figure 1). Settled barnacles were transported to the laboratory at Tjärnö (Sven Lovén Centre for Marine Sciences) and kept in large re-circulating systems at their native mean salinities (6, 15 and 30 PSU) and at a constant temperature of 19°C. Adult barnacles were fed daily with newly hatched Artemia sp. ad libitum. As in other barnacle species, B. improvisus produces free-swimming larvae (six nauplius stages followed by a non-feeding cyprid stage) with pelagic dispersal for up to several weeks before they settle on hard substrates and metamorphose into the adult stage [48]. Nauplius larvae were obtained through natural spawning of the adults in a modified laboratory culture system described previously [49]. Nauplii were cultured (at maximum 500 nauplii L-1) until the cyprid stage in their native salinities, which were created by mixing filtered seawater (0.2 μm) with filtered tap water (0.2 μm) and fed with a mixture of the microalgae Thalassiosira pseudonana and Skeletonema marinoi (at a ratio of 40:60, respectively, and in ad libitum). These microalgae were also cultured at different salinities (10 PSU, 15 PSU and 30 PSU) in order to reduce osmotic disturbances when transferring algae to the larval cultures. We took into account the possible natural variability between larval batches by collecting larvae from the adult barnacle cultures (consisting of several hundred individuals each) on four different occasions during two weeks in September-October 2011, representing four replicate batches of larvae (presumably with different parents, and thus different genotypes). Larval batches were used as level of replication for each population, since natural variability in larval performance between batches has been observed previously [50]. Cyprid larvae from each batch were settled on acrylic plastic (Plexiglas) panels and were transferred to the experimental system (see below) after eight days, leaving enough time for the larvae to attach and metamorphose before the start of the experiment.

Experimental setup

The experiment involved three salinity treatments (6, 15 and 30 PSU), corresponding to the native mean salinities for the three barnacle populations studied (Stockholm, Kiel and Tjärnö, respectively). From each of the four batches obtained per population, two panels with newly settled barnacles (eight days old) were placed in separate aquaria in each of the three salinity treatments, resulting in a split-plot design [51]. Three large independent re-circulating systems (370 l each) were used, each consisting of 24 aquaria (6 l each). The flow through each aquarium was 20-25 l/h. Each system had a biological filter (Fluval filter 405) with coral rubble, activated charcoal and nitrifying bacteria added to maintain water quality and avoid build-up of nitrogenous waste products. In addition, water in each system was completely replaced every two weeks, at which time the culture systems were cleaned, drained and treatment combinations were moved to a new system within the same room to avoid confounding the effects of salinity with effects of aquarium location within the room. The experiment was maintained at 20°C and a light regime of 14:10 h (L:D). The different salinities were obtained by mixing filtered deep saltwater (0.2 μm) from the Kosterfjord (30-34 PSU, total alkalinity, AT, of 2186-2290 μmol l-1) with filtered tap water (AT of 369 μmol l-1). Water quality variables including temperature, salinity, AT, NH3 and NH4 were monitored routinely during the experiment (Additional file 1: Table S1). Salinity was adjusted by adding freshwater whenever needed. Barnacles were fed a mixture of microalgae (Skeletonema marinoi and Chaetoceros calcitrans) at ~ 20,000 cells ml-1 and 30,000 cells ml-1, respectively. Algal cell concentrations were checked regularly using a Multisizer™ 3 Coulter Counter (Beckman Coulter) and levels were adjusted to maintain stable concentrations ad libitum throughout the experiment. The algal cell concentration in the systems never reached below 5,000 cells ml-1. Chlorophyll content of the algal cultures and experimental systems were also checked for possible degradation of food quality, assessed by a spectrophotometric trichromatic method [52], which showed no signs of degradation. After four weeks, the barnacle diet was complemented with newly hatched Artemia (ca. 3,000 Artemia per aquarium) added every second day. Barnacles were cultured in the experimental systems for a total of nine weeks.

Growth modelling and condition index

Barnacle growth was recorded photographically every two weeks, (Olympus E5 DSLR, 50 mm F1.8 Macro lens). To avoid overcrowding, excess barnacles were removed from panels during the first two weeks, and subsequently whenever barnacles came in contact with each other. At the end of the experiment 5 - 96 barnacles remained on each panel (mean 43 ± 15 (SD)). The different barnacle densities were evenly distributed between populations and salinity treatments. Digital images were converted to binary files using ImageJ (version 1.43). An image analysis script in MATLAB (R2012a) was developed to track and measure maximum basal plate diameter (rostro-carinal) for each barnacle at four time points (after 2, 4, 6, and 9 weeks). Only barnacles that remained alive at the end of the experiment were included in the analysis. Two different growth models were fitted to the size data for each individual barnacle. The first model was the von Bertalanffy growth equation [53]:

Y t = S max 1 - b e - kt
(1)

where Y(t) is the shell diameter at time t, S max is the asymptote of the curve, b is the lag phase (initial growth), and k is the growth rate. The second model was the Gompertz growth equation [54]:

Y t =a e b e ct
(2)

where Y(t) is the shell diameter at time t, a is the maximum shell diameter (asymptote), b is the lag phase or early growth phase and c is the maximum growth rate (how quickly the individual reaches the asymptote; note b and c are negative numbers) (Figure 2). The von Bertalanffy growth model is widely used, especially in fisheries studies, but it has also earlier been fitted to barnacle growth data [55]. The Gompertz growth model is commonly used to describe growth in fish [56], but has also been used in studies on marine invertebrates [57]. The von Bertalanffy model assumes that growth rate declines over time, whereas the Gompertz model includes assumptions of slower growth in the beginning and at the end of the growth phase (sigmoid curve). It has been suggested that a sigmoid curve (such as the Gompertz model) is more applicable to data for larval and early juvenile stages [58]. Since barnacle larvae go through a metamorphosis to become adults, including the formation of a calcareous shell, this could potentially result in reduced initial growth rate. We fitted both growth models to our data and compared them using the coefficient of determination, R2, and root mean square error (rmse). We found no statistically significant difference in fit between the two models (t = 0.23, df = 90, P = 0.82) based on rmse, although the Gompertz growth equation provided a better fit than von Bertalanffy in a majority of cases (47/68). We therefore used Gompertz curves in all further growth analysis. Estimates of the three Gompertz model parameters were obtained for each individual barnacle (Figure 2) and the mean value for all barnacles on each replicate panel was used in subsequent statistical analysis. Dry weight (DW) of barnacles from each batch and treatment combination was determined for 20 haphazardly selected individuals from each panel. Samples were dried (80°C, 48 h) and weighed to the nearest 0.0001 g (Sartorius CP analytical balance). Dried samples were then burned (20 h, 500°C) and the remaining inorganic material was weighed (ash weight, AW), and ash-free dry weight (AFDW) calculated. We calculated a condition index (CI) defined as the body-mass to shell-mass ratio (AFDW/AW). This measure has previously been used to evaluate physiological stress under environmental changes in invertebrates, e.g. [59, 60]. Results are presented as norms of reaction, a common way to describe the relation between the genetic background (different populations) and the phenotypic responses (growth parameters) across a range of environments, e.g. different salinities [61].

Figure 2
figure 2

Illustration of the Gompertz growth model. The Gompertz growth model is a sigmoid function described by three parameters; a (asymptote), b (lag phase) and c (maximal growth rate). Model parameters for each replicate panel with barnacles, were obtained by first fitting a model (thin lines) to shell diameter data (circles) for each barnacle on a panel (measurements at four different times) and then calculating an average model from these fittings for each replicate panel (thick line).

Shell strength

To investigate whether salinity affected shell strength, the strength of 20 haphazardly chosen barnacles on each panel was determined using a TAXT2i Texture Analyzer (Stable Micro Systems, 25-1 measuring cell). The pressure (compressive force) required to break the barnacle shell from above (“breaking stress”) was determined according to [59]. Maximum “breaking stress” was measured using a cylinder of 2 mm diameter pushing onto the rostral plate with a speed of 1.0 mm s-1. Measured shell strength (in MPa) was further normalised by dividing by the maximum shell diameter of each individual barnacle.

Reproductive maturity

Twenty barnacles from each batch were inspected under a dissection microscope to determine their state of reproductive maturity. Gonad maturity and reproductive status could easily be determined since B. improvisus broods its offspring for up to several weeks before releasing them into the water column. The proportion of individuals with well-developed gonads (ovaries and/or testes – B. improvisus is an hermaphrodite) or fertilized eggs was recorded at the end of the experiment.

Mortality

Mortality was estimated by comparing the number of living barnacles on each panel after two weeks with the number alive after nine weeks. Mortality estimates for the first two weeks of the experiment could not be obtained due to missing data (missing photographs), although no major mortality events were observed during this period. Mortality estimates were corrected for removal (culling) of barnacles during the experiment.

Statistical analyses

General and generalized linear models (“aov” and “glm” respectively in R) were fitted using R-Studio (version 096.331) and SPSS (version 21). Effects of salinity and population origin on growth, shell strength, DW, AW, AFDW, CI and mortality were tested using a 2-factor ANOVA, with ‘salinity’ and ‘population’ as fixed factors. Batch was used as the level of replication in all analyses. For all the response variables we evaluated the main effects of salinity, population, and their interaction (an indicator of local adaptation). Significant interactions between salinity and population would indicate that populations respond differently to different salinity environments, which could be a result of local adaptations to a specific salinity regime. Assumptions of normality and homogeneity of variances were checked using Q/Q-plots, box plots and Levene’s test. AW was log-transformed and per cent mortality was arcsin-transformed prior to analysis. Tukey’s HSD test was used post-hoc to further resolve significant differences between means. Due to high frequencies of zero count data, fecundity (measured as the proportion of barnacles with fertilized eggs or mature gonads) was analysed using generalized linear models (with a Poisson distribution in R) with salinity and population as the main factors in the models. Significance of these factors was evaluated using log likelihood tests of full versus reduced models.

To test the effects of salinity and population on the overall fitness of barnacles (based on all phenotypic traits measured in this study), PERMANOVA (two-factorial, orthogonal permutational MANOVA) was applied, using Euclidian distance matrices with 9999 permutations. We also ran PERMANOVA pair-wise post-hoc tests for multiple comparisons. Possible trade-offs between individual phenotypic traits in different salinities were visualized and evaluated using a canonical analysis of principal coordinates (CAP) [62, 63]. This type of analysis differs from unconstrained ordinations (e.g. MDS and PCA), because it identifies axes that maximize differences among groups rather than maximizing the variance explained, and thus allows a clearer focus on explicit hypotheses about e.g. salinity effects [63, 64].

Ethics statement

This work has been conducted according to relevant national and international guidelines for ethics and animal welfare, which do not include any specific requirement for barnacles.

Results

For all the response variables we investigated the effects of salinity and population, as well as their interaction (an indicator of local adaptation). None of the interactions were statistically significant, with the exception of shell strength (see below) (Additional file 1: Table S3, S4 and S5).

Growth

Salinity had a small but statistically significant effect on the size of barnacles (maximum shell diameter, as described by parameter a in the growth model; Figure 2) from all three tested populations (F = 3.362, P = 0.0497; Additional file 1: Table S3 and Figure 3a). Barnacles growing in 15 PSU had on average 8% larger shells than those in 30 PSU (Tukey’s HSD: P = 0.039), but were not significantly different from those in 6 PSU (Tukey’s HSD: P = 0.360; Figure 3a). The length of the early lag phase in growth (parameter b in growth model; Figure 2) was clearly influenced by salinity, where the start of growth of barnacles in 6 PSU was significantly delayed compared to the two higher salinities (15 and 30 PSU; F = 16.554, P < 0.001; Additional file 1: Table S3 and Figure 3b). There were no significant differences in maximum growth rate (parameter c; Figure 2) between salinities or between populations (Additional file 1: Table S3 and Figure 3c). Although sample sizes were small (n = 4 panels), the fitting of growth curves was good with average r2 = 0.994 ± 0.001 (SE) (Figure 2). The advantage of fitting a growth model described by several parameters, instead of simply analysing one endpoint, is that it provides a more detailed understanding of how salinity affects complex and dynamic growth patterns.

Figure 3
figure 3

Effects of salinity on growth. Effects of salinity on growth in barnacles from three different populations (Stockholm, Kiel and Tjärnö), described by the growth model parameters (a, b and c, labelled respectively) obtained from the Gompertz growth model fitted to barnacle shell diameter data (± standard error (SE), n = 4). Parameter a estimates the maximum shell diameter of the barnacles, b estimates the lag phase during early growth and c estimates the maximal growth rate, i.e. how rapidly barnacles approach the asymptote.

Ash weight (AW) was significantly influenced by salinity (F = 6.463, P < 0.005; Additional file 1: Table S4), where barnacles in 6 PSU had lower AW after 9 weeks, compared to barnacles in the higher salinities (15 and 30 PSU; Figure 4a). We found no statistically significant difference in AW between populations, although in 30 PSU barnacles from the Stockholm population had lower AW than the other populations (Tjärnö and Kiel; Figure 4a). AFDW did not differ significantly between salinities or populations (Additional file 1: Table S4 and Figure 4b).

Figure 4
figure 4

Effects of salinity on ash weight, ash-free dry weight and condition index. Effects of salinity on a) ash weight, b) ash-free dry weight, c) condition index (ratio between ash-free dry weight and dry weight) in barnacles from three different populations (Stockholm, Kiel and Tjärnö) that were grown in different salinity treatments (mean ± SE, n = 4).

Condition index and shell strength

The highest body-mass to shell-mass ratio (CI) was found in barnacles in low salinity (6 PSU), and CI decreased with increasing salinity (F = 8.658, P = 0.001; Additional file 1: Table S4 and Figure 4c). There were no differences in CI between populations (Additional file 1: Table S4). Furthermore, low salinity (6 PSU) resulted in weaker shells compared to higher salinities (15 and 30 PSU; F = 8.432, P = 0.001; Additional file 1: Table S4 and Figure 5). Importantly, however, there was a significant interaction between the effects of salinity and population on shell strength (F = 2.876, P = 0.042; Additional file 1: Table S4 and Figure 5). Barnacles from the high salinity environment (Tjärnö) had stronger shells in their native salinity (30 PSU) compared to the other two populations (Tukey’s HSD: P = 0.042; Figure 5). This positive relationship between shell strength and native salinity was not observed for the Kiel and Stockholm populations (Figure 5).

Figure 5
figure 5

Effects of salinity on shell strength. Effects of salinity on shell strength in barnacles from different populations (Stockholm, Kiel and Tjärnö) that were grown in different salinity treatments (mean ± SE, n = 4). Shell strength was standardised by dividing the force required to break the shell (in MPa) by the maximum shell diameter of each barnacles.

Reproductive maturity

Barnacles with mature gonads were observed in all three salinity treatments and in all three populations, but not in all combinations (Figure 6). Mature gonads were mostly observed in barnacles growing at intermediate and high salinities (15 and 30 PSU; Figure 6a). Barnacles brooding fertilized eggs, however, were only observed in the two lower salinities (6 and 15 PSU; Figure 6b). Generalized linear modelling revealed that mainly salinity but also to some extent population (but not their interaction) significantly influenced reproductive status (presence of fertilized eggs or mature gonads) (Additional file 1: Table S2). The analyses showed that the number of barnacles with fertilized eggs was best explained by a model including only salinity (P < 0.0001; Additional file 1: Table S2). This result was strongly driven by the absence of individuals with fertilized eggs in 30 PSU, (Figure 6b). For barnacles with mature gonads, but not carrying fertilized eggs, a model containing salinity and population gave the best fit (P < 0.003; Additional file 1: Table S2). For all populations, the greatest numbers of barnacles with mature gonads were observed in 30 PSU, followed by 15 PSU and 6 PSU. Barnacles from Kiel generally had the highest proportion of individuals with mature gonads in all salinities.

Figure 6
figure 6

Effects of salinity on reproduction. Effects of salinity on reproduction in barnacles from three different populations (Stockholm, Kiel and Tjärnö). a) Percentage of barnacles with mature, but not yet fertilized, eggs (mean ± SE, n = 4); b) Percentage of barnacles with fertilized eggs (mean ± SE, n = 4).

Mortality

Mortality was low in all treatments and populations, ranging from 2 to 15% (Figure 7). Mortality was significantly higher in 30 PSU than in 6 PSU (F = 5.597, P = 0.013; Additional file 1: Table S4), however, we detected no population-specific differences in mortality (Additional file 1: Table S4).

Figure 7
figure 7

Effects of salinity on mortality. Effects of salinity on mortality (week 2 to week 9) in barnacles from three different populations (Stockholm, Kiel and Tjärnö) that were grown in different salinity treatments (mean ± SE; n = 2–4).

Trade-offs between phenotypic traits in response to salinity

The CAP analysis that we used to explore effects of salinity and population on trade-offs between phenotypic traits revealed differences between salinity treatments, but not among populations (Figure 8). The first and the most important axis (δ2 = 68.1) separated the barnacles into three groups, based on salinity. This analysis confirmed that the traits that contributed most to this pattern were related mainly to shell strength and reproduction, although these were not 180 degrees opposed to each other as would be expected if the trade-off was complete. In high salinity (30 PSU), barnacles generally had stronger shells but delayed reproduction, whereas in low salinity (6 PSU), barnacles had weaker shells, but more fertilized eggs. In contrast, maximum size of barnacles (parameter a in the growth model) and AFDW did not contribute to the observed effects of salinity.

Figure 8
figure 8

Trade-offs between phenotypic traits in response to salinity. Bi-plot from canonical analysis of principal coordinates (CAP) on barnacles from three populations, exposed to different salinity treatments (grouped by salinity). The first canonical axis explained 68% of the variation whereas the second axis only explained 15%. Vectors (lines) indicate the influence of different phenotypic traits in explaining the differences in responses between salinities (correlations are scaled by multiplying the original value by 0.5 to fit plot size). The phenotypic traits included are: Gomp_a, Gomp_b, Gomp_c (Gompertz growth parameters; see methods), AW (=ash weight), DW (=dry weight), AFDW (=ash-free dry weight), CI (=condition index), shell strength, % Mature and % Fertilized (=proportion of barnacles with mature gonads vs. fertilized eggs).

Discussion

Our study shows that B. improvisus can tolerate a broad range of salinities and thus displays high phenotypic plasticity throughout its post-settlement life. Nonetheless, long-term exposure to different salinities had significant effects on fitness-related phenotypic traits such as shell-strength and fecundity. We also found evidence for a possible trade-off between traits (e.g. shell strength and reproduction), which may have important consequences in coping with changing environmental conditions. Overall responses to different salinities were similar in all populations, although differences in shell strength between populations in high salinity indicate possible local adaptation related to differences in natural predation pressure between the Baltic Sea and North Sea.

The broad salinity tolerance of B. improvisus was first observed by Darwin who found the species living in almost freshwater conditions in the La Plata River (Uruguay) during his voyage with the Beagle in 1831 [38]. Since then, several studies have documented the euryhaline abilities of B. improvisus[41, 43, 45, 65]. Few other barnacle species show similar tolerance to low salinities [66, 67], however most studies have investigated effects of short-term exposure, focusing on behavioural responses in adults, e.g. [45] or larval development, e.g. [42]. Few studies have investigated long-term effects of salinity on the post-larval life stages of barnacles [11, 68]. Our study therefore contributes important new knowledge about long-term fitness-related effects of salinity exposure during the juvenile post-settlement stage until reproductive maturity.

Despite better performance in low and intermediate salinities in terms of e.g. maximum size and survival (Figures 3a and 7), our experiments revealed that long-term exposure to low salinity had negative effects on other fitness-related traits in B. improvisus. Growth during the first weeks after metamorphosis is recognized as a vital phase of a barnacle’s life [69], and therefore our observation of delayed onset of growth in barnacles at 6 PSU (Figure 3b) may have important consequences. Turpayeva and Simkina [43] observed a similar delay in early growth in response to low salinity in B. improvisus. Fitness advantages of rapid initiation of growth include improved competition for space [70], increased tolerance to fluctuating environmental stressors [71, 72] as well as size-based predator avoidance [73]. Reduced growth under environmental stress has been shown in other marine species [74, 75]. The delay in early growth that we observed in 6 PSU may be due to altered physiology (e.g. increased osmoregulatory needs), resulting in re-allocation of energy [41] and/or increased shell production costs in low salinities due to low CaCO3 content in the water [76].

Salinity had a clear effect on reproduction in all populations with the highest reproductive activity (sum of individuals with fertilized eggs and/or mature gonads) in 15 PSU. Although mature gonads were found in barnacles in 30 PSU, no individuals with fertilized eggs were observed even after nine weeks of study, indicating that reproduction may be delayed in high salinity. This may be a result of more energy initially being invested in shell production in high salinities, revealing a trade-off between two important fitness-related traits: reproduction and defence against predators by producing hard shells (i.e. survival). A similar trade-off between shell morphology and reproduction has also been reported by [77] for the barnacle Chthamalus anisopoma. Alternatively, it is possible that reproduction is delayed or disrupted in high salinities, i.e. not representing an adaptive trade-off per se, but rather the result of physiological dysfunction. In either case, delaying reproduction will result in fewer broods per year, and thereby affect recruitment dynamics [78, 79], however the importance of this delay may be overwhelmed by the effects of season on reproductive activity. Studies on other barnacle species have suggested that exposure to low salinities markedly reduces the reproductive output [80, 81]. This does not seem to be the case for B. improvisus, providing further support for the conclusion that B. improvisus is a truly brackish species.

The calcareous outer shell of a barnacle offers protection against predation but also aids in regulating exposure to extreme conditions such as desiccation or salinity changes [82, 83]. The reduction in shell strength we observed at low salinities (Figure 5 and Additional file 1: Table S4) was most likely caused by limited amounts of dissolved CaCO3 at low salinity [84, 85]. Reduced shell strength and ash weight in response to low salinity has also been found in other invertebrate species, including mussels (Mytilus edulis) and oysters [86, 87]. Extreme high salinities (>40 PSU) were also found to reduce shell strength of B. amphitrite[55], indicating that shell strength may not only be limited by available CaCO3 in the water, but may also be associated with other physiological stress in response to changes in salinity.

Interestingly, we found that barnacles from the highest salinity environment (Tjärnö) built stronger shells in 30 PSU compared to the other two populations (Figure 5 and Additional file 1: Table S4). The absence of many of B. improvisus’ main predators (e.g. Nucella lapillus, Carcinus maenas and Asterias rubens) from lower salinity environments in the inner Baltic Sea [87, 88] may have resulted in lower selection pressure for strong shells, especially when CaCO3 is limited. Similar, population-specific differences in the ability to build thick shells have been shown for blue mussels (Mytilus edulis) from the Baltic and North Seas [88]. Interestingly, barnacles from Kiel (15 PSU) only showed a modest increase in shell strength in 30 PSU, compared to barnacles from Tjärnö (Figure 5 and Additional file 1: Table S4) even though the predators A. rubens and C. maenas are abundant in that habitat [89]. We can only speculate with regard to the predation pressure that these species exert on B. improvisus in Kiel, however it is clear that other factors such as maternal effects [90], trade-offs between shell-strength and reproductive output (CAP analysis, Figure 8), and physiological constraints may also have influenced our results. Furthermore, the influence of gene flow from populations further inside the Baltic Sea (with lower predation pressure) to Kiel is estimated to be higher than the gene flow from Tjärnö to Kiel, based on oceanographic modelling, which could partly explain the observed pattern [91, 92].

With the exception of the population-specific responses in shell strength, no evidence of local adaptation was detected in the other measured traits. Several scenarios could explain this result: i) Balanus improvisus was first observed in the Baltic Sea relatively recently (<200 years ago [17, 93], and this may have been insufficient time for new adaptations to arise – especially if standing genetic variation was low in the founding populations, which is often the case after biological invasions [9496]; ii) B. improvisus has free-swimming pelagic larvae that drift with ocean currents for up to several weeks, potentially creating high levels of gene flow that prevent local adaptation from evolving [97, 98]. In addition, B. improvisus is a major fouling organism on ships and human mediated dispersal may enhance gene flow between populations [99]. However, there are several examples in the literature that show that adaptations can evolve rapidly [100, 101] and that adaptations can be maintained despite high gene flow [25, 102, 103]; iii) B. improvisus displays broad phenotypic plasticity, which can facilitate establishment in new environments without the need for strong selection/adaptation of local populations [104].

Phenotypic plasticity has also been suggested to play an important role in organisms’ ability to cope with current and future climate change [105]. Recent modelling of future climate scenarios in the Baltic Sea indicate that increased precipitation may lead to reduced salinity [23], which could potentially result in dramatic shifts in species’ distributions. According to the most extreme climate scenario of Meier et al.[23], the surface salinity of 5 PSU, today situated in the northern Baltic Sea (63° N), would move south to the waters around Bornholm (55° N) by the end of the 21st century (Figure 1) Although hypothetical [3], this type of environmental shift would have dramatic consequences, especially in ecosystems with already low species richness and genetic diversity, such as the Baltic Sea. Our results indicate that it is unlikely that climate-related shifts in salinity will have strong negative effects on B. improvisus populations, but may rather favour the species. It should be remembered, however, that climate change involves multiple environmental parameters, as well as shifts in ecosystem composition, and it is therefore hard to predict the likely effects of salinity changes on the ecology of B. improvisus or the ecosystems in which it lives.

Conclusions

In conclusion, we have shown that B. improvisus is a highly euryhaline species with strong capacity to tolerate a range of salinities through primarily plastic responses. For almost all response variables, B. improvisus performed slightly better at low and intermediate salinities, supporting the idea that B. improvisus is one of few truly brackish species. However, the negative effects of low salinity on early initiation of growth as well as reduced shell strength suggest trade-offs between traits in different environments in order to maximise fitness. Population-specific responses in shell-strength indicated the possibility of local adaptation, perhaps in relation to different predation pressures along the salinity gradient. The existence of these population-specific responses, despite recent colonization, high potential dispersal and broad tolerance, supports earlier work indicating that evolutionary changes can occur rapidly. Projected future climate-driven reductions in salinity in the Baltic Sea will most likely not have major impacts on B. improvisus populations, however, further work is needed to clarify the interactions between salinity tolerance and other stressors such as temperature, acidification and food limitation [106, 107]. Furthermore, selection experiments using multiple generations of B. improvisus are needed to elucidate the respective roles of phenotypic plasticity, trans-generational effects, and adaptations in response to locally strong selection pressures. Finally, population genomic studies in the Baltic Sea could help to elucidate population structure and identify (candidate) genes involved in both phenotypic plasticity and local adaptations.

References

  1. Antonov JI, Levitus S, Boyer TP: Steric sea level variations during 1957–1994: Importance of salinity. J Geophysl Res C. 2002, 107: 8013-

    Google Scholar 

  2. Boyer TP, Levitus S, Antonov JI, Locarnini RA, Garcia HE: Linear trends in salinity for the world ocean, 1955–1998. Geophys Res Lett. 2005, 32: L01604-

    Google Scholar 

  3. Meier HEM, Andersson HC, Arheimer B, Blenckner T, Chubarenko B, Donnelly C, Eilola K, Gustafsson BG, Hansson A, Havenhand J, Höglund A, Kuznetsov I, MacKenzie BR, Müller-Karulis B, Neumann T, Niiranen S, Piwowarczyk J, Raudsepp U, Reckermann M, Ruoho-Airola T, Savchuk OP, Schenk F, Schimanke S, Väli G, Weslawski J-M, Zorita E: Comparing reconstructed past variations and future projections of the Baltic Sea ecosystem-first results from multi-model ensemble simulations. Environ Res Lett. 2012, 7: 034005-

    Google Scholar 

  4. Byrne M: Impacts of Ocean Warming and Ocean Acidification on Marine Invertebrate Life History Stages: Vulnerabilities and Potential For Persistence in a Changing Ocean. Oceanography and Marine Biology: An Annual Review, Volume 49. Edited by: Gibson RN, Atkinson RJA, Gordon JDM. 2011, Boca Raton: Crc Press-Taylor & Francis Group, 1-42.

    Google Scholar 

  5. Doney SC, Fabry VJ, Feely RA, Kleypas JA: Ocean acidification: the other CO2 problem. Ann Rev Mar Sci. 2009, 1: 169-192.

    PubMed  Google Scholar 

  6. Kroeker KJ, Kordas RL, Crim R, Hendriks IE, Ramajo L, Singh GS, Duarte CM, Gattuso JP: Impacts of ocean acidification on marine organisms: quantifying sensitivities and interaction with warming. Glob Chang Biol. 2013, 19: 1884-1896.

    PubMed  PubMed Central  Google Scholar 

  7. Wikner J, Andersson A: Increased freshwater discharge shifts the trophic balance in the coastal zone of the northern Baltic Sea. Glob Chang Biol. 2012, 18: 2509-2519.

    PubMed Central  Google Scholar 

  8. Kinne RKH: The role of organic osmolytes in osmoregulation: From bacteria to mammals. J Exp Zool. 1993, 265: 346-355.

    PubMed  CAS  Google Scholar 

  9. Lapucki T, Nonnant M: Physiological responses to salinity changes of the isopod Idotea chelipes from the Baltic brackish waters. Comp Biochem Physiol A-Mol Integr Physiol. 2008, 149: 299-305.

    PubMed  Google Scholar 

  10. Normant M, Lamprecht I: Does scope for growth change as a result of salinity stress in the amphipod Gammarus oceanicus?. J Exp Mar Biol Ecol. 2006, 334: 158-163.

    Google Scholar 

  11. Qiu JW, Qian PY: Tolerance of the barnacle Balanus amphitrite amphitrite (Darwin) to salinity and temperature stress: effects of previous experience. Mar Ecol-Prog Ser. 1999, 188: 123-132.

    Google Scholar 

  12. Torres G, Gimenez L, Anger K: Growth, tolerance to low salinity, and osmoregulation in decapod crustacean larvae. Aquat Biol. 2011, 12: 249-260.

    Google Scholar 

  13. Bonsdorff E: Zoobenthic diversity-gradients in the Baltic Sea: Continuous post-glacial succession in a stressed ecosystem. J Exp Mar Biol Ecol. 2006, 330: 383-391.

    Google Scholar 

  14. Ojaveer H, Jaanus A, MacKenzie BR, Martin G, Olenin S, Radziejewska T, Telesh I, Zettler ML, Zaiko A: Status of biodiversity in the Baltic Sea. PLoS One. 2010, 5: e12467-

    PubMed  PubMed Central  Google Scholar 

  15. Björck S: A review of the history of the Baltic Sea, 13.0-8.0 KA BP. Quat Int. 1995, 27: 19-40.

    Google Scholar 

  16. Zillén L, Conley DJ, Andren T, Andren E, Björck S: Past occurrences of hypoxia in the Baltic Sea and the role of climate variability, environmental change and human impact. Earth Sci Rev. 2008, 91: 77-92.

    Google Scholar 

  17. Leppäkoski E, Olenin S: Non-native species and rates of spread: Lessons from the brackish Baltic Sea. Biol Invas. 2000, 2: 151-163.

    Google Scholar 

  18. Leppäranta M, Myrberg K: Physical Oceanography of the Baltic Sea. 2009, Chichester, UK: Springer Praxis Publishing Ltd

    Google Scholar 

  19. Omstedt A, Axell LB: Modeling the variations of salinity and temperature in the large Gulfs of the Baltic Sea. Cont Shelf Res. 2003, 23: 265-294.

    Google Scholar 

  20. Bergström L, Kautsky L: Local adaptation of Ceramium tenuicorne (Ceramiales, Rhodophyta) within the Baltic Sea. J Phycol. 2006, 42: 36-42.

    Google Scholar 

  21. Johannesson K, André C: Life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea. Mol Ecol. 2006, 15: 2013-2029.

    PubMed  CAS  Google Scholar 

  22. Johannesson K, Smolarz K, Grahn M, Andre C: The future of Baltic Sea populations: local extinction or evolutionary rescue?. Ambio. 2011, 40: 179-190.

    PubMed  CAS  PubMed Central  Google Scholar 

  23. Meier HEM, Kjellstrom E, Graham LP: Estimating uncertainties of projected Baltic Sea salinity in the late 21st century. Geophys Res Lett. 2006, 33: L15705-

    Google Scholar 

  24. Pigliucci M: Phenotypic Plasticity - Beyond Nature and Nurture. 2001, Baltimore: The John Hopkins University Press

    Google Scholar 

  25. Reusch TBH: Climate change in the oceans: evolutionary versus phenotypically plastic responses of marine animals and plants. Evol Appl. 2014, 7: 104-122.

    PubMed  PubMed Central  Google Scholar 

  26. Crispo E: Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow. J Evol Biol. 2008, 21: 1460-1469.

    PubMed  CAS  Google Scholar 

  27. Ghalambor CK, McKay JK, Carroll SP, Reznick DN: Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007, 21: 394-407.

    Google Scholar 

  28. Waddington CH: Genetic assimilation of an acquired character. Evolution. 1953, 7: 118-126.

    Google Scholar 

  29. Gienapp P, Teplitsky C, Alho JS, Mills JA, Merila J: Climate change and evolution: disentangling environmental and genetic responses. Mol Ecol. 2008, 17: 167-178.

    PubMed  CAS  Google Scholar 

  30. DeWitt TJ, Sih A, Wilson DS: Costs and limits of phenotypic plasticity. Trends Ecol Evol. 1998, 13: 77-81.

    PubMed  CAS  Google Scholar 

  31. Releya RA: Costs of phenotypic plasticity. Am Nat. 2002, 159: 272-282.

    Google Scholar 

  32. Lande R: Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation. J Evol Biol. 2009, 22: 1435-1446.

    PubMed  Google Scholar 

  33. Pigliucci M, Murren CJ, Schlichting CD: Phenotypic plasticity and evolution by genetic assimilation. J Exp Biol. 2006, 209: 2362-2367.

    PubMed  Google Scholar 

  34. Via S: Adaptive phenotypic plasticity - target or by-product of selection in a variable environment. Am Nat. 1993, 142: 352-365.

    PubMed  CAS  Google Scholar 

  35. Sotka EE: Natural selection, larval dispersal, and the geography of phenotype in the sea. Integr Comp Biol. 2012, 52: 538-545.

    PubMed  Google Scholar 

  36. Renborg E, Johannesson K, Havenhand J: Variable salinity tolerance in ascidian larvae is primarily a plastic response to the parental environment. Evol Ecol. 2014, 28: 561-572.

    Google Scholar 

  37. Foster BA: Barnacle Ecology and Adaptation. Barnacle Biology. Edited by: Southward AJ. 1987, Rotterdam: Balkema, [Schram FR (Series Editor): Crustacean Issues]

    Google Scholar 

  38. Darwin C: A Monograph on the Sub-Class Cirripedia: The Balanidae (or sessile cirrepedes); the Verrucidae. 1854, London: The Ray Society

    Google Scholar 

  39. De Rivera CE, Steves BP, Fofonoff PW, Hines AH, Ruiz GM: Potential for high-latitude marine invasions along western North America. Divers Distrib. 2011, 17: 1198-1209.

    Google Scholar 

  40. Blom S-E: Balanus improvisus Darwin on the west coast of Sweden. Zoologiska Bidrag från Uppsala. 1965, 37: 59-76.

    Google Scholar 

  41. Fyhn HJ: Holeuryhalinity and its mechanisms in a cirriped crustacean, Balanus improvisus. Comp Biochem Physiol A Physiol. 1976, 53: 19-30.

    CAS  Google Scholar 

  42. Nasrolahi A, Pansch C, Lenz M, Wahl M: Being young in a changing world: how temperature and salinity changes interactively modify the performance of larval stages of the barnacle Amphibalanus improvisus. Mar Biol. 2012, 159: 331-340.

    Google Scholar 

  43. Turpayeva YP, Simkina RG: The reaction of the Black Sea Balanus improvisus (Darwin) to reduced salinities. T Inst Oceanogr. 1961, XLIX: 205-223.

    Google Scholar 

  44. Lind U, Rosenblad MA, Wrange AL, Sundell KS, Jonsson PR, Andre C, Havenhand J, Blomberg A: Molecular characterization of the alpha-subunit of Na+/K+ ATPase from the euryhaline barnacle Balanus improvisus reveals multiple genes and differential expression of alternative splice variants. PLoS One. 2013, 8: e77069-

    PubMed  CAS  PubMed Central  Google Scholar 

  45. Davenport J: A comparative study of the behaviour of some balanomorph barnacles exposed to fluctuating sea water concentrations. J Mar Biol Ass U K. 1976, 56: 889-907.

    Google Scholar 

  46. Wahl M, Shahnaz L, Dobretsov S, Saha M, Symanowski F, David K, Lachnit T, Vasel M, Weinberger F: Ecology of antifouling resistance in the bladder wrack Fucus vesiculosus: patterns of microfouling and antimicrobial protection. Mar Ecol Progr Ser. 2010, 411: 33-U61.

    Google Scholar 

  47. Berntsson KM, Jonsson PR: Temporal and spatial patterns in recruitment and succession of a temperate marine fouling assemblage: A comparison of static panels and boat hulls during the boating season. Biofouling. 2003, 19: 187-195.

    PubMed  CAS  Google Scholar 

  48. Anil AC, Chiba K, Okamoto K, Kurokura H: Influence of temperature and salinity on larval development of Balanus amphitrite - implications in fouling ecology. Mar Ecol Prog Ser. 1995, 118: 159-166.

    Google Scholar 

  49. Berntsson KM, Jonsson PR, Lejhall M, Gatenholm P: Analysis of behavioural rejection of micro-textured surfaces and implications for recruitment by the barnacle Balanus improvisus. J Exp Mar Biol Ecol. 2000, 251: 59-83.

    PubMed  Google Scholar 

  50. Pansch C, Schlegel P, Havenhand J: Larval development of the barnacle Amphibalanus improvisus responds variably but robustly to near-future ocean acidification. ICES J Mar Sci: Journal du Conseil. 2013, 70: 805-811.

    Google Scholar 

  51. Quinn G, Keough M: Experimental Design and Data Analysis for Biologists. 2002, Cambridge, UK: Cambridge University Press

    Google Scholar 

  52. Aminot A, Rey F: Standard procedure for the determination of chlorophyll a by spectroscopic methods. ICES: Techniques in Marine Environmental Sciences. 2002, 28: 1-15.

    Google Scholar 

  53. von Bertalanffy L: A quantitative theory of organic growth. Hum Biol. 1938, 10: 181-213.

    CAS  Google Scholar 

  54. Winsor C: The Gompertz curve as a growth equation. Proc Natl Acad Sci U S A. 1932, 18: 1-8.

    PubMed  CAS  PubMed Central  Google Scholar 

  55. Simpson E, Hurlbert S: Salinity effects on the growth, mortality and shell strength of Balanus amphitrite from the Salton Sea, California. Hydrobiologia. 1998, 381: 179-190.

    Google Scholar 

  56. Tsang WN, Chaillé PM, Collins PM: Growth and reproductive performance in cultured nearshore rockfish (Sebastes spp.). Aquaculture. 2007, 266: 236-245.

    Google Scholar 

  57. Sköld M, Josefson AB, Loo LO: Sigmoidal growth in the brittle star Amphiura filiformis (Echinodermata : Ophiuroidea). Mar Biol. 2001, 139: 519-526.

    Google Scholar 

  58. Karkach AS: Trajectories and models of individual growth. Demogr Res. 2006, 15: 348-400.

    Google Scholar 

  59. Pansch C, Nasrolahi A, Appelhans YS, Wahl M: Tolerance of juvenile barnacles (Amphibalanus improvisus) to warming and elevated pCO2. Mar Biol. 2013, 160: 2023-2035.

    CAS  Google Scholar 

  60. Range P, Chicharo MA, Ben-Hamadou R, Pilo D, Matias D, Joaquim S, Oliveira AP, Chicharo L: Calcification, growth and mortality of juvenile clams Ruditapes decussatus under increased pCO2 and reduced pH: Variable responses to ocean acidification at local scales?. J Exp Mar Biol Ecol. 2011, 396: 177-184.

    Google Scholar 

  61. Kawecki TJ, Stearns SC: The evolution of life histories in spatially heterogeneous environments - optimal reaction norms revisited. Evol Ecol. 1993, 7: 155-174.

    Google Scholar 

  62. Anderson MJ: CAP: A FORTRAN Computer Program for Canonical Analysis of Principal Coordinates. 2004, New Zealand: Department of Statistics, University of Auckland

    Google Scholar 

  63. Anderson MJ, Willis TJ: Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology. 2003, 84: 511-525.

    Google Scholar 

  64. Meot A, Legendre P, Borcard D: Partialling out the spatial component of ecological variation: questions and propositions in the linear modelling framework. Environ Ecol Stat. 1998, 5: 1-27.

    Google Scholar 

  65. Foster BA: Responses and acclimation to salinity in adults of some balanomorph barnacles. Philos Trans R Soc Lond Ser B-Biol Sci. 1970, 256: 377-400.

    Google Scholar 

  66. Dineen JF, Hines AH: Effects of salinity and adult extract on settlement of the oligohaline barnacle Balanus subalbidus. Mar Biol. 1994, 119: 423-430.

    Google Scholar 

  67. Dineen JF, Hines AH: Larval settlement of the polyhaline barnacle Balanus eburneus (Gould): cue interactions and comparisons with two estuarine congeners. J Exp Mar Biol Ecol. 1994, 179: 223-234.

    Google Scholar 

  68. Nasrolahi A, Pansch C, Lenz M, Wahl M: Temperature and salinity interactively impact early juvenile development: a bottleneck in barnacle ontogeny. Mar Biol. 2013, 160: 1109-1117.

    Google Scholar 

  69. Gosselin LA, Qian PY: Early post-settlement mortality of an intertidal barnacle: A critical period for survival. Mar Ecol-Progr Ser. 1996, 135: 69-75.

    Google Scholar 

  70. Connell JH: Effects of competition, predation by Thais lapillus, and other factors on natural populations of barnacle Balanus balanoides. Ecol Monogr. 1961, 31: 61-104.

    Google Scholar 

  71. Foster BA: Desiccation as a factor in intertidal zonation of barnacles. Mar Biol. 1971, 8: 12-29.

    Google Scholar 

  72. Gosselin LA, Qian PY: Juvenile mortality in benthic marine invertebrates. Mar Ecol Prog Ser. 1997, 146: 265-282.

    Google Scholar 

  73. Miller KM, Carefoot TH: The role of spatial and size refuges in the interaction between juvenile barnacles and grazing limpets. J Exp Mar Biol Ecol. 1989, 134: 157-174.

    Google Scholar 

  74. Anger K, Spivak E, Luppi T: Effects of reduced salinities on development and bioenergetics of early larval shore crab, Carcinus maenas. J Exp Mar Biol Ecol. 1998, 220: 287-304.

    Google Scholar 

  75. Torres G, Gimenez L, Anger K: Cumulative effects of low salinity on larval growth and biochemical composition in an estuarine crab, Neohelice granulata. Aquat Biol. 2008, 2: 37-45.

    Google Scholar 

  76. Malone PG, Dodd JR: Temperature and salinity effects on calcification rate in Mytilus edulis and its paleoecological implications. Limnol Oceanogr. 1967, 12: 432-436.

    CAS  Google Scholar 

  77. Lively CM: Predator-induced shell dimorphism in the acorn barnacle Chthamalus anisopoma. Evolution. 1986, 40: 232-242.

    Google Scholar 

  78. Berger MS: Reproduction of the intertidal barnacle Balanus glandula along an estuarine gradient. Mar Ecol. 2009, 30: 346-353.

    Google Scholar 

  79. Hines AH: Reproduction in three species of intertidal barnacles from central California. Biol Bull. 1978, 154: 262-281.

    Google Scholar 

  80. Barnes H, Barnes M: Egg numbers, metabolic efficiency of egg production and fecundity; Local and regional variations in a number of common cirripedes. J Exp Mar Biol Ecol. 1968, 2: 135-153.

    Google Scholar 

  81. Crisp DJ, Costlow JD: The tolerance of developing cirripede embryos to salinity and temperature. Oikos. 1963, 14: 22-34.

    Google Scholar 

  82. Gosselin LA, Jones LA: Effects of solar radiation on barnacle settlement, early post-settlement mortality and community development in the intertidal zone. Mar Ecol Progr Ser. 2010, 407: 149-158.

    Google Scholar 

  83. Southward AJ: Barnacle Biology. 1987, Rotterdam: A.A. Balkema

    Google Scholar 

  84. Dickinson GH, Ivanina AV, Matoo OB, Portner HO, Lannig G, Bock C, Beniash E, Sokolova IM: Interactive effects of salinity and elevated CO2 levels on juvenile eastern oysters, Crassostrea virginica. J Exp Biol. 2012, 215: 29-43.

    PubMed  CAS  Google Scholar 

  85. Mucci A: The solubility of calcite and aragonite in seawater at various salinities, temperatures and one atmosphere total pressure. Am J Sci. 1983, 283: 780-799.

    CAS  Google Scholar 

  86. Almada-Villela PC: The effects of reduced salinity on the shell growth of small Mytilus edulis. J Mar Biol Assoc UK. 1984, 64: 171-182.

    Google Scholar 

  87. Kautsky N, Johannesson K, Tedengren M: Genotypic and phenotypic differences between Baltic and North-Sea populations of Mytilus edulis evaluated through reciprocal transplantations. 1. Growth and morphology. Mar Ecol Prog Ser. 1990, 59: 203-210.

    Google Scholar 

  88. Reimer O, Harms-Ringdahl S: Predator-inducible changes in blue mussels from the predator-free Baltic Sea. Mar Biol. 2001, 139: 959-965.

    Google Scholar 

  89. Reusch TBH, Chapman ARO: Persistence and space occupancy by subtidal blue mussel patches. Ecol Monogr. 1997, 67: 65-87.

    Google Scholar 

  90. Marshall DJ, Allen RM, Crean AJ: The Ecological and Evolutionary Importance of Maternal Effects in the Sea. Oceanography and Marine Biology: An Annual Review, Volume 46. Edited by: Gibson RN, Atkinson RJA, Gordon JDM. 2008, Boca Raton: Crc Press-Taylor & Francis Group, 203-250.

    Google Scholar 

  91. Teacher AGF, Andre C, Jonsson PR, Merila J: Oceanographic connectivity and environmental correlates of genetic structuring in Atlantic herring in the Baltic Sea. Evol Appl. 2013, 6: 549-567.

    PubMed  PubMed Central  Google Scholar 

  92. Wrange A-L: From genes to geography - evolutionary perspectives on salinity tolerance in the brackish water barnacle Balanus improvisus. PhD thesis. 2014, University of Gothenburg: Department of Biological and Environmental Sciences

    Google Scholar 

  93. Luther A: Om Balanus improvisus i Östersjön. Fauna Flora. 1950, 45: 155-160.

    Google Scholar 

  94. Estoup A, Guillemaud T: Reconstructing routes of invasion using genetic data: why, how and so what?. Mol Ecol. 2010, 19: 4113-4130.

    PubMed  Google Scholar 

  95. Facon B, Genton BJ, Shykoff J, Jarne P, Estoup A, David P: A general eco-evolutionary framework for understanding bioinvasions. Trends Ecol Evol. 2006, 21: 130-135.

    PubMed  Google Scholar 

  96. Lambrinos JG: How interactions between ecology and evolution influence contemporary invasion dynamics. Ecology. 2004, 85: 2061-2070.

    Google Scholar 

  97. Caley MJ, Carr MH, Hixon MA, Hughes TP, Jones GP, Menge BA: Recruitment and the local dynamics of open marine populations. Annu Rev Ecol Syst. 1996, 27: 477-500.

    Google Scholar 

  98. Kawecki TJ, Ebert D: Conceptual issues in local adaptation. Ecol Lett. 2004, 7: 1225-1241.

    Google Scholar 

  99. Seebens H, Gastner MT, Blasius B: The risk of marine bioinvasion caused by global shipping. Ecol Lett. 2013, 16: 782-790.

    PubMed  CAS  Google Scholar 

  100. Bell MA, Aguirre WE, Buck NJ: Twelve years of contemporary armor evolution in a threespine stickleback population. Evolution. 2004, 58: 814-824.

    PubMed  Google Scholar 

  101. Wirgin I, Roy NK, Loftus M, Chambers RC, Franks DG, Hahn ME: Mechanistic basis of resistance to PCBs in Atlantic tomcod from the Hudson River. Science. 2011, 331: 1322-1325.

    PubMed  CAS  PubMed Central  Google Scholar 

  102. Hemmer-Hansen J, Nielsen EE, Frydenberg J, Loeschcke V: Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.). Heredity. 2007, 99: 592-600.

    PubMed  CAS  Google Scholar 

  103. Schmidt PS, Bertness MD, Rand DM: Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides. Proc R Soc Lond Ser B-Biol Sci. 2000, 267: 379-384.

    CAS  Google Scholar 

  104. Lee CE, Remfert JL, Gelembiuk GW: Evolution of physiological tolerance and performance during freshwater invasions. Integr Comp Biol. 2003, 43: 439-449.

    PubMed  Google Scholar 

  105. Munday PL, Warner RR, Monro K, Pandolfi JM, Marshall DJ: Predicting evolutionary responses to climate change in the sea. Ecol Lett. 2013, 16: 1488-1500.

    PubMed  Google Scholar 

  106. Lee CE, Moss WE, Olson N, Chau KF, Chang Y-M, Johnson KE: Feasting in fresh water: impacts of food concentration on freshwater tolerance and the evolution of food × salinity response during the expansion from saline into fresh water habitats. Evol Appl. 2013, 6: 673-689.

    PubMed  PubMed Central  Google Scholar 

  107. Todgham AE, Stillman JH: Physiological responses to shifts in multiple environmental stressors: relevance in a changing world. Integr Comp Biol. 2013, 53: 539-544.

    PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by a Linnaeus-grant from the Swedish Research Councils, VR and Formas (http://www.cemeb.science.gu.se) as well as the Lars Hierta Memorial Foundation. We also want to thank the technical staff and lab assistants at the Lovén Centre - Tjärnö for assisting during the experiment (especially Tomas Johansson, Martin Ogemark, Emma Johansson, Christopher Weller, Daniel Simonsson and Anne Thonig) as well as the Archipelago Foundation in Stockholm (Skärgårdsstiftelsen) for providing animals for the barnacle larval cultures. A special thanks to Gerry Quinn (Deakin University) for statistical advice and to Christian Pansch, GEOMAR, Kiel, for assistance with shell strength measurements as well as giving useful comments on the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Anna-Lisa Wrange.

Additional information

Competing interests

The author’s declare that they have no competing interests.

Authors’ contributions

ALW was involved in planning of the experimental design, carried out the experiment and collected the data, ran the statistical analyses and drafted the manuscript. CA participated in designing the experiment, analysing the data and commenting on the manuscript, TL developed the mathematical script for analysis of growth data and commented on the manuscript, UL and AB contributed to the experimental planning and commented on the manuscript, PJ took part in the planning, analyses and commenting on the manuscript, JH was involved in the planning, experimental setup, statistical analysis and commenting on the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2014_2646_MOESM1_ESM.pdf

Additional file 1: Table S1: Experimental conditions of the common-garden setup. Table S2. Effects of salinity and population on reproduction (frequency of fertilized eggs and mature gonads. Table S3. Results of ANOVA testing effects of salinity and population on Gompertz growth model parameters. Table S4. Results of ANOVA testing effects of salinity and population on fitness-related traits in barnacles; including ash weight, ash-free dry weight, condition index, shell strength, and mortality). Table S5. Results from the PERMANOVA on effects of salinity and population on phenotypic traits in barnacles, including growth, shell strength, condition index and reproductive maturity. (PDF 158 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wrange, AL., André, C., Lundh, T. et al. Importance of plasticity and local adaptation for coping with changing salinity in coastal areas: a test case with barnacles in the Baltic Sea. BMC Evol Biol 14, 156 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2148-14-156

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2148-14-156

Keywords