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

Long-Read Metagenomics of Marine Microbes Reveals Diversely Expressed Secondary Metabolites

Microbial secondary metabolites play crucial roles in microbial competition, communication, resource acquisition, antibiotic production, and a variety of other biotechnological processes. The retrieval of full-length BGC (biosynthetic gene cluster) sequences from uncultivated bacteria is difficult due to the technical constraints of short-read sequencing, making it impossible to determine BGC diversity. Using long-read sequencing and genome mining, 339 mainly full-length BGCs were recovered in this study, illuminating the wide range of BGCs from uncultivated lineages discovered in seawater from Aoshan Bay, Yellow Sea, China. Many extremely diverse BGCs were discovered in bacterial phyla such as Proteobacteria, Bacteroidota, Acidobacteriota, and Verrucomicrobiota as well as the previously uncultured archaeal phylum "Candidatus Thermoplasmatota." The data from metatranscriptomics showed that 30.1% of secondary metabolic genes were being expressed, and they also revealed the expression pattern of BGC core biosynthetic genes and tailoring enzymes. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional expression of BGCs in environmental processes. IMPORTANCE Genome mining of metagenomic data has become the preferred method for the bioprospecting of novel compounds by cataloguing secondary metabolite potential. However, the accurate detection of BGCs requires unfragmented genomic assemblies, which have been technically difficult to obtain from metagenomes until recently with new long-read technologies. We used high-quality metagenome-assembled genomes generated from long-read data to determine the biosynthetic potential of microbes found in the surface water of the Yellow Sea. We recovered 339 highly diverse and mostly full-length BGCs from largely uncultured and underexplored bacterial and archaeal phyla. Additionally, we present long-read metagenomic sequencing combined with metatranscriptomic analysis as a potential method for gaining access to the largely underutilized genetic reservoir of specialized metabolite gene clusters in the majority of microbes that are not cultured. The combination of long-read metagenomic and metatranscriptomic analyses is significant because it can more accurately assess the mechanisms of microbial adaptation to the environment through BGC expression based on metatranscriptomic data.

Microbe-mediated simultaneous nitrogen reduction and sulfamethoxazole/N-acetylsulfamethoxazole removal in lab-scale constructed wetlands

Constructed wetlands (CWs) are increasingly used to treat complex pollution such as nitrogen and emerging organic micropollutants from anthropogenic sources. In this study, the denitrification, anaerobic ammonium oxidation, dissimilatory nitrate reduction to ammonium, and nitrous oxide release rates following exposure to the frequently detected sulfonamides sulfamethoxazole (SMX) and its human metabolite, N-acetylsulfamethoxazole (N-SMX), were investigated in lab-scale CWs. Over a period of 190 d, the denitrification rates were noticeably inhibited in the SMX and N-SMX groups at week 5. Subsequently, the denitrification rates recovered, accompanied by an increase in the relevant nitrogen reduction and antibiotic resistance genes (ARGs). The composition of the microbial community also changed during this process. After the denitrification rates recovered, Burkholderia_Paraburkholderia and Gordonia exhibited a significant positive correlation with SMX exposure, which simultaneously reduced nitrate concentrations and degraded antibiotics. Burkholderia_Paraburkholderia is a key carrier of ARGs. Finally, nitrogen reduction (> 90%) and antibiotic removal (> 80%) also recovered in both SMX- and N-SMX-exposed lab-scale CWs during the operation, which revealed the interaction of SMX or N-SMX removal and nitrogen reduction.

CheckM2: a rapid, scalable and accurate tool for assessing microbial genome quality using machine learning

Advances in sequencing technologies and bioinformatics tools have dramatically increased the recovery rate of microbial genomes from metagenomic data. Assessing the quality of metagenome-assembled genomes (MAGs) is a critical step before downstream analysis. Here, we present CheckM2, an improved method of predicting genome quality of MAGs using machine learning. Using synthetic and experimental data, we demonstrate that CheckM2 outperforms existing tools in both accuracy and computational speed. In addition, CheckM2's database can be rapidly updated with new high-quality reference genomes, including taxa represented only by a single genome. We also show that CheckM2 accurately predicts genome quality for MAGs from novel lineages, even for those with reduced genome size (for example, Patescibacteria and the DPANN superphylum). CheckM2 provides accurate genome quality predictions across bacterial and archaeal lineages, giving increased confidence when inferring biological conclusions from MAGs.

A programmable pAgo nuclease with RNA target-cleavage specificity from the mesophilic bacterium Verrucomicrobia

Argonaute (Ago) proteins are conserved programmable nucleases present in eukaryotes and prokaryotes and provide defense against mobile genetic elements. Almost all characterized pAgos prefer to cleave DNA targets. Here, we describe a novel pAgo from Verrucomicrobia bacterium (VbAgo) that can specifically cleave RNA targets rather than DNA targets at 37°C and function as a multiple-turnover enzyme showing prominent catalytic capacity. VbAgo utilizes DNA guides (gDNAs) to cleave RNA targets at the canonical cleavage site. Meanwhile, the cleavage activity is remarkably strengthened at low concentrations of NaCl. In addition, VbAgo presents a weak tolerance for mismatches between gDNAs and RNA targets, and single-nucleotide mismatches at positions 11‒12 and dinucleotide mismatches at positions 3‒15 dramatically reduce target cleavage. Moreover, VbAgo can efficiently cleave highly structured RNA targets at 37°C. These properties of VbAgo broaden our understanding of Ago proteins and expand the pAgo-based RNA manipulation toolbox.

An efficient method for high molecular weight bacterial DNA extraction suitable for shotgun metagenomics from skin swabs

The human skin microbiome represents a variety of complex microbial ecosystems that play a key role in host health. Molecular methods to study these communities have been developed but have been largely limited to low-throughput quantification and short amplicon-based sequencing, providing limited functional information about the communities present. Shotgun metagenomic sequencing has emerged as a preferred method for microbiome studies as it provides more comprehensive information about the species/strains present in a niche and the genes they encode. However, the relatively low bacterial biomass of skin, in comparison to other areas such as the gut microbiome, makes obtaining sufficient DNA for shotgun metagenomic sequencing challenging. Here we describe an optimised high-throughput method for extraction of high molecular weight DNA suitable for shotgun metagenomic sequencing. We validated the performance of the extraction method, and analysis pipeline on skin swabs collected from both adults and babies. The pipeline effectively characterised the bacterial skin microbiota with a cost and throughput suitable for larger longitudinal sets of samples. Application of this method will allow greater insights into community compositions and functional capabilities of the skin microbiome.

Inter-kingdom interactions and stability of methanogens revealed by machine-learning guided multi-omics analysis of industrial-scale biogas plants

Multi-omics analysis is a powerful tool for the detection and study of inter-kingdom interactions, such as those between bacterial and archaeal members of complex biogas-producing microbial communities. In the present study, the microbiomes of three industrial-scale biogas digesters, each fed with different substrates, were analysed using a machine-learning guided genome-centric metagenomics framework complemented with metatranscriptome data. This data permitted us to elucidate the relationship between abundant core methanogenic communities and their syntrophic bacterial partners. In total, we detected 297 high-quality, non-redundant metagenome-assembled genomes (nrMAGs). Moreover, the assembled 16 S rRNA gene profiles of these nrMAGs showed that the phylum Firmicutes possessed the highest copy number, while the representatives of the archaeal domain had the lowest. Further investigation of the three anaerobic microbial communities showed characteristic alterations over time but remained specific to each industrial-scale biogas plant. The relative abundance of various microorganisms as revealed by metagenome data was independent from corresponding metatranscriptome activity data. Archaea showed considerably higher activity than was expected from their abundance. We detected 51 nrMAGs that were present in all three biogas plant microbiomes with different abundances. The core microbiome correlated with the main chemical fermentation parameters, and no individual parameter emerged as a predominant shaper of community composition. Various interspecies H₂/electron transfer mechanisms were assigned to hydrogenotrophic methanogens in the biogas plants that ran on agricultural biomass and wastewater. Analysis of metatranscriptome data revealed that methanogenesis pathways were the most active of all main metabolic pathways.

Metagenomic next generation sequencing for studying antibiotic resistance genes in the environment

Bacterial antimicrobial resistance (AMR) is a persisting and growing threat to human health. Characterization of antibiotic resistance genes (ARGs) in the environment is important to understand and control ARG-associated microbial risks. Numerous challenges exist in monitoring ARGs in the environment, due to the extraordinary diversity of ARGs, low abundance of ARGs with respect to the complex environmental microbiomes, difficulties in linking ARGs with bacterial hosts by molecular methods, difficulties in achieving quantification and high throughput simultaneously, difficulties in assessing mobility potential of ARGs, and difficulties in determining the specific AMR determinant genes. Advances in the next generation sequencing (NGS) technologies and related computational and bioinformatic tools are facilitating rapid identification and characterization ARGs in genomes and metagenomes from environmental samples. This chapter discusses NGS-based strategies, including amplicon-based sequencing, whole genome sequencing, bacterial population-targeted metagenome sequencing, metagenomic NGS, quantitative metagenomic sequencing, and functional/phenotypic metagenomic sequencing. Current bioinformatic tools for analyzing sequencing data for studying environmental ARGs are also discussed.

CRISPR-Cas phage defense systems and prophages in Candidatus Accumulibacter

Candidatus Accumulibacter plays a major role in enhanced biological phosphorus removal (EBPR) from wastewater. Although bacteriophages have been shown to represent fatal threats to Ca. Accumulibacter organisms and thus interfere with the stability of the EBPR process, little is known about the ability of different Ca. Accumulibacter strains to resist phage infections. We conducted a systematic analysis of the occurrence and characteristics of clustered regularly interspaced short palindromic repeats and associated proteins (CRISPR-Cas) systems and prophages in Ca. Accumulibacter lineage members (43 in total, including 10 newly recovered genomes). Results indicate that 28 Ca. Accumulibacter genomes encode CRISPR-Cas systems. They were likely acquired via horizontal gene transfer, conveying a distinct adaptivity to phage predation to different Ca. Accumulibacter members. Major differences in the number of spacers show the unique phage resistance of these members. A comparison of the spacers in closely related Ca. Accumulibacter members from distinct geographical locations indicates that habitat isolation may have resulted in the acquisition of resistance to different phages by different Ca. Accumulibacter. Long-term operation of three laboratory-scale EBPR bioreactors revealed high relative abundances of Ca. Accumulibacter with CRISPSR-Cas systems. Their specific resistance to phages in these reactors was indicated by spacer analysis. Metatranscriptomic analyses showed the activation of the CRISPR-Cas system under both anaerobic and aerobic conditions. Additionally, 133 prophage regions were identified in 43 Ca. Accumulibacter genomes. Twenty-seven of them (in 19 genomes) were potentially active. Major differences in the occurrence of CRISPR-Cas systems and prophages in Ca. Accumulibacter will lead to distinct responses to phage predation. This study represents the first systematic analysis of CRISPR-Cas systems and prophages in the Ca. Accumulibacter lineage, providing new perspectives on the potential impacts of phages on Ca. Accumulibacter and EBPR systems.

Comparing genomes recovered from time-series metagenomes using long- and short-read sequencing technologies

BACKGROUND

Over the past years, sequencing technologies have expanded our ability to examine novel microbial metabolisms and diversity previously obscured by isolation approaches. Long-read sequencing promises to revolutionize the metagenomic field and recover less fragmented genomes from environmental samples. Nonetheless, how to best benefit from long-read sequencing and whether long-read sequencing can provide recovered genomes of similar characteristics as short-read approaches remains unclear.

RESULTS

We recovered metagenome-assembled genomes (MAGs) from the free-living fraction at four-time points during a spring bloom in the North Sea. The taxonomic composition of all MAGs recovered was comparable between technologies. However, differences consisted of higher sequencing depth for contigs and higher genome population diversity in short-read compared to long-read metagenomes. When pairing population genomes recovered from both sequencing approaches that shared ≥ 99% average nucleotide identity, long-read MAGs were composed of fewer contigs, a higher N50, and a higher number of predicted genes when compared to short-read MAGs. Moreover, 88% of the total long-read MAGs carried a 16S rRNA gene compared to only 23% of MAGs recovered from short-read metagenomes. Relative abundances for population genomes recovered using both technologies were similar, although disagreements were observed for high and low GC content MAGs.

CONCLUSIONS

Our results highlight that short-read technologies recovered more MAGs and a higher number of species than long-read due to an overall higher sequencing depth. Long-read samples produced higher quality MAGs and similar species composition compared to short-read sequencing. Differences in the GC content recovered by each sequencing technology resulted in divergences in the diversity recovered and relative abundance of MAGs within the GC content boundaries.

3D ZnO/Activated Carbon Alginate Beads for the Removal of Antibiotic-Resistant Bacteria and Antibiotic Resistance Genes

The worldwide prevalence of antibiotic-resistant bacteria (ARB) and antibiotic resistance genes (ARGs) have become one of the most urgent issues for public health. Thus, it is critical to explore more sustainable methods with less toxicity for the long-term removal of both ARB and ARGs. In this study, we fabricated a novel material by encapsulating zinc oxide (ZnO) nanoflowers and activated carbon (AC) in an alginate biopolymer. When the dosage of ZnO was 1.0 g (≈2 g/L), the composite beads exhibited higher removal efficiency and a slight release of Zn²⁺ in water treatment. Fixed bed column experiments demonstrated that ZnO/AC alginate beads had excellent removal capacities. When the flow rate was 1 mL/min, and the initial concentration was 10⁷ CFU/mL, the removal efficiency of ARB was 5.69-log, and the absolute abundance of ARGs was decreased by 2.44-2.74-log. Moreover, the mechanism demonstrated that ZnO significantly caused cell lysis, cytoplasmic leakage, and the increase of reactive oxygen species induced subsequent oxidative stress state. These findings suggested that ZnO/AC alginate beads can be a promising material for removing ARB and ARGs from wastewater with eco-friendly and sustainable properties.

Predicting microbial community compositions in wastewater treatment plants using artificial neural networks

BACKGROUND

Activated sludge (AS) of wastewater treatment plants (WWTPs) is one of the world's largest artificial microbial ecosystems and the microbial community of the AS system is closely related to WWTPs' performance. However, how to predict its community structure is still unclear.

RESULTS

Here, we used artificial neural networks (ANN) to predict the microbial compositions of AS systems collected from WWTPs located worldwide. The predictive accuracy R²_1:1 of the Shannon-Wiener index reached 60.42%, and the average R²_1:1 of amplicon sequence variants (ASVs) appearing in at least 10% of samples and core taxa were 35.09% and 42.99%, respectively. We also found that the predictability of ASVs was significantly positively correlated with their relative abundance and occurrence frequency, but significantly negatively correlated with potential migration rate. The typical functional groups such as nitrifiers, denitrifiers, polyphosphate-accumulating organisms (PAOs), glycogen-accumulating organisms (GAOs), and filamentous organisms in AS systems could also be well recovered using ANN models, with R²_1:1 ranging from 32.62% to 56.81%. Furthermore, we found that whether industry wastewater source contained in inflow (IndConInf) had good predictive abilities, although its correlation with ASVs in the Mantel test analysis was weak, which suggested important factors that cannot be identified using traditional methods may be highlighted by the ANN model.

CONCLUSIONS

We demonstrated that the microbial compositions and major functional groups of AS systems are predictable using our approach, and IndConInf has a significant impact on the prediction. Our results provide a better understanding of the factors affecting AS communities through the prediction of the microbial community of AS systems, which could lead to insights for improved operating parameters and control of community structure. Video Abstract.

Viruses of sulfur oxidizing phototrophs encode genes for pigment, carbon, and sulfur metabolisms

Viral infections modulate bacterial metabolism and ecology. Here, we investigated the hypothesis that viruses influence the ecology of purple and green sulfur bacteria in anoxic and sulfidic lakes, analogs of euxinic oceans in the geologic past. By screening metagenomes from lake sediments and water column, in addition to publicly-available genomes of cultured purple and green sulfur bacteria, we identified almost 300 high and medium-quality viral genomes. Viruses carrying the gene psbA, encoding the small subunit of photosystem II protein D1, were ubiquitous, suggesting viral interference with the light reactions of sulfur oxidizing autotrophs. Viruses predicted to infect these autotrophs also encoded auxiliary metabolic genes for reductive sulfur assimilation as cysteine, pigment production, and carbon fixation. These observations show that viruses have the genomic potential to modulate the production of metabolic markers of phototrophic sulfur bacteria that are used to identify photic zone euxinia in the geologic past.

Microbial Community Succession Along a Chronosequence in Constructed Salt Marsh Soils

In this study, we examined the succession of soil microbial communities across a chronosequence of newly constructed salt marshes constructed primarily of fine-grained dredge material, using 16S rRNA amplicon sequences. Alpha diversity in the subsurface horizons was initially low and increased to reference levels within 3 years of marsh construction, while alpha diversity in the newly accumulating organic matter-rich surface soils was initially high and remained unchanged. Microbial community succession was fastest in the surface horizon (~ 24 years to reference equivalency) and became progressively slower with depth in the subsurface horizons (~ 30-67 years). Random forest linear regression analysis was used to identify important taxa driving the trajectories toward reference conditions. In the parent material, putative sulfate-reducers (Desulfobacterota), methanogens (Crenarchaeota, especially Methanosaeta), and fermenters (Chloroflexi and Clostridia) increased over time, suggesting an enrichment of these metabolisms over time, similar to natural marshes. Concurrently in the surface soils, the relative abundances of putative methane-, methyl-, and sulfide oxidizers, especially among Gammaproteobacteria, increased over time, suggesting the co-development of sulfide and methane removal metabolisms in marsh soils. Finally, we observed that the surface soil communities at one of the marshes did not follow the trajectory of the others, exhibiting a greater relative abundance of anaerobic taxa. Uniquely in this dataset, this marsh was developing signs of excessive inundation stress in terms of vegetation coverage and soil geochemistry. Therefore, we suggest that soil microbial community structure may be effective bioindicators of salt marsh inundation and are worthy of further targeted investigation.

Closed genomes uncover a saltwater species of Candidatus Electronema and shed new light on the boundary between marine and freshwater cable bacteria

Cable bacteria of the Desulfobulbaceae family are centimeter-long filamentous bacteria, which are capable of conducting long-distance electron transfer. Currently, all cable bacteria are classified into two candidate genera: Candidatus Electronema, typically found in freshwater environments, and Candidatus Electrothrix, typically found in saltwater environments. This taxonomic framework is based on both 16S rRNA gene sequences and metagenome-assembled genome (MAG) phylogenies. However, most of the currently available MAGs are highly fragmented, incomplete, and thus likely miss key genes essential for deciphering the physiology of cable bacteria. Also, a closed, circular genome of cable bacteria has not been published yet. To address this, we performed Nanopore long-read and Illumina short-read shotgun sequencing of selected environmental samples and a single-strain enrichment of Ca. Electronema aureum. We recovered multiple cable bacteria MAGs, including two circular and one single-contig. Phylogenomic analysis, also confirmed by 16S rRNA gene-based phylogeny, classified one circular MAG and the single-contig MAG as novel species of cable bacteria, which we propose to name Ca. Electronema halotolerans and Ca. Electrothrix laxa, respectively. The Ca. Electronema halotolerans, despite belonging to the previously recognized freshwater genus of cable bacteria, was retrieved from brackish-water sediment. Metabolic predictions showed several adaptations to a high salinity environment, similar to the "saltwater" Ca. Electrothrix species, indicating how Ca. Electronema halotolerans may be the evolutionary link between marine and freshwater cable bacteria lineages.

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.

The Hydroponic Rockwool Root Microbiome: Under Control or Underutilised?

Land plants have an ancient and intimate relationship with microorganisms, which influences the composition of natural ecosystems and the performance of crops. Plants shape the microbiome around their roots by releasing organic nutrients into the soil. Hydroponic horticulture aims to protect crops from damaging soil-borne pathogens by replacing soil with an artificial growing medium, such as rockwool, an inert material made from molten rock spun into fibres. Microorganisms are generally considered a problem to be managed, to keep the glasshouse clean, but the hydroponic root microbiome assembles soon after planting and flourishes with the crop. Hence, microbe–plant interactions play out in an artificial environment that is quite unlike the soil in which they evolved. Plants in a near-ideal environment have little dependency on microbial partners, but our growing appreciation of the role of microbial communities is revealing opportunities to advance practices, especially in agriculture and human health. Hydroponic systems are especially well-suited to active management of the root microbiome because they allow complete control over the root zone environment; however, they receive much less attention than other host–microbiome interactions. Novel techniques for hydroponic horticulture can be identified by extending our understanding of the microbial ecology of this unique environment.

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.

Electricity-Driven Microbial Metabolism of Carbon and Nitrogen: A Waste-to-Resource Solution

Electricity-driven microbial metabolism relies on the extracellular electron transfer (EET) process between microbes and electrodes and provides promise for resource recovery from wastewater and industrial discharges. Over the past decades, tremendous efforts have been dedicated to designing electrocatalysts and microbes, as well as hybrid systems to push this approach toward industrial adoption. This paper summarizes these advances in order to facilitate a better understanding of electricity-driven microbial metabolism as a sustainable waste-to-resource solution. Quantitative comparisons of microbial electrosynthesis and abiotic electrosynthesis are made, and the strategy of electrocatalyst-assisted microbial electrosynthesis is critically discussed. Nitrogen recovery processes including microbial electrochemical N₂ fixation, electrocatalytic N₂ reduction, dissimilatory nitrate reduction to ammonium (DNRA), and abiotic electrochemical nitrate reduction to ammonia (Abio-NRA) are systematically reviewed. Furthermore, the synchronous metabolism of carbon and nitrogen using hybrid inorganic-biological systems is discussed, including advanced physicochemical, microbial, and electrochemical characterizations involved in this field. Finally, perspectives for future trends are presented. The paper provides valuable insights on the potential contribution of electricity-driven microbial valorization of waste carbon and nitrogen toward a green and sustainable society.

Nanopore Sequencing Discloses Compositional Quality of Commercial Probiotic Feed Supplements

The market for the application of probiotics as a livestock health improvement supplement has increased in recent years. However, most of the available products are quality-controlled using low-resolution techniques and un-curated databases, resulting in misidentification and incorrect product labels. In this work, we deployed two workflows and compared results obtained by full-length 16S rRNA genes (16S) and metagenomic (Meta) data to investigate their reliability for the microbial composition of both liquid and solid forms of animal probiotic products using Oxford Nanopore long-read-only (without short-read). Our result revealed that 16S amplicon data permits to detect the bacterial microbiota even with the low abundance in the samples. Moreover, the 16S approach has the potential to provide species-level resolution for prokaryotes but not for assessing yeast communities. Whereas, Meta data has more power to recover of high-quality metagenome-assembled genomes that enables detailed exploration of both bacterial and yeast populations, as well as antimicrobial resistance genes, and functional genes in the population. Our findings clearly demonstrate that implementing these workflows with long-read-only monitoring could be applied to assessing the quality and safety of probiotic products for animals and evaluating the quality of probiotic products on the market. This would benefit the sustained growth of the livestock probiotic industry.

Conservation of Genomic Information in Multiple Displacement Amplified Low-Quantity Metagenomic Material from Marine Invertebrates

Marine invertebrate microbiomes have been a rich source of bioactive compounds and interesting genomic features. In cases where the achievable amounts of metagenomic DNA are too low for direct sequencing, multiple displacement amplification (MDA) can be used for whole genome amplification. However, MDA has known limitations which can affect the quality of the resulting genomes and metagenomes. In this study, we evaluated the conservation of biosynthetic gene clusters (BGCs) and enzymes in MDA products from low numbers of prokaryotic cells (estimated 2–850). Marine invertebrate microbiomes collected from Arctic and sub-Arctic areas served as source material. The cells were separated from the host tissue, lysed, and directly subjected to MDA. The MDA products were sequenced by Illumina sequencing. Corresponding numbers of bacteria from a set of three reference strains were treated the same way. The study demonstrated that useful information on taxonomic, BGC, and enzyme diversities was obtainable from such marginal quantities of metagenomic material. Although high levels of assembly fragmentation resulted in most BGCs being incomplete, we conclude that this genome mining approach has the potential to reveal interesting BGCs and genes from hard-to-reach biological sources.

Towards low carbon demand and highly efficient nutrient removal: Establishing denitrifying phosphorus removal in a biofilm-based system

The combined denitrifying phosphorus removal (DPR) and Anammox process is expected to achieve advanced nutrient removal with low carbon consumption. However, exchanging ammonia/nitrate between them is one limitation. This study investigated the feasibility of conducting DPR in a biofilm reactor to solve that problem. After 46-day anaerobic/aerobic operation, high phosphorus removal efficiency (PRE, 83.15 %) was obtained in the activated sludge (AS) and biofilm co-existed system, in which the AS performed better. Phosphate-accumulating organisms might quickly adapt to the anoxic introduced nitrate, but the following aerobic stage ensured a low effluent orthophosphate (<1.03 mg/L). Because of waste sludge discharging and AS transforming to biofilm, the suspended solids dropped below 60 mg/L on Day 100, resulting in PRE decline (17.17 %) and effluent orthophosphate rise (4.23 mg/L). Metagenomes analysis revealed that Pseudomonas and Thiothrix had genes for denitrification and encoding Pit phosphate transporter, and Candidatus_Competibacter was necessary for biofilm formation.

A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics

BACKGROUND

Lake Baikal, the world's deepest freshwater lake, contains important numbers of Candidatus Patescibacteria (formerly CPR) in its deepest reaches. However, previously obtained CPR metagenome-assembled genomes recruited very poorly indicating the potential of other groups being present. Here, we have applied for the first time a long-read (PacBio CCS) metagenomic approach to analyze in depth the Ca. Patescibacteria living in the bathypelagic water column of Lake Baikal at 1600 m.

RESULTS

The retrieval of nearly complete 16S rRNA genes before assembly has allowed us to detect the presence of a novel and a likely endemic group of Ca. Patescibacteria inhabiting bathypelagic Lake Baikal. This novel group seems to possess extremely high intra-clade diversity, precluding complete genomes' assembly. However, read binning and scaffolding indicate that these microbes are similar to other Ca. Patescibacteria (i.e. parasites or symbionts), although they seem to carry more anabolic pathways, likely reflecting the extremely oligotrophic habitat they inhabit. The novel bins have not been found anywhere, but one of the groups appears in small amounts in an oligotrophic and deep alpine Lake Thun. We propose this novel group be named Baikalibacteria.

CONCLUSION

The recovery of 16S rRNA genes via long-read metagenomics plus the use of long-read binning to uncover highly diverse "hidden" groups of prokaryotes are key strategies to move forward in ecogenomic microbiology. The novel group possesses enormous intraclade diversity akin to what happens with Ca. Patescibacteria at the interclade level, which is remarkable in an environment that has changed little in the last 25 million years.

Active predation, phylogenetic diversity, and global prevalence of myxobacteria in wastewater treatment plants

The operation of modern wastewater treatment plants (WWTPs) is driven by activated sludge microbiota, a complex assemblage of trophically interacting microorganisms. Microbial predation is crucial to fundamental understanding of how biological interactions drive microbiome structuring and functioning of WWTPs. However, predatory bacteria have received little attention regarding their diversity, activity, and ecological function in activated sludge, limiting the exploitation of food web interactions for wastewater microbiome engineering. Here, by using rRNA-stable isotope probing of activated sludge microbiota with ¹³C-labeled prey bacteria, we uncovered diverse as-yet-uncultivated putative predatory bacteria that actively incorporated ¹³C-biomass. Myxobacteria, especially Haliangium and the mle1-27 clade, were found as the dominant active predators, refreshing conventional views based on a few predatory isolates of Bdellovibrionota from WWTPs. The identified predatory bacteria showed more selective predation on prey compared with the protists dominated by ciliates, providing in situ evidence for inter-domain predation behavior divergence in activated sludge. Putative predatory bacteria were tracked over a two-year microbiome monitoring effort at a local WWTP, revealing the predominance of Myxococcota (6.5 ± 1.3%) over Bdellovibrionota (1.0 ± 0.2%) lineages. Phylogenetic analysis unveiled highly diverse myxobacteria inhabiting activated sludge and suggested a habitat filtering effect in global WWTPs. Further mining of a global activated sludge microbiome dataset revealed the prevalence of Myxococcota (5.4 ± 0.1%) species and potential impacts of myxobacterial predation on process performance. Collectively, our findings provided unique insights into the predating activity, diversity, and prevalence of Myxococcota species in activated sludge, highlighting their links with wastewater treatment processes via trophic regulation of enteric and functional bacteria.

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.

Low diversity and microdiversity of comammox bacteria in wastewater systems suggest specific adaptations within the Ca. Nitrospira nitrosa cluster

Studies have found Ca. Nitrospira nitrosa-like bacteria to be the principal or sole comammox bacteria in nitrogen removal systems for wastewater treatment. In contrast, multiple populations of strict ammonia and nitrite oxidizers co-exist in similar systems. This apparent lack of diversity is surprising and could impact the feasibility of leveraging comammox bacteria for nitrogen removal. We used full-length 16S rRNA gene sequencing and genome-resolved metagenomics to compare the species-level diversity of comammox bacteria with that of strict nitrifiers in full-scale wastewater treatment systems and assess whether this comparison is consistent or diverged at the strain-level. Full-length 16S rRNA gene sequencing indicated that Nitrosomonas-like bacteria exhibited higher species-level diversity in comparison with other nitrifying bacteria, while the strain-level diversity (also called microdiversity) of most Nitrospira-like bacteria were higher than Nitrosomonas-like bacteria with few exceptions (one Nitrospira lineage II population). Comammox bacterial metagenome assembled genomes (MAGs) were associated with Ca. Nitrospira nitrosa. The average amino acid identity between principal comammox bacterial MAGs (93% ± 3) across systems was significantly higher than that of the Nitrosomonas-like ammonia oxidizers (73% ± 8), the Nitrospira_A-like nitrite oxidizer (85% ± 4), and the Nitrospira_D-like nitrite oxidizer (83% ± 1). This demonstrated the low species-level diversity of comammox bacteria compared with strict nitrifiers and further suggests that the same comammox population was detected in all systems. Comammox bacteria (Nitrospira lineage II), Nitrosomonas and, Nitrospira_D (Nitrospira lineage II) MAGs were significantly less microdiverse than the Nitrospira_A (lineage I) MAGs. Interestingly, strain-resolved analysis also indicates that different nitrogen removal systems harbor different comammox bacterial strains within the Ca. Nitrospira nitrosa cluster. These results suggest that comammox bacteria associated with Ca. Nitrospira nitrosa have low species- and strain-level diversity in nitrogen removal systems and may thus harbor specific adaptations to the wastewater ecosystem.

Metagenomic Views of Microbial Communities in Sand Sediments Associated with Coral Reefs

Reef sediments, the home for microbes with high abundances, provide an important source of carbonates and nutrients for the growth and maintenance of coral reefs. However, there is a lack of systematic research on the composition of microbial community in sediments of different geographic sites and their potential effect on nutrient recycling and health of the coral reef ecosystem. In combination of biogeochemical measurements with gene- and genome-centric metagenomics, we assessed microbial community compositions and functional diversity, as well as profiles of antibiotic resistance genes in surface sediments of 16 coral reef sites at different depths from the Xisha islands in the South China Sea. Reef sediment microbiomes are diverse and novel at lower taxonomic ranks, dominated by Proteobacteria and Planctomycetota. Most reef sediment bacteria potentially participate in biogeochemical cycling via oxidizing various organic and inorganic compounds as energy sources. High abundances of Proteobacteria (mostly Rhizobiales and Woeseiales) are metabolically flexible and contain rhodopsin genes. Various classes of antibiotic resistance genes, hosted by diverse bacterial lineages, were identified to confer resistance to multidrug, aminoglycoside, and other antibiotics. Overall, our findings expanded the understanding of reef sediment microbial ecology and provided insights for their link to the coral reef ecosystem health.

Genetic potential for exopolysaccharide synthesis in activated sludge bacteria uncovered by genome-resolved metagenomics

A good floc formation of activated sludge (AS) is crucial for solid-liquid separation and production of clean effluent during wastewater treatment. Floc formation is partly controlled by self-produced extracellular polymeric substances (EPS) such as exopolysaccharides, proteins, and nucleic acids. Little is known about the composition, structure, and function of EPS in AS and which bacteria produce them. To address this knowledge gap for the exopolysaccharides, we took advantage of 1083 high-quality metagenome-assembled genomes (MAGs) obtained from 23 Danish wastewater treatment plants. We investigated the genomic potential for exopolysaccharide biosynthesis in bacterial species typical in AS systems based on genome mining and gene synteny analyses. Putative gene clusters associated with the biosynthesis of alginate, cellulose, curdlan, diutan, hyaluronic acids, Pel, poly-β-1,6-N-acetyl-d-glucosamine (PNAG), Psl, S88 capsular polysaccharide, salecan, succinoglycan, and xanthan were identified and linked to individual MAGs, providing a comprehensive overview of the genome-resolved potential for these exopolysaccharides in AS bacteria. The approach and results provide a starting point for a more comprehensive understanding of EPS composition in wastewater treatment systems, which may facilitate a more refined regulation of the activated sludge process for improved stability.

Combining Flow Cytometry and Metagenomics Improves Recovery of Metagenome-Assembled Genomes in a Cell Culture from Activated Sludge

The recovery of metagenome-assembled genomes is biased towards the most abundant species in a given community. To improve the identification of species, even if only dominant species are recovered, we investigated the integration of flow cytometry cell sorting with bioinformatics tools to recover metagenome-assembled genomes. We used a cell culture of a wastewater microbial community as our model system. Cells were separated based on fluorescence signals via flow cytometry cell sorting into sub-communities: dominant gates, low abundant gates, and outer gates into subsets of the original community. Metagenome sequencing was performed for all groups. The unsorted community was used as control. We recovered a total of 24 metagenome-assembled genomes (MAGs) representing 11 species-level genome operational taxonomic units (gOTUs). In addition, 57 ribosomal operational taxonomic units (rOTUs) affiliated with 29 taxa at species level were reconstructed from metagenomic libraries. Our approach suggests a two-fold increase in the resolution when comparing sorted and unsorted communities. Our results also indicate that species abundance is one determinant of genome recovery from metagenomes as we can recover taxa in the sorted libraries that are not present in the unsorted community. In conclusion, a combination of cell sorting and metagenomics allows the recovery of MAGs undetected without cell sorting.

Recovery and Analysis of Long-Read Metagenome-Assembled Genomes

The development of long-read nucleic acid sequencing is beginning to make very substantive impact on the conduct of metagenome analysis, particularly in relation to the problem of recovering the genomes of member species of complex microbial communities. Here we outline bioinformatics workflows for the recovery and characterization of complete genomes from long-read metagenome data and some complementary procedures for comparison of cognate draft genomes and gene quality obtained from short-read sequencing and long-read sequencing.

Long-Read Metagenomics and CAZyme Discovery

Microorganisms play a primary role in regulating biogeochemical cycles and are a valuable source of enzymes that have biotechnological applications, such as carbohydrate-active enzymes (CAZymes). However, the inability to culture the majority of microorganisms that exist in natural ecosystems restricts access to potentially novel bacteria and beneficial CAZymes. While commonplace molecular-based culture-independent methods such as metagenomics enable researchers to study microbial communities directly from environmental samples, recent progress in long-read sequencing technologies are advancing the field. We outline key methodological stages that are required as well as describe specific protocols that are currently used for long-read metagenomic projects dedicated to CAZyme discovery.

Predicting the impact of temperature on metabolic fluxes using resource allocation modelling: Application to polyphosphate accumulating organisms

The understanding of microbial communities and the biological regulation of its members is crucial for implementation of novel technologies using microbial ecology. One poorly understood metabolic principle of microbial communities is resource allocation and biosynthesis. Resource allocation theory in polyphosphate accumulating organisms (PAOs) is limited as a result of their slow imposed growth rate (typical sludge retention times of at least 4 days) and limitations to quantify changes in biomass components over a 6 hours cycle (less than 10% of their growth). As a result, there is no direct evidence supporting that biosynthesis is an exclusive aerobic process in PAOs that alternate continuously between anaerobic and aerobic phases. Here, we apply resource allocation metabolic flux analysis to study the optimal phenotype of PAOs over a temperature range of 4 °C to 20 °C. The model applied in this research allowed to identify optimal metabolic strategies in a core metabolic model with limited constraints based on biological principles. The addition of a constraint limiting biomass synthesis to be an exclusive aerobic process changed the metabolic behaviour and improved the predictability of the model over the studied temperature range by closing the gap between prediction and experimental findings. The results validate the assumption of limited anaerobic biosynthesis in PAOs, specifically "Candidatus Accumulibacter" related species. Interestingly, the predicted growth yield was lower, suggesting that there are mechanistic barriers for anaerobic growth not yet understood nor reflected in the current models of PAOs. Moreover, we identified strategies of resource allocation applied by PAOs at different temperatures as a result of the decreased catalytic efficiencies of their biochemical reactions. Understanding resource allocation is paramount in the study of PAOs and their currently unknown complex metabolic regulation, and metabolic modelling based on biological first principles provides a useful tool to develop a mechanistic understanding.

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.

Occurrence, Diversity, and Genomes of "Candidatus Patescibacteria" along the Early Diagenesis of Marine Sediments

The phylum "Candidatus Patescibacteria" (or Candidate Phyla Radiation [CPR]) accounts for roughly one-quarter of microbial diversity on Earth, but the presence and diversity of these bacteria in marine sediments have been rarely charted. Here, we investigate the abundance, diversity, and metabolic capacities of CPR bacteria in three sediment sites (Mohns Ridge, North Pond, and Costa Rica Margin) with samples covering a wide range of redox zones formed during the early diagenesis of organic matter. Through metagenome sequencing, we found that all investigated sediment horizons contain "Ca. Patescibacteria" (0.4 to 28% of the total communities), which are affiliated with the classes "Ca. Paceibacteria," "Ca. Gracilibacteria," "Ca. Microgenomatia," "Ca. Saccharimonadia," "Ca. ABY1," and "Ca. WWE3." However, only a subset of the diversity of marine sediment "Ca. Patescibacteria," especially the classes "Ca. Paceibacteria" and "Ca. Gracilibacteria," can be captured by 16S rRNA gene amplicon sequencing with commonly used universal primers. We recovered 11 metagenome-assembled genomes (MAGs) of CPR from these sediments, most of which are novel at the family or genus level in the "Ca. Paceibacteria" class and are missed by the amplicon sequencing. While individual MAGs are confined to specific anoxic niches, the lack of capacities to utilize the prevailing terminal electron acceptors indicates that they may not be directly selected by the local redox conditions. These CPR bacteria lack essential biosynthesis pathways and may use a truncated glycolysis pathway to conserve energy as fermentative organotrophs. Our findings suggest that marine sediments harbor some novel yet widespread CPR bacteria during the early diagenesis of organic matter, which needs to be considered in population dynamics assessments in this vast environment. IMPORTANCE Ultrasmall-celled "Ca. Patescibacteria" have been estimated to account for one-quarter of the total microbial diversity on Earth, the parasitic lifestyle of which may exert a profound control on the overall microbial population size of the local ecosystems. However, their diversity and metabolic functions in marine sediments, one of the largest yet understudied ecosystems on Earth, remain virtually uncharacterized. By applying cultivation-independent approaches to a range of sediment redox zones, we reveal that "Ca. Patescibacteria" members are rare but widespread regardless of the prevailing geochemical conditions. These bacteria are affiliated with novel branches of "Ca. Patescibacteria" and have been largely missed in marker gene-based surveys. They do not have respiration capacity but may conserve energy by fermenting organic compounds from their episymbiotic hosts. Our findings suggest that these novel "Ca. Patescibacteria" are among the previously overlooked microbes in diverse marine sediments.

The role of microbial ecology in improving the performance of anaerobic digestion of sewage sludge

The use of next-generation diagnostic tools to optimise the anaerobic digestion of municipal sewage sludge has the potential to increase renewable natural gas recovery, improve the reuse of biosolid fertilisers and help operators expand circular economies globally. This review aims to provide perspectives on the role of microbial ecology in improving digester performance in wastewater treatment plants, highlighting that a systems biology approach is fundamental for monitoring mesophilic anaerobic sewage sludge in continuously stirred reactor tanks. We further highlight the potential applications arising from investigations into sludge ecology. The principal limitation for improvements in methane recoveries or in process stability of anaerobic digestion, especially after pre-treatment or during co-digestion, are ecological knowledge gaps related to the front-end metabolism (hydrolysis and fermentation). Operational problems such as stable biological foaming are a key problem, for which ecological markers are a suitable approach. However, no biomarkers exist yet to assist in monitoring and management of clade-specific foaming potentials along with other risks, such as pollutants and pathogens. Fundamental ecological principles apply to anaerobic digestion, which presents opportunities to predict and manipulate reactor functions. The path ahead for mapping ecological markers on process endpoints and risk factors of anaerobic digestion will involve numerical ecology, an expanding field that employs metrics derived from alpha, beta, phylogenetic, taxonomic, and functional diversity, as well as from phenotypes or life strategies derived from genetic potentials. In contrast to addressing operational issues (as noted above), which are effectively addressed by whole population or individual biomarkers, broad improvement and optimisation of function will require enhancement of hydrolysis and acidogenic processes. This will require a discovery-based approach, which will involve integrative research involving the proteome and metabolome. This will utilise, but overcome current limitations of DNA-centric approaches, and likely have broad application outside the specific field of anaerobic digestion.

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.

Comparative analysis of floc characteristics and microbial communities in anoxic and aerobic suspended growth processes

A fully anoxic suspended growth process is an appealing alternative to conventional activated sludge (AS) due to considerable aeration reduction and improved carbon processing efficiency for biological nutrient removal (BNR). With development of the hybrid membrane aerated biofilm reactor (MABR) technology, implementation of a fully anoxic suspended growth community in BNR facilities became practical. To better understand potential limitations with the elimination of aeration, we carried out microscopic examination and 16S rRNA gene-based microbial community profiling to determine how an anoxic suspended growth would differ from the conventional aerobic process in floc characteristics, microbial diversity, microbial temporal dynamics, and community assembly pattern. Fewer filamentous populations were found in the anoxic mixed liquor, suggesting easily sheared flocs. The anoxic microbial community had distinct composition and structure, but its diversity and temporal dynamics were similar to the conventional aerobic community. A variety of well-studied functional guilds were also identified in the anoxic community. The anoxic microbial community assembly was more stochastic than the conventional aerobic community, but deterministic assembly was still significant with a large core microbiome adapted to the anoxic condition. PRACTITIONER POINTS: Flocs developed under the anoxic conditions had less filamentous backbones, implying reduced flocculation capacity and easily sheared flocs. Knowledge about the ecophysiology of Thauera, Thiothrix, and Trichococcus can help achieve good properties of the anoxic flocs. A diverse microbial community sustainably adapted to the fully anoxic condition, containing a variety of filaments, denitrifiers, and PAOs. The anoxic microbial community displayed a similar degree of diversity and temporal dynamics compared to the aerobic counterpart. The anoxic community's assembly was more stochastic, so it may be less subject to changes in environmental variables.

Development of microbial communities in biofilm and activated sludge in a hybrid reactor

Microorganisms play a key role in biological wastewater treatment. The form in which biomass develops determines the efficiency and mechanisms of organic compound conversion, due to different conditions in various microbial structures. However, the results of studies comparing the microbial communities in biofilm and activated sludge have often conflicted. Therefore, this study compared the composition and development of the bacterial communities in biofilm and activated sludge in a hybrid reactor, employing 16S rRNA sequencing. Statistical analysis of the sequencing data included the identification of taxa characteristic to the biofilm and activated sludge, alpha and beta diversity analysis, and network analysis. These analyses indicated that the biofilm bacterial community was richer and more diverse than the activated sludge community. The mean numbers of OTU were 1614 in the biofilm and 993 in the activated sludge, and the mean values of the Chao1 (1735 vs. 1105) and Shannon (5.3 vs. 4.3) biodiversity indices were significantly higher for the biofilm. The biofilm was a better environment for development of nitrifiers (e.g., Nitrosomonas, Nitrospira) and phosphorus accumulating organisms (Candidatus Accumulibacter). Bacteria in the biofilm co-occurrence network had more connections (based on Spearman's rank correlation coefficient) with each other, indicating that they interact more than those in the activated sludge.

Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Metagenome-assembled genomes (MAGs) have revealed the hidden diversity and functions of uncultivated microbes, but their reconstruction from metagenomes remains a computationally difficult task. Repetitive or exogenous sequences, such as ribosomal RNA and horizontally transferred genes, are frequently absent from MAGs because of misassembly and binning errors. Here, we report that ribosomal protein genes are also often absent from MAGs, although they are neither repetitive nor exogenous. Comprehensive analyses of more than 190,000 MAGs revealed that these genes could be missing in more than 20-40% of near-complete (i.e., with completeness of 90% or higher) MAGs. While some uncultivated environmental microbes intrinsically lack some ribosomal protein genes, we found that this unexpected absence is largely due to special evolutionary patterns of codon usage bias in ribosomal protein genes and algorithmic characteristics of metagenomic binning, which is dependent on tetranucleotide frequencies of contigs. This problem reflects the microbial life-history strategy. Fast-growing microbes tend to have this difficulty, likely because of strong evolutionary pressures on ribosomal protein genes toward the efficient assembly of ribosomes. Our observations caution those who study genomics and phylogeny of uncultivated microbes, the diversity and evolution of microbial genes in the central dogma, and bioinformatics in metagenomics.

Critical assessment of pan-genomic analysis of metagenome-assembled genomes

Pan-genome analyses of metagenome-assembled genomes (MAGs) may suffer from the known issues with MAGs: fragmentation, incompleteness and contamination. Here, we conducted a critical assessment of pan-genomics of MAGs, by comparing pan-genome analysis results of complete bacterial genomes and simulated MAGs. We found that incompleteness led to significant core gene (CG) loss. The CG loss remained when using different pan-genome analysis tools (Roary, BPGA, Anvi'o) and when using a mixture of MAGs and complete genomes. Contamination had little effect on core genome size (except for Roary due to in its gene clustering issue) but had major influence on accessory genomes. Importantly, the CG loss was partially alleviated by lowering the CG threshold and using gene prediction algorithms that consider fragmented genes, but to a less degree when incompleteness was higher than 5%. The CG loss also led to incorrect pan-genome functional predictions and inaccurate phylogenetic trees. Our main findings were supported by a study of real MAG-isolate genome data. We conclude that lowering CG threshold and predicting genes in metagenome mode (as Anvi'o does with Prodigal) are necessary in pan-genome analysis of MAGs. Development of new pan-genome analysis tools specifically for MAGs are needed in future studies.

Chromosome-level assembly and analysis of the Thymus genome provide insights into glandular secretory trichome formation and monoterpenoid biosynthesis in thyme

Thyme has medicinal and aromatic value because of its potent antimicrobial and antioxidant properties. However, the absence of a fully sequenced thyme genome limits functional genomic studies of Chinese native thymes. Thymus quinquecostatus Čelak., which contains large amounts of bioactive monoterpenes such as thymol and carvacrol, is an important wild medicinal and aromatic plant in China. Monoterpenoids are abundant in glandular secretory trichomes. Here, high-fidelity and chromatin conformation capture technologies were used to assemble and annotate the T. quinquecostatus genome at the chromosome level. The 13 chromosomes of T. quinquecostatus had a total length of 528.66 Mb, a contig N50 of 8.06 Mb, and a BUSCO score of 97.34%. We found that T. quinquecostatus had experienced two whole-genome duplications, with the most recent event occurring ∼4.34 million years ago. Deep analyses of the genome, in conjunction with comparative genomic, phylogenetic, transcriptomic, and metabonomic studies, uncovered many regulatory factors and genes related to monoterpenoids and glandular secretory trichome development. Genes encoding terpene synthase (TPS), cytochrome P450 monooxygenases (CYPs), short-chain dehydrogenase/reductase (SDR), R2R3-MYB, and homeodomain-leucine zipper (HD-ZIP) IV were among those present in the T. quinquecostatus genome. Notably, Tq02G002290.1 (TqTPS1) was shown to encode the terpene synthase responsible for catalyzing production of the main monoterpene product γ-terpinene from geranyl diphosphate (GPP). Our study provides significant insight into the mechanisms of glandular secretory trichome formation and monoterpenoid biosynthesis in thyme. This work will facilitate the development of molecular breeding tools to enhance the production of bioactive secondary metabolites in Lamiaceae.

TbasCO: trait-based comparative 'omics identifies ecosystem-level and niche-differentiating adaptations of an engineered microbiome

A grand challenge in microbial ecology is disentangling the traits of individual populations within complex communities. Various cultivation-independent approaches have been used to infer traits based on the presence of marker genes. However, marker genes are not linked to traits with complete fidelity, nor do they capture important attributes, such as the timing of gene expression or coordination among traits. To address this, we present an approach for assessing the trait landscape of microbial communities by statistically defining a trait attribute as a shared transcriptional pattern across multiple organisms. Leveraging the KEGG pathway database as a trait library and the Enhanced Biological Phosphorus Removal (EBPR) model microbial ecosystem, we demonstrate that a majority (65%) of traits present in 10 or more genomes have niche-differentiating expression attributes. For example, while many genomes containing high-affinity phosphorus transporter pstABCS display a canonical attribute (e.g. up-regulation under phosphorus starvation), we identified another attribute shared by many genomes where transcription was highest under high phosphorus conditions. Taken together, we provide a novel framework for unravelling the functional dynamics of uncultivated microorganisms by assigning trait-attributes through genome-resolved time-series metatranscriptomics.

Insights into ecological roles of uncultivated bacteria in Katase hot spring sediment from long-read metagenomics

Diverse yet-uncultivated bacteria and archaea, i.e., microbial dark matter, are present in terrestrial hot spring environments. Numerous metagenome-assembled genomes (MAGs) of these uncultivated prokaryotes by short-read metagenomics have been reported so far, suggesting their metabolic potential. However, more reliable MAGs, i.e., circularized complete MAGs (cMAGs), have been rarely reported from hot spring environments. Here, we report 61 high-quality (HQ)-MAGs, including 14 cMAGs, of diverse uncultivated bacteria and archaea retrieved from hot spring sediment (52°C, pH 7.2) by highly accurate long-read sequencing using PacBio Sequel II. The HQ MAGs were affiliated with one archaeal and 13 bacterial phyla. Notably, nine of the 14 cMAGs were the first reported cMAGs for the family- to class-level clades that these cMAGs belonged to. The genome information suggests that the bacteria represented by MAGs play a significant role in the biogeochemical cycling of carbon, nitrogen, iron, and sulfur at this site. In particular, the genome analysis of six HQ MAGs including two cMAGs of Armatimonadota, of which members are frequently abundant in hot spring environments, predicts that they are aerobic, moderate thermophilic chemoorganoheterotrophs, and potentially oxidize and/or reduce iron. This prediction is consistent with the environmental conditions where they were detected. Our results expand the knowledge regarding the ecological potential of uncultivated bacteria in moderately-high-temperature environments.

Advancing the protection of marine life through genomics

The rapid growth in genomic techniques provides the potential to transform how we protect, manage, and conserve marine life. Further, solutions to boost the resilience of marine species to climate change and other disturbances that characterize the Anthropocene require transformative approaches, made more effective if guided by genomic data. Although genetic techniques have been employed in marine conservation for decades and the availability of genomic data is rapidly expanding, widespread application still lags behind other data types. This Essay reviews how genetics and genomics have been utilized in management initiatives for ocean conservation and restoration, highlights success stories, and presents a pathway forward to enhance the uptake of genomic data for protecting our oceans.

Hybrid metagenome assemblies link carbohydrate structure with function in the human gut microbiome

Complex carbohydrates that escape small intestinal digestion, are broken down in the large intestine by enzymes encoded by the gut microbiome. This is a symbiotic relationship between microbes and host, resulting in metabolic products that influence host health and are exploited by other microbes. However, the role of carbohydrate structure in directing microbiota community composition and the succession of carbohydrate-degrading microbes, is not fully understood. In this study we evaluate species-level compositional variation within a single microbiome in response to six structurally distinct carbohydrates in a controlled model gut using hybrid metagenome assemblies. We identified 509 high-quality metagenome-assembled genomes (MAGs) belonging to ten bacterial classes and 28 bacterial families. Bacterial species identified as carrying genes encoding starch binding modules increased in abundance in response to starches. The use of hybrid metagenomics has allowed identification of several uncultured species with the functional potential to degrade starch substrates for future study.

Longitudinal hybrid metagenomic analyses of a human stool sample reveal compositional and functional variation in response to six structurally-distinct carbohydrates, providing insight into how gut bacteria utilize various carbohydrate sources.

Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes

Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.

Metagenomic binning with assembly graph embeddings

MOTIVATION

Despite recent advancements in sequencing technologies and assembly methods, obtaining high-quality microbial genomes from metagenomic samples is still not a trivial task. Current metagenomic binners do not take full advantage of assembly graphs and are not optimized for long-read assemblies. Deep graph learning algorithms have been proposed in other fields to deal with complex graph data structures. The graph structure generated during the assembly process could be integrated with contig features to obtain better bins with deep learning.

RESULTS

We propose GraphMB, which uses graph neural networks to incorporate the assembly graph into the binning process. We test GraphMB on long-read datasets of different complexities, and compare the performance with other binners in terms of the number of High Quality (HQ) genome bins obtained. With our approach, we were able to obtain unique bins on all real datasets, and obtain more bins on most datasets. In particular, we obtained on average 17.5% more HQ bins when compared with state-of-the-art binners and 13.7% when aggregating the results of our binner with the others. These results indicate that a deep learning model can integrate contig-specific and graph-structure information to improve metagenomic binning.

AVAILABILITY AND IMPLEMENTATION

GraphMB is available from https://github.com/MicrobialDarkMatter/GraphMB.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrating Genome-Resolved Metagenomics with Trait-Based Process Modeling to Determine Biokinetics of Distinct Nitrifying Communities within Activated Sludge

Conventional bioprocess models for wastewater treatment are based on aggregated bulk biomass concentrations and do not incorporate microbial physiological diversity. Such a broad aggregation of microbial functional groups can fail to predict ecosystem dynamics when high levels of physiological diversity exist within trophic guilds. For instance, functional diversity among nitrite-oxidizing bacteria (NOB) can obfuscate engineering strategies for their out-selection in activated sludge (AS), which is desirable to promote energy-efficient nitrogen removal. Here, we hypothesized that different NOB populations within AS can have different physiological traits that drive process performance, which we tested by estimating biokinetic growth parameters using a combination of highly replicated respirometry, genome-resolved metagenomics, and process modeling. A lab-scale AS reactor subjected to a selective pressure for over 90 days experienced resilience of NOB activity. We recovered three coexisting Nitrospira population genomes belonging to two sublineages, which exhibited distinct growth strategies and underwent a compositional shift following the selective pressure. A trait-based process model calibrated at the NOB genus level better predicted nitrite accumulation than a conventional process model calibrated at the NOB guild level. This work demonstrates that trait-based modeling can be leveraged to improve our prediction, control, and design of functionally diverse microbiomes driving key environmental biotechnologies.

Fermentation liquid as a carbon source for wastewater nitrogen removal reduced nitrogenous disinfection byproduct formation potentials of the effluent

Sludge alkaline fermentation liquid (SAFL) is an alternative to sodium acetate (NaAc) in enhancing wastewater nitrogen removal. Upon SAFL addition, dissolved organic nitrogen (DON) can be externally introduced or biologically synthesized during nitrogen removal, which is an important precursor to toxic nitrogenous disinfection by-products (N-DBPs). This study aims to evaluate the effects of different carbon source addition on effluent DON concentration, composition, and N-DBP formation potentials. A lab-scale A²O system treating real municipal wastewater was operated with NaAc or SAFL as external carbon sources. DON molecules and potential N-DBP precursors were identified by Orbitrap mass spectrometry. Subsequently, major microorganisms contributing to DON biosynthesis were suggested based on metagenomics. It was found that effluent DON was higher with SAFL as the carbon source than NaAc (1.51 ± 0.24 v.s. 0.56 ± 0.08 mg N/L, p < 0.05). Nevertheless, dichloroacetonitrile and nitrosamine formation potentials (7.14 ± 1.02 and 1.57 ± 0.07 μg/mg DON-N, respectively) of the effluent with SAFL addition were 42.79 ± 2.42% and 54.89 ± 1.70% lower than those of NaAc. Protein- and lignin-like compounds were the most abundant DON molecules in the effluent, where alanine, glycine and tyrosine were important precursors to N-DBPs. Azonexus and Flavobacterium spp. were positively correlated with these precursors, and possessed key genes involved in precursor synthesis. SAFL is a promising carbon source, not only for achieving efficient inorganic nitrogen and DON removals, but also for reducing N-DBP formation potentials of chlorinated effluent.