East is East and West is West: Population genomics and hierarchical analyses reveal genetic structure and adaptation footprints in the keystone species Paracentrotus lividus (Echinoidea)

Genomic signals of adaptation to a natural CO2 gradient over a striking microgeographic scale

Our study explores genomic signs of adaptation in A. lixula to different water pH conditions. To achieve this, we analysed the genomics variation of A. lixula individuals living across a natural pH gradient in Canary Islands, Spain. We use a 2b-RADseq protocol with 74 samples from sites with varying pH levels (from 7.3 to 7.9 during low tide) and included a control site. We identified 14,883 SNPs, with 432 identified as candidate SNPs under selection to pH variations through redundancy analysis. While all SNPs indicated genomic homogeneity, the 432 candidate SNPs under selection displayed genomic differences among sites and along the pH gradient. Out of these 432 loci, 17 were annotated using published A. lixula transcriptomes, involved in biological functions such as growth. Therefore, our findings suggest local adaptation in A. lixula populations to acidification in CO2 vents, even over short distances of 75 m, underscoring their potential resistance to future Ocean Acidification.

Exploring reported population differences in Norway lobster (Nephrops norvegicus) in the Pomo Pits region of the Adriatic Sea using genome-wide markers

The Norway lobster (Nephrops norvegicus) is one of the most important decapod crustacean seafood species in the Adriatic Sea. Previous research has identified significant differences in growth rates and maturation timing of Nephrops in the Pomo/Jabuka Pits area compared to other subpopulations in Adriatic fishing grounds. Here, we use 1,623 genome-wide single nucleotide polymorphisms (SNPs) to investigate whether the Pomo Pits subpopulation is genetically different from other sites in the Adriatic and neighbouring seas. We found no genetic differentiation among all sampled Adriatic sites, suggesting high gene flow between Pomo Pits Nephrops and those of surrounding areas. We also found genetic homogeneity between the Adriatic sites and single-site samples from the Aegean and Tyrrhenian Seas. However, we detected distinct genetic differentiation between all Mediterranean sites and an Atlantic site in western Scotland, which provides evidence for a phylogenetic break between the Atlantic and the Mediterranean. Our results indicate that Pomo Pits Nephrops are not genetically different from others sampled in the Adriatic and that key biological parameters in Pomo Pits Nephrops could be driven by spatial variation in fishing pressure and/or environmental factors rather than geographic isolation.

Are Microfibers a Threat to Marine Invertebrates? A Sea Urchin Toxicity Assessment

The rise of "fast fashion" has driven up the production of low-cost, short-lived clothing, significantly increasing global textile fiber production and, consequently, exacerbating environmental pollution. This study investigated the ecotoxicological effects of different types of anthropogenic microfibers-cotton, polyester, and mixed fibers (50% cotton: 50% polyester)-on marine organisms, specifically sea urchin embryos. All tested fibers exhibited toxicity, with cotton fibers causing notable effects on embryonic development even at environmentally relevant concentrations. The research also simulated a scenario where microfibers were immersed in seawater for 30 days to assess changes in toxicity over time. The results showed that the toxicity of microfibers increased with both concentration and exposure duration, with polyester being the most toxic among the fibers tested. Although synthetic fibers have been the primary focus of previous research, this study highlights that natural fibers like cotton, which are often overlooked, can also be toxic due to the presence of harmful additives. These natural fibers, despite decomposing faster than synthetic ones, can persist in aquatic environments for extended periods. The findings underline the critical need for further research on both natural and synthetic microfibers to understand their environmental impact and potential threats to marine ecosystems and sea urchin populations.

Population genomic structure of the sea urchin Diadema africanum, a relevant species in the rocky reef systems across the Macaronesian archipelagos

The sea urchin Diadema africanum is a macro-herbivore found in the rocky reef systems of the West African region and Macaronesian archipelagos. Over several decades, high densities of this species have generated marine barrens in certain areas at the Canary Islands. In contrast, more recently, during the last few years, the species has suffered mass mortality events that continue to the present day. In this study, we used 9,109 Single Nucleotide Polymorphisms (SNPs) and a fragment of a mitochondrial gene to evaluate the species' population structure, effects of mass mortalities on its diversity, and potential local adaptation across the Canary Islands and Cabo Verde. Our research provides compelling evidence of low genomic diversity and homogeneity across the studied area for neutral markers, along with recent demographic fluctuations. The high connectivity among distant areas likely allows a rapid recovering of the populations from local mortality events. Interestingly, we also observed genomic sub-structure from 405 SNPs identified as candidate loci under selection for seawater temperature. The lack of divergence among distant sites and the low diversity found can be attributed to the species' divergence from a small ancestral genomic pool, followed by a contemporary demographic expansion, and ongoinggene flow.

Genetic homogeneity and weak signatures of local adaptation in the marine mussel Mytilus chilensis

The natural populations of the marine mussel Mytilus chilensis and the associated aquaculture industry forms a sensitive social-ecological system that relies on the released propagules for cultivation in the highly heterogeneous environment (temperature, productivity, and salinity) of northern Patagonia (42-44 °S). We assessed spatial genetic structure, signals of local adaptation, and population assignment of M. chilensis analyzing 5963 SNPs from 125 individuals across six natural populations sampled over two consecutive years along the southeast Pacific coast (39° 25' to 43° 07' S, ~ 430 km). Neutral and putatively adaptive loci revealed high genetic diversity and low genetic differentiation among populations. Of the whole dataset, less than 1% (50) of loci were identified as putatively adaptive through multiple approaches, with only 0.1% detected in by all of them, and only two loci of them were correlated with environmental variables. No evidence of Isolation by Environment (IBE) was found, albeit a slight differentiation in the southern sampling location (Yaldad). These results suggest that the genetic structure observed is primarily shaped by neutral processes with weak signals of local adaptation. Gene-flow appears to be the main evolutionary force influencing the species' population genetic structure. Because of the importance for the industry, the probability of correct assignment of individuals to their population of origin using allelic frequencies was evaluated. Analyses exhibited relatively low probabilities (< 50% for four out of six sites) of accurately assigning individuals to their geographic origin, with a limited success of SNP markers the for such purposes. Likely, species' high dispersal capacity, seed translocation, and the spill-over effect of mussel aquaculture prevents population genetic differentiation through high effective gene flow, hindering local genetic adaptation.

The past, the recent, and the ongoing evolutionary processes of the worldwide invasive ascidian Styela plicata

Invasive species are one of the main threats to global biodiversity and, within marine ecosystems, tunicates feature some prominent examples. Styela plicata is an ascidian species inhabiting harbours in all temperate oceans and seas, thus being considered a thriving invasive species. However, this species' adaptive mechanisms, introduction history, and population structure have never been completely elucidated. Here, by genotyping 87 S. plicata individuals from 18 localities worldwide with 2b-RADseq, we confirm the global presence of four chromosome inversions, demonstrate population structuring on this species, detect local adaptation signals, and infer historical demographic events. We show that North Carolina individuals constitute an unrelated population, Atlanto-Mediterranean and Pacific localities form their own genetic clusters with substructuring, being the most evident the split between northern and southern Atlantic localities. The locality of South Carolina presents an intermediate genetic position between North Carolina and the other two groups pointing to a hybrid origin with recurrent gene flow. We generate and test demographic models, providing evidence of two independent introduction events to the Atlantic and Pacific, and an admixture that originated the population of South Carolina. Finally, we identify candidate loci for adaptation, with functions involved with cell processes, metabolism, development, and ion transport, among others. Overall, this study highlights the complex historical processes of S. plicata, which have led this species to its current distribution, population structure, and local adaptation footprint in oceans worldwide.

Mitochondrial, nuclear and morphological differentiation in the swimming crab Liocarcinus depurator along the Atlantic-Mediterranean transition

Environmental gradients in the sea may coincide with phenotypic or genetic gradients resulting from an evolutionary balance between selection and dispersal. The population differentiation of the swimming crab, Liocarcinus depurator, an important by-catch species in the Mediterranean Sea and North-East Atlantic, was assessed using both genetic and morphometric approaches. A total of 472 specimens were collected along its distribution area, and 17 morphometric landmarks, one mitochondrial gene (COI) and 11 polymorphic microsatellite markers were scored in 350, 287 and 280 individuals, respectively. Morphometric data lacked significant differences, but genetic analyses showed significant genetic differentiation between Atlantic and Mediterranean populations, with a steeper gradient in COI compared to microsatellite markers. Interestingly, nuclear differentiation was due to an outlier locus with a gradient in the Atlantic-Mediterranean transition area overlapping with the mtDNA gradient. Such overlapping clines are likely to be maintained by natural selection. Our results suggest a scenario of past isolation with local adaptation and secondary contact between the two basins. Local adaptation during the process of vicariance may reinforce genetic differentiation at loci maintained by environmental selection even after secondary contact.

Evaluating restriction enzyme selection for reduced representation sequencing in conservation genomics

Conservation genomic studies in non-model organisms generally rely on reduced representation sequencing techniques based on restriction enzymes to identify population structure as well as candidate loci for local adaptation. While the expectation is that the reduced representation of the genome is randomly distributed, the proportion of the genome sampled might depend on the GC content of the recognition site of the restriction enzyme used. Here, we evaluated the distribution and functional composition of loci obtained after a reduced representation approach using Genotyping-by-Sequencing (GBS). To do so, we compared experimental data from two endemic fish species (Symphodus ocellatus and Symphodus tinca, EcoT22I enzyme) and two ecosystem engineer sea urchins (Paracentrotus lividus and Arbacia lixula, ApeKI enzyme). In brief, we mapped the sequenced loci to the phylogenetically closest reference genome available (Labrus bergylta in the fish and Strongylocentrotus purpuratus in the sea urchin datasets), classified them as exonic, intronic and intergenic, and studied their function by using Gene Ontology (GO) terms. We also simulated the effect of using both enzymes in the two reference genomes. In both simulated and experimental data, we detected an enrichment towards exonic or intergenic regions depending on the restriction enzyme used and failed to detect differences between total loci and candidate loci for adaptation in the empirical dataset. Most of the functions assigned to the mapped loci were shared between the four species and involved a myriad of general functions. Our results highlight the importance of restriction enzyme selection and the need for high-quality annotated genomes in conservation genomic studies.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Occurrence of microfibres in wild specimens of adult sea urchin Paracentrotus lividus (Lamarck, 1816) from a coastal area of the central Mediterranean Sea

This study investigates the occurrence of anthropogenic fibres inside wild Paracentrotus lividus at a Mediterranean coastal area in 2020. From each sea urchin, the coelomic fluid was directly analysed while digestive tracts and gonads were removed, pre-treated with trypsin (0.3%) and digested with H₂O₂ (10%) before analysis. A total of 260 fibres and 1 fragment were found in 100 specimens, with an average of 2.6 items/individual. Fibres were more abundant in the digestive system, less in gonads and in the coelomic fluid, respectively. Fourier transform infrared (FTIR) analysis of representative fibres identified 67% natural (cotton-based) and 33% synthetic polymers (polyester) suggesting their origin from textiles, possibly released from laundry sewages. Overall, these results encourage further in-depth investigations on fibres accumulation and potential transfer through the trophic chain up to humans.

Climate differently influences the genomic patterns of two sympatric marine fish species

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Genome-Wide SNP Discovery and Population Genetic Analysis of Mesocentrotus nudus in China Seas

Mesocentrotus nudus is an important commercially aquatic species because of its high edible and medicinal values. However, wild stocks have dramatically decreased in recent decades. Understanding the population structure and genetic diversity can provide vital information for genetic conservation and improvement. In the present study, the genotyping-by-sequencing (GBS) approach was adopted to identify the genome-wide single-nucleotide polymorphisms (SNPs) from a collection of 80 individuals consisting of five geographical populations (16 individuals from each population), covering the natural habitats of M. nudus in China seas. An average of 0.96-Gb clean reads per sample were sequenced, and a total of 51,738 biallelic SNPs were identified. Based on these SNPs, diversity index analysis showed that all populations have a similar pattern with positive F_is (0.136) and low Ne (724.3). Low genetic differentiation and high genetic connectivity among five geographical populations were detected by pairwise F_st, principal component analysis (PCA), admixture, and phylogenetic analysis. Besides, two YWL individuals originating from an isolated ancestor may imply that there is a genetically differentiated population in the adjacent sea. Overall, the results showed that GBS is an effective method to detect genome-wide SNPs for M. nudus and suggested that the protective measures and the investigation with larger spatial scale and sample size for M. nudus should be carried out in the future.

Genetic homogeneity, lack of larvae recruitment, and clonality in absence of females across western Mediterranean populations of the starfish Coscinasterias tenuispina

We here analysed the populations' genetic structure of Coscinasterias tenuispina, an Atlantic-Mediterranean fissiparous starfish, focusing on the western Mediterranean, to investigate: the distribution and prevalence of genetic variants, the relative importance of asexual reproduction, connectivity across the Atlantic-Mediterranean transition, and the potential recent colonisation of the Mediterranean Sea. Individuals from 11 Atlantic-Mediterranean populations of a previous study added to 172 new samples from five new W Mediterranean sites. Individuals were genotyped at 12 microsatellite loci and their gonads histologically analysed for sex determination. Additionally, four populations were genotyped at two-time points. Results demonstrated genetic homogeneity and low clonal richness within the W Mediterranean, due to the dominance of a superclone, but large genetic divergence with adjacent areas. The lack of new genotypes recruitment over time, and the absence of females, confirmed that W Mediterranean populations were exclusively maintained by fission and reinforced the idea of its recent colonization. The existence of different environmental conditions among basins and/or density-depend processes could explain this lack of recruitment from distant areas. The positive correlation between clonal richness and heterozygote excess suggests that most genetic diversity is retained within individuals in the form of heterozygosity in clonal populations, which might increase their resilience.

Genetic connectivity of two marine gastropods in the Mediterranean Sea: seascape genetics reveals species-specific oceanographic drivers of gene flow

Oceanographic features such as currents, waves, temperature and salinity, together with life history traits, control patterns and rates of gene flow and contribute to shaping the population genetic structure of marine organisms. Seascape genetics is an emerging discipline that adopts a spatially explicit approach to examine biotic and abiotic factors that drive gene flow in marine environments. In this study, we examined factors that contribute to genetic differentiation in two coastal Mediterranean gastropods whose geographical ranges overlap but which inhabit different environments. The two species differ in several life history traits and in their dispersal capabilities. Genetic differentiation was relatively low for the trochid species Gibbula divaricata (F_ST  =0.059), and high for the vermetid species Dendropoma lebeche (F_ST  =0.410). Salinity emerged as the most important variable explaining the genetic structure of both species; sea surface temperature was also important for G. divaricata. For the more sessile D. lebeche, the coastline was predicted to provide important pathways for stepping-stone connectivity and gene flow. Our results provide a greater understanding of the factors influencing marine population connectivity, which may be useful to guide marine conservation and management in the Mediterranean.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

Genetic structure and life history are key factors in species distribution models of spiny lobsters

AIM

We incorporated genetic structure and life history phase in species distribution models (SDMs) constructed for a widespread spiny lobster, to reveal local adaptations specific to individual subspecies and predict future range shifts under the RCP 8.5 climate change scenario.

LOCATION

Indo-West Pacific.

METHODS

MaxEnt was used to construct present-day SDMs for the spiny lobster Panulirus homarus and individually for the three genetically distinct subspecies of which it comprises. SDMs incorporated both sea surface and benthic (seafloor) climate layers to recreate discrete influences of these habitats during the drifting larval and benthic juvenile and adult life history phases. Principle component analysis (PCA) was used to infer environmental variables to which individual subspecies were adapted. SDM projections of present-day habitat suitability were compared with predictions for the year 2,100, under the RCP 8.5 climate change scenario.

RESULTS

In the PCA, salinity best explained P. h. megasculptus habitat suitability, compared with current velocity in P. h. rubellus and sea surface temperature in P. h. homarus. Drifting and benthic life history phases were adapted to different combinations of sea surface and benthic environmental variables considered. Highly suitable habitats for benthic phases were spatially enveloped within more extensive sea surface habitats suitable for drifting larvae. SDMs predicted that present-day highly suitable habitats for P. homarus will decrease by the year 2,100.

MAIN CONCLUSIONS

Incorporating genetic structure in SDMs showed that individual spiny lobster subspecies had unique adaptations, which could not be resolved in species-level models. The use of sea surface and benthic climate layers revealed the relative importance of environmental variables during drifting and benthic life history phases. SDMs that included genetic structure and life history were more informative in predictive models of climate change effects.

Distribution modelling of an introduced species: do adaptive genetic markers affect potential range?

Biological invasions have increased in the last few decades mostly due to anthropogenic causes such as globalization of trade. Because invaders sometimes cause large economic losses and ecological disturbances, estimating their origin and potential geographical ranges is useful. Drosophila subobscura is native to the Old World but was introduced in the New World in the late 1970s and spread widely. We incorporate information on adaptive genetic markers into ecological niche modelling and then estimate the most probable geographical source of colonizers; evaluate whether the genetic bottleneck experienced by founders affects their potential distribution; and finally test whether this species has spread to all its potential suitable habitats worldwide. We find the environmental space occupied by this species in its native and introduced distributions are notably the same, although the introduced niche has shifted slightly towards higher temperature and lower precipitation. The genetic bottleneck of founding individuals was a key factor limiting the spread of this introduced species. We also find that regions in the Mediterranean and north-central Portugal show the highest probability of being the origin of the colonizers. Using genetically informed environmental niche modelling can enhance our understanding of the initial colonization and spread of invasive species, and also elucidate potential areas of future expansions worldwide.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Individual-based population genomics reveal different drivers of adaptation in sympatric fish

Connectivity and local adaptation are two contrasting evolutionary forces highly influencing population structure. To evaluate the impact of early-life traits and environmental conditions on genetic structuring and adaptation, we studied two sympatric fish species in the Western Mediterranean Sea: Symphodus tinca and S. ocellatus. We followed an individual-based approach and measured early-life history traits from otolith readings, gathered information on environmental variables and obtained genome-wide markers from genotyping-by-sequencing (GBS). The two species presented contrasting population structure across the same geographic gradient, with high and significant population differentiation in S. ocellatus, mostly determined by oceanographic fronts, and low differentiation and no front effect in S. tinca. Despite their different levels of genetic differentiation, we identified in both species candidate regions for local adaptation by combining outlier analysis with environmental and phenotypic association analyses. Most candidate loci were associated to temperature and productivity in S. ocellatus and to temperature and turbulence in S. tinca suggesting that different drivers may determine genomic diversity and differentiation in each species. Globally, our study highlights that individual-based approach combining genomic, environmental and phenotypic information is key to identify signals of selection and the processes mediating them.