Assessing the performance and efficiency of environmental DNA / RNA capture methodologies under controlled experimental conditions

Evaluating eDNA and eRNA metabarcoding for aquatic biodiversity assessment: From bacteria to vertebrates

The monitoring and management of aquatic ecosystems depend on precise estimates of biodiversity. Metabarcoding analyses of environmental nucleic acids (eNAs), including environmental DNA (eDNA) and environmental RNA (eRNA), have garnered attention for their cost-effective and non-invasive biomonitoring capabilities. However, the accuracy of biodiversity estimates obtained through eNAs can vary among different organismal groups. Here we evaluate the performance of eDNA and eRNA metabarcoding across nine organismal groups, ranging from bacteria to terrestrial vertebrates, in three cross-sections of the Yangtze River, China. We observe robust complementarity between eDNA and eRNA data. The relative detectability of eNAs was notably influenced by major taxonomic groups and organismal sizes, with eDNA providing more robust signals for larger organisms. Both eDNA and eRNA exhibited similar cross-sectional and longitudinal patterns. However, the detectability of larger organisms declined in eRNA metabarcoding, possibly due to differential RNA release and decay among different organismal groups or sizes. While underscoring the potential of eDNA and eRNA in large river biomonitoring, we emphasize the need for differential interpretation of eDNA versus eRNA data. This highlights the importance of careful method selection and interpretation in biomonitoring studies.

Seasonal dynamics and environmental drivers of tissue and mucus microbiomes in the staghorn coral Acropora pulchra

Background

Rainfall-induced coastal runoff represents an important environmental impact in near-shore coral reefs that may affect coral-associated bacterial microbiomes. Shifts in microbiome community composition and function can stress corals and ultimately cause mortality and reef declines. Impacts of environmental stress may be site specific and differ between coral microbiome compartments (e.g., tissue versus mucus). Coastal runoff and associated water pollution represent a major stressor for near-shore reef-ecosystems in Guam, Micronesia.

Methods

Acropora pulchra colonies growing on the West Hagåtña reef flat in Guam were sampled over a period of 8 months spanning the 2021 wet and dry seasons. To examine bacterial microbiome diversity and composition, samples of A. pulchra tissue and mucus were collected during late April, early July, late September, and at the end of December. Samples were collected from populations in two different habitat zones, near the reef crest (farshore) and close to shore (nearshore). Seawater samples were collected during the same time period to evaluate microbiome dynamics of the waters surrounding coral colonies. Tissue, mucus, and seawater microbiomes were characterized using 16S DNA metabarcoding in conjunction with Illumina sequencing. In addition, water samples were collected to determine fecal indicator bacteria (FIB) concentrations as an indicator of water pollution. Water temperatures were recorded using data loggers and precipitation data obtained from a nearby rain gauge. The correlation structure of environmental parameters (temperature and rainfall), FIB concentrations, and A. pulchra microbiome diversity was evaluated using a structural equation model. Beta diversity analyses were used to investigate spatio-temporal trends of microbiome composition.

Results

Acropora pulchra microbiome diversity differed between tissues and mucus, with mucus microbiome diversity being similar to the surrounding seawater. Rainfall and associated fluctuations of FIB concentrations were correlated with changes in tissue and mucus microbiomes, indicating their role as drivers of A. pulchra microbiome diversity. A. pulchra tissue microbiome composition remained relatively stable throughout dry and wet seasons; tissues were dominated by Endozoicomonadaceae, coral endosymbionts and putative indicators of coral health. In nearshore A. pulchra tissue microbiomes, Simkaniaceae, putative obligate coral endosymbionts, were more abundant than in A. pulchra colonies growing near the reef crest (farshore). A. pulchra mucus microbiomes were more diverse during the wet season than the dry season, a distinction that was also associated with drastic shifts in microbiome composition. This study highlights the seasonal dynamics of coral microbiomes and demonstrates that microbiome diversity and composition may differ between coral tissues and the surface mucus layer.

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

Drop it all: extraction-free detection of targeted marine species through optimized direct droplet digital PCR

Molecular biomonitoring programs increasingly use environmental DNA (eDNA) for detecting targeted species such as marine non-indigenous species (NIS) or endangered species. However, the current molecular detection workflow is cumbersome and time-demanding, and thereby can hinder management efforts and restrict the "opportunity window" for rapid management responses. Here, we describe a direct droplet digital PCR (direct-ddPCR) approach to detect species-specific free-floating extra-cellular eDNA (free-eDNA) signals, i.e., detection of species-specific eDNA without the need for filtration or DNA extraction, with seawater samples. This first proof-of-concept aquarium study was conducted with three distinct marine species: the Mediterranean fanworm Sabella spallanzanii, the ascidian clubbed tunicate Styela clava, and the brown bryozoan Bugula neritina to evaluate the detectability of free-eDNA in seawater. The detectability of targeted free-eDNA was assessed by directly analysing aquarium marine water samples using an optimized species-specific ddPCR assay. The results demonstrated the consistent detection of S. spallanzanii and B. neritina free-eDNA when these organisms were present in high abundance. Once organisms were removed, the free-eDNA signal exponentially declined, noting that free-eDNA persisted between 24-72 h. Results indicate that organism biomass, specimen characteristics (e.g., stress and viability), and species-specific biological differences may influence free-eDNA detectability. This study represents the first step in assessing the feasibility of direct-ddPCR technology for the detection of marine species. Our results provide information that could aid in the development of new technology, such as a field development of ddPCR systems, which could allow for automated continuous monitoring of targeted marine species, enabling point-of-need detection and rapid management responses.

Comparison of concentration and extraction workflows for qPCR quantification of intI1 and vanA in untreated wastewater

Quantitative polymerase chain reaction (qPCR) measurement of antibiotic resistance genes (ARGs) in untreated municipal wastewater may prove useful in combating the antimicrobial resistance crisis. However, harmonizing and optimizing qPCR-based workflows is essential to facilitate comparisons across studies, and includes achieving highly-effective ARG capture through efficient concentration and extraction procedures. In the current study, combinations of sample volume, membrane types and DNA extraction kits within filtration and centrifugation-based workflows were used to quantify 16S ribosomal RNA (16S rRNA), class 1 integron-integrase gene (intI1) and an ARG encoding resistance to vancomycin (vanA) in untreated wastewater sampled from three wastewater treatment plants (WWTPs). Highly abundant 16S rRNA and intI1 were detected in 100 % of samples from all three WWTPs using both 2 and 20 mL sample volumes, while lower prevalence vanA was only detected when using the 20 mL volume. When filtering 2 mL of wastewater, workflows with 0.20-/0.40-μm polycarbonate (PC) membranes generally yielded greater concentrations of the three targets than workflows with 0.22-/0.45-μm mixed cellulose ester (MCE) membranes. The improved performance was diminished when the sample volume was increased to 20 mL. Consistently greater concentrations of 16S rRNA, intI1 and vanA were yielded by filtration-based workflows using PC membranes combined with a DNeasy PowerWater (DPW) Kit, regardless of the sample volume used, and centrifugation-based workflows with DNeasy Blood & Tissue Kit for 2-mL wastewater extractions. Within the filtration-based workflows, the DPW kit yielded more detection and quantifiable results for less abundant vanA than the DNeasy PowerSoil Pro Kit and FastDNA™ SPIN Kit for Soil. These findings indicate that the performance of qPCR-based workflows for surveillance of ARGs in wastewater varies across targets, sample volumes, concentration methods and extraction kits. Workflows must be carefully considered and validated considering the target ARGs to be monitored.

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.

Methodological considerations for aqueous environmental RNA collection, preservation, and extraction

Environmental RNA (eRNA) analysis is expected to infer species' physiological information (health status, developmental stage, and environmental stress response) and their distribution and composition more correctly than environmental DNA (eDNA) analysis. With the prospect of such eRNA applications, there is an increasing need for technological development for efficient eRNA detection because of its physicochemical instability. The present study conducted a series of aquarium experiments using zebrafish (Danio rerio) and validated the methodologies for capture, preservation, and extraction of eRNA in a water sample. In the eRNA extraction experiment, an approximately 1.5-fold increase in lysis buffer volume resulted in a more than sixfold increase in target eRNA concentration. In the eRNA capture experiment, although GF/F and GF/A filters yielded similar eRNA concentrations, a GF/A filter may be capable of passing through more volume of water samples and consequently collecting more eRNA particles, given the time required for water filtration. In the eRNA preservation experiment, the use of RNA stabilization reagent (RNAlater) allowed for stably preserving target eRNA on a filter sample at - 20 and even 4 °C for 6 days at least. Altogether, the findings enable the improvement of eRNA availability from the field and easily preserve eRNA samples without deep-freezing, which will contribute to the refinement of eRNA analysis for biological and physiological monitoring in aquatic ecosystems.

Net overboard: Comparing marine eDNA sampling methodologies at sea to unravel marine biodiversity

Environmental DNA (eDNA) analyses are powerful for describing marine biodiversity but must be optimized for their effective use in routine monitoring. To maximize eDNA detection probabilities of sparsely distributed populations, water samples are usually concentrated from larger volumes and filtered using fine-pore membranes, often a significant cost-time bottleneck in the workflow. This study aimed to streamline eDNA sampling by investigating plankton net versus bucket sampling, direct versus sequential filtration including self-preserving filters. Biodiversity was assessed using metabarcoding of the small ribosomal subunit (18S rRNA) and mitochondrial cytochrome c oxidase I (COI) genes. Multispecies detection probabilities were estimated for each workflow using a probabilistic occupancy modelling approach. Significant workflow-related differences in biodiversity metrics were reported. Highest amplicon sequence variant (ASV) richness was attained by the bucket sampling combined with self-preserving filters, comprising a large portion of microplankton. Less diversity but more metazoan taxa were captured in the net samples combined with 5 μm pore size filters. Prefiltered 1.2 μm samples yielded few or no unique ASVs. The highest average (~32%) metazoan detection probabilities in the 5 μm pore size net samples confirmed the effectiveness of preconcentration plankton for biodiversity screening. These results contribute to streamlining eDNA sampling protocols for uptake and implementation in marine biodiversity research and surveillance.

Capturing marine microbiomes and environmental DNA: A field sampling guide

The expanding interest in marine microbiome and eDNA sequence data has led to a demand for sample collection and preservation standard practices to enable comparative assessments of results across studies and facilitate meta-analyses. We support this effort by providing guidelines based on a review of published methods and field sampling experiences. The major components considered here are environmental and resource considerations, sample processing strategies, sample storage options, and eDNA extraction protocols. It is impossible to provide universal recommendations considering the wide range of eDNA applications; rather, we provide information to design fit-for-purpose protocols. To manage scope, the focus here is on sampling collection and preservation of prokaryotic and microeukaryotic eDNA. Even with a focused view, the practical utility of any approach depends on multiple factors, including habitat type, available resources, and experimental goals. We broadly recommend enacting rigorous decontamination protocols, pilot studies to guide the filtration volume needed to characterize the target(s) of interest and minimize PCR inhibitor collection, and prioritizing sample freezing over (only) the addition of preservation buffer. An annotated list of studies that test these parameters is included for more detailed investigation on specific steps. To illustrate an approach that demonstrates fit-for-purpose methodologies, we provide a protocol for eDNA sampling aboard an oceanographic vessel. These guidelines can aid the decision-making process for scientists interested in sampling and sequencing marine microbiomes and/or eDNA.

Optimal sample type and number vary in small shallow lakes when targeting non-native fish environmental DNA

Non-native fish have been shown to have deleterious impacts on freshwater ecosystems in New Zealand. Early detection is critical for their effective management. Traditional capture-based techniques may not detect newly introduced fish, especially if they are present in low abundance. Molecular techniques that target environmental DNA (eDNA) have been shown, in many instances, to be more sensitive, cost-effective and require lower sampling effort. However, appropriate sampling strategies are needed to ensure robust and interpretable data are obtained. In this study we used droplet digital PCR assays to investigate the presence of two non-native fish in New Zealand, the European perch (Perca fluviatilis) and rudd (Scardinius erythrophthalmus) in three small lakes. Samples were collected from water and surface sediment at near-shore and mid-lake sites. Probabilistic modelling was used to assess the occupancy of fish eDNA and develop guidance on sampling strategies. Based on the detection probability measures from the present study, at least six sites and five replicates per site are needed to reliably detect fish eDNA in sediment samples, and twelve sites with eight replicates per site for water samples. The results highlight the potential of developing monitoring and surveillance programs adapted to lakes, that include the use of assays targeting eDNA. This study focused on small shallow lakes, and it is likely that these recommendations may vary in larger, deeper, and more geomorphologically complex lakes, and this requires further research.

Temperature Controls eDNA Persistence across Physicochemical Conditions in Seawater

Environmental DNA (eDNA) quantification and sequencing are emerging techniques for assessing biodiversity in marine ecosystems. Environmental DNA can be transported by ocean currents and may remain at detectable concentrations far from its source depending on how long it persist. Thus, predicting the persistence time of eDNA is crucial to defining the spatial context of the information derived from it. To investigate the physicochemical controls of eDNA persistence, we performed degradation experiments at temperature, pH, and oxygen conditions relevant to the open ocean and the deep sea. The eDNA degradation process was best explained by a model with two phases with different decay rate constants. During the initial phase, eDNA degraded rapidly, and the rate was independent of physicochemical factors. During the second phase, eDNA degraded slowly, and the rate was strongly controlled by temperature, weakly controlled by pH, and not controlled by dissolved oxygen concentration. We demonstrate that marine eDNA can persist at quantifiable concentrations for over 2 weeks at low temperatures (≤10 °C) but for a week or less at ≥20 °C. The relationship between temperature and eDNA persistence is independent of the source species. We propose a general temperature-dependent model to predict the maximum persistence time of eDNA detectable through single-species eDNA quantification methods.