Local environment-driven adaptive evolution in a marine invasive ascidian (Molgula manhattensis)

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Roles of eyestalk in salinity acclimatization of mud crab (Scylla paramamosain) by transcriptomic analysis

Salinity acclimatization refers to the physiological and behavioral adjustments made by crustaceans to adapt to varying salinity environments. The eyestalk, a neuroendocrine organ in crustaceans, plays a crucial role in salinity acclimatization. To elucidate the molecular mechanisms underlying eyestalk involvement in mud crab (Scylla paramamosain) acclimatization, we employed RNA-seq technology to analyze transcriptomic changes in the eyestalk under low (5 ppt) and standard (23 ppt) salinity conditions. This analysis revealed 5431 differentially expressed genes (DEGs), with 2372 upregulated and 3059 downregulated. Notably, these DEGs were enriched in crucial biological pathways like metabolism, osmoregulation, and signal transduction. To validate the RNA-seq data, we further analyzed 15 DEGs of interest using qRT-PCR. Our results suggest a multifaceted role for the eyestalk: maintaining energy homeostasis, regulating hormone synthesis and release, PKA activity, and downstream signaling, and ensuring proper ion and osmotic balance. Furthermore, our findings indicate that the crustacean hyperglycemic hormone (CHH) may function as a key regulator, modulating carbonic anhydrase expression through the activation of the PKA signaling pathway, thereby influencing cellular osmoregulation, and associated metabolic processes. Overall, our study provides valuable insights into unraveling the molecular mechanisms of mud crab acclimatization to low salinity environments.

Incorporating adaptive genomic variation into predictive models for invasion risk assessment

Global climate change is expected to accelerate biological invasions, necessitating accurate risk forecasting and management strategies. However, current invasion risk assessments often overlook adaptive genomic variation, which plays a significant role in the persistence and expansion of invasive populations. Here we used Molgula manhattensis, a highly invasive ascidian, as a model to assess its invasion risks along Chinese coasts under climate change. Through population genomics analyses, we identified two genetic clusters, the north and south clusters, based on geographic distributions. To predict invasion risks, we employed the gradient forest and species distribution models to calculate genomic offset and species habitat suitability, respectively. These approaches yielded distinct predictions: the gradient forest model suggested a greater genomic offset to future climatic conditions for the north cluster (i.e., lower invasion risks), while the species distribution model indicated higher future habitat suitability for the same cluster (i.e, higher invasion risks). By integrating these models, we found that the south cluster exhibited minor genome-niche disruptions in the future, indicating higher invasion risks. Our study highlights the complementary roles of genomic offset and habitat suitability in assessing invasion risks under climate change. Moreover, incorporating adaptive genomic variation into predictive models can significantly enhance future invasion risk predictions and enable effective management strategies for biological invasions in the future.

(Epi)genomic adaptation driven by fine geographical scale environmental heterogeneity after recent biological invasions

Elucidating processes and mechanisms involved in rapid local adaptation to varied environments is a poorly understood but crucial component in management of invasive species. Recent studies have proposed that genetic and epigenetic variation could both contribute to ecological adaptation, yet it remains unclear on the interplay between these two components underpinning rapid adaptation in wild animal populations. To assess their respective contributions to local adaptation, we explored epigenomic and genomic responses to environmental heterogeneity in eight recently colonized ascidian (Ciona intestinalis) populations at a relatively fine geographical scale. Based on MethylRADseq data, we detected strong patterns of local environment-driven DNA methylation divergence among populations, significant epigenetic isolation by environment (IBE), and a large number of local environment-associated epigenetic loci. Meanwhile, multiple genetic analyses based on single nucleotide polymorphisms (SNPs) showed genomic footprints of divergent selection. In addition, for five genetically similar populations, we detected significant methylation divergence and local environment-driven methylation patterns, indicating the strong effects of local environments on epigenetic variation. From a functional perspective, a majority of functional genes, Gene Ontology (GO) terms, and biological pathways were largely specific to one of these two types of variation, suggesting partial independence between epigenetic and genetic adaptation. The methylation quantitative trait loci (mQTL) analysis showed that the genetic variation explained only 18.67% of methylation variation, further confirming the autonomous relationship between these two types of variation. Altogether, we highlight the complementary interplay of genetic and epigenetic variation involved in local adaptation, which may jointly promote populations' rapid adaptive capacity and successful invasions in different environments. The findings here provide valuable insights into interactions between invaders and local environments to allow invasive species to rapidly spread, thus contributing to better prediction of invasion success and development of management strategies.

Invader at the edge - Genomic origins and physiological differences of round gobies across a steep urban salinity gradient

Species invasions are a global problem of increasing concern, especially in highly connected aquatic environments. Despite this, salinity conditions can pose physiological barriers to their spread, and understanding them is important for management. In Scandinavia's largest cargo port, the invasive round goby (Neogobius melanostomus) is established across a steep salinity gradient. We used 12,937 SNPs to identify the genetic origin and diversity of three sites along the salinity gradient and round goby from western, central and northern Baltic Sea, as well as north European rivers. Fish from two sites from the extreme ends of the gradient were also acclimated to freshwater and seawater, and tested for respiratory and osmoregulatory physiology. Fish from the high-salinity environment in the outer port showed higher genetic diversity, and closer relatedness to the other regions, compared to fish from lower salinity upstream the river. Fish from the high-salinity site also had higher maximum metabolic rate, fewer blood cells and lower blood Ca²⁺. Despite these genotypic and phenotypic differences, salinity acclimation affected fish from both sites in the same way: seawater increased the blood osmolality and Na⁺ levels, and freshwater increased the levels of the stress hormone cortisol. Our results show genotypic and phenotypic differences over short spatial scales across this steep salinity gradient. These patterns of the physiologically robust round goby are likely driven by multiple introductions into the high-salinity site, and a process of sorting, likely based on behaviour or selection, along the gradient. This euryhaline fish risks spreading from this area, and seascape genomics and phenotypic characterization can inform management strategies even within an area as small as a coastal harbour inlet.

Evidence from transcriptome analysis unravelled the roles of eyestalk in salinity adaptation in Pacific white shrimp (Litopenaeus vannamei)

Eyestalk is considered the main neuroendocrine organ in crustaceans. Eyestalk regulates reproduction, molting, and energy metabolism by secreting several neurohormones. However, the role of eyestalk in salinity adaptation in crustaceans remains unclear. To reveal the role of eyestalk in salinity adaptation in Litopenaeus vannamei, we performed RNA-seq to compare the transcriptomic response of the eyestalk under low salinity (salinity 3) with that of the control group (salinity 25) for 8 weeks. A total of 479 mRNAs, including 150 upregulated and 329 downregulated mRNAs, were differentially expressed between the two salinity groups. The majority of the differentially expressed genes (DEGs) were enriched in biological pathways related to osmoregulation, metabolism and energy production, and oxidative stress. The most important DEGs associated with osmoregulation were CA4, ATP1A, ATP2B, ABCB1, ABCC4, PhoA, PhoB, NOS1, ACE, ANPEP, and the V-type H⁺-ATPase E-subunit. The metabolism-related DEGs were divided into three main categories: carbohydrate and energy metabolism (i.e., G6PC, UGT), protein and amino acid metabolism (i.e., SLC15A1, AhcY, GFAT), and lipid and fatty acid metabolism (i.e., GPAT3_4, CYP2J). The key DEGs related to the oxidative stress response were UGT, NDUFB1, QCR7, QCR8, P5CDh, COX6B, and CES1. These results provide evidence for the existence of an eyestalk-salinity adaptation-stress endocrine axis in L. vannamei. These findings provide a better understanding of the molecular mechanism underlying salinity adaptation in L. vannamei.

Hypo-Osmoregulatory Roles of Vasotocinergic and Isotocinergic Systems in the Intestines of Two European Sea Bass Lineages

European sea bass (Dicentrarchus labrax) are a major aquaculture species that live in habitats with fluctuating salinities that are sometimes higher than in seawater (SW). Atlantic and West-Mediterranean genetic lineages were compared regarding intestinal neuropeptide receptor expression in SW (36%) and following a two-week transfer to hypersalinity (HW, 55%). Phylogenetic analysis revealed seven neuropeptide receptors belonging to the arginine vasotocine (AVTR) family and two isotocin receptors (ITR). Among AVTR paralogs, the highest mRNA levels were recorded for v1a2, with a two- to fourfold upregulation in the European sea bass intestinal sections after transfer of fish to HW. Principal component analysis in posterior intestines showed that v1a2 expression grouped together with the expression and activity of main ion transporters and channels involved in solute-coupled water uptake, indicating a possible role of this receptor in triggering water absorption. v1a1 expression, however, was decreased or did not change after transfer to hypersaline water. Among ITR paralogs, itr1 was the most expressed paralog in the intestine and opposite expression patterns were observed following salinity transfer, comparing intestinal sections. Overall, different expression profiles were observed between genetic lineages for several analyzed genes which could contribute to different osmotic stress-related responses in D. labrax lineages.

Shipping voyage simulation reveals abiotic barriers to marine bioinvasions

The shipping industry is considered the main vector of introduction of marine non-indigenous species (NIS). NIS distributions are often a consequence of frequent trade activities that are affected by economic trends. A dominant trend in the shipping industry is the operation of Ultra Large Container Vessels (ULCV), which are over 395 m long and sail mostly on the East-Asia - northern-Europe route. Understanding the risk of NIS introduction by this emerging shipping category is needed for devising strategies for sustainable shipping. Here, we conducted a controlled simulation of key abiotic factors that determine marine bioinvasion success: temperature, salinity, and food availability along selected routes, under two treatments: ULCV and intermediate-size vessels. We tested the effect of each treatment and the varying environmental conditions on the survival of two invasive ascidians (Chordata, Ascidiacea). We used survival analysis methods to locate predictors of ascidian mortality; Environmental conditions at ports with high mortality were used to identify similar major ports on a global scale as potential abiotic barriers. The key factors in ascidian mortality varied between the two species, but for both species, the treatment and salinity were dominant predictors for survival. We identified Port Klang, Rotterdam, and Dammam as ports with high mortality and located several globally distributed major ports that present similar environmental conditions. Our results highlight the potential role of selected major ports as abiotic barriers to fouling organisms during ocean voyages. The tolerance of the tropical-origin Microcosmus exasperatus to the northern-Europe conditions, and of the temperate/sub-tropical origin Styela plicata, to high temperature conditions, point out the urgent need to modify international fouling regulations in view of global change. Further studies on the survival of fouling organisms during a cascade of changing environmental conditions will contribute to the advancement of science-based regulations to reduce the adverse effects of NIS.

Complementary genomic and epigenomic adaptation to environmental heterogeneity

While adaptation is commonly thought to result from selection on DNA sequence-based variation, recent studies have highlighted an analogous epigenetic component as well. However, the relative roles of these mechanisms in facilitating population persistence under environmental heterogeneity remain unclear. To address the underlying genetic and epigenetic mechanisms and their relationship during environmental adaptation, we screened the genomes and epigenomes of nine global populations of a predominately sessile marine invasive tunicate, Botryllus schlosseri, using reduced-representation methods. We detected clear population differentiation at the genetic and epigenetic levels. Patterns of genetic and epigenetic structure were significantly influenced by local environmental variables. Among these variables, minimum annual sea surface temperature was identified as the top explanatory variable for both genetic and epigenetic variation. However, patterns of population structure driven by genetic and epigenetic variation were somewhat distinct, suggesting possible autonomy of epigenetic variation. We found both shared and specific genes and biological pathways among genetic and epigenetic loci associated with environmental factors, consistent with complementary and independent contributions of genetic and epigenetic variation to environmental adaptation in this system. Collectively, these mechanisms may facilitate population persistence under environmental change and sustain successful invasions across novel environments.

Adaptive divergence and underlying mechanisms in response to salinity gradients between two Crassostrea oysters revealed by phenotypic and transcriptomic analyses

Comparing the responses of closely related species to environmental changes is an efficient method to explore adaptive divergence, for a better understanding of the adaptive evolution of marine species under rapidly changing climates. Oysters are keystone species thrive in intertidal and estuarine areas where frequent environmental disturbance occurs including fluctuant salinity. The evolutionary divergence of two sister species of sympatric estuarine oysters, Crassostrea hongkongensis and Crassostrea ariakensis, in response to euryhaline habitats on phenotypes and gene expression, and the relative contribution of species effect, environment effect, and their interaction to the divergence were explored. After a 2-month outplanting at high- and low-salinity locations in the same estuary, the high growth rate, percent survival, and high tolerance indicated by physiological parameters suggested that the fitness of C. ariakensis was higher under high-salinity conditions and that of C. hongkongensis was higher under low-salinity conditions. Moreover, a transcriptomic analysis showed the two species exhibited differentiated transcriptional expression in high- and low-salinity habitats, largely caused by the species effect. Several of the important pathways enriched in divergent genes between species were also salinity-responsive pathways. Specifically, the pyruvate and taurine metabolism pathway and several solute carriers may contribute to the hyperosmotic adaptation of C. ariakensis, and some solute carriers may contribute to the hypoosmotic adaptation of C. hongkongensis. Our findings provide insights into the phenotypic and molecular mechanisms underlying salinity adaptation in marine mollusks, which will facilitate the assessment of the adaptive capacity of marine species in the context of climate change and will also provide practical information for marine resource conservation and aquaculture.