Impacts of climate change on geographical distributions of invasive ascidians

Non-native species in marine protected areas: Global distribution patterns

Marine protected areas (MPAs) across various countries have contributed to safeguarding coastal and marine environments. Despite these efforts, marine non-native species (NNS) continue to threaten biodiversity and ecosystems, even within MPAs. Currently, there is a lack of comprehensive studies on the inventories, distribution patterns, and effect factors of NNS within MPAs. Here we show a database containing over 15,000 occurrence records of 2714 marine NNS across 16,401 national or regional MPAs worldwide. To identify the primary mechanisms driving the occurrence of NNS, we utilize model selection with proxies representing colonization pressure, environmental variables, and MPA characteristics. Among the environmental predictors analyzed, sea surface temperature emerged as the sole factor strongly associated with NNS richness. Higher sea surface temperatures are linked to increased NNS richness, aligning with global marine biodiversity trends. Furthermore, human activities help species overcome geographical barriers and migration constraints. Consequently, this influences the distribution patterns of marine introduced species and associated environmental factors. As global climate change continues to alter sea temperatures, it is crucial to protect marine regions that are increasingly vulnerable to intense human activities and biological invasions.

Incorporating physiological knowledge into correlative species distribution models minimizes bias introduced by the choice of calibration area

Correlative species distribution models (SDMs) are important tools to estimate species' geographic distribution across space and time, but their reliability heavily relies on the availability and quality of occurrence data. Estimations can be biased when occurrences do not fully represent the environmental requirement of a species. We tested to what extent species' physiological knowledge might influence SDM estimations. Focusing on the Japanese sea cucumber Apostichopus japonicus within the coastal ocean of East Asia, we compiled a comprehensive dataset of occurrence records. We then explored the importance of incorporating physiological knowledge into SDMs by calibrating two types of correlative SDMs: a naïve model that solely depends on environmental correlates, and a physiologically informed model that further incorporates physiological information as priors. We further tested the models' sensitivity to calibration area choices by fitting them with different buffered areas around known presences. Compared with naïve models, the physiologically informed models successfully captured the negative influence of high temperature on A. japonicus and were less sensitive to the choice of calibration area. The naïve models resulted in more optimistic prediction of the changes of potential distributions under climate change (i.e., larger range expansion and less contraction) than the physiologically informed models. Our findings highlight benefits from incorporating physiological information into correlative SDMs, namely mitigating the uncertainties associated with the choice of calibration area. Given these promising features, we encourage future SDM studies to consider species physiological information where available.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00226-0.

Geographical distribution and ecological niche dynamics of Crassostrea sikamea (Amemiya, 1928) in China's coastal regions under climate change

Global climate change drives species redistribution, threatening biodiversity and ecosystem heterogeneity. The Kumamoto oyster, Crassostrea sikamea (Amemiya, 1928), one of the most promising aquaculture species because of its delayed reproductive timing, was once prevalent in southern China. In this study, an ensemble species distribution model was employed to analyze the distribution range shift and ecological niche dynamics of C. sikamea along China's coastline under the current and future climate scenarios (RCP 2.6-8.5 covering 2050 s and 2100 s). The model results indicated that the current habitat distribution for C. sikamea consists of a continuous stretch extending from the coastlines of Hainan Province to the northern shores of Jiangsu Province. By the 2050 s, the distribution range will stabilize at its southern end along the coast of Hainan Province, while expanding northward to cover the coastal areas of Shandong Province, showing a more dramatic trend of contraction in the south and invasion in the north by the 2100 s. In RCP8.5, the southern end retracts to the coasts of Guangdong, whereas the northern end covers all of China's coastal areas north of 34°N. C. sikamea can maintain relatively stable ecological niche characteristics, while it may occupy different ecological niche spaces under future climate conditions. Significant niche expansion will occur in lower temperature. We concluded C. sikamea habitats are susceptible to climate change. The rapid northward expansion of C. sikamea may open new possibilities for oyster farming in China, but it will also have important consequences for the ecological balance and biodiversity of receiving areas. It's imperative that we closely examine and strategize to address these repercussions for a win-win situation.

Rolling with the punches: Organism-environment interactions shape spatial pattern of adaptive differentiation in the widespread mantis shrimp Oratosquilla oratoria

Investigating spatial pattern of adaptive variation and its underlying processes can inform the adaptive potential distributed within species ranges, which is increasingly important in the context of a changing climate. A correct interpretation of adaptive variation pattern requires that population history and the ensuing population genetic structure are taken into account. Here we carried out such a study by integrating population genomic analyses, demographic model testing and species distribution modeling to investigate patterns and causes of adaptive differentiation in a widespread mantis shrimp, Oratosquilla oratoria, along a replicated, broad-scale temperature gradient in the northwestern Pacific (NWP). Our results supported a strong hierarchical ecogeographic structure dominated by habitat-linked divergence among O. oratoria populations accompanied with introgressive hybridization. A combined FST outlier and environmental correlation analyses revealed remarkable temperature-associated clines in allele frequency across paired North-South populations on Chinese and Japanese coasts, and identified a suite of loci associated with temperature adaptation. Further demographic model testing revealed the observed clinal variation derived partly from Pleistocene divergence followed by recent secondary contact. More importantly, the likelihood of hybridization is predicted to increase as climate change progresses, which would break barriers to gene flow and enable the spread of adaptive genetic variation. These results support that not only is temperature-driven adaptive differentiation occurs in O. oratoria but that such pattern is likely attributed to ancient adaptive variation, sustained by contemporary ocean conditions and a semi-permeable barrier to gene flow maintained by selection. They moreover provide genomic insights into the distribution of adaptive potential across O. oratoria' s species range. This work can serve as a case study to characterize adaptive diversity of marine species in the NWP by integrating environmental and genetic data at temporal and spatial scales in a population genomic framework, which would improve management and conservation actions under climate change.

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.

Interplays between cis- and trans-Acting Factors for Alternative Splicing in Response to Environmental Changes during Biological Invasions of Ascidians

Alternative splicing (AS), a pivotal biological process contributing to phenotypic plasticity, creates a bridge linking genotypes with phenotypes. Despite its importance, the AS mechanisms underlying environmental response and adaptation have not been well studied, and more importantly, the cis- and trans-acting factors influencing AS variation remain unclear. Using the model invasive congeneric ascidians, Ciona robusta, and Ciona savignyi, we compared their AS responses to environmental changes and explored the potential determinants. Our findings unveiled swift and dynamic AS changes in response to environmental challenges, and differentially alternative spliced genes (DASGs) were functionally enriched in transmembrane transport processes. Interestingly, both the prevalence and level of AS in C. robusta were lower than those observed in C. savignyi. Furthermore, these two indices were higher under temperature stresses compared to salinity stresses in C. savignyi. All the observed patterns underscore the species-specific and environmental context-dependent AS responses to environmental challenges. The dissimilarities in genomic structure and exon/intron size distributions between these two species likely contributed to the observed AS variation. Moreover, we identified a total of 11 and 9 serine/arginine-rich splicing factors (SRSFs) with conserved domains and gene structures in the genomes of C. robusta and C. savignyi, respectively. Intriguingly, our analysis revealed that all detected SRSFs did not exhibit prevalent AS regulations. Instead, we observed AS control over a set of genes related to splicing factors and spliceosome components. Altogether, our results elucidate species-specific and environmental challenge-dependent AS response patterns in closely related invasive ascidians. The identified splicing factors and spliceosome components under AS control offer promising candidates for further investigations into AS-mediated rapid responses to environmental challenges complementary to SRSFs.

Comparison of the Efficiency of Different Eradication Treatments to Minimize the Impacts Caused by the Invasive Tunicate Styela plicata in Mussel Aquaculture

In 2017, aquaculture producers of the Albufeira lagoon, Portugal, reported an invasion of tunicates that was disrupting mussel production, particularly the tunicate Styela plicata (Lesueur, 1823). A totally effective eradication method still does not exist, particularly for S. plicata, and the effects of the eradication treatments on bivalves' performance are also poorly understood. Our study examined the effectiveness of eradication treatments using three laboratory trials and five treatments (air exposure, freshwater immersion, sodium hypochlorite, hypersaline solution and acetic acid) for S. plicata, as well as their effects on survival and growth of blue mussel Mytilus edulis Linnaeus, 1758. While air exposure and freshwater immersion caused a 27% mortality rate in S. plicata, the acetic acid treatment was the most effective in eliminating this species (>90% mortality). However, a 33-40% mortality rate was registered in mussels. Both species were not affected by the hypersaline treatment in the last trial, but the sodium hypochlorite treatment led to a 57% mortality rate in mussels. Differences in mussels' growth rates were not detected. These trials represent a step forward in responding to the needs of aquaculture producers. However, further studies are needed to investigate the susceptibility of tunicates to treatments according to sexual maturation, as well as to ensure minimum mussel mortality in the most effective treatments, and to better understand the effects on mussel physiological performance in the long-term.

Impact of in Situ Simulated Climate Change on Communities and Non-Indigenous Species: Two Climates, Two Responses

Climate change constitutes a major challenge for marine urban ecosystems and ocean warming will likely strongly affect local communities. Non-Indigenous Species (NIS) have been shown to often have higher heat resistance than natives, but studies investigating how forthcoming global warming might affect them in marine urban environments remain scarce, especially in Situ studies. Here we used an in Situ warming experiment in a NW Mediterranean (warm temperate) and a NE Atlantic (cold temperate) marina to see how global warming might affect recruited communities in the near future. In both marinas, warming resulted in significantly different community structure, lower biomass, and more empty space compared to control. However, while in the warm temperate marina, NIS showed an increased surface cover, it was reduced in the cold temperate one. Metabolomic analyses on Bugula neritina in the Atlantic marina revealed potential heat stress experienced by this introduced bryozoan and a potential link between heat stress and the expression of a halogenated alkaloid, Caelestine A. The present results might indicate that the effects of global warming on the prevalence of NIS may differ between geographical provinces, which could be investigated by larger scale studies.

Prediction of the dynamic distribution for Eucheuma denticulatum (Rhodophyta, Solieriaceae) under climate change in the Indo-Pacific Ocean

Eucheuma is one of the most important commercial red seaweeds in Southeast Asia, and plays an important role in the global seaweed aquaculture. It is expected to exhibit great responses to ocean warming. Here, we used maximum entropy species distribution models (SDMs) to estimate the suitable habitat of Eucheuma denticulatum under present conditions, and to predict the future range dynamics under the four representative concentration pathway (RCP) scenarios. The best marine environmental factors for E. denticulatum distribution modeling were distance to shore, sea surface temperature and currents velocity. Our results showed that E. denticulatum' distributions would contract in the Central Indo-Pacific Ocean, especially the regions of the Sunda Shelf, while expanding poleward along the south coast of Australia in 2100. Our study provided important knowledge for the prediction of the tropical seaweed distribution, conservation and sustainable developments of E. denticulatum in the future.

Estimating the impact of climate change on the potential distribution of Indo-Pacific humpback dolphins with species distribution model

As IUCN critically vulnerable species,the Indo-Pacific humpback dolphins (Sousa chinensis) have attracted great public attention in recent years. The threats of human disturbance and environmental pollution to this population have been documented extensively. However, research on the sensitivity of this species to climate change is lacking. To understand the effect of climate change on the potential distribution of Sousa chinensis, we developed a weighted ensemble model based on 82 occurrence records and six predictor variables (e.g., ocean depth, distance to shore, mean temperature, salinity, ice thickness, and current velocity). According to the true skill statistic (TSS) and the area under the receiver operating characteristic curve (AUC), our ensemble model presented higher prediction precision than most of the single-algorithm models. It also indicated that ocean depth and distance to shore were the most important predictors in shaping the distribution patterns. The projections for the 2050s and 2100s from our ensemble model indicated a severe adverse impact of climate change on the Sousa chinensis habitat. Over 75% and 80% of the suitable habitat in the present day will be lost in all representative concentration pathway emission scenarios (RCPS) in the 2050s and 2100s, respectively. With the increased numbers of records of stranding and deaths of Sousa chinensis in recent years, strict management regulations and conservation plans are urgent to safeguard the current suitable habitats. Due to habitat contraction and poleward shift in the future, adaptive management strategies, including designing new reserves and adjusting the location and range of reserves according to the geographical distribution of Sousa chinensis, should be formulated to minimize the impacts of climate change on this species.

Intraspecific genetic variation matters when predicting seagrass distribution under climate change

Seagrasses play a vital role in structuring coastal marine ecosystems, but their distributional range and genetic diversity have declined rapidly in recent decades. To improve conservation of seagrass species, it is important to predict how climate change may impact their ranges. Such predictions are typically made with correlative species distribution models (SDMs), which can estimate a species' potential distribution under present and future climatic scenarios given species' presence data and climatic predictor variables. However, these models are typically constructed with species-level data, and thus ignore intraspecific genetic variability, which can give rise to populations with adaptations to heterogeneous climatic conditions. Here, we explore the link between intraspecific adaptation and niche differentiation in Thalassia hemprichii, a seagrass broadly distributed in the tropical Indo-Pacific Ocean and a crucial provider of habitat for numerous marine species. By retrieving and re-analysing microsatellite data from previous studies, we delimited two distinct phylogeographical lineages within the nominal species and found an intermediate level of differentiation in their multidimensional environmental niches, suggesting the possibility for local adaptation. We then compared projections of the species' habitat suitability under climate change scenarios using species-level and lineage-level SDMs. In the Central Tropical Indo-Pacific region, models for both levels predicted considerable range contraction in the future, but the lineage-level models predicted more severe habitat loss. Importantly, the two modelling approaches predicted opposite patterns of habitat change in the Western Tropical Indo-Pacific region. Our results highlight the necessity of conserving distinct populations and genetic pools to avoid regional extinction due to climate change and have important implications for guiding future management of seagrasses.

Marine macroinvertebrates inhabiting plastic litter in Peru

Marine litter, such as plastic bags, bottles, fabrics, or fishing gear, serve as a shelter for many marine organisms that are likely to colonize artificial substrata. Such assemblages can potentially turn marine litter into vectors of alien invasive species (AIS). Here, we report the abundance and diversity of macroinvertebrates inhabiting marine litter in Peruvian beaches. Results indicate that most of the fouled items found came from land-based sources (81.5%) and Bivalvia was the most abundant class (53.5%), mainly composed of the mussel Semimytilus algosus. No significant differences were found in the abundance and diversity of macroinvertebrates (class level) among sampling sites or sources of litter. Polypropylene and low-density polyethylene items were most frequently found with fouled biota. Although none of the identified species were non-native to the coast of Peru, we discuss marine litter as a potential source of AIS in this region.