A genomic catalog of Earth's microbiomes

National antibiotic consumption is strongly related to the prevalence of antibiotic resistance across bacterial clades

The impact of societal antibiotic consumption on the prevalence of antibiotic resistance across microbial taxa in natural environments has not yet been assessed at global scales. Here, I examine the prevalence of 155 antibiotic resistance genes (ARGs) in 300,209 bacterial genomes, from non-clinical non-human-associated terrestrial environments at over 9,600 locations in 44 countries. I then compare ARG prevalences to nationwide antibiotic consumption rates, distinguishing between different ARG types. I find that depending on country and ARG type, ARG prevalences can be extremely high; for example, the probability that a given quinolone resistance gene is present in a given strain in Thailand was estimated at 42%. Further, I find strong positive correlations between nationwide antibiotic consumption rates and mean ARG prevalences for nearly all ARG types. Thus, national antibiotic consumption leaves a signal on the prevalence of ARGs across the bacterial tree, even in non-clinical environments.

Deciphering the biosynthetic landscape of biofilms in glacier-fed streams

Glacier-fed streams are permanently cold, ultra-oligotrophic, and physically unstable environments, yet microbial life thrives in benthic biofilm communities. Within biofilms, microorganisms rely on secondary metabolites for communication and competition. However, the diversity and genetic potential of secondary metabolites in glacier-fed stream biofilms remain poorly understood. In this study, we present the first large-scale exploration of biosynthetic gene clusters (BGCs) from benthic glacier-fed stream biofilms sampled by the Vanishing Glaciers project from the world's major mountain ranges. We found a remarkable diversity of BGCs, with more than 8,000 of them identified within 2,868 prokaryotic metagenome-assembled genomes, some of them potentially conferring ecological advantages, such as UV protection and quorum sensing. The BGCs were distinct from those sourced from other aquatic microbiomes, with over 40% of them being novel. The glacier-fed stream BGCs exhibited the highest similarity to BGCs from glacier microbiomes. BGC composition displayed geographic patterns and correlated with prokaryotic alpha diversity. We also found that BGC diversity was positively associated with benthic chlorophyll a and prokaryotic diversity, indicative of more biotic interactions in more extensive biofilms. Our study provides new insights into a hitherto poorly explored microbial ecosystem, which is now changing at a rapid pace as glaciers are shrinking due to climate change.

IMPORTANCE

Glacier-fed streams are characterized by low temperatures, high turbidity, and high flow. They host a unique microbiome within biofilms, which form the foundation of the food web and contribute significantly to biogeochemical cycles. Our investigation into secondary metabolites, which likely play an important role in these complex ecosystems, found a unique genetic potential distinct from other aquatic environments. We found the potential to synthesize several secondary metabolites, which may confer ecological advantages, such as UV protection and quorum sensing. This biosynthetic diversity was positively associated with the abundance and complexity of the microbial community, as well as concentrations of chlorophyll a. In the face of climate change, our study offers new insights into a vanishing ecosystem.

A genome-scale metabolic reconstruction resource of 247,092 diverse human microbes spanning multiple continents, age groups, and body sites

Genome-scale modeling of microbiome metabolism enables the simulation of diet-host-microbiome-disease interactions. However, current genome-scale reconstruction resources are limited in scope by computational challenges. We developed an optimized and highly parallelized reconstruction and analysis pipeline to build a resource of 247,092 microbial genome-scale metabolic reconstructions, deemed APOLLO. APOLLO spans 19 phyla, contains >60% of uncharacterized strains, and accounts for strains from 34 countries, all age groups, and multiple body sites. Using machine learning, we predicted with high accuracy the taxonomic assignment of strains based on the computed metabolic features. We then built 14,451 metagenomic sample-specific microbiome community models to systematically interrogate their community-level metabolic capabilities. We show that sample-specific metabolic pathways accurately stratify microbiomes by body site, age, and disease state. APOLLO is freely available, enables the systematic interrogation of the metabolic capabilities of largely still uncultured and unclassified species, and provides unprecedented opportunities for systems-level modeling of personalized host-microbiome co-metabolism.

Peeling off the layers from microbial dark matter (MDM): recent advances, future challenges, and opportunities

Microbes represent the most common organisms on Earth; however, less than 2% of microbial species in the environment can undergo cultivation for study under laboratory conditions, and the rest of the enigmatic, microbial world remains mysterious, constituting a kind of "microbial dark matter" (MDM). In the last two decades, remarkable progress has been made in culture-dependent and culture-independent techniques. More recently, studies of MDM have relied on culture-independent techniques to recover genetic material through either unicellular genomics or shotgun metagenomics to construct single-amplified genomes (SAGs) and metagenome-assembled genomes (MAGs), respectively, which provide information about evolution and metabolism. Despite the remarkable progress made in the past decades, the functional diversity of MDM still remains uncharacterized. This review comprehensively summarizes the recently developed culture-dependent and culture-independent techniques for characterizing MDM, discussing major challenges, opportunities, and potential applications. These activities contribute to expanding our knowledge of the microbial world and have implications for various fields including Biotechnology, Bioprospecting, Functional genomics, Medicine, Evolutionary and Planetary biology. Overall, this review aims to peel off the layers from MDM, shed light on recent advancements, identify future challenges, and illuminate the exciting opportunities that lie ahead in unraveling the secrets of this intriguing microbial realm.

Proteomic diversity in bacteria: Insights and implications for bacterial identification

Mass spectrometry-based proteomics has revolutionized bacterial identification and elucidated many molecular mechanisms underlying bacterial growth, community formation, and drug resistance. However, most research has been focused on a few model bacteria, overlooking bacterial diversity. In this study, we present the most extensive bacterial proteomic resource to date, covering 303 species, 119 genera, and five phyla with over 636,000 unique expressed proteins, confirming the existence of over 38,700 hypothetical proteins. Accessible via the public resource ProteomicsDB, this dataset enables quantitative exploration of proteins within and across species. Additionally, we developed MS2Bac, a bacterial identification algorithm that queries NCBI's bacterial proteome space in two iterations. MS2Bac achieved over 99% species-level and 89% strain-level accuracy, surpassing methods like MALDI-TOF and FTIR, as demonstrated with food-derived bacterial isolates. MS2Bac also effectively identified bacteria in clinical samples, highlighting the potential of MS-based proteomics as a routine diagnostic tool.

Microbial secondary metabolites: advancements to accelerate discovery towards application

Microbial secondary metabolites not only have key roles in microbial processes and relationships but are also valued in various sectors of today's economy, especially in human health and agriculture. The advent of genome sequencing has revealed a previously untapped reservoir of biosynthetic capacity for secondary metabolites indicating that there are new biochemistries, roles and applications of these molecules to be discovered. New predictive tools for biosynthetic gene clusters (BGCs) and their associated pathways have provided insights into this new diversity. Advanced molecular and synthetic biology tools and workflows including cell-based and cell-free expression facilitate the study of previously uncharacterized BGCs, accelerating the discovery of new metabolites and broadening our understanding of biosynthetic enzymology and the regulation of BGCs. These are complemented by new developments in metabolite detection and identification technologies, all of which are important for unlocking new chemistries that are encoded by BGCs. This renaissance of secondary metabolite research and development is catalysing toolbox development to power the bioeconomy.

A metagenomic perspective on the microbial prokaryotic genome census

Following 30 years of sequencing, we assessed the phylogenetic diversity (PD) of >1.5 million microbial genomes in public databases, including metagenome-assembled genomes (MAGs) of uncultivated microbes. As compared to the vast diversity uncovered by metagenomic sequences, cultivated taxa account for a modest portion of the overall diversity, 9.73% in bacteria and 6.55% in archaea, while MAGs contribute 48.54% and 57.05%, respectively. Therefore, a substantial fraction of bacterial (41.73%) and archaeal PD (36.39%) still lacks any genomic representation. This unrepresented diversity manifests primarily at lower taxonomic ranks, exemplified by 134,966 species identified in 18,087 metagenomic samples. Our study exposes diversity hotspots in freshwater, marine subsurface, sediment, soil, and other environments, whereas human samples yielded minimal novelty within the context of existing datasets. These results offer a roadmap for future genome recovery efforts, delineating uncaptured taxa in underexplored environments and underscoring the necessity for renewed isolation and sequencing.

Subspecies phylogeny in the human gut revealed by co-evolutionary constraints across the bacterial kingdom

The human gut microbiome contains many bacterial strains of the same species ("strain-level variants") that shape microbiome function. The tremendous scale and molecular resolution at which microbial communities are being interrogated motivates addressing how to describe strain-level variants. We introduce the "Spectral Tree"-an inferred tree of relatedness built from patterns of co-evolutionary constraint between greater than 7,000 diverse bacteria. Using the Spectral Tree to describe over 600 diverse gut commensal strains that we isolated, whole-genome sequenced, and metabolically profiled revealed (1) widespread phylogenetic structure among strain-level variants, (2) the origins of subspecies phylogeny as a shared history of phage infections across humans, and (3) the key role of inter-human strain variation in predicting strain-level metabolic qualities. Overall, our work demonstrates the existence and metabolic importance of structured phylogeny below the level of species for commensal gut bacteria, motivating a redefinition of individual strains according to their evolutionary context. A record of this paper's transparent peer review process is included in the supplemental information.

Microbial functional diversity and redundancy: moving forward

Microbial functional ecology is expanding as we can now measure the traits of wild microbes that affect ecosystem functioning. Here, we review techniques and advances that could be the bedrock for a unified framework to study microbial functions. These include our newfound access to environmental microbial genomes, collections of microbial traits, but also our ability to study microbes' distribution and expression. We then explore the technical, ecological, and evolutionary processes that could explain environmental patterns of microbial functional diversity and redundancy. Next, we suggest reconciling microbiology with biodiversity-ecosystem functioning studies by experimentally testing the significance of microbial functional diversity and redundancy for the efficiency, resistance, and resilience of ecosystem processes. Such advances will aid in identifying state shifts and tipping points in microbiomes, enhancing our understanding of how and where will microbes guide Earth's biomes in the context of a changing planet.

Expanding the genomic diversity of human anelloviruses

Anelloviruses are a group of small, circular, single-stranded DNA viruses that are found ubiquitously across mammalian hosts. Here, we explored a large number of publicly available human microbiome datasets and retrieved a total of 829 anellovirus genomes, substantially expanding the known diversity of these viruses. The majority of new genomes fall within the three major human anellovirus genera: Alphatorquevirus, Betatorquevirus, and Gammatorquevirus, while we also present new genomes of the under-sampled Hetorquevirus, Memtorquevirus, and Samektorquevirus genera. We performed recombination analysis and show evidence of extensive recombination across all human anelloviruses. Interestingly, more than 95% of the detected events are between members of the same genus and only 15 inter-genus recombination events were detected. The breakpoints of recombination cluster in hotspots at the ends and outside of the ORF1 gene, while a recombination coldspot was detected within the gene. Our analysis suggests that anellovirus evolution is governed by homologous recombination; however, events between distant viruses or ones producing chimaeric ORF1s likely lead to nonviable recombinants. The large number of genomes further allowed us to examine how essential genomic features vary across anelloviruses. These include functional domains in the ORF1 protein and the nucleotide motif of the replication loop region, required for the viruses' rolling-circle replication. A subset of the genomes assembled in both this and previous studies are completely lacking these essential elements, opening up the possibility that anellovirus intracellular populations contain nonstandard viral genomes. However, low-read depth of the metagenomically assembled contigs may partly explain the lack of some features. Overall, our study highlights key features of anellovirus genomics and evolution, a largely understudied group of viruses whose potential in virus-based therapeutics is recently being explored.

DFAST_QC: quality assessment and taxonomic identification tool for prokaryotic Genomes

BACKGROUND

Accurate taxonomic classification in genome databases is essential for reliable biological research and effective data sharing. Mislabeling or inaccuracies in genome annotations can lead to incorrect scientific conclusions and hinder the reproducibility of research findings. Despite advances in genome analysis techniques, challenges persist in ensuring precise and reliable taxonomic assignments. Existing tools for genome verification often involve extensive computational resources or lengthy processing times, which can limit their accessibility and scalability for large-scale projects. There is a need for more efficient, user-friendly solutions that can handle diverse datasets and provide accurate results with minimal computational demands. This work aimed to address these challenges by introducing a novel tool that enhances taxonomic accuracy, offers a user-friendly interface, and supports large-scale analyses.

RESULTS

We introduce a novel tool for the quality control and taxonomic classification tool of prokaryotic genomes, called DFAST_QC, which is available as both a command-line tool and a web service. DFAST_QC can quickly identify species based on NCBI and GTDB taxonomies by combining genome-distance calculations using MASH with ANI calculations using Skani. We evaluated DFAST_QC's performance in species identification and found it to be highly consistent with existing taxonomic standards, successfully identifying species across diverse datasets. In several cases, DFAST_QC identified potential mislabeling of species names in public databases and highlighted discrepancies in current classifications, demonstrating its capability to uncover errors and enhance taxonomic accuracy. Additionally, the tool's efficient design allows it to operate smoothly on local machines with minimal computational requirements, making it a practical choice for large-scale genome projects.

CONCLUSIONS

DFAST_QC is a reliable and efficient tool for accurate taxonomic identification and genome quality control, well-suited for large-scale genomic studies. Its compatibility with limited-resource environments, combined with its user-friendly design, ensures seamless integration into existing workflows. DFAST_QC's ability to refine species assignments in public databases highlights its value as a complementary tool for maintaining and enhancing the accuracy of taxonomic data in genomic research. The web version is available at https://dfast.ddbj.nig.ac.jp/dqc/submit/ , and the source code for local use can be found at https://github.com/nigyta/dfast_qc .

Fish gut microbiome and its application in aquaculture and biological conservation

Understanding the diversity and function of fish gut microbiomes has advanced substantially, yet many aspects remain poorly understood, particularly the interplay among microbiota, host species, and environmental factors in the context of conservation. This review explores the composition and abundance of gut bacterial communities in key aquaculture fish groups-cyprinids, ictalurids (catfish), salmonids, and cichlids (tilapia)-alongside the model organism zebrafish, across diverse geographic regions. The findings highlight environmental habitats and host species as primary determinants of gut microbiome structure, offering a global perspective on these microbial communities. Across all fish groups, the phyla Firmicutes, Fusobacteria, and Proteobacteria consistently dominated, while temperate, sub-equatorial, and sub-tropical regions exhibited the highest microbiome diversity, underscoring the contribution of taxonomic and environmental factors. The gut bacterial diversity of farm-raised fish shows a significant divergence from that of wild-caught fish, reflecting the impacts of ecological and management differences. Understanding the dynamic responses of fish gut microbiota is vital for guiding conservation efforts, safeguarding aquatic biodiversity, and advancing sustainable aquaculture practices. Future research should leverage innovative techniques and integrative approaches, both experimental and theoretical, to uncover the functional roles of microbiomes and predict their responses to environmental changes. Expanding geographic and taxonomic coverage will be critical for creating a comprehensive framework to inform global aquaculture and conservation strategies. Collectively, this perspective highlights the transformative potential of microbiome research in addressing global challenges in aquaculture and conservation biology.

The Gut Microbiota-Related Antihyperglycemic Effect of Metformin

It is critical to sustain the diversity of the microbiota to maintain host homeostasis and health. Growing evidence indicates that changes in gut microbial biodiversity may be associated with the development of several pathologies, including type 2 diabetes mellitus (T2DM). Metformin is still the first-line drug for treatment of T2DM unless there are contra-indications. The drug primarily inhibits hepatic gluconeogenesis and increases the sensitivity of target cells (hepatocytes, adipocytes and myocytes) to insulin; however, increasing evidence suggests that it may also influence the gut. As T2DM patients exhibit gut dysbiosis, the intestinal microbiome has gained interest as a key target for metabolic diseases. Interestingly, changes in the gut microbiome were also observed in T2DM patients treated with metformin compared to those who were not. Therefore, the aim of this review is to present the current state of knowledge regarding the association of the gut microbiome with the antihyperglycemic effect of metformin. Numerous studies indicate that the reduction in glucose concentration observed in T2DM patients treated with metformin is due in part to changes in the biodiversity of the gut microbiota. These changes contribute to improved intestinal barrier integrity, increased production of short-chain fatty acids (SCFAs), regulation of bile acid metabolism, and enhanced glucose absorption. Therefore, in addition to the well-recognized reduction of gluconeogenesis, metformin also appears to exert its glucose-lowering effect by influencing gut microbiome biodiversity. However, we are only beginning to understand how metformin acts on specific microorganisms in the intestine, and further research is needed to understand its role in regulating glucose metabolism, including the impact of this remarkable drug on specific microorganisms in the gut.

A functional microbiome catalogue crowdsourced from North American rivers

Predicting elemental cycles and maintaining water quality under increasing anthropogenic influence requires knowledge of the spatial drivers of river microbiomes. However, understanding of the core microbial processes governing river biogeochemistry is hindered by a lack of genome-resolved functional insights and sampling across multiple rivers. Here we used a community science effort to accelerate the sampling, sequencing and genome-resolved analyses of river microbiomes to create the Genome Resolved Open Watersheds database (GROWdb). GROWdb profiles the identity, distribution, function and expression of microbial genomes across river surface waters covering 90% of United States watersheds. Specifically, GROWdb encompasses microbial lineages from 27 phyla, including novel members from 10 families and 128 genera, and defines the core river microbiome at the genome level. GROWdb analyses coupled to extensive geospatial information reveals local and regional drivers of microbial community structuring, while also presenting foundational hypotheses about ecosystem function. Building on the previously conceived River Continuum Concept1, we layer on microbial functional trait expression, which suggests that the structure and function of river microbiomes is predictable. We make GROWdb available through various collaborative cyberinfrastructures2,3, so that it can be widely accessed across disciplines for watershed predictive modelling and microbiome-based management practices.

A Global Relationship Between Genome Size and Encoded Carbon Metabolic Strategies of Soil Bacteria

Microbial traits are critical for carbon sequestration and degradation in terrestrial ecosystems. Yet, our understanding of the relationship between carbon metabolic strategies and genomic traits like genome size remains limited. To address this knowledge gap, we conducted a global-scale meta-analysis of 2650 genomes, integrated whole-genome sequencing data, and performed a continental-scale metagenomic field study. We found that genome size was tightly associated with an increase in the ratio between genes encoding for polysaccharide decomposition and biomass synthesis that we defined as the carbon acquisition-to-biomass yield ratio (A/Y). We also show that horizontal gene transfer played a major evolutionary role in the expanded bacterial capacities in carbon acquisition. Our continental-scale field study further revealed a significantly negative relationship between the A/Y ratio and soil organic carbon stocks. Our work demonstrates a global relationship between genome size and the encoded carbon metabolic strategies of soil bacteria across terrestrial microbiomes.

Microbial diversity and biogeochemical interactions in the seismically active and CO2- rich Eger Rift ecosystem

The Eger Rift subsurface is characterized by frequent seismic activity and consistently high CO2 concentrations, making it a unique deep biosphere ecosystem and a suitable site to study the interactions between volcanism, tectonics, and microbiological activity. Pulses of geogenic H2 during earthquakes may provide substrates for methanogenic and chemolithoautotrophic processes, but very little is currently known about the role of subsurface microorganisms and their cellular processes in this type of environment. To assess the impact of geologic activity on microbial life, we analyzed the geological, geochemical, and microbiological composition of rock and sediment samples from a 238 m deep drill core, running across six lithostratigraphic zones. We evaluated the diversity and distribution of bacterial and archaeal communities. Our investigation revealed a distinct low-biomass community, with a surprisingly diverse archaeal population, providing strong support that methanogenic archaea reside in the Eger subsurface. Geochemical analysis demonstrated that ion concentrations (mostly sodium and sulfate) were highest in sediments from 50 to 100 m depth and in weathered rock below 200 m, indicating an elevated potential for ion solution in these areas. Microbial communities were dominated by common soil and water bacteria. Together with the occurrence of freshwater cyanobacteria at specific depths, these observations emphasize the heterogenous character of the sediments and are indicators for vertical groundwater movement across the Eger Rift subsurface. Our investigations also found evidence for anaerobic, autotrophic, and acidophilic communities in Eger Rift sediments, as sulfur-cycling taxa like Thiohalophilus and Desulfosporosinus were specifically enriched at depths below 100 m. The detection of methanogenic, halophilic, and ammonia-oxidizing archaeal populations demonstrate that the unique features of the Eger Rift subsurface environment provide the foundation for diverse types of microbial life, including the microbial utilization of geologically derived CO2 and, when available, H2, as a primary energy source.

A global comparison of surface and subsurface microbiomes reveals large-scale biodiversity gradients, and a marine-terrestrial divide

Subsurface environments are among Earth's largest habitats for microbial life. Yet, until recently, we lacked adequate data to accurately differentiate between globally distributed marine and terrestrial surface and subsurface microbiomes. Here, we analyzed 478 archaeal and 964 bacterial metabarcoding datasets and 147 metagenomes from diverse and widely distributed environments. Microbial diversity is similar in marine and terrestrial microbiomes at local to global scales. However, community composition greatly differs between sea and land, corroborating a phylogenetic divide that mirrors patterns in plant and animal diversity. In contrast, community composition overlaps between surface to subsurface environments supporting a diversity continuum rather than a discrete subsurface biosphere. Differences in microbial life thus seem greater between land and sea than between surface and subsurface. Diversity of terrestrial microbiomes decreases with depth, while marine subsurface diversity and phylogenetic distance to cultured isolates rivals or exceeds that of surface environments. We identify distinct microbial community compositions but similar microbial diversity for Earth's subsurface and surface environments.

Identification of diversity-generating retroelements in host-associated and environmental genomes: prevalence, diversity, and roles

BACKGROUND

The diversity-generating retroelements (DGRs) are a family of genetic elements that can produce mutations in target genes often related to ligand-binding functions, which possess a C-type lectin (CLec) domain that tolerates massive variations. They were first identified in viruses, then in bacteria and archaea from human-associated and environmental genomes. This DGR mechanism represents a fast adaptation of organisms to ever- changing environments. However, their existence, phylogenetic and structural diversity, and functions in a wide range of environments are largely unknown.

RESULTS

Here we present a study of DGR systems based on metagenome-assembled genomes (MAGs) from host-associated, aquatic, terrestrial and engineered environments. In total, we identified 861 non-redundant DGR-RTs and ~ 5.7% are new. We found that microbes associated with human hosts harbor the highest number of DGRs and also exhibit a higher prevalence of DGRs. After normalizing with genome size and including more genome data, we found that DGRs occur more frequently in organisms with smaller genomes. Overall, we identified nine main clades in the phylogenetic tree of reverse transcriptases (RTs), some comprising specific phyla and cassette architectures. We identified 38 different cassette patterns and 6 of them were shown in at least 10 DGRs, showing differences in terms of the numbers, arrangements, and orientations of their components. Finally, most of the target genes were related to ligand-binding and signaling functions, but we discovered a few cases in which the VRs were situated in domains different from the CLec.

CONCLUSIONS

Our research sheds light on the widespread prevalence of DGRs within environments and taxa, and supports the DGR phylogenetic divergence in different organisms. These variations might also occur in their structures since some cassette architectures were common in specific underrepresented phyla. In addition, we suggest that VRs could be found in domains different to the CLec, which should be further explored for organisms in scarcely studied environments.

Microbial communities associated with marine sponges from diverse geographic locations harbor biosynthetic novelty

Marine sponges are a prolific source of biologically active small molecules, many of which originate from sponge-associated bacteria. Identifying the producing bacteria is a key step in developing sustainable routes for the production of these metabolites. To facilitate the required computational analyses, we developed MetaSing, a reproducible singularity-based pipeline for assembly, identification of high-quality metagenome-assembled genomes (MAGs), and analysis of biosynthetic gene clusters (BGCs) from metagenomic short-read data. We applied this pipeline to metagenomic sequencing data from 16 marine sponges collected from New Zealand, Tonga, and the Mediterranean Sea. This analysis yielded 643 MAGs representing 510 species. Of the 2,670 BGCs identified across all samples, 70.8% were linked to a MAG. Comparison of BGCs to those identified from previously sequenced bacteria revealed high biosynthetic novelty in variety of underexplored phyla, including Poribacteria, Acidobacteriota, and Dadabacteria. Alongside the observation that each sample contains unique biosynthetic potential, this holds great promise for natural product discovery and for furthering the understanding of different sponge holobionts.IMPORTANCEDiscovery of new chemical compounds such as natural products is a crucial endeavor to combat the increasing resistance to antibiotics and other drugs. This manuscript demonstrates that microbial communities associated with marine sponges investigated in this work encode the potential to produce novel chemistry. Lesser studied bacterial taxa that are often difficult to cultivate are particularly rich in potential.

Recovery of metagenome-assembled genomes from the rumen and fecal microbiomes of Bos indicus beef cattle

Nelore is a Bos indicus beef breed that is well-adapted to tropical environments and constitutes most of the world's largest commercial cattle herd: the Brazilian bovine herd. Despite its significance, microbial genome recovery from ruminant microbiomes has largely excluded representatives from Brazilian Nelore cattle. To address this gap, this study presents a comprehensive dataset of microbial genomes recovered from the rumen and feces of 52 Brazilian Nelore bulls. A total of 1,526 non-redundant metagenome-assembled genomes (MAGs) were recovered from their gastrointestinal tract, with 497 ruminal and 486 fecal classified as high-quality. Phylogenetic analysis revealed that the bacterial MAGs fall into 12 phyla, with Firmicutes and Bacteroidota being the most predominant, while all archaeal MAGs belong to the genus Methanobrevibacter. The exploration of these microbial genomes will provide valuable insights into the metabolic potential and functional roles of individual microorganisms within host-microbiome interactions, contributing to a better understanding of the microbiome's roles in bovine performance.

Long-read sequencing sheds light on key bacteria contributing to deadwood decomposition processes

BACKGROUND

Deadwood decomposition is an essential ecological process in forest ecosystems, playing a key role in nutrient cycling and carbon sequestration by enriching soils with organic matter. This process is driven by diverse microbial communities encompassing specialized functions in breaking down organic matter, but the specific roles of individual microorganisms in this process are still not fully understood.

RESULTS

Here, we characterized the deadwood microbiome in a natural mixed temperate forest in Central Europe using PacBio HiFi long-read sequencing and a genome-resolved transcriptomics approach in order to uncover key microbial contributors to wood decomposition. We obtained high quality assemblies, which allowed attribution of complex microbial functions such as nitrogen fixation to individual microbial taxa and enabled the recovery of metagenome-assembled genomes (MAGs) from both abundant and rare deadwood bacteria. We successfully assembled 69 MAGs (including 14 high-quality and 7 single-contig genomes) from 4 samples, representing most of the abundant bacterial phyla in deadwood. The MAGs exhibited a rich diversity of carbohydrate-active enzymes (CAZymes), with Myxococcota encoding the highest number of CAZymes and the full complement of enzymes required for cellulose decomposition. For the first time we observed active nitrogen fixation by Steroidobacteraceae, as well as hemicellulose degradation and chitin recycling by Patescibacteria. Furthermore, PacBio HiFi sequencing identified over 1000 biosynthetic gene clusters, highlighting a vast potential for secondary metabolite production in deadwood, particularly in Pseudomonadota and Myxococcota.

CONCLUSIONS

PacBio HiFi long-read sequencing offers comprehensive insights into deadwood decomposition processes by advancing the identification of functional features involving multiple genes. It represents a robust tool for unraveling novel microbial genomes in complex ecosystems and allows the identification of key microorganisms contributing to deadwood decomposition.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

The way to uncovering and utilizing marine microbial resources

Recently, Chen et al. published their breakthrough results on a marine microbial genomic catalog and genetic potentials in bioprospecting in Nature, providing unprecedented opportunities for development and utilization of genetic resources of marine microorganisms. To highlight this article, we summarized and highlighted their breakthroughs seriatim.

Ecological significance of Candidatus ARS69 and Gemmatimonadota in the Arctic glacier foreland ecosystems

The Gemmatimonadota phylum has been widely detected in diverse natural environments, yet their specific ecological roles in many habitats remain poorly investigated. Similarly, the Candidatus ARS69 phylum has been identified only in a few habitats, and literature on their metabolic functions is relatively scarce. In the present study, we investigated the ecological significance of phyla Ca. ARS69 and Gemmatimonadota in the Arctic glacier foreland (GF) ecosystems through genome-resolved metagenomics. We have reconstructed the first high-quality metagenome-assembled genome (MAG) belonging to Ca. ARS69 and 12 other MAGs belonging to phylum Gemmatimonadota from the three different Arctic GF samples. We further elucidated these two groups phylogenetic lineage and their metabolic function through phylogenomic and pangenomic analysis. The analysis showed that all the reconstructed MAGs potentially belonged to novel species. The MAGs belonged to Ca. ARS69 consist about 8296 gene clusters, of which only about 8% of single-copy core genes (n = 980) were shared among them. The study also revealed the potential ecological role of Ca. ARS69 is associated with carbon fixation, denitrification, sulfite oxidation, and reduction biochemical processes in the GF ecosystems. Similarly, the study demonstrates the widespread distribution of different classes of Gemmatimonadota across wide ranges of ecosystems and their metabolic functions, including in the polar region. KEY POINTS: • Glacier foreland ecosystems act as a natural laboratory to study microbial community structure. • We have reconstructed 13 metagenome-assembled genomes from the soil samples. • All the reconstructed MAGs belonged to novel species with different metabolic processes. • Ca. ARS69 and Gemmatimonadota MAGs were found to participate in carbon fixation and denitrification processes.

Metagenomic characterization of viruses and mobile genetic elements associated with the DPANN archaeal superphylum

The archaeal superphylum DPANN (an acronym formed from the initials of the first five phyla discovered: Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanohaloarchaeota and Nanoarchaeota) is a group of ultrasmall symbionts able to survive in extreme ecosystems. The diversity and dynamics between DPANN archaea and their virome remain largely unknown. Here we use a metagenomic clustered regularly interspaced short palindromic repeats (CRISPR) screening approach to identify 97 globally distributed, non-redundant viruses and unclassified mobile genetic elements predicted to infect hosts across 8 DPANN phyla, including 7 viral groups not previously characterized. Genomic analysis suggests a diversity of viral morphologies including head-tailed, tailless icosahedral and spindle-shaped viruses with the potential to establish lytic, chronic or lysogenic infections. We also find evidence of a virally encoded Cas12f1 protein (probably originating from uncultured DPANN archaea) and a mini-CRISPR array, which could play a role in modulating host metabolism. Many metagenomes have virus-to-host ratios >10, indicating that DPANN viruses play an important role in controlling host populations. Overall, our study illuminates the underexplored diversity, functional repertoires and host interactions of the DPANN virome.

Synoptic Variation Drives Genetic Diversity and Transmission Mode of Airborne DNA Viruses in Urban Space

Airborne viruses are ubiquitous and play critical roles in maintaining ecosystem balance, however, they remain unexplored. Here, it is aimed to demonstrate that highly diverse airborne viromes carry out specific metabolic functions and use different transmission modes under different air quality conditions. A total of 263.5-Gb data are collected from 13 air samples for viral metagenomic analysis. After assembly and curation, a total of 12 484 viral contigs (1.5-184.2 kb) are assigned to 221 genus-level clades belonging to 47 families, 19 orders, and 15 classes. The composition of viral communities is influenced by weather conditions, with the main biomarker being Caudoviricetes. The most dominant viruses in these air samples belong to the dsDNA Caudoviricetes (54.0%) and ssDNA Repensiviricetes (31.2%) classes. Twelve novel candidate viruses are identified at the order/family/genus levels by alignment of complete genomes and core genes. Notably, Caudoviricetes are highly prevalent in cloudy and smoggy air, whereas Repensiviricetes are highly dominant in sunny and rainy air. Diverse auxiliary metabolic genes of airborne viruses are mainly involved in deoxynucleotide synthesis, implying their unique roles in atmosphere ecosystem. These findings deepen the understanding of the meteorological impacts on viral composition, transmission mode, and ecological roles in the air that we breathe.

Metagenome-Assembled Genomes of Pig Fecal Samples in Nine European Countries: Insights into Antibiotic Resistance Genes and Viruses

Gut microbiota plays a crucial role in the health and productivity of pigs. However, the spread of antibiotic resistance genes (ARGs) and viruses within the pig intestinal microbiota poses significant threats to animal and public health. This study utilized 181 pig samples from nine European countries and employed metagenomic assembly methods to investigate the dynamics and distribution of ARGs and viruses within the pig intestinal microbiota, aiming to observing their associations with potential bacterial hosts. We identified 4605 metagenome-assembled genomes (MAGs), corresponding to 19 bacterial phyla, 97 families, 309 genera, and a total of 449 species. Additionally, 44 MAGs were classified as archaea. Analysis of ARGs revealed 276 ARG types across 21 ARG classes, with Glycopeptide being the most abundant ARG class, followed by the class of Multidrug. Treponema D sp016293915 was identified as a primary potential bacterial host for Glycopeptide. Aligning nucleotide sequences with a viral database, we identified 1044 viruses. Among the viral genome families, Peduoviridae and Intestiviridae were the most prevalent, with CAG-914 sp000437895 being the most common potential host species for both. These findings highlight the importance of MAGs in enhancing our understanding of the gut microbiome, revealing microbial diversity, antibiotic resistance, and virus-bacteria interactions. The data analysis for the article was based on the public dataset PRJEB22062 in the European Nucleotide Archive.

Interpretation of bacterial composition patterns and community assembly processes in the rhizosphere soil of tea trees in karst areas

RESEARCH BACKGROUND AND PURPOSE

Soil microorganisms that are closely related to plants are important factors affecting plant health. Therefore, elucidating the abundant and rare bacterial species in soil associated with plant diseases is crucial for understanding ecological processes, maintaining the stability of microecological environments, and formulating microbial strategies that are consistent with modern agricultural development.

RESULTS

Tea leaf blight leads to an increase in bacterial diversity in the rhizosphere. Random processes dominate the assembly of abundant and rare taxa, while abundant taxa are also influenced by deterministic processes. In the co-occurrence network, the increase in bacterial community diversity mediated by tea cloud leaf blight enhances the stability of the network. Meanwhile, the proportion of positive correlation between rare taxa is relatively high, and the relationship between rare taxa and intermediate taxa is closer. This highly diverse bacteria community maintained the structure and stability of the community to a certain extent.

CONCLUSION

The rare taxa in the rhizosphere and the rhizosphere bacterial community mediated by tea leaf blight have high diversity, which is of great significance for maintaining the stability of the rhizosphere bacterial ecological network. In the future, we will further explore the dynamic changes and interaction patterns of the species in the rhizosphere soil affected by tea tree diseases, and their ecological functions and importance in areas of habitat fragmentation. Overall, there are many microbial resources in the rhizosphere microbiota under the influence of plant diseases that can be used for agricultural practice. The results of this study will enrich the insights into ecodynamics of bacteria in karst areas, especially in karst tea gardens.

Exploring the biosynthesis potential of permafrost microbiomes

BACKGROUND

Permafrost microbiomes are of paramount importance for the biogeochemistry of high latitude soils and while endemic biosynthetic domain sequences involved in secondary metabolism have been found in polar surface soils, the biosynthetic potential of permafrost microbiomes remains unexplored. Moreover, the nature of these ecosystems facilitates the unique opportunity to study the distribution and diversity of biosynthetic genes in relic DNA from ancient microbiomes. To explore the biosynthesis potential in permafrost, we used adenylation (AD) domain sequencing to evaluate non-ribosomal peptide (NRP) production in permafrost cores housing microbiomes separated at kilometer and kiloyear scales.

RESULTS

Permafrost microbiomes represented NRP repertoires significantly different from that of temperate soil microbiomes, but as for temperate soils, the estimated domain richness and diversity was strongly correlated to the bacterial taxonomic diversity across locations. Furthermore, we found significant differences in both community composition and AD domain composition across geographical and temporal distances. Overall, the vast majority of biosynthetic domains showed below 90% amino acid similarity to characterized BGCs, confirming the high degree of novelty of NRPs inherent to permafrost microbiomes. Using available metagenomic sequences, we further identified a high biosynthetic diversity beyond NRPs throughout arctic surface soils down to deep and ancient (megayear old) permafrost microbiomes.

CONCLUSION

We have shown that arctic permafrost microbiomes harbor a unique biosynthetic repertoire rich in hitherto undescribed NRPs. This diversity is driven by geographic separation across kilometer scales and by the bacterial taxonomic diversity between microbiomes confined in separate permafrost layers. Hence the permafrost biome represents a unique resource for studying secondary metabolism, and potentially for the discovery of novel drug leads.

kMetaShot: a fast and reliable taxonomy classifier for metagenome-assembled genomes

The advent of high-throughput sequencing (HTS) technologies unlocked the complexity of the microbial world through the development of metagenomics, which now provides an unprecedented and comprehensive overview of its taxonomic and functional contribution in a huge variety of macro- and micro-ecosystems. In particular, shotgun metagenomics allows the reconstruction of microbial genomes, through the assembly of reads into MAGs (metagenome-assembled genomes). In fact, MAGs represent an information-rich proxy for inferring the taxonomic composition and the functional contribution of microbiomes, even if the relevant analytical approaches are not trivial and still improvable. In this regard, tools like CAMITAX and GTDBtk have implemented complex approaches, relying on marker gene identification and sequence alignments, requiring a large processing time. With the aim of deploying an effective tool for fast and reliable MAG taxonomic classification, we present here kMetaShot, a taxonomy classifier based on k-mer/minimizer counting. We benchmarked kMetaShot against CAMITAX and GTDBtk by using both in silico and real mock communities and demonstrated how, while implementing a fast and concise algorithm, it outperforms the other tools in terms of classification accuracy. Additionally, kMetaShot is an easy-to-install and easy-to-use bioinformatic tool that is also suitable for researchers with few command-line skills. It is available and documented at https://github.com/gdefazio/kMetaShot.

Marine biofilms: cyanobacteria factories for the global oceans

Marine biofilms were newly revealed as a giant microbial diversity pool for global oceans. However, the cyanobacterial diversity in marine biofilms within the upper seawater column and its ecological and evolutionary implications remains undetermined. Here, we reconstructed a full picture of modern marine cyanobacteria habitats by re-analyzing 9.3 terabyte metagenomic data sets and 2,648 metagenome-assembled genomes (MAGs). The abundances of cyanobacteria lineages exclusively detected in marine biofilms were up to ninefold higher than those in seawater at similar sample size. Analyses revealed that cyanobacteria in marine biofilms are specialists with strong geographical and environmental constraints on their genome and functional adaption, which is in stark contrast to the generalistic features of seawater-derived cyanobacteria. Molecular dating suggests that the important diversifications in biofilm-forming cyanobacteria appear to coincide with the Great Oxidation Event (GOE), "boring billion" middle Proterozoic, and the Neoproterozoic Oxidation Event (NOE). These new insights suggest that marine biofilms are large and important cyanobacterial factories for the global oceans.

IMPORTANCE

Cyanobacteria, highly diverse microbial organisms, play a crucial role in Earth's oxygenation and biogeochemical cycling. However, their connection to these processes remains unclear, partly due to incomplete surveys of oceanic niches. Our study uncovered significant cyanobacterial diversity in marine biofilms, showing distinct niche differentiation compared to seawater counterparts. These patterns reflect three key stages of marine cyanobacterial diversification, coinciding with major geological events in the Earth's history.

Discovering CRISPR-Cas system with self-processing pre-crRNA capability by foundation models

The discovery of CRISPR-Cas systems has paved the way for advanced gene editing tools. However, traditional Cas discovery methods relying on sequence similarity may miss distant homologs and aren't suitable for functional recognition. With protein large language models (LLMs) evolving, there is potential for Cas system modeling without extensive training data. Here, we introduce CHOOSER (Cas HOmlog Observing and SElf-processing scReening), an AI framework for alignment-free discovery of CRISPR-Cas systems with self-processing pre-crRNA capability using protein foundation models. By using CHOOSER, we identify 11 Casλ homologs, nearly doubling the known catalog. Notably, one homolog, EphcCasλ, is experimentally validated for self-processing pre-crRNA, DNA cleavage, and trans-cleavage, showing promise for CRISPR-based pathogen detection. This study highlights an innovative approach for discovering CRISPR-Cas systems with specific functions, emphasizing their potential in gene editing.

Exploring Microbial Dark Matter for the Discovery of Novel Natural Products: Characteristics, Abundance Challenges and Methods

The objective of this review is to investigate microbial dark matter (MDM) with a focus on its potential for discovering novel natural products (NPs). This first part will examine the characteristics and abundance of these previously unexplored microbial communities, as well as the challenges faced in identifying and harnessing their unique biochemical properties and novel methods in this field. MDMs are thought to hold great potential for the discovery of novel NPs, which could have significant applications in medicine, agriculture, and industry. In recent years, there has been a growing interest in exploring MDM to unlock its potential. In fact, developments in genome-sequencing technologies and sophisticated phylogenetic procedures and metagenomic techniques have contributed to drastically make important changes in our sights on the diversity of microbial life, including the very outline of the tree of life. This has led to the development of novel technologies and methodologies for studying these elusive microorganisms, such as single-cell genomics, metagenomics, and culturomics. These approaches enable researchers to isolate and analyze individual microbial cells, as well as entire communities, providing insights into their genetic and metabolic potential. By delving into the MDM, scientists hope to uncover new compounds and biotechnological advancements that could have far-reaching impacts on various fields.

Acidic proteomes are linked to microbial alkaline preference in African lakes

Microbial amino acid composition (AA) reflects adaptive strategies of cellular and molecular regulations such as a high proportion of acidic AAs, including glutamic and aspartic acids in alkaliphiles. It remains understudied how microbial AA content is linked to their pH adaptation especially in natural environments. Here we examined prokaryotic communities and their AA composition of genes with metagenomics for 39 water and sediments of East African lakes along a gradient of pH spanning from 7.2 to 10.1. We found that Shannon diversity declined with the increasing pH and that species abundance were either positively or negatively associated with pH, indicating their distinct habitat preference in lakes. Microbial communities showed higher acidic proteomes in alkaline than neutral lakes. Species acidic proteomes were also positively correlated with their pH preference, which was consistent across major bacterial lineages. These results suggest selective pressure associated with high pH likely shape microbial amino acid composition both at the species and community levels. Comparative genome analyses further revealed that alkaliphilic microbes contained more functional genes with higher acidic AAs when compared to those in neutral conditions. These traits included genes encoding diverse classes of cation transmembrane transporters, antiporters, and compatible solute transporters, which are involved in cytoplasmic pH homeostasis and osmotic stress defense under high pH conditions. Our results provide the field evidence for the strong relationship between prokaryotic AA composition and their habitat preference and highlight amino acid optimization as strategies for environmental adaptation.

Genetic potential for aerobic respiration and denitrification in globally distributed respiratory endosymbionts

The endosymbiont Candidatus Azoamicus ciliaticola was proposed to generate ATP for its eukaryotic host, an anaerobic ciliate of the Plagiopylea class, fulfilling a function analogous to mitochondria in other eukaryotic cells. The discovery of this respiratory endosymbiosis has major implications for both evolutionary history and ecology of microbial eukaryotes. However, with only a single species described, knowledge of its environmental distribution and diversity is limited. Here we report four complete, circular metagenome assembled genomes (cMAGs) representing respiratory endosymbionts inhabiting groundwater in California, Ohio, and Germany. These cMAGs form two lineages comprising a monophyletic clade within the uncharacterized gammaproteobacterial order UBA6186, enabling evolutionary analysis of their key protein complexes. Strikingly, all four cMAGs encode a cytochrome cbb3 oxidase, which indicates that these endosymbionts have the capacity for aerobic respiration. Accordingly, we detect these respiratory endosymbionts in diverse habitats worldwide, thus further expanding the ecological scope of this respiratory symbiosis.

MAGqual: a stand-alone pipeline to assess the quality of metagenome-assembled genomes

BACKGROUND

Metagenomics, the whole genome sequencing of microbial communities, has provided insight into complex ecosystems. It has facilitated the discovery of novel microorganisms, explained community interactions and found applications in various fields. Advances in high-throughput and third-generation sequencing technologies have further fuelled its popularity. Nevertheless, managing the vast data produced and addressing variable dataset quality remain ongoing challenges. Another challenge arises from the number of assembly and binning strategies used across studies. Comparing datasets and analysis tools is complex as it requires the quantitative assessment of metagenome quality. The inherent limitations of metagenomic sequencing, which often involves sequencing complex communities, mean community members are challenging to interrogate with traditional culturing methods leading to many lacking reference sequences. MIMAG standards aim to provide a method to assess metagenome quality for comparison but have not been widely adopted.

RESULTS

To address the need for simple and quick metagenome quality assignation, here we introduce the pipeline MAGqual (Metagenome-Assembled Genome qualifier) and demonstrate its effectiveness at determining metagenomic dataset quality in the context of the MIMAG standards.

CONCLUSIONS

The MAGqual pipeline offers an accessible way to evaluate metagenome quality and generate metadata on a large scale. MAGqual is built in Snakemake to ensure readability and scalability, and its open-source nature promotes accessibility, community development, and ease of updates. MAGqual is built in Snakemake, R, and Python and is available under the MIT license on GitHub at https://github.com/ac1513/MAGqual . Video Abstract.

Genome-guided isolation of the hyperthermophilic aerobe Fervidibacter sacchari reveals conserved polysaccharide metabolism in the Armatimonadota

Few aerobic hyperthermophilic microorganisms degrade polysaccharides. Here, we describe the genome-enabled enrichment and optical tweezer-based isolation of an aerobic polysaccharide-degrading hyperthermophile, Fervidibacter sacchari, previously ascribed to candidate phylum Fervidibacteria. F. sacchari uses polysaccharides and monosaccharides for growth at 65-87.5 °C and expresses 191 carbohydrate-active enzymes (CAZymes) according to RNA-Seq and proteomics, including 31 with unusual glycoside hydrolase domains (GH109, GH177, GH179). Fluorescence in-situ hybridization and nanoscale secondary ion mass spectrometry confirmed rapid assimilation of 13C-starch in spring sediments. Purified GHs were optimally active at 80-100 °C on ten different polysaccharides. Finally, we propose reassigning Fervidibacteria as a class within phylum Armatimonadota, along with 18 other species, and show that a high number and diversity of CAZymes is a hallmark of the phylum, in both aerobic and anaerobic lineages. Our study establishes Fervidibacteria as hyperthermophilic polysaccharide degraders in terrestrial geothermal springs and suggests a broad role for Armatimonadota in polysaccharide catabolism.

Genomic diversity of phages infecting the globally widespread genus Sulfurimonas

The widespread bacterial genus Sulfurimonas is metabolically versatile and occupies a key ecological niche in different habitats, but its interaction with bacteriophages remains unexplored. Here we systematically investigated the genetic diversity, taxonomy and interaction patterns of Sulfurimonas-associated phages based on sequenced microbial genomes and metagenomes. High-confidence phage contigs related to Sulfurimonas were retrieved from various ecosystems, clustered into 61 viral operational taxonomic units across three viral realms, including Duplodnaviria, Monodnaviria and Varidnaviria. Head-tail phages of Caudoviricetes were assigned to 19 genus-level viral clusters, the majority of which were distantly related to known viruses. Notably, diverse double jelly-roll viruses and inoviruses were also linked to Sulfurimonas, representing two commonly overlooked phage groups. Historical and current phage infections were revealed, implying viral impact on the evolution of host adaptive immunity. Additionally, phages carrying auxiliary metabolic genes might benefit hosts by compensating or augmenting sulfur metabolism. This study highlights the diversity and novelty of Sulfurimonas-associated phages with divergent tailless lineages, providing basis for further investigation of phage-host interactions within this genus.

pan-Draft: automated reconstruction of species-representative metabolic models from multiple genomes

The accurate reconstruction of genome-scale metabolic models (GEMs) for unculturable species poses challenges due to the incomplete and fragmented genetic information typical of metagenome-assembled genomes (MAGs). While existing tools leverage sequence homology from single genomes, this study introduces pan-Draft, a pan-reactome-based approach exploiting recurrent genetic evidence to determine the solid core structure of species-level GEMs. By comparing MAGs clustered at the species-level, pan-Draft addresses the issues due to the incompleteness and contamination of individual genomes, providing high-quality draft models and an accessory reactions catalog supporting the gapfilling step. This approach will improve our comprehension of metabolic functions of uncultured species.

Mini review: antimicrobial compounds produced by bacteria associated with marine invertebrates

The marine environment is considered one of the most important ecosystems with high biodiversity. Microorganisms in this environment are variable and coexist with other marine organisms. The microbes associated with other marine organisms produce compounds with biological activity that may help the host's defense against invading organisms. The symbiotic association of bacteria with marine invertebrates is of ecological and biotechnological importance. Biologically active metabolites isolated from bacteria associated with marine invertebrates are considered potential sources of natural antimicrobial molecules for treating infectious diseases. Many studies have been conducted to screen the antimicrobial activity of metabolites produced by bacteria associated with marine invertebrates. This work provides an overview of the advancements in antimicrobial compound research on bacteria associated with marine invertebrates.

Annotation-free prediction of microbial dioxygen utilization

Aerobes require dioxygen (O2) to grow; anaerobes do not. However, nearly all microbes-aerobes, anaerobes, and facultative organisms alike-express enzymes whose substrates include O2, if only for detoxification. This presents a challenge when trying to assess which organisms are aerobic from genomic data alone. This challenge can be overcome by noting that O2 utilization has wide-ranging effects on microbes: aerobes typically have larger genomes encoding distinctive O2-utilizing enzymes, for example. These effects permit high-quality prediction of O2 utilization from annotated genome sequences, with several models displaying ≈80% accuracy on a ternary classification task for which blind guessing is only 33% accurate. Since genome annotation is compute-intensive and relies on many assumptions, we asked if annotation-free methods also perform well. We discovered that simple and efficient models based entirely on genomic sequence content-e.g., triplets of amino acids-perform as well as intensive annotation-based classifiers, enabling rapid processing of genomes. We further show that amino acid trimers are useful because they encode information about protein composition and phylogeny. To showcase the utility of rapid prediction, we estimated the prevalence of aerobes and anaerobes in diverse natural environments cataloged in the Earth Microbiome Project. Focusing on a well-studied O2 gradient in the Black Sea, we found quantitative correspondence between local chemistry (O2:sulfide concentration ratio) and the composition of microbial communities. We, therefore, suggest that statistical methods like ours might be used to estimate, or "sense," pivotal features of the chemical environment using DNA sequencing data.IMPORTANCEWe now have access to sequence data from a wide variety of natural environments. These data document a bewildering diversity of microbes, many known only from their genomes. Physiology-an organism's capacity to engage metabolically with its environment-may provide a more useful lens than taxonomy for understanding microbial communities. As an example of this broader principle, we developed algorithms that accurately predict microbial dioxygen utilization directly from genome sequences without annotating genes, e.g., by considering only the amino acids in protein sequences. Annotation-free algorithms enable rapid characterization of natural samples, highlighting quantitative correspondence between sequences and local O2 levels in a data set from the Black Sea. This example suggests that DNA sequencing might be repurposed as a multi-pronged chemical sensor, estimating concentrations of O2 and other key facets of complex natural settings.

Diverse and specialized metabolic capabilities of microbes in oligotrophic built environments

BACKGROUND

Built environments (BEs) are typically considered to be oligotrophic and harsh environments for microbial communities under normal, non-damp conditions. However, the metabolic functions of microbial inhabitants in BEs remain poorly understood. This study aimed to shed light on the functional capabilities of microbes in BEs by analyzing 860 representative metagenome-assembled genomes (rMAGs) reconstructed from 738 samples collected from BEs across the city of Hong Kong and from the skin surfaces of human occupants. The study specifically focused on the metabolic functions of rMAGs that are either phylogenetically novel or prevalent in BEs.

RESULTS

The diversity and composition of BE microbiomes were primarily shaped by the sample type, with Micrococcus luteus and Cutibacterium acnes being prevalent. The metabolic functions of rMAGs varied significantly based on taxonomy, even at the strain level. A novel strain affiliated with the Candidatus class Xenobia in the Candidatus phylum Eremiobacterota and two novel strains affiliated with the superphylum Patescibacteria exhibited unique functions compared with their close relatives, potentially aiding their survival in BEs and on human skins. The novel strains in the class Xenobia possessed genes for transporting nitrate and nitrite as nitrogen sources and nitrosative stress mitigation induced by nitric oxide during denitrification. The two novel Patescibacteria strains both possessed a broad array of genes for amino acid and trace element transport, while one of them carried genes for carotenoid and ubiquinone biosynthesis. The globally prevalent M. luteus in BEs displayed a large and open pangenome, with high infraspecific genomic diversity contributed by 11 conspecific strains recovered from BEs in a single geographic region. The versatile metabolic functions encoded in the large accessory genomes of M. luteus may contribute to its global ubiquity and specialization in BEs.

CONCLUSIONS

This study illustrates that the microbial inhabitants of BEs possess metabolic potentials that enable them to tolerate and counter different biotic and abiotic conditions. Additionally, these microbes can efficiently utilize various limited residual resources from occupant activities, potentially enhancing their survival and persistence within BEs. A better understanding of the metabolic functions of BE microbes will ultimately facilitate the development of strategies to create a healthy indoor microbiome. Video Abstract.

Hidden Links Between Skin Microbiome and Skin Imaging Phenome

Despite the skin microbiome has been linked to skin health and diseases, its role in modulating human skin appearance remains understudied. Using a total of 1244 face imaging phenomes and 246 cheek metagenomes, we first established three skin age indices by machine learning, including skin phenotype age (SPA), skin microbiota age (SMA), and skin integration age (SIA) as surrogates of phenotypic aging, microbial aging, and their combination, respectively. Moreover, we found that besides aging and gender as intrinsic factors, skin microbiome might also play a role in shaping skin imaging phenotypes (SIPs). Skin taxonomic and functional α diversity was positively linked to melanin, pore, pigment, and ultraviolet spot levels, but negatively linked to sebum, lightening, and porphyrin levels. Furthermore, certain species were correlated with specific SIPs, such as sebum and lightening levels negatively correlated with Corynebacterium matruchotii, Staphylococcus capitis, and Streptococcus sanguinis. Notably, we demonstrated skin microbial potential in predicting SIPs, among which the lightening level presented the least error of 1.8%. Lastly, we provided a reservoir of potential mechanisms through which skin microbiome adjusted the SIPs, including the modulation of pore, wrinkle, and sebum levels by cobalamin and heme synthesis pathways, predominantly driven by Cutibacterium acnes. This pioneering study unveils the paradigm for the hidden links between skin microbiome and skin imaging phenome, providing novel insights into how skin microbiome shapes skin appearance and its healthy aging.

Small amounts of misassembly can have disproportionate effects on pangenome-based metagenomic analyses

Individual genes from microbiomes can drive host-level phenotypes. To help identify such candidate genes, several recent tools estimate microbial gene copy numbers directly from metagenomes. These tools rely on alignments to pangenomes, which in turn are derived from the set of all individual genomes from one species. While large-scale metagenomic assembly efforts have made pangenome estimates more complete, mixed communities can also introduce contamination into assemblies, and it is unknown how robust pangenome-based metagenomic analyses are to these errors. To gain insight into this problem, we re-analyzed a case-control study of the gut microbiome in cirrhosis, focusing on commensal Clostridia previously implicated in this disease. We tested for differentially prevalent genes in the Lachnospiraceae, then investigated which were likely to be contaminants using sequence similarity searches. Out of 86 differentially prevalent genes, we found that 33 (38%) were probably contaminants originating in taxa such as Veillonella and Haemophilus, unrelated genera that were independently correlated with disease status. Our results demonstrate that even small amounts of contamination in metagenome assemblies, below typical quality thresholds, can threaten to overwhelm gene-level metagenomic analyses. However, we also show that such contaminants can be accurately identified using a method based on gene-to-species correlation. After removing these contaminants, we observe that several flagellar motility gene clusters in the Lachnospira eligens pangenome are associated with cirrhosis status. We have integrated our analyses into an analysis and visualization pipeline, PanSweep, that can automatically identify cases where pangenome contamination may bias the results of gene-resolved analyses.

Importance

Metagenome-assembled genomes, or MAGs, can be constructed without pure cultures of microbes. Large scale efforts to build MAGs have yielded more complete pangenomes (i.e., sets of all genes found in one species). Pangenomes allow us to measure strain variation in gene content, which can strongly affect phenotype. However, because MAGs come from mixed communities, they can contaminate pangenomes with unrelated DNA, and how much this impacts downstream analyses has not been studied. Using a metagenomic study of gut microbes in cirrhosis as our test case, we investigate how contamination affects analyses of microbial gene content. Surprisingly, even small, typical amounts of MAG contamination (<5%) result in disproportionately high levels of false positive associations (38%). Fortunately, we show that most contaminants can be automatically flagged, and provide a simple method for doing so. Furthermore, applying this method reveals a new association between cirrhosis and gut microbial motility.

Rapid species-level metagenome profiling and containment estimation with sylph

Profiling metagenomes against databases allows for the detection and quantification of microorganisms, even at low abundances where assembly is not possible. We introduce sylph, a species-level metagenome profiler that estimates genome-to-metagenome containment average nucleotide identity (ANI) through zero-inflated Poisson k-mer statistics, enabling ANI-based taxa detection. On the Critical Assessment of Metagenome Interpretation II (CAMI2) Marine dataset, sylph was the most accurate profiling method of seven tested. For multisample profiling, sylph took >10-fold less central processing unit time compared to Kraken2 and used 30-fold less memory. Sylph's ANI estimates provided an orthogonal signal to abundance, allowing for an ANI-based metagenome-wide association study for Parkinson disease (PD) against 289,232 genomes while confirming known butyrate-PD associations at the strain level. Sylph took <1 min and 16 GB of random-access memory to profile metagenomes against 85,205 prokaryotic and 2,917,516 viral genomes, detecting 30-fold more viral sequences in the human gut compared to RefSeq. Sylph offers precise, efficient profiling with accurate containment ANI estimation even for low-coverage genomes.

Biodiversity of mudflat intertidal viromes along the Chinese coasts

Viruses constitute the most diverse and abundant biological entities on Earth. However, our understanding of this tiniest life form in complex ecosystems remains limited. Here, we recover 20,102 viral OTUs from twelve intertidal zones along the Chinese coasts. Our analysis demonstrates high viral diversity and functional potential in intertidal zones, encoding important functional genes that can be potentially transferred to microbial hosts and mediate elemental biogeochemical cycles, especially carbon, phosphate and sulfur. Virus-host abundance dynamics vary among different microbial lineages. Viral community composition is closely associated with environmental conditions, including dissolved organic matter. Concordant biogeographic patterns are observed for viruses and microbes. Viral communities are generally habitat specific with low overlaps between intertidal and other habitats. Environmental factors and geographic distance dominate the compositional variation of intertidal viromes. Overall, these findings expand our understanding of intertidal viromes within an ecological framework, providing insights into the virus-host coevolutionary arms race.

Comprehensive evaluation of primer pairs targeting the ammonia monooxygenase subunit A gene of complete ammonia-oxidizing Nitrospira

Since the discovery of complete ammonia oxidizers (comammox) within the genus Nitrospira, their distribution and abundance across habitats have been intensively studied to better understand their ecological significance. Many primers targeting their ammonia monooxygenase subunit A gene (amoA) have been designed to detect and quantify comammox bacteria and to describe their community structure. We identified 38 published primers, but only few had high coverage and specificity for all known comammox Nitrospira or one of the two described subclades. For each target group, we comprehensively evaluated selected primer pairs using in silico analyses, endpoint PCRs, qPCRs, and amplicon sequencing on samples from various environments. Endpoint PCRs and qPCRs showed that the most commonly used primer pairs (comaA-244F/659R, comaB-244F/659R, and Ntsp-amoA162F/359R) produced several bands, which likely inflated quantifications via qPCR. In contrast, the recently published primer combinations CA377F/C576R, CB377F/C576R, and CA-CB377F/C576R resulted mostly in a single band. Furthermore, amplicon sequencing demonstrated that these primer combinations also captured the highest richness of comammox Nitrospira. Taken together, our results indicate that few existing comammox amoA primer combinations have both high specificity and coverage and that the choice of these high-specificity and high-coverage primer pairs substantially impacts the accurate detection, quantification, and community description of comammox bacteria. We, therefore, recommend using the CA377F/C576R, CB377F/C576R, and CA-CB377F/C576R primer pairs.IMPORTANCEBacteria that can fully convert ammonia via nitrite to nitrate, the complete ammonia oxidizers (comammox), were recently discovered and are found in many natural and engineered environments. PCR-based tools to study their abundance and diversity were rapidly developed, resulting in a plethora of primers available, many of which are widely used. The presence of comammox bacteria in an environment can, however, only be correctly determined if the used primers detect all members of this group while not detecting any other guilds. This study assesses the coverage and specificity of existing primers targeting comammox bacteria using both computational and standard molecular techniques, revealing large differences in their performance. The uniform usage of well-performing primers across studies could aid in generating comparable and generalizable data to better understand the importance of comammox bacteria in the environment.

A single amplified genome catalog reveals the dynamics of mobilome and resistome in the human microbiome

BACKGROUND

The increase in metagenome-assembled genomes (MAGs) has advanced our understanding of the functional characterization and taxonomic assignment within the human microbiome. However, MAGs, as population consensus genomes, often aggregate heterogeneity among species and strains, thereby obfuscating the precise relationships between microbial hosts and mobile genetic elements (MGEs). In contrast, single amplified genomes (SAGs) derived via single-cell genome sequencing can capture individual genomic content, including MGEs.

RESULTS

We introduce the first substantial SAG dataset (bbsag20) from the human oral and gut microbiome, comprising 17,202 SAGs above medium-quality without co-assembly. This collection unveils a diversity of bacterial lineages across 312 oral and 647 gut species, demonstrating different taxonomic compositions from MAGs. Moreover, the SAGs showed cellular-level evidence of the translocation of oral bacteria to the gut. We also identified broad-host-range MGEs harboring antibiotic resistance genes (ARGs), which were not detected in the MAGs.

CONCLUSIONS

The difference in taxonomic composition between SAGs and MAGs indicates that combining both methods would be effective in expanding the genome catalog. By connecting mobilomes and resistomes in individual samples, SAGs could meticulously chart a dynamic network of ARGs on MGEs, pinpointing potential ARG reservoirs and their spreading patterns in the microbial community. Video Abstract.

Plant growth promoting endophyte modulates soil ecological characteristics during the enhancement process of cadmium phytoremediation

Endophyte assisted phytoremediation of cadmium (Cd) contaminated soil represents a promising strategy. However, the precise soil ecological regulatory mechanisms by which endophyte enhance the Cd phytoextraction remain unclear. Here, we employed the plant growth promoting endophyte (PGPE) Pseudomonas sp. E3, which has been validated to effectively enhance Cd extraction in Solanum nigrum L., to investigate its regulatory mechanism on soil ecology. The results demonstrated that while PGPE inoculation resulted in minimal alterations to the physicochemical properties of the bulk soil, it led to a notable increase in acid phosphatase activity by 17.86% and urease activity by 24.85% in the rhizosphere soil. This, in turn, significantly raised the available nitrogen and phosphorus contents by 16.93% and 21.27%, respectively, in the rhizosphere soil. Additionally, PGPE inoculation effectively replenished the bioavailable fractions of Fe and Cd, which had been depleted due to root uptake. Importantly, the inoculation specifically augmented the abundance of biomarkers p_Patescibacteria, f_Saccharimonadales, and g_Saccharimonadales in the rhizosphere soil. These biomarkers exhibited a significant positive correlation with the available nutrient and metal element contents. Moreover, the co-occurrence network analysis demonstrated that the inoculation resulted in a simplified bacterial community network, which may have facilitated community synergism by displacing bacteria with a negative association. This regulation appears to occur independently of PGPE colonization. Overall, our findings suggested that PGPE also exerts a regulatory influence on soil ecological features, significantly aiding hyperaccumulators in nutrient acquisition and heavy metal accumulation.

Taxometer: Improving taxonomic classification of metagenomics contigs

For taxonomy based classification of metagenomics assembled contigs, current methods use sequence similarity to identify their most likely taxonomy. However, in the related field of metagenomic binning, contigs are routinely clustered using information from both the contig sequences and their abundance. We introduce Taxometer, a neural network based method that improves the annotations and estimates the quality of any taxonomic classifier using contig abundance profiles and tetra-nucleotide frequencies. We apply Taxometer to five short-read CAMI2 datasets and find that it increases the average share of correct species-level contig annotations of the MMSeqs2 tool from 66.6% to 86.2%. Additionally, it reduce the share of wrong species-level annotations in the CAMI2 Rhizosphere dataset by an average of two-fold for Metabuli, Centrifuge, and Kraken2. Futhermore, we use Taxometer for benchmarking taxonomic classifiers on two complex long-read metagenomics data sets where ground truth is not known. Taxometer is available as open-source software and can enhance any taxonomic annotation of metagenomic contigs.