Distinct interspecific and intraspecific vulnerability of coastal species to global change

Predictive modelling reveals Australian continental risk hotspots for marine debris interactions with key threatened species

Anthropogenic debris is a global threat that impacts threatened species through various lethal and sub-lethal consequences, as well as overall ecosystem health. This study used a database of over 24,000 beach surveys of marine debris collated by the Australian Marine Debris Initiative from 2012 to 2021, with two key objectives: (1) identify variables that most influence the occurrence of debris hotspots on a continental scale and (2) use these findings to identify likely hotspots of interaction between threatened species and marine debris. The number of particles found in each beach survey was modelled alongside fifteen biological, social, and physical spatial variables including land use, physical oceanography, population, rainfall, distance to waste facilities, ports, and mangroves to identify the significant drivers of debris deposition. The model of best fit for predicting debris particle abundance was calculated using a generalized additive model. Overall, debris was more abundant at sites near catchments with high annual rainfall (mm), intensive land use (km2), and that were nearer to ports (km) and mangroves (km). These results support previous studies which state that mangroves are a significant sink for marine debris, and that large ports and urbanized catchments are significant sources for marine debris. We illustrate the applicability of these models by quantifying significant overlap between debris hotspots and the distributions for four internationally listed threatened species that exhibit debris interactions; green turtle (26,868 km2), dugong (16,164 km2), Australian sea lion (2903 km2) and Flesh-footed Shearwater (2413 km2). This equates to less than 1% (Flesh-footed Shearwater, Australian sea lion), over 2% (green sea turtle) and over 5% (dugong) of their habitat being identified as areas of high risk for marine debris interactions. The results of this study hold practical value, informing decision-making processes, managing debris pollution at continental scales, as well as identifying gaps in species monitoring.

Local adaptation of Pinus leiophylla under climate and land use change models in the Avocado Belt of Michoacán

Climate change and land use change are two main drivers of global biodiversity decline, decreasing the genetic diversity that populations harbour and altering patterns of local adaptation. Landscape genomics allows measuring the effect of these anthropogenic disturbances on the adaptation of populations. However, both factors have rarely been considered simultaneously. Based on a set of 3660 SNPs from which 130 were identified as outliers by a genome-environment association analysis (LFMM), we modelled the spatial turnover of allele frequencies in 19 localities of Pinus leiophylla across the Avocado Belt in Michoacán state, Mexico. Then, we evaluated the effect of climate change and land use change scenarios, in addition to evaluating assisted gene flow strategies and connectivity metrics across the landscape to identify priority conservation areas for the species. We found that localities in the centre-east of the Avocado Belt would be more vulnerable to climate change, while localities in the western area are more threatened by land conversion to avocado orchards. Assisted gene flow actions could aid in mitigating both threats. Connectivity patterns among forest patches will also be modified by future habitat loss, with central and eastern parts of the Avocado Belt maintaining the highest connectivity. These results suggest that areas with the highest priority for conservation are in the eastern part of the Avocado Belt, including the Monarch Butterfly Biosphere Reserve. This work is useful as a framework that incorporates distinct layers of information to provide a more robust representation of the response of tree populations to anthropogenic disturbances.

Variation in genomic vulnerability to climate change across temperate populations of eelgrass (Zostera marina)

A global decline in seagrass populations has led to renewed calls for their conservation as important providers of biogenic and foraging habitat, shoreline stabilization and carbon storage. Eelgrass (Zostera marina) occupies the largest geographic range among seagrass species spanning a commensurately broad spectrum of environmental conditions. In Canada, eelgrass is managed as a single phylogroup despite occurring across three oceans and a range of ocean temperatures and salinity gradients. Previous research has focused on applying relatively few markers to reveal population structure of eelgrass, whereas a whole-genome approach is warranted to investigate cryptic structure among populations inhabiting different ocean basins and localized environmental conditions. We used a pooled whole-genome re-sequencing approach to characterize population structure, gene flow and environmental associations of 23 eelgrass populations ranging from the Northeast United States to Atlantic, subarctic and Pacific Canada. We identified over 500,000 SNPs, which when mapped to a chromosome-level genome assembly revealed six broad clades of eelgrass across the study area, with pairwise F ST ranging from 0 among neighbouring populations to 0.54 between Pacific and Atlantic coasts. Genetic diversity was highest in the Pacific and lowest in the subarctic, consistent with colonization of the Arctic and Atlantic oceans from the Pacific less than 300 kya. Using redundancy analyses and two climate change projection scenarios, we found that subarctic populations are predicted to be potentially more vulnerable to climate change through genomic offset predictions. Conservation planning in Canada should thus ensure that representative populations from each identified clade are included within a national network so that latent genetic diversity is protected, and gene flow is maintained. Northern populations, in particular, may require additional mitigation measures given their potential susceptibility to a rapidly changing climate.

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.

Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species

Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.

Research on the Changes in Distribution and Habitat Suitability of the Chinese Red Panda Population

The study of the dynamics of species habitat is of great significance for maintaining or adjusting the current habitat protection management strategy. However, the current research on the Chinese red panda's habitat is limited to the analysis of a single period, which makes it difficult to quantify the changes in its habitat on a temporal scale and greatly hinders the formulation of the overall protection and management strategies that are to be used for the Chinese red panda. This study simulated habitat suitability at different temporal scales to quantify the trend of changes in habitat quality and analyzed the reasons for the changes in habitat suitability in certain regions. The results showed that the overall suitability of the Chinese red panda's habitat increased and that the area of suitable habitats expanded. Suitable Chinese red panda habitats in the mountains of Qionglai (1662.73 km2), Daxiangling (230.30 km2), Xiaoxiangling (549.47 km2), and Liangshan (50.39 km2) increased by a total of 2452.89 km2. The suitability of habitats in the central part of the Liangshan Mountains has declined significantly, which is positively correlated with changes in temperature seasonality (BIO4, R = 0.18) and negatively correlated with changes in annual average temperature (BIO1, R = -0.03) as well as changes in the proportion of farmland (FARMLAND, R = -0.14). The local extinction of isolated populations of Chinese red pandas in the Minshan Mountains is the main factor leading to their distribution retreat rather than a decrease in habitat quality. The research results help us to provide a scientific basis for the formulation of conservation and management strategies for Chinese red pandas at different scales.

Assessing the global vulnerability of dryland birds to heatwaves

As global average surface temperature increases, extreme climatic events such as heatwaves are becoming more frequent and intense, which can drive biodiversity responses such as rapid population declines and/or shifts in species distributions and even local extirpations. However, the impacts of extreme climatic events are largely ignored in conservation plans. Birds are known to be susceptible to heatwaves, especially in dryland ecosystems. Understanding which birds are most vulnerable to heatwaves, and where these birds occur, can offer a scientific basis for adaptive management and conservation. We assessed the relative vulnerability of 1196 dryland bird species to heatwaves using a trait-based approach. Among them, 888 bird species are estimated to be vulnerable to heatwaves (170 highly vulnerable, eight extremely vulnerable), of which ~91% are currently considered non-threatened by the IUCN, which suggests that many species will likely become newly threatened with intensifying climate change. We identified the top three hotspot areas of heatwave-vulnerable species in Australia (208 species), Southern Africa (125 species) and Eastern Africa (99 species). Populations of vulnerable species recorded in the Living Planet Database were found to be declining significantly faster than those of non-vulnerable species (p = .048) after heatwaves occurred. In contrast, no significant difference in population trends between vulnerable and non-vulnerable species was detected when no heatwave occurred (p = .34). This suggests that our vulnerability framework correctly identified vulnerable species and that heatwaves are already impacting the population trends of these species. Our findings will help prioritize heatwave-vulnerable birds in dryland ecosystems in risk mitigation and adaptation management as the frequency of heatwaves accelerates in the coming decades.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Evaluation of species thermal sensitivity with individual-based physiological performance

Ignoring intraspecific variations can prevent us from accurately assessing species' thermal sensitivity to global warming. Individual-based physiological performance provides a feasible solution to depict species' thermal sensitivity using a bottom-up approach. We measured the cardiac performance of intertidal bivalves (1159 individuals from multiple populations of six bivalves), determined the upper thermal limit of each individual, calculated the proportions of individuals suffering sublethal/lethal heat stress, and mapped sensitive regions to high temperatures. Results showed that high inter-individual variations of physiological performance existed in levels of populations and species, and species' thermal sensitivity was positively related to the intraspecific variations of heat tolerance. This bottom-up approach scaled up from individual, population to species emphasizes the importance of individual-based physiology performance in assessing thermal sensitivity across different hierarchical levels and enables better evaluating and forecasting of species responses to global warming.

Responses of population structure and genomic diversity to climate change and fishing pressure in a pelagic fish

The responses of marine species to environmental changes and anthropogenic pressures (e.g., fishing) interact with ecological and evolutionary processes that are not well understood. Knowledge of changes in the distribution range and genetic diversity of species and their populations into the future is essential for the conservation and sustainable management of resources. Almaco jack (Seriola rivoliana) is a pelagic fish with high importance to fisheries and aquaculture in the Pacific Ocean. In this study, we assessed contemporary genomic diversity and structure in loci that are putatively under selection (outlier loci) and determined their potential functions. Using a combination of genotype-environment association, spatial distribution models, and demogenetic simulations, we modeled the effects of climate change (under three different RCP scenarios) and fishing pressure on the species' geographic distribution and genomic diversity and structure to 2050 and 2100. Our results show that most of the outlier loci identified were related to biological and metabolic processes that may be associated with temperature and salinity. The contemporary genomic structure showed three populations-two in the Eastern Pacific (Cabo San Lucas and Eastern Pacific) and one in the Central Pacific (Hawaii). Future projections suggest a loss of suitable habitat and potential range contractions for most scenarios, while fishing pressure decreased population connectivity. Our results suggest that future climate change scenarios and fishing pressure will affect the genomic structure and genotypic composition of S. rivoliana and lead to loss of genomic diversity in populations distributed in the eastern-central Pacific Ocean, which could have profound effects on fisheries that depend on this resource.

The effects of intraspecific variation on forecasts of species range shifts under climate change

As global climate change is altering the distribution range of macroalgae across the globe, it is critical to assess its impact on species range shifts to inform the biodiversity conservation of macroalgae. Latitude/environmental gradients could cause intraspecific variability, which may result in distinct responses to climate change. It remains unclear whether geographical variation occurs in the response of species' populations to climate change. We tested this assumption using the brown alga Sargassum thunbergii, a habitat-forming macroalgae encompassing multiple divergent lineages along the Northwest Pacific. Previous studies revealed a distinct lineage of S. thunbergii in rear-edge populations. Given the phylogeographic structure and temperature gradients, we divided these populations into the southern and northern groups. We assessed the physiological responses of the two groups to temperature changes and estimated their niche differences using n-dimensional hypervolumes. A higher photosynthetic rate and antioxidative abilities were detected in the southern group of S. thunbergii than in the northern group. In addition, significant niche differentiation was detected between the two groups, suggesting the possibility for local adaptation. Given these results, we inferred that the southern group (rear-edge populations) may be more resilient to climate change. To examine climate-driven range shifts of S. thunbergii, we constructed species- and lineage-level species distribution models (SDMs). Predictions of both levels showed considerable distribution contracts along the Chinese coasts in the future. For the southern group, the lineage-level model predicted less habitat loss than the species-level model. Our results highlight the importance of considering intraspecific variation in climate change vulnerability assessments for coastal species.

Climate adaptation and vulnerability of foundation species in a global change hotspot

Climate change is altering species ranges, and relative abundances within ranges, as populations become differentially adapted and vulnerable to the climates they face. Understanding present species ranges, whether species harbour and exchange adaptive variants, and how variants are distributed across landscapes undergoing rapid change, is therefore crucial to predicting responses to future climates and informing conservation strategies. Such insights are nonetheless lacking for most species of conservation concern. We assess genomic patterns of neutral variation, climate adaptation and climate vulnerability (offsets in predicted distributions of putatively adaptive variants across present and future landscapes) for sister foundation species, the marine tubeworms Galeolaria caespitosa and Galeolaria gemineoa, in a sentinel region for climate change impacts. We find that species are genetically isolated despite uncovering sympatry in their ranges, show parallel and nonparallel signals of thermal adaptation on spatial scales smaller than gene flow across their ranges, and are predicted to face different risks of maladaptation under future temperatures across their ranges. Our findings have implications for understanding local adaptation in the face of gene flow, and generate spatially explicit predictions for climatic disruption of adaptation and species distributions in coastal ecosystems that could guide experimental validation and conservation planning.

Ecological and genomic vulnerability to climate change across native populations of Robusta coffee (Coffea canephora)

The assessment of population vulnerability under climate change is crucial for planning conservation as well as for ensuring food security. Coffea canephora is, in its native habitat, an understorey tree that is mainly distributed in the lowland rainforests of tropical Africa. Also known as Robusta, its commercial value constitutes a significant revenue for many human populations in tropical countries. Comparing ecological and genomic vulnerabilities within the species' native range can provide valuable insights about habitat loss and the species' adaptive potential, allowing to identify genotypes that may act as a resource for varietal improvement. By applying species distribution models, we assessed ecological vulnerability as the decrease in climatic suitability under future climatic conditions from 492 occurrences. We then quantified genomic vulnerability (or risk of maladaptation) as the allelic composition change required to keep pace with predicted climate change. Genomic vulnerability was estimated from genomic environmental correlations throughout the native range. Suitable habitat was predicted to diminish to half its size by 2050, with populations near coastlines and around the Congo River being the most vulnerable. Whole-genome sequencing revealed 165 candidate SNPs associated with climatic adaptation in C. canephora, which were located in genes involved in plant response to biotic and abiotic stressors. Genomic vulnerability was higher for populations in West Africa and in the region at the border between DRC and Uganda. Despite an overall low correlation between genomic and ecological vulnerability at broad scale, these two components of vulnerability overlap spatially in ways that may become damaging. Genomic vulnerability was estimated to be 23% higher in populations where habitat will be lost in 2050 compared to regions where habitat will remain suitable. These results highlight how ecological and genomic vulnerabilities are relevant when planning on how to cope with climate change regarding an economically important species.

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.