Marine reserves are being established worldwide to conserve and manage marine resources, and are also often assumed to conserve the evolutionary potential of marine populations. However, comparisons of genomic patterns inside and outside marine reserves in multiple species are scarce. Here, we aim to fill this gap by 1) comparing genomic variation in protected versus unprotected areas in a network of eight reserves in the Mediterranean Sea, and 2) disentangling the effect of seascape, space and marine reserves on this genomic variation in 1297 individuals from four species genotyped at more than 10 000 SNPs each. We report homogeneous patterns of local genetic diversity within and outside reserves in three of the four species (the white seabream Diplodus sargus , the striped red mullet Mullus surmuletus and the European spiny lobster Palinurus elephas ). Surprisingly, the comber Serranus cabrilla shows lower genetic diversity in protected areas, highlighting an apparent conservation paradox that may be due to increased predation by protected fishes in the reserves. Multivariate analyses reveal that seascape factors (salinity, temperature and chlorophyll) and space exert a significant influence on genomic variation in some cases, while protection status has no detectable effect. Nevertheless, four marine reserves (Cabo de Palos, Cerbère‐Banyuls, Cap de Creus and Illes Columbretes) harbor singular genetic variation in one or two species. This pattern, observed in the four species, suggests that reserves contribute to preserve genetic variation locally. Our study illustrates that we can document key aspects of the genetic variation of marine species and their interaction with seascape factors, space and protection status in multiple species. It stresses the need for a multi‐species approach to inform marine conservation planning, opening up new perspectives at the community level.

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 alpha-and beta-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 alpha-and beta-diversity measured from our barcodes, and phylogenetic but also functional diversity. Results: We show that fish eDNA MOTU and sequence alpha-and beta-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.

Environmental DNA is mainly not only used at the interspecific level, to quantify species diversity in ecosystems, but can also be used to quantify intraspecific genetic variability, thus avoiding the need to sample individual tissue. However, errors in the amplification and sequencing of eDNA samples can blur this intraspecific signal and strongly over-estimate genetic diversity. Existing bioinformatics pipelines therefore need to be tested to evaluate whether reliable levels of intraspecific genetic variability can be derived from eDNA samples. Here, we compare the ability of twelve metabarcoding pipelines to detect intraspecific genetic variability combining five programs. All pipelines have common pre-processing steps, a processing data step using programs among obiclean ; DADA2; SWARM; and LULU. An additional chimera removal step is also investigated based on two programs (VSEARCH or DADA2). The case study was the natural intraspecific variation within Mullus surmuletus in experimental settings. We developed specific primers for this species, located on the mitochondrial D-loop fragment (barcode MS-DL06). Thirty-nine individuals were collected from the Mediterranean Sea, placed into four aquariums, and their DNA was sequenced on this marker to build an intraspecific reference database. After filtering the aquarium water, DNA was extracted, amplified, and sequenced using the primer pair developed. We then quantified the number of true haplotypes returned by each pipeline and its capacity to eliminate most of the erroneous sequences. We show that the program DADA2 with a two-parent chimeric sequence removal step is the best tool to estimate intraspecific diversity from eDNA. Furthermore, our approach was also able to detect true M. surmuletus haplotypes in two eDNA samples collected in the Mediterranean Sea. We conclude that the combination of an appropriate intrapopulation barcode and a denoising pipeline like DADA2 with a chimeric sequence removal step is promising to make population-level inference using environmental DNA possible.

Increasing speed and magnitude of global change threaten the world's biodiversity and particularly coral reef fishes. A better understanding of large-scale patterns and processes on coral reefs is essential to prevent fish biodiversity decline but it requires new monitoring approaches. Here, we use environmental DNA metabarcoding to reconstruct well-known patterns of fish biodiversity on coral reefs and uncover hidden patterns on these highly diverse and threatened ecosystems. We analysed 226 environmental DNA (eDNA) seawater samples from 100 stations in five tropical regions (Caribbean, Central and Southwest Pacific, Coral Triangle and Western Indian Ocean) and compared those to 2047 underwater visual censuses from the Reef Life Survey in 1224 stations. Environmental DNA reveals a higher (16%) fish biodiversity, with 2650 taxa, and 25% more families than underwater visual surveys. By identifying more pelagic, reef-associated and crypto-benthic species, eDNA offers a fresh view on assembly rules across spatial scales. Nevertheless, the reef life survey identified more species than eDNA in 47 shared families, which can be due to incomplete sequence assignment, possibly combined with incomplete detection in the environment, for some species. Combining eDNA metabarcoding and extensive visual census offers novel insights on the spatial organization of the richest marine ecosystems.

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.

Assessing the impact of global changes and protection effectiveness is a key step in monitoring marine fishes. Most traditional census methods are demanding or destructive. Nondisturbing and nonlethal approaches based on video and environmental DNA are alternatives to underwater visual census or fishing. However, their ability to detect multiple biodiversity factors beyond traditional taxonomic diversity is still unknown. For bony fishes and elasmobranchs, we compared the performance of eDNA metabarcoding and long‐term remote video to assess species’ phylogenetic and functional diversity. We used 10 eDNA samples from 30 L of water each and 25 hr of underwater videos over 4 days on Malpelo Island (pacific coast of Colombia), a remote marine protected area. Metabarcoding of eDNA detected 66% more molecular operational taxonomic units (MOTUs) than species on video. We found 66 and 43 functional entities with a single eDNA marker and videos, respectively, and higher functional richness for eDNA than videos. Despite gaps in genetic reference databases, eDNA also detected a higher fish phylogenetic diversity than videos; accumulation curves showed how 1 eDNA transect detected as much phylogenetic diversity as 25 hr of video. Environmental DNA metabarcoding can be used to affordably, efficiently, and accurately census biodiversity factors in marine systems. Although taxonomic assignments are still limited by species coverage in genetic reference databases, use of MOTUs highlights the potential of eDNA metabarcoding once reference databases have expanded.

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F‐measure, root‐mean‐square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F‐measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study.

Coral reefs host the highest fish diversity on Earth despite covering less than 0.1% of the ocean’s seafloor. At the same time they are also extremely threatened. Data syntheses over decades of surveys estimate the total number of coral reef fishes to vary from 2,400 to 8,000 species distributed among roughly 100 families. But this diversity remains largely unknown. Here, we investigated how environmental DNA (eDNA) could describe the distribution of fish diversity in coral reefs. We generated 504,457,267 raw 12S ribosomal DNA (rDNA) sequence reads from 251 samples (2,693 PCR replicates) collected at 25 sites in 145 stations covering five regions across the Indian, Pacific and Atlantic Oceans. Bioinformatic analysis clustered these sequences into 2,160 molecular operational taxonomic units (MOTUs) corresponding to distinct species (Marques 2020) We compared our results, with visual census surveys from Reef Life Survey, on 2,813 transects in tropical regions. Our outcomes demonstrate the capacity of eDNA metabarcoding from water samples to reconstruct well-known biogeographic patterns of fish diversity on coral reefs, such as species richness gradients towards the coral triangle, and family proportion stability across sites (Bellwood and Hughes 2001). Additionally, eDNA survey data documented a higher fish species (16%) and family (50%) diversity than estimates obtained with underwater visual surveys carried out on 20 times more sites. MOTU richness per family retrieved with eDNA closely matched fish species richness within families recorded in visual census data. However, eDNA revealed higher richness of reef-associated species and species from adjacent habitats, of cryptobenthic or nocturnal species, but also of pelagic and wide-ranging species. eDNA survey data showed that fish diversity is characterized by spatially heterogeneous species assemblages among regions, with more dissimilarity among adjacent coral reefs than detected with visual survey data. Unlike visual surveys, eDNA metabarcoding revealed the same prevalence of rarity as expected under the neutral theory of biodiversity, suggesting the predominance of random processes and ecological equivalence within trophic groups at large scale to explain fish biodiversity patterns on coral reefs. Our study demonstrates how sequencing eDNA from water provides a rapid and effective approach to characterize and assess coral reef diversity across large spatial scales, thereby also uncovering hidden biodiversity patterns.

How far do marine larvae disperse in the ocean? Decades of population genetic studies have revealed generally low levels of genetic structure at large spatial scales (hundreds of kilometres). Yet this result, typically based on discrete sampling designs, does not necessarily imply extensive dispersal. Here, we adopt a continuous sampling strategy along 950 km of coast in the northwestern Mediterranean Sea to address this question in four species. In line with expectations, we observe weak genetic structure at a large spatial scale. Nevertheless, our continuous sampling strategy uncovers a pattern of isolation by distance at small spatial scales (few tens of kilometres) in two species. Individual-based simulations indicate that this signal is an expected signature of restricted dispersal. At the other extreme of the connectivity spectrum, two pairs of individuals that are closely related genetically were found more than 290 km apart, indicating long-distance dispersal. Such a combination of restricted dispersal with rare long-distance dispersal events is supported by a high-resolution biophysical model of larval dispersal in the study area, and we posit that it may be common in marine species. Our results bridge population genetic studies with direct dispersal studies and have implications for the design of marine reserve networks.

Extensive fishing has led to fish stock declines throughout the last decades. While clear stock identification is required for designing management schemes, stock delineation is problematic due to generally low levels of genetic structure in marine species. The development of genomic resources can help to solve this issue. Here, we present the first mitochondrial and nuclear draft genome assemblies of three economically important Mediterranean fishes, the white seabream, the striped red mullet, and the comber. The assemblies are between 613 and 785 Mbp long and contain between 27,222 and 32,375 predicted genes. They were used as references to map Restriction-site Associated DNA markers, which were developed with a single-digest approach. This approach provided between 15,710 and 21,101 Single Nucleotide Polymorphism markers per species. These genomic resources will allow uncovering subtle genetic structure, identifying stocks, assigning catches to populations and assessing connectivity. Furthermore, the annotated genomes will help to characterize adaptive divergence.

Genetic diversity is estimated to be declining faster than species diversity under escalating threats, but its spatial distribution remains poorly documented at the global scale. Theory predicts that similar processes should foster congruent spatial patterns of genetic and species diversity, but empirical studies are scarce. Using a mined database of 50,588 georeferenced mitochondrial DNA barcode sequences (COI) for 3,815 marine and 1,611 freshwater fish species respectively, we examined the correlation between genetic diversity and species diversity and their global distributions in relation to climate and geography. Genetic diversity showed a clear spatial organisation, but a weak association with species diversity for both marine and freshwater species. We found a predominantly positive relationship between genetic diversity and sea surface temperature for marine species. Genetic diversity of freshwater species varied primarily across the regional basins and was negatively correlated with average river slope. The detection of genetic diversity patterns suggests that conservation measures should consider mismatching spatial signals across multiple facets of biodiversity.

Human activities impact all ecosystems on Earth, which urges scientists to better understand biodiversity changes across temporal and spatial scales. Environmental DNA (eDNA) metabarcoding is a promising non‐invasive method to assess species composition in a wide range of ecosystems. Yet, this method requires the completeness of a reference database, i.e. a list of DNA sequences attached to each species of the regional pool, which is rarely met. As an alternative, molecular operational taxonomic units (MOTUs) can be extracted as clusters of sequences. However, the extent to which the diversity of MOTUs can predict the diversity of species across spatial scales is unknown. Here, we used 196 samples along the Rhone river (France) for which the reference database is complete to assess whether a blind eDNA approach can reliably predict the ground‐truth number of species at different spatial scales. Using the 12S rDNA teleo primer, we curated and clustered 60 million sequences into MOTUs using a new assembled bioinformatic pipeline. We show that stringent quality filters were necessary to remove artefact noise, notably MOTUs present in a single PCR replicate, which represented 55% of MOTUs (103). Post‐clustering cleaning also removed 19 additional erroneous MOTUs and only discarded one truly present species. We then show that the diversity of retained fish MOTUs accurately predicted the local (α, r = 0.98) and regional (γ) ground‐truth species diversity (67 MOTUs versus 63 species), but also the species dissimilarity between samples (β‐diversity, r = 0.98). This work paves the way towards extending the use of eDNA metabarcoding in community ecology and biogeography despite major gaps in genetic reference databases.

Genome-environment association methods aim to detect genetic markers associated with environmental variables. The detected associations are usually analysed separately to identify the genomic regions involved in local adaptation. However, a recent study suggests that single-locus associations can be combined and used in a predictive way to estimate environmental variables for new individuals on the basis of their genotypes. Here, we introduce an original approach to predict the environmental range (values and upper and lower limits) of species genotypes from the genetic markers significantly associated with those environmental variables in an independent set of individuals. We illustrate this approach to predict aridity in a database constituted of 950 individuals of wild beets and 299 individuals of cultivated beets genotyped at 14,409 random single nucleotide polymorphisms (SNPs). We detected 66 alleles associated with aridity and used them to calculate the fraction (I) of aridity-associated alleles in each individual. The fraction I correctly predicted the values of aridity in an independent validation set of wild individuals and was then used to predict aridity in the 299 cultivated individuals. Wild individuals had higher median values and a wider range of values of aridity than the cultivated individuals, suggesting that wild individuals have higher ability to resist to stress-aridity conditions and could be used to improve the resistance of cultivated varieties to aridity.

Protein structure prediction based on comparative modeling is the most efficient way to produce structural models when it can be performed. ORION is a dedicated webserver based on a new strategy that performs this task. The identification by ORION of suitable templates is performed using an original profile-profile approach that combines sequence and structure evolution information. Structure evolution information is encoded into profiles using structural features, such as solvent accessibility and local conformation-with Protein Blocks-, which give an accurate description of the local protein structure. ORION has recently been improved, increasing by 5% the quality of its results. The ORION web server accepts a single protein sequence as input and searches homologous protein structures within minutes. Various databases such as PDB, SCOP and HOMSTRAD can be mined to find an appropriate structural template. For the modeling step, a protein 3D structure can be directly obtained from the selected template by MODELLER and displayed with global and local quality model estimation measures. The sequence and the predicted structure of 4 examples from the CAMEO server and a recent CASP11 target from the 'Hard' category (T0818-D1) are shown as pertinent examples. Our web server is accessible at http://www.dsimb.inserm.fr/ORION/.