Metagenomic and Resistome Analysis of a Full-Scale Municipal Wastewater Treatment Plant in Singapore Containing Membrane Bioreactors.

Reclaimed water provides a water supply alternative to address problems of scarcity in urbanized cities with high living densities and limited natural water resources. In this study, wastewater metagenomes from 6 stages of a wastewater treatment plant (WWTP) integrating conventional and membrane bioreactor (MBR) treatment were evaluated for diversity of antibiotic resistance genes (ARGs) and bacteria, and relative abundance of class 1 integron integrases (intl1). ARGs confering resistance to 12 classes of antibiotics (ARG types) persisted through the treatment stages, which included genes that confer resistance to aminoglycoside [aadA, aph(6)-I, aph(3')-I, aac(6')-I, aac(6')-II, ant(2″)-I], beta-lactams [class A, class C, class D beta-lactamases (bla OXA)], chloramphenicol (acetyltransferase, exporters, floR, cmIA), fosmidomycin (rosAB), macrolide-lincosamide-streptogramin (macAB, ereA, ermFB), multidrug resistance (subunits of transporters), polymyxin (arnA), quinolone (qnrS), rifamycin (arr), sulfonamide (sul1, sul2), and tetracycline (tetM, tetG, tetE, tet36, tet39, tetR, tet43, tetQ, tetX). Although the ARG subtypes in sludge and MBR effluents reduced in diversity relative to the influent, clinically relevant beta lactamases (i.e., bla KPC, bla OXA) were detected, casting light on other potential point sources of ARG dissemination within the wastewater treatment process. To gain a deeper insight into the types of bacteria that may survive the MBR removal process, genome bins were recovered from metagenomic data of MBR effluents. A total of 101 close to complete draft genomes were assembled and annotated to reveal a variety of bacteria bearing metal resistance genes and ARGs in the MBR effluent. Three bins in particular were affiliated to Mycobacterium smegmatis, Acinetobacter Iwoffii, and Flavobacterium psychrophila, and carried aquired ARGs aac(2')-Ib, bla OXA-278, and tet36 respectively. In terms of indicator organisms, cumulative log removal values (LRV) of Escherichia coli, Enterococci, and P. aeruginosa from influent to conventional treated effluent was lower (0-2.4), compared to MBR effluent (5.3-7.4). We conclude that MBR is an effective treatment method for reducing fecal indicators and ARGs; however, incomplete removal of P. aeruginosa in MBR treated effluents (<8 MPN/100 mL) and the presence of ARGs and intl1 underscores the need to establish if further treatment should be applied prior to reuse.

Draft Genome Sequences of Four Multidrug-Resistant Pseudomonas aeruginosa Isolates from Hospital Wastewater in Singapore.

Four multidrug-resistant Pseudomonas aeruginosa isolates were cultured from intensive care unit wastewater. All isolates exhibited resistance to carbapenem and extended-spectrum beta-lactam antibiotics. Genome characterization revealed the presence of beta-lactamase resistance genes (bla PAO and bla OXA), and three out of the four isolates carried the bla NDM-1 gene encoding resistance against carbapenems.

Corrigendum: Characterization of Metagenomes in Urban Aquatic Compartments Reveals High Prevalence of Clinically Relevant Antibiotic Resistance Genes in Wastewaters.

[This corrects the article on p. 2200 in vol. 8, PMID: 29201017.].

Characterization of Metagenomes in Urban Aquatic Compartments Reveals High Prevalence of Clinically Relevant Antibiotic Resistance Genes in Wastewaters.

The dissemination of antimicrobial resistance (AMR) is an escalating problem and a threat to public health. Comparative metagenomics was used to investigate the occurrence of antibiotic resistant genes (ARGs) in wastewater and urban surface water environments in Singapore. Hospital and municipal wastewater (n = 6) were found to have higher diversity and average abundance of ARGs (303 ARG subtypes, 197,816 x/Gb) compared to treated wastewater effluent (n = 2, 58 ARG subtypes, 2,692 x/Gb) and surface water (n = 5, 35 subtypes, 7,985 x/Gb). A cluster analysis showed that the taxonomic composition of wastewaters was highly similar and had a bacterial community composition enriched in gut bacteria (Bacteroides, Faecalibacterium, Bifidobacterium, Blautia, Roseburia, Ruminococcus), the Enterobacteriaceae group (Klebsiella, Aeromonas, Enterobacter) and opportunistic pathogens (Prevotella, Comamonas, Neisseria). Wastewater, treated effluents and surface waters had a shared resistome of 21 ARGs encoding multidrug resistant efflux pumps or resistance to aminoglycoside, macrolide-lincosamide-streptogramins (MLS), quinolones, sulfonamide, and tetracycline resistance which suggests that these genes are wide spread across different environments. Wastewater had a distinctively higher average abundance of clinically relevant, class A beta-lactamase resistant genes (i.e., blaKPC, blaCTX-M, blaSHV, blaTEM). The wastewaters from clinical isolation wards, in particular, had a exceedingly high levels of blaKPC-2 genes (142,200 x/Gb), encoding for carbapenem resistance. Assembled scaffolds (16 and 30 kbp) from isolation ward wastewater samples indicated this gene was located on a Tn3-based transposon (Tn4401), a mobilization element found in Klebsiella pneumonia plasmids. In the longer scaffold, transposable elements were flanked by a toxin-antitoxin (TA) system and other metal resistant genes that likely increase the persistence, fitness and propagation of the plasmid in the bacterial host under conditions of stress. A few bacterial species (Enterobacter cloacae, Klebsiella pneumoniae, Citrobacter freundii, Pseudomonas aeruginosa) that were cultured from the isolation ward wastewaters on CHROMagar media harbored the blaKPC-2 gene. This suggests that hospital wastewaters derived from clinical specialty wards are hotspots for the spread of AMR. Assembled scaffolds of other mobile genetic elements such as IncQ and IncF plasmids bearing quinolone resistance genes (qnrS1, qnrS2) and the class A beta-lactamase gene (blaTEM-1) were recovered in wastewater samples which may aid the transfer of AMR.

Vitamin and Amino Acid Auxotrophy in Anaerobic Consortia Operating under Methanogenic Conditions.

Syntrophy among Archaea and Bacteria facilitates the anaerobic degradation of organic compounds to CH4 and CO2. Particularly during aliphatic and aromatic hydrocarbon mineralization, as in the case of crude oil reservoirs and petroleum-contaminated sediments, metabolic interactions between obligate mutualistic microbial partners are of central importance. Using micromanipulation combined with shotgun metagenomic approaches, we describe the genomes of complex consortia within short-chain alkane-degrading cultures operating under methanogenic conditions. Metabolic reconstruction revealed that only a small fraction of genes in the metagenome-assembled genomes encode the capacity for fermentation of alkanes facilitated by energy conservation linked to H2 metabolism. Instead, the presence of inferred lifestyles based on scavenging anabolic products and intermediate fermentation products derived from detrital biomass was a common feature. Additionally, inferred auxotrophy for vitamins and amino acids suggests that the hydrocarbon-degrading microbial assemblages are structured and maintained by multiple interactions beyond the canonical H2-producing and syntrophic alkane degrader-methanogen partnership. Compared to previous work, our report points to a higher order of complexity in microbial consortia engaged in anaerobic hydrocarbon transformation. IMPORTANCE Microbial interactions between Archaea and Bacteria mediate many important chemical transformations in the biosphere from degrading abundant polymers to synthesis of toxic compounds. Two of the most pressing issues in microbial interactions are how consortia are established and how we can modulate these microbial communities to express desirable functions. Here, we propose that public goods (i.e., metabolites of high energy demand in biosynthesis) facilitate energy conservation for life under energy-limited conditions and determine the assembly and function of the consortia. Our report suggests that an understanding of public good dynamics could result in new ways to improve microbial pollutant degradation in anaerobic systems.

Microbial potential for carbon and nutrient cycling in a geogenic supercritical carbon dioxide reservoir.

Microorganisms catalyze carbon cycling and biogeochemical reactions in the deep subsurface and thus may be expected to influence the fate of injected supercritical (sc) CO2 following geological carbon sequestration (GCS). We hypothesized that natural subsurface scCO2 reservoirs, which serve as analogs for the long-term fate of sequestered scCO2 , harbor a 'deep carbonated biosphere' with carbon cycling potential. We sampled subsurface fluids from scCO2 -water separators at a natural scCO2 reservoir at McElmo Dome, Colorado for analysis of 16S rRNA gene diversity and metagenome content. Sequence annotations indicated dominance of Sulfurospirillum, Rhizobium, Desulfovibrio and four members of the Clostridiales family. Genomes extracted from metagenomes using homology and compositional approaches revealed diverse mechanisms for growth and nutrient cycling, including pathways for CO2 and N2 fixation, anaerobic respiration, sulfur oxidation, fermentation and potential for metabolic syntrophy. Differences in biogeochemical potential between two production well communities were consistent with differences in fluid chemical profiles, suggesting a potential link between microbial activity and geochemistry. The existence of a microbial ecosystem associated with the McElmo Dome scCO2 reservoir indicates that potential impacts of the deep biosphere on CO2 fate and transport should be taken into consideration as a component of GCS planning and modelling.

Next-Generation Sequencing Assessment of Eukaryotic Diversity in Oil Sands Tailings Ponds Sediments and Surface Water.

Tailings ponds in the Athabasca oil sands (Canada) contain fluid wastes, generated by the extraction of bitumen from oil sands ores. Although the autochthonous prokaryotic communities have been relatively well characterized, almost nothing is known about microbial eukaryotes living in the anoxic soft sediments of tailings ponds or in the thin oxic layer of water that covers them. We carried out the first next-generation sequencing study of microbial eukaryotic diversity in oil sands tailings ponds. In metagenomes prepared from tailings sediment and surface water, we detected very low numbers of sequences encoding eukaryotic small subunit ribosomal RNA representing seven major taxonomic groups of protists. We also produced and analysed three amplicon-based 18S rRNA libraries prepared from sediment samples. These revealed a more diverse set of taxa, 169 different OTUs encompassing up to eleven higher order groups of eukaryotes, according to detailed classification using homology searching and phylogenetic methods. The 10 most abundant OTUs accounted for > 90% of the total of reads, vs. large numbers of rare OTUs (< 1% abundance). Despite the anoxic and hydrocarbon-enriched nature of the environment, the tailings ponds harbour complex communities of microbial eukaryotes indicating that these organisms should be taken into account when studying the microbiology of the oil sands.

Next-generation sequencing (NGS) for assessment of microbial water quality: current progress, challenges, and future opportunities.

Water quality is an emergent property of a complex system comprised of interacting microbial populations and introduced microbial and chemical contaminants. Studies leveraging next-generation sequencing (NGS) technologies are providing new insights into the ecology of microbially mediated processes that influence fresh water quality such as algal blooms, contaminant biodegradation, and pathogen dissemination. In addition, sequencing methods targeting small subunit (SSU) rRNA hypervariable regions have allowed identification of signature microbial species that serve as bioindicators for sewage contamination in these environments. Beyond amplicon sequencing, metagenomic and metatranscriptomic analyses of microbial communities in fresh water environments reveal the genetic capabilities and interplay of waterborne microorganisms, shedding light on the mechanisms for production and biodegradation of toxins and other contaminants. This review discusses the challenges and benefits of applying NGS-based methods to water quality research and assessment. We will consider the suitability and biases inherent in the application of NGS as a screening tool for assessment of biological risks and discuss the potential and limitations for direct quantitative interpretation of NGS data. Secondly, we will examine case studies from recent literature where NGS based methods have been applied to topics in water quality assessment, including development of bioindicators for sewage pollution and microbial source tracking, characterizing the distribution of toxin and antibiotic resistance genes in water samples, and investigating mechanisms of biodegradation of harmful pollutants that threaten water quality. Finally, we provide a short review of emerging NGS platforms and their potential applications to the next generation of water quality assessment tools.

Anaerobic alkane biodegradation by cultures enriched from oil sands tailings ponds involves multiple species capable of fumarate addition.

A methanogenic short-chain alkane-degrading culture (SCADC) was enriched from oil sands tailings and transferred several times with a mixture of C6, C7, C8 and C10 n-alkanes as the predominant organic carbon source, plus 2-methylpentane, 3-methylpentane and methylcyclopentane as minor components. Cultures produced ∼40% of the maximum theoretical methane during 18 months incubation while depleting the n-alkanes, 2-methylpentane and methylcyclopentane. Substrate depletion correlated with detection of metabolites characteristic of fumarate activation of 2-methylpentane and methylcyclopentane, but not n-alkane metabolites. During active methanogenesis with the mixed alkanes, reverse-transcription PCR confirmed the expression of functional genes (assA and bssA) associated with hydrocarbon addition to fumarate. Pyrosequencing of 16S rRNA genes amplified during active alkane degradation revealed enrichment of Clostridia (particularly Peptococcaceae) and methanogenic Archaea (Methanosaetaceae and Methanomicrobiaceae). Methanogenic cultures transferred into medium containing sulphate produced sulphide, depleted n-alkanes and produced the corresponding succinylated alkane metabolites, but were slow to degrade 2-methylpentane and methylcyclopentane; these cultures were enriched in Deltaproteobacteria rather than Clostridia. 3-Methylpentane was not degraded by any cultures. Thus, nominally methanogenic oil sands tailings harbour dynamic and versatile hydrocarbon-degrading fermentative syntrophs and sulphate reducers capable of degrading n-, iso- and cyclo-alkanes by addition to fumarate.

Comparative analysis of metagenomes from three methanogenic hydrocarbon-degrading enrichment cultures with 41 environmental samples.

Methanogenic hydrocarbon metabolism is a key process in subsurface oil reservoirs and hydrocarbon-contaminated environments and thus warrants greater understanding to improve current technologies for fossil fuel extraction and bioremediation. In this study, three hydrocarbon-degrading methanogenic cultures established from two geographically distinct environments and incubated with different hydrocarbon substrates (added as single hydrocarbons or as mixtures) were subjected to metagenomic and 16S rRNA gene pyrosequencing to test whether these differences affect the genetic potential and composition of the communities. Enrichment of different putative hydrocarbon-degrading bacteria in each culture appeared to be substrate dependent, though all cultures contained both acetate- and H2-utilizing methanogens. Despite differing hydrocarbon substrates and inoculum sources, all three cultures harbored genes for hydrocarbon activation by fumarate addition (bssA, assA, nmsA) and carboxylation (abcA, ancA), along with those for associated downstream pathways (bbs, bcr, bam), though the cultures incubated with hydrocarbon mixtures contained a broader diversity of fumarate addition genes. A comparative metagenomic analysis of the three cultures showed that they were functionally redundant despite their enrichment backgrounds, sharing multiple features associated with syntrophic hydrocarbon conversion to methane. In addition, a comparative analysis of the culture metagenomes with those of 41 environmental samples (containing varying proportions of methanogens) showed that the three cultures were functionally most similar to each other but distinct from other environments, including hydrocarbon-impacted environments (for example, oil sands tailings ponds and oil-affected marine sediments). This study provides a basis for understanding key functions and environmental selection in methanogenic hydrocarbon-associated communities.

Draft Genome Sequence of Uncultivated Desulfosporosinus sp. Strain Tol-M, Obtained by Stable Isotope Probing Using [13C6]Toluene.

A draft Desulfosporosinus genome was assembled from the metagenome of a methanogenic [(13)C6]toluene-degrading community. The Desulfosporosinus sp. strain Tol-M genome is distinguished from that of previously published Desulfosporosinus strain by containing bss, bbs, and bam genes encoding enzymes for anaerobic biodegradation of monoaromatic hydrocarbons and lacking dsrAB genes for dissimilatory sulfate reduction.

Draft Genome Sequence of Uncultivated Toluene-Degrading Desulfobulbaceae Bacterium Tol-SR, Obtained by Stable Isotope Probing Using [13C6]Toluene.

The draft genome of a member of the bacterial family Desulfobulbaceae (phylum Deltaproteobacteria) was assembled from the metagenome of a sulfidogenic [(13)C6]toluene-degrading enrichment culture. The "Desulfobulbaceae bacterium Tol-SR" genome is distinguished from related, previously sequenced genomes by suites of genes associated with anaerobic toluene metabolism, including bss, bbs, and bam.

Draft genome sequences of campylobacterales (epsilonproteobacteria) obtained from methanogenic oil sands tailings pond metagenomes.

Draft genome sequences of two Campylobacterales (Sulfurospirillum sp. strain SCADC and Sulfuricurvum sp. strain MLSB [Mildred Lake Settling Basin]) were obtained by taxonomic binning of metagenomes originating from an oil sands tailings pond. Both genomes contain soxABXYZ genes involved in sulfur oxidation, highlighting their potential roles in sulfur cycling in oil sands tailings ponds.

Draft Genome Sequences of Three Smithella spp. Obtained from a Methanogenic Alkane-Degrading Culture and Oil Field Produced Water.

Two draft genomes affiliated with Smithella spp. were obtained from a methanogenic alkane-degrading enrichment culture by single-cell sorting and metagenome contig binning, and a third was obtained by single-cell sorting of oil field produced water. Two genomes contained putative assABC genes encoding alkylsuccinate synthase, indicating genetic potential for fumarate activation of alkanes.

Draft Genome Sequence of Uncultivated Firmicutes (Peptococcaceae SCADC) Single Cells Sorted from Methanogenic Alkane-Degrading Cultures.

The draft genome of an uncultivated bacterium affiliated with the Peptococcaceae was reconstructed by co-assembling Illumina MiSeq sequences from three single cells sorted by microfluidics from two methanogenic alkane-degrading cultures. Peptococcaceae SCADC (short-chain alkane-degrading culture) may be genetically capable of anaerobic alkane activation by fumarate addition in the absence of sulfate.

Metagenomic analysis of an anaerobic alkane-degrading microbial culture: potential hydrocarbon-activating pathways and inferred roles of community members.

A microbial community (short-chain alkane-degrading culture, SCADC) enriched from an oil sands tailings pond was shown to degrade C6-C10 alkanes under methanogenic conditions. Total genomic DNA from SCADC was subjected to 454 pyrosequencing, Illumina paired-end sequencing, and 16S rRNA amplicon pyrotag sequencing; the latter revealed 320 operational taxonomic units at 5% distance. Metagenomic sequences were subjected to in-house quality control and co-assembly, yielding 984 086 contigs, and annotation using MG-Rast and IMG. Substantial nucleotide and protein recruitment to Methanosaeta concilii, Syntrophus aciditrophicus, and Desulfobulbus propionicus reference genomes suggested the presence of closely related strains in SCADC; other genomes were not well mapped, reflecting the paucity of suitable reference sequences for such communities. Nonetheless, we detected numerous homologues of putative hydrocarbon succinate synthase genes (e.g., assA, bssA, and nmsA) implicated in anaerobic hydrocarbon degradation, suggesting the ability of the SCADC microbial community to initiate methanogenic alkane degradation by addition to fumarate. Annotation of a large contig revealed analogues of the ass operon 1 in the alkane-degrading sulphate-reducing bacterium Desulfatibacillum alkenivorans AK-01. Despite being enriched under methanogenic-fermentative conditions, additional metabolic functions inferred by COG profiling indicated multiple CO(2) fixation pathways, organic acid utilization, hydrogenase activity, and sulphate reduction.