Biodiversity assessment of marine hard-bottom communities is hindered by the high diversity and size-ranges of the organisms present. We developed a DNA metabarcoding protocol for biodiversity characterization of structurally complex natural marine hard-bottom communities. We used two molecular markers: the “Leray fragment” of mitochondrial cytochrome c oxidase (COI), for which a novel primer set was developed, and the V7 region of the nuclear small subunit ribosomal RNA (18S). Eight different shallow marine littoral communities from two National Parks in Spain (one in the Atlantic Ocean and another in the Mediterranean Sea) were studied. Samples were sieved into three size fractions from where DNA was extracted separately. Bayesian clustering was used for delimiting molecular operational taxonomic units (MOTUs) and custom reference databases were constructed for taxonomic assignment. Despite applying stringent filters, we found high values for MOTU richness (2,510 and 9,679 MOTUs with 18S and COI, respectively), suggesting that these communities host a large amount of yet undescribed eukaryotic biodiversity. Significant gaps are still found in sequence reference databases, which currently prevent the complete taxonomic assignment of the detected sequences. In our dataset, 85% of 18S MOTUs and 64% of COI MOTUs could be identified to phylum or lower taxonomic level. Nevertheless, those unassigned were mostly rare MOTUs with low numbers of reads, and assigned MOTUs comprised over 90% of the total sequence reads. The identification rate might be significantly improved in the future, as reference databases are further completed. Our results show that marine metabarcoding, currently applied mostly to plankton or sediments, can be adapted to structurally complex hard bottom samples. Thus, eukaryotic metabarcoding emerges as a robust, fast, objective and affordable method to comprehensively characterize the diversity of marine benthic communities dominated by macroscopic seaweeds and colonial or modular sessile metazoans. The 18S marker lacks species-level resolution and thus cannot be recommended to assess the detailed taxonomic composition of these communities. Our new universal primers for COI can potentially be used for biodiversity assessment with high taxonomic resolution in a wide array of marine, terrestrial or freshwater eukaryotic communities.

Metabarcoding extra‐organismal DNA from environmental samples is now a key technique in aquatic biomonitoring and ecosystem health assessment. Of critical consideration when designing experiments, and especially so when developing community standards and legislative frameworks, is the choice of genetic marker and primer set. Mitochondrial cytochrome c oxidase subunit I (COI), the standard DNA barcode marker for animals, with its extensive reference library, taxonomic discriminatory power and predictable sequence variation, is the natural choice for many metabarcoding applications. However, for targeting specific taxonomic groups in environmental samples, the utility of COI has yet to be fully scrutinized. Here, by using a case study of marine and freshwater fishes from the British Isles, we quantify the in silico performance of twelve primer pairs from four mitochondrial loci – COI, cytochrome b, 12S and 16S – in terms of reference library coverage, taxonomic discriminatory power and primer universality. We subsequently test in vitro four primer pairs – three COI and one 12S – for their specificity, reproducibility, and congruence with independent datasets derived from traditional survey methods at five estuarine and coastal sites around the English Channel and North Sea. Our results show that for aqueous extra‐organismal DNA at low template concentrations, both metazoan‐targeted and fish‐targeted COI primers perform poorly in comparison to 12S, exhibiting low levels of reproducibility due to non‐specific amplification of prokaryotic and non‐target eukaryotic DNAs. An ideal metabarcode would have an extensive reference library upon which custom primers could be designed, either for broad assessments of biodiversity, or taxon specific surveys. Such a database is available for COI, but low primer specificity hinders practical application, while conversely, 12S primers offer high specificity, but lack adequate references. The latter, however, can be mitigated by expanding the concept of DNA barcodes to include whole mitochondrial genomes generated by genome‐skimming existing tissue collections.

As environmental DNA (eDNA) becomes an increasingly valuable resource for marine ecosystem monitoring, understanding variation in its persistence across contrasting environments is critical. Here, we quantify the breakdown of macrobial eDNA over a spatio-temporal axis of locally extreme conditions, varying from ocean-influenced offshore to urban-inshore, and between winter and summer. We report that eDNA degrades 1.6 times faster in the inshore environment than the offshore environment, but contrary to expectation we find no difference over season. Analysis of environmental covariables show a spatial gradient of salinity and a temporal gradient of pH, with salinity—or the biotic correlates thereof—most important. Based on our estimated inshore eDNA half-life and naturally occurring eDNA concentrations, we estimate that eDNA may be detected for around 48 h, offering potential to collect ecological community data of high local fidelity. We conclude by placing these results in the context of previously published eDNA decay rates.

Environmental DNA reveals unsuspected shark diversity and calls for monitoring and protection of residual populations.

Marine sediments are home to one of the richest species pools on Earth, but logistics and a dearth of taxonomic work-force hinders the knowledge of their biodiversity. We characterized α- and β-diversity of deep-sea assemblages from submarine canyons in the western Mediterranean using an environmental DNA metabarcoding. We used a new primer set targeting a short eukaryotic 18S sequence (ca. 110 bp). We applied a protocol designed to obtain extractions enriched in extracellular DNA from replicated sediment corers. With this strategy we captured information from DNA (local or deposited from the water column) that persists adsorbed to inorganic particles and buffered short-term spatial and temporal heterogeneity. We analysed replicated samples from 20 localities including 2 deep-sea canyons, 1 shallower canal, and two open slopes (depth range 100–2,250 m). We identified 1,629 MOTUs, among which the dominant groups were Metazoa (with representatives of 19 phyla), Alveolata, Stramenopiles, and Rhizaria. There was a marked small-scale heterogeneity as shown by differences in replicates within corers and within localities. The spatial variability between canyons was significant, as was the depth component in one of the canyons where it was tested. Likewise, the composition of the first layer (1 cm) of sediment was significantly different from deeper layers. We found that qualitative (presence-absence) and quantitative (relative number of reads) data showed consistent trends of differentiation between samples and geographic areas. The subset of exclusively benthic MOTUs showed similar patterns of β-diversity and community structure as the whole dataset. Separate analyses of the main metazoan phyla (in number of MOTUs) showed some differences in distribution attributable to different lifestyles. Our results highlight the differentiation that can be found even between geographically close assemblages, and sets the ground for future monitoring and conservation efforts on these bottoms of ecological and economic importance.

Metabarcoding is by now a well‐established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α‐ and β‐diversity estimates, that is, interspecies (or inter‐MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra‐MOTU) information—mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra‐MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co‐occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two‐step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra‐ and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.

Environmental DNA (eDNA) metabarcoding has revolutionized biomonitoring in both marine and freshwater ecosystems. However, for semi‐aquatic and terrestrial animals, the application of this technique remains relatively untested. We first assess the efficiency of eDNA metabarcoding in detecting semi‐aquatic and terrestrial mammals in natural lotic ecosystems in the UK by comparing sequence data recovered from water and sediment samples to the mammalian communities expected from historical data. Secondly, using occupancy modelling we compared the detection efficiency of eDNA metabarcoding to multiple conventional non‐invasive survey methods (latrine surveys and camera trapping). eDNA metabarcoding detected a large proportion of the expected mammalian community within each area. Common species in the areas were detected at the majority of sites. Several key species of conservation concern in the UK were detected by eDNA sampling in areas where authenticated records do not currently exist, but potential false positives were also identified. Water‐based eDNA metabarcoding provided comparable results to conventional survey methods in per unit of survey effort for three species (water vole, field vole and red deer) using occupancy models. The comparison between survey ‘effort’ to reach a detection probability of ≥.95 revealed that 3–6 water replicates would be equivalent to 3–5 latrine surveys and 5–30 weeks of single camera deployment, depending on the species. Synthesis and applications. eDNA metabarcoding can be used to generate an initial ‘distribution map’ of mammalian diversity at the landscape level. If conducted during times of peak abundance, carefully chosen sampling points along multiple river courses provide a reliable snapshot of the species that are present in a catchment area. In order to fully capture solitary, rare and invasive species, we would currently recommend the use of eDNA metabarcoding alongside other non‐invasive surveying methods (i.e. camera traps) to maximize monitoring efforts.

Sharks are charismatic predators that play a key role in most marine food webs. Their demonstrated vulnerability to exploitation has recently turned them into flagship species in ocean conservation. Yet, the assessment and monitoring of the distribution and abundance of such mobile species in marine environments remain challenging, often invasive and resource-intensive. Here we pilot a novel, rapid and non-invasive environmental DNA (eDNA) metabarcoding approach specifically targeted to infer shark presence, diversity and eDNA read abundance in tropical habitats. We identified at least 21 shark species, from both Caribbean and Pacific Coral Sea water samples, whose geographical patterns of diversity and read abundance coincide with geographical differences in levels of anthropogenic pressure and conservation effort. We demonstrate that eDNA metabarcoding can be effectively employed to study shark diversity. Further developments in this field have the potential to drastically enhance our ability to assess and monitor elusive oceanic predators, and lead to improved conservation strategies.