MicrobeAnnotator: a user-friendly, comprehensive functional annotation pipeline for microbial genomes.

BACKGROUND: High-throughput sequencing has increased the number of available microbial genomes recovered from isolates, single cells, and metagenomes. Accordingly, fast and comprehensive functional gene annotation pipelines are needed to analyze and compare these genomes. Although several approaches exist for genome annotation, these are typically not designed for easy incorporation into analysis pipelines, do not combine results from different annotation databases or offer easy-to-use summaries of metabolic reconstructions, and typically require large amounts of computing power for high-throughput analysis not available to the average user. RESULTS: Here, we introduce MicrobeAnnotator, a fully automated, easy-to-use pipeline for the comprehensive functional annotation of microbial genomes that combines results from several reference protein databases and returns the matching annotations together with key metadata such as the interlinked identifiers of matching reference proteins from multiple databases [KEGG Orthology (KO), Enzyme Commission (E.C.), Gene Ontology (GO), Pfam, and InterPro]. Further, the functional annotations are summarized into Kyoto Encyclopedia of Genes and Genomes (KEGG) modules as part of a graphical output (heatmap) that allows the user to quickly detect differences among (multiple) query genomes and cluster the genomes based on their metabolic similarity. MicrobeAnnotator is implemented in Python 3 and is freely available under an open-source Artistic License 2.0 from https://github.com/cruizperez/MicrobeAnnotator . CONCLUSIONS: We demonstrated the capabilities of MicrobeAnnotator by annotating 100 Escherichia coli and 78 environmental Candidate Phyla Radiation (CPR) bacterial genomes and comparing the results to those of other popular tools. We showed that the use of multiple annotation databases allows MicrobeAnnotator to recover more annotations per genome compared to faster tools that use reduced databases and is computationally efficient for use in personal computers. The output of MicrobeAnnotator can be easily incorporated into other analysis pipelines while the results of other annotation tools can be seemingly incorporated into MicrobeAnnotator to generate summary plots.

Does age matter? Sexual event-level analysis of age-disparate sexual partners among gay, bisexual and other men who have sex with men (GBM) in Vancouver, Canada.

BACKGROUND: To determine factors associated with age-disparate sexual partners among Vancouver gay, bisexual and other men who have sex with men (GBM). METHODS: Sexually active GBM aged ≥16 years were recruited from February 2012 to February 2014. Participants self-completed a questionnaire on demographics, attitudes and sexual behaviour and substance use at last sexual event with five most recent partners. Two generalised linear mixed models identified factors associated with: (1) 'same-age' (referent), 'younger' or 'much-younger' and (2) 'same-age' (referent), 'older' or 'much-older' partners. Statistical interactions between age and HIV status were tested. RESULTS: Participants (n=719) were predominantly gay (85.1%), White (75.0%), HIV-negative/unknown status (72.9%) with median age of 33 years (Q1,Q3: 26,47). A minority of sexual events were reported with much-older/much-younger partners (13.7%). In the multivariable models, GBM reporting older partners were more likely to be Asian or Latino, have greater Escape Motivation scores, report their partner used erectile dysfunction drugs (EDDs) and have received something for sex; compared with condom-protected insertive anal sex, participants with older partners were more likely to report condomless insertive anal sex with a serodiscordant or unknown status partner or no insertive anal sex. GBM reporting older partners were less likely to be bisexual-identified, have given something for sex and report event-level alcohol and EDD use. GBM reporting younger partners were more likely to have annual incomes >$30 000 and have met their partner online. As per significant statistical interactions, age-disparate relations were more common for younger HIV-positive and older HIV-negative GBM. CONCLUSIONS: Differences among age-disparate partners highlight important targets for health promotion and future research.

Awareness and use of nonoccupational post-exposure prophylaxis among men who have sex with men in Vancouver, Canada.

OBJECTIVES: Nonoccupational post-exposure prophylaxis (nPEP) is a strategy to reduce the risk of HIV infection in those with high-risk exposure. This study characterized nPEP awareness among gay, bisexual and other men who have sex with men (MSM) in Metro Vancouver, British Columbia, Canada after a pilot nPEP programme established in 2012. METHODS: Momentum Health Study participants were MSM aged ≥16 years recruited via respondent-driven sampling (RDS) who completed a computer-assisted self-interview. Stratifying patients by HIV status, we used multivariable logistic regression with backward selection to identify factors associated with nPEP awareness. All analyses were RDS-adjusted. RESULTS: A total of 51.9% (112 of 173) of HIV-positive and 48.5% (272 of 500) of HIV-negative participants had heard of nPEP. Only 3% (five of 106) of HIV-negative participants who reported recent high-risk sex used nPEP. Generally, nPEP awareness was higher for participants who engaged in sexual activities with increased HIV transmission potential. Factors associated with greater awareness among HIV-negative participants included recent alcohol use, higher communal sexual altruism, previous sexually transmitted infection diagnosis, and greater perceived condom use self-efficacy. Other factors associated with greater awareness among HIV-negative participants included white race/ethnicity, gay sexual identity, more formal education, lower personal sexual altruism, and Vancouver residence. Greater nPEP awareness among HIV-positive participants was associated with greater perceived agency to ask sexual partners' HIV status and more frequently reporting doing so, a higher number of lifetime receptive sex partners, and greater access to condoms. CONCLUSIONS: Following implementation of an nPEP pilot programme, nPEP awareness among HIV-negative MSM was 51% and use was 3%. These data support the need to expand access to and actively promote nPEP services.

Branched isomeric 1,2,3-triazolium-based ionic liquids: new insight into structure-property relationships.

A series of four isomeric 1,2,3-triazolium-based ionic liquids (ILs) with vary degree of branching were synthesized and characterized to investigate the effect of ion branching on thermal and physical properties of the resulting IL. It was found that increased branching led to a higher ionicity and higher viscosity. The thermal properties were also altered significantly and spectral changes in the near edge X-ray absorption fine structure (NEXAFS) spectra show that branching affects intermolecular interaction. While the ionicity and viscosity varying linearly with branching, the MDSC and NEXAFS measurements show that the cation shape has a stronger influence on the melting temperature and absorptive properties than the number of branched alkyl substituents.

A natural ANI gap that can define intra-species units of bacteriophages and other viruses.

Despite the importance of intra-species variants of viruses for causing disease and/or disrupting ecosystem functioning, there is no universally applicable standard to define these. A (natural) gap in whole-genome average nucleotide identity (ANI) values around 95% is commonly used to define species, especially for bacteriophages, but whether a similar gap exists within species that can be used to define intra-species units has not been evaluated yet. Whole-genome comparisons among members of 1,016 bacteriophage (Caudoviricetes) species revealed a region of low frequency of ANI values around 99.2%-99.8%, showing threefold or fewer pairs than expected for an even distribution. This second gap is prevalent in viruses infecting various cultured or uncultured hosts from a variety of environments, although a few exceptions to this pattern were also observed (3.7% of total species) and are likely attributed to cultivation biases or other factors. Similar results were observed for a limited set of eukaryotic viruses that are adequately sampled, including SARS-CoV-2, whose ANI-based clusters matched well with the WHO-defined variants of concern, indicating that our findings from bacteriophages might be more broadly applicable and the ANI-based clusters may represent functionally and/or ecologically distinct units. These units appear to be predominantly driven by (high) ecological cohesiveness coupled to either frequent recombination for bacteriophages or selection and clonal evolution for other viruses such as SARS-CoV-2, indicating that fundamentally different underlying mechanisms could lead to similar diversity patterns. Accordingly, we propose the ANI gap approach outlined above for defining viral intra-species units, for which we propose the term genomovars. IMPORTANCE: Viral species are composed of an ensemble of intra-species variants whose individual dynamics may have major implications for human and animal health and/or ecosystem functioning. However, the lack of universally accepted standards to define these intra-species variants has led researchers to use different approaches for this task, creating inconsistent intra-species units across different viral families and confusion in communication. By comparing hundreds of mostly bacteriophage genomes, we show that there is a widely distributed natural gap in whole-genome average nucleotide identity values in most, but not all, of these species that can be used to define intra-species units. Therefore, these results advance the molecular toolbox for tracking viral intra-species units and should facilitate future epidemiological and environmental studies.

Potential routes of plastics biotransformation involving novel plastizymes revealed by global multi-omic analysis of plastic associated microbes.

Despite increasing efforts across various disciplines, the fate, transport, and impact of synthetic plastics on the environment and public health remain poorly understood. To better elucidate the microbial ecology of plastic waste and its potential for biotransformation, we conducted a large-scale analysis of all publicly available meta-omic studies investigating plastics (n = 27) in the environment. Notably, we observed low prevalence of known plastic degraders throughout most environments, except for substantial enrichment in riverine systems. This indicates rivers may be a highly promising environment for discovery of novel plastic bioremediation products. Ocean samples associated with degrading plastics showed clear differentiation from non-degrading polymers, showing enrichment of novel putative biodegrading taxa in the degraded samples. Regarding plastisphere pathogenicity, we observed significant enrichment of antimicrobial resistance genes on plastics but not of virulence factors. Additionally, we report a co-occurrence network analysis of 10 + million proteins associated with the plastisphere. This analysis revealed a localized sub-region enriched with known and putative plastizymes-these may be useful for deeper investigation of nature's ability to biodegrade man-made plastics. Finally, the combined data from our meta-analysis was used to construct a publicly available database, the Plastics Meta-omic Database (PMDB)-accessible at plasticmdb.org. These data should aid in the integrated exploration of the microbial plastisphere and facilitate research efforts investigating the fate and bioremediation potential of environmental plastic waste.

An ANI gap within bacterial species that advances the definitions of intra-species units.

IMPORTANCE: Bacterial strains and clonal complexes are two cornerstone concepts for microbiology that remain loosely defined, which confuses communication and research. Here we identify a natural gap in genome sequence comparisons among isolate genomes of all well-sequenced species that has gone unnoticed so far and could be used to more accurately and precisely define these and related concepts compared to current methods. These findings advance the molecular toolbox for accurately delineating and following the important units of diversity within prokaryotic species and thus should greatly facilitate future epidemiological and micro-diversity studies across clinical and environmental settings.

Towards estimating the number of strains that make up a natural bacterial population.

What a strain is and how many strains make up a natural bacterial population remain elusive concepts despite their apparent importance for assessing the role of intra-population diversity in disease emergence or response to environmental perturbations. To advance these concepts, we sequenced 138 randomly selected Salinibacter ruber isolates from two solar salterns and assessed these genomes against companion short-read metagenomes from the same samples. The distribution of genome-aggregate average nucleotide identity (ANI) values among these isolates revealed a bimodal distribution, with four-fold lower occurrence of values between 99.2% and 99.8% relative to ANI >99.8% or <99.2%, revealing a natural "gap" in the sequence space within species. Accordingly, we used this ANI gap to define genomovars and a higher ANI value of >99.99% and shared gene-content >99.0% to define strains. Using these thresholds and extrapolating from how many metagenomic reads each genomovar uniquely recruited, we estimated that -although our 138 isolates represented about 80% of the Sal. ruber population- the total population in one saltern pond is composed of 5,500 to 11,000 genomovars, the great majority of which appear to be rare in-situ. These data also revealed that the most frequently recovered isolate in lab media was often not the most abundant genomovar in-situ, suggesting that cultivation biases are significant, even in cases that cultivation procedures are thought to be robust. The methodology and ANI thresholds outlined here should represent a useful guide for future microdiversity surveys of additional microbial species.

Radiosensitization by histone deacetylase inhibition in an osteosarcoma mouse model.

BACKGROUND: Osteosarcomas (OS) are highly malignant and radioresistant tumors. Histone deacetylase inhibitors (HDACi) constitute a novel class of anticancer agents. We sought to investigate the effect of combined treatment with suberoylanilide hydroxamic acid (SAHA) and radiotherapy in OS in vivo. METHODS: Clonogenic survival of human OS cell lines as well as tumor growth delay of OS xenografts were tested after treatment with either vehicle, radiotherapy (XRT), SAHA, or XRT and SAHA. Tumor proliferation, necrosis, microvascular density, apoptosis, and p53/p21 were monitored by immunohistochemistry. The CD95 pathway was performed by flow cytometry, caspase (3/7/8) activity measurements, and functional inhibition of CD95 death signaling. RESULTS: Combined treatment with SAHA and XRT markedly reduced the surviving fraction of OS cells as compared to XRT alone. Likewise, dual therapy significantly inhibited OS tumor growth in vivo as compared to XRT alone, reflected by reduced tumor proliferation, impaired angiogenesis, and increased apoptosis. Addition of HDACi to XRT led to elevated p53, p21, CD95, and CD95L expression. Inhibition of CD95 signaling reduced HDACi- and XRT-induced apoptosis. CONCLUSION: Our data show that HDACi increases the radiosensitivity of osteosarcoma cells at least in part via ligand-induced apoptosis. HDACi thus emerge as potentially useful treatment components of OS.

Identification and characterization of the gene causing type 1 spinocerebellar ataxia.

Spinocerebellar ataxia type 1 (SCA1) is a neurodegenerative disorder caused by expansion of a CAG trinucleotide repeat. In this study, we describe the identification and characterization of the gene harbouring this repeat. The SCA1 transcript is 10,660 bases and is transcribed from both the wild type and SCA1 alleles. The CAG repeat, coding for a polyglutamine tract, lies within the coding region. The gene spans 450 kb of genomic DNA and is organized in nine exons. The first seven fall in the 5' untranslated region and the last two contain the coding region, and a 7,277 basepairs 3' untranslated region. The first four non-coding exons undergo alternative splicing in several tissues. These features suggest that the transcriptional and translational regulation of ataxin-1, the SCA1 encoded protein, may be complex.

A high resolution deletion map of human chromosome Xp22.

We have developed a 32-interval deletion panel for human chromosome Xp22 spanning about 30 megabases of genomic DNA. DNA samples from 50 patients with chromosomal rearrangements involving Xp22 were tested with 60 markers using a polymerase chain reaction strategy. The ensuing deletion map allowed us to confirm and refine the order of previously isolated and newly developed markers. Our mapping panel will provide the framework for mapping new sequences, for orienting chromosome walks in the region and for projects aimed at isolating genes responsible for diseases mapping to Xp22.

Early urodynamic findings after complete primary repair of exstrophy.

INTRODUCTION: Detrusor contraction in bladder exstrophy (BE) patients following reconstruction is poorly understood as there are few published studies assessing urodynamic findings in this population. Understanding the ability of the detrusor to contract in BE patients early after closure may be able to inform the longer-term management and potential for the development of future continence in this population. OBJECTIVE: We sought to evaluate early detrusor contraction using urodynamic studies (UDS) in children who had previously undergone complete primary repair of bladder exstrophy (CPRE). We hypothesized that a majority of children with BE would display the presence of normal detrusor contractile function after CPRE. STUDY DESIGN: A retrospective review of our prospectively collected database was performed for all patients with a diagnosis of classic BE who underwent primary CPRE between 2013 and 2017. From this cohort we identified patients with at least one post-operative UDS at 3 years of age or older who had undergone an initial CPRE. Our primary outcome was the presence of a detrusor contraction demonstrated on UDS. RESULTS: There were 50 children (31 male, 19 female) with CBE who underwent CPRE between 2013 and 2017.There were 26 (13 male, 13 female) who met inclusion criteria. Median age was 3.5 (IQR: 3.2-4.7) years at the time of UDS Sixteen of the 26 (61.5%) generated a sustained detrusor contraction generating a void, with a median peak voiding pressure of 38 cm H20 (IQR: 28-51). The median bladder capacity reached was 48 ml, which represented a median of 30% of expected bladder capacity. The median post void residual (PVR) for the entire cohort was 26 ml (IQR: 9, 47) or 51% (IQR: 20%-98%) of their actual bladder capacity, while the median PVR for those children with a sustained detrusor contraction was 18 ml (IQR: 5, 46) or 33% (IQR: 27%, 98%) of their actual bladder capacity. Intraoperative bladder width and bladder dome to bladder neck length did not correlate with the presence of voiding via a detrusor contraction (p = 0.64). DISCUSSION: We present the first study assessing early UDS finding of detrusor contraction in BE patients after CPRE. In our cohort, 61.5% of patients were able to generate a sustained detrusor contraction on UDS which is a higher percentage than has been reported in previous series. A difference in initial surgical management may account for these findings. CONCLUSION: At short term follow up, the majority of children in our cohort were able to produce sustained detrusor contractions sufficient to generate a void per urethra with a modest post void residual volume. Long-term follow-up and repeated UDS will be needed to track detrusor contractility rates, bladder capacities, compliance, post void residuals and ultimately continence rates over time.

Description of two cultivated and two uncultivated new Salinibacter species, one named following the rules of the bacteriological code: Salinibacter grassmerensis sp. nov.; and three named following the rules of the SeqCode: Salinibacter pepae sp. nov., Salinibacter abyssi sp. nov., and Salinibacter pampae sp. nov.

Current -omics methods allow the collection of a large amount of information that helps in describing the microbial diversity in nature. Here, and as a result of a culturomic approach that rendered the collection of thousands of isolates from 5 different hypersaline sites (in Spain, USA and New Zealand), we obtained 21 strains that represent two new Salinibacter species. For these species we propose the names Salinibacter pepae sp. nov. and Salinibacter grassmerensis sp. nov. (showing average nucleotide identity (ANI) values < 95.09% and 87.08% with Sal. ruber M31T, respectively). Metabolomics revealed species-specific discriminative profiles. Sal. ruber strains were distinguished by a higher percentage of polyunsaturated fatty acids and specific N-functionalized fatty acids; and Sal. altiplanensis was distinguished by an increased number of glycosylated molecules. Based on sequence characteristics and inferred phenotype of metagenome-assembled genomes (MAGs), we describe two new members of the genus Salinibacter. These species dominated in different sites and always coexisted with Sal. ruber and Sal. pepae. Based on the MAGs from three Argentinian lakes in the Pampa region of Argentina and the MAG of the Romanian lake Fara Fund, we describe the species Salinibacter pampae sp. nov. and Salinibacter abyssi sp. nov. respectively (showing ANI values 90.94% and 91.48% with Sal. ruber M31T, respectively). Sal. grassmerensis sp. nov. name was formed according to the rules of the International Code for Nomenclature of Prokaryotes (ICNP), and Sal. pepae, Sal. pampae sp. nov. and Sal. abyssi sp. nov. are proposed following the rules of the newly published Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). This work constitutes an example on how classification under ICNP and SeqCode can coexist, and how the official naming a cultivated organism for which the deposit in public repositories is difficult finds an intermediate solution.

Fecal microbiota transplantation promotes reduction of antimicrobial resistance by strain replacement.

Multidrug-resistant organism (MDRO) colonization is a fundamental challenge in antimicrobial resistance. Limited studies have shown that fecal microbiota transplantation (FMT) can reduce MDRO colonization, but its mechanisms are poorly understood. We conducted a randomized, controlled trial of FMT for MDRO decolonization in renal transplant recipients called PREMIX (NCT02922816). Eleven participants were enrolled and randomized 1:1 to FMT or an observation period followed by delayed FMT if stool cultures were MDRO positive at day 36. Participants who were MDRO positive after one FMT were treated with a second FMT. At last visit, eight of nine patients who completed all treatments were MDRO culture negative. FMT-treated participants had longer time to recurrent MDRO infection versus PREMIX-eligible controls who were not treated with FMT. Key taxa (Akkermansia muciniphila, Alistipes putredinis, Phocaeicola dorei, Phascolarctobacterium faecium, Alistipes species, Mesosutterella massiliensis, Barnesiella intestinihominis, and Faecalibacterium prausnitzii) from the single feces donor used in the study that engrafted in recipients and metabolites such as short-chain fatty acids and bile acids in FMT-responding participants uncovered leads for rational microbiome therapeutic and diagnostic development. Metagenomic analyses revealed a previously unobserved mechanism of MDRO eradication by conspecific strain competition in an FMT-treated subset. Susceptible Enterobacterales strains that replaced baseline extended-spectrum ß-lactamase-producing strains were not detectable in donor microbiota manufactured as FMT doses but in one case were detectable in the recipient before FMT. These data suggest that FMT may provide a path to exploit strain competition to reduce MDRO colonization.

Bone mineral density in patients with stroke: relationship with motor impairment and functional mobility.

PURPOSE: Patients with stroke have a 2- to 4-fold increased risk of hip fracture compared to the general population, because of decreased bone mineral density (BMD) on the paretic side and the high incidence of accidental falls. However, the relationship between BMD and stroke-related motor impairment and functional mobility is not known. The purpose of this study was to investigate these relationships. METHOD: A convenience sample of 87 patients with stroke was recruited from an outpatient rehabilitation clinic. Demographics and clinical history were collected, and patients answered questionnaires regarding functional status. Motor impairment was assessed using motor items of the National Institutes of Health Stroke Scale (NIHSS), and BMD was measured using dual energy X-ray absorptiometry. RESULTS: Mean BMD measured at the total hip was lower on the paretic side relative to the contralateral side (0.883 ± 0.148 g/cm2 vs 0.923 ± 0.136 g/cm2;P < .001). Compared to patients without limitations in walking, those reporting limitations had lower BMD at the paretic total hip (0.808 ± 0.141 g/cm2 vs 0.917 ± 0.139 g/cm2;P = .001) and lower BMDZ scores (-0.282 ± 0.888 vs -0.028 ± 0.813;P = .035). A significant correlation was found between mean BMD and the BMDZ score at the total hip on the paretic side and motor impairment in that lower extremity (r = -0.326,P = .003;r = -0.312,P = .004, respectively). CONCLUSION: In patients with stroke, BMD at the paretic hip correlated with motor impairment. Furthermore, ability to ambulate was shown to be a simple yet useful test to determine which individuals had increased bone loss at the paretic versus nonparetichip.

The Dark Side:Stigma in purpose-built senior environments.

This paper focuses on stigma in collective living environments for older adults, specifically multi-level campuses. We contrast two design profiles, a purpose-built campus which opened in 1997, and an older setting that grew by accretion over decades. Purpose-built housing is used as originally intended, designed, and constructed; housing built by accretion has been modified over time to meet changing needs and uses. The separation by care levels in both sites is reflected in their cultures as residents and staff relate to physical levels of care through a vocabulary of fear and stigma. Residents of the independent living building on the purpose-built campus refer to the assisted living building, accessed only through a second floor link, as "the other side" or "the dark side." In this setting we observe stigma assigned to a place in the built environment. By contrast, in the older setting built by accretion over the last century, levels of care feature a less-structured clustering of independent living and assisted living, and common areas were shared by residents from multiple care levels. We have observed less stigma associated with levels of care in this older building. Grounding our analysis in data drawn from ongoing ethnographic research, we focus on the built environment as it relates to stigma in the social environment. The paper concludes with a discussion of the importance and pervasiveness of stigma in senior environments.

Solar salterns as model systems to study the units of bacterial diversity that matter for ecosystem functioning.

Microbial communities often harbor overwhelming species and gene diversity, making it challenging to determine the important units to study this diversity. We argue that the reduced, and thus tractable, microbial diversity of manmade salterns provides an ideal system to advance this cornerstone issue. We review recent time-series genomic and metagenomic studies of the saltern-dominating bacterial and archaeal taxa to show that these taxa form persistent, sequence-discrete, species-like populations. While these populations harbor extensive intra-population gene diversity, even within a single saltern site, only a small minority of these genes appear to be functionally important during environmental perturbations. We outline an approach to detect and track such populations and their ecologically important genes that should be broadly applicable.

Toward quantifying the adaptive role of bacterial pangenomes during environmental perturbations.

Metagenomic surveys have revealed that natural microbial communities are predominantly composed of sequence-discrete, species-like populations but the genetic and/or ecological processes that maintain such populations remain speculative, limiting our understanding of population speciation and adaptation to perturbations. To address this knowledge gap, we sequenced 112 Salinibacter ruber isolates and 12 companion metagenomes from four adjacent saltern ponds in Mallorca, Spain that were experimentally manipulated to dramatically alter salinity and light intensity, the two major drivers of this ecosystem. Our analyses showed that the pangenome of the local Sal. ruber population is open and similar in size (~15,000 genes) to that of randomly sampled Escherichia coli genomes. While most of the accessory (noncore) genes were isolate-specific and showed low in situ abundances based on the metagenomes compared to the core genes, indicating that they were functionally unimportant and/or transient, 3.5% of them became abundant when salinity (but not light) conditions changed and encoded for functions related to osmoregulation. Nonetheless, the ecological advantage of these genes, while significant, was apparently not strong enough to purge diversity within the population. Collectively, our results provide an explanation for how this immense intrapopulation gene diversity is maintained, which has implications for the prokaryotic species concept.

Toward shotgun metagenomic approaches for microbial source tracking sewage spills based on laboratory mesocosms.

Little is known about the genomic diversity of the microbial communities associated with raw municipal wastewater (sewage), including whether microbial populations specific to sewage exist and how such populations could be used to improve source attribution and apportioning in contaminated waters. Herein, we used the influent of three wastewater treatment plants in Atlanta, Georgia (USA) to perturb laboratory freshwater mesocosms, simulating sewage contamination events, and followed these mesocosms with shotgun metagenomics over a 7-day observational period. We describe 15 abundant non-redundant bacterial metagenome-assembled genomes (MAGs) ubiquitous within all sewage inocula yet absent from the unperturbed freshwater control at our analytical limit of detection. Tracking the dynamics of the populations represented by these MAGs revealed varied decay kinetics, depending on (inferred) phenotypes, e.g., anaerobes decayed faster than aerobes under the well-aerated incubation conditions. Notably, a portion of these populations showed decay patterns similar to those of common markers, Enterococcus and HF183. Despite the apparent decay of these populations, the abundance of ß-lactamase encoding genes remained high throughout incubation relative to the control. Lastly, we constructed genomic libraries representing several different fecal sources and outline a bioinformatic approach which leverages these libraries for identifying and apportioning contamination signal among multiple probable sources using shotgun metagenomic data.

ROCker Models for Reliable Detection and Typing of Short-Read Sequences Carrying ß-Lactamase Genes.

Identification of genes encoding ß-lactamases (BLs) from short-read sequences remains challenging due to the high frequency of shared amino acid functional domains and motifs in proteins encoded by BL genes and related non-BL gene sequences. Divergent BL homologs can be frequently missed during similarity searches, which has important practical consequences for monitoring antibiotic resistance. To address this limitation, we built ROCker models that targeted broad classes (e.g., class A, B, C, and D) and individual families (e.g., TEM) of BLs and challenged them with mock 150-bp- and 250-bp-read data sets of known composition. ROCker identifies most-discriminant bit score thresholds in sliding windows along the sequence of the target protein sequence and hence can account for nondiscriminative domains shared by unrelated proteins. BL ROCker models showed a 0% false-positive rate (FPR), a 0% to 4% false-negative rate (FNR), and an up-to-50-fold-higher F1 score [2 × precision × recall/(precision + recall)] compared to alternative methods, such as similarity searches using BLASTx with various e-value thresholds and BL hidden Markov models, or tools like DeepARG, ShortBRED, and AMRFinder. The ROCker models and the underlying protein sequence reference data sets and phylogenetic trees for read placement are freely available through http://enve-omics.ce.gatech.edu/data/rocker-bla. Application of these BL ROCker models to metagenomics, metatranscriptomics, and high-throughput PCR gene amplicon data should facilitate the reliable detection and quantification of BL variants encoded by environmental or clinical isolates and microbiomes and more accurate assessment of the associated public health risk, compared to the current practice. IMPORTANCE Resistance genes encoding ß-lactamases (BLs) confer resistance to the widely prescribed antibiotic class ß-lactams. Therefore, it is important to assess the prevalence of BL genes in clinical or environmental samples for monitoring the spreading of these genes into pathogens and estimating public health risk. However, detecting BLs in short-read sequence data is technically challenging. Our ROCker model-based bioinformatics approach showcases the reliable detection and typing of BLs in complex data sets and thus contributes toward solving an important problem in antibiotic resistance surveillance. The ROCker models developed substantially expand the toolbox for monitoring antibiotic resistance in clinical or environmental settings.