Connecting structure to function with the recovery of over 1000 high-quality metagenome-assembled genomes from activated sludge using long-read sequencing

Magnetite-mediated shifts in denitrifying consortia in bioelectrochemical system: Insights into species selection and metabolic dynamics

Conductive materials, such as magnetite, are recognized for their ability to enhance electron transfer and stimulate microbial metabolic activities. This study aimed to elucidate the metabolic potential and species interactions of dominant microbial species within complex communities influenced by magnetite. It indicated that the optimal dosage of magnetite at 4.5 mg/cm², would significantly improve denitrification efficiency and then reduce the time for removing 50 mg/L nitrate by 24.33 %. This enhancement was attributed to the reduced charge transfer resistance and the promoted formation of extracellular polymeric substances (EPS) facilitated by magnetite. Metagenomic analysis revealed that magnetite addition mitigated the competition among truncated denitrifiers for downstream nitrogen species, diminished the contribution of bacteria with complete nitrogen metabolism pathways to denitrification, and fostered a transition towards co-denitrification through interspecies cooperation, consequently leading to decreased nitrite accumulation and increased tolerance to nitrate shock loads. Furthermore, an in-depth study on a key species, Geobacter anodireducens JN93 within the bioelectrochemical system revealed that while magnetite with varying Fe(II) and Fe(III) ratios improved denitrification performance, the metabolic potential of Geobacter sp. varied for different nitrogen metabolism pathways. Collectively, this research provides insights into the microecological effects of magnetite on denitrifying consortia by shifting interspecific interactions via enhanced electron transfer.

Growth of complete ammonia oxidizers on guanidine

Guanidine is a chemically stable nitrogen compound that is excreted in human urine and is widely used in manufacturing of plastics, as a flame retardant and as a component of propellants, and is well known as a protein denaturant in biochemistry1-3. Guanidine occurs widely in nature and is used by several microorganisms as a nitrogen source, but microorganisms growing on guanidine as the only substrate have not yet been identified. Here we show that the complete ammonia oxidizer (comammox) Nitrospira inopinata and probably most other comammox microorganisms can grow on guanidine as the sole source of energy, reductant and nitrogen. Proteomics, enzyme kinetics and the crystal structure of a N. inopinata guanidinase homologue demonstrated that it is a bona fide guanidinase. Incubation experiments with comammox-containing agricultural soil and wastewater treatment plant microbiomes suggested that guanidine serves as substrate for nitrification in the environment. The identification of guanidine as a growth substrate for comammox shows an unexpected niche of these globally important nitrifiers and offers opportunities for their isolation.

Non-redundant metagenome-assembled genomes of activated sludge reactors at different disturbances and scales

Metagenome-assembled genomes (MAGs) are microbial genomes reconstructed from metagenomic data and can be assigned to known taxa or lead to uncovering novel ones. MAGs can provide insights into how microbes interact with the environment. Here, we performed genome-resolved metagenomics on sequencing data from four studies using sequencing batch reactors at microcosm (~25 mL) and mesocosm (~4 L) scales inoculated with sludge from full-scale wastewater treatment plants. These studies investigated how microbial communities in such plants respond to two environmental disturbances: the presence of toxic 3-chloroaniline and changes in organic loading rate. We report 839 non-redundant MAGs with at least 50% completeness and 10% contamination (MIMAG medium-quality criteria). From these, 399 are of putative high-quality, while sixty-seven meet the MIMAG high-quality criteria. MAGs in this catalogue represent the microbial communities in sixty-eight laboratory-scale reactors used for the disturbance experiments, and in the full-scale wastewater treatment plant which provided the source sludge. This dataset can aid meta-studies aimed at understanding the responses of microbial communities to disturbances, particularly as ecosystems confront rapid environmental changes.

From soil to sequence: filling the critical gap in genome-resolved metagenomics is essential to the future of soil microbial ecology

Soil microbiomes are heterogeneous, complex microbial communities. Metagenomic analysis is generating vast amounts of data, creating immense challenges in sequence assembly and analysis. Although advances in technology have resulted in the ability to easily collect large amounts of sequence data, soil samples containing thousands of unique taxa are often poorly characterized. These challenges reduce the usefulness of genome-resolved metagenomic (GRM) analysis seen in other fields of microbiology, such as the creation of high quality metagenomic assembled genomes and the adoption of genome scale modeling approaches. The absence of these resources restricts the scale of future research, limiting hypothesis generation and the predictive modeling of microbial communities. Creating publicly available databases of soil MAGs, similar to databases produced for other microbiomes, has the potential to transform scientific insights about soil microbiomes without requiring the computational resources and domain expertise for assembly and binning.

Advancing environmental biotechnology with microbial community modeling rooted in functional 'omics

Emerging biotechnologies that solve pressing environmental and climate emergencies will require harnessing the vast functional diversity of the underlying microbiomes driving such engineered processes. Modeling is a critical aspect of process engineering that informs system design as well as aids diagnostic optimization of performance. 'Conventional' bioprocess models assume homogenous biomass within functional guilds and thus fail to predict emergent properties of diverse microbial physiologies, such as product specificity and community interactions. Yet, recent advances in functional 'omics-based approaches can provide a 'lens' through which we can probe and measure in situ ecophysiologies of environmental microbiomes. Here, we overview microbial community modeling approaches that incorporate functional 'omics data, which we posit can advance our ability to design and control new environmental biotechnologies going forward.

Metagenome-assembled genomes (MAGs) suggest an acetate-driven protective role in gut microbiota disrupted by Clostridioides difficile

Clostridioides difficile may have a negative impact on gut microbiota composition in terms of diversity and abundance, thereby triggering functional changes supported by the differential presence of genes involved in significant metabolic pathways, such as short-chain fatty acids (SCFA). This work has evaluated shotgun metagenomics data regarding 48 samples from four groups classified according to diarrhea acquisition site (community- and healthcare facility-onset) and positive or negative Clostridioides difficile infection (CDI) result. The metagenomic-assembled genomes (MAGs) obtained from each sample were taxonomically assigned for preliminary comparative analysis concerning differences in composition among groups. The predicted genes involved in metabolism, transport, and signaling remained constant in microbiota members; characteristic patterns were observed in MAGs and genes involved in SCFA butyrate and acetate metabolic pathways for each study group. A decrease in genera and species, as well as relative MAG abundance with the presence of the acetate metabolism-related gene, was evident in the HCFO/- group. Increased antibiotic resistance markers (ARM) were observed in MAGs along with the genes involved in acetate metabolism. The results highlight the need to explore the role of acetate in greater depth as a potential protector of the imbalances produced by CDI, as occurs in other inflammatory intestinal diseases.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs.

IMPORTANCE

Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Distributions, interactions, and dynamics of prokaryotes and phages in a hybrid biological wastewater treatment system

BACKGROUND

Understanding the interactions and dynamics of microbiotas within biological wastewater treatment systems is essential for ensuring their stability and long-term sustainability. In this study, we developed a systematic framework employing multi-omics and Hi-C sequencing to extensively investigate prokaryotic and phage communities within a hybrid biofilm and activated sludge system.

RESULTS

We uncovered distinct distribution patterns, metabolic capabilities, and activities of functional prokaryotes through the analysis of 454 reconstructed prokaryotic genomes. Additionally, we reconstructed a phage catalog comprising 18,645 viral operational taxonomic units (vOTUs) with high length and contiguity using hybrid assembly, and a distinct distribution of phages was depicted between activated sludge (AS) and biofilm. Importantly, 1340 host-phage pairs were established using Hi-C and conventional in silico methods, unveiling the host-determined phage prevalence. The majority of predicted hosts were found to be involved in various crucial metabolic processes, highlighting the potential vital roles of phages in influencing substance metabolism within this system. Moreover, auxiliary metabolic genes (AMGs) related to various categories (e.g., carbohydrate degradation, sulfur metabolism, transporter) were predicted. Subsequent activity analysis emphasized their potential ability to mediate host metabolism during infection. We also profiled the temporal dynamics of phages and their associated hosts using 13-month time-series metagenomic data, further demonstrating their tight interactions. Notably, we observed lineage-specific infection patterns, such as potentially host abundance- or phage/host ratio-driven phage population changes.

CONCLUSIONS

The insights gained from this research contribute to the growing body of knowledge surrounding interactions and dynamics of host-phage and pave the way for further exploration and potential applications in the field of microbial ecology. Video Abstract.

METASEED: a novel approach to full-length 16S rRNA gene reconstruction from short read data

BACKGROUND

With the emergence of Oxford Nanopore technology, now the on-site sequencing of 16S rRNA from environments is available. Due to the error level and structure, the analysis of such data demands some database of reference sequences. However, many taxa from complex and diverse environments, have poor representation in publicly available databases. In this paper, we propose the METASEED pipeline for the reconstruction of full-length 16S sequences from such environments, in order to improve the reference for the subsequent use of on-site sequencing.

RESULTS

We show that combining high-precision short-read sequencing of both 16S and full metagenome from the same samples allow us to reconstruct high-quality 16S sequences from the more abundant taxa. A significant novelty is the carefully designed collection of metagenome reads that matches the 16S amplicons, based on a combination of uniqueness and abundance. Compared to alternative approaches this produces superior results.

CONCLUSION

Our pipeline will facilitate numerous studies associated with various unknown microorganisms, thus allowing the comprehension of the diverse environments. The pipeline is a potential tool in generating a full length 16S rRNA gene database for any environment.

Unveiling genomic features linked to traits of plant growth-promoting bacterial communities from sugarcane

Microorganisms are ubiquitous, and those inhabiting plants have been the subject of several studies. Plant-associated bacteria exhibit various biological mechanisms that enable them to colonize host plants and, in some cases, enhance their fitness. In this study, we describe the genomic features predicted to be associated with plant growth-promoting traits in six bacterial communities isolated from sugarcane. The use of highly accurate single-molecule real-time sequencing technology for metagenomic samples from these bacterial communities allowed us to recover 17 genomes. The taxonomic assignments for the binned genomes were performed, revealing taxa distributed across three main phyla: Bacillota, Bacteroidota, and Pseudomonadota, with the latter being the most representative. Subsequently, we functionally annotated the metagenome-assembled genomes (MAGs) to characterize their metabolic pathways related to plant growth-promoting traits. Our study successfully identified the enrichment of important functions related to phosphate and potassium acquisition, modulation of phytohormones, and mechanisms for coping with abiotic stress. These findings could be linked to the robust colonization of these sugarcane endophytes.

Metatranscriptomics-guided genome-scale metabolic reconstruction reveals the carbon flux and trophic interaction in methanogenic communities

BACKGROUND

Despite rapid advances in genomic-resolved metagenomics and remarkable explosion of metagenome-assembled genomes (MAGs), the function of uncultivated anaerobic lineages and their interactions in carbon mineralization remain largely uncertain, which has profound implications in biotechnology and biogeochemistry.

RESULTS

In this study, we combined long-read sequencing and metatranscriptomics-guided metabolic reconstruction to provide a genome-wide perspective of carbon mineralization flow from polymers to methane in an anaerobic bioreactor. Our results showed that incorporating long reads resulted in a substantial improvement in the quality of metagenomic assemblies, enabling the effective recovery of 132 high-quality genomes meeting stringent criteria of minimum information about a metagenome-assembled genome (MIMAG). In addition, hybrid assembly obtained 51% more prokaryotic genes in comparison to the short-read-only assembly. Metatranscriptomics-guided metabolic reconstruction unveiled the remarkable metabolic flexibility of several novel Bacteroidales-affiliated bacteria and populations from Mesotoga sp. in scavenging amino acids and sugars. In addition to recovering two circular genomes of previously known but fragmented syntrophic bacteria, two newly identified bacteria within Syntrophales were found to be highly engaged in fatty acid oxidation through syntrophic relationships with dominant methanogens Methanoregulaceae bin.74 and Methanothrix sp. bin.206. The activity of bin.206 preferring acetate as substrate exceeded that of bin.74 with increasing loading, reinforcing the substrate determinantal role.

CONCLUSION

Overall, our study uncovered some key active anaerobic lineages and their metabolic functions in this complex anaerobic ecosystem, offering a framework for understanding carbon transformations in anaerobic digestion. These findings advance the understanding of metabolic activities and trophic interactions between anaerobic guilds, providing foundational insights into carbon flux within both engineered and natural ecosystems. Video Abstract.

Effects of different assembly strategies on gene annotation in activated sludge

Activated sludge comprises diverse bacteria, fungi, and other microorganisms, featuring a rich repertoire of genes involved in antibiotic resistance, pollutant degradation, and elemental cycling. In this regard, hybrid assembly technology can revolutionize metagenomics by detecting greater gene diversity in environmental samples. Nonetheless, the optimal utilization and comparability of genomic information between hybrid assembly and short- or long-read technology remain unclear. To address this gap, we compared the performance of the hybrid assembly, short- and long-read technologies, abundance and diversity of annotated genes, and taxonomic diversity by analysing 46, 161, and 45 activated sludge metagenomic datasets, respectively. The results revealed that hybrid assembly technology exhibited the best performance, generating the most contiguous and longest contigs but with a lower proportion of high-quality metagenome-assembled genomes than short-read technology. Compared with short- or long-read technologies, hybrid assembly technology can detect a greater diversity of microbiota and antibiotic resistance genes, as well as a wider range of potential hosts. However, this approach may yield lower gene abundance and pathogen detection. Our study revealed the specific advantages and disadvantages of hybrid assembly and short- and long-read applications in wastewater treatment plants, and our approach could serve as a blueprint to be extended to terrestrial environments.

MiDAS 5: Global diversity of bacteria and archaea in anaerobic digesters

Anaerobic digestion of organic waste into methane and carbon dioxide (biogas) is carried out by complex microbial communities. Here, we use full-length 16S rRNA gene sequencing of 285 full-scale anaerobic digesters (ADs) to expand our knowledge about diversity and function of the bacteria and archaea in ADs worldwide. The sequences are processed into full-length 16S rRNA amplicon sequence variants (FL-ASVs) and are used to expand the MiDAS 4 database for bacteria and archaea in wastewater treatment systems, creating MiDAS 5. The expansion of the MiDAS database increases the coverage for bacteria and archaea in ADs worldwide, leading to improved genus- and species-level classification. Using MiDAS 5, we carry out an amplicon-based, global-scale microbial community profiling of the sampled ADs using three common sets of primers targeting different regions of the 16S rRNA gene in bacteria and/or archaea. We reveal how environmental conditions and biogeography shape the AD microbiota. We also identify core and conditionally rare or abundant taxa, encompassing 692 genera and 1013 species. These represent 84-99% and 18-61% of the accumulated read abundance, respectively, across samples depending on the amplicon primers used. Finally, we examine the global diversity of functional groups with known importance for the anaerobic digestion process.

Advancing metagenome-assembled genome-based pathogen identification: unraveling the power of long-read assembly algorithms in Oxford Nanopore sequencing

Oxford Nanopore sequencing is one of the high-throughput sequencing technologies that facilitates the reconstruction of metagenome-assembled genomes (MAGs). This study aimed to assess the potential of long-read assembly algorithms in Oxford Nanopore sequencing to enhance the MAG-based identification of bacterial pathogens using both simulated and mock communities. Simulated communities were generated to mimic those on fresh spinach and in surface water. Long reads were produced using R9.4.1+SQK-LSK109 and R10.4 + SQK-LSK112, with 0.5, 1, and 2 million reads. The simulated bacterial communities included multidrug-resistant Salmonella enterica serotypes Heidelberg, Montevideo, and Typhimurium in the fresh spinach community individually or in combination, as well as multidrug-resistant Pseudomonas aeruginosa in the surface water community. Real data sets of the ZymoBIOMICS HMW DNA Standard were also studied. A bioinformatic pipeline (MAGenie, freely available at https://github.com/jackchen129/MAGenie) that combines metagenome assembly, taxonomic classification, and sequence extraction was developed to reconstruct draft MAGs from metagenome assemblies. Five assemblers were evaluated based on a series of genomic analyses. Overall, Flye outperformed the other assemblers, followed by Shasta, Raven, and Unicycler, while Canu performed least effectively. In some instances, the extracted sequences resulted in draft MAGs and provided the locations and structures of antimicrobial resistance genes and mobile genetic elements. Our study showcases the viability of utilizing the extracted sequences for precise phylogenetic inference, as demonstrated by the consistent alignment of phylogenetic topology between the reference genome and the extracted sequences. R9.4.1+SQK-LSK109 was more effective in most cases than R10.4+SQK-LSK112, and greater sequencing depths generally led to more accurate results.IMPORTANCEBy examining diverse bacterial communities, particularly those housing multiple Salmonella enterica serotypes, this study holds significance in uncovering the potential of long-read assembly algorithms to improve metagenome-assembled genome (MAG)-based pathogen identification through Oxford Nanopore sequencing. Our research demonstrates that long-read assembly stands out as a promising avenue for boosting precision in MAG-based pathogen identification, thus advancing the development of more robust surveillance measures. The findings also support ongoing endeavors to fine-tune a bioinformatic pipeline for accurate pathogen identification within complex metagenomic samples.

New avenues for potentially seeking microbial responses to climate change beneath Antarctic ice shelves

The signs of climate change are undeniable, and the impact of these changes on ecosystem function heavily depends on the response of microbes that underpin the food web. Antarctic ice shelf is a massive mass of floating ice that extends from the continent into the ocean, exerting a profound influence on global carbon cycles. Beneath Antarctic ice shelves, marine ice stores valuable genetic information, where marine microbial communities before the industrial revolution are archived. Here, in this proof-of-concept, by employing a combination of single-cell technologiesand metagenomics, we have been able to sequence frozen microbial DNA (≈300 years old) stored in the marine ice core B15 collected from the Filchnner-Ronne Ice Shelf. Metagenomic data indicated that Proteobacteria and Thaumarchaeota (e.g., Nitrosopumilus spp.), followed by Actinobacteria (e.g., Actinomarinales), were abundant. Remarkably, our data allow us to "travel to the past" and calibrate genomic and genetic evolutionary changes for ecologically relevant microbes and functions, such as Nitrosopumilus spp., preserved in the marine ice (≈300 years old) with those collected recently in seawater under an ice shelf (year 2017). The evolutionary divergence for the ammonia monooxygenase gene amoA involved in chemolithoautotrophy was about 0.88 amino acid and 2.8 nucleotide substitution rate per 100 sites in a century, while the accumulated rate of genomic SNPs was 2,467 per 1 Mb of genome and 100 years. Whether these evolutionary changes remained constant over the last 300 years or accelerated during post-industrial periods remains an open question that will be further elucidated.

IMPORTANCE

Several efforts have been undertaken to predict the response of microbes under climate change, mainly based on short-term microcosm experiments under forced conditions. A common concern is that manipulative experiments cannot properly simulate the response of microbes to climate change, which is a long-term evolutionary process. In this proof-of-concept study with a limited sample size, we demonstrate a novel approach yet to be fully explored in science for accessing genetic information from putative past marine microbes preserved under Antarctic ice shelves before the industrial revolution. This potentially allows us estimating evolutionary changes as exemplified in our study. We advocate for gathering a more comprehensive Antarctic marine ice core data sets across various periods and sites. Such a data set would enable the establishment of a robust baseline, facilitating a better assessment of the potential effects of climate change on key genetic signatures of microbes.

Improving Bacterial Metagenomic Research through Long-Read Sequencing

Metagenomic sequencing analysis is central to investigating microbial communities in clinical and environmental studies. Short-read sequencing remains the primary approach for metagenomic research; however, long-read sequencing may offer advantages of improved metagenomic assembly and resolved taxonomic identification. To compare the relative performance for metagenomic studies, we simulated short- and long-read datasets using increasingly complex metagenomes comprising 10, 20, and 50 microbial taxa. Additionally, we used an empirical dataset of paired short- and long-read data generated from mouse fecal pellets to assess real-world performance. We compared metagenomic assembly quality, taxonomic classification, and metagenome-assembled genome (MAG) recovery rates. We show that long-read sequencing data significantly improve taxonomic classification and assembly quality. Metagenomic assemblies using simulated long reads were more complete and more contiguous with higher rates of MAG recovery. This resulted in more precise taxonomic classifications. Principal component analysis of empirical data demonstrated that sequencing technology affects compositional results as samples clustered by sequence type, not sample type. Overall, we highlight strengths of long-read metagenomic sequencing for microbiome studies, including improving the accuracy of classification and relative abundance estimates. These results will aid researchers when considering which sequencing approaches to use for metagenomic projects.

Putative genome contamination has minimal impact on the GTDB taxonomy

The Genome Taxonomy Database (GTDB) provides a species to domain classification of publicly available genomes based on average nucleotide identity (ANI) (for species) and a concatenated gene phylogeny normalized by evolutionary rates (for genus to phylum), which has been widely adopted by the scientific community. Here, we use the Genome UNClutterer (GUNC) software to identify putatively contaminated genomes in GTDB release 07-RS207. We found that GUNC reported 35,723 genomes as putatively contaminated, comprising 11.25 % of the 317,542 genomes in GTDB release 07-RS207. To assess the impact of this high level of inferred contamination on the delineation of taxa, we created 'clean' versions of the 34,846 putatively contaminated bacterial genomes by removing the most contaminated half. For each clean half, we re-calculated the ANI and concatenated gene phylogeny and found that only 77 (0.22 %) of the genomes were not consistent with their original classification. We conclude that the delineation of taxa in GTDB is robust to the putative contamination detected by GUNC.

A multifaceted investigation of lactococcal strain diversity in undefined mesophilic starter cultures

The dairy fermentation industry relies on the activity of lactic acid bacteria in robust starter cultures to accomplish milk acidification. Maintenance of the composition of these starter cultures, whether defined or undefined, is essential to ensure consistent and high-quality fermentation end products. To date, limited information exists regarding the microbial composition of undefined starter culture systems. Here, we describe a culture-based analysis combined with a metagenomics approach to evaluate the composition of two undefined mesophilic starter cultures. In addition, we describe a qPCR-based genotype detection assay, which is capable of discerning nine distinct lactococcal genotypes to characterize these undefined starter cultures, and which can be applied to monitor compositional changes in an undefined starter culture during a fermentation.

IMPORTANCE

This study reports on the development of a combined culture-based analysis and metagenomics approach to evaluate the composition of two undefined mesophilic starter cultures. In addition, a novel qPCR-based genotype detection assay, capable of discerning nine distinct lactococcal genotypes (based on lactococcal cell wall polysaccharide biosynthesis gene clusters), was used to monitor compositional changes in an undefined starter culture following phage attack. These analytical approaches facilitate a multifaceted assessment of starter culture compositional stability during milk fermentation, which has become an important QC aspect due to the increasing demand for consistent and high-quality dairy products.

Global abundance patterns, diversity, and ecology of Patescibacteria in wastewater treatment plants

BACKGROUND

Microorganisms are responsible for nutrient removal and resource recovery in wastewater treatment plants (WWTPs), and their diversity is often studied by 16S rRNA gene amplicon sequencing. However, this approach underestimates the abundance and diversity of Patescibacteria due to the low coverage of commonly used PCR primers for this highly divergent bacterial phylum. Therefore, our current understanding of the global diversity, distribution, and ecological role of Patescibacteria in WWTPs is very incomplete. This is particularly relevant as Patescibacteria are considered to be associated with microbial host cells and can therefore influence the abundance and temporal variability of other microbial groups that are important for WWTP functioning.

RESULTS

Here, we evaluated the in silico coverage of widely used 16S rRNA gene-targeted primer pairs and redesigned a primer pair targeting the V4 region of bacterial and archaeal 16S rRNA genes to expand its coverage for Patescibacteria. We then experimentally evaluated and compared the performance of the original and modified V4-targeted primers on 565 WWTP samples from the MiDAS global sample collection. Using the modified primer pair, the percentage of ASVs classified as Patescibacteria increased from 5.9 to 23.8%, and the number of detected patescibacterial genera increased from 560 to 1576, while the detected diversity of the remaining microbial community remained similar. Due to this significantly improved coverage of Patescibacteria, we identified 23 core genera of Patescibacteria in WWTPs and described the global distribution pattern of these unusual microbes in these systems. Finally, correlation network analysis revealed potential host organisms that might be associated with Patescibacteria in WWTPs. Interestingly, strong indications were found for an association between Patescibacteria of the Saccharimonadia and globally abundant polyphosphate-accumulating organisms of the genus Ca. Phosphoribacter.

CONCLUSIONS

Our study (i) provides an improved 16S rRNA gene V4 region-targeted amplicon primer pair inclusive of Patescibacteria with little impact on the detection of other taxa, (ii) reveals the diversity and distribution patterns of Patescibacteria in WWTPs on a global scale, and (iii) provides new insights into the ecological role and potential hosts of Patescibacteria in WWTPs. Video Abstract.

Discovery of cytosine deaminases enables base-resolution methylome mapping using a single enzyme

Deaminases have important uses in modification detection and genome editing. However, the range of applications is limited by the small number of characterized enzymes. To expand the toolkit of deaminases, we developed an in vitro approach that bypasses a major hurdle with their toxicity in cells. We assayed 175 putative cytosine deaminases on a variety of substrates and found a broad range of activity on double- and single-stranded DNA in various sequence contexts, including CpG-specific deaminases and enzymes without sequence preference. We also characterized enzyme selectivity across six DNA modifications and reported enzymes that do not deaminate modified cytosines. The detailed analysis of diverse deaminases opens new avenues for biotechnological and medical applications. As a demonstration, we developed SEM-seq, a non-destructive single-enzyme methylation sequencing method using a modification-sensitive double-stranded DNA deaminase. The streamlined protocol enables accurate, base-resolution methylome mapping of scarce biological material, including cell-free DNA and 10 pg input DNA.

Genome-centric analyses of 165 metagenomes show that mobile genetic elements are crucial for the transmission of antimicrobial resistance genes to pathogens in activated sludge and wastewater

Wastewater is considered a reservoir of antimicrobial resistance genes (ARGs), where the abundant antimicrobial-resistant bacteria and mobile genetic elements facilitate horizontal gene transfer. However, the prevalence and extent of these phenomena in different taxonomic groups that inhabit wastewater are still not fully understood. Here, we determined the presence of ARGs in metagenome-assembled genomes (MAGs) and evaluated the risks of MAG-carrying ARGs in potential human pathogens. The potential of these ARGs to be transmitted horizontally or vertically was also determined. A total of 5,916 MAGs (completeness >50%, contamination <10%) were recovered, covering 68 phyla and 279 genera. MAGs were dereplicated into 1,204 genome operational taxonomic units (gOTUs) as a proxy for species ( average nucleotide identity >0.95). The dominant ARG classes detected were bacitracin, multi-drug, macrolide-lincosamide-streptogramin (MLS), glycopeptide, and aminoglycoside, and 10.26% of them were located on plasmids. The main hosts of ARGs belonged to Escherichia, Klebsiella, Acinetobacter, Gresbergeria, Mycobacterium, and Thauera. Our data showed that 253 MAGs carried virulence factor genes (VFGs) divided into 44 gOTUs, of which 45 MAGs were carriers of ARGs, indicating that potential human pathogens carried ARGs. Alarmingly, the MAG assigned as Escherichia coli contained 159 VFGs, of which 95 were located on chromosomes and 10 on plasmids. In addition to shedding light on the prevalence of ARGs in individual genomes recovered from activated sludge and wastewater, our study demonstrates a workflow that can identify antimicrobial-resistant pathogens in complex microbial communities.

IMPORTANCE

Antimicrobial resistance (AMR) threatens the health of humans, animals, and natural ecosystems. In our study, an analysis of 165 metagenomes from wastewater revealed antibiotic-targeted alteration, efflux, and inactivation as the most prevalent AMR mechanisms. We identified several genera correlated with multiple ARGs, including Klebsiella, Escherichia, Acinetobacter, Nitrospira, Ottowia, Pseudomonas, and Thauera, which could have significant implications for AMR transmission. The abundance of bacA, mexL, and aph(3")-I in the genomes calls for their urgent management in wastewater. Our approach could be applied to different ecosystems to assess the risk of potential pathogens containing ARGs. Our findings highlight the importance of managing AMR in wastewater and can help design measures to reduce the transmission and evolution of AMR in these systems.

Horizontal transfer potential of antibiotic resistance genes in wastewater treatment plants unraveled by microfluidic-based mini-metagenomics

Wastewater treatment plants (WWTPs) are known to harbor antibiotic resistance genes (ARGs), which can potentially spread to the environment and human populations. However, the extent and mechanisms of ARG transfer in WWTPs are not well understood due to the high microbial diversity and limitations of molecular techniques. In this study, we used a microfluidic-based mini-metagenomics approach to investigate the transfer potential and mechanisms of ARGs in activated sludge from WWTPs. Our results show that while diverse ARGs are present in activated sludge, only a few highly similar ARGs are observed across different taxa, indicating limited transfer potential. We identified two ARGs, ermF and tla-1, which occur in a variety of bacterial taxa and may have high transfer potential facilitated by mobile genetic elements. Interestingly, genes that are highly similar to the sequences of these two ARGs, as identified in this study, display varying patterns of abundance across geographic regions. Genes similar to ermF found are widely found in Asia and the Americas, while genes resembling tla-1 are primarily detected in Asia. Genes similar to both genes are barely detected in European WWTPs. These findings shed light on the limited horizontal transfer potential of ARGs in WWTPs and highlight the importance of monitoring specific ARGs in different regions to mitigate the spread of antibiotic resistance.

Functional traits and health implications of the global household drinking-water microbiome retrieved using an integrative genome-centric approach

The biological safety of drinking water plays a crucial role in public health protection. However, research on the drinking water microbiome remains in its infancy, especially little is known about the potentially pathogenic bacteria in and functional characteristics of the microbiome in household tap water that people are directly exposed to. In this study, we used a genomic-centric approach to construct a genetic catalogue of the drinking water microbiome by analysing 116 metagenomic datasets of household tap water worldwide, spanning nine countries/regions on five continents. We reconstructed 859 high-quality metagenome-assembled genomes (MAGs) spanning 27 bacterial and 2 archaeal phyla, and found that the core MAGs belonging to the phylum Proteobacteria encoded the highest metabolic functional diversity of the 33 key complete metabolic modules. In particular, we found that two core MAGs of Brevibacillus and Methylomona encoded genes for methane metabolism, which may support the growth of heterotrophic organisms observed in the oligotrophic ecosystem. Four MAGs of complete ammonia oxidation (comammox) Nitrospira were identified and functional metabolic analysis suggested these may enable mixotrophic growth and encode genes for reactive oxygen stress defence and arsenite reduction that could aid survival in the environment of oligotrophic drinking water systems. Four MAGs were annotated as potentially pathogenic bacteria (PPB) and thus represented a possible public health concern. They belonged to the genera Acinetobacter (n = 3) and Mycobacterium (n = 1), with a total relative abundance of 1.06 % in all samples. The genomes of PPB A. junii and A. ursingii were discovered to contain antibiotic resistance genes and mobile genetic elements that could contribute to antimicrobial dissemination in drinking water. Further network analysis suggested that symbiotic microbes which support the growth of pathogenic bacteria can be targets for future surveillance and removal.

From metagenomes to metabolism: Systematically assessing the metabolic flux feasibilities for "Candidatus Accumulibacter" species during anaerobic substrate uptake

With the rapid growing availability of metagenome assembled genomes (MAGs) and associated metabolic models, the identification of metabolic potential in individual community members has become possible. However, the field still lacks an unbiassed systematic evaluation of the generated metagenomic information to uncover not only metabolic potential, but also feasibilities of these models under specific environmental conditions. In this study, we present a systematic analysis of the metabolic potential in species of "Candidatus Accumulibacter", a group of polyphosphate-accumulating organisms (PAOs). We constructed a metabolic model of the central carbon metabolism and compared the metabolic potential among available MAGs for "Ca. Accumulibacter" species. By combining Elementary Flux Modes Analysis (EFMA) with max-min driving force (MDF) optimization, we obtained all possible flux distributions of the metabolic network and calculated their individual thermodynamic feasibility. Our findings reveal significant variations in the metabolic potential among "Ca. Accumulibacter" MAGs, particularly in the presence of anaplerotic reactions. EFMA revealed 700 unique flux distributions in the complete metabolic model that enable the anaerobic uptake of acetate and its conversion into polyhydroxyalkanoates (PHAs), a well-known phenotype of "Ca. Accumulibacter". However, thermodynamic constraints narrowed down this solution space to 146 models that were stoichiometrically and thermodynamically feasible (MDF > 0 kJ/mol), of which only 8 were strongly feasible (MDF > 7 kJ/mol). Notably, several novel flux distributions for the metabolic model were identified, suggesting putative, yet unreported, functions within the PAO communities. Overall, this work provides valuable insights into the metabolic variability among "Ca. Accumulibacter" species and redefines the anaerobic metabolic potential in the context of phosphate removal. More generally, the integrated workflow presented in this paper can be applied to any metabolic model obtained from a MAG generated from microbial communities to objectively narrow the expected phenotypes from community members.

Membrane aerated biofilm reactor system driven by pure oxygen for wastewater treatment

Pure oxygen is proposed for wastewater treatment due to its advantages over conventional air aeration. This study investigates a Pure Oxygen-based Membrane Aerated Biofilm Reactor (PO-MABR) for the first time under various operating conditions. The PO-MABR employs a gas-permeable membrane for direct diffusion of low-pressurized pure oxygen to the biofilm, ensuring exceptional carbon and nitrogen removal. The effectiveness of PO-MABR was investigated by varying operational conditions, including temperature, carbon-to-nitrogen ratio, gas pressure, and flow rate. Results indicate superior performance, with a 97% chemical oxygen demand removal and 19% higher total nitrogen removal than Air-Ventilated MABR (A-MABR) due to thicker biofilm and unique microbial structures in PO-MABR. Also, PO-MABR demonstrated resilience to low temperatures and effectively treated both high and low-strength wastewater. The findings emphasize the efficiency of PO-MABR in wastewater treatment, advocating for its adoption due to superior carbon and nitrogen removal across diverse operational conditions.

Closing the genome of unculturable cable bacteria using a combined metagenomic assembly of long and short sequencing reads

Many environmentally relevant micro-organisms cannot be cultured, and even with the latest metagenomic approaches, achieving complete genomes for specific target organisms of interest remains a challenge. Cable bacteria provide a prominent example of a microbial ecosystem engineer that is currently unculturable. They occur in low abundance in natural sediments, but due to their capability for long-distance electron transport, they exert a disproportionately large impact on the biogeochemistry of their environment. Current available genomes of marine cable bacteria are highly fragmented and incomplete, hampering the elucidation of their unique electrogenic physiology. Here, we present a metagenomic pipeline that combines Nanopore long-read and Illumina short-read shotgun sequencing. Starting from a clonal enrichment of a cable bacterium, we recovered a circular metagenome-assembled genome (5.09 Mbp in size), which represents a novel cable bacterium species with the proposed name Candidatus Electrothrix scaldis. The closed genome contains 1109 novel identified genes, including key metabolic enzymes not previously described in incomplete genomes of cable bacteria. We examined in detail the factors leading to genome closure. Foremost, native, non-amplified long reads are crucial to resolve the many repetitive regions within the genome of cable bacteria, and by analysing the whole metagenomic assembly, we found that low strain diversity is key for achieving genome closure. The insights and approaches presented here could help achieve genome closure for other keystone micro-organisms present in complex environmental samples at low abundance.

Comparative genomic analyses of Cutibacterium granulosum provide insights into genomic diversity

Cutibacterium granulosum, a commensal bacterium found on human skin, formerly known as Propionibacterium granulosum, rarely causes infections and is generally considered non-pathogenic. Recent research has revealed the transferability of the multidrug-resistant plasmid pTZC1 between C. granulosum and Cutibacterium acnes, the latter being an opportunistic pathogen in surgical site infections. However, there is a noticeable lack of research on the genome of C. granulosum, and the genetic landscape of this species remains largely uncharted. We investigated the genomic features and evolutionary structure of C. granulosum by analyzing a total of 30 Metagenome-Assembled Genomes (MAGs) and isolate genomes retrieved from public databases, as well as those generated in this study. A pan-genome of 6,077 genes was identified for C. granulosum. Remarkably, the 'cloud genes' constituted 62.38% of the pan-genome. Genes associated with mobilome: prophages, transposons [X], defense mechanisms [V] and replication, recombination and repair [L] were enriched in the cloud genome. Phylogenomic analysis revealed two distinct mono-clades, highlighting the genomic diversity of C. granulosum. The genomic diversity was further confirmed by the distribution of Average Nucleotide Identity (ANI) values. The functional profiles analysis of C. granulosum unveiled a wide range of potential Antibiotic Resistance Genes (ARGs) and virulence factors, suggesting its potential tolerance to various environmental challenges. Subtype I-E of the CRISPR-Cas system was the most abundant in these genomes, a feature also detected in C. acnes genomes. Given the widespread distribution of C. granulosum strains within skin microbiome, our findings make a substantial contribution to our broader understanding of the genetic diversity, which may open new avenues for investigating the mechanisms and treatment of conditions such as acne vulgaris.

Niche differentiation in microbial communities with stable genomic traits over time in engineered systems

Microbial communities in full-scale engineered systems undergo dynamic compositional changes. However, mechanisms governing assembly of such microbes and succession of their functioning and genomic traits under various environmental conditions are unclear. In this study, we used the activated sludge and anaerobic treatment systems of four full-scale industrial wastewater treatment plants as models to investigate the niches of microbes in communities and the temporal succession patterns of community compositions. High-quality representative metagenome-assembled genomes revealed that taxonomic, functional, and trait-based compositions were strongly shaped by environmental selection, with replacement processes primarily driving variations in taxonomic and functional compositions. Plant-specific indicators were associated with system environmental conditions and exhibited strong determinism and trajectory directionality over time. The partitioning of microbes in a co-abundance network according to groups of plant-specific indicators, together with significant between-group differences in genomic traits, indicated the occurrence of niche differentiation. The indicators of the treatment plant with rich nutrient input and high substrate removal efficiency exhibited a faster predicted growth rate, lower guanine-cytosine content, smaller genome size, and higher codon usage bias than the indicators of the other plants. In individual plants, taxonomic composition displayed a more rapid temporal succession than functional and trait-based compositions. The succession of taxonomic, functional, and trait-based compositions was correlated with the kinetics of treatment processes in the activated sludge systems. This study provides insights into ecological niches of microbes in engineered systems and succession patterns of their functions and traits, which will aid microbial community management to improve treatment performance.

Two new clades recovered at high temperatures provide novel phylogenetic and genomic insights into Candidatus Accumulibacter

Candidatus Accumulibacter, a key genus of polyphosphate-accumulating organisms, plays key roles in lab- and full-scale enhanced biological phosphorus removal (EBPR) systems. A total of 10 high-quality Ca. Accumulibacter genomes were recovered from EBPR systems operated at high temperatures, providing significantly updated phylogenetic and genomic insights into the Ca. Accumulibacter lineage. Among these genomes, clade IIF members SCELSE-3, SCELSE-4, and SCELSE-6 represent the to-date known genomes encoding a complete denitrification pathway, suggesting that Ca. Accumulibacter alone could achieve complete denitrification. Clade IIC members SSA1, SCUT-1, SCELCE-2, and SCELSE-8 lack the entire set of denitrifying genes, representing to-date known non-denitrifying Ca. Accumulibacter. A pan-genomic analysis with other Ca. Accumulibacter members suggested that all Ca. Accumulibacter likely has the potential to use dicarboxylic amino acids. Ca. Accumulibacter aalborgensis AALB and Ca. Accumulibacter affinis BAT3C720 seemed to be the only two members capable of using glucose for EBPR. A heat shock protein Hsp20 encoding gene was found exclusively in genomes recovered at high temperatures, which was absent in clades IA, IC, IG, IIA, IIB, IID, IIG, and II-I members. High transcription of this gene in clade IIC members SCUT-2 and SCUT-3 suggested its role in surviving high temperatures for Ca. Accumulibacter. Ambiguous clade identity was observed for newly recovered genomes (SCELSE-9 and SCELSE-10). Five machine learning models were developed using orthogroups as input features. Prediction results suggested that they belong to a new clade (IIK). The phylogeny of Ca. Accumulibacter was re-evaluated based on the laterally derived polyphosphokinase 2 gene, showing improved resolution in differentiating different clades.

Disentangling the effects of sulfate and other seawater ions on microbial communities and greenhouse gas emissions in a coastal forested wetland

Seawater intrusion into freshwater wetlands causes changes in microbial communities and biogeochemistry, but the exact mechanisms driving these changes remain unclear. Here we use a manipulative laboratory microcosm experiment, combined with DNA sequencing and biogeochemical measurements, to tease apart the effects of sulfate from other seawater ions. We examined changes in microbial taxonomy and function as well as emissions of carbon dioxide, methane, and nitrous oxide in response to changes in ion concentrations. Greenhouse gas emissions and microbial richness and composition were altered by artificial seawater regardless of whether sulfate was present, whereas sulfate alone did not alter emissions or communities. Surprisingly, addition of sulfate alone did not lead to increases in the abundance of sulfate reducing bacteria or sulfur cycling genes. Similarly, genes involved in carbon, nitrogen, and phosphorus cycling responded more strongly to artificial seawater than to sulfate. These results suggest that other ions present in seawater, not sulfate, drive ecological and biogeochemical responses to seawater intrusion and may be drivers of increased methane emissions in soils that received artificial seawater addition. A better understanding of how the different components of salt water alter microbial community composition and function is necessary to forecast the consequences of coastal wetland salinization.

Amplicon-guided isolation and cultivation of previously uncultured microbial species from activated sludge

IMPORTANCE

Biological wastewater treatment relies on complex microbial communities that assimilate nutrients and break down pollutants in the wastewater. Knowledge about the physiology and metabolism of bacteria in wastewater treatment plants (WWTPs) may therefore be used to improve the efficacy and economy of wastewater treatment. Our current knowledge is largely based on 16S rRNA gene amplicon profiling, fluorescence in situ hybridization studies, and predictions based on metagenome-assembled genomes. Bacterial isolates are often required to validate genome-based predictions as they allow researchers to analyze a specific species without interference from other bacteria and with simple bulk measurements. Unfortunately, there are currently very few pure cultures representing the microbes commonly found in WWTPs. To address this, we introduce an isolation strategy that takes advantage of state-of-the-art microbial profiling techniques to uncover suitable growth conditions for key WWTP microbes. We furthermore demonstrate that this information can be used to isolate key organisms representing global WWTPs.

A comprehensive overview of the Chloroflexota community in wastewater treatment plants worldwide

IMPORTANCE

Chloroflexota are often abundant members of the biomass in wastewater treatment plants (WWTPs) worldwide, typically with a filamentous morphology, forming the backbones of the activated sludge floc. However, their overgrowth can often cause operational issues connected to poor settling or foaming, impairing effluent quality and increasing operational costs. Despite their importance, few Chloroflexota genera have been characterized so far. Here, we present a comprehensive overview of Chloroflexota abundant in WWTPs worldwide and an in-depth characterization of their morphology, phylogeny, and ecophysiology, obtaining a broad understanding of their ecological role in activated sludge.

Dancing the Nanopore limbo - Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BACKGROUND

While genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.

RESULTS

We generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.

CONCLUSION

Our results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

First shotgun metagenomics study of Juan de Fuca deep-sea sediments reveals distinct microbial communities above, within, between, and below sulfate methane transition zones

The marine deep subsurface is home to a vast microbial ecosystem, affecting biogeochemical cycles on a global scale. One of the better-studied deep biospheres is the Juan de Fuca (JdF) Ridge, where hydrothermal fluid introduces oxidants into the sediment from below, resulting in two sulfate methane transition zones (SMTZs). In this study, we present the first shotgun metagenomics study of unamplified DNA from sediment samples from different depths in this stratified environment. Bioinformatic analyses showed a shift from a heterotrophic, Chloroflexota-dominated community above the upper SMTZ to a chemolithoautotrophic Proteobacteria-dominated community below the secondary SMTZ. The reintroduction of sulfate likely enables respiration and boosts active cells that oxidize acetate, iron, and complex carbohydrates to degrade dead biomass in this low-abundance, low-diversity environment. In addition, analyses showed many proteins of unknown function as well as novel metagenome-assembled genomes (MAGs). The study provides new insights into microbial communities in this habitat, enabled by an improved DNA extraction protocol that allows a less biased view of taxonomic composition and metabolic activities, as well as uncovering novel taxa. Our approach presents the first successful attempt at unamplified shotgun sequencing samples from beyond 50 meters below the seafloor and opens new ways for capturing the true diversity and functional potential of deep-sea sediments.

Meta-omic profiling reveals ubiquity of genes encoding for the nitrogen-rich biopolymer cyanophycin in activated sludge microbiomes

Recovering nitrogen (N) from municipal wastewater is a promising approach to prevent nutrient pollution, reduce energy use, and transition toward a circular N bioeconomy, but remains a technologically challenging endeavor. Existing N recovery techniques are optimized for high-strength, low-volume wastewater. Therefore, developing methods to concentrate dilute N from mainstream wastewater will bridge the gap between existing technologies and practical implementation. The N-rich biopolymer cyanophycin is a promising candidate for N bioconcentration due to its pH-tunable solubility characteristics and potential for high levels of accumulation. However, the cyanophycin synthesis pathway is poorly explored in engineered microbiomes. In this study, we analyzed over 3,700 publicly available metagenome assembled genomes (MAGs) and found that the cyanophycin synthesis gene cphA was ubiquitous across common activated sludge bacteria. We found that cphA was present in common phosphorus accumulating organisms (PAO) Ca. 'Accumulibacter' and Tetrasphaera, suggesting potential for simultaneous N and P bioconcentration in the same organisms. Using metatranscriptomic data, we confirmed the expression of cphA in lab-scale bioreactors enriched with PAO. Our findings suggest that cyanophycin synthesis is a ubiquitous metabolic activity in activated sludge microbiomes. The possibility of combined N and P bioconcentration could lower barriers to entry for N recovery, since P concentration by PAO is already a widespread biotechnology in municipal wastewater treatment. We anticipate this work to be a starting point for future evaluations of combined N and P bioaccumulation, with the ultimate goal of advancing widespread adoption of N recovery from municipal wastewater.

Cross-Feeding between Filamentous Cyanobacteria and Symbiotic Bacteria Favors Rapid Photogranulation

Photogranules are dense algal-bacterial aggregates used in aeration-free and carbon-negative wastewater treatment, wherein filamentous cyanobacteria (FC) are essential components. However, little is known about the functional role of symbiotic bacteria in photogranulation. Herein, we combined cyanobacterial isolation, reactor operation, and multiomics analysis to investigate the cyanobacterial-bacterial interaction during photogranulation. The addition of FC to the inoculated sludge achieved a 1.4-fold higher granule size than the control, and the aggregation capacity of FC-dominant photogranules was closely related to the extracellular polysaccharide (PS) concentration (R = 0.86). Importantly, we found that cross-feeding between FC and symbiotic bacteria for macromolecular PS synthesis is at the heart of photogranulation and substantially enhanced the granular stability. Chloroflexi-affiliated bacteria intertwined with FC throughout the photogranules and promoted PS biosynthesis using the partial nucleotide sugars produced by FC. Proteobacteria-affiliated bacteria were spatially close to FC, and highly expressed genes for vitamin B1 and B12 synthesis, contributing the necessary cofactors to promote FC proliferation. In addition, Bacteroidetes-affiliated bacteria degraded FC-derived carbohydrates and influenced granules development. Our metabolic characterization identified the functional role of symbiotic bacteria of FC during photogranulation and shed light on the critical cyanobacterial-bacterial interactions in photogranules from the viewpoint of cross-feeding.

Metaproteomics, metagenomics and 16S rRNA sequencing provide different perspectives on the aerobic granular sludge microbiome

The tremendous progress in sequencing technologies has made DNA sequencing routine for microbiome studies. Additionally, advances in mass spectrometric techniques have extended conventional proteomics into the field of microbial ecology. However, systematic studies that provide a better understanding of the complementary nature of these 'omics' approaches, particularly for complex environments such as wastewater treatment sludge, are urgently needed. Here, we describe a comparative metaomics study on aerobic granular sludge from three different wastewater treatment plants. For this, we employed metaproteomics, whole metagenome, and 16S rRNA amplicon sequencing to study the same granule material with uniform size. We furthermore compare the taxonomic profiles using the Genome Taxonomy Database (GTDB) to enhance the comparability between the different approaches. Though the major taxonomies were consistently identified in the different aerobic granular sludge samples, the taxonomic composition obtained by the different omics techniques varied significantly at the lower taxonomic levels, which impacts the interpretation of the nutrient removal processes. Nevertheless, as demonstrated by metaproteomics, the genera that were consistently identified in all techniques cover the majority of the protein biomass. The established metaomics data and the contig classification pipeline are publicly available, which provides a valuable resource for further studies on metabolic processes in aerobic granular sludge.

Adversarial and variational autoencoders improve metagenomic binning

Assembly of reads from metagenomic samples is a hard problem, often resulting in highly fragmented genome assemblies. Metagenomic binning allows us to reconstruct genomes by re-grouping the sequences by their organism of origin, thus representing a crucial processing step when exploring the biological diversity of metagenomic samples. Here we present Adversarial Autoencoders for Metagenomics Binning (AAMB), an ensemble deep learning approach that integrates sequence co-abundances and tetranucleotide frequencies into a common denoised space that enables precise clustering of sequences into microbial genomes. When benchmarked, AAMB presented similar or better results compared with the state-of-the-art reference-free binner VAMB, reconstructing ~7% more near-complete (NC) genomes across simulated and real data. In addition, genomes reconstructed using AAMB had higher completeness and greater taxonomic diversity compared with VAMB. Finally, we implemented a pipeline Integrating VAMB and AAMB that enabled improved binning, recovering 20% and 29% more simulated and real NC genomes, respectively, compared to VAMB, with moderate additional runtime.

Metagenome sequencing and recovery of 444 metagenome-assembled genomes from the biofloc aquaculture system

Biofloc technology is increasingly recognised as a sustainable aquaculture method. In this technique, bioflocs are generated as microbial aggregates that play pivotal roles in assimilating toxic nitrogenous substances, thereby ensuring high water quality. Despite the crucial roles of the floc-associated bacterial (FAB) community in pathogen control and animal health, earlier microbiota studies have primarily relied on the metataxonomic approaches. Here, we employed shotgun sequencing on eight biofloc metagenomes from a commercial aquaculture system. This resulted in the generation of 106.6 Gbp, and the reconstruction of 444 metagenome-assembled genomes (MAGs). Among the recovered MAGs, 230 were high-quality (≥90% completeness, ≤5% contamination), and 214 were medium-quality (≥50% completeness, ≤10% contamination). Phylogenetic analysis unveiled Rhodobacteraceae as dominant members of the FAB community. The reported metagenomes and MAGs are crucial for elucidating the roles of diverse microorganisms and their functional genes in key processes such as nitrification, denitrification, and remineralization. This study will contribute to scientific understanding of phylogenetic diversity and metabolic capabilities of microbial taxa in aquaculture environments.

Unraveling the rate-limiting step in microorganisms' mediation of denitrification and phosphorus absorption/transport processes in a highly regulated river-lake system

River-lake ecosystems are indispensable hubs for water transfers and flow regulation engineering, which have frequent and complex artificial hydrological regulation processes, and the water quality is often unstable. Microorganisms usually affect these systems by driving the nutrient cycling process. Thus, understanding the key biochemical rate-limiting steps under highly regulated conditions was critical for the water quality stability of river-lake ecosystems. This study investigated how the key microorganisms and genes involving nitrogen and phosphorus cycling contributed to the stability of water by combining 16S rRNA and metagenomic sequencing using the Dongping river-lake system as the case study. The results showed that nitrogen and phosphorus concentrations were significantly lower in lake zones than in river inflow and outflow zones (p < 0.05). Pseudomonas, Acinetobacter, and Microbacterium were the key microorganisms associated with nitrate and phosphate removal. These microorganisms contributed to key genes that promote denitrification (nirB/narG/narH/nasA) and phosphorus absorption and transport (pstA/pstB/pstC/pstS). Partial least squares path modeling (PLS-PM) revealed that environmental factors (especially flow velocity and COD concentration) have a significant negative effect on the key microbial abundance (p < 0.001). Our study provides theoretical support for the effective management and protection of water transfer and the regulation function of the river-lake system.

Primary settling changes the microbial community of influent wastewater to wastewater treatment plants

The continuous immigration of bacteria in influent wastewater strongly impacts the microbial community of activated sludge (AS) in wastewater treatment plants (WWTP), both in terms of species composition and their abundance. Therefore, it is of interest to elucidate the route of immigrating bacteria into the biological tanks, including the effect of primary settlers. These are commonly used pretreatment units that can possibly selectively increase or reduce the relative abundance of certain bacteria. Species-level identification of the microbial composition of influent wastewater before and after primary settling was carried out in four full-scale municipal WWTPs biweekly over one year by 16S rRNA gene amplicon sequencing. Overall, 37-49% of incoming COD was removed in the primary settlers. Most genera and species were present in the wastewater to all four plants and the trend of these were investigated across the primary settlers. Approximately 50% of the genera had the same trend across at least three WWTPs. Few genera significantly increased in relative read abundance (3.7%) after settling, while 22.3% showed a significant reduction in relative abundance. We investigated process-critical species in AS, such as known nitrifiers, polyphosphate-accumulating organisms, and filamentous bacteria. Most taxa were affected similarly in all WWTPs including multiple genera involved in bulking in AS. However, some genera, e.g., important polyphosphate-accumulating bacteria, had inconsistent trends across WWTPs, suggesting that the characteristics of the wastewater are important for the trend of some bacteria through primary settling. In all cases, primary settling changed the microbial community of the influent wastewater, posing an obvious candidate for upstream control to optimize the assembly of the microbial communities in activated sludge.

Treatment performance and microbial community structure in an aerobic granular sludge sequencing batch reactor amended with diclofenac, erythromycin, and gemfibrozil

This study characterizes the effects of three commonly detected pharmaceuticals-diclofenac, erythromycin, and gemfibrozil-on aerobic granular sludge. Approximately 150 μg/L of each pharmaceutical was fed in the influent to a sequencing batch reactor for 80 days, and the performance of the test reactor was compared with that of a control reactor. Wastewater treatment efficacy in the test reactor dropped by approximately 30-40%, and ammonia oxidation was particularly inhibited. The relative abundance of active Rhodocyclaceae, Nitrosomonadaceae, and Nitrospiraceae families declined throughout exposure, likely explaining reductions in wastewater treatment performance. Pharmaceuticals were temporarily removed in the first 12 days of the test via both sorption and degradation; both removal processes declined sharply thereafter. This study demonstrates that aerobic granular sludge may successfully remove pharmaceuticals in the short term, but long-term tests are necessary to confirm if pharmaceutical removal is sustainable.

ACR: metagenome-assembled prokaryotic and eukaryotic genome refinement tool

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k-means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics.

Microaerophilic Activated Sludge System for Ammonia Retention toward Recovery from High-Strength Nitrogenous Wastewater: Performance and Microbial Communities

A transition to ammonia recovery from wastewater has started; however, a technology for sustainable nitrogen retention in the form of ammonia and organic carbon removal is still in development. This study validated a microaerophilic activated sludge (MAS) system to efficiently retain ammonia from high-strength nitrogenous wastewater. The MAS is based on conventional activated sludge (CAS) with aerobic and settling compartments. Low dissolved oxygen (DO) concentrations (<0.2 mg/L) and short solids retention times (SRTs) (<5 days) eliminated nitrifying bacteria. The two parallel MASs were successfully operated for 300 days and had ammonia retention of 101.7 ± 24.9% and organic carbon removal of 85.5 ± 8.9%. The MASs mitigated N2O emissions with an emission factor of <0.23%, much lower than the default value of CAS (1.6%). A short-term step-change test demonstrated that N2O indicated the initiation of nitrification and the completion of denitrification in the MAS. The parallel MASs had comparable microbial diversity, promoting organic carbon oxidation while inhibiting ammonia-oxidizing microorganisms (AOMs), as revealed by 16S rRNA gene amplicon sequencing, the quantitative polymerase chain reaction of functional genes, and fluorescence in situ hybridization of β-proteobacteria AOB. The microbial analyses also uncovered that filamentous bacteria were positively correlated with effluent turbidity. Together, controlling DO and SRT achieved organic carbon removal and successful ammonia retention, mainly by suppressing AOM activity. This process represents a new nitrogen management paradigm.

Microbiome assembly mechanism and functional potential in enhanced biological phosphorus removal system enriched with Tetrasphaera-related polyphosphate accumulating organisms

Tetrasphaera-related polyphosphate accumulating organisms (PAOs) are the key functional guilds for enhanced biological phosphorus removal (EBPR) systems. Their biomass enrichment can be enhanced by the nitrification inhibitor allylthiourea (ATU). However, the underlying assembly mechanism and the functional potential of the EBPR microbiome regulated by ATU are unclear. This study investigates the effect of ATU on microbiome assembly and functional potential by closely following the microbiota dynamics in an EBPR system enriched with Tetrasphaera-related PAOs for 288-days before, during and after ATU addition. The results showed that ATU addition increased microbiota structural similarity and compositional convergence, and enhanced determinism in the assembly of EBPR microbiome. During exposure to ATU, Tetrasphaera-related PAOs were governed by homogeneous selection and the dominant species revealed by 16S rRNA gene-based phylogenetic analysis shifted from clade III to clade I. Meanwhile, ATU supply promoted significant enrichment of functional genes involved in phosphate transport (pit) and polyphosphate synthesis and degradation (ppk1 and ppk2), whereas both Nitrosomonas and ammonia monooxygenase-encoding genes (amoA/B/C) assignable to this group of nitrifying bacteria decreased. Moreover, ATU addition relieved the significant abundance correlation between filamentous bacteria Ca. Promineofilum and denitrifying Brevundimonas (FDR-adjusted P < 0.01), damaging their potential synergic or cooperative interactions, thus weakening their competitiveness against Tetrasphaera-related PAOs. Notably, ATU withdrawn created opportunistic conditions for the unexpected explosive growth and predominance of Thiothrix filaments, leading to a serious bulking event. Our study provides new insights into the microbial ecology of Tetrasphaera-related PAOs in EBPR system, which could guide the establishment of an efficient microbiota for EBPR.

Identification of key degraders for controlling toxicity risks of disguised toxic pollutants with division of labor mechanisms in activated sludge microbiomes: Using nonylphenol ethoxylate as an example

Efficient management of disguised toxic pollutants (DTPs), which can undergo microbial degradation and convert into more toxic substances, necessitates the collaboration of diverse microbial populations in wastewater treatment plants. However, the identification of key bacterial degraders capable of controlling the toxicity risks of DTPs through division of labor mechanisms in activated sludge microbiomes has received limited attention. In this study, we investigated the key degraders capable of controlling the risk of estrogenicity associated with nonylphenol ethoxylate (NPEO), a representative DTP, in textile activated sludge microbiomes. The results of our batch experiments revealed that the transformation of NPEO into NP and subsequent NP degradation were the rate-limiting processes for controlling the risk of estrogenicity, resulting in an inverted V-shaped curve of estrogenicity in water samples during the biodegradation of NPEO by textile activated sludge. By utilizing enrichment sludge microbiomes treated with NPEO or NP as the sole carbon and energy source, a total of 15 bacterial degraders, including Sphingbium, Pseudomonas, Dokdonella, Comamonas, and Hyphomicrobium, were identified as capable of participating in these processes, Among them, Sphingobium and Pseudomonas were the two key degraders that could cooperatively interact in the degradation of NPEO with division of labor mechanisms. Co-culturing Sphingobium and Pseudomonas isolates exhibited a synergistic effect in degrading NPEO and reducing estrogenicity. Our study underscores the potential of the identified functional bacteria for controlling estrogenicity associated with NPEO and provides a methodological framework for identifying key cooperators engaged in labor division, contributing to the management of risks associated with DTPs by leveraging intrinsic microbial metabolic interactions.

Metabolic implications for predatory and parasitic bacterial lineages in activated sludge wastewater treatment systems

Deciphering unclear microbial interactions is key to improving biological wastewater treatment processes. Microbial predation and parasitism in wastewater treatment ecosystems are unexplored survival strategies that have long been known and have recently attracted attention because these interspecies interactions may contribute to the reduction of excess sludge. Here, microbial community profiling of 600 activated sludge samples taken from six industrial and one municipal wastewater treatment processes (WWTPs) was conducted. To identify the shared lineages in the WWTPs, the shared microbial constituents were defined as the family level taxa that had ≥ 0.1% average relative abundance and detected in all processes. The microbial community analysis assigned 106 families as the shared microbial constituents in the WWTPs. Correlation analysis showed that 98 of the 106 shared families were significantly correlated with total carbon (TC) and/or total nitrogen (TN) concentrations, suggesting that they may contribute to wastewater remediation. Most possible predatory or parasitic bacteria belonging to the phyla Bdellovibrionota, Myxococcota, and Candidatus Patescibacteria were found to be the shared families and negatively correlated with TC/TN; thus, they were frequently present in the WWTPs and could be involved in the removal of carbon/nitrogen derived from cell components. Shotgun metagenome-resolved metabolic reconstructions indicated that gene homologs associated with predation or parasitism are conserved in the Bdellovibrionota, Myxococcota, and Ca. Patescibacteria genomes (e.g., host interaction (hit) locus, Tad-like secretion complexes, and type IV pilus assembly proteins). This study provides insights into the complex microbial interactions potentially linked to the reduction of excess sludge biomass in these processes.

Biochar Enhances Soil Resource Availability and Suppresses Microbial Metabolism Genes in the Rhizosphere of Wheat

Despite the well-documented role of biochar in promoting soil quality and crop productivity, the underlying biological mechanisms remain poorly understood. Here, we explored the effects of straw biochar on soil microbiome in the rhizosphere from wheat using metagenomic sequencing. Our results showed that straw return decreased the yields of wheat, while the straw biochar return increased the wheat yields. Further, both the richness and community composition confirmed different effects of the straw return and straw biochar return. The straw biochar return also resulted in greater rhizosphere effects from wheat, represented by resource availability, including soil organic carbon, soil total nitrogen, available phosphorus, and available potassium. The rhizosphere effects from wheat, represented by microbial metabolism genes involved in carbon, nitrogen, phosphorus, and potassium cycling, however, were decreased by straw biochar returning. In addition, the rhizosphere effects from nitrogen content and the nitrogen cycling genes showed negative relationships with wheat yields. Together, these results revealed that straw biochar enhanced soil resource availability but suppressed microbial metabolism genes in the rhizosphere from wheat, supporting the idea that straw biochar serves as a nutrient pool for crops.

Expanding the Diversity of Accumulibacter with a Novel Type and Deciphering the Transcriptional and Morphological Features among Co-Occurring Strains

"Candidatus Accumulibacter" is the major polyphosphate-accumulating organism (PAO) in global wastewater treatment systems, and its phylogenetic and functional diversity have expanded in recent years. In addition to the widely recognized type I and II sublineages, we discovered a novel type enriched in laboratory bioreactors. Core gene and machine learning-based gene feature profiling supported the assertion that type III "Ca. Accumulibacter" is a potential PAO with the unique function of using dimethyl sulfoxide as an electron acceptor. Based on the correlation between ppk1 and genome similarity, the species-level richness of Accumulibacter was estimated to be over 100, suggesting that the currently recognized species are only the tip of the iceberg. Meanwhile, the interstrain transcriptional and morphological features of multiple "Ca. Accumulibacter" strains co-occurring in a bioreactor were investigated. Metatranscriptomics of seven co-occurring strains indicated that the expression level and interphasic dynamics of PAO phenotype-related genes had minimal correlation with their phylogeny. In particular, the expression of denitrifying and polyphosphate (poly-P) metabolism genes exhibited higher interstrain and interphasic divergence than expression of glycogen and polyhydroxyalkanoate metabolic genes. A strategy of cloning rRNA genes from different strains based on similar genomic synteny was successfully applied to differentiate their morphology via fluorescence in situ hybridization. Our study further expands the phylogenetic and functional diversity of "Ca. Accumulibacter" and proposes that deciphering the function and capability of certain "Ca. Accumulibacter" should be tailored to the environment and population in question. IMPORTANCE In the last 2 decades, "Ca. Accumulibacter" has garnered significant attention as the core functional but uncultured taxon for enhanced biological phosphorus removal due to its phylogenetic and functional diversity and intragenus niche differentiation. Since 2002, it has been widely known that this genus has two sublineages (type I and II). However, in this study, a metagenomic approach led to the discovery of a novel type (type III) with proposed novel functional features. By comparing the average nucleotide identity of "Ca. Accumulibacter" genomes and the similarity of ppk1, a phylogenetic biomarker largely deposited in databases, the global species-level richness of "Ca. Accumulibacter" was estimated for the first time to be over 100. Furthermore, we observed the co-occurrence of multiple "Ca. Accumulibacter" strains in a single bioreactor and found the simultaneous transcriptional divergence of these strains intriguing with regard to their niche differentiation within a single community. Our results indicated a decoupling feature between transcriptional pattern and phylogeny for co-occurring strains.