Evaluating detection limits of next-generation sequencing for the surveillance and monitoring of international marine pests

Development and Testing of Species-Specific Primers for Detecting the Presence of the Northern Pacific Sea Star (Asterias amurensis) from Environmental DNA

The starfish Asterias amurensis, a well-known predator of molluscan species in intertidal ecosystems, has caused substantial ecological and economic losses in North China such as offshore Qingdao. Effective monitoring and prevention measures are urged to minimize its negative impacts. Compared with traditional biomonitoring methods, environmental DNA technology has emerged as a powerful and cost-efficient tool for inferring species' presence and abundance. In this study, we developed a pair of species-specific primers (i.e., Ast-F and Ast-R) for the A. amurensis mitochondrial COI gene and tested its utility in amplifying and quantifying the DNA fragments from environmental samples under both laboratory and field conditions. The results of controlled water tank experiments demonstrated that the amount of eDNA released by A. amurensis was positively related to its biomass; after the removal of the starfish, the eDNA degraded significantly in 24 h and remained detectable for 8 days. The number of eDNA copies enriched tended to increase with smaller pore size of filter membrane and larger volume of filtered water. For field tests, we confirmed the validation of our approach in six locations in Qingdao by filtering 1000 ml water per sample with a 0.45-µm pore size filtration. All the amplification products generated a single and bright band via gel electrophoresis, and the quantitative PCR results unveiled significant differences in eDNA copies. This study provided an eDNA-based approach for investigating the distribution and biomass of A. amurensis, which may help to formulate early warning and management strategies in coastal Qingdao and other regions.

How to identify biofouling species in marine and freshwater

The identification and management of biofouling remain pressing challenges in marine and freshwater ecosystems, with significant implications for environmental sustainability and industrial operations. This comprehensive review synthesizes the current state-of-the-art in biofouling identification technologies, examining eight prominent methodologies: Microscopy Examination, Molecular Biology, Remote Sensing, Community Involvement, Ecological Methods, Artificial Intelligence, Chemical Analysis, and Macro Photography. Each method is evaluated for its respective advantages and disadvantages, considering factors such as precision, scalability, cost, and data quality. Furthermore, the review identifies current obstacles that inhibit the optimal utilization of these technologies, ranging from technical limitations and high operational costs to issues of data inconsistency and subjectivity. Finally, the review posits a future outlook, advocating for the development of integrated, standardized systems that amalgamate the strengths of individual approaches. Such advancement will pave the way for more effective and sustainable strategies for biofouling identification and management.

Assessing the Seasonal and Spatial Dynamics of Zooplankton through DNA Metabarcoding in a Temperate Estuary

Zooplankton are key components of estuarine trophic networks. However, routine monitoring is hindered by the difficulty of morphology-based identification. DNA-based methods allow us to circumvent some of these hurdles, providing precise species identifications regardless of the taxonomic expertise of the investigator or the developmental stage of the specimens. However, the process is dependent on the completeness of the reference libraries. In this study, we sought to evaluate the potential of DNA metabarcoding to assess the seasonal (summer, autumn, and early spring) and spatial dynamics of zooplankton (four locations spanning ca. 6 km) in the Lima estuary (NW Portugal). Two genetic markers were used: the cytochrome c oxidase subunit I and the V4 hypervariable region of the ribosomal 18S rRNA genes. Overall, 327 species were recovered, and both markers displayed minute overlap (7% were detected with both markers). Species richness, composition, and taxonomic distinctness were majorly influenced by the season, with a declining tendency from summer (highest number of exclusive species, n = 74) to spring. Second to season, the taxa composition was influenced by spatial variation where the most downstream site displayed the highest number of exclusive species, n = 53. A total of 16 non-indigenous species were detected using metabarcoding, but only one (Austrominus modestus) has been documented out in the estuary. In conclusion, both the seasonal and spatial gradients influenced the recovered richness, composition, and taxonomic distinctness, confirming the great aptitude of DNA metabarcoding for providing higher density monitoring and shedding new light on the composition and dynamics of complex zooplankton communities.

Chromosome-level genome assembly of the northern Pacific seastar Asterias amurensis

Asterias amurensis has attracted widespread concern because of its population outbreaks, which has impacted fisheries and aquaculture, as well as disrupting local ecosystems. A high-quality reference genome is necessary to better investigate mechanisms of outbreak and adaptive changes. Combining PacBio HiFi and Hi-C sequencing data, we generated a chromosome-level A. amurensis genome with a size of 491.53 Mb. The contig N50 and scaffold N50 were 8.05 and 23.75 Mb, respectively. The result of BUSCO analysis revealed a completeness score of 98.85%. A total of 16,531 protein-coding genes were predicted in the genome, of which 94.63% were functionally annotated. The high-quality genome assembly resulting from this study will provide a valuable genetic resource for future research on the mechanism of population outbreaks and invasion ecology.

An integrative protocol for one-step PCR amplicon library construction and accurate demultiplexing of pooled sequencing data

High-throughput sequencing of amplicons has been widely used to precisely and efficiently identify species compositions and analyze community structures, greatly promoting biological studies involving large amounts of complex samples, especially those involving environmental and pathogen-monitoring ones. Commercial library preparation kits for amplicon sequencing, which generally require multiple steps, including adapter ligation and indexing, are expensive and time-consuming, especially for applications at a large scale. To overcome these limitations, a "one-step PCR approach" has been previously proposed for constructions of amplicon libraries using long fusion primers. However, efficient amplifications of target genes and accurate demultiplexing of pooled sequencing data remain to be addressed. To tackle these, we present an integrative protocol for one-step PCR amplicon library construction (OSPALC). High-quality reads have been generated by this approach to reliably identify species compositions of mock bacterial communities and environmental samples. With this protocol, the amplicon library is constructed through one regular PCR with long primers, and the total cost per DNA/cDNA sample decreases to just 7% of the typical cost via the multi-step PCR approach. Empirically tested primers and optimized PCR conditions to construct OSPALC libraries for 16S rDNA V4 regions are demonstrated as a case study. Tools to design primers targeting at any genomic regions are also presented. In principle, OSPALC can be readily applied to construct amplicon libraries of any target genes using DNA or RNA samples, and will facilitate research in numerous fields.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00182-1.

Harnessing environmental DNA to reveal biogeographical patterns of non-indigenous species for improved co-governance of the marine environment in Aotearoa New Zealand

Aotearoa New Zealand's Northern region is a major gateway for the incursion and establishment of non-indigenous species (NIS) populations due to high numbers of recreational and commercial vessels. This region also holds a unique marine ecosystem, home to many taonga (treasured) species of cultural and economic importance. Regular surveillance, eradication plans and public information sharing are undertaken by local communities and governmental organizations to protect these ecosystems from the impact of NIS. Recently, considerable investments went into environmental DNA (eDNA) research, a promising approach for the early detection of NIS for complementing existing biosecurity systems. We applied eDNA metabarcoding for elucidating bioregional patterns of NIS distributions across a gradient from harbors (NIS hotspots) to open seas (spreading areas). Samples were collected during a research cruise sailing across three Aotearoa New Zealand harbors, Waitematā, Whangārei and Pēwhairangi (Bay of Islands), and their adjacent coastal waters. The small-ribosomal subunit (18S rRNA) and mitochondrial cytochrome c oxidase I (COI) genes were screened using the online Pest Alert Tool for automated detection of putative NIS sequences. Using a probabilistic modelling approach, location-dependent occupancies of NIS were investigated and related to the current information on species distribution from biosecurity surveillance programs. This study was collaboratively designed with Māori partners to initiate a model of co-governance within the existing science system.

Plumbing the depths with environmental DNA (eDNA): Metabarcoding reveals biodiversity zonation at 45-60 m on mesophotic coral reefs

Mesophotic coral ecosystems (MCEs) are tropical reefs found at depths of ~30-150 m, below the region most heavily impacted by heat stress and other disturbances. Hence, MCEs may serve as potential refugia for threatened shallow reefs, but they also harbour depth-endemic fauna distinct from shallow reefs. Previous studies have characterized biodiversity patterns along depth gradients, but focussed primarily on conspicuous taxa (fishes, corals, etc.). Environmental DNA (eDNA) metabarcoding offers a more holistic approach to assess biodiversity patterns across the tree of life. Here, we use three metabarcoding assays targeting fishes (16S rRNA), eukaryotes (18S rDNA) and metazoans (COI) to assess biodiversity change from the surface to ~90 m depth across 15-m intervals at three sites within the Hawaiian Archipelago. We observed significant community differences between most depth zones, with distinct zonation centred at 45-60 m for eukaryotes and metazoans, but not for fishes. This finding may be attributable to the higher mobility of reef fishes, although methodological limitations are likely a contributing factor. The possibility for MCEs to serve as refugia is not excluded for fishes, but invertebrate communities >45 m are distinct, indicating limited connectivity for the majority of reef fauna. This study provides a new approach for surveying biodiversity on MCEs, revealing patterns in a much broader context than the limited-taxon studies that comprise the bulk of our present knowledge.

Gut content metabarcoding of specialized feeders is not a replacement for environmental DNA assays of seawater in reef environments

In tropical marine ecosystems, the coral-based diet of benthic-feeding reef fishes provides a window into the composition and health of coral reefs. In this study, for the first time, we compare multi-assay metabarcoding sequences of environmental DNA (eDNA) isolated from seawater and partially digested gut items from an obligate corallivore butterflyfish (Chaetodon lunulatus) resident to coral reef sites in the South China Sea. We specifically tested the proportional and statistical overlap of the different approaches (seawater vs gut content metabarcoding) in characterizing eukaryotic community composition on coral reefs. Based on 18S and ITS2 sequence data, which differed in their taxonomic sensitivity, we found that gut content detections were only partially representative of the eukaryotic communities detected in the seawater based on low levels of taxonomic overlap (3 to 21%) and significant differences between the sampling approaches. Overall, our results indicate that dietary metabarcoding of specialized feeders can be complimentary to, but is no replacement for, more comprehensive environmental DNA assays of reef environments that might include the processing of different substrates (seawater, sediment, plankton) or traditional observational surveys. These molecular assays, in tandem, might be best suited to highly productive but cryptic oceanic environments (kelp forests, seagrass meadows) that contain an abundance of organisms that are often small, epiphytic, symbiotic, or cryptic.

eDNA metabarcoding reveals shifts in sediment eukaryote communities in a metal contaminated estuary

Metal contamination is a global issue impacting biodiversity in urbanised estuaries. Traditional methods to assess biodiversity are time consuming, costly and often exclude small or cryptic organisms due to difficulties with morphological identification. Metabarcoding approaches have been increasingly recognised for their utility in monitoring, however studies have focused on freshwater and marine systems despite the ecological significance of estuaries. We targeted estuarine eukaryote communities within the sediments of Australia's largest urbanised estuary, where a history of industrial activity has resulted in a metal contamination gradient. We identified specific eukaryote families with significant correlations with bioavailable metal concentrations, indicating sensitivity or tolerance to specific metals. While polychaete families Terebellidae and Syllidae demonstrated tolerance to the contamination gradient, members of the meio- and microfaunal communities including diatoms, dinoflagellates and nematodes displayed sensitivities. These may have high value as indicators but are frequently missed in traditional surveys due to sampling limitations.

Deciphering the diversity and distribution of chromophytic phytoplankton in the Bohai Sea and the Yellow Sea via RuBisCO genes (rbcL)

Present study investigated composition and distribution of chromophytic phytoplankton in the Bohai Sea (BS) and the Yellow Sea (YS) by using rbcL genes. Bacillariophyceae, Haptophyceae and Pelagophyceae were the most abundant phytoplankton groups. Distinct phytoplankton communities were observed in the BS and the YS: offshore stations were dominated by bloom forming genera Thalassiosira and Skeletonema, while brown tide-forming species including Chrysochromulina spp. and Aureococcus anophagefferens were commonly found in the nearshore areas. Redundancy analysis showed that phosphate, temperature and silicic acid play key roles in structuring chromophytic phytoplankton, such as phytoplankton at nearshore stations were affected by nutrient runoff from adjacent rivers (Yellow River). Anthropogenic activities in the Bohai Sea and seasonal circulation of ocean currents may also contribute to shaping chromophytic phytoplankton communities. This study provides data support and foundational observations of chromophytic phytoplankton in the BS and the YS, and their responses to environmental gradients and human activities.

Skimming for barcodes: rapid production of mitochondrial genome and nuclear ribosomal repeat reference markers through shallow shotgun sequencing

DNA barcoding is critical to conservation and biodiversity research, yet public reference databases are incomplete. Existing barcode databases are biased toward cytochrome oxidase subunit I (COI) and frequently lack associated voucher specimens or geospatial metadata, which can hinder reliable species assignments. The emergence of metabarcoding approaches such as environmental DNA (eDNA) has necessitated multiple marker techniques combined with barcode reference databases backed by voucher specimens. Reference barcodes have traditionally been generated by Sanger sequencing, however sequencing multiple markers is costly for large numbers of specimens, requires multiple separate PCR reactions, and limits resulting sequences to targeted regions. High-throughput sequencing techniques such as genome skimming enable assembly of complete mitogenomes, which contain the most commonly used barcoding loci (e.g., COI, 12S, 16S), as well as nuclear ribosomal repeat regions (e.g., ITS1&2, 18S). We evaluated the feasibility of genome skimming to generate barcode references databases for marine fishes by assembling complete mitogenomes and nuclear ribosomal repeats. We tested genome skimming across a taxonomically diverse selection of 12 marine fish species from the collections of the National Museum of Natural History, Smithsonian Institution. We generated two sequencing libraries per species to test the impact of shearing method (enzymatic or mechanical), extraction method (kit-based or automated), and input DNA concentration. We produced complete mitogenomes for all non-chondrichthyans (11/12 species) and assembled nuclear ribosomal repeats (18S-ITS1-5.8S-ITS2-28S) for all taxa. The quality and completeness of mitogenome assemblies was not impacted by shearing method, extraction method or input DNA concentration. Our results reaffirm that genome skimming is an efficient and (at scale) cost-effective method to generate all mitochondrial and common nuclear DNA barcoding loci for multiple species simultaneously, which has great potential to scale for future projects and facilitate completing barcode reference databases for marine fishes.

The use of environmental DNA to monitor impacted coastal estuaries

Environmental DNA (eDNA) metabarcoding is increasingly being used to assess community composition in coastal ecosystems. In this study, we chose to examine temporal and spatial changes in the aquatic community of Manly Lagoon - one of the most heavily developed and polluted estuaries in eastern Australia. Based on metabarcoding of the 16S mitochondrial gene (for fish) and the 18S nuclear gene (for macroinvertebrates), we identified seasonal differences in fish and macroinvertebrate community composition as well as species richness, which correlated, in some cases, with the environmental parameters of sea surface temperature and freshwater input. Moreover, given the greater taxonomic resolution of fish versus macroinvertebrate assignments, we identified several known migratory fish species of management importance that contributed significantly to the overall patterns observed. Overall, our data support the use of eDNA metabarcoding to track fish assemblages shifting in response to environmental drivers in polluted estuaries with increased sampling and consultation with historical data.

Sea Snacks from DNA Tracks: Using DNA Metabarcoding to Characterize The Diet of Green Turtles (Chelonia Mydas)

The green turtle (Chelonia mydas) is a circumglobal species with a wide dietary breadth that varies among regions and life history stages. Comprehensive understanding of foraging ecology over space and time is critical to inform conservation and management of this species and its habitats. Here, we used DNA metabarcoding to test candidate primer sets with 39 gut content homogenates from stranded green turtles (FL, USA) to identify primer sets that maximize detection of food items and specificity of taxonomic classifications. We tested six existing universal primer sets to detect plants, animals, and eukaryotes more broadly (CO1, 18SV1-V3, 18SV4, rbcL, UPA, ITS). The CO1 and 18SV4 primer sets produced the greatest number of dietary amplicon sequence variants (ASVs) and unique taxonomic classifications, and they were the only primer sets to amplify taxa from all three kingdoms relevant to green turtle diet (Animalia, Chromista, and Plantae). Even though the majority of CO1-derived reads were of host origin (>90%), this primer set still produced the largest number of dietary ASVs classified to species among the six primer sets. However, because the CO1 primer set failed to detect both vascular plants and green algae, we do not recommend the use of this primer set on its own to characterize green turtle diet. Instead, our findings support previous research highlighting the utility of using multiple primer sets, specifically targeting CO1 and the V4 region of the 18S gene, as doing so will provide the most comprehensive understanding of green turtle diet. More generally, our results highlight the importance of primer and loci selection and the need to validate primer sets against the study system of interest. The addition of DNA metabarcoding with optimized primer sets to the sea turtle researcher's toolbox will both increase our understanding of foraging ecology and better inform science-based conservation and ecosystem management.

Strategies for sample labelling and library preparation in DNA metabarcoding studies

Metabarcoding of DNA extracted from environmental or bulk specimen samples is increasingly used to profile biota in basic and applied biodiversity research because of its targeted nature that allows sequencing of genetic markers from many samples in parallel. To achieve this, PCR amplification is carried out with primers designed to target a taxonomically informative marker within a taxonomic group, and sample-specific nucleotide identifiers are added to the amplicons prior to sequencing. The latter enables assignment of the sequences back to the samples they originated from. Nucleotide identifiers can be added during the metabarcoding PCR and during "library preparation", that is, when amplicons are prepared for sequencing. Different strategies to achieve this labelling exist. All have advantages, challenges and limitations, some of which can lead to misleading results, and in the worst case compromise the fidelity of the metabarcoding data. Given the range of questions addressed using metabarcoding, ensuring that data generation is robust and fit for the chosen purpose is critically important for practitioners seeking to employ metabarcoding for biodiversity assessments. Here, we present an overview of the three main workflows for sample-specific labelling and library preparation in metabarcoding studies on Illumina sequencing platforms; one-step PCR, two-step PCR, and tagged PCR. Further, we distill the key considerations for researchers seeking to select an appropriate metabarcoding strategy for their specific study. Ultimately, by gaining insights into the consequences of different metabarcoding workflows, we hope to further consolidate the power of metabarcoding as a tool to assess biodiversity across a range of applications.

Shifts in DNA yield and biological community composition in stored sediment: implications for paleogenomic studies

Lake sediments hold a wealth of information from past environments that is highly valuable for paleolimnological reconstructions. These studies increasingly apply modern molecular tools targeting sedimentary DNA (sedDNA). However, sediment core sampling can be logistically difficult, making immediate subsampling for sedDNA challenging. Sediment cores are often refrigerated (4 °C) for weeks or months before subsampling. We investigated the impact of storage time on changes in DNA (purified or as cell lysate) concentrations and shifts in biological communities following storage of lake surface sediment at 4 °C for up to 24 weeks. Sediment samples (~ 0.22 g, in triplicate per time point) were spiked with purified DNA (100 or 200 ng) or lysate from a brackish water cyanobacterium that produces the cyanotoxin nodularin or non-spiked. Samples were analysed every 1–4 weeks over a 24-week period. Droplet digital PCR showed no significant decrease in the target gene (nodularin synthetase – subunit F; ndaF) over the 24-week period for samples spiked with purified DNA, while copy number decreased by more than half in cell lysate-spiked samples. There was significant change over time in bacteria and eukaryotic community composition assessed using metabarcoding. Amongst bacteria, the cyanobacterial signal became negligible after 5 weeks while Proteobacteria increased. In the eukaryotic community, Cercozoa became dominant after 6 weeks. These data demonstrate that DNA yields and community composition data shift significantly when sediments are stored chilled for more than 5 weeks. This highlights the need for rapid subsampling and appropriate storage of sediment core samples for paleogenomic studies.

Environmental DNA sampling protocols for the surveillance of marine non-indigenous species in Irish coastal waters

Understanding the spread and distribution of Non-Indigenous Species (NIS) is key when implementing legislation to maintain good ecosystem health. Environmental DNA (eDNA) has shown great potential to detect aquatic organisms in a rapid and cost-effective way, however their applicability to new environments must be validated prior to their implementation. Here, we tested different field sampling methods in combination with eDNA metabarcoding to develop a tool to detect NIS. Large and small volumes of seawater were filtered, in addition to the collection of sediment and horizontal tow net samples at 12 locations across four distinct geographic areas in Ireland. The biggest dissimilarity in the species recovered was found between sediment and town net samples. Tow nets showed to be the most efficient. A total of 357 taxa were identified, including 16 NIS. Fine mesh tow nets were identified as the most cost-efficient for large-scale monitoring and surveillance of NIS.

Pathogens co-transported with invasive non-native aquatic species: implications for risk analysis and legislation

Invasive Non-Native Species (INNS) can co-transport externally and internally other organisms including viruses, bacteria and other eukaryotes (including metazoan parasites), collectively referred to as the symbiome. These symbiotic organisms include pathogens, a small minority of which are subject to surveillance and regulatory control, but most of which are currently unscrutinized and/or unknown. These putatively pathogenetic symbionts can potentially pose diverse risks to other species, with implications for increased epidemiological risk to agriculture and aquaculture, wildlife/ecosystems, and human health (zoonotic diseases). The risks and impacts arising from co-transported known pathogens and other symbionts of unknown pathogenic virulence, remain largely unexplored, unlegislated, and difficult to identify and quantify. Here, we propose a workflow using PubMed and Google Scholar to systematically search existing literature to determine any known and potential pathogens of aquatic INNS. This workflow acts as a prerequisite for assessing the nature and risk posed by co-transported pathogens of INNS; of which a better understanding is necessary to inform policy and INNS risk assessments. Addressing this evidence gap will be instrumental to devise an appropriate set of statutory responsibilities with respect to these symbionts, and to underpin new and more effective legislative processes relating to the disease screening and risk assessment of INNS.

Pathogens co-transported with invasive non-native aquatic species: implications for risk analysis and legislation

Invasive Non-Native Species (INNS) can co-transport externally and internally other organisms including viruses, bacteria and other eukaryotes (including metazoan parasites), collectively referred to as the symbiome. These symbiotic organisms include pathogens, a small minority of which are subject to surveillance and regulatory control, but most of which are currently unscrutinized and/or unknown. These putatively pathogenetic symbionts can potentially pose diverse risks to other species, with implications for increased epidemiological risk to agriculture and aquaculture, wildlife/ecosystems, and human health (zoonotic diseases). The risks and impacts arising from co-transported known pathogens and other symbionts of unknown pathogenic virulence, remain largely unexplored, unlegislated, and difficult to identify and quantify. Here, we propose a workflow using PubMed and Google Scholar to systematically search existing literature to determine any known and potential pathogens of aquatic INNS. This workflow acts as a prerequisite for assessing the nature and risk posed by co-transported pathogens of INNS; of which a better understanding is necessary to inform policy and INNS risk assessments. Addressing this evidence gap will be instrumental to devise an appropriate set of statutory responsibilities with respect to these symbionts, and to underpin new and more effective legislative processes relating to the disease screening and risk assessment of INNS.

Development of droplet digital Polymerase Chain Reaction assays for the detection of long-finned (Anguilla dieffenbachii) and short-finned (Anguilla australis) eels in environmental samples

Freshwater eels are ecologically, and culturally important worldwide. The New Zealand long-finned eel (Anguilla dieffenbachii) and short-finned eel (Anguilla australis) are apex predators, playing an important role in ecosystem functioning of rivers and lakes. Recently, there has been a national decline in their populations due to habitat destruction and commercial harvest. The emergence of targeted environmental DNA detection methodologies provides an opportunity to enhance information about their past and present distributions. In this study we successfully developed species-specific droplet digital Polymerase Chain Reaction (ddPCR) assays to detect A. dieffenbachii and A. australis DNA in water and sediment samples. Assays utilized primers and probes designed for regions of the mitochondrial cytochrome b and 16S ribosomal RNA genes in A. dieffenbachii and A. australis, respectively. River water samples (n = 27) were analyzed using metabarcoding of fish taxa and were compared with the ddPCR assays. The presence of A. dieffenbachii and A. australis DNA was detected in a greater number of water samples using ddPCR in comparison to metabarcoding. There was a strong and positive correlation between gene copies (ddPCR analyses) and relative eel sequence reads (metabarcoding analyses) when compared to eel biomass. These ddPCR assays provide a new method for assessing spatial distributions of A. dieffenbachii and A. australis in a range of environments and sample types.

Assessing the Use of Molecular Barcoding and qPCR for Investigating the Ecology of Prorocentrum minimum (Dinophyceae), a Harmful Algal Species

Prorocentrum minimum is a species of marine dinoflagellate that occurs worldwide and can be responsible for harmful algal blooms (HABs). Some studies have reported it to produce tetrodotoxin; however, results have been inconsistent. qPCR and molecular barcoding (amplicon sequencing) using high-throughput sequencing have been increasingly applied to quantify HAB species for ecological analyses and monitoring. Here, we isolated a strain of P. minimum from eastern Australian waters, where it commonly occurs, and developed and validated a qPCR assay for this species based on a region of ITS rRNA in relation to abundance estimates from the cultured strain as determined using light microscopy. We used this tool to quantify and examine ecological drivers of P. minimum in Botany Bay, an estuary in southeast Australia, for over ~14 months in 2016–2017. We compared abundance estimates using qPCR with those obtained using molecular barcoding based on an 18S rRNA amplicon. There was a significant correlation between the abundance estimates from amplicon sequencing and qPCR, but the estimates from light microscopy were not significantly correlated, likely due to the counting method applied. Using amplicon sequencing, ~600 unique actual sequence variants (ASVs) were found, much larger than the known phytoplankton diversity from this region. P. minimum abundance in Botany Bay was found to be significantly associated with lower salinities and higher dissolved CO₂ levels.

Comparison of Larval Fish Detections Using Morphology-Based Taxonomy versus High-Throughput Sequencing for Invasive Species Early Detection

When first introduced, invasive species typically evade detection; DNA barcoding coupled with high-throughput sequencing (HTS) may be more sensitive and accurate than morphology-based taxonomy, and thereby improve invasive (or rare) species detection. We quantified the relative error of species detection between morphology-based and HTS-based taxonomic identification of ichthyoplankton collections from the Port of Duluth, Minnesota, an aquatic non-native species introduction 'hot-spot' in the Laurentian Great Lakes. We found HTS-based taxonomy identified 28 species and morphology-based taxonomy 30 species, of which 27 were common to both. Among samples, 76% of family-level taxonomic assignments agreed; however, only 42% of species assignments agreed. Most errors were attributed to morphology-based taxonomy, whereas HTS-based taxonomy error was low. For this study system, for most non-native fishes, the detection probability by randomized survey for larvae was similar to that by a survey that is optimized for non-native species early detection of juveniles and adults. We conclude that classifying taxonomic errors by comparing HTS results against morphology-based taxonomy is an important step toward incorporating HTS-based taxonomy into biodiversity surveys.

Status and prospects of marine NIS detection and monitoring through (e)DNA metabarcoding

In coastal ecosystems, non-indigenous species (NIS) are recognized as a major threat to biodiversity, ecosystem functioning and socio-economic activities. Here we present a systematic review on the use of metabarcoding for NIS surveillance in marine and coastal ecosystems, through the analysis of 42 publications. Metabarcoding has been mainly applied to environmental DNA (eDNA) from water samples, but also to DNA extracted from bulk organismal samples. DNA extraction kits have been widely used and the 18S rRNA and the COI genes the most employed markers, but less than half of the studies targeted more than one marker loci. The Illumina MiSeq platform has been used in >50% of the publications. Current weaknesses include potential occurrence of false negatives due to the primer-biased or faulty DNA amplification and the incompleteness of reference libraries. This is particularly concerning in the case of NIS surveillance, where proficiency in species level detection is critical. Until these weaknesses are resolved, ideally NIS metabarcoding should be supported by complementary approaches, such as morphological analysis or more targeted molecular approaches (e.g. qPCR, ddPCR). Even so, metabarcoding has already proved to be a highly sensitive tool to detect small organisms or undifferentiated life stages across a wide taxonomic range. In addition, it also seems to be very effective in ballast water management and to improve the spatial and temporal sampling frequency of NIS surveillance in marine and coastal ecosystems. Although specific protocols may be required for species-specific NIS detection, for general monitoring it would be vital to settle on a standard protocol able to generate comparable results among surveillance campaigns and regions of the globe, seeking the best approach for detecting the broadest range of species, while minimizing the chances of a false positive or negative detection.

Biosecurity implications of drifting marine plastic debris: Current knowledge and future research

The introduction and spread of marine non-indigenous species (NIS) and pathogens into new habitats are a major threat to biodiversity, ecosystem services, human health, and can have substantial economic consequences. Shipping is considered the main vector for marine biological invasions; less well understood is the increased spread of marine NIS and pathogens rafting on marine plastic debris (MPD). Despite an increasing research interest and recent progress in characterizing the plastisphere, this manuscript highlights critical knowledge gaps and research priorities towards a better understanding of the biosecurity implications of MPD. We advocate for future research to (i) investigate plastisphere community succession and the factors influencing NIS propagules and pathogens recruitment through robust experimental investigations; (ii) combine microscopy and molecular approaches to effectively assess the presence of specific taxa; (iii) include additional genetic markers to thoroughly characterize the biodiversity associated with MPD and explore the presence of specific marine pests.

Under the karst: detecting hidden subterranean assemblages using eDNA metabarcoding in the caves of Christmas Island, Australia

Subterranean ecosystems are understudied and challenging to conventionally survey given the inaccessibility of underground voids and networks. In this study, we conducted a eukaryotic environmental DNA (eDNA) metabarcoding survey across the karst landscape of Christmas Island, (Indian Ocean, Australia) to evaluate the utility of this non-invasive technique to detect subterranean aquatic 'stygofauna' assemblages. Three metabarcoding assays targeting the mitochondrial 16S rRNA and nuclear 18S genes were applied to 159 water and sediment samples collected from 23 caves and springs across the island. Taken together, our assays detected a wide diversity of chordates, cnidarians, porifera, arthropods, molluscs, annelids and bryozoans from 71 families across 60 orders. We report a high level of variation between cave and spring subterranean community compositions which are significantly influenced by varying levels of salinity. Additionally, we show that dissolved oxygen and longitudinal gradients significantly affect biotic assemblages within cave communities. Lastly, we combined eDNA-derived community composition and environmental (water quality) data to predict potential underground interconnectivity across Christmas Island. We identified three cave and spring groups that showed a high degree of biotic and abiotic similarity indicating likely local connectivity. This study demonstrates the applicability of eDNA metabarcoding to detect subterranean eukaryotic communities and explore underground interconnectivity.

Translational Molecular Ecology in practice: Linking DNA-based methods to actionable marine environmental management

Molecular-based approaches can provide timely biodiversity assessments, showing an immense potential to facilitate decision-making in marine environmental management. However, the uptake of molecular data into environmental policy remains minimal. Here, we showcase a selection of local to global scale studies applying molecular-based methodologies for environmental management at various stages of implementation. Drawing upon lessons learned from these case-studies, we provide a roadmap to facilitate applications of DNA-based methods to marine policies and to overcome the existing challenges. The main impediment identified is the need for standardized protocols to guarantee data comparison across spatial and temporal scales. Adoption of Translational Molecular Ecology - the sustained collaboration between molecular ecologists and stakeholders, will enhance consensus with regards to the objectives, methods, and outcomes of environmental management projects. Establishing a sustained dialogue among stakeholders is key to accelerating the adoption of molecular-based approaches for marine monitoring and assessment.

Marine biomonitoring with eDNA: Can metabarcoding of water samples cut it as a tool for surveying benthic communities?

In the marine realm, biomonitoring using environmental DNA (eDNA) of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at four distances (0, 0.5, 1.5, and 20 m) from the benthic habitat. Using broad range primers to amplify a ca. 313 bp fragment of the cytochrome oxidase subunit I gene, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs). The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida and Stramenopiles being the most diverse groups in benthic samples, and Hacrobia, Metazoa and Alveolata in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples and, of these, 180 (5.1%) were identified as benthic taxa that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex communities via metabarcoding.

Metabarcoding of native and invasive species in stomach contents of Great Lakes fishes

As aquatic invasive species (AIS) proliferate worldwide, a better understanding of their roles in invaded habitats is needed to inform management and introduction prevention strategies and priorities. Metabarcoding of stomach content DNA (scDNA) shows considerable promise in such regard. We thus metabarcoded scDNA from two non-native fish species (alewife (Alosa pseudoharengus) and rainbow smelt (Osmerus mordax)), and three native ones (bloater (Coregonus hoyi), ninespine stickleback (Pungitius pungitius), and slimy sculpin (Cottus cognatus)). Fishes (N = 376) were sampled in spring 2009 and 2010 from 73–128 m depths at three Lake Michigan sites. Four mitochondrial cytochrome oxidase 1 (CO1) primer sets designed to target five potential AIS prey, and a universal aquatic invertebrate CO1 primer set targeting both native and AIS prey were used. Quality controlled prey amplicons were matched to three AIS prey: Bythotrephes longimanus (mean percent frequency occurrence, all samples = 7%), Cercopagis pengoi (5%), and Dreissena rostriformis bugensis (11%). Neither invasive prey Dreissena polymorpha nor Hemimysis anomala were detected. Native prey Leptodiaptomus sicilis, Limnocalanus macrurus, and Mysis diluviana were relatively common in scDNA (respective mean percent occurrences, all samples: 48%, 25%, 42%). Analysis of variation in prey occurrences for sample site, predator species, sample year, sample depth, and predator total length (TL) indicated site and predator species were most important. However, B. longimanus occurrence in scDNA depended upon predator TL, perhaps indicative of its unique defensive spine limiting susceptibility to predation until fishes exceed species-specific gape-based limitations. Our analysis of native and invasive prey species indicated possible indirect AIS impacts such as native predators switching their diet due to AIS-driven losses of preferred native prey. Metabarcoding demonstrated that AIS are integrated components of the offshore Lake Michigan food web, with both native and non-native predators, and both invasive and native prey are affecting species interactions across multiple trophic levels.

Monitoring of biofouling communities in a Portuguese port using a combined morphological and metabarcoding approach

Marine biofouling remains an unsolved problem with a serious economic impact on several marine associated industries and constitutes a major vector for the spread of non-indigenous species (NIS). The implementation of biofouling monitoring programs allows for better fouling management and also for the early identification of NIS. However, few monitoring studies have used recent methods, such as metabarcoding, that can significantly enhance the detection of those species. Here, we employed monthly monitoring of biofouling growth on stainless steel plates in the Atlantic Port of Leixões (Northern Portugal), over one year to test the effect of commercial anti-corrosion paint in the communities. Fouling organisms were identified by combining morpho-taxonomy identification with community DNA metabarcoding using multiple markers (16S rRNA, 18S rRNA, 23S rRNA, and COI genes). The dominant colonizers found at this location were hard foulers, namely barnacles and mussels, while other groups of organisms such as cnidarians, bryozoans, and ascidians were also abundant. Regarding the temporal dynamics of the fouling communities, there was a progressive increase in the colonization of cyanobacteria, green algae, and red algae during the sampled period with the replacement of less abundant groups. The tested anticorrosion paint demonstrated to have a significant prevention effect against the biofouling community resulting in a biomass reduction. Our study also reports, for the first time, 29 NIS in this port, substantiating the need for the implementation of recurring biofouling monitoring programs in ports and harbours.

Environmental DNA can act as a biodiversity barometer of anthropogenic pressures in coastal ecosystems

Loss of biodiversity from lower to upper trophic levels reduces overall productivity and stability of coastal ecosystems in our oceans, but rarely are these changes documented across both time and space. The characterisation of environmental DNA (eDNA) from sediment and seawater using metabarcoding offers a powerful molecular lens to observe marine biota and provides a series of ‘snapshots’ across a broad spectrum of eukaryotic organisms. Using these next-generation tools and downstream analytical innovations including machine learning sequence assignment algorithms and co-occurrence network analyses, we examined how anthropogenic pressures may have impacted marine biodiversity on subtropical coral reefs in Okinawa, Japan. Based on 18 S ribosomal RNA, but not ITS2 sequence data due to inconsistent amplification for this marker, as well as proxies for anthropogenic disturbance, we show that eukaryotic richness at the family level significantly increases with medium and high levels of disturbance. This change in richness coincides with compositional changes, a decrease in connectedness among taxa, an increase in fragmentation of taxon co-occurrence networks, and a shift in indicator taxa. Taken together, these findings demonstrate the ability of eDNA to act as a barometer of disturbance and provide an exemplar of how biotic networks and coral reefs may be impacted by anthropogenic activities.

Metabarcoding reveals distinct microbiotypes in the giant clam Tridacna maxima

BACKGROUND

Giant clams and scleractinian (reef-building) corals are keystone species of coral reef ecosystems. The basis of their ecological success is a complex and fine-tuned symbiotic relationship with microbes. While the effect of environmental change on the composition of the coral microbiome has been heavily studied, we know very little about the composition and sensitivity of the microbiome associated with clams. Here, we explore the influence of increasing temperature on the microbial community (bacteria and dinoflagellates from the family Symbiodiniaceae) harbored by giant clams, maintained either in isolation or exposed to other reef species. We created artificial benthic assemblages using two coral species (Pocillopora damicornis and Acropora cytherea) and one giant clam species (Tridacna maxima) and studied the microbial community in the latter using metagenomics.

RESULTS

Our results led to three major conclusions. First, the health status of giant clams depended on the composition of the benthic species assemblages. Second, we discovered distinct microbiotypes in the studied T. maxima population, one of which was disproportionately dominated by Vibrionaceae and directly linked to clam mortality. Third, neither the increase in water temperature nor the composition of the benthic assemblage had a significant effect on the composition of the Symbiodiniaceae and bacterial communities of T. maxima.

CONCLUSIONS

Altogether, our results suggest that at least three microbiotypes naturally exist in the studied clam populations, regardless of water temperature. These microbiotypes plausibly provide similar functions to the clam host via alternate molecular pathways as well as microbiotype-specific functions. This redundancy in functions among microbiotypes together with their specificities provides hope that giant clam populations can tolerate some levels of environmental variation such as increased temperature. Importantly, the composition of the benthic assemblage could make clams susceptible to infections by Vibrionaceae, especially when water temperature increases. Video abstract.

Partitioning of diet between species and life history stages of sympatric and cryptic snappers (Lutjanidae) based on DNA metabarcoding

Lutjanus erythropterus and L. malabaricus are sympatric, sister taxa that are important to fisheries throughout the Indo-Pacific. Their juveniles are morphologically indistinguishable (i.e. cryptic). A DNA metabarcoding dietary study was undertaken to assess the diet composition and partitioning between the juvenile and adult life history stages of these two lutjanids. Major prey taxa were comprised of teleosts and crustaceans for all groups except adult L. erythropterus, which instead consumed soft bodied invertebrates (e.g. tunicates, comb jellies and medusae) as well as teleosts, with crustaceans being notably absent. Diet composition was significantly different among life history stages and species, which may be associated with niche habitat partitioning or differences in mouth morphology within adult life stages. This study provides the first evidence of diet partitioning between cryptic juveniles of overlapping lutjanid species, thus providing new insights into the ecological interactions, habitat associations, and the specialised adaptations required for the coexistence of closely related species. This study has improved our understanding of the differential contributions of the juvenile and adult diets of these sympatric species within food webs. The diet partitioning reported in this study was only revealed by the taxonomic resolution provided by the DNA metabarcoding approach and highlights the potential utility of this method to refine the dietary components of reef fishes more generally.

eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate-related threats and as such there is a pressing need for cost-effective whole-ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high-resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)-a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km² area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

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.

Considerations for metabarcoding-based port biological baseline surveys aimed at marine nonindigenous species monitoring and risk assessments

Monitoring introduction and spread of nonindigenous species via maritime transport and performing risk assessments require port biological baseline surveys. Yet, the comprehensiveness of these surveys is often compromised by the large number of habitats present in a port, the seasonal variability, and the time-consuming morphological approach used for taxonomic identification. Metabarcoding represents a promising alternative for rapid comprehensive port biological baseline surveys, but its application in this context requires further assessments.We applied metabarcoding (based on barcodes of the cytochrome c oxidase subunit I and of the 18S ribosomal RNA gene) to 192 port samples collected (a) from diverse habitats (water column-including environmental DNA and zooplankton, sediment, and fouling structures), (b) at different sites (from inner to outer estuary), and iii) during the four seasons of the year.By comparing the biodiversity metrics derived from each sample group, we show that each sampling method resulted in a distinct community profile and that environmental DNA alone cannot substitute for organismal sampling, and that, although sampling at different seasons and locations resulted in higher observed biodiversity, operational results can be obtained by sampling selected locations and seasons.By assessing the taxonomic composition of the samples, we show that metabarcoding data allowed the detection of previously recorded nonindigenous species as well as to reveal presence of new ones, even if in low abundance. Synthesis and application. Our comprehensive assessment of metabarcoding for port biological baseline surveys sets the basics for cost-effective, standardized, and comprehensive monitoring of nonindigenous species and for performing risk assessments in ports. This development will contribute to the implementation of the recently entered into force International Convention for the Control and Management of Ships' Ballast Water and Sediments.

Clinical metagenomics for infectious corneal ulcers: Rags to riches?

The emergence of clinical metagenomics as an unbiased, hypothesis-free approach to diagnostic testing is set to fundamentally alter the way infectious diseases are detected. Long envisioned as the solution to the limitations of culture-based conventional microbiology, next generation sequencing methods will soon mature, and our attention will inevitably turn to how they can be applied to areas of medicine which need it most urgently. In ophthalmology, the demand for this technology is particularly pressing for the care of infectious corneal ulcers, where current diagnostic tests may fail to identify a causative organism in over half of cases. However, the optimism found in the budding discourse surrounding clinical metagenomics belies the reality that clinicians and scientists will soon be inundated by oppressive volumes of sequencing data, much of which will be foreign and unfamiliar. Therefore, our success in translating clinical metagenomics is likely to hinge on how we make sense of these data, and understanding its implications for the interpretation and implementation of sequencing into routine clinical care. In this consortium-led review, we provide an outline of these data-related issues and how they may be used to inform technical workflows, with the hope that we may edge closer to realizing the potential of clinical metagenomics for this important unmet need.

Diversity and Spatial Distribution of Chromophytic Phytoplankton in the Bay of Bengal Revealed by RuBisCO Genes (rbcL)

Phytoplankton are the basis of primary production and play important roles in regulating energy export in marine ecosystems. Compared to other regions, chromophytic phytoplankton are considerably understudied in the Bay of Bengal (BOB). Here, we investigated community structure and spatial distribution of chromophytic phytoplankton in the BOB by using RuBisCO genes (Form ID rbcL). High throughput sequencing of rbcL genes revealed that diatoms, cyanobacteria (Cyanophyceae), Pelagophyceae, Haptophyceae, Chrysophyceae, Eustigamatophyceae, Xanthophyceae, Cryptophyceae, Dictyochophyceae, and Pinguiophyceae were the most abundant groups recovered in the BOB. Abundances and distribution of diatoms and Pelagophyceae were further verified using quantitative PCR analyses which showed the dominance of these groups near the Equator region (p < 0.01) where upwelling was likely the source of nutrients. Further, redundancy analysis demonstrated that temperature was an important environmental driver in structuring distributions of Cyanophyceae and dominant chromophytic phytoplankton. Morphological identification and quantification confirmed the dominance of diatoms, and also detected other cyanobacteria and dinoflagellates that were missing in our molecular characterizations. Pearson’s correlations of these morphologically identified phytoplankton with environmental gradients also indicated that nutrients and temperature were key variables shaping community structure. Combination of molecular characterization and morphological identification provided a comprehensive overview of chromophytic phytoplankton. This is the first molecular study of chromophytic phytoplankton accomplished in the BOB, and our results highlight a combination of molecular analysis targeting rbcL genes and microscopic detection in examining phytoplankton composition and diversity.

Universal human papillomavirus typing by whole genome sequencing following target enrichment: evaluation of assay reproducibility and limit of detection

BACKGROUND

We recently described a method for unbiased detection of all known human papillomaviruses (HPV) types with the potential for the determination of their variant and integration from the resulting whole genome sequence data. Considering the complex workflow for target-enriched next generation sequencing (NGS), we focused on the reproducibility and limit of detection (LOD) of this new universal HPV typing assay in this study.

RESULTS

We evaluated the reproducibility and LOD for HPV genotyping based on our recently published method that used RNA-baits targeting whole genomes of 191 HPV types, Agilent SureSelect protocol for target enrichment and Illumina HiSeq 2500 for sequencing (eWGS, enriched whole genome sequencing). Two libraries, prepared from pooled plasmids representing 9 vaccine HPV types at varying input (1-625 copies/reaction), were sequenced twice giving four replicates for evaluating reproducibility and LOD. eWGS showed high correlation in the number of reads mapped to HPV reference genomes between the two flow-cell lanes within (R² = 1) and between experiments (R² = 0.99). The number of mapped reads was positively correlated to copy number (β = 13.9, p < 0.0001). The limit of blank (LOB) could be calculated based on mapped reads to HPV types not included in each sample. HPV genotyping was reproducible for all 9 types at 625 copies using multiple cut-off criteria but LOD was 25 copies based on number of reads above LOB even when multiple types were present. eWGS showed no bias for HPV genotyping under single or multiple infection (p = 0.16-0.99).

CONCLUSIONS

The universal eWGS method for HPV genotyping has sensitivity, competitive with widely used consensus PCR methods with reduced type competition, and with the potential for determination of variant and integration status. The protocol used in this study, using defined samples varying in complexity and copy number, analyzed in replicate and duplicate assays, is applicable to most WGS methods.

Deciphering host-parasitoid interactions and parasitism rates of crop pests using DNA metabarcoding

An accurate estimation of parasitism rates and diversity of parasitoids of crop insect pests is a prerequisite for exploring processes leading to efficient natural biocontrol. Traditional methods such as rearing have been often limited by taxonomic identification, insect mortality and intensive work, but the advent of high-throughput sequencing (HTS) techniques, such as DNA metabarcoding, is increasingly seen as a reliable and powerful alternative approach. Little has been done to explore the benefits of such an approach for estimating parasitism rates and parasitoid diversity in an agricultural context. In this study, we compared the composition of parasitoid species and parasitism rates between rearing and DNA metabarcoding of host eggs and larvae of the millet head miner, Heliocheilus albipunctella De Joannis (Lepidoptera, Noctuidae), collected from millet fields in Senegal. We first assessed the detection threshold for the main ten endoparasitoids, by sequencing PCR products obtained from artificial dilution gradients of the parasitoid DNAs in the host moth. We then assessed the potential of DNA metabarcoding for diagnosing parasitism rates in samples collected from the field. Under controlled conditions, our results showed that relatively small quantities of parasitoid DNA (0.07 ng) were successfully detected within an eight-fold larger quantity of host DNA. Parasitoid diversity and parasitism rate estimates were always higher for DNA metabarcoding than for host rearing. Furthermore, metabarcoding detected multi-parasitism, cryptic parasitoid species and differences in parasitism rates between two different sampling sites. Metabarcoding shows promise for gaining a clearer understanding of the importance and complexity of host-parasitoid interactions in agro-ecosystems, with a view to improving pest biocontrol strategies.

DNA metabarcoding assays reveal a diverse prey assemblage for Mobula rays in the Bohol Sea, Philippines

Diet studies provide base understanding of trophic structure and are a valuable initial step for many fields of marine ecology, including conservation and fisheries biology. Considerable complexity in marine trophic structure can exist due to the presence of highly mobile species with long life spans. Mobula rays are highly mobile, large, planktivorous elasmobranchs that are frequently caught either directly or as bycatch in fisheries, which, combined with their conservative life history strategy, makes their populations susceptible to decline in intensely fished regions. Effective management of these iconic and vulnerable species requires an understanding of the diets that sustain them, which can be difficult to determine using conventional sampling methods. We use three DNA metabarcode assays to identify 44 distinct taxa from the stomachs (n = 101) of four sympatric Mobula ray species (Mobula birostris, Mobula tarapacana, Mobula japanica, and Mobula thurstoni) caught over 3 years (2013-2015) in a direct fishery off Bohol in the Philippines. The diversity and incidence of bony fishes observed in ray diets were unprecedented. Nevertheless, rays showed dietary overlap, with krill (Euphausia) dominating their diet. Our results provide a more detailed assessment of sympatric ray diets than was previously described and reveal the complexity that can exist in food webs at critical foraging habitats.

Differential effects of coral-giant clam assemblages on biofouling formation

To prevent the settlement and/or the growth of fouling organisms (i.e. bacteria, fungi or microalgae), benthic sessile species have developed various defense mechanisms among which the production of chemical molecules. While studies have mostly focused on the release of chemical compounds by single species, there exist limited data on multi-species assemblages. We used an integrative approach to explore the potential interactive effects of distinct assemblages of two corals species and one giant clam species on biofouling appearance and composition. Remarkably, we found distinct biofouling communities suggesting the importance of benthic sessile assemblages in biofouling control. Moreover, the assemblage of 3 species led to an inhibition of biofouling, likely through a complex of secondary metabolites. Our results highlight that through their different effect on their near environment, species assemblages might be of upmost importance for their survival and therefore, should now be taken into account for sustainable management of coral reefs.

Digging for DNA at depth: rapid universal metabarcoding surveys (RUMS) as a tool to detect coral reef biodiversity across a depth gradient

Background

Effective biodiversity monitoring is fundamental in tracking changes in ecosystems as it relates to commercial, recreational, and conservation interests. Current approaches to survey coral reef ecosystems center on the use of indicator species and repeat surveying at specific sites. However, such approaches are often limited by the narrow snapshot of total marine biodiversity that they describe and are thus hindered in their ability to contribute to holistic ecosystem-based monitoring. In tandem, environmental DNA (eDNA) and next-generation sequencing metabarcoding methods provide a new opportunity to rapidly assess the presence of a broad spectrum of eukaryotic organisms within our oceans, ranging from microbes to macrofauna.

Methods

We here investigate the potential for rapid universal metabarcoding surveys (RUMS) of eDNA in sediment samples to provide snapshots of eukaryotic subtropical biodiversity along a depth gradient at two coral reefs in Okinawa, Japan based on 18S rRNA.

Results

Using 18S rRNA metabarcoding, we found that there were significant separations in eukaryotic community assemblages (at the family level) detected in sediments when compared across different depths ranging from 10 to 40 m (p = 0.001). Significant depth zonation was observed across operational taxonomic units assigned to the class Demospongiae (sponges), the most diverse class (contributing 81% of species) within the phylum Porifera; the oldest metazoan phylum on the planet. However, zonation was not observed across the class Anthozoa (i.e., anemones, stony corals, soft corals, and octocorals), suggesting that the former may serve as a better source of indicator species based on sampling over fine spatial scales and using this universal assay. Furthermore, despite their abundance on the examined coral reefs, we did not detect any octocoral DNA, which may be due to low cellular shedding rates, assay sensitivities, or primer biases.

Discussion

Overall, our pilot study demonstrates the importance of exploring depth effects in eDNA and suggest that RUMS may be applied to provide a baseline of information on eukaryotic marine taxa at coastal sites of economic and conservation importance.

Species-level biodiversity assessment using marine environmental DNA metabarcoding requires protocol optimization and standardization

DNA extraction from environmental samples (environmental DNA; eDNA) for metabarcoding-based biodiversity studies is gaining popularity as a noninvasive, time-efficient, and cost-effective monitoring tool. The potential benefits are promising for marine conservation, as the marine biome is frequently under-surveyed due to its inaccessibility and the consequent high costs involved. With increasing numbers of eDNA-related publications have come a wide array of capture and extraction methods. Without visual species confirmation, inconsistent use of laboratory protocols hinders comparability between studies because the efficiency of target DNA isolation may vary. We determined an optimal protocol (capture and extraction) for marine eDNA research based on total DNA yield measurements by comparing commonly employed methods of seawater filtering and DNA isolation. We compared metabarcoding results of both targeted (small taxonomic group with species-level assignment) and universal (broad taxonomic group with genus/family-level assignment) approaches obtained from replicates treated with the optimal and a low-performance capture and extraction protocol to determine the impact of protocol choice and DNA yield on biodiversity detection. Filtration through cellulose-nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit outperformed other combinations of capture and extraction methods, showing a ninefold improvement in DNA yield over the poorest performing methods. Use of optimized protocols resulted in a significant increase in OTU and species richness for targeted metabarcoding assays. However, changing protocols made little difference to the OTU and taxon richness obtained using universal metabarcoding assays. Our results demonstrate an increased risk of false-negative species detection for targeted eDNA approaches when protocols with poor DNA isolation efficacy are employed. Appropriate optimization is therefore essential for eDNA monitoring to remain a powerful, efficient, and relatively cheap method for biodiversity assessments. For seawater, we advocate filtration through cellulose-nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit or phenol-chloroform-isoamyl for successful implementation of eDNA multi-marker metabarcoding surveys.

High-throughput identification and diagnostics of pathogens and pests: Overview and practical recommendations

High-throughput identification technologies provide efficient tools for understanding the ecology and functioning of microorganisms. Yet, these methods have been only rarely used for monitoring and testing ecological hypotheses in plant pathogens and pests in spite of their immense importance in agriculture, forestry and plant community dynamics. The main objectives of this manuscript are the following: (a) to provide a comprehensive overview about the state-of-the-art high-throughput quantification and molecular identification methods used to address population dynamics, community ecology and host associations of microorganisms, with a specific focus on antagonists such as pathogens, viruses and pests; (b) to compile available information and provide recommendations about specific protocols and workable primers for bacteria, fungi, oomycetes and insect pests; and (c) to provide examples of novel methods used in other microbiological disciplines that are of great potential use for testing specific biological hypotheses related to pathology. Finally, we evaluate the overall perspectives of the state-of-the-art and still evolving methods for diagnostics and population- and community-level ecological research of pathogens and pests.

How quantitative is metabarcoding: A meta-analytical approach

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta‐analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.

Evaluating metabarcoding to analyse diet composition of species foraging in anthropogenic landscapes using Ion Torrent and Illumina sequencing

DNA metabarcoding of faecal samples is being successfully used to study the foraging niche of species. We assessed the ability of two benchtop high-throughput sequencing (HTS) platforms, to identify a large taxonomic array of food items from domestic cats Felis silvestris catus, including prey and human-related food taxa (pet food and leftovers leaving undetectable solid remains in faeces). Scats from a captive feeding trial (n = 41) and from free-ranging individuals (n = 326) were collected and analysed using a cytb mini-barcode in independent PCR replicates on the Ion PGM and the MiSeq platforms. Outputs from MiSeq were more sensitive and reproducible than those from Ion PGM due to a higher sequencing depth and sequence quality on MiSeq. DNA from intact prey taxa was detected more often (82% of the expected occurrences) than DNA from pet food (54%) and raw fish and meat (31%). We assumed that this variability was linked to different degree of DNA degradation: The Ion PGM detected significantly less human-linked food, birds, field voles, murids and shrews in the field-collected samples than the MiSeq platform. Pooling the replicates from both platforms and filtering the data allowed identification of at least one food item in 87.4% of the field-collected samples. Our DNA metabarcoding approach identified 29 prey taxa, of which 25 to species level (90% of items) including 9 rodents, 3 insectivores, 12 birds and 1 reptile and 33 human-related food taxa of which 23 were identified to genus level (75% of items). Our results demonstrate that using HTS platforms such as MiSeq, which provide reads of sufficiently high quantity and quality, with sufficient numbers of technical replicates, is a robust and non-invasive approach for further dietary studies on animals foraging on a wide range of food items in anthropogenic landscapes.