A genomic catalog of Earth's microbiomes

A mixed blessing of influent leachate microbes in downstream biotreatment systems of a full-scale landfill leachate treatment plant

In landfill leachate treatment plants (LLTPs), the microbiome plays a pivotal role in the decomposition of organic compounds, reduction in nutrient levels, and elimination of toxins. However, the effects of microbes in landfill leachate influents on downstream treatment systems remain poorly understood. To address this knowledge gap, we collected 23 metagenomic and 12 metatranscriptomic samples from landfill leachate and activated sludge from various treatment units in a full-scale LLTP. We successfully recovered 1,152 non-redundant metagenome-assembled genomes (MAGs), encompassing a wide taxonomic range, including 48 phyla, 95 classes, 166 orders, 247 families, 238 genera, and 1,152 species. More diverse microbes were observed in the influent leachate than in the downstream biotreatment systems, among which, an unprecedented ∼30 % of microbes with transcriptional expression migrated from the influent to the biological treatment units. Network analysis revealed that 399 shared MAGs across the four units exhibited high node centrality and degree, thus supporting enhanced interactions and increased stability of microbial communities. Functional reconstruction and genome characterization of MAGs indicated that these shared MAGs possessed greater capabilities for carbon, nitrogen, sulfur, and arsenic metabolism compared to non-shared MAGs. We further identified a novel species of Zixibacteria in the leachate influent with discrete lineages from those in other environments that accounted for up to 17 % of the abundance of the shared microbial community and exhibited notable metabolic versatility. Meanwhile, we presented groundbreaking evidence of the involvement of Zixibacteria-encoded genes in the production of harmful gas emissions, such as N2O and H2S, at the transcriptional level, thus suggesting that influent microbes may pose safety risks to downstream treatment systems. In summary, this study revealed the complex impact of the influent microbiome on LLTP and emphasizes the need to consider these microbial characteristics when designing treatment technologies and strategies for landfill leachate management.

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

Biosynthetic potential of the sediment microbial subcommunities of an unexplored karst ecosystem and its ecological implications

Microbial communities from various environments have been studied in the quest for new natural products with a broad range of applications in medicine and biotechnology. We employed an enrichment method and genome mining tools to examine the biosynthetic potential of microbial communities in the sediments of a coastal sinkhole within the karst ecosystem of the Yucatán Peninsula, Mexico. Our investigation led to the detection of 203 biosynthetic gene clusters (BGCs) and 55 secondary metabolites (SMs) within 35 high-quality metagenome-assembled genomes (MAGs) derived from these subcommunities. The most abundant types of BGCs were Terpene, Nonribosomal peptide-synthetase, and Type III polyketide synthase. Some of the in silico identified BGCs and SMs have been previously reported to exhibit biological activities against pathogenic bacteria and fungi. Others could play significant roles in the sinkhole ecosystem, such as iron solubilization and osmotic stress protection. Interestingly, 75% of the BGCs showed no sequence homology with bacterial BGCs previously reported in the MiBIG database. This suggests that the microbial communities in this environment could be an untapped source of genes encoding novel specialized compounds. The majority of the BGCs were identified in pathways found in the genus Virgibacillus, followed by Sporosarcina, Siminovitchia, Rhodococcus, and Halomonas. The latter, along with Paraclostridium and Lysinibacillus, had the highest number of identified BGC types. This study offers fresh insights into the potential ecological role of SMs from sediment microbial communities in an unexplored environment, underscoring their value as a source of novel natural products.

An in-depth evaluation of metagenomic classifiers for soil microbiomes

BACKGROUND

Recent endeavours in metagenomics, exemplified by projects such as the human microbiome project and TARA Oceans, have illuminated the complexities of microbial biomes. A robust bioinformatic pipeline and meticulous evaluation of their methodology have contributed to the success of these projects. The soil environment, however, with its unique challenges, requires a specialized methodological exploration to maximize microbial insights. A notable limitation in soil microbiome studies is the dearth of soil-specific reference databases available to classifiers that emulate the complexity of soil communities. There is also a lack of in-vitro mock communities derived from soil strains that can be assessed for taxonomic classification accuracy.

RESULTS

In this study, we generated a custom in-silico mock community containing microbial genomes commonly observed in the soil microbiome. Using this mock community, we simulated shotgun sequencing data to evaluate the performance of three leading metagenomic classifiers: Kraken2 (supplemented with Bracken, using a custom database derived from GTDB-TK genomes along with its own default database), Kaiju, and MetaPhlAn, utilizing their respective default databases for a robust analysis. Our results highlight the importance of optimizing taxonomic classification parameters, database selection, as well as analysing trimmed reads and contigs. Our study showed that classifiers tailored to the specific taxa present in our samples led to fewer errors compared to broader databases including microbial eukaryotes, protozoa, or human genomes, highlighting the effectiveness of targeted taxonomic classification. Notably, an optimal classifier performance was achieved when applying a relative abundance threshold of 0.001% or 0.005%. The Kraken2 supplemented with bracken, with a custom database demonstrated superior precision, sensitivity, F1 score, and overall sequence classification. Using a custom database, this classifier classified 99% of in-silico reads and 58% of real-world soil shotgun reads, with the latter identifying previously overlooked phyla using a custom database.

CONCLUSION

This study underscores the potential advantages of in-silico methodological optimization in metagenomic analyses, especially when deciphering the complexities of soil microbiomes. We demonstrate that the choice of classifier and database significantly impacts microbial taxonomic profiling. Our findings suggest that employing Kraken2 with Bracken, coupled with a custom database of GTDB-TK genomes and fungal genomes at a relative abundance threshold of 0.001% provides optimal accuracy in soil shotgun metagenome analysis.

Sequence-structure-function characterization of the emerging tetracycline destructase family of antibiotic resistance enzymes

Tetracycline destructases (TDases) are flavin monooxygenases which can confer resistance to all generations of tetracycline antibiotics. The recent increase in the number and diversity of reported TDase sequences enables a deep investigation of the TDase sequence-structure-function landscape. Here, we evaluate the sequence determinants of TDase function through two complementary approaches: (1) constructing profile hidden Markov models to predict new TDases, and (2) using multiple sequence alignments to identify conserved positions important to protein function. Using the HMM-based approach we screened 50 high-scoring candidate sequences in Escherichia coli, leading to the discovery of 13 new TDases. The X-ray crystal structures of two new enzymes from Legionella species were determined, and the ability of anhydrotetracycline to inhibit their tetracycline-inactivating activity was confirmed. Using the MSA-based approach we identified 31 amino acid positions 100% conserved across all known TDase sequences. The roles of these positions were analyzed by alanine-scanning mutagenesis in two TDases, to study the impact on cell and in vitro activity, structure, and stability. These results expand the diversity of TDase sequences and provide valuable insights into the roles of important residues in TDases, and flavin monooxygenases more broadly.

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.

The defensome of complex bacterial communities

Bacteria have developed various defense mechanisms to avoid infection and killing in response to the fast evolution and turnover of viruses and other genetic parasites. Such pan-immune system (defensome) encompasses a growing number of defense lines that include well-studied innate and adaptive systems such as restriction-modification, CRISPR-Cas and abortive infection, but also newly found ones whose mechanisms are still poorly understood. While the abundance and distribution of defense systems is well-known in complete and culturable genomes, there is a void in our understanding of their diversity and richness in complex microbial communities. Here we performed a large-scale in-depth analysis of the defensomes of 7759 high-quality bacterial population genomes reconstructed from soil, marine, and human gut environments. We observed a wide variation in the frequency and nature of the defensome among large phyla, which correlated with lifestyle, genome size, habitat, and geographic background. The defensome's genetic mobility, its clustering in defense islands, and genetic variability was found to be system-specific and shaped by the bacterial environment. Hence, our results provide a detailed picture of the multiple immune barriers present in environmentally distinct bacterial communities and set the stage for subsequent identification of novel and ingenious strategies of diversification among uncultivated microbes.

Extracellular vesicle formation in Euryarchaeota is driven by a small GTPase

Since their discovery, extracellular vesicles (EVs) have changed our view on how organisms interact with their extracellular world. EVs are able to traffic a diverse array of molecules across different species and even domains, facilitating numerous functions. In this study, we investigate EV production in Euryarchaeota, using the model organism Haloferax volcanii. We uncover that EVs enclose RNA, with specific transcripts preferentially enriched, including those with regulatory potential, and conclude that EVs can act as an RNA communication system between haloarchaea. We demonstrate the key role of an EV-associated small GTPase for EV formation in H. volcanii that is also present across other diverse evolutionary branches of Archaea. We propose the name, ArvA, for the identified family of archaeal vesiculating GTPases. Additionally, we show that two genes in the same operon with arvA (arvB and arvC) are also involved in EV formation. Both, arvB and arvC, are closely associated with arvA in the majority of other archaea encoding ArvA. Our work demonstrates that small GTPases involved in membrane deformation and vesiculation, ubiquitous in Eukaryotes, are also present in Archaea and are widely distributed across diverse archaeal phyla.

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

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

IMPORTANCE

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

An optimized culturomics strategy for isolation of human milk microbiota

Viable microorganisms and a diverse microbial ecosystem found in human milk play a crucial role in promoting healthy immune system and shaping the microbial community in the infant's gut. Culturomics is a method to obtain a comprehensive repertoire of human milk microbiota. However, culturomics is an onerous procedure, and needs expertise, making it difficult to be widely implemented. Currently, there is no efficient and feasible culturomics method specifically designed for human milk microbiota yet. Therefore, the aim of this study was to develop a more efficient and feasible culturomics method specifically designed for human milk microbiota. We obtained fresh samples of human milk from healthy Chinese mothers and conducted a 27-day enrichment process using blood culture bottles. Bacterial isolates were harvested at different time intervals and cultured on four different types of media. Using matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) analysis, we identified a total of 6601 colonies and successfully obtained 865 strains, representing 4 phyla, 21 genera, and 54 species. By combining CBA and MRS media, we were able to cultivate over 94.4% of bacterial species with high diversity, including species-specific microorganisms. Prolonged pre-incubation in blood culture bottles significantly increased the number of bacterial species by about 33% and improved the isolation efficiency of beneficial bacteria with low abundance in human milk. After optimization, we reduced the pre-incubation time in blood culture bottles and selected optimal picking time-points (0, 3, and 6 days) at 37°C. By testing 6601 colonies using MALDI-TOF MS, we estimated that this new protocol could obtain more than 90% of bacterial species, reducing the workload by 57.0%. In conclusion, our new culturomics strategy, which involves the combination of CBA and MRS media, extended pre-incubation enrichment, and optimized picking time-points, is a feasible method for studying the human milk microbiota. This protocol significantly improves the efficiency of culturomics and allows for the establishment of a comprehensive repertoire of bacterial species and strains in human milk.

Why and how to use the SeqCode

The SeqCode, formally called the Code of Nomenclature of Prokaryotes Described from Sequence Data, is a new code of nomenclature in which genome sequences are the nomenclatural types for the names of prokaryotic species. While similar to the International Code of Nomenclature of Prokaryotes (ICNP) in structure and rules of priority, it does not require the deposition of type strains in international culture collections. Thus, it allows for the formation of permanent names for uncultured prokaryotes whose nearly complete genome sequences have been obtained directly from environmental DNA as well as other prokaryotes that cannot be deposited in culture collections. Because the diversity of uncultured prokaryotes greatly exceeds that of readily culturable prokaryotes, the SeqCode is the only code suitable for naming the majority of prokaryotic species. The start date of the SeqCode was January 1, 2022, and the online Registry (https://seqco.de/) was created to ensure valid publication of names. The SeqCode recognizes all names validly published under the ICNP before 2022. After that date, names validly published under the SeqCode compete with ICNP names for priority. As a result, species can have only one name, either from the SeqCode or ICNP, enabling effective communication and the creation of unified taxonomies of uncultured and cultured prokaryotes. The SeqCode is administered by the SeqCode Committee, which is comprised of the SeqCode Community and elected administrative components. Anyone with an interest in the systematics of prokaryotes is encouraged to join the SeqCode Community and participate in the development of this resource.

A genome and gene catalog of the aquatic microbiomes of the Tibetan Plateau

The Tibetan Plateau supplies water to nearly 2 billion people in Asia, but climate change poses threats to its aquatic microbial resources. Here, we construct the Tibetan Plateau Microbial Catalog by sequencing 498 metagenomes from six water ecosystems (saline lakes, freshwater lakes, rivers, hot springs, wetlands and glaciers). Our catalog expands knowledge of regional genomic diversity by presenting 32,355 metagenome-assembled genomes that de-replicated into 10,723 representative genome-based species, of which 88% were unannotated. The catalog contains nearly 300 million non-redundant gene clusters, of which 15% novel, and 73,864 biosynthetic gene clusters, of which 50% novel, thus expanding known functional diversity. Using these data, we investigate the Tibetan Plateau aquatic microbiome's biogeography along a distance of 2,500 km and >5 km in altitude. Microbial compositional similarity and the shared gene count with the Tibetan Plateau microbiome decline along with distance and altitude difference, suggesting a dispersal pattern. The Tibetan Plateau Microbial Catalog stands as a substantial repository for high-altitude aquatic microbiome resources, providing potential for discovering novel lineages and functions, and bridging knowledge gaps in microbiome biogeography.

A metagenomic alpha-diversity index for microbial functional biodiversity

Alpha-diversity indices are an essential tool for describing and comparing biodiversity. Microbial ecologists apply indices originally intended for, or adopted by, macroecology to address questions relating to taxonomy (conserved marker) and function (metagenome-based data). In this Perspective piece, I begin by discussing the nature and mathematical quirks important for interpreting routinely employed alpha-diversity indices. Secondly, I propose a metagenomic alpha-diversity index (MD) that measures the (dis)similarity of protein-encoding genes within a community. MD has defined limits, whereby a community comprised mostly of similar, poorly diverse protein-encoding genes pulls the index to the lower limit, while a community rich in divergent homologs and unique genes drives it toward the upper limit. With data acquired from an in silico and three in situ metagenome studies, I derive MD and typical alpha-diversity indices applied to taxonomic (ribosomal rRNA) and functional (all protein-encoding) genes, and discuss their relationships with each other. Not all alpha-diversity indices detect biological trends, and taxonomic does not necessarily follow functional biodiversity. Throughout, I explain that protein Richness and MD provide complementary and easily interpreted information, while probability-based indices do not. Finally, considerations regarding the unique nature of microbial metagenomic data and its relevance for describing functional biodiversity are discussed.

Metagenomic analysis of carbohydrate-active enzymes and their contribution to marine sediment biodiversity

Marine sediments constitute the world's most substantial long-term carbon repository. The microorganisms dwelling in these sediments mediate the transformation of fixed oceanic carbon, but their contribution to the carbon cycle is not fully understood. Previous culture-independent investigations into sedimentary microorganisms have underscored the significance of carbohydrates in the carbon cycle. In this study, we employ a metagenomic methodology to investigate the distribution and abundance of carbohydrate-active enzymes (CAZymes) in 37 marine sediments sites. These sediments exhibit varying oxygen availability and were isolated in diverse regions worldwide. Our comparative analysis is based on the metabolic potential for oxygen utilisation, derived from genes present in both oxic and anoxic environments. We found that extracellular CAZyme modules targeting the degradation of plant and algal detritus, necromass, and host glycans were abundant across all metagenomic samples. The analysis of these results indicates that the oxic/anoxic conditions not only influence the taxonomic composition of the microbial communities, but also affect the occurrence of CAZyme modules involved in the transformation of necromass, algae and plant detritus. To gain insight into the sediment microbial taxa, we reconstructed metagenome assembled genomes (MAG) and examined the presence of primary extracellular carbohydrate active enzyme (CAZyme) modules. Our findings reveal that the primary CAZyme modules and the CAZyme gene clusters discovered in our metagenomes were prevalent in the Bacteroidia, Gammaproteobacteria, and Alphaproteobacteria classes. We compared those MAGs to organisms from the same taxonomic classes found in soil, and we found that they were similar in its CAZyme repertoire, but the soil MAG contained a more abundant and diverse CAZyme content. Furthermore, the data indicate that abundant classes in our metagenomic samples, namely Alphaproteobacteria, Bacteroidia and Gammaproteobacteria, play a pivotal role in carbohydrate transformation within the initial few metres of the sediments.

Hydrolytic endonucleolytic ribozyme (HYER) is programmable for sequence-specific DNA cleavage

Ribozymes are catalytic RNAs with diverse functions including self-splicing and polymerization. This work aims to discover natural ribozymes that behave as hydrolytic and sequence-specific DNA endonucleases, which could be repurposed as DNA manipulation tools. Focused on bacterial group II-C introns, we found that many systems without intron-encoded protein propagate multiple copies in their resident genomes. These introns, named HYdrolytic Endonucleolytic Ribozymes (HYERs), cleaved RNA, single-stranded DNA, bubbled double-stranded DNA (dsDNA), and plasmids in vitro. HYER1 generated dsDNA breaks in the mammalian genome. Cryo-electron microscopy analysis revealed a homodimer structure for HYER1, where each monomer contains a Mg2+-dependent hydrolysis pocket and captures DNA complementary to the target recognition site (TRS). Rational designs including TRS extension, recruiting sequence insertion, and heterodimerization yielded engineered HYERs showing improved specificity and flexibility for DNA manipulation.

Tools for microbial single-cell genomics for obtaining uncultured microbial genomes

The advent of next-generation sequencing technologies has facilitated the acquisition of large amounts of DNA sequence data at a relatively low cost, leading to numerous breakthroughs in decoding microbial genomes. Among the various genome sequencing activities, metagenomic analysis, which entails the direct analysis of uncultured microbial DNA, has had a profound impact on microbiome research and has emerged as an indispensable technology in this field. Despite its valuable contributions, metagenomic analysis is a "bulk analysis" technique that analyzes samples containing a wide diversity of microbes, such as bacteria, yielding information that is averaged across the entire microbial population. In order to gain a deeper understanding of the heterogeneous nature of the microbial world, there is a growing need for single-cell analysis, similar to its use in human cell biology. With this paradigm shift in mind, comprehensive single-cell genomics technology has become a much-anticipated innovation that is now poised to revolutionize microbiome research. It has the potential to enable the discovery of differences at the strain level and to facilitate a more comprehensive examination of microbial ecosystems. In this review, we summarize the current state-of-the-art in microbial single-cell genomics, highlighting the potential impact of this technology on our understanding of the microbial world. The successful implementation of this technology is expected to have a profound impact in the field, leading to new discoveries and insights into the diversity and evolution of microbes.

A genomic view of environmental and life history controls on microbial nitrogen acquisition strategies

Microorganisms have evolved diverse strategies to acquire the vital element nitrogen (N) from the environment. Ecological and physiological controls on the distribution of these strategies among microbes remain unclear. In this study, we examine the distribution of 10 major N acquisition strategies in taxonomically and metabolically diverse microbial genomes, including those from the Genomic Catalogue of Earth's Microbiomes dataset. We utilize a marker gene-based approach to assess relationships between N acquisition strategy prevalence and microbial life history strategies. Our results underscore energetic costs of assimilation as a broad control on strategy distribution. The most prevalent strategies are the uptake of ammonium and simple amino acids, which have relatively low energetic costs, while energy-intensive biological nitrogen fixation is the least common. Deviations from the energy-based framework include the higher-than-expected prevalence of the assimilatory pathway for chitin, a large organic polymer. Energy availability is also important, with aerobic chemoorganotrophs  and oxygenic phototrophs notably possessing ~2-fold higher numbers of total strategies compared to anaerobic microbes. Environmental controls are evidenced by the enrichment of inorganic N assimilation strategies among free-living taxa compared to host-associated taxa. Physiological constraints such as pathway incompatibility add complexity to N acquisition strategy distributions. Finally, we discuss the necessity for microbially-relevant spatiotemporal environmental metadata for improving mechanistic and prediction-oriented analyses of genomic data.

Functional and evolutionary significance of unknown genes from uncultivated taxa

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

Recent advances in the culture-independent discovery of natural products using metagenomic approaches

Natural products derived from bacterial sources have long been pivotal in the discovery of drug leads. However, the cultivation of only about 1% of bacteria in laboratory settings has left a significant portion of biosynthetic diversity hidden within the genomes of uncultured bacteria. Advances in sequencing technologies now enable the exploration of genetic material from these metagenomes through culture-independent methods. This approach involves extracting genetic sequences from environmental DNA and applying a hybrid methodology that combines functional screening, sequence tag-based homology screening, and bioinformatic-assisted chemical synthesis. Through this process, numerous valuable natural products have been identified and synthesized from previously uncharted metagenomic territories. This paper provides an overview of the recent advancements in the utilization of culture-independent techniques for the discovery of novel biosynthetic gene clusters and bioactive small molecules within metagenomic libraries.

scMicrobe PTA: Near Complete Genomes from Single Bacterial Cells

Microbial genomes produced by single-cell amplification are largely incomplete. Here, we show that primary template amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard amplification approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

Structure-guided discovery of anti-CRISPR and anti-phage defense proteins

Bacteria use a variety of defense systems to protect themselves from phage infection. In turn, phages have evolved diverse counter-defense measures to overcome host defenses. Here, we use protein structural similarity and gene co-occurrence analyses to screen >66 million viral protein sequences and >330,000 metagenome-assembled genomes for the identification of anti-phage and counter-defense systems. We predict structures for ~300,000 proteins and perform large-scale, pairwise comparison to known anti-CRISPR (Acr) and anti-phage proteins to identify structural homologs that otherwise may not be uncovered using primary sequence search. This way, we identify a Bacteroidota phage Acr protein that inhibits Cas12a, and an Akkermansia muciniphila anti-phage defense protein, termed BxaP. Gene bxaP is found in loci encoding Bacteriophage Exclusion (BREX) and restriction-modification defense systems, but confers immunity independently. Our work highlights the advantage of combining protein structural features and gene co-localization information in studying host-phage interactions.

Microbial metabolites as modulators of host physiology

The gut microbiota is increasingly recognised as a key player in influencing human health and changes in the gut microbiota have been strongly linked with many non-communicable conditions in humans such as type 2 diabetes, obesity and cardiovascular disease. However, characterising the molecular mechanisms that underpin these associations remains an important challenge for researchers. The gut microbiota is a complex microbial community that acts as a metabolic interface to transform ingested food (and other xenobiotics) into metabolites that are detected in the host faeces, urine and blood. Many of these metabolites are only produced by microbes and there is accumulating evidence to suggest that these microbe-specific metabolites do act as effectors to influence human physiology. For example, the gut microbiota can digest dietary complex polysaccharides (such as fibre) into short-chain fatty acids (SCFA) such as acetate, propionate and butyrate that have a pervasive role in host physiology from nutrition to immune function. In this review we will outline our current understanding of the role of some key microbial metabolites, such as SCFA, indole and bile acids, in human health. Whilst many studies linking microbial metabolites with human health are correlative we will try to highlight examples where genetic evidence is available to support a specific role for a microbial metabolite in host health and well-being.

Genomic profiles of four novel cyanobacteria MAGs from Lake Vanda, Antarctica: insights into photosynthesis, cold tolerance, and the circadian clock

Cyanobacteria in polar environments face environmental challenges, including cold temperatures and extreme light seasonality with small diurnal variation, which has implications for polar circadian clocks. However, polar cyanobacteria remain underrepresented in available genomic data, and there are limited opportunities to study their genetic adaptations to these challenges. This paper presents four new Antarctic cyanobacteria metagenome-assembled genomes (MAGs) from microbial mats in Lake Vanda in the McMurdo Dry Valleys in Antarctica. The four MAGs were classified as Leptolyngbya sp. BulkMat.35, Pseudanabaenaceae cyanobacterium MP8IB2.15, Microcoleus sp. MP8IB2.171, and Leptolyngbyaceae cyanobacterium MP9P1.79. The MAGs contain 2.76 Mbp - 6.07 Mbp, and the bin completion ranges from 74.2-92.57%. Furthermore, the four cyanobacteria MAGs have average nucleotide identities (ANIs) under 90% with each other and under 77% with six existing polar cyanobacteria MAGs and genomes. This suggests that they are novel cyanobacteria and demonstrates that polar cyanobacteria genomes are underrepresented in reference databases and there is continued need for genome sequencing of polar cyanobacteria. Analyses of the four novel and six existing polar cyanobacteria MAGs and genomes demonstrate they have genes coding for various cold tolerance mechanisms and most standard circadian rhythm genes with the Leptolyngbya sp. BulkMat.35 and Leptolyngbyaceae cyanobacterium MP9P1.79 contained kaiB3, a divergent homolog of kaiB.

Precision enzyme discovery through targeted mining of metagenomic data

Metagenomics has opened new avenues for exploring the genetic potential of uncultured microorganisms, which may serve as promising sources of enzymes and natural products for industrial applications. Identifying enzymes with improved catalytic properties from the vast amount of available metagenomic data poses a significant challenge that demands the development of novel computational and functional screening tools. The catalytic properties of all enzymes are primarily dictated by their structures, which are predominantly determined by their amino acid sequences. However, this aspect has not been fully considered in the enzyme bioprospecting processes. With the accumulating number of available enzyme sequences and the increasing demand for discovering novel biocatalysts, structural and functional modeling can be employed to identify potential enzymes with novel catalytic properties. Recent efforts to discover new polysaccharide-degrading enzymes from rumen metagenome data using homology-based searches and machine learning-based models have shown significant promise. Here, we will explore various computational approaches that can be employed to screen and shortlist metagenome-derived enzymes as potential biocatalyst candidates, in conjunction with the wet lab analytical methods traditionally used for enzyme characterization.

SeqCode in the golden age of prokaryotic systematics

The SeqCode is a new code of prokaryotic nomenclature that was developed to validate taxon names using genome sequences as the type material. The present article provides an independent view about the SeqCode, highlighting its history, current status, basic features, pros and cons, and use to date. We also discuss important topics to consider for validation of novel prokaryotic taxon names using genomes as the type material. Owing to significant advances in metagenomics and cultivation methods, hundreds of novel prokaryotic species are expected to be discovered in the coming years. This manuscript aims to stimulate and enrich the debate around the use of the SeqCode in the upcoming golden age of prokaryotic taxon discovery and systematics.

Sheaths are diverse and abundant cell surface layers in archaea

Prokaryotic cells employ multiple protective layers crucial for defense, structural integrity, and cellular interactions in the environment. Archaea often feature an S-layer, with some species possessing additional and remarkably resistant sheaths. The archaeal sheath has been studied in Methanothrix and Methanospirillum, revealing a complex structure consisting of amyloid proteins organized into rings. Here, we conducted a comprehensive survey of sheath-forming proteins (SH proteins) across archaeal genomes. Structural modeling reveals a rich diversity of SH proteins, indicating the presence of a sheath in members of the TACK superphylum (Thermoprotei), as well as in the methanotrophic ANME-1. SH proteins are present in up to 40 copies per genome and display diverse domain arrangements suggesting multifunctional roles within the sheath, and potential involvement in cell-cell interaction with syntrophic partners. We uncover a complex evolutionary dynamic, indicating active exchange of SH proteins in archaeal communities. We find that viruses infecting sheathed archaea encode a diversity of SH-like proteins and we use them as markers to identify 580 vOTUs potentially associated with sheathed archaea. Structural modeling suggests that viral SH proteins can form complexes with the host SH proteins. We propose a previously unreported egress strategy where the expression of viral SH-like proteins may disrupt the integrity of the host sheath and facilitate viral exit during lysis. Together, our results significantly expand knowledge of the diversity and evolution of the archaeal sheath, which has been largely understudied but might have an important role in shaping microbial communities.

Globally distributed bacteriophage genomes reveal mechanisms of tripartite phage-bacteria-coral interactions

Reef-building corals depend on an intricate community of microorganisms for functioning and resilience. The infection of coral-associated bacteria by bacteriophages can modify bacterial ecological interactions, yet very little is known about phage functions in the holobiont. This gap stems from methodological limitations that have prevented the recovery of high-quality viral genomes and bacterial host assignment from coral samples. Here, we introduce a size fractionation approach that increased bacterial and viral recovery in coral metagenomes by 9-fold and 2-fold, respectively, and enabled the assembly and binning of bacterial and viral genomes at relatively low sequencing coverage. We combined these viral genomes with those derived from 677 publicly available metagenomes, viromes, and bacterial isolates from stony corals to build a global coral virus database of over 20,000 viral genomic sequences spanning four viral realms. The tailed bacteriophage families Kyanoviridae and Autographiviridae were the most abundant, replacing groups formerly referred to as Myoviridae and Podoviridae, respectively. Prophage and CRISPR spacer linkages between these viruses and 626 bacterial metagenome-assembled genomes and bacterial isolates showed that most viruses infected Alphaproteobacteria, the most abundant class, and less abundant taxa like Halanaerobiia and Bacteroidia. A host-phage-gene network identified keystone viruses with the genomic capacity to modulate bacterial metabolic pathways and direct molecular interactions with eukaryotic cells. This study reveals the genomic basis of nested symbioses between bacteriophage, bacteria, and the coral host and its endosymbiotic algae.

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

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

Genomic basis of environmental adaptation in the widespread poly-extremophilic Exiguobacterium group

Delineating cohesive ecological units and determining the genetic basis for their environmental adaptation are among the most important objectives in microbiology. In the last decade, many studies have been devoted to characterizing the genetic diversity in microbial populations to address these issues. However, the impact of extreme environmental conditions, such as temperature and salinity, on microbial ecology and evolution remains unclear so far. In order to better understand the mechanisms of adaptation, we studied the (pan)genome of Exiguobacterium, a poly-extremophile bacterium able to grow in a wide range of environments, from permafrost to hot springs. To have the genome for all known Exiguobacterium type strains, we first sequenced those that were not yet available. Using a reverse-ecology approach, we showed how the integration of phylogenomic information, genomic features, gene and pathway enrichment data, regulatory element analyses, protein amino acid composition, and protein structure analyses of the entire Exiguobacterium pangenome allows to sharply delineate ecological units consisting of mesophilic, psychrophilic, halophilic-mesophilic, and halophilic-thermophilic ecotypes. This in-depth study clarified the genetic basis of the defined ecotypes and identified some key mechanisms driving the environmental adaptation to extreme environments. Our study points the way to organizing the vast microbial diversity into meaningful ecologically units, which, in turn, provides insight into how microbial communities adapt and respond to different environmental conditions in a changing world.

IMG/PR: a database of plasmids from genomes and metagenomes with rich annotations and metadata

Plasmids are mobile genetic elements found in many clades of Archaea and Bacteria. They drive horizontal gene transfer, impacting ecological and evolutionary processes within microbial communities, and hold substantial importance in human health and biotechnology. To support plasmid research and provide scientists with data of an unprecedented diversity of plasmid sequences, we introduce the IMG/PR database, a new resource encompassing 699 973 plasmid sequences derived from genomes, metagenomes and metatranscriptomes. IMG/PR is the first database to provide data of plasmid that were systematically identified from diverse microbiome samples. IMG/PR plasmids are associated with rich metadata that includes geographical and ecosystem information, host taxonomy, similarity to other plasmids, functional annotation, presence of genes involved in conjugation and antibiotic resistance. The database offers diverse methods for exploring its extensive plasmid collection, enabling users to navigate plasmids through metadata-centric queries, plasmid comparisons and BLAST searches. The web interface for IMG/PR is accessible at https://img.jgi.doe.gov/pr. Plasmid metadata and sequences can be downloaded from https://genome.jgi.doe.gov/portal/IMG_PR.

SPIRE: a Searchable, Planetary-scale mIcrobiome REsource

Meta'omic data on microbial diversity and function accrue exponentially in public repositories, but derived information is often siloed according to data type, study or sampled microbial environment. Here we present SPIRE, a Searchable Planetary-scale mIcrobiome REsource that integrates various consistently processed metagenome-derived microbial data modalities across habitats, geography and phylogeny. SPIRE encompasses 99 146 metagenomic samples from 739 studies covering a wide array of microbial environments and augmented with manually-curated contextual data. Across a total metagenomic assembly of 16 Tbp, SPIRE comprises 35 billion predicted protein sequences and 1.16 million newly constructed metagenome-assembled genomes (MAGs) of medium or high quality. Beyond mapping to the high-quality genome reference provided by proGenomes3 (http://progenomes.embl.de), these novel MAGs form 92 134 novel species-level clusters, the majority of which are unclassified at species level using current tools. SPIRE enables taxonomic profiling of these species clusters via an updated, custom mOTUs database (https://motu-tool.org/) and includes several layers of functional annotation, as well as crosslinks to several (micro-)biological databases. The resource is accessible, searchable and browsable via http://spire.embl.de.

ABC-HuMi: the Atlas of Biosynthetic Gene Clusters in the Human Microbiome

The human microbiome has emerged as a rich source of diverse and bioactive natural products, harboring immense potential for therapeutic applications. To facilitate systematic exploration and analysis of its biosynthetic landscape, we present ABC-HuMi: the Atlas of Biosynthetic Gene Clusters (BGCs) in the Human Microbiome. ABC-HuMi integrates data from major human microbiome sequence databases and provides an expansive repository of BGCs compared to the limited coverage offered by existing resources. Employing state-of-the-art BGC prediction and analysis tools, our database ensures accurate annotation and enhanced prediction capabilities. ABC-HuMi empowers researchers with advanced browsing, filtering, and search functionality, enabling efficient exploration of the resource. At present, ABC-HuMi boasts a catalog of 19 218 representative BGCs derived from the human gut, oral, skin, respiratory and urogenital systems. By capturing the intricate biosynthetic potential across diverse human body sites, our database fosters profound insights into the molecular repertoire encoded within the human microbiome and offers a comprehensive resource for the discovery and characterization of novel bioactive compounds. The database is freely accessible at https://www.ccb.uni-saarland.de/abc_humi/.

Challenges in estimating effective population sizes from metagenome-assembled genomes

Effective population size (Ne) plays a critical role in shaping the relative efficiency between natural selection and genetic drift, thereby serving as a cornerstone for understanding microbial ecological dynamics. Direct Ne estimation relies on neutral genetic diversity within closely related genomes, which is, however, often constrained by the culturing difficulties for the vast majority of prokaryotic lineages. Metagenome-assembled genomes (MAGs) offer a high-throughput alternative for genomic data acquisition, yet their accuracy in Ne estimation has not been fully verified. This study examines the Thermococcus genus, comprising 66 isolated strains and 29 MAGs, to evaluate the reliability of MAGs in Ne estimation. Despite the even distribution across the Thermococcus phylogeny and the comparable internal average nucleotide identity (ANI) between isolate populations and MAG populations, our results reveal consistently lower Ne estimates from MAG populations. This trend of underestimation is also observed in various MAG populations across three other bacterial genera. The underrepresentation of genetic variation in MAGs, including loss of allele frequency data and variable genomic segments, likely contributes to the underestimation of Ne. Our findings underscore the necessity for caution when employing MAGs for evolutionary studies, which often depend on high-quality genome assemblies and nucleotide-level diversity.

CoCoPyE: feature engineering for learning and prediction of genome quality indices

BACKGROUND

The exploration of the microbial world has been greatly advanced by the reconstruction of genomes from metagenomic sequence data. However, the rapidly increasing number of metagenome-assembled genomes has also resulted in a wide variation in data quality. It is therefore essential to quantify the achieved completeness and possible contamination of a reconstructed genome before it is used in subsequent analyses. The classical approach for the estimation of quality indices solely relies on a relatively small number of universal single-copy genes. Recent tools try to extend the genomic coverage of estimates for an increased accuracy.

RESULTS

We developed CoCoPyE, a fast tool based on a novel 2-stage feature extraction and transformation scheme. First, it identifies genomic markers and then refines the marker-based estimates with a machine learning approach. In our simulation studies, CoCoPyE showed a more accurate prediction of quality indices than the existing tools. While the CoCoPyE web server offers an easy way to try out the tool, the freely available Python implementation enables integration into existing genome reconstruction pipelines.

CONCLUSIONS

CoCoPyE provides a new approach to assess the quality of genome data. It complements and improves existing tools and may help researchers to better distinguish between low-quality draft and high-quality genome assemblies in metagenome sequencing projects.

A unified compendium of prokaryotic and viral genomes from over 300 anaerobic digestion microbiomes

BACKGROUND

The anaerobic digestion process degrades organic matter into simpler compounds and occurs in strictly anaerobic and microaerophilic environments. The process is carried out by a diverse community of microorganisms where each species has a unique role and it has relevant biotechnological applications since it is used for biogas production. Some aspects of the microbiome, including its interaction with phages, remains still unclear: a better comprehension of the community composition and role of each species is crucial for a cured understanding of the carbon cycle in anaerobic systems and improving biogas production.

RESULTS

The primary objective of this study was to expand our understanding on the anaerobic digestion microbiome by jointly analyzing its prokaryotic and viral components. By integrating 192 additional datasets into a previous metagenomic database, the binning process generated 11,831 metagenome-assembled genomes from 314 metagenome samples published between 2014 and 2022, belonging to 4,568 non-redundant species based on ANI calculation and quality verification. CRISPR analysis on these genomes identified 76 archaeal genomes with active phage interactions. Moreover, single-nucleotide variants further pointed to archaea as the most critical members of the community. Among the MAGs, two methanogenic archaea, Methanothrix sp. 43zhSC_152 and Methanoculleus sp. 52maCN_3230, had the highest number of SNVs, with the latter having almost double the density of most other MAGs.

CONCLUSIONS

This study offers a more comprehensive understanding of microbial community structures that thrive at different temperatures. The findings revealed that the fraction of archaeal species characterized at the genome level and reported in public databases is higher than that of bacteria, although still quite limited. The identification of shared spacers between phages and microbes implies a history of phage-bacterial interactions, and specifically lysogenic infections. A significant number of SNVs were identified, primarily comprising synonymous and nonsynonymous variants. Together, the findings indicate that methanogenic archaea are subject to intense selective pressure and suggest that genomic variants play a critical role in the anaerobic digestion process. Overall, this study provides a more balanced and diverse representation of the anaerobic digestion microbiota in terms of geographic location, temperature range and feedstock utilization.

How to Obtain and Compare Metagenome-Assembled Genomes

Metagenome-assembled genomes, or MAGs, are genomes retrieved from metagenome datasets. In the vast majority of cases, MAGs are genomes from prokaryotic species that have not been isolated or cultivated in the lab. They, therefore, provide us with information on these species that are impossible to obtain otherwise, at least until new cultivation methods are devised. Thanks to improvements and cost reductions of DNA sequencing technologies and growing interest in microbial ecology, the rise in number of MAGs in genome repositories has been exponential. This chapter covers the basics of MAG retrieval and processing and provides a practical step-by-step guide using a real dataset and state-of-the-art tools for MAG analysis and comparison.

scMicrobe PTA: near complete genomes from single bacterial cells

Microbial genomes produced by standard single-cell amplification methods are largely incomplete. Here, we show that primary template-directed amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard multiple displacement amplification-based approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

Phylogenetics and environmental distribution of nitric oxide-forming nitrite reductases reveal their distinct functional and ecological roles

The two evolutionarily unrelated nitric oxide-producing nitrite reductases, NirK and NirS, are best known for their redundant role in denitrification. They are also often found in organisms that do not perform denitrification. To assess the functional roles of the two enzymes and to address the sequence and structural variation within each, we reconstructed robust phylogenies of both proteins with sequences recovered from 6973 isolate and metagenome-assembled genomes and identified 32 well-supported clades of structurally distinct protein lineages. We then inferred the potential niche of each clade by considering other functional genes of the organisms carrying them as well as the relative abundances of each nir gene in 4082 environmental metagenomes across diverse aquatic, terrestrial, host-associated, and engineered biomes. We demonstrate that Nir phylogenies recapitulate ecology distinctly from the corresponding organismal phylogeny. While some clades of the nitrite reductase were equally prevalent across biomes, others had more restricted ranges. Nitrifiers make up a sizeable proportion of the nitrite-reducing community, especially for NirK in marine waters and dry soils. Furthermore, the two reductases showed distinct associations with genes involved in oxidizing and reducing other compounds, indicating that the NirS and NirK activities may be linked to different elemental cycles. Accordingly, the relative abundance and diversity of NirS versus NirK vary between biomes. Our results show the divergent ecological roles NirK and NirS-encoding organisms may play in the environment and provide a phylogenetic framework to distinguish the traits associated with organisms encoding the different lineages of nitrite reductases.

Functional prediction of proteins from the human gut archaeome

The human gastrointestinal tract contains diverse microbial communities, including archaea. Among them, Methanobrevibacter smithii represents a highly active and clinically relevant methanogenic archaeon, being involved in gastrointestinal disorders, such as inflammatory bowel disease and obesity. Herein, we present an integrated approach using sequence and structure information to improve the annotation of M. smithii proteins using advanced protein structure prediction and annotation tools, such as AlphaFold2, trRosetta, ProFunc, and DeepFri. Of an initial set of 873 481 archaeal proteins, we found 707 754 proteins exclusively present in the human gut. Having analysed archaeal proteins together with 87 282 994 bacterial proteins, we identified unique archaeal proteins and archaeal-bacterial homologs. We then predicted and characterized functional domains and structures of 73 unique and homologous archaeal protein clusters linked the human gut and M. smithii. We refined annotations based on the predicted structures, extending existing sequence similarity-based annotations. We identified gut-specific archaeal proteins that may be involved in defense mechanisms, virulence, adhesion, and the degradation of toxic substances. Interestingly, we identified potential glycosyltransferases that could be associated with N-linked and O-glycosylation. Additionally, we found preliminary evidence for interdomain horizontal gene transfer between Clostridia species and M. smithii, which includes sporulation Stage V proteins AE and AD. Our study broadens the understanding of archaeal biology, particularly M. smithii, and highlights the importance of considering both sequence and structure for the prediction of protein function.

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

With the explosion of available genomic information, comparative genomics has become a central approach to understanding microbial ecology and evolution. We developed DiGAlign (https://www.genome.jp/digalign/), a web server that provides versatile functionality for comparative genomics with an intuitive interface. It allows the user to perform the highly customizable visualization of a synteny map by simply uploading nucleotide sequences of interest, ranging from a specific region to the whole genome landscape of microorganisms and viruses. DiGAlign will serve a wide range of biological researchers, particularly experimental biologists, with multifaceted features that allow the rapid characterization of genomic sequences of interest and the generation of a publication-ready figure.

Metabolic and phylogenetic diversity in the phylum Nitrospinota revealed by comparative genome analyses

The most abundant known nitrite-oxidizing bacteria in the marine water column belong to the phylum Nitrospinota. Despite their importance in marine nitrogen cycling and primary production, there are only few cultured representatives that all belong to the class Nitrospinia. Moreover, although Nitrospinota were traditionally thought to be restricted to marine environments, metagenome-assembled genomes have also been recovered from groundwater. Over the recent years, metagenomic sequencing has led to the discovery of several novel classes of Nitrospinota (UBA9942, UBA7883, 2-12-FULL-45-22, JACRGO01, JADGAW01), which remain uncultivated and have not been analyzed in detail. Here, we analyzed a nonredundant set of 98 Nitrospinota genomes with focus on these understudied Nitrospinota classes and compared their metabolic profiles to get insights into their potential role in biogeochemical element cycling. Based on phylogenomic analysis and average amino acid identities, the highly diverse phylum Nitrospinota could be divided into at least 33 different genera, partly with quite distinct metabolic capacities. Our analysis shows that not all Nitrospinota are nitrite oxidizers and that members of this phylum have the genomic potential to use sulfide and hydrogen for energy conservation. This study expands our knowledge of the phylogeny and potential ecophysiology of the phylum Nitrospinota and offers new avenues for the isolation and cultivation of these elusive bacteria.

Putative novel hydrogen- and iron-oxidizing sheath-producing Zetaproteobacteria thrive at the Fåvne deep-sea hydrothermal vent field

IMPORTANCE

Knowledge on microbial iron oxidation is important for understanding the cycling of iron, carbon, nitrogen, nutrients, and metals. The current study yields important insights into the niche sharing, diversification, and Fe(III) oxyhydroxide morphology of Ghiorsea, an iron- and hydrogen-oxidizing Zetaproteobacteria representative belonging to Zetaproteobacteria operational taxonomic unit 9. The study proposes that Ghiorsea exhibits a more extensive morphology of Fe(III) oxyhydroxide than previously observed. Overall, the results increase our knowledge on potential drivers of Zetaproteobacteria diversity in iron microbial mats and can eventually be used to develop strategies for the cultivation of sheath-forming Zetaproteobacteria.

Diverse and abundant phages exploit conjugative plasmids

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally-encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 64 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

Integration of molecular and computational approaches paints a holistic portrait of obscure metabolisms

Microorganisms are essential drivers of earth's geochemical cycles. However, the significance of elemental redox cycling mediated by microorganisms is often underestimated beyond the most well-studied nutrient cycles. Phosphite, (per)chlorate, and iodate are each considered esoteric substrates metabolized by microorganisms. However, recent investigations have indicated that these metabolisms are widespread and ubiquitous, affirming a need to continue studying the underlying microbiology to understand their biogeochemical effects and their interface with each other and our biosphere. This review focuses on combining canonical techniques of culturing microorganisms with modern omic approaches to further our understanding of obscure metabolic pathways and elucidate their importance in global biogeochemical cycles. Using these approaches, marker genes of interest have already been identified for phosphite, (per)chlorate, and iodate using traditional microbial physiology and genetics. Subsequently, their presence was queried to reveal the distribution of metabolic pathways in the environment using publicly available databases. In conjunction with each other, computational and experimental techniques provide a more comprehensive understanding of the location of these microorganisms, their underlying biochemistry and genetics, and how they tie into our planet's geochemical cycles.

The undiscovered biosynthetic potential of the Greenland Ice Sheet microbiome

The Greenland Ice Sheet is a biome which is mainly microbially driven. Several different niches can be found within the glacial biome for those microbes able to withstand the harsh conditions, e.g., low temperatures, low nutrient conditions, high UV radiation in summer, and contrasting long and dark winters. Eukaryotic algae can form blooms during the summer on the ice surface, interacting with communities of bacteria, fungi, and viruses. Cryoconite holes and snow are also habitats with their own microbial community. Nevertheless, the microbiome of supraglacial habitats remains poorly studied, leading to a lack of representative genomes from these environments. Under-investigated extremophiles, like those living on the Greenland Ice Sheet, may provide an untapped reservoir of chemical diversity that is yet to be discovered. In this study, an inventory of the biosynthetic potential of these organisms is made, through cataloging the presence of biosynthetic gene clusters in their genomes. There were 133 high-quality metagenome-assembled genomes (MAGs) and 28 whole genomes of bacteria obtained from samples of the ice sheet surface, cryoconite, biofilm, and snow using culturing-dependent and -independent approaches. AntiSMASH and BiG-SCAPE were used to mine these genomes and subsequently analyze the resulting predicted gene clusters. Extensive sets of predicted Biosynthetic Gene Clusters (BGCs) were collected from the genome collection, with limited overlap between isolates and MAGs. Additionally, little overlap was found in the biosynthetic potential among different environments, suggesting specialization of organisms in specific habitats. The median number of BGCs per genome was significantly higher for the isolates compared to the MAGs. The most talented producers were found among Proteobacteria. We found evidence for the capacity of these microbes to produce antimicrobials, carotenoid pigments, siderophores, and osmoprotectants, indicating potential survival mechanisms to cope with extreme conditions. The majority of identified BGCs, including those in the most prevalent gene cluster families, have unknown functions, presenting a substantial potential for bioprospecting. This study underscores the diverse biosynthetic potential in Greenland Ice Sheet genomes, revealing insights into survival strategies and highlighting the need for further exploration and characterization of these untapped resources.

An evolution-based framework for describing human gut bacteria

The human gut microbiome contains many bacterial strains of the same species ('strain-level variants'). Describing strains in a biologically meaningful way rather than purely taxonomically is an important goal but challenging due to the genetic complexity of strain-level variation. Here, we measured patterns of co-evolution across >7,000 strains spanning the bacterial tree-of-life. Using these patterns as a prior for studying hundreds of gut commensal strains that we isolated, sequenced, and metabolically profiled revealed widespread structure beneath the phylogenetic level of species. Defining strains by their co-evolutionary signatures enabled predicting their metabolic phenotypes and engineering consortia from strain genome content alone. Our findings demonstrate a biologically relevant organization to strain-level variation and motivate a new schema for describing bacterial strains based on their evolutionary history.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Potential for homoacetogenesis via the Wood-Ljungdahl pathway in Korarchaeia lineages from marine hydrothermal vents

The Wood-Ljungdahl pathway (WLP) is a key metabolic component of acetogenic bacteria where it acts as an electron sink. In Archaea, despite traditionally being linked to methanogenesis, the pathway has been found in several Thermoproteota and Asgardarchaeota lineages. In Bathyarchaeia and Lokiarchaeia, its presence has been linked to a homoacetogenic metabolism. Genomic evidence from marine hydrothermal genomes suggests that lineages of Korarchaeia could also encode the WLP. In this study, we reconstructed 50 Korarchaeia genomes from marine hydrothermal vents along the Arctic Mid-Ocean Ridge, substantially expanding the Korarchaeia class with several taxonomically novel genomes. We identified a complete WLP in several deep-branching lineages, showing that the presence of the WLP is conserved at the root of the Korarchaeia. No methyl-CoM reductases were encoded by genomes with the WLP, indicating that the WLP is not linked to methanogenesis. By assessing the distribution of hydrogenases and membrane complexes for energy conservation, we show that the WLP is likely used as an electron sink in a fermentative homoacetogenic metabolism. Our study confirms previous hypotheses that the WLP has evolved independently from the methanogenic metabolism in Archaea, perhaps due to its propensity to be combined with heterotrophic fermentative metabolisms.

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

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