Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data.

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic 'novel' lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

Diverse and abundant phages exploit conjugative plasmids.

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 65 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the evolutionarily promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

Changing fitness effects of mutations through long-term bacterial evolution.

The distribution of fitness effects of new mutations shapes evolution, but it is challenging to observe how it changes as organisms adapt. Using Escherichia coli lineages spanning 50,000 generations of evolution, we quantify the fitness effects of insertion mutations in every gene. Macroscopically, the fraction of deleterious mutations changed little over time whereas the beneficial tail declined sharply, approaching an exponential distribution. Microscopically, changes in individual gene essentiality and deleterious effects often occurred in parallel; altered essentiality is only partly explained by structural variation. The identity and effect sizes of beneficial mutations changed rapidly over time, but many targets of selection remained predictable because of the importance of loss-of-function mutations. Taken together, these results reveal the dynamic-but statistically predictable-nature of mutational fitness effects.

Efficient and Robust Search of Microbial Genomes via Phylogenetic Compression.

Comprehensive collections approaching millions of sequenced genomes have become central information sources in the life sciences. However, the rapid growth of these collections makes it effectively impossible to search these data using tools such as BLAST and its successors. Here, we present a technique called phylogenetic compression, which uses evolutionary history to guide compression and efficiently search large collections of microbial genomes using existing algorithms and data structures. We show that, when applied to modern diverse collections approaching millions of genomes, lossless phylogenetic compression improves the compression ratios of assemblies, de Bruijn graphs, and k-mer indexes by one to two orders of magnitude. Additionally, we develop a pipeline for a BLAST-like search over these phylogeny-compressed reference data, and demonstrate it can align genes, plasmids, or entire sequencing experiments against all sequenced bacteria until 2019 on ordinary desktop computers within a few hours. Phylogenetic compression has broad applications in computational biology and may provide a fundamental design principle for future genomics infrastructure.

Resolving Deleterious and Near-Neutral Effects Requires Different Pooled Fitness Assay Designs.

Pooled sequencing-based fitness assays are a powerful and widely used approach to quantifying fitness of thousands of genetic variants in parallel. Despite the throughput of such assays, they are prone to biases in fitness estimates, and errors in measurements are typically larger for deleterious fitness effects, relative to neutral effects. In practice, designing pooled fitness assays involves tradeoffs between the number of timepoints, the sequencing depth, and other parameters to gain as much information as possible within a feasible experiment. Here, we combined simulations and reanalysis of an existing experimental dataset to explore how assay parameters impact measurements of near-neutral and deleterious fitness effects using a standard fitness estimator. We found that sequencing multiple timepoints at relatively modest depth improved estimates of near-neutral fitness effects, but systematically biased measurements of deleterious effects. We showed that a fixed total number of reads, deeper sequencing at fewer timepoints improved resolution of deleterious fitness effects. Our results highlight a tradeoff between measurement of deleterious and near-neutral effect sizes for a fixed amount of data and suggest that fitness assay design should be tuned for fitness effects that are relevant to the specific biological question.

Diverse and abundant phages exploit conjugative plasmids.

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally-encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 64 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

A swapped genetic code prevents viral infections and gene transfer.

Engineering the genetic code of an organism has been proposed to provide a firewall from natural ecosystems by preventing viral infections and gene transfer1-6. However, numerous viruses and mobile genetic elements encode parts of the translational apparatus7-9, potentially rendering a genetic-code-based firewall ineffective. Here we show that such mobile transfer RNAs (tRNAs) enable gene transfer and allow viral replication in Escherichia coli despite the genome-wide removal of 3 of the 64 codons and the previously essential cognate tRNA and release factor genes. We then establish a genetic firewall by discovering viral tRNAs that provide exceptionally efficient codon reassignment allowing us to develop cells bearing an amino acid-swapped genetic code that reassigns two of the six serine codons to leucine during translation. This amino acid-swapped genetic code renders cells resistant to viral infections by mistranslating viral proteomes and prevents the escape of synthetic genetic information by engineered reliance on serine codons to produce leucine-requiring proteins. As these cells may have a selective advantage over wild organisms due to virus resistance, we also repurpose a third codon to biocontain this virus-resistant host through dependence on an amino acid not found in nature10. Our results may provide the basis for a general strategy to make any organism safely resistant to all natural viruses and prevent genetic information flow into and out of genetically modified organisms.

Mycobacterial nucleoid-associated protein Lsr2 is required for productive mycobacteriophage infection.

Mycobacteriophages are a diverse group of viruses infecting Mycobacterium with substantial therapeutic potential. However, as this potential becomes realized, the molecular details of phage infection and mechanisms of resistance remain ill-defined. Here we use live-cell fluorescence microscopy to visualize the spatiotemporal dynamics of mycobacteriophage infection in single cells and populations, showing that infection is dependent on the host nucleoid-associated Lsr2 protein. Mycobacteriophages preferentially adsorb at Mycobacterium smegmatis sites of new cell wall synthesis and following DNA injection, Lsr2 reorganizes away from host replication foci to establish zones of phage DNA replication (ZOPR). Cells lacking Lsr2 proceed through to cell lysis when infected but fail to generate consecutive phage bursts that trigger epidemic spread of phage particles to neighbouring cells. Many mycobacteriophages code for their own Lsr2-related proteins, and although their roles are unknown, they do not rescue the loss of host Lsr2.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Decreased thermal niche breadth as a trade-off of antibiotic resistance.

Evolutionary theory predicts that adaptations, including antibiotic resistance, should come with associated fitness costs; yet, many resistance mutations seemingly contradict this prediction by inducing no growth rate deficit. However, most growth assays comparing sensitive and resistant strains have been performed under a narrow range of environmental conditions, which do not reflect the variety of contexts that a pathogenic bacterium might encounter when causing infection. We hypothesized that reduced niche breadth, defined as diminished growth across a diversity of environments, can be a cost of antibiotic resistance. Specifically, we test whether chloramphenicol-resistant Escherichia coli incur disproportionate growth deficits in novel thermal conditions. Here we show that chloramphenicol-resistant bacteria have greater fitness costs at novel temperatures than their antibiotic-sensitive ancestors. In several cases, we observed no resistance cost in growth rate at the historic temperature but saw diminished growth at warmer and colder temperatures. These results were consistent across various genetic mechanisms of resistance. Thus, we propose that decreased thermal niche breadth is an under-documented fitness cost of antibiotic resistance. Furthermore, these results demonstrate that the cost of antibiotic resistance shifts rapidly as the environment changes; these context-dependent resistance costs should select for the rapid gain and loss of resistance as an evolutionary strategy.

The bacterial toxin colibactin triggers prophage induction.

Colibactin is a chemically unstable small-molecule genotoxin that is produced by several different bacteria, including members of the human gut microbiome1,2. Although the biological activity of colibactin has been extensively investigated in mammalian systems3, little is known about its effects on other microorganisms. Here we show that colibactin targets bacteria that contain prophages, and induces lytic development through the bacterial SOS response. DNA, added exogenously, protects bacteria from colibactin, as does expressing a colibactin resistance protein (ClbS) in non-colibactin-producing cells. The prophage-inducing effects that we observe apply broadly across different phage-bacteria systems and in complex communities. Finally, we identify bacteria that have colibactin resistance genes but lack colibactin biosynthetic genes. Many of these bacteria are infected with predicted prophages, and we show that the expression of their ClbS homologues provides immunity from colibactin-triggered induction. Our study reveals a mechanism by which colibactin production could affect microbiomes and highlights a role for microbial natural products in influencing population-level events such as phage outbreaks.

Prophages encode phage-defense systems with cognate self-immunity.

Temperate phages are pervasive in bacterial genomes, existing as vertically inherited islands termed prophages. Prophages are vulnerable to predation of their host bacterium by exogenous phages. Here, we identify BstA, a family of prophage-encoded phage-defense proteins in diverse Gram-negative bacteria. BstA localizes to sites of exogenous phage DNA replication and mediates abortive infection, suppressing the competing phage epidemic. During lytic replication, the BstA-encoding prophage is not itself inhibited by BstA due to self-immunity conferred by the anti-BstA (aba) element, a short stretch of DNA within the bstA locus. Inhibition of phage replication by distinct BstA proteins from Salmonella, Klebsiella, and Escherichia prophages is generally interchangeable, but each possesses a cognate aba element. The specificity of the aba element ensures that immunity is exclusive to the replicating prophage, preventing exploitation by variant BstA-encoding phages. The BstA protein allows prophages to defend host cells against exogenous phage attack without sacrificing the ability to replicate lytically.

Variant abundance estimation for SARS-CoV-2 in wastewater using RNA-Seq quantification.

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.

Identification of a pathway for electron uptake in Shewanella oneidensis.

Extracellular electron transfer (EET) could enable electron uptake into microbial metabolism for the synthesis of complex, energy dense organic molecules from CO2 and renewable electricity1-6. Theoretically EET could do this with an efficiency comparable to H2-oxidation7,8 but without the need for a volatile intermediate and the problems it causes for scale up9. However, significant gaps remain in understanding the mechanism and genetics of electron uptake. For example, studies of electron uptake in electroactive microbes have shown a role for the Mtr EET complex in the electroactive microbe Shewanella oneidensis MR-110-14, though there is substantial variation in the magnitude of effect deletion of these genes has depending on the terminal electron acceptor used. This speaks to the potential for previously uncharacterized and/or differentially utilized genes involved in electron uptake. To address this, we screened gene disruption mutants for 3667 genes, representing ≈99% of all nonessential genes, from the S. oneidensis whole genome knockout collection using a redox dye oxidation assay. Confirmation of electron uptake using electrochemical testing allowed us to identify five genes from S. oneidensis that are indispensable for electron uptake from a cathode. Knockout of each gene eliminates extracellular electron uptake, yet in four of the five cases produces no significant defect in electron donation to an anode. This result highlights both distinct electron uptake components and an electronic connection between aerobic and anaerobic electron transport chains that allow electrons from the reversible EET machinery to be coupled to different respiratory processes in S. oneidensis. Homologs to these genes across many different genera suggesting that electron uptake by EET coupled to respiration could be widespread. These gene discoveries provide a foundation for: studying this phenotype in exotic metal-oxidizing microbes, genetic optimization of electron uptake in S. oneidensis; and genetically engineering electron uptake into a highly tractable host like E. coli to complement recent advances in synthetic CO2 fixation15.

Simplitigs as an efficient and scalable representation of de Bruijn graphs.

de Bruijn graphs play an essential role in bioinformatics, yet they lack a universal scalable representation. Here, we introduce simplitigs as a compact, efficient, and scalable representation, and ProphAsm, a fast algorithm for their computation. For the example of assemblies of model organisms and two bacterial pan-genomes, we compare simplitigs to unitigs, the best existing representation, and demonstrate that simplitigs provide a substantial improvement in the cumulative sequence length and their number. When combined with the commonly used Burrows-Wheeler Transform index, simplitigs reduce memory, and index loading and query times, as demonstrated with large-scale examples of GenBank bacterial pan-genomes.

Perfect as the enemy of good: tracing transmissions with low-sensitivity tests to mitigate SARS-CoV-2 outbreaks.

Throughout the COVID-19 pandemic, governments and individuals have attempted a wide variety of strategies to limit the damage of the pandemic on human lives, population health, and economies. Contact tracing has been a commonly used strategy, and various approaches have been proposed and attempted. We summarise some methods of contact tracing and testing, considering the resources demanded by each and how features of SARS-CoV-2 transmission affect their effectiveness. We also propose an approach focusing on tracing transmission events, which can be particularly effective when superspreading events play a large role in transmission. Accounting for the best available evidence on a pathogen and for the availability of resources can make control strategies more effective, even if they are not perfect.

Metastable hybridization-based DNA information storage to allow rapid and permanent erasure.

The potential of DNA as an information storage medium is rapidly growing due to advances in DNA synthesis and sequencing. However, the chemical stability of DNA challenges the complete erasure of information encoded in DNA sequences. Here, we encode information in a DNA information solution, a mixture of true message- and false message-encoded oligonucleotides, and enables rapid and permanent erasure of information. True messages are differentiated by their hybridization to a "truth marker" oligonucleotide, and only true messages can be read; binding of the truth marker can be effectively randomized even with a brief exposure to the elevated temperature. We show 8 separate bitmap images can be stably encoded and read after storage at 25 °C for 65 days with an average of over 99% correct information recall, which extrapolates to a half-life of over 15 years at 25 °C. Heating to 95 °C for 5 minutes, however, permanently erases the message.

Barcoded microbial system for high-resolution object provenance.

Determining where an object has been is a fundamental challenge for human health, commerce, and food safety. Location-specific microbes in principle offer a cheap and sensitive way to determine object provenance. We created a synthetic, scalable microbial spore system that identifies object provenance in under 1 hour at meter-scale resolution and near single-spore sensitivity and can be safely introduced into and recovered from the environment. This system solves the key challenges in object provenance: persistence in the environment, scalability, rapid and facile decoding, and biocontainment. Our system is compatible with SHERLOCK, a Cas13a RNA-guided nucleic acid detection assay, facilitating its implementation in a wide range of applications.

Escape mutations circumvent a tradeoff between resistance to a beta-lactam and resistance to a beta-lactamase inhibitor.

Beta-lactamase inhibitors are increasingly used to counteract antibiotic resistance mediated by beta-lactamase enzymes. These inhibitors compete with the beta-lactam antibiotic for the same binding site on the beta-lactamase, thus generating an evolutionary tradeoff: mutations that increase the enzyme's beta-lactamase activity tend to increase also its susceptibility to the inhibitor. Here, we investigate how common and accessible are mutants that escape this adaptive tradeoff. Screening a deep mutant library of the blaampC beta-lactamase gene of Escherichia coli, we identified mutations that allow growth at beta-lactam concentrations far exceeding those inhibiting growth of the wildtype strain, even in the presence of the enzyme inhibitor (avibactam). These escape mutations are rare and drug-specific, and some combinations of avibactam with beta-lactam drugs appear to prevent such escape phenotypes. Our results, showing differential adaptive potential of blaampC to combinations of avibactam and different beta-lactam antibiotics, suggest that it may be possible to identify treatments that are more resilient to evolution of resistance.

Rapid inference of antibiotic resistance and susceptibility by genomic neighbour typing.

Surveillance of drug-resistant bacteria is essential for healthcare providers to deliver effective empirical antibiotic therapy. However, traditional molecular epidemiology does not typically occur on a timescale that could affect patient treatment and outcomes. Here, we present a method called 'genomic neighbour typing' for inferring the phenotype of a bacterial sample by identifying its closest relatives in a database of genomes with metadata. We show that this technique can infer antibiotic susceptibility and resistance for both Streptococcus pneumoniae and Neisseria gonorrhoeae. We implemented this with rapid k-mer matching, which, when used on Oxford Nanopore MinION data, can run in real time. This resulted in the determination of resistance within 10 min (91% sensitivity and 100% specificity for S. pneumoniae and 81% sensitivity and 100% specificity for N. gonorrhoeae from isolates with a representative database) of starting sequencing, and within 4 h of sample collection (75% sensitivity and 100% specificity for S. pneumoniae) for clinical metagenomic sputum samples. This flexible approach has wide application for pathogen surveillance and may be used to greatly accelerate appropriate empirical antibiotic treatment.