Metabarcoding of marine environmental DNA based on mitochondrial and nuclear genes

A belly full of jelly? DNA metabarcoding shows evidence for gelatinous zooplankton predation by several fish species in Greenland waters

The waters of Greenland harbour a high species richness and biomass of gelatinous zooplankton (GZP); however, their role in the diet of the many fish species, including commercially exploited species, has not yet been verified. Traditionally, GZP was considered to be a trophic dead end, i.e. with a limited contribution as prey for higher trophic levels. We applied DNA metabarcoding of two gene fragments (COI, 18S V1-V2) to the stomach contents of seven pelagic and demersal fish species in Greenland waters, to identify their prey composition as well as the occurrence of GZP predation. We detected GZP DNA reads in the stomachs of all investigated fish species, with frequency of occurrences ranging from 12.5% (for Melanogrammus aeglefinus) to 50% (for Argentina silus). GZP predation had not yet been reported for several of these species. GZP were found to majorly contribute to the diet of A. silus and Anarhichas denticulatus, particularly, the siphonophore Nanomia cara and the scyphozoan Atolla were of a high importance as prey, respectively. The use of multiple genetic markers enabled us to detect a total of 59 GZP taxa in the fish stomachs with several GZP species being detected only by one of the markers.

Mare-MAGE curated reference database of fish mitochondrial genes

Biodiversity assessment approaches based on molecular biology techniques such as metabarcoding, RAD-seq, or SnaPshot sequencing are increasingly applied in assessing marine and aquatic ecosystems. Here we present a new reference database for fish meta-barcoding based on mitochondrial genes. The Mare-MAGE database contains quality-checked sequences of the mitochondrial 12S ribosomal RNA and Cytochrome c Oxidase I gene. All sequences were obtained from the National Center for Biotechnology Information- GenBank (NBCI-GenBank), the European Nucleotide Archive (ENA), AquaGene Database and BOLD database, and have undergone intensive processing. They were checked for false annotations and non-target anomalies, according to the Integrated Taxonomic Information System (ITIS) and FishBase. The dataset is compiled in ARB-Home, FASTA and Qiime2 formats, and is publicly available from the Mare-MAGE database website ( http://mare-mage.weebly.com/ ). It includes altogether 231,333 COI and 12S rRNA gene sequences of fish, covering 19,506 species of 4,058 genera and 586 families.

Molecular based identification and phylogenetic relationship by using cytochrome b gene of Pangasius pangasius

Molecular appraoch for identification of unknown species by using Cytochrome b gene is an effective and reliable as compared with morphological based identification. For DNA barcoding universal molecular genes were used to identify the species. Cytochrome b is a specific gene used for identification purpose. DNA barcoding is a reliable and effective method compared to the different traditional morphological methods of specie identification. So,in the present study which was conducted to identify the species, a total of 50 fish samples were collected from five different sites. DNA was extracted by using the Phenol Chloroform method from muscle tissue. Five sequences were sequenced (one from each site), analyzed, and identified specific species as Pangasius pangasius. Identified sequences were variable in length from 369 bp (Site 1), 364 bp (Site 2), 364 bp (Site 3), 352 bp (Site 4), and 334 bp (Site 5). Identity matches on the NCBI database confirmed the specific specie as P. pangasius. A distancing tree was drawn to show maximum likelihood among the same and different species. Yet, in many cases fishes on diverse development stages are difficult to identify by morphological characters. DNA-based identification methods offer an analytically powerful addition or even an alternative tool for species identification and phylogenetic study. This work intends to provide an updated and extensive overview on the DNA based methods for fish species identification by using Cytochrome b gene as targeted markers for identification purpose.

Benthic invertebrates in Svalbard fjords-when metabarcoding does not outperform traditional biodiversity assessment

To protect and restore ecosystems and biodiversity is one of the 10 challenges identified by the United Nations's Decade of the Ocean Science. In this study we used eDNA from sediments collected in two fjords of the Svalbard archipelago and compared the taxonomic composition with traditional methods through metabarcoding, targeting mitochondrial CO1, to survey benthos. Clustering of 21.6 mill sequence reads with a d value of 13 in swarm, returned about 25 K OTU reads. An identification search with the BOLD database returned 12,000 taxonomy annotated sequences spanning a similarity range of 50% to 100%. Using an acceptance filter of minimum 90% similarity to the CO1 reference sequence, we found that 74% of the ca 100 taxon identified sequence reads were Polychaeta and 22% Nematoda. Relatively few other benthic invertebrate species were detected. Many of the identified sequence reads were extra-organismal DNA from terrestrial, planktonic, and photic zone sources. For the species rich Polychaeta, we found that, on average, only 20.6% of the species identified from morphology were also detected with DNA. This discrepancy was not due to missing reference sequences in the search database, because 90-100% (mean 96.7%) of the visually identified species at each station were represented with barcodes in Boldsystems. The volume of DNA samples is small compared with the volume searched in visual sorting, and the replicate DNA-samples in sum covered only about 2% of the surface area of a grab. This may considerably reduce the detection rate of species that are not uniformly distributed in the sediments. Along with PCR amplification bias and primer mismatch, this may be an important reason for the limited congruence of species identified with the two approaches. However, metabarcoding also identified 69 additional species that are usually overlooked in visual sample sorting, demonstrating how metabarcoding can complement traditional methodology by detecting additional, less conspicuous groups of organisms.

Using DNA Metabarcoding to Characterize the Prey Spectrum of Two Co-Occurring Themisto Amphipods in the Rapidly Changing Atlantic-Arctic Gateway Fram Strait

The two congeneric hyperiids Themisto libellula and T. abyssorum provide an important trophic link between lower and higher trophic levels in the rapidly changing Arctic marine ecosystem. These amphipods are characterized by distinct hydrographic affinities and are hence anticipated to be impacted differently by environmental changes, with major consequences for the Arctic food web. In this study, we applied DNA metabarcoding to the stomach contents of these Themisto species, to comprehensively reveal their prey spectra at an unprecedented-high-taxonomic-resolution and assess the regional variation in their diet across the Fram Strait. Both species feed on a wide variety of prey but their diet strongly differed in the investigated summer season, showing overlap for only a few prey taxa, such as calanoid copepods. The spatially structured prey field of T. libellula clearly differentiated it from T. abyssorum, of which the diet was mainly dominated by chaetognaths. Our approach also allowed the detection of previously overlooked prey in the diet of T. libellula, such as fish species and gelatinous zooplankton. We discuss the reasons for the differences in prey spectra and which consequences these may have in the light of ongoing environmental changes.

Environmental RNA as a Tool for Marine Community Biodiversity Assessments

Microscopic organisms are often overlooked in traditional diversity assessments due to the difficulty of identifying them based on morphology. Metabarcoding is a method for rapidly identifying organisms where Environmental DNA (eDNA) is used as a template. However, legacy DNA is problematically detected from organisms no longer in the environment during sampling. Environmental RNA (eRNA), which is only produced by living organisms, can also be collected from environmental samples and used for metabarcoding. The aim of this study was to determine differences in community composition and diversity between eRNA and eDNA templates for metabarcoding. Using mesocosms containing field-collected communities from an estuary, RNA and DNA were co-extracted from sediment, libraries were prepared for two loci (18S and COI), and sequenced using an Illumina MiSeq. Results show a higher number of unique sequences detected from eRNA in both markers and higher α-diversity compared to eDNA. Significant differences between eRNA and eDNA for all β-diversity metrics were also detected. This study is the first to demonstrate community differences detected with eRNA compared to eDNA from an estuarine system and illustrates the broad applications of eRNA as a tool for assessing benthic community diversity, particularly for environmental conservation and management applications.

Monitoring of environmental DNA from nonindigenous species of algae, dinoflagellates and animals in the North East Atlantic

Monitoring the distribution of marine nonindigenous species is a challenging task. To support this monitoring, we developed and validated the specificity of 12 primer-probe assays for detection of environmental DNA (eDNA) from marine species, all nonindigenous to Europe. The species include sturgeons, a Pacific red algae, oyster thief, a freshwater hydroid from the Black Sea, Chinese mitten crab, Pacific oyster, warty comb jelly, sand gaper, round goby, pink salmon, rainbow trout and North American mud crab. We tested all assays in the laboratory, on DNA extracted from both the target and non-target species to ensure that they only amplified DNA from the intended species. Subsequently, all assays were used to analyse water samples collected at 16 different harbours across two different seasons during 2017. We also included six previously published assays targeting eDNA from goldfish, European carp, two species of dinoflagellates of the genera Karenia and Prorocentrum, two species of the heterokont flagellate genus Pseudochattonella. Conventional monitoring was carried out alongside eDNA sampling but with only one sampling event over the one year. Because eDNA was relatively fast and easy to collect compared to conventional sampling, we sampled eDNA twice during 2017, which showed seasonal changes in the distribution of nonindigenous species. Comparing eDNA levels with salinity gradients did not show any correlation. A significant correlation was observed between number of species detected with conventional monitoring methods and number of species found using eDNA at each location. This supports the use of eDNA for surveillance of the distribution of marine nonindigenous species, where the speed and relative easy sampling in the field combined with fast molecular analysis may provide advantages compared to conventional monitoring methods. Prior validation of assays increases taxonomic precision, and laboratorial setup facilitates analysis of multiple samples simultaneously. The specific eDNA assays presented here can be implemented directly in monitoring programmes across Europe and potentially worldwide to infer a more precise picture of the dynamics in the distribution of marine nonindigenous species.

Hidden diversity of Ctenophora revealed by new mitochondrial COI primers and sequences

The mitochondrial gene cytochrome-c-oxidase subunit 1 (COI) is useful in many taxa for phylogenetics, population genetics, metabarcoding, and rapid species identifications. However, the phylum Ctenophora (comb jellies) has historically been difficult to study due to divergent mitochondrial sequences and the corresponding inability to amplify COI with degenerate and standard COI "barcoding" primers. As a result, there are very few COI sequences available for ctenophores, despite over 200 described species in the phylum. Here, we designed new primers and amplified the COI fragment from members of all major groups of ctenophores, including many undescribed species. Phylogenetic analyses of the resulting COI sequences revealed high diversity within many groups that was not evident from more conserved 18S rDNA sequences, in particular among the Lobata (Ctenophora; Tentaculata; Lobata). The COI phylogenetic results also revealed unexpected community structure within the genus Bolinopsis, suggested new species within the genus Bathocyroe, and supported the ecological and morphological differences of some species such as Lampocteis cruentiventer and similar undescribed lobates (Lampocteis sp. "V" stratified by depth, and "A" differentiated by colour). The newly designed primers reported herein provide important tools to enable researchers to illuminate the diversity of ctenophores worldwide via quick molecular identifications, improve the ability to analyse environmental DNA by improving reference libraries and amplifications, and enable a new breadth of population genetic studies.

Environmental conditions drive zooplankton community structure in the epipelagic oceanic water of the southern Gulf of Mexico: A molecular approach

Zooplankton plays a pivotal role in sustaining the majority of marine ecosystems. The distribution patterns and diversity of zooplankton provide key information for understanding the functioning of these ecosystems. Nevertheless, due to the numerous cryptic and sibling species and the lack of diagnostic characteristics for early developmental stages, the identification of the global-to-local patterns of zooplankton biodiversity and biogeography remains challenging in different research fields. The spatial and temporal changes in the zooplankton community in the open waters of the southern Gulf of Mexico were assessed using metabarcoding analysis of the V9 region of 18S rRNA and mitochondrial cytochrome oxidase c subunit I (COI). Additionally, a multiscale analysis was implemented to evaluate which environmental predictors may explain the variability in the structure of the zooplankton community. Our findings suggest that the synergistic effects of dissolved oxygen concentration, temperature, and longitude (intended as a proxy for still unidentified predictors) may explain both spatial and temporal zooplankton variability even with low contribution. Furthermore, the zooplankton distribution probably reflects the coexistence of three heterogeneous ecoregions and a bio-physical partitioning of the studied area. Finally, some taxa were either exclusive or predominant with either 18S or COI markers. This may suggest that comprehensive assessments of the zooplankton community may be more accurately met by the use of multilocus approaches.

Metabarcoding confirms the opportunistic foraging behaviour of Atlantic bluefin tuna and reveals the importance of gelatinous prey

Studies of the diet, feeding habits and trophic activity of top marine predators are essential for understanding their trophodynamics. The main direct method used for such studies thus far has been morphological inventories of stomach contents. This approach presents limitations such as missing gelatinous prey, which are usually digested too quickly to be detectable. Here, we analysed the stomachs of 48 Atlantic bluefin tuna (Thunnus thynnus, approximately 15 to 60 kg, including juveniles and adult fishes) collected from the Mediterranean Sea through the metabarcoding of two gene regions (cytochrome c oxidase subunit I (COI) and the ribosomal 18S-V1V2 region). The identified prey taxa and their relative read abundances (RRAs) estimated using COI results were in line with the findings of morphologically based inventories simultaneously performed on the same set of tuna samples. In both cases (and with the same rankings), the prey taxa included anchovy (Engraulis encrasicolus, here detected in more than 80% of samples, RRA = 43%), sardine (Sardina pilchardus, also approximately 80%, RRA = 30%), sprat (Sprattus sprattus, approximately 66%, RRA = 8%), mackerel (Scomber colias, approximately 44%, RRA = 7%) and cephalopods (approximately 15%, RRA = 1.4%). Another striking result was the detection, based on 18S (with which vertebrates were detected as the most abundant group, RRA = 61.6%), of a high prevalence and diversity of gelatinous organisms (RRA = 27.1%), including cnidarians (6.7%), salps (11.7%), and ctenophores (8.7%), the latter increasing with the size of the predator. These results thus support the hypothesis of the role of gelatinous prey in the diet of Atlantic bluefin tuna, suggesting that this species is even more generalist and opportunistic than previously thought. This study further confirms that DNA metabarcoding can be a powerful tool for assessing the diet and trophodynamics of top marine predators.

Suitability of a dual COI marker for marine zooplankton DNA metabarcoding

As DNA metabarcoding has become an emerging tool for surveying biodiversity, including its application in legally binding assessments, reliable and efficient barcodes are requested, especially for the highly diverse group of zooplankton. This study focuses on comparing the efficiency of two mitochondrial COI barcodes based on the internal primers mlCOIintF and mlCOIintR utilizing mesozooplankton samples collected in a Mediterranean lagoon. Our results indicate that after a slight adjustment, the mlCOIintR primer performs in combination with jdgLCO1490 (herein) very comparably to the much more widely used primer system mlCOIintF/jgHCO2198+dgHCO2198, in terms of level of taxonomic resolution, species detection and their relative abundance in terms of numbers of reads. As for some groups, like Ctenophora, this barcode is not suitable; a combination of them may be the best option to rely on the Folmer region in its entirety without the risk of losing information for a limited primer match.

High-throughput sequencing outperforms traditional morphological methods in Blue Catfish diet analysis and reveals novel insights into diet ecology

Blue Catfish Ictalurus furcatus are an invasive, yet economically important species in the Chesapeake Bay. However, their impact on the trophic ecology of this system is not well understood. In order to provide in-depth analysis of predation by Blue Catfish, we identified prey items using high-throughput DNA sequencing (HTS) of entire gastrointestinal tracts from 134 samples using two genetic markers, mitochondrial cytochrome c oxidase I (COI) and the nuclear 18S ribosomal RNA gene. We compared our HTS results to a more traditional "hybrid" approach that coupled morphological identification with DNA barcoding. The hybrid study was conducted on additional Blue Catfish samples (n = 617 stomachs) collected from the same location and season in the previous year. Taxonomic representation with HTS vastly surpassed that achieved with the hybrid methodology in Blue Catfish. Significantly, our HTS study identified several instances of at-risk and invasive species consumption not identified using the hybrid method, supporting the hypothesis that previous studies using morphological methods may greatly underestimate consumption of critical species. Finally, we report the novel finding that Blue Catfish diet diversity inversely correlates to daily flow rates, perhaps due to higher mobility and prey-seeking behaviors exhibited during lower flow.

eDNA metabarcoding of small plankton samples to detect fish larvae and their preys from Atlantic and Pacific waters

Zooplankton community inventories are the basis of fisheries management for containing fish larvae and their preys; however, the visual identification of early-stage larvae (the "missing biomass") is difficult and laborious. Here, eDNA metabarcoding was employed to detect zooplankton species of interest for fisheries from open and coastal waters. High-Throughput sequencing (HTS) from environmental samples using small water volumes has been proposed to detect species of interest whose DNA is the most abundant. We analyzed 6-L water samples taken from subtropical and tropical waters using Cytochrome oxidase I (COI) gene as metabarcode. In the open ocean, several commercial fish larvae and invertebrate species important in fish diet were found from metabarcodes and confirmed from individual barcoding. Comparing Atlantic, Mediterranean, Red Sea, and Pacific samples we found a lower taxonomic depth of OTU assignments in samples from tropical waters than in those from temperate ones, suggesting large gaps in reference databases for those areas; thus a higher effort of zooplankton barcoding in tropical oceans is highly recommended. This and similar simplified sampling protocols could be applied in early detection of species important for fisheries.

Integrating eDNA metabarcoding and simultaneous underwater visual surveys to describe complex fish communities in a marine biodiversity hotspot

Marine biodiversity can be surveyed using underwater visual censuses and recently with eDNA metabarcoding. Although a promising tool, eDNA studies have shown contrasting results related to its detection scale and the number of species identified compared to other survey methods. Also, its accuracy relies on complete reference databases used for taxonomic assignment and, as other survey methods, species detection may show false-negative and false-positive errors. Here, we compared results from underwater visual censuses and simultaneous eDNA metabarcoding fish surveys in terms of observed species and community composition. We also assess the effect of a custom reference database in the taxonomic assignment, and evaluate occupancy, capture and detection probabilities, as well as error rates of eDNA survey data. We amplified a 12S rRNA fish barcode from 24 sampling sites in the gulf of California. More species were detected with eDNA metabarcoding than with UVC. Because each survey method largely detected different sets of species, the combined approach doubled the number of species registered. Both survey methods recovered a known biodiversity gradient and a biogeographic break, but eDNA captured diversity over a broader geographic and bathymetric scale. Furthermore, the use of a modest-sized custom reference database significantly increased taxonomic assignment. In a subset of species, occupancy models revealed eDNA surveys provided similar or higher detection probabilities compared to UVC. The occupancy value of each species had a large influence on eDNA detectability, and in the false positive and negative error. Overall, these results highlight the potential of eDNA metabarcoding in complementing other established ecological methods for studies of marine fishes.

Marine health of the Arabian Gulf: Drivers of pollution and assessment approaches focusing on desalination activities

The Arabian Gulf is one of the most adversely affected marine environments worldwide, which results from combined pollution drivers including climate change, oil and gas activities, and coastal anthropogenic disturbances. Desalination activities are one of the major marine pollution drivers regionally and internationally. Arabian Gulf countries represent a hotspot of desalination activities as they are responsible for nearly 50% of the global desalination capacity. Building desalination plants, up-taking seawater, and discharging untreated brine back into the sea adversely affects the biodiversity of the marine ecosystems. The present review attempted to reveal the potential negative effects of desalination plants on the Gulf's marine environments. We emphasised different conventional and innovative assessment tools used to assess the health of marine environments and evaluate the damage exerted by desalination activity in the Gulf. Finally, we suggested effective management approaches to tackle the issue including the significance of national regulations and regional cooperation.

Early detection of marine non-indigenous species on Svalbard by DNA metabarcoding of sediment

Non-indigenous species (NIS) in the Arctic have an increased likelihood of arrival from ship traffic in the region, while the survival potential of the species becomes more likely in a warming environment. Monitoring is essential to detect the rate and magnitude of the establishment of NIS. In this study, a list of 123 potential marine NIS for Svalbard was drafted and the presence of marine NIS in soft sediment of Kongsfjorden in Svalbard was assessed using molecular metabarcoding techniques. For 37 species, including eight potential Arctic NIS, we generated new 18S and/or COI barcode sequences to improve the available online reference databases. In total, 299 species were identified in the sediment samples, including seven potential NIS. Three of these potential NIS have not been reported before in Svalbard: the harpacticoid copepod Euterpina acutifrons, and the ascidians Botrylloides violaceus and Molgula manhattensis. Another novel observation for Svalbard was the polychaete Chone mollis. Additional studies are needed to assess whether the NIS have been established on Svalbard and what their potential impact on the local system may be. Metabarcoding proved to be an effective monitoring tool to detect the presence of new species in Svalbard marine waters. We advise its use to set up a baseline record for the presence of NIS at points of entry, especially harbours. This approach is also valuable for biodiversity monitoring, in particular the detection of small organisms and life stages that are hard to identify using current visual techniques.

A novel metabarcoding primer pair for environmental DNA analysis of Cephalopoda (Mollusca) targeting the nuclear 18S rRNA region

Cephalopods are pivotal components of marine food webs, but biodiversity studies are hampered by challenges to sample these agile marine molluscs. Metabarcoding of environmental DNA (eDNA) is a potentially powerful technique to study oceanic cephalopod biodiversity and distribution but has not been applied thus far. We present a novel universal primer pair for metabarcoding cephalopods from eDNA, Ceph18S (Forward: 5'-CGC GGC GCT ACA TAT TAG AC-3', Reverse: 5'-GCA CTT AAC CGA CCG TCG AC-3'). The primer pair targets the hypervariable region V2 of the nuclear 18S rRNA gene and amplifies a relatively short target sequence of approximately 200 bp in order to allow the amplification of degraded DNA. In silico tests on a reference database and empirical tests on DNA extracts from cephalopod tissue estimate that 44-66% of cephalopod species, corresponding to about 310-460 species, can be amplified and identified with this primer pair. A multi-marker approach with the novel Ceph18S and two previously published cephalopod mitochondrial 16S rRNA primer sets targeting the same region (Jarman et al. 2006 Mol. Ecol. Notes. 6, 268-271; Peters et al. 2015 Mar. Ecol. 36, 1428-1439) is estimated to amplify and identify 89% of all cephalopod species, of which an estimated 19% can only be identified by Ceph18S. All sequences obtained with Ceph18S were submitted to GenBank, resulting in new 18S rRNA sequences for 13 cephalopod taxa.

Field application of an improved protocol for environmental DNA extraction, purification, and measurement using Sterivex filter

Environmental DNA (eDNA) is increasingly popular as a useful non-invasive method to monitor and study biodiversity and community structure in freshwater and marine environments. To effectively extract eDNA from the filter surface is a fundamental factor determining the representativeness of the samples. We improved the eDNA extraction efficiency of an established Sterivex method by 12- to 16-fold using a larger volume of lysis buffer mix coupled with backflushing the cartridges. The DNeasy extraction column could be overloaded when the environmental sample input is high, possibly due to a higher nonspecific binding present in environmental samples, thus resulting in a relatively lower quantity measured. Therefore, we included an internal control DNA in the extraction to monitor the extraction and purification efficiencies in field samples, which is crucial for quantification of original eDNA concentration. The use of Takara Probe qPCR Mix supplemented with protein-based additives improved the robustness of the real time PCR assay on inhibitor-rich environmental samples, but prior purification by Qiagen PowerClean Pro Cleanup kit could be essential for inhibitor-rich water samples, even though the recovery rate was unexpectedly low (average 33.0%). The improved extraction and quantification complement the qualitative analyses including metabarcoding and metagenomics in field application.

The crossover from microscopy to genes in marine diversity: from species to assemblages in marine pelagic copepods

An accurate identification of species and communities is a prerequisite for analysing and recording biodiversity and community shifts. In the context of marine biodiversity conservation and management, this review outlines past, present and forward-looking perspectives on identifying and recording planktonic diversity by illustrating the transition from traditional species identification based on morphological diagnostic characters to full molecular genetic identification of marine assemblages. In this process, the article presents the methodological advancements by discussing progress and critical aspects of the crossover from traditional to novel and future molecular genetic identifications and it outlines the advantages of integrative approaches using the strengths of both morphological and molecular techniques to identify species and assemblages. We demonstrate this process of identifying and recording marine biodiversity on pelagic copepods as model taxon. Copepods are known for their high taxonomic and ecological diversity and comprise a huge variety of behaviours, forms and life histories, making them a highly interesting and well-studied group in terms of biodiversity and ecosystem functioning. Furthermore, their short life cycles and rapid responses to changing environments make them good indicators and core research components for ecosystem health and status in the light of environmental change. This article is part of the theme issue 'Integrative research perspectives on marine conservation'.

Comparing fish prey diversity for a critically endangered aquatic mammal in a reserve and the wild using eDNA metabarcoding

Using environmental DNA (eDNA) metabarcoding, we compared fish diversity in two distinct water bodies within the Yangtze River Basin with known populations of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis; YFP): the Tian-e-Zhou Reserve and Poyang Lake. We aimed to create a fish surveying tool for use in the Yangtze River Basin, while also gaining a better understanding of the prey distribution and diversity within two of the remaining strongholds of YFP. 16S rRNA universal primers were developed to amplify fish eDNA. After high-throughput sequencing and stringent data filtering, we identified a total of 75 fish species (6 orders, 9 families, 57 genera) across seasons and regions. Nine of the 75 fish species were among the 28 known YFP prey species, three of which were detected in all water samples. Our eDNA metabarcoding identified many species that had been previously captured using traditional netting practices, but also numerous species not previously collected in these water bodies. Fish diversity was higher in Poyang Lake than in Tian-e-Zhou Reserve, as well as higher in the spring than in summer. These methods provide a broadly applicable tool to quantify fish diversity and distributions throughout the Yangtze River Basin, and to inform conservation strategies of YFP.

Marine water environmental DNA metabarcoding provides a comprehensive fish diversity assessment and reveals spatial patterns in a large oceanic area

Current methods for monitoring marine fish (including bony fishes and elasmobranchs) diversity mostly rely on trawling surveys, which are invasive, costly, and time-consuming. Moreover, these methods are selective, targeting a subset of species at the time, and can be inaccessible to certain areas. Here, we used environmental DNA (eDNA), the DNA present in the water column as part of shed cells, tissues, or mucus, to provide comprehensive information about fish diversity in a large marine area. Further, eDNA results were compared to the fish diversity obtained in pelagic trawls. A total of 44 5 L-water samples were collected onboard a wide-scale oceanographic survey covering about 120,000 square kilometers in Northeast Atlantic Ocean. A short region of the 12S rRNA gene was amplified and sequenced through metabarcoding generating almost 3.5 million quality-filtered reads. Trawl and eDNA samples resulted in the same most abundant species (European anchovy, European pilchard, Atlantic mackerel, and blue whiting), but eDNA metabarcoding resulted in more detected bony fish and elasmobranch species (116) than trawling (16). Although an overall correlation between fishes biomass and number of reads was observed, some species deviated from the common trend, which could be explained by inherent biases of each of the methods. Species distribution patterns inferred from eDNA metabarcoding data coincided with current ecological knowledge of the species, suggesting that eDNA has the potential to draw sound ecological conclusions that can contribute to fish surveillance programs. Our results support eDNA metabarcoding for broad-scale marine fish diversity monitoring in the context of Directives such as the Common Fisheries Policy or the Marine Strategy Framework Directive.

From metabarcoding to metaphylogeography: separating the wheat from the chaff

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.

Biodiversity assessment of tropical shelf eukaryotic communities via pelagic eDNA metabarcoding

Our understanding of marine communities and their functions in an ecosystem relies on the ability to detect and monitor species distributions and abundances. Currently, the use of environmental DNA (eDNA) metabarcoding is increasingly being applied for the rapid assessment and monitoring of aquatic species. Most eDNA metabarcoding studies have either focussed on the simultaneous identification of a few specific taxa/groups or have been limited in geographical scope. Here, we employed eDNA metabarcoding to compare beta diversity patterns of complex pelagic marine communities in tropical coastal shelf habitats spanning the whole Caribbean Sea. We screened 68 water samples using a universal eukaryotic COI barcode region and detected highly diverse communities, which varied significantly among locations, and proved good descriptors of habitat type and environmental conditions. Less than 15% of eukaryotic taxa were assigned to metazoans, most DNA sequences belonged to a variety of planktonic "protists," with over 50% of taxa unassigned at the phylum level, suggesting that the sampled communities host an astonishing amount of micro-eukaryotic diversity yet undescribed or absent from COI reference databases. Although such a predominance of micro-eukaryotes severely reduces the efficiency of universal COI markers to investigate vertebrate and other metazoans from aqueous eDNA, the study contributes to the advancement of rapid biomonitoring methods and brings us closer to a full inventory of extant marine biodiversity.

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.