Spatial genetic structure in a crustacean herbivore highlights the need for local considerations in Baltic Sea biodiversity management

A decade of progress in marine evolutionary biology

This article summarizes the Evolutionary Applications Special Issue, "A decade of progress in Marine Evolutionary Biology." The globally connected ocean, from its pelagic depths to its highly varied coastlines, inspired Charles Darwin to develop the theory of evolution during the voyage of the Beagle. As technology has developed, there has been a dramatic increase in our knowledge about life on our blue planet. This Special Issue, composed of 19 original papers and seven reviews, represents a small contribution to the larger picture of recent research in evolutionary biology, and how such advancements come about through the connection of researchers, their fields, and their knowledge. The first European network for marine evolutionary biology, the Linnaeus Centre for Marine Evolutionary Biology (CeMEB), was developed to study evolutionary processes in the marine environment under global change. Though hosted by the University of Gothenburg in Sweden, the network quickly grew to encompass researchers throughout Europe and beyond. Today, more than a decade after its foundation, CeMEB's focus on the evolutionary consequences of global change is more relevant than ever, and knowledge gained from marine evolution research is urgently needed in management and conservation. This Special Issue, organized and developed through the CeMEB network, contains contributions from all over the world and provides a snapshot of the current state of the field, thus forming an important basis for future research directions.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

Adaptive, maladaptive, neutral, or absent plasticity: Hidden caveats of reaction norms

Adaptive phenotypic plasticity may improve the response of individuals when faced with new environmental conditions. Typically, empirical evidence for plasticity is based on phenotypic reaction norms obtained in reciprocal transplant experiments. In such experiments, individuals from their native environment are transplanted into a different environment, and a number of trait values, potentially implicated in individuals' response to the new environment, are measured. However, the interpretations of reaction norms may differ depending on the nature of the assessed traits, which may not be known beforehand. For example, for traits that contribute to local adaptation, adaptive plasticity implies nonzero slopes of reaction norms. By contrast, for traits that are correlated to fitness, high tolerance to different environments (possibly due to adaptive plasticity in traits that contribute to adaptation) may, instead, result in flat reaction norms. Here we investigate reaction norms for adaptive versus fitness-correlated traits and how they may affect the conclusions regarding the contribution of plasticity. To this end, we first simulate range expansion along an environmental gradient where plasticity evolves to different values locally and then perform reciprocal transplant experiments in silico. We show that reaction norms alone cannot inform us whether the assessed trait exhibits locally adaptive, maladaptive, neutral, or no plasticity, without any additional knowledge of the traits assessed and species' biology. We use the insights from the model to analyse and interpret empirical data from reciprocal transplant experiments involving the marine isopod Idotea balthica sampled from two geographical locations with different salinities, concluding that the low-salinity population likely has reduced adaptive plasticity relative to the high-salinity population. Overall, we conclude that, when interpreting results from reciprocal transplant experiments, it is necessary to consider whether traits assessed are locally adaptive with respect to the environmental variable accounted for in the experiments or correlated to fitness.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Ecological connectivity of the marine protected area network in the Baltic Sea, Kattegat and Skagerrak: Current knowledge and management needs

Marine protected areas (MPAs) have become a key component of conservation and fisheries management to alleviate anthropogenic pressures. For MPA networks to efficiently promote persistence and recovery of populations, ecological connectivity, i.e. dispersal and movement of organisms and material across ecosystems, needs to be taken into account. To improve the ecological coherence of MPA networks, there is hence a need to evaluate the connectivity of species spreading through active migration and passive dispersal. We reviewed knowledge on ecological connectivity in the Baltic Sea, Kattegat and Skagerrak in the northeast Atlantic and present available information on species-specific dispersal and migration distances. Studies on genetic connectivity are summarised and discussed in relation to dispersal-based analyses. Threats to ecological connectivity, limiting dispersal of populations and lowering the resilience to environmental change, were examined. Additionally, a review of studies evaluating the ecological coherence of MPA networks in the Baltic Sea, Kattegat and Skagerrak was performed, and suggestions for future evaluations to meet management needs are presented.

Biophysical models of dispersal contribute to seascape genetic analyses

Dispersal is generally difficult to directly observe. Instead, dispersal is often inferred from genetic markers and biophysical modelling where a correspondence indicates that dispersal routes and barriers explain a significant part of population genetic differentiation. Biophysical models are used for wind-driven dispersal in terrestrial environments and for propagules drifting with ocean currents in the sea. In the ocean, such seascape genetic or seascape genomic studies provide promising tools in applied sciences, as actions within management and conservation rely on an understanding of population structure, genetic diversity and presence of local adaptations, all dependent on dispersal within the metapopulation. Here, we surveyed 87 studies that combine population genetics and biophysical models of dispersal. Our aim was to understand if biophysical dispersal models can generally explain genetic differentiation. Our analysis shows that genetic differentiation and lack of genetic differentiation can often be explained by dispersal, but the realism of the biophysical model, as well as local geomorphology and species biology also play a role. The review supports the use of a combination of both methods, and we discuss our findings in terms of recommendations for future studies and pinpoint areas where further development is necessary, particularly on how to compare both approaches. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

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

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

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

BACKGROUND

As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis).

RESULTS

Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus.

CONCLUSION

The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

Trans-Atlantic Distribution and Introgression as Inferred from Single Nucleotide Polymorphism: Mussels Mytilus and Environmental Factors

Large-scale climate changes influence the geographic distribution of biodiversity. Many taxa have been reported to extend or reduce their geographic range, move poleward or displace other species. However, for closely related species that can hybridize in the natural environment, displacement is not the only effect of changes of environmental variables. Another option is subtler, hidden expansion, which can be found using genetic methods only. The marine blue mussels Mytilus are known to change their geographic distribution despite being sessile animals. In addition to natural dissemination at larval phase—enhanced by intentional or accidental introductions and rafting—they can spread through hybridization and introgression with local congeners, which can create mixed populations sustaining in environmental conditions that are marginal for pure taxa. The Mytilus species have a wide distribution in coastal regions of the Northern and Southern Hemisphere. In this study, we investigated the inter-regional genetic differentiation of the Mytilus species complex at 53 locations in the North Atlantic and adjacent Arctic waters and linked this genetic variability to key local environmental drivers. Of seventy-nine candidate single nucleotide polymorphisms (SNPs), all samples were successfully genotyped with a subset of 54 SNPs. There was a clear interregional separation of Mytilus species. However, all three Mytilus species hybridized in the contact area and created hybrid zones with mixed populations. Boosted regression trees (BRT) models showed that inter-regional variability was important in many allele models but did not prevail over variability in local environmental factors. Local environmental variables described over 40% of variability in about 30% of the allele frequencies of Mytilus spp. For the 30% of alleles, variability in their frequencies was only weakly coupled with local environmental conditions. For most studied alleles the linkages between environmental drivers and the genetic variability of Mytilus spp. were random in respect to “coding” and “non-coding” regions. An analysis of the subset of data involving functional genes only showed that two SNPs at Hsp70 and ATPase genes correlated with environmental variables. Total predictive ability of the highest performing models (r² between 0.550 and 0.801) were for alleles that discriminated most effectively M. trossulus from M. edulis and M. galloprovincialis, whereas the best performing allele model (BM101A) did the best at discriminating M. galloprovincialis from M. edulis and M. trossulus. Among the local environmental variables, salinity, water temperature, ice cover and chlorophyll a concentration were by far the greatest predictors, but their predictive performance varied among different allele models. In most cases changes in the allele frequencies along these environmental gradients were abrupt and occurred at a very narrow range of environmental variables. In general, regions of change in allele frequencies for M. trossulus occurred at 8–11 psu, 0–10 °C, 60%–70% of ice cover and 0–2 mg m⁻³ of chlorophyll a, M. edulis at 8–11 and 30–35 psu, 10–14 °C and 60%–70% of ice cover and for M. galloprovincialis at 30–35 psu, 14–20 °C.