Virginie Marques

Habitat configuration governs the movement of organisms across landscapes, thereby shaping both population structure and community assembly. While theoretical and empirical studies have assessed how habitat connectivity simultaneously influences intra‐ and interspecific diversity, direct comparisons across contrasting biogeographic regions remain limited. Here, we investigate patterns of genetic and species β‐diversity in tropical reef fishes across two ocean basins with distinct spatial configurations: the Caribbean Sea and the Western Indian Ocean. Using a comparative framework based on species occurrence data from five fish families and single nucleotide polymorphism (SNP) data from 19 species, we detected significant isolation by distance at both population and community levels in the Western Indian Ocean, but only at the community level in the Caribbean Sea. Additionally, genetic and species β‐diversity were positively correlated among species in the Western Indian Ocean, but not in the Caribbean Sea. Together, these results suggest that the shorter inter‐reef distances of the Caribbean Sea promote higher connectivity, leading to a decoupling of intra‐ and interspecific β‐diversity patterns.

Understanding the intricate dynamics of biodiversity within and across riverine ecosystems, influenced by geological history and environmental factors, is crucial for effective conservation and management strategies. Italy, particularly the Ligurian region, harbors diverse freshwater fish communities and populations shaped by unique geological and hydrological conditions. Here, we investigated the suitability of environmental DNA (eDNA) metabarcoding to identify inter‐ and intraspecific diversity patterns of riverine fish populations in drainage basins on both sides of the main drainage divide (MDD) between the Adriatic and Ligurian river basins in Northern Italy. We collected 96 aquatic eDNA samples across 48 riverine sites, amplified them using a cytochrome b primer pair, and denoised the sequences to retrieve amplicon sequence variants (ASVs). We calculated communities' phylogenetic distance with betaMPD based on genetic distances derived from the ASVs, combined them with conductance‐based landscape metrics, and applied generalized dissimilarity models to assess spatial genetic structure. Our results reveal genetic differentiation among populations of several fish species, with some displaying clustering patterns across the main drainage divide and isolation by distance patterns. Overall, taxon richness was higher on the Adriatic side (12.82 ± 3.57; 26 unique taxa) than on the Ligurian side of the MDD (8.35 ± 3.66 SD, 25 unique taxa), but the other way around for ASV richness across all species (Ligurian side: 51.94 ± 25.79, 308 unique ASV, Adriatic side: 68.00 ± 32.97, 274 unique ASV). Our findings highlight the effectiveness of eDNA metabarcoding in uncovering various facets of diversity, shedding light on hidden genetic diversity among ASVs, and revealing significant spatial genetic structuring in freshwater fish populations across multiple species.

Arctic marine ecosystems are rapidly transforming due to climate change. Warming temperatures and shrinking sea ice are enabling boreal fish to expand northward, possibly disturbing cold‐adapted Arctic species assemblages. Species range shifts have been documented in the Bering and Barents Seas, raising concerns about ecosystem restructuring. Range shifts are especially difficult to detect in the Arctic due to sparse and inconsistent data. Here, we studied fish composition from eDNA water samples taken in East Greenland, Svalbard, the Barents Sea, and the Kara Sea during the TOPtoTOP and Arctic Century expeditions. We examined the environmental drivers of fish community structure using global dissimilarity models. We calculated the decadal rate of temperature change to identify the fastest‐changing areas. We compared fish detections from eDNA with published historical records for the Kara Sea to assess possible range expansions. We found that temperature was the main factor influencing the taxa turnover of fish communities, with Gadidae and Liparis sp. driving the greatest compositional differences. Over the past 30 years, temperatures increased by 0.2 to 0.6°C per decade at our study sites, with the highest increases in western Svalbard and the lowest in the eastern Kara Sea. Despite the apparent dependence on temperature, we identified only one species detected outside its known latitudinal range, and five species in the Kara Sea with recent occurrences or representing an extended distribution. Our study suggests that temperature, the main driver of fish community assembly, is increasing rapidly in the Arctic, and a few species have likely already shifted recently, or at least their detections are new in some areas. While these detections cannot be definitively linked to range shifts, our results highlight the need to improve monitoring of high‐latitude fish communities to detect and predict future ecosystem changes.

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.

Spatial and temporal monitoring of species threatened with extinction is of critical importance for conservation and ecosystem management. In the Mediterranean coast, the fan mussel (Pinna nobilis) is listed as critically endangered after suffering from a mass mortality event since 2016, leading to 100% mortality in most marine populations. Conventional monitoring for this macroinvertebrate is done using scuba, which is challenging in dense meadows or with low visibility. Here we developed an environmental DNA assay targeting the fan mussel and assessed the influence of several environmental parameters on the species detectability in situ. We developed and tested an eDNA molecular marker and collected 48 water samples in two sites at the Thau lagoon (France) with distinct fan mussel density, depths and during two seasons (summer and autumn). Our marker can amplify fan mussel DNA but lacks specificity since it also amplifies a conspecific species (Pinna rudis). We successfully amplified fan mussel DNA from in situ samples with 46 positive samples (out of 48) using ddPCR, although the DNA concentrations measured were low over almost all samples. Deeper sampling depth slightly increased DNA concentrations, but no seasonal effect was found. We highlight a putative spawning event on a single summer day with much higher DNA concentration compared to all other samples. We present an eDNA molecular assay able to detect the endangered fan mussel and provide guidelines to optimize the sampling protocol to maximize detectability. Effective and non‐invasive monitoring tools for endangered species are promising to monitor remaining populations and have the potential of ecological restoration or habitat recolonization following a mass mortality event.

Arctic fjords are experiencing rapid environmental shifts due to climate change, which may have significant impacts on marine biodiversity and ecosystem functioning. However, the impact of climate change on fjord biodiversity is difficult to quantify given the low accessibility and high cost to sample these areas. In this study, we sampled locations from inside to outside an Arctic fjord and used environmental DNA metabarcoding to assess how the biodiversity of fish and eukaryotic plankton communities relate with environmental conditions. We detected a total of 12 fish taxa and 872 MOTUs for eukaryotes and found marked gradients of temperature and salinity driven by the distance to the glacier at the terminal part of the fjord and depth. Eukaryotic richness was mainly associated with lower temperature and chlorophyll a. Moreover, co-inertia analyses showed a shared structure between eukaryotes assemblages and the environmental gradients. A partial association between eukaryote and fishes suggest a potential shared effect of environmental gradients among these taxa. Our findings provide a baseline for future studies to assess how these assemblages may be impacted by ongoing environmental changes and highlight how fjord ice loss and warming might shift environmental gradients and species distribution under climate change.

Aim: Coastal fishes have a fundamental role in marine ecosystem functioning and contributions to people, but face increasing threats due to climate change, habitat degradation and overexploitation. The extent to which human pressures are impacting coastal fish biodiversity in comparison with geographic and environmental factors at large spatial scale is still under scrutiny. Here, we took advantage of environmental DNA (eDNA) metabarcoding to investigate the relationship between fish biodiversity, including taxonomic and genetic components, and environmental but also socio‐economic factors.Location: Tropical, temperate and polar coastal areas.Time period: Present day.Major taxa studied: Marine fishes.Methods: We analysed fish eDNA in 263 stations (samples) in 68 sites distributed across polar, temperate and tropical regions. We modelled the effect of environmental, geographic and socio‐economic factors on α‐ and β‐diversity. We then computed the partial effect of each factor on several fish biodiversity components using taxonomic molecular units (MOTU) and genetic sequences. We also investigated the relationship between fish genetic α‐ and β‐diversity measured from our barcodes, and phylogenetic but also functional diversity.Results: We show that fish eDNA MOTU and sequence α‐ and β‐diversity have the strongest correlation with environmental factors on coastal ecosystems worldwide. However, our models also reveal a negative correlation between biodiversity and human dependence on marine ecosystems. In areas with high dependence, diversity of all fish, cryptobenthic fish and large fish MOTUs declined steeply. Finally, we show that a sequence diversity index, accounting for genetic distance between pairs of MOTUs, within and between communities, is a reliable proxy of phylogenetic and functional diversity.Main conclusions: Together, our results demonstrate that short eDNA sequences can be used to assess climate and direct human impacts on marine biodiversity at large scale in the Anthropocene and can further be extended to investigate biodiversity in its phylogenetic and functional dimensions.

To mitigate the ongoing threats to coastal ecosystems, and the biodiversity erosion they are causing, marine‐protected areas (MPAs) have emerged as powerful and widespread conservation tools. Strictly no‐take MPAs, also called marine reserves, undeniably promote fish biomass and density, but it remains unclear how biodiversity responds to protection. Identifying which facets of biodiversity respond to protection is critical for the management of MPAs and the development of relevant conservation strategies towards the achievement of biodiversity targets. We collected 99 environmental DNA (eDNA) samples inside and outside nine marine reserves in the Mediterranean Sea to assess the effect of protection on 11 biodiversity indicators based on fish traits, phylogeny and vulnerability to fishing. We controlled for the effect of environmental heterogeneity (habitat, bathymetry, productivity, temperature and accessibility) using a principal component analysis, and for spatial autocorrelation due to potential unmeasured factors. We found a positive and significant effect of protection on only 3 out of 11 indicators: functional and phylogenic diversity but also the ratio between demerso‐pelagic and benthic species richness. Rather, total fish richness responded significantly and negatively to protection. We did not detect any significant effect of protection on threatened and elasmobranch species richness, probably due to their large home range compared to the size of Mediterranean marine reserves. Synthesis and applications: Our findings highlight the importance of looking beyond the mere number of species to fully depict and understand the effect of marine reserves on biodiversity and evaluate the effectiveness of conservation measures. Rather, we propose a dashboard of three eDNA‐based indicators that can provide an early signal of ecosystem deterioration or recovery. eDNA metabarcoding offers a powerful tool to supply site‐specific and standardized taxonomic‐, phylogenetic‐ and trait‐based biodiversity assessments, in complement to other classical techniques, such as visual censuses or video surveys, able to estimate species abundance but also individual life‐stage and size.

Species inventories are the building blocks of our assessment of biodiversity patterns and human impact. Yet, historical inventories based on visual observations are often incomplete, impairing subsequent analyses of ecological mechanisms, extinction risk and management success. Environmental DNA (eDNA) metabarcoding is an emerging tool that can provide wider biodiversity assessments than classical visual‐based surveys. However, eDNA‐based inventories remain limited by sampling effort and reference database incompleteness. In this study, we propose a new framework coupling eDNA surveys and sampling‐theory methods to estimate species richness in under‐sampled and hyper‐diverse regions where some species remain absent from the checklist or undetected by visual surveys. We applied this framework to the coastal fish diversity in the heart of the coral triangle, the richest marine biodiversity hotspot worldwide. Combining data from 279 underwater visual censuses, 92 eDNA samples and an extensive custom genetic reference database, we show that eDNA metabarcoding recorded 196 putative species not detected by underwater visual census including 37 species absent from the regional checklist. We provide an updated checklist of marine fishes in the ‘Raja Ampat Bird's Head Peninsula' ecoregion with 2534 species including 1761 confirmed and 773 highly probable presences. The Chao lower‐bound diversity estimator, based on the incidence of rare species, shows that the region potentially hosts an additional 123 fish species, including pelagic, cryptobenthic and vulnerable species. The extended and hidden biodiversity along with their asymptotic estimates highlight the ability of eDNA to expand regional inventories and species distributions to better guide conservation strategies.

Quantifying fish species diversity in rich tropical marine environments remains challenging. Environmental DNA (eDNA) metabarcoding is a promising tool to face this challenge through the filtering, amplification, and sequencing of DNA traces from water samples. However, because eDNA concentration is low in marine environments, the reliability of eDNA to detect species diversity can be limited. Using an eDNA metabarcoding approach to identify fish Molecular Taxonomic Units (MOTUs) with a single 12S marker, we aimed to assess how the number of sampling replicates and filtered water volume affect biodiversity estimates. We used a paired sampling design of 30 L per replicate on 68 reef transects from 8 sites in 3 tropical regions. We quantified local and regional sampling variability by comparing MOTU richness, compositional turnover, and compositional nestedness. We found strong turnover of MOTUs between replicated pairs of samples undertaken in the same location, time, and conditions. Paired samples contained non‐overlapping assemblages rather than subsets of one another. As a result, non‐saturated localized diversity accumulation curves suggest that even 6 replicates (180 L) in the same location can underestimate local diversity (for an area <1 km). However, sampling regional diversity using ~25 replicates in variable locations (often covering 10 s of km) often saturated biodiversity accumulation curves. Our results demonstrate variability of diversity estimates possibly arising from heterogeneous distribution of eDNA in seawater, highly skewed frequencies of eDNA traces per MOTU, in addition to variability in eDNA processing. This high compositional variability has consequences for using eDNA to monitor temporal and spatial biodiversity changes in local assemblages. Avoiding false‐negative detections in future biomonitoring efforts requires increasing replicates or sampled water volume to better inform management of marine biodiversity using eDNA.