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

Challenges for assessing vertebrate diversity in turbid Saharan water-bodies using environmental DNA

The Sahara desert is the largest warm desert in the world and a poorly explored area. Small water-bodies occur across the desert and are crucial habitats for vertebrate biodiversity. Environmental DNA (eDNA) is a powerful tool for species detection and is being increasingly used to conduct biodiversity assessments. However, there are a number of difficulties with sampling eDNA from such turbid water-bodies and it is often not feasible to rely on electrical tools in remote desert environments. We trialled a manually powered filtering method in Mauritania, using pre-filtration to circumvent problems posed by turbid water in remote arid areas. From nine vertebrate species expected in the water-bodies, four were detected visually, two via metabarcoding, and one via both methods. Difficulties filtering turbid water led to severe constraints, limiting the sampling protocol to only one sampling point per study site, which alone may largely explain why many of the expected vertebrate species were not detected. The amplification of human DNA using general vertebrate primers is also likely to have contributed to the low number of taxa identified. Here we highlight a number of challenges that need to be overcome to successfully conduct metabarcoding eDNA studies for vertebrates in desert environments in Africa.

Historical baselines in marine bioinvasions: Implications for policy and management

The human-mediated introduction of marine non-indigenous species is a centuries- if not millennia-old phenomenon, but was only recently acknowledged as a potent driver of change in the sea. We provide a synopsis of key historical milestones for marine bioinvasions, including timelines of (a) discovery and understanding of the invasion process, focusing on transfer mechanisms and outcomes, (b) methodologies used for detection and monitoring, (c) approaches to ecological impacts research, and (d) management and policy responses. Early (until the mid-1900s) marine bioinvasions were given little attention, and in a number of cases actively and routinely facilitated. Beginning in the second half of the 20th century, several conspicuous non-indigenous species outbreaks with strong environmental, economic, and public health impacts raised widespread concerns and initiated shifts in public and scientific perceptions. These high-profile invasions led to policy documents and strategies to reduce the introduction and spread of non-indigenous species, although with significant time lags and limited success and focused on only a subset of transfer mechanisms. Integrated, multi-vector management within an ecosystem-based marine management context is urgently needed to address the complex interactions of natural and human pressures that drive invasions in marine ecosystems.

A Rapid and Precise Mutation-Activated Fluorescence Reporter for Analyzing Acute Mutagenesis Frequency

Mutagenesis reporters are critical for quantifying genome stability. However, current methods rely on cell survival/death to report mutation, which takes weeks and prevents evaluation of acute or time-dependent changes. Existing methods also have other limitations, such as cell type restrictions. Using our discovery that mCherryFP fluorescence depends on residue Trp98, we replaced this codon with a stop codon to generate a mutation biosensor (termed CherryOFF), with a green fluorescence protein (GFP) as an internal control. We found that the red fluorescence of this biosensor is activated by a specific A/T-G/C nucleotide transition. Compared with the established hypoxanthine phosphoribosyl transferase assay, our reporter has similar or better ability to detect changes of mutation frequency induced by physical/chemical mutagens or manipulation of mutation-related genes. Furthermore, CherryOFF-GFP can report mutagenesis independently of cell-death events, can be adapted to many cell types, and can generate readouts within 1 day for the measurement of acute or time-dependent events.

Metabarcoding and metabolomics offer complementarity in deciphering marine eukaryotic biofouling community shifts

Metabarcoding and metabolomics were used to explore the taxonomic composition and functional diversity of eukaryotic biofouling communities on plates with antifouling paints at two French coastal sites: Lorient (North Eastern Atlantic Ocean; temperate and eutrophic) and Toulon (North-Western Mediterranean Sea; mesotrophic but highly contaminated). Four distinct coatings were tested at each site and season for one month. Metabarcoding showed biocidal coatings had less impact on eukaryotic assemblages compared to spatial and temporal effects. Ciliophora, Chlorophyceae or Cnidaria (mainly hydrozoans) were abundant at Lorient, whereas Arthropoda (especially crustaceans), Nematoda, and Ochrophyta dominated less diversified assemblages at Toulon. Seasonal shifts were observed at Lorient, but not Toulon. Metabolomics also showed clear site discrimination, but these were associated with a coating and not season dependent clustering. The meta-omics analysis enabled identifications of some associative patterns between metabolomic profiles and specific taxa, in particular those colonizing the plates with biocidal coatings at Lorient.

Nucleic acids-based tools for ballast water surveillance, monitoring, and research

Understanding the risks of biological invasion posed by ballast water-whether in the context of compliance testing, routine monitoring, or basic research-is fundamentally an exercise in biodiversity assessment, and as such should take advantage of the best tools available for tackling that problem. The past several decades have seen growing application of genetic methods for the study of biodiversity, driven in large part by dramatic technological advances in nucleic acids analysis. Monitoring approaches based on such methods have the potential to increase dramatically sampling throughput for biodiversity assessments, and to improve on the sensitivity, specificity, and taxonomic accuracy of traditional approaches. The application of targeted detection tools (largely focused on PCR but increasingly incorporating novel probe-based methodologies) has led to a paradigm shift in rare species monitoring, and such tools have already been applied for early detection in the context of ballast water surveillance. Rapid improvements in community profiling approaches based on high throughput sequencing (HTS) could similarly impact broader efforts to catalogue biodiversity present in ballast tanks, and could provide novel opportunities to better understand the risks of biotic exchange posed by ballast water transport-and the effectiveness of attempts to mitigate those risks. These various approaches still face considerable challenges to effective implementation, depending on particular management or research needs. Compliance testing, for instance, remains dependent on accurate quantification of viable target organisms; while tools based on RNA detection show promise in this context, the demands of such testing require considerable additional investment in methods development. In general surveillance and research contexts, both targeted and community-based approaches are still limited by various factors: quantification remains a challenge (especially for taxa in larger size classes), gaps in nucleic acids reference databases are still considerable, uncertainties in taxonomic assignment methods persist, and many applications have not yet matured sufficiently to offer standardized methods capable of meeting rigorous quality assurance standards. Nevertheless, the potential value of these tools, their growing utilization in biodiversity monitoring, and the rapid methodological advances over the past decade all suggest that they should be seriously considered for inclusion in the ballast water surveillance toolkit.

A cross-taxa study using environmental DNA/RNA metabarcoding to measure biological impacts of offshore oil and gas drilling and production operations

Standardized ecosystem-based monitoring surveys are critical for providing information on marine ecosystem health. Environmental DNA/RNA (eDNA/eRNA) metabarcoding may facilitate such surveys by quickly and effectively characterizing multi-trophic levels. In this study, we assessed the suitability of eDNA/eRNA metabarcoding to evaluate changes in benthic assemblages of bacteria, Foraminifera and other eukaryotes along transects at three offshore oil and gas (O&G) drilling and production sites, and compared these to morphologically characterized macro-faunal assemblages. Bacterial communities were the most responsive to O&G activities, followed by Foraminifera, and macro-fauna (the latter assessed by morphology). The molecular approach enabled detection of hydrocarbon degrading taxa such as the bacteria Alcanivorax and Microbulbifer at petroleum impacted stations. Most identified indicator taxa, notably among macro-fauna, were highly specific to site conditions. Based on our results we suggest that eDNA/eRNA metabarcoding can be used as a stand-alone method for biodiversity assessment or as a complement to morphology-based monitoring approaches.

Censusing marine eukaryotic diversity in the twenty-first century

The ocean constitutes one of the vastest and richest biomes on our planet. Most recent estimations, all based on indirect approaches, suggest that there are millions of marine eukaryotic species. Moreover, a large majority of these are small (less than 1 mm), cryptic and still unknown to science. However, this knowledge gap, caused by the lack of diagnostic morphological features in small organisms and the limited sampling of the global ocean, is currently being filled, thanks to new DNA-based approaches. The molecular technique of PCR amplification of homologous gene regions combined with high-throughput sequencing, routinely used to census unculturable prokaryotes, is now also being used to characterize whole communities of marine eukaryotes. Here, we review how this methodological advancement has helped to better quantify the magnitude and patterns of marine eukaryotic diversity, with an emphasis on taxonomic groups previously largely overlooked. We then discuss obstacles remaining to achieve a global understanding of marine eukaryotic diversity. In particular, we argue that 18S variable regions do not provide sufficient taxonomic resolution to census marine life, and suggest combining broad eukaryotic surveys targeting the 18S rRNA region with more taxon-focused analyses of hypervariable regions to improve our understanding of the diversity of species, the functional units of marine ecosystems.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

DNA metabarcoding for high-throughput monitoring of estuarine macrobenthic communities

Morphology-based profiling of benthic communities has been extensively applied to aquatic ecosystems' health assessment. However, it remains a low-throughput, and sometimes ambiguous, procedure. Despite DNA metabarcoding has been applied to marine benthos, a comprehensive approach providing species-level identifications for estuarine macrobenthos is still lacking. Here we report a combination of experimental and field studies to assess the aptitude of COI metabarcoding to provide robust species-level identifications for high-throughput monitoring of estuarine macrobenthos. To investigate the ability of metabarcoding to detect all species present in bulk DNA extracts, we contrived three phylogenetically diverse communities, and applied four different primer pairs to generate PCR products within the COI barcode region. Between 78-83% of the species in the contrived communities were recovered through HTS. Subsequently, we compared morphology and metabarcoding-based approaches to determine the species composition from four distinct estuarine sites. Our results indicate that species richness would be considerably underestimated if only morphological methods were used: globally 27 species identified through morphology versus 61 detected by metabarcoding. Although further refinement is required to improve efficiency and output of this approach, here we show the great aptitude of COI metabarcoding to provide high quality and auditable species identifications in estuarine macrobenthos monitoring.

Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications

High-throughput sequencing metabarcoding studies in marine biosecurity have largely focused on targeting environmental DNA (eDNA). DNA can persist extracellularly in the environment, making discrimination of living organisms difficult. In this study, bilge water samples (i.e., water accumulating on-board a vessel during transit) were collected from 15 small recreational and commercial vessels. eDNA and eRNA molecules were co-extracted and the V4 region of the 18S ribosomal RNA gene targeted for metabarcoding. In total, 62.7% of the Operational Taxonomic Units (OTUs) were identified at least once in the corresponding eDNA and eRNA reads, with 19.5% unique to eDNA and 17.7% to eRNA. There were substantial differences in diversity between molecular compartments; 57% of sequences from eDNA-only OTUs belonged to fungi, likely originating from legacy DNA. In contrast, there was a higher percentage of metazoan (50.2%) and ciliate (31.7%) sequences in the eRNA-only OTUs. Our data suggest that the presence of eRNA-only OTUs could be due to increased cellular activities of some rare taxa that were not identified in the eDNA datasets, unusually high numbers of rRNA transcripts in ciliates, and/or artefacts produced during the reverse transcriptase, PCR and sequencing steps. The proportions of eDNA/eRNA shared and unshared OTUs were highly heterogeneous within individual bilge water samples. Multiple factors including boat type and the activities performed on-board, such as washing of scientific equipment, may play a major role in contributing to this variability. For some marine biosecurity applications analysis, eDNA-only data may be sufficient, however there are an increasing number of instances where distinguishing the living portion of a community is essential. For these circumstances, we suggest only including OTUs that are present in both eDNA and eRNA data. OTUs found only in the eRNA data need to be interpreted with caution until further research provides conclusive evidence for their origin.

Early detection monitoring for aquatic non-indigenous species: Optimizing surveillance, incorporating advanced technologies, and identifying research needs

Following decades of ecologic and economic impacts from a growing list of nonindigenous and invasive species, government and management entities are committing to systematic early- detection monitoring (EDM). This has reinvigorated investment in the science underpinning such monitoring, as well as the need to convey that science in practical terms to those tasked with EDM implementation. Using the context of nonindigenous species in the North American Great Lakes, this article summarizes the current scientific tools and knowledge - including limitations, research needs, and likely future developments - relevant to various aspects of planning and conducting comprehensive EDM. We begin with the scope of the effort, contrasting target-species with broad-spectrum monitoring, reviewing information to support prioritization based on species and locations, and exploring the challenge of moving beyond individual surveys towards a coordinated monitoring network. Next, we discuss survey design, including effort to expend and its allocation over space and time. A section on sample collection and analysis overviews the merits of collecting actual organisms versus shed DNA, reviews the capabilities and limitations of identification by morphology, DNA target markers, or DNA barcoding, and examines best practices for sample handling and data verification. We end with a section addressing the analysis of monitoring data, including methods to evaluate survey performance and characterize and communicate uncertainty. Although the body of science supporting EDM implementation is already substantial, research and information needs (many already actively being addressed) include: better data to support risk assessments that guide choice of taxa and locations to monitor; improved understanding of spatiotemporal scales for sample collection; further development of DNA target markers, reference barcodes, genomic workflows, and synergies between DNA-based and morphology-based taxonomy; and tools and information management systems for better evaluating and communicating survey outcomes and uncertainty.

Next-Generation Sequencing in Diagnostic Pathology

Interrogation of tissue informs on patient management through delivery of a diagnosis together with associated clinically relevant data. The diagnostic pathologist will usually evaluate the morphological appearances of a tissue sample and, occasionally, the pattern of expression of a limited number of biomarkers. Recent developments in sequencing technology mean that DNA and RNA from tissue samples can now be interrogated in great detail. These new technologies, collectively known as next-generation sequencing (NGS), generate huge amounts of data which can be used to support patient management. In order to maximize the utility of tissue interrogation, the molecular data need to be interpreted and integrated with the morphological data. However, in order to interpret the molecular data, the pathologist must understand the utility and the limitations of NGS data. In this review, the principles behind NGS technologies are described. In addition, the caveats in the interpretation of the data are discussed, and a scheme is presented to "classify" the types of data which are generated. Finally, a glossary of new terminology is included to help pathologists become familiar with the lexicon of NGS-derived molecular data.

Ecosystem biomonitoring with eDNA: metabarcoding across the tree of life in a tropical marine environment

Effective marine management requires comprehensive data on the status of marine biodiversity. However, efficient methods that can document biodiversity in our oceans are currently lacking. Environmental DNA (eDNA) sourced from seawater offers a new avenue for investigating the biota in marine ecosystems. Here, we investigated the potential of eDNA to inform on the breadth of biodiversity present in a tropical marine environment. Directly sequencing eDNA from seawater using a shotgun approach resulted in only 0.34% of 22.3 million reads assigning to eukaryotes, highlighting the inefficiency of this method for assessing eukaryotic diversity. In contrast, using 'tree of life' (ToL) metabarcoding and 20-fold fewer sequencing reads, we could detect 287 families across the major divisions of eukaryotes. Our data also show that the best performing 'universal' PCR assay recovered only 44% of the eukaryotes identified across all assays, highlighting the need for multiple metabarcoding assays to catalogue biodiversity. Lastly, focusing on the fish genus Lethrinus, we recovered intra- and inter-specific haplotypes from seawater samples, illustrating that eDNA can be used to explore diversity beyond taxon identifications. Given the sensitivity and low cost of eDNA metabarcoding we advocate this approach be rapidly integrated into biomonitoring programs.

DNA in a bottle-Rapid metabarcoding survey for early alerts of invasive species in ports

Biota monitoring in ports is increasingly needed for biosecurity reasons and safeguarding marine biodiversity from biological invasion. Present and future international biosecurity directives can be accomplished only if the biota acquired by maritime traffic in ports is controlled. Methodologies for biota inventory are diverse and now rely principally on extensive and labor-intensive sampling along with taxonomic identification by experts. In this study, we employed an extremely simplified environmental DNA (eDNA) sampling methodology from only three 1-L bottles of water per port, followed by metabarcoding (high-throughput sequencing and DNA-based species identification) using 18S rDNA and Cytochrome oxidase I as genetic barcodes. Eight Bay of Biscay ports with available inventory of fouling invertebrates were employed as a case study. Despite minimal sampling efforts, three invasive invertebrates were detected: the barnacle Austrominius modestus, the tubeworm Ficopomatus enigmaticus and the polychaete Polydora triglanda. The same species have been previously found from visual and DNA barcoding (genetic identification of individuals) surveys in the same ports. The current costs of visual surveys, conventional DNA barcoding and this simplified metabarcoding protocol were compared. The results encourage the use of metabarcoding for early biosecurity alerts.

Adaptation of Surface-Associated Bacteria to the Open Ocean: A Genomically Distinct Subpopulation of Phaeobacter gallaeciensis Colonizes Pacific Mesozooplankton

The marine Roseobacter group encompasses numerous species which occupy a large variety of ecological niches. However, members of the genus Phaeobacter are specifically adapted to a surface-associated lifestyle and have so far been found nearly exclusively in disjunct, man-made environments including shellfish and fish aquacultures, as well as harbors. Therefore, the possible natural habitats, dispersal and evolution of Phaeobacter spp. have largely remained obscure. Applying a high-throughput cultivation strategy along a longitudinal Pacific transect, the present study revealed for the first time a widespread natural occurrence of Phaeobacter in the marine pelagial. These bacteria were found to be specifically associated to mesoplankton where they constitute a small but detectable proportion of the bacterial community. The 16S rRNA gene sequences of 18 isolated strains were identical to that of Phaeobacter gallaeciensis DSM26640^T but sequences of internal transcribed spacer and selected genomes revealed that the strains form a distinct clade within P. gallaeciensis. The genomes of the Pacific and the aquaculture strains were highly conserved and had a fraction of the core genome of 89.6%, 80 synteny breakpoints, and differed 2.2% in their nucleotide sequences. Diversification likely occurred through neutral mutations. However, the Pacific strains exclusively contained two active Type I restriction modification systems which is commensurate with a reduced acquisition of mobile elements in the Pacific clade. The Pacific clade of P. gallaeciensis also acquired a second, homolog phosphonate transport system compared to all other P. gallaeciensis. Our data indicate that a previously unknown, distinct clade of P. gallaeciensis acquired a limited number of clade-specific genes that were relevant for its association with mesozooplankton and for colonization of the marine pelagial. The divergence of the Pacific clade most likely was driven by the adaptation to this novel ecological niche rather than by geographic isolation.

Development of a real-time polymerase chain reaction assay for the detection of the invasive Mediterranean fanworm, Sabella spallanzanii, in environmental samples

The Mediterranean fanworm, Sabella spallanzanii Gmelin 1791, was first detected in the Southern Hemisphere in the 1990s and is now abundant in many parts of southern Australia and in several locations around northern New Zealand. Once established, it can proliferate rapidly, reaching high densities with potential ecological and economic impacts. Early detection of new S. spallanzanii incursions is important to prevent its spread, guide eradication or control efforts and to increase knowledge on the species' dispersal pathways. In this study, we developed a TaqMan probe real-time polymerase chain reaction assay targeting a region of the mitochondrial cytochrome oxidase I gene. The assay was validated in silico and in vitro using DNA from New Zealand and Australian Sabellidae with no cross-reactivity detected. The assay has a linear range of detection over seven orders of magnitude with a limit of detection reached at 12.4 × 10^-4 ng/μL of DNA. We analysed 145 environmental (water, sediment and biofouling) samples and obtained positive detections only from spiked samples and those collected at a port where S. spallanzanii is known to be established. This assay has the potential to enhance current morphological and molecular-based methods, through its ability to rapidly and accurately identify S. spallanzanii in environmental samples.

Metabarcoding monitoring analysis: the pros and cons of using co-extracted environmental DNA and RNA data to assess offshore oil production impacts on benthic communities

Sequencing environmental DNA (eDNA) is increasingly being used as an alternative to traditional morphological-based identification to characterize biological assemblages and monitor anthropogenic impacts in marine environments. Most studies only assess eDNA which, compared to eRNA, can persist longer in the environment after cell death. Therefore, eRNA may provide a more immediate census of the environment due to its relatively weaker stability, leading some researchers to advocate for the use of eRNA as an additional, or perhaps superior proxy for portraying ecological changes. A variety of pre-treatment techniques for screening eDNA and eRNA derived operational taxonomic units (OTUs) have been employed prior to statistical analyses, including removing singleton taxa (i.e., OTUs found only once) and discarding those not present in both eDNA and eRNA datasets. In this study, we used bacterial (16S ribosomal RNA gene) and eukaryotic (18S ribosomal RNA gene) eDNA- and eRNA-derived data from benthic communities collected at increasing distances along a transect from an oil production platform (Taranaki, New Zealand). Macro-infauna (visual classification of benthic invertebrates) and physico-chemical data were analyzed in parallel. We tested the effect of removing singleton taxa, and removing taxa not present in the eDNA and eRNA libraries from the same environmental sample (trimmed by shared OTUs), by comparing the impact of the oil production platform on alpha- and beta-diversity of the eDNA/eRNA-based biological assemblages, and by correlating these to the morphologically identified macro-faunal communities and the physico-chemical data. When trimmed by singletons, presence/absence information from eRNA data represented the best proxy to detect changes on species diversity for both bacteria and eukaryotes. However, assessment of quantitative beta-diversity from read abundance information of bacteria eRNA did not, contrary to eDNA, reveal any impact from the oil production activity. Overall, the data appeared more robust when trimmed by shared OTUs, showing a greater effect of the platform on alpha- and beta-diversity. Trimming by shared OTUs likely removes taxa derived from legacy DNA and technical artefacts introduced through reverse transcriptase, polymerase-chain-reaction and sequencing. Findings from our scoping study suggest that metabarcoding-based biomonitoring surveys should, if funds, time and expertise allow, be assessed using both eDNA and eRNA products.

Community barcoding reveals little effect of ocean acidification on the composition of coastal plankton communities: Evidence from a long-term mesocosm study in the Gullmar Fjord, Skagerrak

The acidification of the oceans could potentially alter marine plankton communities with consequences for ecosystem functioning. While several studies have investigated effects of ocean acidification on communities using traditional methods, few have used genetic analyses. Here, we use community barcoding to assess the impact of ocean acidification on the composition of a coastal plankton community in a large scale, in situ, long-term mesocosm experiment. High-throughput sequencing resulted in the identification of a wide range of planktonic taxa (Alveolata, Cryptophyta, Haptophyceae, Fungi, Metazoa, Hydrozoa, Rhizaria, Straminipila, Chlorophyta). Analyses based on predicted operational taxonomical units as well as taxonomical compositions revealed no differences between communities in high CO2 mesocosms (~ 760 μatm) and those exposed to present-day CO2 conditions. Observed shifts in the planktonic community composition were mainly related to seasonal changes in temperature and nutrients. Furthermore, based on our investigations, the elevated CO2 did not affect the intraspecific diversity of the most common mesozooplankter, the calanoid copepod Pseudocalanus acuspes. Nevertheless, accompanying studies found temporary effects attributed to a raise in CO2. Differences in taxa composition between the CO2 treatments could, however, only be observed in a specific period of the experiment. Based on our genetic investigations, no compositional long-term shifts of the plankton communities exposed to elevated CO2 conditions were observed. Thus, we conclude that the compositions of planktonic communities, especially those in coastal areas, remain rather unaffected by increased CO2.

Sensitivity and accuracy of high-throughput metabarcoding methods for early detection of invasive fish species

High-throughput DNA metabarcoding has gained recognition as a potentially powerful tool for biomonitoring, including early detection of aquatic invasive species (AIS). DNA based techniques are advancing, but our understanding of the limits to detection for metabarcoding complex samples is inadequate. For detecting AIS at an early stage of invasion when the species is rare, accuracy at low detection limits is key. To evaluate the utility of metabarcoding in future fish community monitoring programs, we conducted several experiments to determine the sensitivity and accuracy of routine metabarcoding methods. Experimental mixes used larval fish tissue from multiple “common” species spiked with varying proportions of tissue from an additional “rare” species. Pyrosequencing of genetic marker, COI (cytochrome c oxidase subunit I) and subsequent sequence data analysis provided experimental evidence of low-level detection of the target “rare” species at biomass percentages as low as 0.02% of total sample biomass. Limits to detection varied interspecifically and were susceptible to amplification bias. Moreover, results showed some data processing methods can skew sequence-based biodiversity measurements from corresponding relative biomass abundances and increase false absences. We suggest caution in interpreting presence/absence and relative abundance in larval fish assemblages until metabarcoding methods are optimized for accuracy and precision.

Marine invasions enter the genomic era: three lessons from the past, and the way forward

The expanding scale and increasing rate of marine biological invasions have been documented since the early 20th century. Besides their global ecological and economic impacts, non-indigenous species (NIS) also have attracted much attention as opportunities to explore important eco-evolutionary processes such as rapid adaptation, long-distance dispersal and range expansion, and secondary contacts between divergent evolutionary lineages. In this context, genetic tools have been extensively used in the past 20 years. Three important issues appear to have emerged from such studies. First, the study of NIS has revealed unexpected cryptic diversity in what had previously been assumed homogeneous entities. Second, there has been surprisingly little evidence of strong founder events accompanying marine introductions, a pattern possibly driven by large propagule loads. Third, the evolutionary processes leading to successful invasion have been difficult to ascertain due to faint genetic signals. Here we explore the potential of novel tools associated with high-throughput sequencing (HTS) to address these still pressing issues. Dramatic increase in the number of loci accessible via HTS has the potential to radically increase the power of analyses aimed at species delineation, exploring the population genomic consequences of range expansions, and examining evolutionary processes such as admixture, introgression, and adaptation. Nevertheless, the value of this new wealth of genomic data will ultimately depend on the ability to couple it with expanded "traditional" efforts, including exhaustive sampling of marine populations over large geographic scales, integrated taxonomic analyses, and population level exploration of quantitative trait differentiation through common-garden and other laboratory experiments.

Priming of a DNA metabarcoding approach for species identification and inventory in marine macrobenthic communities

In marine and estuarine benthic communities, the inventory and estimation of species richness are often hampered by the need for broad taxonomic expertise across several phyla. The use of DNA metabarcoding has emerged as a powerful tool for the fast assessment of species composition in a diversity of ecological communities. Here, we tested the amplification success of five primer sets targeting different COI-5P regions by 454 pyrosequencing to maximize the recovery of two simulated macrobenthic communities containing 21 species (SimCom1 and SimCom 2). Species identification was first performed against a compiled reference library of macrobenthic species. Reads with similarity results to reference sequences between 70% and 97% were then submitted to GenBank and BOLD to attempt the identification of concealed species in the bulk sample. The combination of at least three primer sets was able to recover more species than any primer set alone, achieving 85% of represented species in SimCom1 and 76% in SimCom2. Our approach was successful to detect low-frequency specimens, as well as concealed species, in the bulk sample, indicating the potential for the application of this approach on marine bioassessment and inventory, including the detection of "hidden" biodiversity that would hardly be possible based on morphology only.

DNA barcodes for bio-surveillance: regulated and economically important arthropod plant pests

Many of the arthropod species that are important pests of agriculture and forestry are impossible to discriminate morphologically throughout all of their life stages. Some cannot be differentiated at any life stage. Over the past decade, DNA barcoding has gained increasing adoption as a tool to both identify known species and to reveal cryptic taxa. Although there has not been a focused effort to develop a barcode library for them, reference sequences are now available for 77% of the 409 species of arthropods documented on major pest databases. Aside from developing the reference library needed to guide specimen identifications, past barcode studies have revealed that a significant fraction of arthropod pests are a complex of allied taxa. Because of their importance as pests and disease vectors impacting global agriculture and forestry, DNA barcode results on these arthropods have significant implications for quarantine detection, regulation, and management. The current review discusses these implications in light of the presence of cryptic species in plant pests exposed by DNA barcoding.

Real-time PCR detection of Didemnum perlucidum (Monniot, 1983) and Didemnum vexillum (Kott, 2002) in an applied routine marine biosecurity context

Prevention and early detection are well recognized as the best strategies for minimizing the risks posed by nonindigenous species (NIS) that have the potential to become marine pests. Central to this is the ability to rapidly and accurately identify the presence of NIS, often from complex environmental samples like biofouling and ballast water. Molecular tools have been increasingly applied to assist with the identification of NIS and can prove particularly useful for taxonomically difficult groups like ascidians. In this study, we have developed real-time PCR assays suited to the specific identification of the ascidians Didemnum perlucidum and Didemnum vexillum. Despite being recognized as important global pests, this is the first time specific molecular detection methods have been developed that can support the early identification and detection of these species from a broad range of environmental sample types. These fast, robust and high-throughput assays represent powerful tools for routine marine biosecurity surveillance, as detection and confirmation of the early presence of species could assist in the timely establishment of emergency responses and control strategies. This study applied the developed assays to confirm the ability to detect Didemnid eDNA in water samples. While previous work has focused on detection of marine larvae from water samples, the development of real-time PCR assays specifically aimed at detecting eDNA of sessile invertebrate species in the marine environment represents a world first and a significant step forwards in applied marine biosecurity surveillance. Demonstrated success in the detection of D. perlucidum eDNA from water samples at sites where it could not be visually identified suggests value in incorporating such assays into biosecurity survey designs targeting Didemnid species.

Metabarcoding improves detection of eukaryotes from early biofouling communities: implications for pest monitoring and pathway management

In this experimental study the patterns in early marine biofouling communities and possible implications for surveillance and environmental management were explored using metabarcoding, viz. 18S ribosomal RNA gene barcoding in combination with high-throughput sequencing. The community structure of eukaryotic assemblages and the patterns of initial succession were assessed from settlement plates deployed in a busy port for one, five and 15 days. The metabarcoding results were verified with traditional morphological identification of taxa from selected experimental plates. Metabarcoding analysis identified > 400 taxa at a comparatively low taxonomic level and morphological analysis resulted in the detection of 25 taxa at varying levels of resolution. Despite the differences in resolution, data from both methods were consistent at high taxonomic levels and similar patterns in community shifts were observed. A high percentage of sequences belonging to genera known to contain non-indigenous species (NIS) were detected after exposure for only one day.

18S rRNA V9 metabarcoding for diet characterization: a critical evaluation with two sympatric zooplanktivorous fish species

The potential of the 18S rRNA V9 metabarcoding approach for diet assessment was explored using MiSeq paired‐end (PE; 2 × 150 bp) technology. To critically evaluate the method′s performance with degraded/digested DNA, the diets of two zooplanktivorous fish species from the Bay of Biscay, European sardine (Sardina pilchardus) and European sprat (Sprattus sprattus), were analysed. The taxonomic resolution and quantitative potential of the 18S V9 metabarcoding was first assessed both in silico and with mock and field plankton samples. Our method was capable of discriminating species within the reference database in a reliable way providing there was at least one variable position in the 18S V9 region. Furthermore, it successfully discriminated diet between both fish species, including habitat and diel differences among sardines, overcoming some of the limitations of traditional visual‐based diet analysis methods. The high sensitivity and semi‐quantitative nature of the 18S V9 metabarcoding approach was supported by both visual microscopy and qPCR‐based results. This molecular approach provides an alternative cost and time effective tool for food‐web analysis.

eDNA and specific primers for early detection of invasive species--A case study on the bivalve Rangia cuneata, currently spreading in Europe

Intense human activities facilitate the successful spread and establishment of non-indigenous aquatic organisms in marine and freshwater ecosystems. In some cases such intrusions result in noticeable and adverse changes in the recipient environments. In the Baltic Sea, the discovery and rapid initial spread of the North American wedge clam Rangia cuneata represents a new wave of invasion which may trigger unpredictable changes of the local benthic communities. In this study we present a species-specific DNA-based marker developed in silico and experimentally tested on environmental samples. Marker specificity and sensitivity were assessed in vitro from water samples containing different mixtures of the target species and other five bivalves currently present in the region: the native Cerastoderma glaucum, Macoma balthica and Mytilus trossulus, the invasive Dreissena polymorpha and the cryptogenic Mya arenaria. Cross-species amplification was not found in any case. The method allows to detecting at least 0.4 ng of R. cuneata DNA per μl, and 0.1 g of tissue per liter of water. Finally, the marker performance was assessed in water samples from the Baltic Sea and Vistula Lagoon. The coincidence between independent visual observations of R. cuneata and positive PCR amplification of the marker from the water samples confirmed the efficiency of this highly reproducible, fast, and technically easy method. R. cuneata traces can be detected from environmental DNA even when the population is sparse and small, enabling rapid management responses and allowing to track the invasion dynamics.

Detecting nuisance species using NGST: Methodology shortcomings and possible application in ballast water monitoring

Detecting the presence of potential invasive species in ballast water is a priority for preventing their spread into new environments. Next generation sequencing technologies are being increasingly used for exploring and assessing biodiversity from environmental samples. Here we apply high throughput sequencing from DNA extracted from ballast water (BW) samples employing two different platforms, Ion Torrent and 454, and compare the putative species catalogues from the resulting Operational Taxonomic Units (OTU). Water samples were taken from the RV Polastern ballast tank in five different days between the second and the twentieth navigation day. Pronounced decrease of oxygen concentration and increase of temperature occurred in the BW during this time, coincident with a progressively higher proportion of unassigned OTU and short reads indicating DNA degradation. Discrepancy between platforms for species catalogues was consistent with previously published bias in AT-rich sequences for Ion Torrent platform. Some putative species detected from the two platforms increased in frequency during the Polarstern travel, which suggests they were alive and therefore tolerant to adverse conditions. OTU assigned to the highly invasive red alga Polysiphonia have been detected at low but increasing frequency from the two platforms. Although in this moment NGST could not replace current methods of sampling, sorting and individual taxonomic identification of BW biota, it has potential as an exploratory methodology especially for detecting scarce species.

Accurate assessment of the impact of salmon farming on benthic sediment enrichment using foraminiferal metabarcoding

Assessing the environmental impact of salmon farms on benthic systems is traditionally undertaken using biotic indices derived from microscopic analyses of macrobenthic infaunal (MI) communities. In this study, we tested the applicability of using foraminiferal-specific high-throughput sequencing (HTS) metabarcoding for monitoring these habitats. Sediment samples and physico-chemical data were collected along an enrichment gradient radiating out from three Chinook salmon (Oncorhynchus tshawytscha) farms in New Zealand. HTS of environmental DNA and RNA (eDNA/eRNA) resulted in 1,875,300 sequences that clustered into 349 Operational Taxonomic Units. Strong correlations were observed among various biotic indices calculated from MI data and normalized fourth-root transformed HTS data. Correlations were stronger using eRNA compared to eDNA data. Quantile regression spline analyses identified 12 key foraminiferal taxa that have potential to be used as bioindicator species. This study demonstrates the huge potential for using this method for biomonitoring of fish-farming and other marine industrial activities.

Metabarcoding approach for nonindigenous species surveillance in marine coastal waters

In this study, high-throughput sequencing (HTS) metabarcoding was applied for the surveillance of plankton communities within the southeastern (SE) Baltic Sea coastal zone. These results were compared with those from routine monitoring survey and morphological analyses. Four of five nonindigenous species found in the samples were identified exclusively by metabarcoding. All of them are considered as invasive in the Baltic Sea with reported impact on the ecosystem and biodiversity. This study indicates that, despite some current limitations, HTS metabarcoding can provide information on the presence of exotic species and advantageously complement conventional approaches, only requiring the same monitoring effort as before. Even in the currently immature status of HTS, this combination of HTS metabarcoding and observational records is recommended in the early detection of marine pests and delivery of the environmental status metrics of nonindigenous species.

Deep-Sea, Deep-Sequencing: Metabarcoding Extracellular DNA from Sediments of Marine Canyons

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

NGS metabarcoding proves successful for quantitative assessment of symbiont abundance: the case of feather mites on birds

Understanding the ecological function of species and the structure of communities is crucial in the study of ecological interactions among species. For this purpose, not only the occurrence of particular species but also their abundance in ecological communities is required. However, abundance quantification of species through morphological characters is often difficult or time/money consuming when dealing with elusive or small taxa. Here we tested the use of next-generation sequencing (NGS) for abundance estimation of two species of feather mites (Proctophyllodes stylifer and Pteronyssoides parinus) under five proportions (16:1, 16:4, 16:16, 16:64, and 16:256 mites) against a mock community composed by Proctophyllodes clavatus and Proctophyllodes sylviae. In all mixtures, we retrieved sequence reads from all species. We found a strong linear relationship between 454 reads and the real proportion of individuals in the mixture for both focal species. The slope for Pr. stylifer was close to one (0.904), and the intercept close to zero (-0.007), thus showing an almost perfect correspondence between real and estimated proportions. The slope for Pt. parinus was 0.351 and the intercept 0.307, showing that while the estimated proportion increased linearly relative to real proportions of individuals in the samples, proportions were overestimated at low real proportions and underestimated at larger ones. Additionally, pyrosequencing replicates from each DNA extraction were highly repeatable (R = 0.920 and 0.972, respectively), showing that the quantification method is highly consistent given a DNA extract. Our study suggests that NGS is a promising tool for abundance estimation of feather mites' communities in birds.

Metabarcoding approach for the ballast water surveillance--an advantageous solution or an awkward challenge?

Transfer of organisms with ships' ballast water is recognized as a major pathway of non-indigenous species introduction and addressed in a few recent legislative initiatives. Among other they imply scientific and technical research and monitoring to be conducted in a efficient and reliable way. The recent development of DNA barcoding and metabarcoding technologies opens new opportunities for biodiversity and biosecurity surveillance. In the current study, the performance of metabarcoding approach was assessed in comparison to the conventional (visual) observations, during the en route experimental ballast water survey. Opportunities and limitations of the molecular method were identified from taxonomical datasets rendered by two molecular markers of different degree of universality - the universal cytochrome oxydase sub-unit I gene and a fragment of RuBisCO gene. The cost-efficacy and possible improvements of these methods are discussed for the further successful development and implementation of the approach in ballast water control and NIS surveillance.

Early detection of eukaryotic communities from marine biofilm using high-throughput sequencing: an assessment of different sampling devices

Marine biofilms are precursors for colonization by larger fouling organisms, including non-indigenous species (NIS). In this study, high-throughput sequencing (HTS) of 18S rRNA metabarcodes was used to investigate four sampling methods (modified syringe, sterilized sponge, underwater tape and sterilized swab) for characterizing eukaryotic communities in marine biofilms. Perspex™ plates were sampled in and out of water. DNA collected with tape did not amplify. Otherwise, there were no statistical differences in communities among the remaining three sampling devices or between the two environments. Sterilized sponges are recommended for ease of use underwater. In-depth HTS analysis identified diverse eukaryotic communities, dominated by Metazoa and Chromoalveolata. Among the latter, diatoms (Bacillariophyceae) were particularly abundant (33% of reads assigned to Chromalveolata). The NIS Ciona savignyi was detected in all samples. The application of HTS in marine biofilm surveillance could facilitate early detection of NIS, improving the probability of successful eradication.

First identification of tetrodotoxin (TTX) in the flatworm Stylochoplana sp.; a source of TTX for the sea slug Pleurobranchaea maculata

High concentrations of the neurotoxin tetrodotoxin (TTX) were detected by liquid chromatography-mass spectrometry (LC-MS) in the Platyhelminthes Stylochoplana sp. from Pilot Bay (Tauranga, New Zealand). This is the first detection of TTX in this genus. Concentrations were monitored from March to November (2013) and found to significantly decrease from a peak in July (avg. 551 mg kg(-1)) to November (avg. 140 mg kg(-1)). Stylochoplana sp. co-occurred with TTX-containing Pleurobranchaea maculata (Opisthobranchia). A Stylochoplana sp.-specific real-time PCR assay was developed targeting the mitochondrial cytochrome c oxidase subunit I gene to determine if P. maculata consumed Stylochoplana sp. Positive Stylochoplana sp. signals were obtained for 7 of 19 P. maculata tested. Mass calculations indicate Stylochoplana sp. could supply Pilot Bay P. maculata with the TTX required to account for the concentrations reported in previous studies (ca. 1.04 mg TTX per individual) based on an ingestion rate of one individual every 2-3 days throughout their lifetime. However, due to the lack of Stylochoplana sp. in areas with dense P. maculata populations, and high concentration (ca. 1400 mg kg(-1)) of TTX detected in some individuals, it is unlikely that Stylochoplana sp. represent the sole source of TTX in P. maculata.

Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods

The quantification of the biological diversity in environmental samples using high-throughput DNA sequencing is hindered by the PCR bias caused by variable primer-template mismatches of the individual species. In some dietary studies, there is the added problem that samples are enriched with predator DNA, so often a predator-specific blocking oligonucleotide is used to alleviate the problem. However, specific blocking oligonucleotides could coblock nontarget species to some degree. Here, we accurately estimate the extent of the PCR biases induced by universal and blocking primers on a mock community prepared with DNA of twelve species of terrestrial arthropods. We also compare universal and blocking primer biases with those induced by variable annealing temperature and number of PCR cycles. The results show that reads of all species were recovered after PCR enrichment at our control conditions (no blocking oligonucleotide, 45 °C annealing temperature and 40 cycles) and high-throughput sequencing. They also show that the four factors considered biased the final proportions of the species to some degree. Among these factors, the number of primer-template mismatches of each species had a disproportionate effect (up to five orders of magnitude) on the amplification efficiency. In particular, the number of primer-template mismatches explained most of the variation (~3/4) in the amplification efficiency of the species. The effect of blocking oligonucleotide concentration on nontarget species relative abundance was also significant, but less important (below one order of magnitude). Considering the results reported here, the quantitative potential of the technique is limited, and only qualitative results (the species list) are reliable, at least when targeting the barcoding COI region.

Performance comparison of genetic markers for high-throughput sequencing-based biodiversity assessment in complex communities

Metabarcode surveys of DNA extracted from environmental samples are increasingly popular for biodiversity assessment in natural communities. Such surveys rely heavily on robust genetic markers. Therefore, analysis of PCR efficiency and subsequent biodiversity estimation for different types of genetic markers and their corresponding primers is important. Here, we test the PCR efficiency and biodiversity recovery potential of three commonly used genetic markers - nuclear small subunit ribosomal DNA (18S), mitochondrial cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (mt16S) - using 454 pyrosequencing of a zooplankton community collected from Hamilton Harbour, Ontario. We found that biodiversity detection power and PCR efficiency varied widely among these markers. All tested primers for COI failed to provide high-quality PCR products for pyrosequencing, but newly designed primers for 18S and 16S passed all tests. Furthermore, multiple analyses based on large-scale pyrosequencing (i.e. 1/2 PicoTiter plate for each marker) showed that primers for 18S recover more (38 orders) groups than 16S (10 orders) across all taxa, and four vs. two orders and nine vs. six families for Crustacea. Our results showed that 18S, using newly designed primers, is an efficient and powerful tool for profiling biodiversity in largely unexplored communities, especially when amplification difficulties exist for mitochondrial markers such as COI. Universal primers for higher resolution markers such as COI are still needed to address the possible low resolution of 18S for species-level identification.

Comparison of three next-generation sequencing platforms for metagenomic sequencing and identification of pathogens in blood

Background

The introduction of benchtop sequencers has made adoption of whole genome sequencing possible for a broader community of researchers than ever before. Concurrently, metagenomic sequencing (MGS) is rapidly emerging as a tool for interrogating complex samples that defy conventional analyses. In addition, next-generation sequencers are increasingly being used in clinical or related settings, for instance to track outbreaks. However, information regarding the analytical sensitivity or limit of detection (LoD) of benchtop sequencers is currently lacking. Furthermore, the specificity of sequence information at or near the LoD is unknown.

Results

In the present study, we assess the ability of three next-generation sequencing platforms to identify a pathogen (viral or bacterial) present in low titers in a clinically relevant sample (blood). Our results indicate that the Roche-454 Titanium platform is capable of detecting Dengue virus at titers as low as 1X10^2.5 pfu/mL, corresponding to an estimated 5.4X10⁴ genome copies/ml maximum. The increased throughput of the benchtop sequencers, the Ion Torrent PGM and Illumina MiSeq platforms, enabled detection of viral genomes at concentrations as low as 1X10⁴ genome copies/mL. Platform-specific biases were evident in sequence read distributions as well as viral genome coverage. For bacterial samples, only the MiSeq platform was able to provide sequencing reads that could be unambiguously classified as originating from Bacillus anthracis.

Conclusion

The analytical sensitivity of all three platforms approaches that of standard qPCR assays. Although all platforms were able to detect pathogens at the levels tested, there were several noteworthy differences. The Roche-454 Titanium platform produced consistently longer reads, even when compared with the latest chemistry updates for the PGM platform. The MiSeq platform produced consistently greater depth and breadth of coverage, while the Ion Torrent was unequaled for speed of sequencing. None of the platforms were able to verify a single nucleotide polymorphism responsible for antiviral resistance in an Influenza A strain isolated from the 2009 H1N1 pandemic. Overall, the benchtop platforms perform well for identification of pathogens from a representative clinical sample. However, unlike identification, characterization of pathogens is likely to require higher titers, multiple libraries and/or multiple sequencing runs.