Environmental DNA recovers fish composition turnover of the coral reefs of West Indian Ocean islands.

Islands have been used as model systems to study ecological and evolutionary processes, and they provide an ideal set-up for validating new biodiversity monitoring methods. The application of environmental DNA metabarcoding for monitoring marine biodiversity requires an understanding of the spatial scale of the eDNA signal, which is best tested in island systems. Here, we investigated the variation in Actinopterygii and Elasmobranchii species composition recovered from eDNA metabarcoding along a gradient of distance-to-reef in four of the five French Scattered Islands in the Western Indian Ocean. We collected surface water samples at an increasing distance from reefs (0 m, 250 m, 500 m, 750 m). We used a metabarcoding protocol based on the 'teleo' primers to target marine reef fishes and classified taxa according to their habitat types (benthic or pelagic). We investigated the effect of distance-to-reef on ß diversity variation using generalised linear mixed models and estimated species-specific distance-to-reef effects using a model-based approach for community data. Environmental DNA metabarcoding analyses recovered distinct fish species compositions across the four inventoried islands and variations along the distance-to-reef gradient. The analysis of ß-diversity variation showed significant taxa turnover between the eDNA samples on and away from the reefs. In agreement with a spatially localised signal from eDNA, benthic species were distributed closer to the reef than pelagic ones. Our findings demonstrate that the combination of eDNA inventories and spatial modelling can provide insights into species habitat preferences related to distance-to-reef gradients at a small scale. As such, eDNA can not only recover large compositional differences among islands but also help understand habitat selection and distribution of marine species at a finer spatial scale.

The Tree of Life eDNA metabarcoding reveals a similar taxonomic richness but dissimilar evolutionary lineages between seaports and marine reserves.

Coastal areas host a major part of marine biodiversity but are seriously threatened by ever-increasing human pressures. Transforming natural coastlines into urban seascapes through habitat artificialization may result in loss of biodiversity and key ecosystem functions. Yet, the extent to which seaports differ from nearby natural habitats and marine reserves across the whole Tree of Life is still unknown. This study aimed to assess the level of α and ß-diversity between seaports and reserves, and whether these biodiversity patterns are conserved across taxa and evolutionary lineages. For that, we used environmental DNA (eDNA) metabarcoding to survey six seaports on the French Mediterranean coast and four strictly no-take marine reserves nearby. By targeting four different groups-prokaryotes, eukaryotes, metazoans and fish-with appropriate markers, we provide a holistic view of biodiversity on contrasted habitats. In the absence of comprehensive reference databases, we used bioinformatic pipelines to gather similar sequences into molecular operational taxonomic units (MOTUs). In contrast to our expectations, we obtained no difference in MOTU richness (α-diversity) between habitats except for prokaryotes and threatened fishes with higher diversity in reserves than in seaports. However, we observed a marked dissimilarity (ß-diversity) between seaports and reserves for all taxa. Surprisingly, this biodiversity signature of seaports was preserved across the Tree of Life, up to the order. This result reveals that seaports and nearby marine reserves share few taxa and evolutionary lineages along urbanized coasts and suggests major differences in terms of ecosystem functioning between both habitats.

Early signals of Posidonia oceanica meadows recovery in a context of wastewater treatment improvements.

Natural ecological restoration is a cornerstone of modern conservation science and managers need more documented "success stories" to lead the way. In French mediterranean sea, we monitored Posidonia oceanica lower limit using acoustic telemetry and photogrammetry and investigated the descriptors driving its variations, at a national scale and over more than a decade. We showed significant effects of environmental descriptors (region, sea surface temperature and bottom temperature) but also of wastewater treatment plant (WWTP) effluents proxies (size of WWTP, time since conformity, and distance to the closest effluent) on the meadows lower limit progression. This work indicates a possible positive response of P. oceanica meadows to improvements in wastewater treatment and a negative effect of high temperatures. While more data is needed, the example of French wastewater policy should inspire stakeholders and coastal managers in their efforts to limit anthropogenic pressures on vulnerable ecosystems.

Rarity mediates species-specific responses of tropical reef fishes to protection.

Marine protected areas (MPAs) are the most widely applied tool for marine biodiversity conservation, yet many gaps remain in our understanding of their species-specific effects, partly because the socio-environmental context and spatial autocorrelation may blur and bias perceived conservation outcomes. Based on a large data set of nearly 3000 marine fish surveys spanning all tropical regions of the world, we build spatially explicit models for 658 fish species to estimate species-specific responses to protection while controlling for the environmental, habitat and socio-economic contexts experienced across their geographic ranges. We show that the species responses are highly variable, with ~40% of fishes not benefitting from protection. When investigating how traits influence species' responses, we find that rare top-predators and small herbivores benefit the most from MPAs while mid-trophic level species benefit to a lesser extent, and rare large herbivores experience adverse effects, indicating potential trophic cascades.

Divergent responses of pelagic and benthic fish body-size structure to remoteness and protection from humans.

Animal body-size variation influences multiple processes in marine ecosystems, but habitat heterogeneity has prevented a comprehensive assessment of size across pelagic (midwater) and benthic (seabed) systems along anthropic gradients. In this work, we derive fish size indicators from 17,411 stereo baited-video deployments to test for differences between pelagic and benthic responses to remoteness from human pressures and effectiveness of marine protected areas (MPAs). From records of 823,849 individual fish, we report divergent responses between systems, with pelagic size structure more profoundly eroded near human markets than benthic size structure, signifying greater vulnerability of pelagic systems to human pressure. Effective protection of benthic size structure can be achieved through MPAs placed near markets, thereby contributing to benthic habitat restoration and the recovery of associated fishes. By contrast, recovery of the world's largest and most endangered fishes in pelagic systems requires the creation of highly protected areas in remote locations, including on the High Seas, where protection efforts lag.

Functional diversity of sharks and rays is highly vulnerable and supported by unique species and locations worldwide.

Elasmobranchs (sharks, rays and skates) are among the most threatened marine vertebrates, yet their global functional diversity remains largely unknown. Here, we use a trait dataset of >1000 species to assess elasmobranch functional diversity and compare it against other previously studied biodiversity facets (taxonomic and phylogenetic), to identify species- and spatial- conservation priorities. We show that threatened species encompass the full extent of functional space and disproportionately include functionally distinct species. Applying the conservation metric FUSE (Functionally Unique, Specialised, and Endangered) reveals that most top-ranking species differ from the top Evolutionarily Distinct and Globally Endangered (EDGE) list. Spatial analyses further show that elasmobranch functional richness is concentrated along continental shelves and around oceanic islands, with 18 distinguishable hotspots. These hotspots only marginally overlap with those of other biodiversity facets, reflecting a distinct spatial fingerprint of functional diversity. Elasmobranch biodiversity facets converge with fishing pressure along the coast of China, which emerges as a critical frontier in conservation. Meanwhile, several components of elasmobranch functional diversity fall in high seas and/or outside the global network of marine protected areas. Overall, our results highlight acute vulnerability of the world's elasmobranchs' functional diversity and reveal global priorities for elasmobranch functional biodiversity previously overlooked.

Sustainable reference points for multispecies coral reef fisheries.

Sustainably managing fisheries requires regular and reliable evaluation of stock status. However, most multispecies reef fisheries around the globe tend to lack research and monitoring capacity, preventing the estimation of sustainable reference points against which stocks can be assessed. Here, combining fish biomass data for >2000 coral reefs, we estimate site-specific sustainable reference points for coral reef fisheries and use these and available catch estimates to assess the status of global coral reef fish stocks. We reveal that >50% of sites and jurisdictions with available information have stocks of conservation concern, having failed at least one fisheries sustainability benchmark. We quantify the trade-offs between biodiversity, fish length, and ecosystem functions relative to key benchmarks and highlight the ecological benefits of increasing sustainability. Our approach yields multispecies sustainable reference points for coral reef fisheries using environmental conditions, a promising means for enhancing the sustainability of the world's coral reef fisheries.

New deep learning-based methods for visualizing ecosystem properties using environmental DNA metabarcoding data.

Environmental DNA (eDNA) metabarcoding provides an efficient approach for documenting biodiversity patterns in marine and terrestrial ecosystems. The complexity of these data prevents current methods from extracting and analyzing all the relevant ecological information they contain, and new methods may provide better dimensionality reduction and clustering. Here we present two new deep learning-based methods that combine different types of neural networks (NNs) to ordinate eDNA samples and visualize ecosystem properties in a two-dimensional space: the first is based on variational autoencoders and the second on deep metric learning. The strength of our new methods lies in the combination of two inputs: the number of sequences found for each molecular operational taxonomic unit (MOTU) detected and their corresponding nucleotide sequence. Using three different datasets, we show that our methods accurately represent several biodiversity indicators in a two-dimensional latent space: MOTU richness per sample, sequence α-diversity per sample, Jaccard's and sequence ß-diversity between samples. We show that our nonlinear methods are better at extracting features from eDNA datasets while avoiding the major biases associated with eDNA. Our methods outperform traditional dimension reduction methods such as Principal Component Analysis, t-distributed Stochastic Neighbour Embedding, Nonmetric Multidimensional Scaling and Uniform Manifold Approximation and Projection for dimension reduction. Our results suggest that NNs provide a more efficient way of extracting structure from eDNA metabarcoding data, thereby improving their ecological interpretation and thus biodiversity monitoring.

Industrial fisheries have reversed the carbon sequestration by tuna carcasses into emissions.

To limit climate warming to 2°C above preindustrial levels, most economic sectors will need a rapid transformation toward a net zero emission of CO2 . Tuna fisheries is a key food production sector that burns fossil fuel to operate but also reduces the deadfall of large-bodied fish so the capacity of this natural carbon pump to deep sea. Yet, the carbon balance of tuna populations, so the net difference between CO2 emission due to industrial exploitation and CO2 sequestration by fish deadfall after natural mortality, is still unknown. Here, by considering the dynamics of two main contrasting tuna species (Katsuwonus pelamis and Thunnus obesus) across the Pacific since the 1980s, we show that most tuna populations became CO2 sources instead of remaining natural sinks. Without considering the supply chain, the main factors associated with this shift are exploitation rate, transshipment intensity, fuel consumption, and climate change. Our study urges for a better global ocean stewardship, by curbing subsidies and limiting transshipment in remote international waters, to quickly rebuild most pelagic fish stocks above their target management reference points and reactivate a neglected carbon pump toward the deep sea as an additional Nature Climate Solution in our portfolio. Even if this potential carbon sequestration by surface unit may appear low compared to that of coastal ecosystems or tropical forests, the ocean covers a vast area and the sinking biomass of dead vertebrates can sequester carbon for around 1000 years in the deep sea. We also highlight the multiple co-benefits and trade-offs from engaging the industrial fisheries sector with carbon neutrality.

A marine protected area network does not confer community structure resilience to a marine heatwave across coastal ecosystems.

Marine protected areas (MPAs) have gained attention as a conservation tool for enhancing ecosystem resilience to climate change. However, empirical evidence explicitly linking MPAs to enhanced ecological resilience is limited and mixed. To better understand whether MPAs can buffer climate impacts, we tested the resistance and recovery of marine communities to the 2014-2016 Northeast Pacific heatwave in the largest scientifically designed MPA network in the world off the coast of California, United States. The network consists of 124 MPAs (48 no-take state marine reserves, and 76 partial-take or special regulation conservation areas) implemented at different times, with full implementation completed in 2012. We compared fish, benthic invertebrate, and macroalgal community structure inside and outside of 13 no-take MPAs across rocky intertidal, kelp forest, shallow reef, and deep reef nearshore habitats in California's Central Coast region from 2007 to 2020. We also explored whether MPA features, including age, size, depth, proportion rock, historic fishing pressure, habitat diversity and richness, connectivity, and fish biomass response ratios (proxy for ecological performance), conferred climate resilience for kelp forest and rocky intertidal habitats spanning 28 MPAs across the full network. Ecological communities dramatically shifted due to the marine heatwave across all four nearshore habitats, and MPAs did not facilitate habitat-wide resistance or recovery. Only in protected rocky intertidal habitats did community structure significantly resist marine heatwave impacts. Community shifts were associated with a pronounced decline in the relative proportion of cold water species and an increase in warm water species. MPA features did not explain resistance or recovery to the marine heatwave. Collectively, our findings suggest that MPAs have limited ability to mitigate the impacts of marine heatwaves on community structure. Given that mechanisms of resilience to climate perturbations are complex, there is a clear need to expand assessments of ecosystem-wide consequences resulting from acute climate-driven perturbations, and the potential role of regulatory protection in mitigating community structure changes.

The ecological causes of functional distinctiveness in communities.

Recent work has shown that evaluating functional trait distinctiveness, the average trait distance of a species to other species in a community offers promising insights into biodiversity dynamics and ecosystem functioning. However, the ecological mechanisms underlying the emergence and persistence of functionally distinct species are poorly understood. Here, we address the issue by considering a heterogeneous fitness landscape whereby functional dimensions encompass peaks representing trait combinations yielding positive population growth rates in a community. We identify four ecological cases contributing to the emergence and persistence of functionally distinct species. First, environmental heterogeneity or alternative phenotypic designs can drive positive population growth of functionally distinct species. Second, sink populations with negative population growth can deviate from local fitness peaks and be functionally distinct. Third, species found at the margin of the fitness landscape can persist but be functionally distinct. Fourth, biotic interactions (positive or negative) can dynamically alter the fitness landscape. We offer examples of these four cases and guidelines to distinguish between them. In addition to these deterministic processes, we explore how stochastic dispersal limitation can yield functional distinctiveness. Our framework offers a novel perspective on the relationship between fitness landscape heterogeneity and the functional composition of ecological assemblages.

Temperature, species identity and morphological traits predict carbonate excretion and mineralogy in tropical reef fishes.

Anthropogenic pressures are restructuring coral reefs globally. Sound predictions of the expected changes in key reef functions require adequate knowledge of their drivers. Here we investigate the determinants of a poorly-studied yet relevant biogeochemical function sustained by marine bony fishes: the excretion of intestinal carbonates. Compiling carbonate excretion rates and mineralogical composition from 382 individual coral reef fishes (85 species and 35 families), we identify the environmental factors and fish traits that predict them. We find that body mass and relative intestinal length (RIL) are the strongest predictors of carbonate excretion. Larger fishes and those with longer intestines excrete disproportionately less carbonate per unit mass than smaller fishes and those with shorter intestines. The mineralogical composition of excreted carbonates is highly conserved within families, but also controlled by RIL and temperature. These results fundamentally advance our understanding of the role of fishes in inorganic carbon cycling and how this contribution will change as community composition shifts under increasing anthropogenic pressures.

The global depth range of marine fishes and their genetic coverage for environmental DNA metabarcoding.

The bathymetric and geographical distribution of marine species represent a key information in biodiversity conservation. Yet, deep-sea ecosystems are among the least explored on Earth and are increasingly impacted by human activities. Environmental DNA (eDNA) metabarcoding has emerged as a promising method to study fish biodiversity but applications to the deep-sea are still scarce. A major limitation in the application of eDNA metabarcoding is the incompleteness of species sequences available in public genetic databases which reduces the extent of detected species. This incompleteness by depth is still unknown. Here, we built the global bathymetric and geographical distribution of 10,826 actinopterygian and 960 chondrichthyan fish species. We assessed their genetic coverage by depth and by ocean for three main metabarcoding markers used in the literature: teleo and MiFish-U/E. We also estimated the number of primer mismatches per species amplified by in silico polymerase chain reaction which influence the probability of species detection. Actinopterygians show a stronger decrease in species richness with depth than Chondrichthyans. These richness gradients are accompanied by a continuous species turnover between depths. Fish species coverage with the MiFish-U/E markers is higher than with teleo while threatened species are more sequenced than the others. "Deep-endemic" species, those not ascending to the shallow depth layer, are less sequenced than not threatened species. The number of primer mismatches is not higher for deep-sea species than for shallower ones. eDNA metabarcoding is promising for species detection in the deep-sea to better account for the 3-dimensional structure of the ocean in marine biodiversity monitoring and conservation. However, we argue that sequencing efforts on "deep-endemic" species are needed.

Recovery of carbon benefits by overharvested baleen whale populations is threatened by climate change.

Despite the importance of marine megafauna on ecosystem functioning, their contribution to the oceanic carbon cycle is still poorly known. Here, we explored the role of baleen whales in the biological carbon pump across the southern hemisphere based on the historical and forecasted abundance of five baleen whale species. We modelled whale-mediated carbon sequestration through the sinking of their carcasses after natural death. We provide the first temporal dynamics of this carbon pump from 1890 to 2100, considering both the effects of exploitation and climate change on whale populations. We reveal that at their pre-exploitation abundance, the five species of southern whales could sequester 4.0 × 105 tonnes of carbon per year (tC yr-1). This estimate dropped to 0.6 × 105 tC yr-1 by 1972 following commercial whaling. However, with the projected restoration of whale populations under a RCP8.5 climate scenario, the sequestration would reach 1.7 × 105 tC yr-1 by 2100, while without climate change, recovered whale populations could sequester nearly twice as much (3.2 × 105 tC yr-1) by 2100. This highlights the persistence of whaling damages on whale populations and associated services as well as the predicted harmful impacts of climate change on whale ecosystem services.

Safeguarding nutrients from coral reefs under climate change.

The sustainability of coral reef fisheries is jeopardized by complex and interacting socio-ecological stressors that undermine their contribution to food and nutrition security. Climate change has emerged as one of the key stressors threatening coral reefs and their fish-associated services. How fish nutrient concentrations respond to warming oceans remains unclear but these responses are probably affected by both direct (metabolism and trophodynamics) and indirect (habitat and species range shifts) effects. Climate-driven coral habitat loss can cause changes in fish abundance and biomass, revealing potential winners and losers among major fisheries targets that can be predicted using ecological indicators and biological traits. A critical next step is to extend research focused on the quantity of available food (fish biomass) to also consider its nutritional quality, which is relevant to progress in the fields of food security and malnutrition. Biological traits are robust predictors of fish nutrient content and thus potentially indicate how climate-driven changes are expected to impact nutrient availability within future food webs on coral reefs. Here, we outline future research priorities and an anticipatory framework towards sustainable reef fisheries contributing to nutrition-sensitive food systems in a warming ocean.

A functional vulnerability framework for biodiversity conservation.

Setting appropriate conservation strategies in a multi-threat world is a challenging goal, especially because of natural complexity and budget limitations that prevent effective management of all ecosystems. Safeguarding the most threatened ecosystems requires accurate and integrative quantification of their vulnerability and their functioning, particularly the potential loss of species trait diversity which imperils their functioning. However, the magnitude of threats and associated biological responses both have high uncertainties. Additionally, a major difficulty is the recurrent lack of reference conditions for a fair and operational measurement of vulnerability. Here, we present a functional vulnerability framework that incorporates uncertainty and reference conditions into a generalizable tool. Through in silico simulations of disturbances, our framework allows us to quantify the vulnerability of communities to a wide range of threats. We demonstrate the relevance and operationality of our framework, and its global, scalable and quantitative comparability, through three case studies on marine fishes and mammals. We show that functional vulnerability has marked geographic and temporal patterns. We underline contrasting contributions of species richness and functional redundancy to the level of vulnerability among case studies, indicating that our integrative assessment can also identify the drivers of vulnerability in a world where uncertainty is omnipresent.

Applying convolutional neural networks to speed up environmental DNA annotation in a highly diverse ecosystem.

High-throughput DNA sequencing is becoming an increasingly important tool to monitor and better understand biodiversity responses to environmental changes in a standardized and reproducible way. Environmental DNA (eDNA) from organisms can be captured in ecosystem samples and sequenced using metabarcoding, but processing large volumes of eDNA data and annotating sequences to recognized taxa remains computationally expensive. Speed and accuracy are two major bottlenecks in this critical step. Here, we evaluated the ability of convolutional neural networks (CNNs) to process short eDNA sequences and associate them with taxonomic labels. Using a unique eDNA data set collected in highly diverse Tropical South America, we compared the speed and accuracy of CNNs with that of a well-known bioinformatic pipeline (OBITools) in processing a small region (60 bp) of the 12S ribosomal DNA targeting freshwater fishes. We found that the taxonomic labels from the CNNs were comparable to those from OBITools, with high correlation levels for the composition of the regional fish fauna. The CNNs enabled the processing of raw fastq files at a rate of approximately 1 million sequences per minute, which was about 150 times faster than with OBITools. Given the good performance of CNNs in the highly diverse ecosystem considered here, the development of more elaborate CNNs promises fast deployment for future biodiversity inventories using eDNA.