Critical considerations for communicating environmental DNA science.

The economic and methodological efficiencies of environmental DNA (eDNA) based survey approaches provide an unprecedented opportunity to assess and monitor aquatic environments. However, instances of inadequate communication from the scientific community about confidence levels, knowledge gaps, reliability, and appropriate parameters of eDNA-based methods have hindered their uptake in environmental monitoring programs and, in some cases, has created misperceptions or doubts in the management community. To help remedy this situation, scientists convened a session at the Second National Marine eDNA Workshop to discuss strategies for improving communications with managers. These include articulating the readiness of different eDNA applications, highlighting the strengths and limitations of eDNA tools for various applications or use cases, communicating uncertainties associated with specified uses transparently, and avoiding the exaggeration of exploratory and preliminary findings. Several key messages regarding implementation, limitations, and relationship to existing methods were prioritized. To be inclusive of the diverse managers, practitioners, and researchers, we and the other workshop participants propose the development of communication workflow plans, using RACI (Responsible, Accountable, Consulted, Informed) charts to clarify the roles of all pertinent individuals and parties and to minimize the chance for miscommunications. We also propose developing decision support tools such as Structured Decision-Making (SDM) to help balance the benefits of eDNA sampling with the inherent uncertainty, and developing an eDNA readiness scale to articulate the technological readiness of eDNA approaches for specific applications. These strategies will increase clarity and consistency regarding our understanding of the utility of eDNA-based methods, improve transparency, foster a common vision for confidently applying eDNA approaches, and enhance their benefit to the monitoring and assessment community.

Tourmaline: A containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake.

BACKGROUND: Amplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results. FINDINGS: We developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files-a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences-it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm; assigns taxonomy to the resulting representative sequences; performs analyses of taxonomic, alpha, and beta diversity; and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S ribosomal RNA amplicon data set from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks. CONCLUSION: Automated workflows like Tourmaline enable rapid analysis of environmental amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline.

Dominance of Sulfurospirillum in Metagenomes Associated with the Methane Ice Worm (Sirsoe methanicola).

Sirsoe methanicola, commonly known as the methane ice worm, is the only macrofaunal species known to inhabit the Gulf of Mexico methane hydrates. Little is known about this elusive marine polychaete that can colonize rich carbon and energy reserves. Metagenomic analysis of gut contents and worm fragments predicted diverse metabolic capabilities with the ability to utilize a range of nitrogen, sulfur, and organic carbon compounds through microbial taxa affiliated with Campylobacterales, Desulfobacterales, Enterobacterales, SAR324, Alphaproteobacteria, and Mycoplasmatales. Entomoplasmatales and Chitinivibrionales were additionally identified from extracted full-length 16S rRNA sequences, and read analysis identified 196 bacterial families. Overall, the microbial community appeared dominated by uncultured Sulfurospirillum, a taxon previously considered free-living rather than host-associated. Metagenome-assembled genomes (MAGs) classified as uncultured Sulfurospirillum predicted thiosulfate disproportionation and the reduction of tetrathionate, sulfate, sulfide/polysulfide, and nitrate. Microbial amino acid and vitamin B12 biosynthesis genes were identified in multiple MAGs, suggesting nutritional value to the host. Reads assigned to aerobic or anaerobic methanotrophic taxa were rare. IMPORTANCE Methane hydrates represent vast reserves of natural gas with roles in global carbon cycling and climate change. This study provided the first analysis of metagenomes associated with Sirsoe methanicola, the only polychaete species known to colonize methane hydrates. Previously unrecognized participation of Sulfurospirillum in a gut microbiome is provided, and the role of sulfur compound redox reactions within this community is highlighted. The comparative biology of S. methanicola is of general interest given research into the adverse effects of sulfide production in human gut microbiomes. In addition, taxonomic assignments are provided for nearly 200 bacterial families, expanding our knowledge of microbiomes in the deep sea.

Influence of nutrient supply on plankton microbiome biodiversity and distribution in a coastal upwelling region.

The ecological and oceanographic processes that drive the response of pelagic ocean microbiomes to environmental changes remain poorly understood, particularly in coastal upwelling ecosystems. Here we show that seasonal and interannual variability in coastal upwelling predicts pelagic ocean microbiome diversity and community structure in the Southern California Current region. Ribosomal RNA gene sequencing, targeting prokaryotic and eukaryotic microbes, from samples collected seasonally during 2014-2020 indicate that nitracline depth is the most robust predictor of spatial microbial community structure and biodiversity in this region. Striking ecological changes occurred due to the transition from a warm anomaly during 2014-2016, characterized by intense stratification, to cooler conditions in 2017-2018, representative of more typical upwelling conditions, with photosynthetic eukaryotes, especially diatoms, changing most strongly. The regional slope of nitracline depth exerts strong control on the relative proportion of highly diverse offshore communities and low biodiversity, but highly productive nearshore communities.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed.

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

The founding charter of the Omic Biodiversity Observation Network (Omic BON).

Omic BON is a thematic Biodiversity Observation Network under the Group on Earth Observations Biodiversity Observation Network (GEO BON), focused on coordinating the observation of biomolecules in organisms and the environment. Our founding partners include representatives from national, regional, and global observing systems; standards organizations; and data and sample management infrastructures. By coordinating observing strategies, methods, and data flows, Omic BON will facilitate the co-creation of a global omics meta-observatory to generate actionable knowledge. Here, we present key elements of Omic BON's founding charter and first activities.

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.

Improving metabarcoding taxonomic assignment: A case study of fishes in a large marine ecosystem.

DNA metabarcoding is an important tool for molecular ecology. However, its effectiveness hinges on the quality of reference sequence databases and classification parameters employed. Here we evaluate the performance of MiFish 12S taxonomic assignments using a case study of California Current Large Marine Ecosystem fishes to determine best practices for metabarcoding. Specifically, we use a taxonomy cross-validation by identity framework to compare classification performance between a global database comprised of all available sequences and a curated database that only includes sequences of fishes from the California Current Large Marine Ecosystem. We demonstrate that the regional database provides higher assignment accuracy than the comprehensive global database. We also document a tradeoff between accuracy and misclassification across a range of taxonomic cutoff scores, highlighting the importance of parameter selection for taxonomic classification. Furthermore, we compared assignment accuracy with and without the inclusion of additionally generated reference sequences. To this end, we sequenced tissue from 597 species using the MiFish 12S primers, adding 252 species to GenBank's existing 550 California Current Large Marine Ecosystem fish sequences. We then compared species and reads identified from seawater environmental DNA samples using global databases with and without our generated references, and the regional database. The addition of new references allowed for the identification of 16 additional native taxa representing 17.0% of total reads from eDNA samples, including species with vast ecological and economic value. Together these results demonstrate the importance of comprehensive and curated reference databases for effective metabarcoding and the need for locus-specific validation efforts.

Red Sea SAR11 and Prochlorococcus Single-Cell Genomes Reflect Globally Distributed Pangenomes.

Evidence suggests many marine bacteria are cosmopolitan, with widespread but sparse strains poised to seed abundant populations under conducive growth conditions. However, studies supporting this "microbial seed bank" hypothesis have analyzed taxonomic marker genes rather than whole genomes/metagenomes, leaving open the possibility that disparate ocean regions harbor endemic gene content. The Red Sea is isolated geographically from the rest of the ocean and has a combination of high irradiance, high temperature, and high salinity that is unique among the oceans; we therefore asked whether it harbors endemic gene content. We sequenced and assembled single-cell genomes of 21 SAR11 (subclades Ia, Ib, Id, and II) and 5 Prochlorococcus (ecotype HLII) samples from the Red Sea and combined them with globally sourced reference genomes to cluster genes into ortholog groups (OGs). Ordination of OG composition could distinguish clades, including phylogenetically cryptic Prochlorococcus ecotypes LLII and LLIII. Compared with reference genomes, 1% of Prochlorococcus and 17% of SAR11 OGs were unique to the Red Sea genomes (RS-OGs). Most (83%) RS-OGs had no annotated function, but 65% of RS-OGs were expressed in diel Red Sea metatranscriptomes, suggesting they are functional. Searching Tara Oceans metagenomes, RS-OGs were as likely to be found as non-RS-OGs; nevertheless, Red Sea and other warm samples could be distinguished from cooler samples using the relative abundances of OGs. The results suggest that the prevalence of OGs in these surface ocean bacteria is largely cosmopolitan, with differences in population metagenomes manifested by differences in relative abundance rather than complete presence/absence of OGs.IMPORTANCE Studies have shown that as we sequence seawater from a selected environment deeper and deeper, we approach finding every bacterial taxon known for the ocean as a whole. However, such studies have focused on taxonomic marker genes rather than on whole genomes, raising the possibility that the lack of endemism results from the method of investigation. We took a geographically isolated water body, the Red Sea, and sequenced single cells from it. We compared those single-cell genomes to available genomes from around the ocean and to ocean-spanning metagenomes. We showed that gene ortholog groups found in Red Sea genomes but not in other genomes are nevertheless common across global ocean metagenomes. These results suggest that Baas Becking's hypothesis "everything is everywhere, but the environment selects" also applies to gene ortholog groups. This widely dispersed functional diversity may give oceanic microbial communities the functional capacity to respond rapidly to changing conditions.

A communal catalogue reveals Earth's multiscale microbial diversity.

Our growing awareness of the microbial world's importance and diversity contrasts starkly with our limited understanding of its fundamental structure. Despite recent advances in DNA sequencing, a lack of standardized protocols and common analytical frameworks impedes comparisons among studies, hindering the development of global inferences about microbial life on Earth. Here we present a meta-analysis of microbial community samples collected by hundreds of researchers for the Earth Microbiome Project. Coordinated protocols and new analytical methods, particularly the use of exact sequences instead of clustered operational taxonomic units, enable bacterial and archaeal ribosomal RNA gene sequences to be followed across multiple studies and allow us to explore patterns of diversity at an unprecedented scale. The result is both a reference database giving global context to DNA sequence data and a framework for incorporating data from future studies, fostering increasingly complete characterization of Earth's microbial diversity.

Regional Assessment of Human Fecal Contamination in Southern California Coastal Drainages.

Host-associated genetic markers that allow for fecal source identification have been used extensively as a diagnostic tool to determine fecal sources within watersheds, but have not been used in routine monitoring to prioritize remediation actions among watersheds. Here, we present a regional assessment of human marker prevalence among drainages that discharge to the U.S. southern California coast. Approximately 50 samples were analyzed for the HF183 human marker from each of 22 southern California coastal drainages under summer dry weather conditions, and another 50 samples were targeted from each of 23 drainages during wet weather. The HF183 marker was ubiquitous, detected in all but two sites in dry weather and at all sites during wet weather. However, there was considerable difference in the extent of human fecal contamination among sites. Similar site ranking was produced regardless of whether the assessment was based on frequency of HF183 detection or site average HF183 concentration. However, site ranking differed greatly between dry and wet weather. Site ranking also differed greatly when based on enterococci, which do not distinguish between pollution sources, vs. HF183, which distinguishes higher risk human fecal sources from other sources, indicating the additional value of the human-associated marker as a routine monitoring tool.

Consideration of Natural Sources in a Bacteria TMDL-Lines of Evidence, Including Beach Microbial Source Tracking.

Total Maximum Daily Load (TMDL) stipulations remained unmet at a southern California beach despite a suite of management actions carried out since 2001, prompting exploration of a Natural Sources Exclusion (NSE) provision within the TMDL. Quantitative Microbial Source Tracking (MST) was employed from 2012 to 2015 to inventory sources of natural and anthropogenic fecal indicator bacteria (FIB). Data suggested FIB exceedances could be traced to gulls based on gull marker prevalence and correlations with FIB concentrations in seawater, sand, and eelgrass. In contrast, human marker concentrations and a tracer dye test did not indicate prevalent human sources. Exponential decay of gull marker in sand amended with live Catellicoccus marimammalium suggested that measured marker reflected fecal inputs versus growth outside the host. Improved water quality was coincident with a 2013 bird exclusion structure, consistent with NSE. However, load allocation needed for TMDL reconsideration was hampered by variable ratios of FIB, MST markers, and pathogens measured in seawater and in gull, cat, and raccoon feces. Quantitative Microbial Risk Assessment is a suggested path forward because such models can incorporate distributions from a combination of FIB sources and communicate criteria in terms of human health risk.

Differential Impacts of Land-Based Sources of Pollution on the Microbiota of Southeast Florida Coral Reefs.

Coral reefs are dynamic ecosystems known for decades to be endangered due, in large part, to anthropogenic impacts from land-based sources of pollution (LBSP). In this study, we utilized an Illumina-based next-generation sequencing approach to characterize prokaryotic and fungal communities from samples collected off the southeast coast of Florida. Water samples from coastal inlet discharges, oceanic outfalls of municipal wastewater treatment plants, treated wastewater effluent before discharge, open ocean samples, and coral tissue samples (mucus and polyps) were characterized to determine the relationships between microbial communities in these matrices and those in reef water and coral tissues. Significant differences in microbial communities were noted among all sample types but varied between sampling areas. Contamination from outfalls was found to be the greatest potential source of LBSP influencing native microbial community structure among all reef samples, although pollution from inlets was also noted. Notably, reef water and coral tissue communities were found to be more greatly impacted by LBSP at southern reefs, which also experienced the most degradation during the course of the study. The results of this study provide new insights into how microbial communities from LBSP can impact coral reefs in southeast Florida and suggest that wastewater outfalls may have a greater influence on the microbial diversity and structure of these reef communities than do contaminants carried in runoff, although the influences of runoff and coastal inlet discharge on coral reefs are still substantial.IMPORTANCE Coral reefs are known to be endangered due to sewage discharge and to runoff of nutrients, pesticides, and other substances associated with anthropogenic activity. Here, we used next-generation sequencing to characterize the microbial communities of potential contaminant sources in order to determine how environmental discharges of microbiota and their genetic material may influence the microbiomes of coral reef communities and coastal receiving waters. Runoff delivered through inlet discharges impacted coral microbial communities, but impacts from oceanic outfalls carrying treated wastewater were greater. Geographic differences in the degree of impact suggest that coral microbiomes may be influenced by the microbiological quality of treated wastewater.

Watershed Assessment with Beach Microbial Source Tracking and Outcomes of Resulting Gull Management.

Total maximum daily load (TMDL) implementation at a southern California beach involved ultraviolet treatment of watershed drainage that provided >97% reduction in fecal indicator bacteria (FIB) concentrations. However, this pollutant control measure did not provide sufficient improvement of beach water quality, prompting further assessment. Investigation included microbial source tracking (MST) for human, gull, and canine fecal sources, monitoring of enterococci and fecal coliform, and measurement of chemical and physical water quality parameters for samples collected from watershed, groundwater, and beach sites, including a beach scour pond and tidal creek. FIB variability remained poorly modeled in regression analysis. However, MST revealed correlations between FIB and gull source tracking markers, leading to recommendations to manage gulls as a pollutant source. Beach conditions were followed for three years after implementation of a best management practice (BMP) to abate gulls using a falconry program for the beach and an upland landfill. The gull abatement BMP was associated with improved beach water quality, and this appears to be the first report of falconry in the context of TMDL implementation. Overall, MST data enabled management action despite an inability to fully model FIB dynamics in the coupled watershed-beach system.

Multi-laboratory evaluations of the performance of Catellicoccus marimammalium PCR assays developed to target gull fecal sources.

Here we report results from a multi-laboratory (n = 11) evaluation of four different PCR methods targeting the 16S rRNA gene of Catellicoccus marimammalium originally developed to detect gull fecal contamination in coastal environments. The methods included a conventional end-point PCR method, a SYBR(®) Green qPCR method, and two TaqMan(®) qPCR methods. Different techniques for data normalization and analysis were tested. Data analysis methods had a pronounced impact on assay sensitivity and specificity calculations. Across-laboratory standardization of metrics including the lower limit of quantification (LLOQ), target detected but not quantifiable (DNQ), and target not detected (ND) significantly improved results compared to results submitted by individual laboratories prior to definition standardization. The unit of measure used for data normalization also had a pronounced effect on measured assay performance. Data normalization to DNA mass improved quantitative method performance as compared to enterococcus normalization. The MST methods tested here were originally designed for gulls but were found in this study to also detect feces from other birds, particularly feces composited from pigeons. Sequencing efforts showed that some pigeon feces from California contained sequences similar to C. marimammalium found in gull feces. These data suggest that the prevalence, geographic scope, and ecology of C. marimammalium in host birds other than gulls require further investigation. This study represents an important first step in the multi-laboratory assessment of these methods and highlights the need to broaden and standardize additional evaluations, including environmentally relevant target concentrations in ambient waters from diverse geographic regions.

Evaluation of the repeatability and reproducibility of a suite of qPCR-based microbial source tracking methods.

Many PCR-based methods for microbial source tracking (MST) have been developed and validated within individual research laboratories. Inter-laboratory validation of these methods, however, has been minimal, and the effects of protocol standardization regimes have not been thoroughly evaluated. Knowledge of factors influencing PCR in different laboratories is vital to future technology transfer for use of MST methods as a tool for water quality management. In this study, a blinded set of 64 filters (containing 32 duplicate samples generated from 12 composite fecal sources) were analyzed by three to five core laboratories with a suite of PCR-based methods utilizing standardized reagents and protocols. Repeatability (intra-laboratory variability) and reproducibility (inter-laboratory variability) of observed results were assessed. When standardized methodologies were used, intra- and inter-laboratory %CVs were generally low (median %CV 0.1-3.3% and 1.9-7.1%, respectively) and comparable to those observed in similar inter-laboratory validation studies performed on other methods of quantifying fecal indicator bacteria (FIB) in environmental samples. ANOVA of %CV values found three human-associated methods (BsteriF1, BacHum, and HF183Taqman) to be similarly reproducible (p > 0.05) and significantly more reproducible (p < 0.05) than HumM2. This was attributed to the increased variability associated with low target concentrations detected by HumM2 (approximately 1-2 log10copies/filter lower) compared to other human-associated methods. Cow-associated methods (BacCow and CowM2) were similarly reproducible (p > 0.05). When using standardized protocols, variance component analysis indicated sample type (fecal source and concentration) to be the major contributor to total variability with that from replicate filters and inter-laboratory analysis to be within the same order of magnitude but larger than inherent intra-laboratory variability. However, when reagents and protocols were not standardized, inter-laboratory %CV generally increased with a corresponding decline in reproducibility. Overall, these findings verify the repeatability and reproducibility of these MST methods and highlight the need for standardization of protocols and consumables prior to implementation of larger scale MST studies involving multiple laboratories.

Performance evaluation of canine-associated Bacteroidales assays in a multi-laboratory comparison study.

The contribution of fecal pollution from dogs in urbanized areas can be significant and is an often underestimated problem. Microbial source tracking methods (MST) utilizing quantitative PCR of dog-associated gene sequences encoding 16S rRNA of Bacteroidales are a useful tool to estimate these contributions. However, data about the performance of available assays are scarce. The results of a multi-laboratory study testing two assays for the determination of dog-associated Bacteroidales (DogBact and BacCan-UCD) on 64 single and mixed fecal source samples created from pooled fecal samples collected in California are presented here. Standardization of qPCR data treatment lowered inter-laboratory variability of sensitivity and specificity results. Both assays exhibited 100% sensitivity. Normalization methods are presented that eliminated random and confirmed non-target responses. The combination of standardized qPCR data treatment, use of normalization via a non-target specific Bacteroidales assay (GenBac3), and application of threshold criteria improved the calculated specificity significantly for both assays. Such measures would reasonably improve MST data interpretation not only for canine-associated assays, but for all qPCR assays used in identifying and monitoring fecal pollution in the environment.

Recommendations following a multi-laboratory comparison of microbial source tracking methods.

Microbial source tracking (MST) methods were evaluated in the Source Identification Protocol Project (SIPP), in which 27 laboratories compared methods to identify host sources of fecal pollution from blinded water samples containing either one or two different fecal types collected from California. This paper details lessons learned from the SIPP study and makes recommendations to further advance the field of MST. Overall, results from the SIPP study demonstrated that methods are available that can correctly identify whether particular host sources including humans, cows and birds have contributed to contamination in a body of water. However, differences between laboratory protocols and data processing affected results and complicated interpretation of MST method performance in some cases. This was an issue particularly for samples that tested positive (non-zero Ct values) but below the limits of quantification or detection of a PCR assay. Although false positives were observed, such samples in the SIPP study often contained the fecal pollution source that was being targeted, i.e., the samples were true positives. Given these results, and the fact that MST often requires detection of targets present in low concentrations, we propose that such samples be reported and identified in a unique category to facilitate data analysis and method comparisons. Important data can be lost when such samples are simply reported as positive or negative. Actionable thresholds were not derived in the SIPP study due to limitations that included geographic scope, age of samples, and difficulties interpreting low concentrations of target in environmental samples. Nevertheless, the results of the study support the use of MST for water management, especially to prioritize impaired waters in need of remediation. Future integration of MST data into quantitative microbial risk assessments and other models could allow managers to more efficiently protect public health based on site conditions.