Effects of Copper on Legionella pneumophila Revealed via Viability Assays and Proteomics.

Cu is an antimicrobial that is commonly applied to premise (i.e., building) plumbing systems for Legionella control, but the precise mechanisms of inactivation are not well defined. Here, we applied a suite of viability assays and mass spectrometry-based proteomics to assess the mechanistic effects of Cu on L. pneumophila. Although a five- to six-log reduction in culturability was observed with 5 mg/L Cu2+ exposure, cell membrane integrity only indicated a <50% reduction. Whole-cell proteomic analysis revealed that AhpD, a protein related to oxidative stress, was elevated in Cu-exposed Legionella relative to culturable cells. Other proteins related to cell membrane synthesis and motility were also higher for the Cu-exposed cells relative to controls without Cu. While the proteins related to primary metabolism decreased for the Cu-exposed cells, no significant differences in the abundance of proteins related to virulence or infectivity were found, which was consistent with the ability of VBNC cells to cause infections. Whereas the cell-membrane integrity assay provided an upper-bound measurement of viability, an amoebae co-culture assay provided a lower-bound limit. The findings have important implications for assessing Legionella risk following its exposure to copper in engineered water systems.

Selection and horizontal gene transfer underlie microdiversity-level heterogeneity in resistance gene fate during wastewater treatment.

Activated sludge is the centerpiece of biological wastewater treatment, as it facilitates removal of sewage-associated pollutants, fecal bacteria, and pathogens from wastewater through semi-controlled microbial ecology. It has been hypothesized that horizontal gene transfer facilitates the spread of antibiotic resistance genes within the wastewater treatment plant, in part because of the presence of residual antibiotics in sewage. However, there has been surprisingly little evidence to suggest that sewage-associated antibiotics select for resistance at wastewater treatment plants via horizontal gene transfer or otherwise. We addressed the role of sewage-associated antibiotics in promoting antibiotic resistance using lab-scale sequencing batch reactors fed field-collected wastewater, metagenomic sequencing, and our recently developed bioinformatic tool Kairos. Here, we found confirmatory evidence that fluctuating levels of antibiotics in sewage are associated with horizontal gene transfer of antibiotic resistance genes, microbial ecology, and microdiversity-level differences in resistance gene fate in activated sludge.

Highly Multiplexed Reverse-Transcription Loop-Mediated Isothermal Amplification and Nanopore Sequencing (LAMPore) for Wastewater-Based Surveillance.

Wastewater-based surveillance (WBS) has gained attention as a strategy to monitor and provide an early warning for disease outbreaks. Here, we applied an isothermal gene amplification technique, reverse-transcription loop-mediated isothermal amplification (RT-LAMP), coupled with nanopore sequencing (LAMPore) as a means to detect SARS-CoV-2. Specifically, we combined barcoding using both an RT-LAMP primer and the nanopore rapid barcoding kit to achieve highly multiplexed detection of SARS-CoV-2 in wastewater. RT-LAMP targeting the SARS-CoV-2 N region was conducted on 96 reactions including wastewater RNA extracts and positive and no-target controls. The resulting amplicons were pooled and subjected to nanopore sequencing, followed by demultiplexing based on barcodes that differentiate the source of each SARS-CoV-2 N amplicon derived from the 96 RT-LAMP products. The criteria developed and applied to establish whether SARS-CoV-2 was detected by the LAMPore assay indicated high consistency with polymerase chain reaction-based detection of the SARS-CoV-2 N gene, with a sensitivity of 89% and a specificity of 83%. We further profiled sequence variations on the SARS-CoV-2 N amplicons, revealing a number of mutations on a sample collected after viral variants had emerged. The results demonstrate the potential of the LAMPore assay to facilitate WBS for SARS-CoV-2 and the emergence of viral variants in wastewater.

mobileOG-db: a Manually Curated Database of Protein Families Mediating the Life Cycle of Bacterial Mobile Genetic Elements.

Bacterial mobile genetic elements (MGEs) encode functional modules that perform both core and accessory functions for the element, the latter of which are often only transiently associated with the element. The presence of these accessory genes, which are often close homologs to primarily immobile genes, incur high rates of false positives and, therefore, limits the usability of these databases for MGE annotation. To overcome this limitation, we analyzed 10,776,849 protein sequences derived from eight MGE databases to compile a comprehensive set of 6,140 manually curated protein families that are linked to the "life cycle" (integration/excision, replication/recombination/repair, transfer, stability/transfer/defense, and phage-specific processes) of plasmids, phages, integrative, transposable, and conjugative elements. We overlay experimental information where available to create a tiered annotation scheme of high-quality annotations and annotations inferred exclusively through bioinformatic evidence. We additionally provide an MGE-class label for each entry (e.g., plasmid or integrative element), and assign to each entry a major and minor category. The resulting database, mobileOG-db (for mobile orthologous groups), comprises over 700,000 deduplicated sequences encompassing five major mobileOG categories and more than 50 minor categories, providing a structured language and interpretable basis for an array of MGE-centered analyses. mobileOG-db can be accessed at mobileogdb.flsi.cloud.vt.edu/, where users can select, refine, and analyze custom subsets of the dynamic mobilome. IMPORTANCE The analysis of bacterial mobile genetic elements (MGEs) in genomic data is a critical step toward profiling the root causes of antibiotic resistance, phenotypic or metabolic diversity, and the evolution of bacterial genera. Existing methods for MGE annotation pose high barriers of biological and computational expertise to properly harness. To bridge this gap, we systematically analyzed 10,776,849 proteins derived from eight databases of MGEs to identify 6,140 MGE protein families that can serve as candidate hallmarks, i.e., proteins that can be used as "signatures" of MGEs to aid annotation. The resulting resource, mobileOG-db, provides a multilevel classification scheme that encompasses plasmid, phage, integrative, and transposable element protein families categorized into five major mobileOG categories and more than 50 minor categories. mobileOG-db thus provides a rich resource for simple and intuitive element annotation that can be integrated seamlessly into existing MGE detection pipelines and colocalization analyses.

Critical evaluation of short, long, and hybrid assembly for contextual analysis of antibiotic resistance genes in complex environmental metagenomes.

In the fight to limit the global spread of antibiotic resistance, the assembly of environmental metagenomes has the potential to provide rich contextual information (e.g., taxonomic hosts, carriage on mobile genetic elements) about antibiotic resistance genes (ARG) in the environment. However, computational challenges associated with assembly can impact the accuracy of downstream analyses. This work critically evaluates the impact of assembly leveraging short reads, nanopore MinION long-reads, and a combination of the two (hybrid) on ARG contextualization for ten environmental metagenomes using seven prominent assemblers (IDBA-UD, MEGAHIT, Canu, Flye, Opera-MS, metaSpades and HybridSpades). While short-read and hybrid assemblies produced similar patterns of ARG contextualization, raw or assembled long nanopore reads produced distinct patterns. Based on an in-silico spike-in experiment using real and simulated reads, we show that low to intermediate coverage species are more likely to be incorporated into chimeric contigs across all assemblers and sequencing technologies, while more abundant species produce assemblies with a greater frequency of inversions and insertion/deletions (indels). In sum, our analyses support hybrid assembly as a valuable technique for boosting the reliability and accuracy of assembly-based analyses of ARGs and neighboring genes at environmentally-relevant coverages, provided that sufficient short-read sequencing depth is achieved.

Whole genome sequence analysis reveals the broad distribution of the RtxA type 1 secretion system and four novel putative type 1 secretion systems throughout the Legionella genus.

Type 1 secretion systems (T1SSs) are broadly distributed among bacteria and translocate effectors with diverse function across the bacterial cell membrane. Legionella pneumophila, the species most commonly associated with Legionellosis, encodes a T1SS at the lssXYZABD locus which is responsible for the secretion of the virulence factor RtxA. Many investigations have failed to detect lssD, the gene encoding the membrane fusion protein of the RtxA T1SS, in non-pneumophila Legionella, which has led to the assumption that this system is a virulence factor exclusively possessed by L. pneumophila. Here we discovered RtxA and its associated T1SS in a novel Legionella taurinensis strain, leading us to question whether this system may be more widespread than previously thought. Through a bioinformatic analysis of publicly available data, we classified and determined the distribution of four T1SSs including the RtxA T1SS and four novel T1SSs among diverse Legionella spp. The ABC transporter of the novel Legionella T1SS Legionella repeat protein secretion system shares structural similarity to those of diverse T1SS families, including the alkaline protease T1SS in Pseudomonas aeruginosa. The Legionella bacteriocin (1-3) secretion systems T1SSs are novel putative bacteriocin transporting T1SSs as their ABC transporters include C-39 peptidase domains in their N-terminal regions, with LB2SS and LB3SS likely constituting a nitrile hydratase leader peptide transport T1SSs. The LB1SS is more closely related to the colicin V T1SS in Escherichia coli. Of 45 Legionella spp. whole genomes examined, 19 (42%) were determined to possess lssB and lssD homologs. Of these 19, only 7 (37%) are known pathogens. There was no difference in the proportions of disease associated and non-disease associated species that possessed the RtxA T1SS (p = 0.4), contrary to the current consensus regarding the RtxA T1SS. These results draw into question the nature of RtxA and its T1SS as a singular virulence factor. Future studies should investigate mechanistic explanations for the association of RtxA with virulence.

Comparison of Whole-Genome Sequences of Legionella pneumophila in Tap Water and in Clinical Strains, Flint, Michigan, USA, 2016.

During the water crisis in Flint, Michigan, USA (2014-2015), 2 outbreaks of Legionnaires' disease occurred in Genesee County, Michigan. We compared whole-genome sequences of 10 clinical Legionella pneumophila isolates submitted to a laboratory in Genesee County during the second outbreak with 103 water isolates collected the following year. We documented a genetically diverse range of L. pneumophila strains across clinical and water isolates. Isolates belonging to 1 clade (3 clinical isolates, 3 water isolates from a Flint hospital, 1 water isolate from a Flint residence, and the reference Paris strain) had a high degree of similarity (2-1,062 single-nucleotide polymorphisms), all L. pneumophila sequence type 1, serogroup 1. Serogroup 6 isolates belonging to sequence type 2518 were widespread in Flint hospital water samples but bore no resemblance to available clinical isolates. L. pneumophila strains in Flint tap water after the outbreaks were diverse and similar to some disease-causing strains.