Search for: All records

ABSTRACT Targeted amplicon sequencing is widely used in microbial ecology studies. However, sequencing artifacts and amplification biases are of great concern. To identify sources of these artifacts, a systematic analysis was performed using mock communities comprised of 16S rRNA genes from 33 bacterial strains. Our results indicated that while sequencing errors were generally isolated to low-abundance operational taxonomic units, chimeric sequences were a major source of artifacts. Singleton and doubleton sequences were primarily chimeras. Formation of chimeric sequences was significantly correlated with the GC content of the targeted sequences. Low-GC-content mock community members exhibited lower rates of chimeric sequence formation. GC content also had a large impact on sequence recovery. The quantitative capacity was notably limited, with substantial recovery variations and weak correlation between anticipated and observed strain abundances. The mock community strains with higher GC content had higher recovery rates than strains with lower GC content. Amplification bias was also observed due to the differences in primer affinity. A two-step PCR strategy reduced the number of chimeric sequences by half. In addition, comparative analyses based on the mock communities showed that several widely used sequence processing pipelines/methods, including DADA2, Deblur, UCLUST, UNOISE, and UPARSE, had different advantages and disadvantages in artifact removal and rare species detection. These results are important for improving sequencing quality and reliability and developing new algorithms to process targeted amplicon sequences. IMPORTANCEAmplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches. 

Abstract Complete ammonia oxidizing bacteria coexist with canonical ammonia and nitrite oxidizing bacteria in a wide range of environments. Whether this is due to competitive or cooperative interactions, or a result of niche separation is not yet clear. Understanding the factors driving coexistence of nitrifiers is critical to manage nitrification processes occurring in engineered and natural ecosystems. In this study, microcosm-based experiments were used to investigate the impact of nitrogen source and loading on the population dynamics of nitrifiers in drinking water biofilter media. Shotgun sequencing of DNA followed by co-assembly and reconstruction of metagenome assembled genomes revealed clade A2 comammox bacteria were likely the primary nitrifiers within microcosms and increased in abundance over Nitrosomonas-like ammonia and Nitrospira-like nitrite oxidizing bacteria irrespective of nitrogen source type or loading. Changes in comammox bacterial abundance did not correlate with either ammonia or nitrite oxidizing bacterial abundance in urea-amended systems, where metabolic reconstruction indicated potential for cross-feeding between strict ammonia and nitrite oxidizers. In contrast, comammox bacterial abundance demonstrated a negative correlation with nitrite oxidizers in ammonia-amended systems. This suggests potentially weaker synergistic relationships between strict ammonia and nitrite oxidizers might enable comammox bacteria to displace strict nitrite oxidizers from complex nitrifying communities. 

ABSTRACT Reconstructing microbial genomes from metagenomic short-read data can be challenging due to the unknown and uneven complexity of microbial communities. This complexity encompasses highly diverse populations, which often includes strain variants. Reconstructing high-quality genomes is a crucial part of the metagenomic workflow, as subsequent ecological and metabolic inferences depend on their accuracy, quality, and completeness. In contrast to microbial communities in other ecosystems, there has been no systematic assessment of genome-centric metagenomic workflows for drinking water microbiomes. In this study, we assessed the performance of a combination of assembly and binning strategies for time series drinking water metagenomes that were collected over 6 months. The goal of this study was to identify the combination of assembly and binning approaches that result in high-quality and -quantity metagenome-assembled genomes (MAGs), representing most of the sequenced metagenome. Our findings suggest that the metaSPAdes coassembly strategies had the best performance, as they resulted in larger and less fragmented assemblies, with at least 85% of the sequence data mapping to contigs greater than 1 kbp. Furthermore, a combination of metaSPAdes coassembly strategies and MetaBAT2 produced the highest number of medium-quality MAGs while capturing at least 70% of the metagenomes based on read recruitment. Utilizing different assembly/binning approaches also assists in the reconstruction of unique MAGs from closely related species that would have otherwise collapsed into a single MAG using a single workflow. Overall, our study suggests that leveraging multiple binning approaches with different metaSPAdes coassembly strategies may be required to maximize the recovery of good-quality MAGs. IMPORTANCE Drinking water contains phylogenetic diverse groups of bacteria, archaea, and eukarya that affect the esthetic quality of water, water infrastructure, and public health. Taxonomic, metabolic, and ecological inferences of the drinking water microbiome depend on the accuracy, quality, and completeness of genomes that are reconstructed through the application of genome-resolved metagenomics. Using time series metagenomic data, we present reproducible genome-centric metagenomic workflows that result in high-quality and -quantity genomes, which more accurately signifies the sequenced drinking water microbiome. These genome-centric metagenomic workflows will allow for improved taxonomic and functional potential analysis that offers enhanced insights into the stability and dynamics of drinking water microbial communities. 

ABSTRACT Ammonia availability due to chloramination can promote the growth of nitrifying organisms, which can deplete chloramine residuals and result in operational problems for drinking water utilities. In this study, we used a metagenomic approach to determine the identity and functional potential of microorganisms involved in nitrogen biotransformation within chloraminated drinking water reservoirs. Spatial changes in the nitrogen species included an increase in nitrate concentrations accompanied by a decrease in ammonium concentrations with increasing distance from the site of chloramination. This nitrifying activity was likely driven by canonical ammonia-oxidizing bacteria (i.e., Nitrosomonas ) and nitrite-oxidizing bacteria (i.e., Nitrospira ) as well as by complete-ammonia-oxidizing (i.e., comammox) Nitrospira -like bacteria. Functional annotation was used to evaluate genes associated with nitrogen metabolism, and the community gene catalogue contained mostly genes involved in nitrification, nitrate and nitrite reduction, and nitric oxide reduction. Furthermore, we assembled 47 high-quality metagenome-assembled genomes (MAGs) representing a highly diverse assemblage of bacteria. Of these, five MAGs showed high coverage across all samples, which included two Nitrosomonas, Nitrospira, Sphingomonas , and Rhizobiales -like MAGs. Systematic genome-level analyses of these MAGs in relation to nitrogen metabolism suggest that under ammonia-limited conditions, nitrate may be also reduced back to ammonia for assimilation. Alternatively, nitrate may be reduced to nitric oxide and may potentially play a role in regulating biofilm formation. Overall, this study provides insight into the microbial communities and their nitrogen metabolism and, together with the water chemistry data, improves our understanding of nitrogen biotransformation in chloraminated drinking water distribution systems. IMPORTANCE Chloramines are often used as a secondary disinfectant when free chlorine residuals are difficult to maintain. However, chloramination is often associated with the undesirable effect of nitrification, which results in operational problems for many drinking water utilities. The introduction of ammonia during chloramination provides a potential source of nitrogen either through the addition of excess ammonia or through chloramine decay. This promotes the growth of nitrifying microorganisms and provides a nitrogen source (i.e., nitrate) for the growth for other organisms. While the roles of canonical ammonia-oxidizing and nitrite-oxidizing bacteria in chloraminated drinking water systems have been extensively investigated, those studies have largely adopted a targeted gene-centered approach. Further, little is known about the potential long-term cooccurrence of complete-ammonia-oxidizing (i.e., comammox) bacteria and the potential metabolic synergies of nitrifying organisms with their heterotrophic counterparts that are capable of denitrification and nitrogen assimilation. This study leveraged data obtained for genome-resolved metagenomics over a time series to show that while nitrifying bacteria are dominant and likely to play a major role in nitrification, their cooccurrence with heterotrophic organisms suggests that nitric oxide production and nitrate reduction to ammonia may also occur in chloraminated drinking water systems. 

Microcystins produced during harmful cyanobacterial blooms are a public health concern. Although patterns are emerging, the environmental cues that stimulate production of microcystin remain confusing, hindering our ability to predict fluctuations in bloom toxicity. In earlier work, growth at cool temperatures relative to optimum (18°C vs. 26°C) was confirmed to increase microcystin quota in batch cultures of Microcystis aeruginosa NIES-843. Here, we tested this response in M. aeruginosa PCC 7806 using continuous cultures to examine temporal dynamics and using RNA-sequencing to investigate the physiological nature of the response. A temperature reduction from 26 to 19°C increased microcystin quota ∼2-fold, from an average of ∼464 ag μm –3 cell volume to ∼891 ag μm –3 over a 7–9 d period. Reverting the temperature to 26°C returned the cellular microcystin quota to ∼489 ag μm –3 . Long periods (31–42 d) at 19°C did not increase or decrease microcystin quota beyond that observed at 7–9 d. Nitrogen concentration had little effect on the overall response. RNA sequencing indicated that the decrease in temperature to 19°C induced a classic cold-stress response in M. aeruginosa PCC 7806, but this operated on a different timescale than the increased microcystin production. Microcystin quota showed a strong 48- to 72-h time-lag correlation to mcy gene expression, but no correlation to concurrent mcy expression. This work confirms an effect of temperature on microcystin quota and extends our understanding of the physiological nature of the response.