A Pleistocene legacy structures variation in modern seagrass ecosystems.

Distribution of Earth's biomes is structured by the match between climate and plant traits, which in turn shape associated communities and ecosystem processes and services. However, that climate-trait match can be disrupted by historical events, with lasting ecosystem impacts. As Earth's environment changes faster than at any time in human history, critical questions are whether and how organismal traits and ecosystems can adjust to altered conditions. We quantified the relative importance of current environmental forcing versus evolutionary history in shaping the growth form (stature and biomass) and associated community of eelgrass (Zostera marina), a widespread foundation plant of marine ecosystems along Northern Hemisphere coastlines, which experienced major shifts in distribution and genetic composition during the Pleistocene. We found that eelgrass stature and biomass retain a legacy of the Pleistocene colonization of the Atlantic from the ancestral Pacific range and of more recent within-basin bottlenecks and genetic differentiation. This evolutionary legacy in turn influences the biomass of associated algae and invertebrates that fuel coastal food webs, with effects comparable to or stronger than effects of current environmental forcing. Such historical lags in phenotypic acclimatization may constrain ecosystem adjustments to rapid anthropogenic climate change, thus altering predictions about the future functioning of ecosystems.

Comparative landscape genomics has arrived with a splash.

Comparative methods are central for understanding the distribution of biodiversity. Assessing spatial variation of multiple species can identify biodiversity drivers across landscapes, including on the genetic level. Molecular ecology approaches have expanded in spatial precision, from phylogeography to landscape genetics, as have molecular tools that now allow genome-scale inferences for virtually any organism. These welcome increases in genomic data richness, however, have not been matched by growth in the multispecies dimension, and empirical studies continue to focus predominantly on single species. In this issue of Molecular Ecology, Zbinden et al. (Molecular Ecology, 2022, 32) present a landmark achievement for comparative landscape genomics, surveying an impressive 31 species of fishes, sampled from 75 locations in the White River Basin and genotyping each species for thousands of SNPs. Zbinden et al. (Molecular Ecology, 2022, 32) draw upon the statistical tools of landscape genetics to comprehensively interrogate four hypotheses-that populations of multiple freshwater fish species are isolated by river distance, barriers to dispersal, stream hierarchy or environment. Stream hierarchy overwhelmingly predicts population structure within co-distributed fish species pointing to commonalities among species at the large landscape scale (100s of km). The tight alignment between intraspecific genetic spatial structure and stream hierarchies, moreover, provides clear validation for conservation and fisheries management to use watershed divisions as distinct management units.

Seascape genomics identify adaptive barriers correlated to tidal amplitude in the shore crab Carcinus maenas.

Most marine invertebrates disperse during a planktonic larval stage that may drift for weeks with ocean currents. A challenge for larvae of coastal species is to return to coastal nursery habitats. Shore crab (Carcinus maenas L.) larvae are known to show tidal rhythmicity in vertical migration in tidal areas and circadian rhythmicity in microtidal areas, which seems to increase successful coastal settlement. We studied genome-wide differentiation based on 24,000 single nucleotide polymorphisms of 12 native populations of shore crab sampled from a large tidal amplitude gradient from macrotidal (~8 m) to microtidal (~0.2 m). Dispersal and recruitment success of larvae was assessed with a Lagrangian biophysical model, which showed a strong effect of larval behaviour on long-term connectivity, and dispersal barriers that partly coincided with different tidal environments. The genetic population structure showed a subdivision of the samples into three clusters, which represent micro-, meso- and macrotidal areas. The genetic differentiation was mostly driven by 0.5% outlier loci, which showed strong allelic clines located at the limits between the three tidal areas. Demographic modelling suggested that the two genetic barriers have different origins. Differential gene expression of two clock genes (cyc and pdp1) further highlighted phenotypic differences among genetic clusters that are potentially linked to the differences in larval behaviour. Taken together, our seascape genomic study suggests that tidal regime acts as a strong selection force on shore crab population structure, consistent with larval behaviour affecting dispersal and recruitment success.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species.

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double-digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long-term (~30-year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation-by-distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Integrating genetics, biophysical, and demographic insights identifies critical sites for seagrass conservation.

The eelgrass Zostera marina is an important foundation species of coastal areas in the Northern Hemisphere, but is continuing to decline, despite management actions. The development of new management tools is therefore urgent in order to prioritize limited resources for protecting meadows most vulnerable to local extinctions and identifying most valuable present and historic meadows to protect and restore, respectively. We assessed 377 eelgrass meadows along the complex coastlines of two fjord regions on the Swedish west coast-one is currently healthy and the other is substantially degraded. Shoot dispersal for all meadows was assessed with Lagrangian biophysical modeling (scale: 100-1,000 m) and used for barrier analysis and clustering; a subset (n = 22) was also assessed with population genetic methods (20 microsatellites) including diversity, structure, and network connectivity. Both approaches were in very good agreement, resulting in seven subpopulation groupings or management units (MUs). The MUs correspond to a spatial scale appropriate for coastal management of "waterbodies" used in the European Water Framework Directive. Adding demographic modeling based on the genetic and biophysical data as a third approach, we are able to assess past, present, and future metapopulation dynamics to identify especially vulnerable and valuable meadows. In a further application, we show how the biophysical approach, using eigenvalue perturbation theory (EPT) and distribution records from the 1980s, can be used to identify lost meadows where restoration would best benefit the present metapopulation. The combination of methods, presented here as a toolbox, allows the assessment of different temporal and spatial scales at the same time, as well as ranking of specific meadows according to key genetic, demographic and ecological metrics. It could be applied to any species or region, and we exemplify its versatility as a management guide for eelgrass along the Swedish west coast.

m6A RNA Methylation in Marine Plants: First Insights and Relevance for Biological Rhythms.

Circadian regulations are essential for enabling organisms to synchronize physiology with environmental light-dark cycles. Post-transcriptional RNA modifications still represent an understudied level of gene expression regulation in plants, although they could play crucial roles in environmental adaptation. N6-methyl-adenosine (m6A) is the most prevalent mRNA modification, established by "writer" and "eraser" proteins. It influences the clockwork in several taxa, but only few studies have been conducted in plants and none in marine plants. Here, we provided a first inventory of m6A-related genes in seagrasses and investigated daily changes in the global RNA methylation and transcript levels of writers and erasers in Cymodocea nodosa and Zostera marina. Both species showed methylation peaks during the dark period under the same photoperiod, despite exhibiting asynchronous changes in the m6A profile and related gene expression during a 24-h cycle. At contrasting latitudes, Z. marina populations displayed overlapping daily patterns of the m6A level and related gene expression. The observed rhythms are characteristic for each species and similar in populations of the same species with different photoperiods, suggesting the existence of an endogenous circadian control. Globally, our results indicate that m6A RNA methylation could widely contribute to circadian regulation in seagrasses, potentially affecting the photo-biological behaviour of these plants.

Adaptive responses along a depth and a latitudinal gradient in the endemic seagrass Posidonia oceanica.

Seagrass meadows provide important ecosystem services and are critical for the survival of the associated invertebrate community. However, they are threatened worldwide by human-driven environmental change. Understanding the seagrasses' potential for adaptation is critical to assess not only their ability to persist under future global change scenarios, but also to assess the persistence of the associated communities. Here we screened a wild population of Posidonia oceanica, an endemic long-lived seagrass in the Mediterranean Sea, for genes that may be target of environmental selection, using an outlier and a genome-wide transcriptome analysis. We identified loci where polymorphism or differential expression was associated with either a latitudinal or a bathymetric gradient, as well as with both gradients in an effort to identify loci associated with temperature and light. We found the candidate genes underlying growth and immunity to be divergent between populations adapted to different latitudes and/or depths, providing evidence for local adaptation. Furthermore, we found evidence of reduced gene flow among populations including adjacent populations. Reduced gene flow, combined with low sexual recombination, small effective population size, and long generation time of P. oceanica raises concerns for the long-term persistence of this species, especially in the face of rapid environmental change driven by human activities.

A meta-analysis reveals a positive correlation between genetic diversity metrics and environmental status in the long-lived seagrass Posidonia oceanica.

The seagrass Posidonia oceanica is a key engineering species structuring coastal marine systems throughout much of the Mediterranean basin. Its decline is of concern, leading to the search for short- and long-term indicators of seagrass health. Using ArcGIS maps from a recent, high-resolution (1-4 km) modelling study of 18 disturbance factors affecting coastal marine systems across the Mediterranean (Micheli et al. 2013, http://globalmarine.nceas.ucsb.edu/mediterranean/), we tested for correlations with genetic diversity metrics (allelic diversity, genotypic/clonal diversity and heterozygosity) in a meta-analysis of 56 meadows. Contrary to initial predictions, weak but significantly positive correlations were found for commercial shipping, organic pollution (pesticides) and cumulative impact. This counterintuitive finding suggests greater resistance and resilience of individuals with higher genetic and genotypic diversity under disturbance (at least for a time) and/or increased sexual reproduction under an intermediate disturbance model. We interpret the absence of low and medium levels of genetic variation at impacted locations as probable local extinctions of individuals that already exceeded their resistance capacity. Alternatively, high diversity at high-impact sites is likely a temporal artefact, reflecting the mismatch with pre-environmental impact conditions, especially because flowering and sexual recruitment are seldom observed. While genetic diversity metrics are a valuable tool for restoration and mitigation, caution must be exercised in the interpretation of correlative patterns as found in this study, because the exceptional longevity of individuals creates a temporal mismatch that may falsely suggest good meadow health status, while gradual deterioration of allelic diversity might go unnoticed.

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.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina).

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

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

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.

Evolution and phylogeography of the nonpathogenic calicivirus RCV-A1 in wild rabbits in Australia.

Despite its potential importance for the biological control of European rabbits, relatively little is known about the evolution and molecular epidemiology of rabbit calicivirus Australia 1 (RCV-A1). To address this issue we undertook an extensive evolutionary analysis of 36 RCV-A1 samples collected from wild rabbit populations in southeast Australia between 2007 and 2009. Based on phylogenetic analysis of the entire capsid sequence, six clades of RCV-A1 were defined, each exhibiting strong population subdivision. Strikingly, our estimates of the time to the most recent common ancestor of RCV-A1 coincide with the introduction of rabbits to Australia in the mid-19th century. Subsequent divergence events visible in the RCV-A1 phylogenies likely reflect key moments in the history of the European rabbit in Australia, most notably the bottlenecks in rabbit populations induced by the two viral biocontrol agents used on the Australian continent, myxoma virus and rabbit hemorrhagic disease virus (RHDV). RCV-A1 strains therefore exhibit strong phylogeographic separation and may constitute a useful tool to study recent host population dynamics and migration patterns, which in turn could be used to monitor rabbit control in Australia.

Population genetic structure and connectivity of the seagrass Thalassia hemprichii in the Western Indian Ocean is influenced by predominant ocean currents.

This study is the first large-scale genetic population study of a widespread climax species of seagrass, Thalassia hemprichii, in the Western Indian Ocean (WIO). The aim was to understand genetic population structure and connectivity of T. hemprichii in relation to hydrodynamic features. We genotyped 205 individual seagrass shoots from 11 sites across the WIO, spanning over a distance of ~2,700 km, with twelve microsatellite markers. Seagrass shoots were sampled in Kenya, Tanzania (mainland and Zanzibar), Mozambique, and Madagascar: 4-26°S and 33-48°E. We assessed clonality and visualized genetic diversity and genetic population differentiation. We used Bayesian clustering approaches (TESS) to trace spatial ancestry of populations and used directional migration rates (DivMigrate) to identify sources of gene flow. We identified four genetically differentiated groups: (a) samples from the Zanzibar channel; (b) Mozambique; (c) Madagascar; and (d) the east coast of Zanzibar and Kenya. Significant pairwise population genetic differentiation was found among many sites. Isolation by distance was detected for the estimated magnitude of divergence (D EST), but the three predominant ocean current systems (i.e., East African Coastal Current, North East Madagascar Current, and the South Equatorial Current) also determine genetic connectivity and genetic structure. Directional migration rates indicate that Madagascar acts as an important source population. Overall, clonality was moderate to high with large differences among sampling sites, indicating relatively low, but spatially variable sexual reproduction rates. The strongest genetic break was identified for three sites in the Zanzibar channel. Although isolation by distance is present, this study suggests that the three regionally predominant ocean current systems (i.e., East African Coastal Current, North East Madagascar Current, and the South Equatorial Current) rather than distance determine genetic connectivity and structure of T. hemprichii in the WIO. If the goal is to maintain genetic connectivity of T. hemprichii within the WIO, conservation planning and implementation of marine protection should be considered at the regional scale-across national borders.

Seascape genetics and biophysical connectivity modelling support conservation of the seagrass Zostera marina in the Skagerrak-Kattegat region of the eastern North Sea.

Maintaining and enabling evolutionary processes within meta-populations are critical to resistance, resilience and adaptive potential. Knowledge about which populations act as sources or sinks, and the direction of gene flow, can help to focus conservation efforts more effectively and forecast how populations might respond to future anthropogenic and environmental pressures. As a foundation species and habitat provider, Zostera marina (eelgrass) is of critical importance to ecosystem functions including fisheries. Here, we estimate connectivity of Z. marina in the Skagerrak-Kattegat region of the North Sea based on genetic and biophysical modelling. Genetic diversity, population structure and migration were analysed at 23 locations using 20 microsatellite loci and a suite of analytical approaches. Oceanographic connectivity was analysed using Lagrangian dispersal simulations based on contemporary and historical distribution data dating back to the late 19th century. Population clusters, barriers and networks of connectivity were found to be very similar based on either genetic or oceanographic analyses. A single-generation model of dispersal was not realistic, whereas multigeneration models that integrate stepping-stone dispersal and extant and historic distribution data were able to capture and model genetic connectivity patterns well. Passive rafting of flowering shoots along oceanographic currents is the main driver of gene flow at this spatial-temporal scale, and extant genetic connectivity strongly reflects the "ghost of dispersal past" sensu Benzie, 1999. The identification of distinct clusters, connectivity hotspots and areas where connectivity has become limited over the last century is critical information for spatial management, conservation and restoration of eelgrass.

Tissue-specific transcriptomic profiling provides new insights into the reproductive ecology and biology of the iconic seagrass species Posidonia oceanica.

Seagrasses form extensive meadows in shallow coastal waters and are among the world's most productive ecosystems. Seagrasses can produce both clonally and sexually, and flowering has long been considered infrequent, but important for maintaining genetically diverse stands. Here we investigate the molecular mechanisms involved in flowering of the seagrass Posidonia oceanica, an iconic species endemic to the Mediterranean. We generated a de novo transcriptome of this non-model species for leaf, male and female flower tissue of three individuals, and present molecular evidence for genes that may be involved in the flowering process and on the reproductive biology of the species. We present evidence that suggests that P. oceanica exhibits a strategy of protogyny, where the female part of the hermaphroditic flower develops before the male part, in order to avoid self-fertilization. We found photosynthetic genes to be up-regulated in the female flower tissues, indicating that this may be capable of photosynthesis. Finally, we detected a number of interesting genes, previously known to be involved in flowering pathways responding to light and temperature cues and in pathways involved in anthocyanin and exine synthesis. This first comparative transcriptomic approach of leaf, male and female tissue provides a basis for functional genomics research on flower development in P. oceanica and other seagrass species.

Patterns and mechanisms of dispersal in a keystone seagrass species.

Mechanisms and vectors of long-distance dispersal remain unknown for many coastal benthic species, including plants. Indications for the possibility for long-distance dispersal come from dispersal modelling and from genetic assessments, but have rarely been assessed with both methods. To this end, we assessed dispersal of the seagrass Zostera noltei, an important foundation species of the coastal zone. We investigate whether small scale seed dispersal and long-distance propagule dispersal do play a role for meta-population dynamics, using both genetic assessments based on eight microsatellite markers and physical modelling of ocean currents. Such assessments enhance our understanding of the biology and population dynamics of an important coastal foundation species. They are relevant for large scale conservation strategies as they give insights in the maintenance of genetic diversity and connectivity that may enhance resilience and resistance to stresses associated with seagrass loss.

Adaptation and acclimatization to ocean acidification in marine ectotherms: an in situ transplant experiment with polychaetes at a shallow CO2 vent system.

Metabolic rate determines the physiological and life-history performances of ectotherms. Thus, the extent to which such rates are sensitive and plastic to environmental perturbation is central to an organism's ability to function in a changing environment. Little is known of long-term metabolic plasticity and potential for metabolic adaptation in marine ectotherms exposed to elevated pCO2. Consequently, we carried out a series of in situ transplant experiments using a number of tolerant and sensitive polychaete species living around a natural CO2 vent system. Here, we show that a marine metazoan (i.e. Platynereis dumerilii) was able to adapt to chronic and elevated levels of pCO2. The vent population of P. dumerilii was physiologically and genetically different from nearby populations that experience low pCO2, as well as smaller in body size. By contrast, different populations of Amphiglena mediterranea showed marked physiological plasticity indicating that adaptation or acclimatization are both viable strategies for the successful colonization of elevated pCO2 environments. In addition, sensitive species showed either a reduced or increased metabolism when exposed acutely to elevated pCO2. Our findings may help explain, from a metabolic perspective, the occurrence of past mass extinction, as well as shed light on alternative pathways of resilience in species facing ongoing ocean acidification.