Charting the complexity of the activated sludge microbiome through a hybrid sequencing strategy

Stochastic and deterministic mechanisms jointly drive the assembly of microbial communities in cold-rolling wastewater across China

Microorganisms play a pivotal role in industrial wastewater treatment, serving as a critical barrier to water purification and safeguarding human and environmental health. Despite their importance, the biogeographic distribution and assembly mechanisms of microbial communities in cold-rolling wastewater treatment systems remain poorly understood. This study analyzed 101 microbial samples from nine regions using high-throughput sequencing, revealing rich microbial diversity and distinct regional aggregation patterns. Random forest analysis identified key biomarkers, often low-abundance species, while a unique core microbial community was strongly correlated with pollutant removal efficiencies, including chemical oxygen demand (COD), total organic carbon (TOC), and total nitrogen (TN). Neutral community model analysis demonstrated that microbial community assembly is driven by both stochastic and deterministic processes. Co-occurrence network analysis further highlighted o__1-20 and g__Ellin6067 as pivotal taxa influencing community structure. Among environmental factors, nitrite nitrogen (NO₂-N) and COD were identified as critical drivers of community assembly. This study provides the first comprehensive characterization of microbial biogeographic patterns in cold-rolling wastewater treatment plants across China. The findings deepen our understanding of microbial diversity, distribution, and community dynamics in industrial wastewater systems, offering valuable insights for optimizing treatment processes.

Nanopore Environmental Analysis

As global pollution continues to escalate, timely and accurate monitoring is essential for guiding pollution governance and safeguarding public health. The increasing diversity of pollutants across environmental matrices poses a significant challenge for instrumental analysis methods, which often require labor-intensive and time-consuming sample pretreatment. Nanopore technology, an emerging single-molecule technique, presents a promising solution by enabling the rapid identification of multiple targets within complex mixtures with minimal sample preparation. A wide range of pollutants have been characterized using natural biological nanopores or artificial solid-state nanopores, and their distinct advantages include simple sample preparation, high sensitivity, and rapid onsite analysis. In particular, long-read nanopore sequencing has led to dramatic improvements in the analyses of environmental microbial communities, allows species-level taxonomic assignment using amplicon sequencing, and simplifies the assembly of metagenomes. In this Perspective, we review the latest advancements in analyzing chemical and biological pollutants through nanopore sensing and sequencing techniques. We also explore the challenges that remain in this rapidly evolving field and provide an outlook on the potential for nanopore environmental analysis to transform pollution monitoring, risk assessment, and public health protection.

How small is big enough? Big data-driven machine learning predictions for a full-scale wastewater treatment plant

Wastewater treatment plants (WWTPs) generate vast amounts of water quality, operational, and biological data. The potential of these big data, particularly through machine learning (ML), to improve WWTP management is increasingly recognized. However, the costs associated with data collection and processing can rise sharply as datasets grow larger, and research on determining the optimal data volume for effective ML application remains limited. In this study, we comprehensively analyzed water quality, operational, and biological data collected from a full-scale WWTP over 970 days. Our results demonstrate that ML models can predict not only operational and water quality parameters (concentrations of dissolved oxygen and effluent chemical oxygen demand) but also the abundances of functional bacteria. Notably, we discovered that increasing data volume does not always improve model performance, and that data collection intervals do not need to be excessively small, as moderate intervals can still yield reliable predictions. These findings suggest that excessively large datasets may not be necessary for effective ML predictions in WWTPs. Overall, this study underscores the importance of optimizing dataset size to balance computation efficiency and prediction accuracy, providing valuable insights into data management strategies that can enhance the operational efficiency and sustainability of WWTPs.

Exploring phylogenetic diversity of antibiotic resistance genes in activated sludge: A host and genomic location perspective

Antibiotic resistance has emerged as a significant global public health issue. The environmental behaviors of antibiotic resistance genes (ARGs), such as their persistence and horizontal transfer, have been extensively investigated. However, the genetic diversity characteristics of ARGs remain underexplored, which limits a comprehensive analysis of their roles in the environment. In this study, we examined the genetic diversity of ARGs in activated sludge from 44 wastewater treatment plants in five countries. Most ARGs detected in activated sludge possessed multiple variants, with a median of 48. The number of variants of gd-ARGs varied among different resistance mechanisms and ARG types. The number of potential variants of ARGs was strongly correlated with host diversity. Pseudomonas spp. and Klebsiella pneumoniae, identified as pathogenic bacteria, harbored multiple ARGs and had the most variants. Most ARG subtypes on plasmids and chromosomes showed divergent evolution. Molecular docking of AdeH proteins revealed that genomic location affects tetracycline binding energy. The findings underscore the intricate interplay between genetic variation and environmental adaptation in ARGs, offering a novel perspective on the spread of antibiotic resistance.

Description of new diplonemids (Diplonemea, Euglenozoa) and their endosymbionts: Charting the morphological diversity of these poorly known heterotrophic flagellates

Diplonemids are a hyperdiverse group of flagellated protists, but with less than two dozen formally described representatives. Here, we describe four new species of cultured diplonemids, identified on the basis of their 18S rRNA sequences, light-, fluorescence-, scanning- and transmission electron microscopy. Three new species belong to the genus Rhynchopus (R. asiaticus sp.n., R. granulatus sp.n., and R. valaseki sp.n.), while the fourth species is an unusual representative of the genus Lacrimia (L. aflagellata sp.n.). The latter organism is the first diplonemid outside the genus Rhynchopus (as defined previously) to show a gliding trophic stage with flagellar stubs concealed inside the flagellar pocket and a highly motile dispersive swimming stage. Since this character is thus no longer a genus-specific apomorphy, we provide a taxonomic revision of the genus Rhynchopus with separation of the new genus Natarhynchopus gen. n. We also identify bacterial endosymbionts of L. aflagellata and R. asiaticus as Ca. Syngnamydia medusae (Chlamydiales, Simkaniaceae) and Ca. Cytomitobacter rhynchopi sp. n. (Alphaproteobacteria, Holosporaceae), respectively, and discuss their potential functions. This is the first report of a chlamydial symbiont within a diplonemid host. We also propose that diplonemids may serve as vectors for chlamydial pathogens of marine fish.

Soil warming increases the active antibiotic resistome in the gut of invasive giant African snails

BACKGROUND

Global warming is redrawing the map for invasive species, spotlighting the globally harmful giant African snail as a major ecological disruptor and public health threat. Known for harboring extensive antibiotic resistance genes (ARGs) and human pathogens, it remains uncertain whether global warming exacerbates these associated health risks.

METHODS

We use phenotype-based single-cell Raman with D2O labeling (Raman-D2O) and genotype-based metagenomic sequencing to investigate whether soil warming increases active antibiotic-resistant bacteria (ARBs) in the gut microbiome of giant African snails.

RESULTS

We show a significant increase in beta-lactam phenotypic resistance of active ARBs with rising soil temperatures, mirrored by a surge in beta-lactamase genes such as SHV, TEM, OCH, OKP, and LEN subtypes. Through a correlation analysis between the abundance of phenotypically active ARBs and genotypically ARG-carrying gut microbes, we identify species that contribute to the increased activity of antibiotic resistome under soil warming. Among 299 high-quality ARG-carrying metagenome-assembled genomes (MAGs), we further revealed that the soil warming enhances the abundance of "supercarriers" including human pathogens with multiple ARGs and virulence factors. Furthermore, we identified elevated biosynthetic gene clusters (BGCs) within these ARG-carrying MAGs, with a third encoding at least one BGC. This suggests a link between active ARBs and secondary metabolism, enhancing the environmental adaptability and competitive advantage of these organisms in warmer environments.

CONCLUSIONS

The study underscores the complex interactions between soil warming and antibiotic resistance in the gut microbiome of the giant African snail, highlighting a potential escalation in environmental health risks due to global warming. These findings emphasize the urgent need for integrated environmental and health strategies to manage the rising threat of antibiotic resistance in the context of global climate change. Video Abstract.

Longitudinal dynamics and cross-domain interactions of eukaryotic populations in wastewater treatment plants

Activated sludge is a large reservoir of novel microorganisms and microbial genetic diversity. While much attention has been given to the profile and functions of prokaryotes, the eukaryotic diversity remains largely unexplored. In this study, we analysed longitudinal activated sludge samples spanning 13 years from the largest secondary wastewater treatment plants in Hong Kong, unveiling a wealth of eukaryotic taxa and 681 856 non-redundant protein-coding genes, the majority (416 044) of which appeared novel. Ciliophora was the most dominant phylum with a significant increase after a transient intervention (bleaching event). Our metagenomic analysis reveals close linkage and covariation of eukaryotes, prokaryotes, and prokaryotic viruses (phages), indicating common responses to environmental changes such as transient intervention and intermittent fluctuations. Furthermore, high-resolution cross-domain relationships were interpreted by S-map, demonstrating a predatory role of Arthropoda, Ascomycota, Mucoromycota, and Rotifera. This high-resolution profile of microbial dynamics expands our knowledge on yet-to-be-cultured populations and their cross-domain interactions and highlights the ecological importance of eukaryotes in the activated sludge ecosystem.

High-throughput single-cell sequencing of activated sludge microbiome

Wastewater treatment plants (WWTPs) represent one of biotechnology's largest and most critical applications, playing a pivotal role in environmental protection and public health. In WWTPs, activated sludge (AS) plays a major role in removing contaminants and pathogens from wastewater. While metagenomics has advanced our understanding of microbial communities, it still faces challenges in revealing the genomic heterogeneity of cells, uncovering the microbial dark matter, and establishing precise links between genetic elements and their host cells as a bulk method. These issues could be largely resolved by single-cell sequencing, which can offer unprecedented resolution to show the unique genetic information. Here we show the high-throughput single-cell sequencing to the AS microbiome. The single-amplified genomes (SAGs) of 15,110 individual cells were clustered into 2,454 SAG bins. We find that 27.5% of the genomes in the AS microbial community represent potential novel species, highlighting the presence of microbial dark matter. Furthermore, we identified 1,137 antibiotic resistance genes (ARGs), 10,450 plasmid fragments, and 1,343 phage contigs, with shared plasmid and phage groups broadly distributed among hosts, indicating a high frequency of horizontal gene transfer (HGT) within the AS microbiome. Complementary analysis using 1,529 metagenome-assembled genomes from the AS samples allowed for the taxonomic classification of 98 SAG bins, which were previously unclassified. Our study establishes the feasibility of single-cell sequencing in characterizing the AS microbiome, providing novel insights into its ecological dynamics, and deepening our understanding of HGT processes, particularly those involving ARGs. Additionally, this valuable tool could monitor the distribution, spread, and pathogenic hosts of ARGs both within AS environments and between AS and other environments, which will ultimately contribute to developing a health risk evaluation system for diverse environments within a One Health framework.

Multi-omics illuminates the functional significance of previously unknown species in a full-scale landfill leachate treatment plant

Landfill leachate treatment plants (LLTPs) harbor a vast reservoir of uncultured microbes, yet limited studies have systematically unraveled their functional potentials within LLTPs. Combining 36 metagenomic and 18 metatranscriptomic datasets from a full-scale LLTP, we unveiled a double-edged sword role of unknown species in leachate biotreatment and environmental implication. We identified 655 species-level genome bins (SGBs) spanning 47 bacterial and 3 archaeal phyla, with 75.9 % unassigned to any known species. Over 90 % of up-regulated functional genes in biotreatment units, compared to the leachate influent, were carried by unknown species and actively participated in carbon, nitrogen, and sulfur cycles. Approximately 79 % of the 37,366 carbohydrate active enzymes (CAZymes), with ∼90 % novelty and high expression, were encoded by unknown species, exhibiting great potential in biodegrading carbohydrate compounds linked to human meat-rich diets. Unknown species offered a valuable genetic resource of thousands of versatile, abundant, and actively expressed metabolic gene clusters (MGCs) and biosynthetic gene clusters (BGCs) for enhancing leachate treatment. However, unknown species may contribute to the emission of hazardous N2O/H2S and represented significant reservoirs for antibiotic-resistant pathogens that posed environmental safety risks. This study highlighted the significance of considering both positive and adverse effects of LLTP microbes to optimize LLTP performance.

Advances in lasso peptide discovery, biosynthesis, and function

Lasso peptides are a large and sequence-diverse class of ribosomally synthesized and post-translationally modified peptide (RiPP) natural products characterized by their slip knot-like shape. These unique, highly stable peptides are produced by bacteria for various purposes. Their stability and sequence diversity make them a potentially useful scaffold for biomedically relevant folded peptides. However, many questions remain about lasso peptide biosynthesis, ecological function, and diversification potential for biomedical and agricultural applications. This review discusses new insights and open questions about lasso peptide biosynthesis and biological function. The role that genome mining has played in the development of new methodologies for discovering and diversifying lasso peptides is also discussed.

The increase of particle size shifts the biogeochemical cycle functions of mineral-associated microorganisms and weakens the mineral-associated organic carbon sink in mangrove soils

Mineral-organic matter-microbe interactions greatly impact the biogeochemical processes and biodiversity in soils. An increasing trend of particle size (PS) in mangrove soils has been observed because of the relative sea level rise. However, the impacts of PS increase on the microbial biogeochemical functions and carbon sink in the mineral-associated microcosms are exceedingly nebulous. This work showed a remarkable difference in the communities of mineral-associated microorganisms (MMOs) in various PS fractions. Heavy metal contents and urease activity were the factors that mostly driven the MMO community variation in different PS fractions. Large PS fraction attenuated the stability of MMO communities according to the co-occurrence network characteristics. The PS increase significantly (P < 0.05) lowered the gene abundances for carbon input (e.g., carbon fixation) and raised the gene abundances for carbon loss (e.g., aerobic respiration). Combined with the significant decrease of mineral-associated organic carbon (MOC) in large PS fraction (P < 0.05), this work inferred that the PS increase could weaken the MOC sink partially due to the MMO function shift for carbon cycle. The current work indicated unhealthy changes of MMO communities and MOC storage in mangrove soils, and PS was of significance as an indicator for predicting the carbon sink function, especially for the stable form, such as MOC, in the soils of mangroves under the ecological background of climate migration.

IMPORTANCE

Carbon with stable forms, such as mineral-associated organic carbon (MOC), is crucial for the sink capabilities in mangrove soils, and mineral-associated microorganisms (MMOs) are important players for the formation and metabolism of MOC. Therefore, the future successions of the MMO functions and MOC contents under the background of climate change are of value for a deeper understanding of mangrove ecology. The relative sea level rise caused by the global warming results in the increase of mangrove soil particle size (PS), which provides distinct microcosms for MMOs and MOC. However, the responses of MMO functions and MOC content to the PS increase of mangrove soils are unknown. The current study aims to reveal the succession regulations of MMO functions and their potential ecological impacts for the storages of MOC in different PS fractions, therefore widening our knowledge of future function migration and promoting the research development of mangrove.

Metagenomic insights into ecological risk of antibiotic resistome and mobilome in riverine plastisphere under impact of urbanization

Microplastics (MPs) are of increasing concern due to their role as reservoirs for antibiotic resistance genes (ARGs) and pathogens. To date, few studies have explored the influence of anthropogenic activities on ARGs and mobile genetic elements (MGEs) within various riverine MPs, in comparison to their natural counterparts. Here an in-situ incubation was conducted along heavily anthropogenically-impacted Houxi River to characterize the geographical pattern of antibiotic resistome, mobilome and pathogens inhabiting MPs- and leaf-biofilms. The metagenomics result showed a clear urbanization-driven profile in the distribution of ARGs, MGEs and pathogens, with their abundances sharply increasing 4.77 to 19.90 times from sparsely to densely populated regions. The significant correlation between human fecal marker crAssphage and ARG (R2 = 0.67, P=0.003) indicated the influence of anthropogenic activity on ARG proliferation in plastisphere and natural leaf surfaces. And mantel tests and random forest analysis revealed the impact of 17 socio-environmental factors, e.g., population density, antibiotic concentrations, and pore volume of materials, on the dissemination of ARGs. Partial least squares-path modeling further unveiled that intensifying human activities not only directly boosted ARGs abundance but also exerted a comparable indirect impact on ARGs propagation. Furthermore, the polyvinylchloride plastisphere created a pathogen-friendly habitat, harboring higher abundances of ARGs and MGEs, while polylactic acid are not likely to serve as vectors for pathogens in river, with a lower resistome risk score than that in leaf-biofilms. This study highlights the diverse ecological risks associated with the dissemination of ARGs and pathogens in varied MPs, offering insights for the policymaking of usage and control of plastics within urbanization.

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.

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.

Dirammox-dominated microbial community for biological nitrogen removal from wastewater

Direct ammonia oxidation (Dirammox) might be of great significance to advance the innovation of biological nitrogen removal process in wastewater treatment systems. However, it remains unknown whether Dirammox bacteria can be selectively enriched in activated sludge. In this study, a lab-scale bioreactor was established and operated for 2 months to treat synthetic wastewater with hydroxylamine as a selection pressure. Three Dirammox strains (Alcaligenes aquatilis SDU_AA1, Alcaligenes aquatilis SDU_AA2, and Alcaligenes sp. SDU_A2) were isolated from the activated sludge, and their capability to perform Dirammox process was confirmed. Although these three Dirammox bacteria were undetectable in the seed sludge (0%), their relative abundances rapidly increased after a month of operation, reaching 12.65%, 0.69%, and 0.69% for SDU_A2, SDU_AA1, and SDU_AA2, respectively. Among them, the most dominant Dirammox (SDU_A2) exhibited higher nitrogen removal rate (32.35%) than the other two strains (13.57% of SDU_AA1 and 14.52% of SDU_AA2). Comparative genomic analysis demonstrated that the most dominant Dirammox bacterium (SDU_A2) possesses fewer complete metabolic modules compared to the other two less abundant Alcaligenes strains. Our findings expanded the understanding of the application of Dirammox bacteria as key functional microorganisms in a novel biological nitrogen and carbon removal process if they could be well stabilized. KEY POINTS: • Dirammox-dominated microbial community was enriched in activated sludge bioreactor. • The addition of hydroxylamine played a role in Dirammox enrichment. • Three Dirammox bacterial strains, including one novel species, were isolated.

Corallococcus caeni sp. nov., a novel myxobacterium isolated from activated sludge

Two myxobacterial strains (KH5-1T and NO1) were isolated from the activated sludge tanks treating municipal sewage wastewater in Japan. These strains were recognised as myxobacteria based on their phenotypic characteristics of swarming colonies and fruiting bodies. Phylogenetic analyses using the 16S rRNA gene revealed that strains KH5-1T and NO1 were affiliated with the genus Corallococcus, with the closest neighbours being Corallococcus exercitus AB043AT (99.77% and 99.84%, respectively). Genome comparisons using orthologous average nucleotide identity (orthoANI) and digital DNA-DNA hybridisation similarity (dDDH) with strains KH5-1T and NO1 and their phylogenetically close relatives in Corallococcus spp. were below the thresholds. The major cellular fatty acids of strains KH5-1T and NO1 were iso-C15:0 (31.9%, 30.0%), summed feature 3 (comprising C16:1ω7c and/or C16:1ω6c) (20.2%, 17.7%), and iso-C17:0 (12.1%, 14.8%), and the major respiratory quinone was found to be menaquinone (MK)-8. Based on the phenotypic, chemotaxonomic, and phylogenetic evidence, strains KH5-1T and NO1 represent a new species in the genus Corallococcus, for which the proposed name is Corallococcus caeni sp. nov. The type strain is KH5-1T (= NCIMB 15510T = JCM 36609T).

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.

Comprehensive analysis of genomic variation, pan-genome and biosynthetic potential of Corynebacterium glutamicum strains

Corynebacterium glutamicum is a non-pathogenic species of the Corynebacteriaceae family. It has been broadly used in industrial biotechnology for the production of valuable products. Though it is widely accepted at the industrial level, knowledge about the genomic diversity of the strains is limited. Here, we investigated the comparative genomic features of the strains and pan-genomic characteristics. We also observed phylogenetic relationships among the strains based on average nucleotide identity (ANI). We found diversity between strains at the genomic and pan-genomic levels. Less than one-third of the C. glutamicum pan-genome consists of core genes and soft-core genes. Whereas, a large number of strain-specific genes covered about half of the total pan-genome. Besides, C. glutamicum pan-genome is open and expanding, which indicates the possible addition of new gene families to the pan-genome. We also investigated the distribution of biosynthetic gene clusters (BGCs) among the strains. We discovered slight variations of BGCs at the strain level. Several BGCs with the potential to express novel bioactive secondary metabolites have been identified. Therefore, by utilizing the characteristic advantages of C. glutamicum, different strains can be potential applicants for natural drug discovery.

Critical steps in an environmental metaproteomics workflow

Environmental metaproteomics is a rapidly advancing field that provides insights into the structure, dynamics, and metabolic activity of microbial communities. As the field is still maturing, it lacks consistent workflows, making it challenging for non-expert researchers to navigate. This review aims to introduce the workflow of environmental metaproteomics. It outlines the standard practices for sample collection, processing, and analysis, and offers strategies to overcome the unique challenges presented by common environmental matrices such as soil, freshwater, marine environments, biofilms, sludge, and symbionts. The review also highlights the bottlenecks in data analysis that are specific to metaproteomics samples and provides suggestions for researchers to obtain high-quality datasets. It includes recent benchmarking studies and descriptions of software packages specifically built for metaproteomics analysis. The article is written without assuming the reader's familiarity with single-organism proteomic workflows, making it accessible to those new to proteomics or mass spectrometry in general. This primer for environmental metaproteomics aims to improve accessibility to this exciting technology and empower researchers to tackle challenging and ambitious research questions. While it is primarily a resource for those new to the field, it should also be useful for established researchers looking to streamline or troubleshoot their metaproteomics experiments.

Effects of operational parameters on bacterial communities in Hong Kong and global wastewater treatment plants

Wastewater treatment plants (WWTPs) are indispensable biotechnology facilities for modern cities and play an essential role in modern urban infrastructure by employing microorganisms to remove pollutants in wastewater, thus protecting public health and the environment. This study conducted a 13-month bacterial community survey of six full-scale WWTPs in Hong Kong with samples of influent, activated sludge (AS), and effluent to explore their synchronism and asynchronism of bacterial community. Besides, we compared AS results of six Hong Kong WWTPs with data from 1,186 AS amplicon data in 269 global WWTPs and a 9-year metagenomic sequencing survey of a Hong Kong WWTP. Our results showed the compositions of bacterial communities varied and the bacterial community structure of AS had obvious differences across Hong Kong WWTPs. The co-occurrence analysis identified 40 pairs of relationships that existed among Hong Kong WWTPs to show solid associations between two species and stochastic processes took large proportions for the bacterial community assembly of six WWTPs. The abundance and distribution of the functional bacteria in worldwide and Hong Kong WWTPs were examined and compared, and we found that ammonia-oxidizing bacteria had more diversity than nitrite-oxidizing bacteria. Besides, Hong Kong WWTPs could make great contributions to the genome mining of microbial dark matter in the global "wanted list." Operational parameters had important effects on OTUs' abundance, such as the temperature to the genera of Tetrasphaera, Gordonia and Nitrospira. All these results obtained from this study can deepen our understanding of the microbial ecology in WWTPs and provide foundations for further studies.

IMPORTANCE

Wastewater treatment plants (WWTPs) are an indispensable component of modern cities, as they can remove pollutants in wastewater to prevent anthropogenic activities. Activated sludge (AS) is a fundamental wastewater treatment process and it harbors a highly complex microbial community that forms the main components and contains functional groups. Unveiling "who is there" is a long-term goal of the research on AS microbiology. High-throughput sequencing provides insights into the inventory diversity of microbial communities to an unprecedented level of detail. At present, the analysis of communities in WWTPs usually comes from a specific WWTP and lacks comparisons and verification among different WWTPs. The wide-scale and long-term sampling project and research in this study could help us evaluate the AS community more accurately to find the similarities and different results for different WWTPs in Hong Kong and other regions of the world.

BASALT refines binning from metagenomic data and increases resolution of genome-resolved metagenomic analysis

Metagenomic binning is an essential technique for genome-resolved characterization of uncultured microorganisms in various ecosystems but hampered by the low efficiency of binning tools in adequately recovering metagenome-assembled genomes (MAGs). Here, we introduce BASALT (Binning Across a Series of Assemblies Toolkit) for binning and refinement of short- and long-read sequencing data. BASALT employs multiple binners with multiple thresholds to produce initial bins, then utilizes neural networks to identify core sequences to remove redundant bins and refine non-redundant bins. Using the same assemblies generated from Critical Assessment of Metagenome Interpretation (CAMI) datasets, BASALT produces up to twice as many MAGs as VAMB, DASTool, or metaWRAP. Processing assemblies from a lake sediment dataset, BASALT produces ~30% more MAGs than metaWRAP, including 21 unique class-level prokaryotic lineages. Functional annotations reveal that BASALT can retrieve 47.6% more non-redundant opening-reading frames than metaWRAP. These results highlight the robust handling of metagenomic sequencing data of BASALT.

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.

Draft genome sequences of Corallococcus strains KH5-1 and NO1, isolated from activated sludge

Myxobacteria are known as prolific producers of secondary metabolites with a unique and wide spectrum of bioactivities. Here, we report draft genome sequences of KH5-1 and NO1, myxobacteria isolated from activated sludge, which consist of 9.89 and 9.86 Mb, both of which have G + C contents of 70.7%.

Genome-Resolved Metagenomics and Denitrifying Strain Isolation Reveal New Insights into Microbial Denitrification in the Deep Vadose Zone

The heavy use of nitrogen fertilizer in intensive agricultural areas often leads to nitrate accumulation in subsurface soil and nitrate contamination in groundwater, which poses a serious risk to public health. Denitrifying microorganisms in the subsoil convert nitrate to gaseous forms of nitrogen, thereby mitigating the leaching of nitrate into groundwater. Here, we investigated denitrifying microorganisms in the deep vadose zone of a typical intensive agricultural area in China through microcosm enrichment, genome-resolved metagenomic analysis, and denitrifying bacteria isolation. A total of 1000 metagenome-assembled genomes (MAGs) were reconstructed, resulting in 98 high-quality, dereplicated MAGs that contained denitrification genes. Among them, 32 MAGs could not be taxonomically classified at the genus or species level, indicating that a broader spectrum of taxonomic groups is involved in subsoil denitrification than previously recognized. A denitrifier isolate library was constructed by using a strategy combining high-throughput and conventional cultivation techniques. Assessment of the denitrification characteristics of both the MAGs and isolates demonstrated the dominance of truncated denitrification. Functional screening revealed the highest denitrification activity in two complete denitrifiers belonging to the genus Pseudomonas. These findings greatly expand the current knowledge of the composition and function of denitrifying microorganisms in subsoils. The constructed isolate library provided the first pool of subsoil-denitrifying microorganisms that could facilitate the development of microbe-based technologies for nitrate attenuation in groundwater.

Strain-level diversity in sulfonamide biodegradation: adaptation of Paenarthrobacter to sulfonamides

The widespread occurrence of sulfonamides raises significant concerns about the evolution and spread of antibiotic resistance genes. Biodegradation represents not only a resistance mechanism but also a clean-up strategy. Meanwhile, dynamic and diverse environments could influence the cellular function of individual sulfonamide-degrading strains. Here, we present Paenarthrobacter from different origins that demonstrated diverse growth patterns and sulfonamide-degrading abilities. Generally, the degradation performance was largely associated with the number of sadA gene copies and also relied on its genotype. Based on the survey of sad genes in the public database, an independent mobilization of transposon-borne genes between chromosome and plasmid was observed. Insertions of multiple sadA genes could greatly enhance sulfonamide-degrading performance. Moreover, the sad gene cluster and sadA transposable element showed phylogenetic conservation currently, being identified only in two genera of Paenarthrobacter (Micrococcaceae) and Microbacterium (Microbacteriaceae). Meanwhile, Paenarthrobacter exhibited a high capacity for genome editing to adapt to the specific environmental niche, opening up new opportunities for bioremediation applications.

Comprehensive Insights into Potential Metabolic Functions of Myxococcota in Activated Sludge Systems

Myxobacteria, belonging to the phylum Myxococcota, are ubiquitous in soil, marine, and other environments. A recent metagenomic sequencing ana-lysis showed that Myxococcota are predominant in activated sludge systems; however, their metabolic traits remain unclear. In the present study, we exami-ned the potential biological functions of 46 metagenomic bins of Myxococcota reconstructed from activated sludge samples from four municipal sewage treatment plants. The results obtained showed that most Myxococcota bins had an almost complete set of genes associated with glycolysis and the TCA cycle. The Palsa-1104 and Polyangiales bins contained the glycoside hydrolase GH5 and peptidase M23, which are presumably involved in lysis of the cell wall and cellular cytoplasm, suggesting that some Myxococcota from activated sludge prey on other microorganisms. The cell contact-dependent predatory functions of Myxococcus xanthus are conserved in the family Myxococcaceae, but not in other families. Two bins belonging to Palsa-1104 had phototrophic gene clusters, indicating the potential for heterotrophic and autotrophic metabolism by these microbes. In assessments of the social behavior of Myxococcota in activated sludge, the FruA gene and C-signal gene, which are involved in the regulation of fruiting body formation, were lacking in Myxococcota bins, suggesting their inability to form fruiting bodies. In addition, multiple bins of Myxococcota had novel secondary metabolite biosynthesis gene clusters that may be used for the predation of other bacteria in activated sludge. Our metagenome-based ana-lyses provide novel insights into the microbial interactions associated with Myxococcota in activated sludge ecosystems.

A small step to discover candidate biological control agents from preexisting bioresources by using novel nonribosomal peptide synthetases hidden in activated sludge metagenomes

Biological control agents (BCAs), beneficial organisms that reduce the incidence or severity of plant disease, have been expected to be alternatives to replace chemical pesticides worldwide. To date, BCAs have been screened by culture-dependent methods from various environments. However, previously unknown BCA candidates may be buried and overlooked because this approach preferentially selects only easy-to-culture microbial lineages. To overcome this limitation, as a small-scale test case, we attempted to explore novel BCA candidates by employing the shotgun metagenomic information of the activated sludge (AS) microbiome, which is thought to contain unutilized biological resources. We first performed genome-resolved metagenomics for AS taken from a municipal sewage treatment plant and obtained 97 nonribosomal peptide synthetase (NRPS)/polyketide synthase (PKS)-related gene sequences from 43 metagenomic assembled bins, most of which were assigned to the phyla Proteobacteria and Myxococcota. Furthermore, these NRPS/PKS-related genes are predicted to be novel because they were genetically dissimilar to known NRPS/PKS gene clusters. Of these, the condensation domain of the syringomycin-related NRPS gene cluster was detected in Rhodoferax- and Rhodocyclaceae-related bins, and its homolog was found in previously reported AS metagenomes as well as the genomes of three strains available from the microbial culture collections, implying their potential BCA ability. Then, we tested the antimicrobial activity of these strains against phytopathogenic fungi to investigate the potential ability of BCA by in vitro cultivation and successfully confirmed the actual antifungal activity of three strains harboring a possibly novel NRPS gene cluster. Our findings provide a possible strategy for discovering novel BCAs buried in the environment using genome-resolved metagenomics.

A genomic catalogue of soil microbiomes boosts mining of biodiversity and genetic resources

Soil harbors a vast expanse of unidentified microbes, termed as microbial dark matter, presenting an untapped reservo)ir of microbial biodiversity and genetic resources, but has yet to be fully explored. In this study, we conduct a large-scale excavation of soil microbial dark matter by reconstructing 40,039 metagenome-assembled genome bins (the SMAG catalogue) from 3304 soil metagenomes. We identify 16,530 of 21,077 species-level genome bins (SGBs) as unknown SGBs (uSGBs), which expand archaeal and bacterial diversity across the tree of life. We also illustrate the pivotal role of uSGBs in augmenting soil microbiome's functional landscape and intra-species genome diversity, providing large proportions of the 43,169 biosynthetic gene clusters and 8545 CRISPR-Cas genes. Additionally, we determine that uSGBs contributed 84.6% of previously unexplored viral-host associations from the SMAG catalogue. The SMAG catalogue provides an useful genomic resource for further studies investigating soil microbial biodiversity and genetic resources.

MetaSVs: A pipeline combining long and short reads for analysis and visualization of structural variants in metagenomes

Structural variants (SVs, including large-scale insertions, deletions, inversions, and translocations) significantly impact the functions of genes in the microbial genome, and SVs in the microbiome are associated with diverse biological processes and human diseases. With the advancements in sequencing and bioinformatics technologies, increasingly, sequencing data and analysis tools are already being extensively utilized for microbiome SV analyses, leading to a higher demand for more dedicated SV analysis workflows. Moreover, due to the unique detection biases of various sequencing technologies, including short-read sequencing (such as Illumina platforms) and long-read sequencing (e.g., Oxford Nanopore and PacBio), SV discovery based on multiple platforms is necessary to comprehensively identify the wide variety of SVs. Here, we establish an integrated pipeline MetaSVs combining Nanopore long reads and Illumina short reads to analyze SVs in the microbial genomes from gut microbiome and further identify differential SVs that can be reflective of metabolic differences. Our pipeline provides researchers easy access to SVs and relevant metabolites in the microbial genomes without the requirement of specific technical expertise, which is particularly useful to researchers interested in metagenomic SVs but lacking sophisticated bioinformatic knowledge.

Establishing reference material for the quest towards standardization in environmental microbial metagenomic studies

Breakthroughs in DNA-based technologies, especially in metagenomic sequencing, have drastically enhanced researchers' ability to explore environmental microbiome and the associated interplays within. However, as new methodologies are being actively developed for improvements in different aspects, metagenomic workflows become diversified and heterogeneous. Through a single-variable control approach, we quantified the microbial profiling variations arising from 6 common technical variables associated with metagenomic workflows for both simple and complex samples. The incurred variations were constantly the lowest in replicates of DNA isolation and DNA sequencing library construction. Different DNA extraction kits often caused the highest variation among all the tested variables. Additionally, sequencing run batch was an important source of variability for targeted platforms. As such, the development of an environmental reference material for complex environmental samples could be beneficial in benchmarking accrued non-biological variability within and between protocols and insuring reliable and reproducible sequencing outputs immediately upstream of bioinformatic analysis. To develop an environment reference material, sequencing of a well-homogenized environmental sample composed of activated sludge was performed using different pre-analytical assays in replications. In parallel, a certified mock community was processed and sequenced. Assays were ranked based on the reconstruction of the theoretical mock community profile. The reproducibility of the best-performing assay and the microbial profile of the reference material were further ascertained. We propose the adoption of our complex environmental reference material, which could reflect the degree of diversity in environmental microbiome studies, to facilitate accurate, reproducible, and comparable environmental metagenomics-based studies.

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.

The microbial dark matter and "wanted list" in worldwide wastewater treatment plants

BACKGROUND

Wastewater treatment plants (WWTPs) are one of the largest biotechnology applications in the world and are of critical importance to modern urban societies. An accurate evaluation of the microbial dark matter (MDM, microorganisms whose genomes remain uncharacterized) proportions in WWTPs is of great value, while there is no such research yet. This study conducted a global meta-analysis of MDM in WWTPs with 317,542 prokaryotic genomes from the Genome Taxonomy Database and proposed a "wanted list" for priority targets in further investigations of activated sludge.

RESULTS

Compared with the Earth Microbiome Project data, WWTPs had relatively lower genome-sequenced proportions of prokaryotes than other ecosystems, such as the animal related environments. Analysis showed that the median proportions of the genome-sequenced cells and taxa (100% identity and 100% coverage in 16S rRNA gene region) in WWTPs reached 56.3% and 34.5% for activated sludge, 48.6% and 28.5% for aerobic biofilm, and 48.3% and 28.5% for anaerobic digestion sludge, respectively. This result meant MDM had high proportions in WWTPs. Besides, all of the samples were occupied by a few predominant taxa, and the majority of the sequenced genomes were from pure cultures. The global-scale "wanted list" for activated sludge contained four phyla that have few representatives and 71 operational taxonomic units with the majority of them having no genome or isolate yet. Finally, several genome mining methods were verified to successfully recover genomes from activated sludge such as hybrid assembly of the second- and third-generation sequencing.

CONCLUSIONS

This work elucidated the proportion of MDM in WWTPs, defined the "wanted list" of activated sludge for future investigations, and certified potential genome recovery methods. The proposed methodology of this study can be applied to other ecosystems and improve understanding of ecosystem structure across diverse habitats. Video Abstract.

Genome-centric metagenomics reveals the host-driven dynamics and ecological role of CPR bacteria in an activated sludge system

BACKGROUND

Candidate phyla radiation (CPR) constitutes highly diverse bacteria with small cell sizes and are likely obligate intracellular symbionts. Given their distribution and complex associations with bacterial hosts, genetic and biological features of CPR bacteria in low-nutrient environments have received increasing attention. However, CPR bacteria in wastewater treatment systems remain poorly understood. We utilized genome-centric metagenomics to answer how CPR communities shift over 11 years and what kind of ecological roles they act in an activated sludge system.

RESULTS

We found that approximately 9% (135) of the 1,526 non-redundant bacterial and archaeal metagenome-assembled genomes were affiliated with CPR. CPR bacteria were consistently abundant with a relative abundance of up to 7.5% in the studied activated sludge system. The observed striking fluctuations in CPR community compositions and the limited metabolic and biosynthetic capabilities in CPR bacteria collectively revealed the nature that CPR dynamics may be directly determined by the available hosts. Similarity-based network analysis further confirmed the broad bacterial hosts of CPR lineages. The proteome contents of activated sludge-associated CPR had a higher similarity to those of environmental-associated CPR than to those of human-associated ones. Comparative genomic analysis observed significant enrichment of genes for oxygen stress resistance in activated sludge-associated CPR bacteria. Furthermore, genes for carbon cycling and horizontal gene transfer were extensively identified in activated sludge-associated CPR genomes.

CONCLUSIONS

These findings highlight the presence of specific host interactions among CPR lineages in activated sludge systems. Despite the lack of key metabolic pathways, these small, yet abundant bacteria may have significant involvements in biogeochemical cycling and bacterial evolution in activated sludge systems. Video Abstract.

Metagenomic Analysis of Anaerobic Microbial Communities Degrading Short-Chain Fatty Acids as Sole Carbon Sources

Analyzing microbial communities using metagenomes is a powerful approach to understand compositional structures and functional connections in anaerobic digestion (AD) microbiomes. Whereas short-read sequencing approaches based on the Illumina platform result in highly fragmented metagenomes, long-read sequencing leads to more contiguous assemblies. To evaluate the performance of a hybrid approach of these two sequencing approaches we compared the metagenome-assembled genomes (MAGs) resulting from five AD microbiome samples. The samples were taken from reactors fed with short-chain fatty acids at different feeding regimes (continuous and discontinuous) and organic loading rates (OLR). Methanothrix showed a high relative abundance at all feeding regimes but was strongly reduced in abundance at higher OLR, when Methanosarcina took over. The bacterial community composition differed strongly between reactors of different feeding regimes and OLRs. However, the functional potential was similar regardless of feeding regime and OLR. The hybrid sequencing approach using Nanopore long-reads and Illumina MiSeq reads improved assembly statistics, including an increase of the N50 value (on average from 32 to 1740 kbp) and an increased length of the longest contig (on average from 94 to 1898 kbp). The hybrid approach also led to a higher share of high-quality MAGs and generated five potentially circular genomes while none were generated using MiSeq-based contigs only. Finally, 27 hybrid MAGs were reconstructed of which 18 represent potentially new species-15 of them bacterial species. During pathway analysis, selected MAGs revealed similar gene patterns of butyrate degradation and might represent new butyrate-degrading bacteria. The demonstrated advantages of adding long reads to metagenomic analyses make the hybrid approach the preferable option when dealing with complex microbiomes.

Evaluation of multiplex nanopore sequencing for Salmonella serotype prediction and antimicrobial resistance gene and virulence gene detection

In a previous study, Multiplex-nanopore-sequencing based whole genome sequencing (WGS) allowed for accurate in silico serotype prediction of Salmonella within one day for five multiplexed isolates, using both SISTR and SeqSero2. Since only ten serotypes were tested in our previous study, the conclusions above were yet to be evaluated in a larger scale test. In the current study we evaluated this workflow with 69 Salmonella serotypes and also explored the feasibility of using multiplex-nanopore-sequencing based WGS for antimicrobial resistance gene (AMR) and virulence gene detection. We found that accurate in silico serotype prediction with nanopore-WGS data was achieved within about five hours of sequencing at a minimum of 30× Salmonella genome coverage, with SeqSero2 as the serotype prediction tool. For each tested isolate, small variations were observed between the AMR/virulence gene profiles from the Illumina and Nanopore sequencing platforms. Taking results generated using Illumina data as the benchmark, the average precision value per isolate was 0.99 for both AMR and virulence gene detection. We found that the resistance gene identifier - RGI identified AMR genes with nanopore data at a much lower accuracy compared to Abricate, possibly due to RGI's less stringent minimum similarity and coverage by default for database matching. This study is an evaluation of multiplex-nanopore-sequencing based WGS as a cost-efficient and rapid Salmonella classification method, and a starting point for future validation and verification of using it as a AMR/virulence gene profiling tool for the food industry. This study paves the way for the application of nanopore sequencing in surveillance, tracking, and risk assessment of Salmonella across the food supply chain.

Strategies and tools in illumina and nanopore-integrated metagenomic analysis of microbiome data

Metagenomic strategy serves as the foundation for the ecological exploration of novel bioresources (e.g., industrial enzymes and bioactive molecules) and biohazards (e.g., pathogens and antibiotic resistance genes) in natural and engineered microbial systems across multiple disciplines. Recent advancements in sequencing technology have fostered rapid development in the field of microbiome research where an increasing number of studies have applied both illumina short reads (SRs) and nanopore long reads (LRs) sequencing in their metagenomic workflow. However, given the high complexity of an environmental microbiome data set and the bioinformatic challenges caused by the unique features of these sequencing technologies, integrating SRs and LRs is not as straightforward as one might assume. The fast renewal of existing tools and growing diversity of new algorithms make access to this field even more difficult. Therefore, here we systematically summarized the complete workflow from DNA extraction to data processing strategies for applying illumina and nanopore-integrated metagenomics in the investigation in environmental microbiomes. Overall, this review aims to provide a timely knowledge framework for researchers that are interested in or are struggling with the SRs and LRs integration in their metagenomic analysis. The discussions presented will facilitate improved ecological understanding of community functionalities and assembly of natural, engineered, and human microbiomes, benefiting researchers from multiple disciplines.

Comprehensive analysis of 84 Faecalibacterium prausnitzii strains uncovers their genetic diversity, functional characteristics, and potential risks

Faecalibacterium prausnitzii is a beneficial human gut microbe and a candidate for next-generation probiotics. With probiotics now being used in clinical treatments, concerns about their safety and side effects need to be considered. Therefore, it is essential to obtain a comprehensive understanding of the genetic diversity, functional characteristics, and potential risks of different F. prausnitzii strains. In this study, we collected the genetic information of 84 F . prausnitzii strains to conduct a pan-genome analysis with multiple perspectives. Based on single-copy genes and the sequences of 16S rRNA and the compositions of the pan-genome, different phylogenetic analyses of F. prausnitzii strains were performed, which showed the genetic diversity among them. Among the proteins of the pan-genome, we found that the accessory clusters made a greater contribution to the primary genetic functions of F. prausnitzii strains than the core and specific clusters. The functional annotations of F. prausnitzii showed that only a very small number of proteins were related to human diseases and there were no secondary metabolic gene clusters encoding harmful products. At the same time, complete fatty acid metabolism was detected in F. prausnitzii. In addition, we detected harmful elements, including antibiotic resistance genes, virulence factors, and pathogenic genes, and proposed the probiotic potential risk index (PPRI) and probiotic potential risk score (PPRS) to classify these 84 strains into low-, medium-, and high-risk groups. Finally, 15 strains were identified as low-risk strains and prioritized for clinical application. Undoubtedly, our results provide a comprehensive understanding and insight into F. prausnitzii, and PPRI and PPRS can be applied to evaluate the potential risks of probiotics in general and to guide the application of probiotics in clinical application.

Long-Read Metagenome-Assembled Genomes Improve Identification of Novel Complete Biosynthetic Gene Clusters in a Complex Microbial Activated Sludge Ecosystem

Microorganisms produce a wide variety of secondary/specialized metabolites (SMs), the majority of which are yet to be discovered. These natural products play multiple roles in microbiomes and are important for microbial competition, communication, and success in the environment. SMs have been our major source of antibiotics and are used in a range of biotechnological applications. In silico mining for biosynthetic gene clusters (BGCs) encoding the production of SMs is commonly used to assess the genetic potential of organisms. However, as BGCs span tens to over 200 kb, identifying complete BGCs requires genome data that has minimal assembly gaps within the BGCs, a prerequisite that was previously only met by individually sequenced genomes. Here, we assess the performance of the currently available genome mining platform antiSMASH on 1,080 high-quality metagenome-assembled bacterial genomes (HQ MAGs) previously produced from wastewater treatment plants (WWTPs) using a combination of long-read (Oxford Nanopore) and short-read (Illumina) sequencing technologies. More than 4,200 different BGCs were identified, with 88% of these being complete. Sequence similarity clustering of the BGCs implies that the majority of this biosynthetic potential likely encodes novel compounds, and few BGCs are shared between genera. We identify BGCs in abundant and functionally relevant genera in WWTPs, suggesting a role of secondary metabolism in this ecosystem. We find that the assembly of HQ MAGs using long-read sequencing is vital to explore the genetic potential for SM production among the uncultured members of microbial communities. IMPORTANCE Cataloguing secondary metabolite (SM) potential using genome mining of metagenomic data has become the method of choice in bioprospecting for novel compounds. However, accurate biosynthetic gene cluster (BGC) detection requires unfragmented genomic assemblies, which have been technically difficult to obtain from metagenomes until very recently with new long-read technologies. Here, we determined the biosynthetic potential of activated sludge (AS), the microbial community used in resource recovery and wastewater treatment, by mining high-quality metagenome-assembled genomes generated from long-read data. We found over 4,000 BGCs, including BGCs in abundant process-critical bacteria, with no similarity to the BGCs of characterized products. We show how long-read MAGs are required to confidently assemble complete BGCs, and we determined that the AS BGCs from different studies have very little overlap, suggesting that AS is a rich source of biosynthetic potential and new bioactive compounds.

Nanopore long-read-only metagenomics enables complete and high-quality genome reconstruction from mock and complex metagenomes

BACKGROUND

The accurate and comprehensive analyses of genome-resolved metagenomics largely depend on the reconstruction of reference-quality (complete and high-quality) genomes from diverse microbiomes. Closing gaps in draft genomes have been approaching with the inclusion of Nanopore long reads; however, genome quality improvement requires extensive and time-consuming high-accuracy short-read polishing.

RESULTS

Here, we introduce NanoPhase, an open-source tool to reconstruct reference-quality genomes from complex metagenomes using only Nanopore long reads. Using Kit 9 and Q20+ chemistries, we first evaluated the feasibility of NanoPhase using a ZymoBIOMICS gut microbiome standard (including 21 strains), then sequenced the complex activated sludge microbiome and reconstructed 275 MAGs with median completeness of ~ 90%. As a result, NanoPhase improved the MAG contiguity (median MAG N50: 735 Kb, 44-86X compared to conventional short-read-based methods) while maintaining high accuracy, allowing for a full and accurate investigation of target microbiomes. Additionally, leveraging these high-contiguity reference-quality genomes, we identified 165 prophages within 111 MAGs, with 5 as active prophages, indicating the prophage was a neglected source of genetic diversity within microbial populations and influencer in shaping microbial composition in the activated sludge microbiome.

CONCLUSIONS

Our results demonstrated that NanoPhase enables reference-quality genome reconstruction from complex metagenomes directly using only Nanopore long reads. Furthermore, besides the 16S rRNA genes and biosynthetic gene clusters, the generated high-accuracy and high-contiguity MAGs improved the host identification of critical mobile genetic elements, e.g., prophage, serving as a genomic blueprint to investigate the microbial potential and ecology in the activated sludge ecosystem. Video Abstract.

Niche differentiation and symbiotic association among ammonia/nitrite oxidizers in a full-scale rotating biological contactor

Although the distribution of ammonia/nitrite oxidizers had been profiled in different habitats, current understanding is still limited regarding their niche differentiation in the integrated biofilm reactors, the symbiotic associations of ammonia/nitrite oxidizers, as well as the parasitic interaction between viruses and those functional organisms involved in the nitrogen cycle. Here, the integrated metagenomics and metatranscriptomics are applied to profile the ammonia/nitrite oxidizers communities and transcriptional activities changes along the flowpath of a concatenated full-scale rotating biological contactor (RBC) (frontend Stage-A and backend Stage-B). 19 metagenome-assembled genomes (MAGs) of ammonia/nitrite oxidizers were recovered by using a hybrid assembly approach, including four ammonia-oxidizing bacteria (AOB), two ammonia-oxidizing archaea (AOA), two complete ammonia oxidation bacteria (comammox), eight nitrite-oxidizing bacteria (NOB), and three anaerobic ammonium oxidation bacteria (anammox). Diverse AOB and anammox dominated Stage-A and collectively contributed to nitrogen conversion. With the decline of ammonia concentration along the flowpath, comammox and AOA appeared and increased in relative abundance in Stage-B, accounting for 8.8% of the entire community at the end of this reactor, and their dominating role in nitrogen turnover was indicated by the high transcription activity of their corresponding function genes. Moreover, the variation in the abundance of viruses infecting ammonia and nitrite oxidizers suggests that viruses likely act as a biotic factor mediating ammonia/nitrite oxidizer populations. This study demonstrates that complex factors shaped niche differentiation and symbiotic associations of ammonia/nitrite oxidizers in the RBC and highlights the importance of RBCs as model systems for the investigation of biotic and abiotic factors affecting the composition of microbiomes.