Olivar: automated variant aware primer design for multiplex tiled amplicon sequencing of pathogens

Tiled amplicon sequencing has served as an essential tool for tracking the spread and evolution of pathogens. Over 2 million complete SARS-CoV-2 genomes are now publicly available, most sequenced and assembled via tiled amplicon sequencing. While computational tools for tiled amplicon design exist, they require downstream manual optimization both computationally and experimentally, which is slow and costly. Here we present Olivar, a first step towards a fully automated, variant-aware design of tiled amplicons for pathogen genomes. Olivar converts each nucleotide of the target genome into a numeric risk score, capturing undesired sequence features that should be avoided. In a direct comparison with PrimalScheme, we show that Olivar has fewer SNPs overlapping with primers and predicted PCR byproducts. We also compared Olivar head-to-head with ARTIC v4.1, the most widely used primer set for SARS-CoV-2 sequencing, and show Olivar yields similar read mapping rates (~90%) and better coverage to the manually designed ARTIC v4.1 amplicons. We also evaluated Olivar on real wastewater samples and found that Olivar had up to 3-fold higher mapping rates while retaining similar coverage. In summary, Olivar automates and accelerates the generation of tiled amplicons, even in situations of high mutation frequency and/or density. Olivar is available as a web application at https://olivar.rice.edu. Olivar can also be installed locally as a command line tool with Bioconda. Source code, installation guide and usage are available at https://github.com/treangenlab/Olivar.

Wastewater surveillance of SARS-CoV-2 and influenza in preK-12 schools shows school, community, and citywide infections

Wastewater surveillance is a passive and efficient way to monitor the spread of infectious diseases in large populations and high transmission areas such as preK-12 schools. Infections caused by respiratory viruses in school-aged children are likely underreported, particularly because many children may be asymptomatic or mildly symptomatic. Wastewater monitoring of SARS-CoV-2 has been studied extensively and primarily by sampling at centralized wastewater treatment plants, and there are limited studies on SARS-CoV-2 in preK-12 school wastewater. Similarly, wastewater detections of influenza have only been reported in wastewater treatment plant and university manhole samples. Here, we present the results of a 17-month wastewater monitoring program for SARS-CoV-2 (n = 2176 samples) and influenza A and B (n = 1217 samples) in 51 preK-12 schools. We show that school wastewater concentrations of SARS-CoV-2 RNA were strongly associated with COVID-19 cases in schools and community positivity rates, and that influenza detections in school wastewater were significantly associated with citywide influenza diagnosis rates. Results were communicated back to schools and local communities to enable mitigation strategies to stop the spread, and direct resources such as testing and vaccination clinics. This study demonstrates that school wastewater surveillance is reflective of local infections at several population levels and plays a crucial role in the detection and mitigation of outbreaks.

Direct comparison of RT-ddPCR and targeted amplicon sequencing for SARS-CoV-2 mutation monitoring in wastewater

Over the course of the COVID-19 pandemic, variants of SARS-CoV-2 have emerged that are more contagious and more likely to cause breakthrough infections. Targeted amplicon sequencing approach is a gold standard for identification and analysis of variants. However, when applied to environmental samples such as wastewater, it remains unclear how sensitive this method is for detecting variant-associated mutations in environmental samples. Here we directly compare a targeted amplicon sequencing approach (using ARTIC v3; hereafter referred to as sequencing) with RT-ddPCR quantification for the detection of five mutations that are characteristic of variants of concern (VoCs) in wastewater samples. In total, 547 wastewater samples were analyzed using both methods in parallel. When we observed positive mutation detections by RT-ddPCR, 42.6% of the detection events were missed by sequencing, due to negative detection or the limited read coverage at the mutation position. Further, when sequencing reported negative or depth-limited mutation detections, 26.7% of those events were instead positive detections by RT-ddPCR, highlighting the relatively poor sensitivity of sequencing. No or weak associations were observed between quantitative measurements of target mutations determined by RT-ddPCR and sequencing. These findings caution the use of quantitative measurements of SARS-CoV-2 variants in wastewater samples determined solely based on sequencing.

Direct comparison of RT-ddPCR and targeted amplicon sequencing for SARS-CoV-2 mutation monitoring in wastewater

Over the course of the COVID-19 pandemic, variants of SARS-CoV-2 have emerged that are more contagious and more likely to cause breakthrough infections. Targeted amplicon sequencing approach is a gold standard for identification and analysis of variants. However, when applied to environmental samples such as wastewater, it remains unclear how sensitive this method is for detecting variant-associated mutations in environmental samples. Here we directly compare a targeted amplicon sequencing approach (using ARTIC v3; hereafter referred to as sequencing) with RT-ddPCR quantification for the detection of five mutations that are characteristic of variants of concern (VoCs) in wastewater samples. In total, 547 wastewater samples were analyzed using both methods in parallel. When we observed positive mutation detections by RT-ddPCR, 42.6% of the detection events were missed by sequencing, due to negative detection or the limited read coverage at the mutation position. Further, when sequencing reported negative or depth-limited mutation detections, 26.7% of those events were instead positive detections by RT-ddPCR, highlighting the relatively poor sensitivity of sequencing. No or weak associations were observed between quantitative measurements of target mutations determined by RT-ddPCR and sequencing. These findings caution the use of quantitative measurements of SARS-CoV-2 variants in wastewater samples determined solely based on sequencing.

Evaluating recovery, cost, and throughput of different concentration methods for SARS-CoV-2 wastewater-based epidemiology

As the COVID-19 pandemic continues to affect communities across the globe, the need to contain the spread of the outbreaks is of paramount importance. Wastewater monitoring of the SARS-CoV-2 virus, the causative agent responsible for COVID-19, has emerged as a promising tool for health officials to anticipate outbreaks. As interest in wastewater monitoring continues to grow and municipalities begin to implement this approach, there is a need to further identify and evaluate methods used to concentrate SARS-CoV-2 virus RNA from wastewater samples. Here we evaluate the recovery, cost, and throughput of five different concentration methods for quantifying SARS-CoV-2 virus RNA in wastewater samples. We tested the five methods on six different wastewater samples. We also evaluated the use of a bovine coronavirus vaccine as a process control and pepper mild mottle virus as a normalization factor. Of the five methods we tested head-to-head, we found that HA filtration with bead beating performed the best in terms of sensitivity and cost. This evaluation can serve as a guide for laboratories establishing a protocol to perform wastewater monitoring of SARS-CoV-2.

Translating New Synthetic Biology Advances for Biosensing Into the Earth and Environmental Sciences

The rapid diversification of synthetic biology tools holds promise in making some classically hard-to-solve environmental problems tractable. Here we review longstanding problems in the Earth and environmental sciences that could be addressed using engineered microbes as micron-scale sensors (biosensors). Biosensors can offer new perspectives on open questions, including understanding microbial behaviors in heterogeneous matrices like soils, sediments, and wastewater systems, tracking cryptic element cycling in the Earth system, and establishing the dynamics of microbe-microbe, microbe-plant, and microbe-material interactions. Before these new tools can reach their potential, however, a suite of biological parts and microbial chassis appropriate for environmental conditions must be developed by the synthetic biology community. This includes diversifying sensing modules to obtain information relevant to environmental questions, creating output signals that allow dynamic reporting from hard-to-image environmental materials, and tuning these sensors so that they reliably function long enough to be useful for environmental studies. Finally, ethical questions related to the use of synthetic biosensors in environmental applications are discussed.

Using cellular fitness to map the structure and function of a major facilitator superfamily effluxer

The major facilitator superfamily (MFS) effluxers are prominent mediators of antimicrobial resistance. The biochemical characterization of MFS proteins is hindered by their complex membrane environment that makes in vitro biochemical analysis challenging. Since the physicochemical properties of proteins drive the fitness of an organism, we posed the question of whether we could reverse that relationship and derive meaningful biochemical parameters for a single protein simply from fitness changes it confers under varying strengths of selection. Here, we present a physiological model that uses cellular fitness as a proxy to predict the biochemical properties of the MFS tetracycline efflux pump, TetB, and a family of single amino acid variants. We determined two lumped biochemical parameters roughly describing K_m and V_max for TetB and variants. Including in vivo protein levels into our model allowed for more specified prediction of pump parameters relating to substrate binding affinity and pumping efficiency for TetB and variants. We further demonstrated the general utility of our model by solely using fitness to assay a library of tet(B) variants and estimate their biochemical properties.