Towards the biogeography of prokaryotic genes

Genomic representativeness and chimerism in large collections of SAGs and MAGs of marine prokaryoplankton

BACKGROUND

Single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs) are the predominant sources of information about the coding potential of uncultured microbial lineages, but their strengths and limitations remain poorly understood. Here, we performed a direct comparison of two previously published collections of thousands of SAGs and MAGs obtained from the same, global environment.

RESULTS

We found that SAGs were less prone to chimerism and more accurately reflected the relative abundance and the pangenome content of microbial lineages inhabiting the epipelagic of the tropical and subtropical ocean, as compared to MAGs. SAGs were also better suited to link genome information with taxa discovered through 16S rRNA amplicon analyses. Meanwhile, MAGs had the advantage of more readily recovering genomes of rare lineages.

CONCLUSIONS

Our analyses revealed the relative strengths and weaknesses of the two most commonly used genome recovery approaches in environmental microbiology. These considerations, as well as the need for better tools for genome quality assessment, should be taken into account when designing studies and interpreting data that involve SAGs or MAGs. Video Abstract.

Discovery of antimicrobial peptides in the global microbiome with machine learning

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine-learning-based approach to predict antimicrobial peptides (AMPs) within the global microbiome and leverage a vast dataset of 63,410 metagenomes and 87,920 prokaryotic genomes from environmental and host-associated habitats to create the AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, few of which match existing databases. AMPSphere provides insights into the evolutionary origins of peptides, including by duplication or gene truncation of longer sequences, and we observed that AMP production varies by habitat. To validate our predictions, we synthesized and tested 100 AMPs against clinically relevant drug-resistant pathogens and human gut commensals both in vitro and in vivo. A total of 79 peptides were active, with 63 targeting pathogens. These active AMPs exhibited antibacterial activity by disrupting bacterial membranes. In conclusion, our approach identified nearly one million prokaryotic AMP sequences, an open-access resource for antibiotic discovery.

Quantifying functional redundancy in polysaccharide-degrading prokaryotic communities

BACKGROUND

Functional redundancy (FR) is widely present, but there is no consensus on its formation process and influencing factors. Taxonomically distinct microorganisms possessing genes for the same function in a community lead to within-community FR, and distinct assemblies of microorganisms in different communities playing the same functional roles are termed between-community FR. We proposed two formulas to respectively quantify the degree of functional redundancy within and between communities and analyzed the FR degrees of carbohydrate degradation functions in global environment samples using the genetic information of glycoside hydrolases (GHs) encoded by prokaryotes.

RESULTS

Our results revealed that GHs are each encoded by multiple taxonomically distinct prokaryotes within a community, and the enzyme-encoding prokaryotes are further distinct between almost any community pairs. The within- and between-FR degrees are primarily affected by the alpha and beta community diversities, respectively, and are also affected by environmental factors (e.g., pH, temperature, and salinity). The FR degree of the prokaryotic community is determined by deterministic factors.

CONCLUSIONS

We conclude that the functional redundancy of GHs is a stabilized community characteristic. This study helps to determine the FR formation process and influencing factors and provides new insights into the relationships between prokaryotic community biodiversity and ecosystem functions. Video Abstract.

A systematic comparison of natural product potential, with an emphasis on RiPPs, by mining of bacteria of three large ecosystems

The implementation of several global microbiome studies has yielded extensive insights into the biosynthetic potential of natural microbial communities. However, studies on the distribution of several classes of ribosomally synthesized and post-translationally modified peptides (RiPPs), non-ribosomal peptides (NRPs) and polyketides (PKs) in different large microbial ecosystems have been very limited. Here, we collected a large set of metagenome-assembled bacterial genomes from marine, freshwater and terrestrial ecosystems to investigate the biosynthetic potential of these bacteria. We demonstrate the utility of public dataset collections for revealing the different secondary metabolite biosynthetic potentials among these different living environments. We show that there is a higher occurrence of RiPPs in terrestrial systems, while in marine systems, we found relatively more terpene-, NRP-, and PK encoding gene clusters. Among the many new biosynthetic gene clusters (BGCs) identified, we analyzed various Nif-11-like and nitrile hydratase leader peptide (NHLP) containing gene clusters that would merit further study, including promising products, such as mersacidin-, LAP- and proteusin analogs. This research highlights the significance of public datasets in elucidating the biosynthetic potential of microbes in different living environments and underscores the wide bioengineering opportunities within the RiPP family.

Longitudinal multi-omics analysis of host microbiome architecture and immune responses during short-term spaceflight

Maintenance of astronaut health during spaceflight will require monitoring and potentially modulating their microbiomes. However, documenting microbial shifts during spaceflight has been difficult due to mission constraints that lead to limited sampling and profiling. Here we executed a six-month longitudinal study to quantify the high-resolution human microbiome response to three days in orbit for four individuals. Using paired metagenomics and metatranscriptomics alongside single-nuclei immune cell profiling, we characterized time-dependent, multikingdom microbiome changes across 750 samples and 10 body sites before, during and after spaceflight at eight timepoints. We found that most alterations were transient across body sites; for example, viruses increased in skin sites mostly during flight. However, longer-term shifts were observed in the oral microbiome, including increased plaque-associated bacteria (for example, Fusobacteriota), which correlated with immune cell gene expression. Further, microbial genes associated with phage activity, toxin-antitoxin systems and stress response were enriched across multiple body sites. In total, this study reveals in-depth characterization of microbiome and immune response shifts experienced by astronauts during short-term spaceflight and the associated changes to the living environment, which can help guide future missions, spacecraft design and space habitat planning.

Challenges in computational discovery of bioactive peptides in 'omics data

Peptides have a plethora of activities in biological systems that can potentially be exploited biotechnologically. Several peptides are used clinically, as well as in industry and agriculture. The increase in available 'omics data has recently provided a large opportunity for mining novel enzymes, biosynthetic gene clusters, and molecules. While these data primarily consist of DNA sequences, other types of data provide important complementary information. Due to their size, the approaches proven successful at discovering novel proteins of canonical size cannot be naïvely applied to the discovery of peptides. Peptides can be encoded directly in the genome as short open reading frames (smORFs), or they can be derived from larger proteins by proteolysis. Both of these peptide classes pose challenges as simple methods for their prediction result in large numbers of false positives. Similarly, functional annotation of larger proteins, traditionally based on sequence similarity to infer orthology and then transferring functions between characterized proteins and uncharacterized ones, cannot be applied for short sequences. The use of these techniques is much more limited and alternative approaches based on machine learning are used instead. Here, we review the limitations of traditional methods as well as the alternative methods that have recently been developed for discovering novel bioactive peptides with a focus on prokaryotic genomes and metagenomes.

Functional Genomics of a Collection of Gammaproteobacteria Isolated from Antarctica

Antarctica, one of the most extreme environments on Earth, hosts diverse microbial communities. These microbes have evolved and adapted to survive in these hostile conditions, but knowledge on the molecular mechanisms underlying this process remains limited. The Italian Collection of Antarctic Bacteria (Collezione Italiana Batteri Antartici (CIBAN)), managed by the University of Messina, represents a valuable repository of cold-adapted bacterial strains isolated from various Antarctic environments. In this study, we sequenced and analyzed the genomes of 58 marine Gammaproteobacteria strains from the CIBAN collection, which were isolated during Italian expeditions from 1990 to 2005. By employing genome-scale metrics, we taxonomically characterized these strains and assigned them to four distinct genera: Pseudomonas, Pseudoalteromonas, Shewanella, and Psychrobacter. Genome annotation revealed a previously untapped functional potential, including secondary metabolite biosynthetic gene clusters and antibiotic resistance genes. Phylogenomic analyses provided evolutionary insights, while assessment of cold-shock protein presence shed light on adaptation mechanisms. Our study emphasizes the significance of CIBAN as a resource for understanding Antarctic microbial life and its biotechnological potential. The genomic data unveil new horizons for insight into bacterial existence in Antarctica.

Pre-Bleaching Coral Microbiome Is Enriched in Beneficial Taxa and Functions

Coral reef health is tightly connected to the coral holobiont, which is the association between the coral animal and a diverse microbiome functioning as a unit. The coral holobiont depends on key services such as nitrogen and sulfur cycling mediated by the associated bacteria. However, these microbial services may be impaired in response to environmental changes, such as thermal stress. A perturbed microbiome may lead to coral bleaching and disease outbreaks, which have caused an unprecedented loss in coral cover worldwide, particularly correlated to a warming ocean. The response mechanisms of the coral holobiont under high temperatures are not completely understood, but the associated microbial community is a potential source of acquired heat-tolerance. Here we investigate the effects of increased temperature on the taxonomic and functional profiles of coral surface mucous layer (SML) microbiomes in relationship to coral-algal physiology. We used shotgun metagenomics in an experimental setting to understand the dynamics of microbial taxa and genes in the SML microbiome of the coral Pseudodiploria strigosa under heat treatment. The metagenomes of corals exposed to heat showed high similarity at the level of bacterial genera and functional genes related to nitrogen and sulfur metabolism and stress response. The coral SML microbiome responded to heat with an increase in the relative abundance of taxa with probiotic potential, and functional genes for nitrogen and sulfur acquisition. Coral-algal physiology significantly explained the variation in the microbiome at taxonomic and functional levels. These consistent and specific microbial taxa and gene functions that significantly increased in proportional abundance in corals exposed to heat are potentially beneficial to coral health and thermal resistance.

Competitive fungal commensalism mitigates candidiasis pathology

The mycobiota are a critical part of the gut microbiome, but host-fungal interactions and specific functional contributions of commensal fungi to host fitness remain incompletely understood. Here, we report the identification of a new fungal commensal, Kazachstania heterogenica var. weizmannii, isolated from murine intestines. K. weizmannii exposure prevented Candida albicans colonization and significantly reduced the commensal C. albicans burden in colonized animals. Following immunosuppression of C. albicans colonized mice, competitive fungal commensalism thereby mitigated fatal candidiasis. Metagenome analysis revealed K. heterogenica or K. weizmannii presence among human commensals. Our results reveal competitive fungal commensalism within the intestinal microbiota, independent of bacteria and immune responses, that could bear potential therapeutic value for the management of C. albicans-mediated diseases.

A global survey of prokaryotic genomes reveals the eco-evolutionary pressures driving horizontal gene transfer

Horizontal gene transfer, the exchange of genetic material through means other than reproduction, is a fundamental force in prokaryotic genome evolution. Genomic persistence of horizontally transferred genes has been shown to be influenced by both ecological and evolutionary factors. However, there is limited availability of ecological information about species other than the habitats from which they were isolated, which has prevented a deeper exploration of ecological contributions to horizontal gene transfer. Here we focus on transfers detected through comparison of individual gene trees to the species tree, assessing the distribution of gene-exchanging prokaryotes across over a million environmental sequencing samples. By analysing detected horizontal gene transfer events, we show distinct functional profiles for recent versus old events. Although most genes transferred are part of the accessory genome, genes transferred earlier in evolution tend to be more ubiquitous within present-day species. We find that co-occurring, interacting and high-abundance species tend to exchange more genes. Finally, we show that host-associated specialist species are most likely to exchange genes with other host-associated specialist species, whereas species found across different habitats have similar gene exchange rates irrespective of their preferred habitat. Our study covers an unprecedented scale of integrated horizontal gene transfer and environmental information, highlighting broad eco-evolutionary trends.

Interplay of xenobiotic-degrading and antibiotic-resistant microorganisms among the microbiome found in the air, handrail, and floor of the subway station

Investigating the quality of the subway environment, especially regarding antibiotic resistance genes (ARGs) and xenobiotics, conveys ecological and health impacts. In this study, compositions and relations of microorganisms harboring ARGs and xenobiotic degradation and metabolism genes (XDGs) in the Sukhumvit subway station (MRT-SKV) in Bangkok was assessed by analyzing the taxonomic and genetic diversity of the microbiome in the air and on the surfaces of floor and handrail. The major bacteria in the MRT-SKV (including Moraxella, which was abundant in the bioaerosol and handrail samples, and Staphylococcus, which was abundant in the bioaerosol samples) were found to contain both ARGs and XDGs. The co-abundance correlation network revealed notable relationships among bacteria harboring antibiotic resistance genes (ARGs) and xenobiotic degradation genes (XDGs). Significant associations were observed between ARGs linked to glycopeptide and fluoroquinolone resistance and genes associated with benzoate, styrene, and atrazine degradation pathways, as well as between ARGs related to cephamycin, cephalosporin, and MLS resistance and XDGs associated with the cytochrome P450-dependent drug metabolism pathway. These correlations suggested that selective pressure exerted by certain xenobiotics and antibiotics can simultaneously affect both ARGs and XDGs in the environment and should favor correlations and co-survival among ARG- and XDG-containing bacteria in the environments. The correlations may occur via shared mechanisms of resistance to both xenobiotics and antibiotics. Finally, different correlation pairs were seen in different niches (air, handrail, floor) of the subway environment or different geolocations. Thus, the relationship between ARG and XDG pairs most likely depends on the unique characteristics of the niches and on the prominent types of xenobiotics and antibiotics in the subway environment. The results indicated that interactions and connections between microbial communities can impact how they function. These microorganisms can have profound effects on accumulation of xenobiotics and ARGs in the MRT-SKV.

A compendium of ruminant gastrointestinal phage genomes revealed a higher proportion of lytic phages than in any other environments

BACKGROUND

Ruminants are important livestock animals that have a unique digestive system comprising multiple stomach compartments. Despite significant progress in the study of microbiome in the gastrointestinal tract (GIT) sites of ruminants, we still lack an understanding of the viral community of ruminants. Here, we surveyed its viral ecology using 2333 samples from 10 sites along the GIT of 8 ruminant species.

RESULTS

We present the Unified Ruminant Phage Catalogue (URPC), a comprehensive survey of phages in the GITs of ruminants including 64,922 non-redundant phage genomes. We characterized the distributions of the phage genomes in different ruminants and GIT sites and found that most phages were organism-specific. We revealed that ~ 60% of the ruminant phages were lytic, which was the highest as compared with those in all other environments and certainly will facilitate their applications in microbial interventions. To further facilitate the future applications of the phages, we also constructed a comprehensive virus-bacteria/archaea interaction network and identified dozens of phages that may have lytic effects on methanogenic archaea.

CONCLUSIONS

The URPC dataset represents a useful resource for future microbial interventions to improve ruminant production and ecological environmental qualities. Phages have great potential for controlling pathogenic bacterial/archaeal species and reducing methane emissions. Our findings provide insights into the virome ecology research of the ruminant GIT and offer a starting point for future research on phage therapy in ruminants. Video Abstract.

Inter-phylum circulation of a beta-lactamase-encoding gene: a rare but observable event

Beta-lactamase-mediated degradation of beta-lactams is the most common mechanism of beta-lactam resistance in Gram-negative bacteria. Beta-lactamase-encoding genes can be transferred between closely related bacteria, but spontaneous inter-phylum transfers (between distantly related bacteria) have never been reported. Here, we describe an extended-spectrum beta-lactamase (ESBL)-encoding gene (blaMUN-1) shared between the Pseudomonadota and Bacteroidota phyla. An Escherichia coli strain was isolated from a patient in Münster (Germany). Its genome was sequenced. The ESBL-encoding gene (named blaMUN-1) was cloned, and the corresponding enzyme was characterized. The distribution of the gene among bacteria was investigated using the RefSeq Genomes database. The frequency and relative abundance of its closest homolog in the global microbial gene catalog (GMGC) were analyzed. The E. coli strain exhibited two distinct morphotypes. Each morphotype possessed two chromosomal copies of the blaMUN-1 gene, with one morphotype having two additional copies located on a phage-plasmid p0111. Each copy was located within a 7.6-kb genomic island associated with mobility. blaMUN-1 encoded for an extended-spectrum Ambler subclass A2 beta-lactamase with 43.0% amino acid identity to TLA-1. blaMUN-1 was found in species among the Bacteroidales order and in Sutterella wadsworthensis (Pseudomonadota). Its closest homolog in GMGC was detected frequently in human fecal samples. This is, to our knowledge, the first reported instance of inter-phylum transfer of an ESBL-encoding gene, between the Bacteroidota and Pseudomonadota phyla. Although the gene was frequently detected in the human gut, inter-phylum transfer was rare, indicating that inter-phylum barriers are effective in impeding the spread of ESBL-encoding genes, but not entirely impenetrable.

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.

Infant age inversely correlates with gut carriage of resistance genes, reflecting modifications in microbial carbohydrate metabolism during early life

The infant gut microbiome is increasingly recognized as a reservoir of antibiotic resistance genes, yet the assembly of gut resistome in infants and its influencing factors remain largely unknown. We characterized resistome in 4132 metagenomes from 963 infants in six countries and 4285 resistance genes were observed. The inherent resistome pattern of healthy infants (N = 272) could be distinguished by two stages: a multicompound resistance phase (Months 0-7) and a tetracycline-mupirocin-β-lactam-dominant phase (Months 8-14). Microbial taxonomy explained 40.7% of the gut resistome of healthy infants, with Escherichia (25.5%) harboring the most resistance genes. In a further analysis with all available infants (N = 963), we found age was the strongest influencer on the resistome and was negatively correlated with the overall resistance during the first 3 years (p < 0.001). Using a random-forest approach, a set of 34 resistance genes could be used to predict age (R 2 = 68.0%). Leveraging microbial host inference analyses, we inferred the age-dependent assembly of infant resistome was a result of shifts in the gut microbiome, primarily driven by changes in taxa that disproportionately harbor resistance genes across taxa (e.g., Escherichia coli more frequently harbored resistance genes than other taxa). We performed metagenomic functional profiling and metagenomic assembled genome analyses whose results indicate that the development of gut resistome was driven by changes in microbial carbohydrate metabolism, with an increasing need for carbohydrate-active enzymes from Bacteroidota and a decreasing need for Pseudomonadota during infancy. Importantly, we observed increased acquired resistance genes over time, which was related to increased horizontal gene transfer in the developing infant gut microbiome. In summary, infant age was negatively correlated with antimicrobial resistance gene levels, reflecting a composition shift in the gut microbiome, likely driven by the changing need for microbial carbohydrate metabolism during early life.

Deep learning methods in metagenomics: a review

The ever-decreasing cost of sequencing and the growing potential applications of metagenomics have led to an unprecedented surge in data generation. One of the most prevalent applications of metagenomics is the study of microbial environments, such as the human gut. The gut microbiome plays a crucial role in human health, providing vital information for patient diagnosis and prognosis. However, analysing metagenomic data remains challenging due to several factors, including reference catalogues, sparsity and compositionality. Deep learning (DL) enables novel and promising approaches that complement state-of-the-art microbiome pipelines. DL-based methods can address almost all aspects of microbiome analysis, including novel pathogen detection, sequence classification, patient stratification and disease prediction. Beyond generating predictive models, a key aspect of these methods is also their interpretability. This article reviews DL approaches in metagenomics, including convolutional networks, autoencoders and attention-based models. These methods aggregate contextualized data and pave the way for improved patient care and a better understanding of the microbiome's key role in our health.

Clinically used broad-spectrum antibiotics compromise inflammatory monocyte-dependent antibacterial defense in the lung

Hospital-acquired pneumonia (HAP) is associated with high mortality and costs, and frequently caused by multidrug-resistant (MDR) bacteria. Although prior antimicrobial therapy is a major risk factor for HAP, the underlying mechanism remains incompletely understood. Here, we demonstrate that antibiotic therapy in hospitalized patients is associated with decreased diversity of the gut microbiome and depletion of short-chain fatty acid (SCFA) producers. Infection experiments with mice transplanted with patient fecal material reveal that these antibiotic-induced microbiota perturbations impair pulmonary defense against MDR Klebsiella pneumoniae. This is dependent on inflammatory monocytes (IMs), whose fatty acid receptor (FFAR)2/3-controlled and phagolysosome-dependent antibacterial activity is compromized in mice transplanted with antibiotic-associated patient microbiota. Collectively, we characterize how clinically relevant antibiotics affect antimicrobial defense in the context of human microbiota, and reveal a critical impairment of IM´s antimicrobial activity. Our study provides additional arguments for the rational use of antibiotics and offers mechanistic insights for the development of novel prophylactic strategies to protect high-risk patients from HAP.

Reduced PaxillinB localization to cell-substrate adhesions promotes cell migration in Dictyostelium

Many cells adhere to extracellular matrix for efficient cell migration. This adhesion is mediated by focal adhesions, a protein complex linking the extracellular matrix to the intracellular cytoskeleton. Focal adhesions have been studied extensively in mesenchymal cells, but recent research in physiological contexts and amoeboid cells suggest focal adhesion regulation differs from the mesenchymal focal adhesion paradigm. We used Dictyostelium discoideum to uncover new mechanisms of focal adhesion regulation, as Dictyostelium are amoeboid cells that form focal adhesion-like structures for migration. We show that PaxillinB, the Dictyostelium homologue of Paxillin, localizes to dynamic focal adhesion-like structures during Dictyostelium migration. Unexpectedly, reduced PaxillinB recruitment to these structures increases Dictyostelium cell migration. Quantitative analysis of focal adhesion size and dynamics show that lack of PaxillinB recruitment to focal adhesions does not alter focal adhesion size, but rather increases focal adhesion turnover. These findings are in direct contrast to Paxillin function at focal adhesions during mesenchymal migration, challenging the established focal adhesion model.

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.

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.

A bacterial sensor taxonomy across earth ecosystems for machine learning applications

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and "transduce" signals to adjust internal processes. We hypothesized that an ecosystem's unique stimuli leave a sensor "fingerprint," able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection's nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient's disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized model's feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

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.

Low-level resource partitioning supports coexistence among functionally redundant bacteria during successional dynamics

Members of microbial communities can substantially overlap in substrate use. However, what enables functionally redundant microorganisms to coassemble or even stably coexist remains poorly understood. Here, we show that during unstable successional dynamics on complex, natural organic matter, functionally redundant bacteria can coexist by partitioning low-concentration substrates even though they compete for one simple, dominant substrate. We allowed ocean microbial communities to self-assemble on leachates of the brown seaweed Fucus vesiculosus and then analyzed the competition among 10 taxonomically diverse isolates representing two distinct stages of the succession. All, but two isolates, exhibited an average of 90% ± 6% pairwise overlap in resource use, and functional redundancy of isolates from the same assembly stage was higher than that from between assembly stages, leading us to construct a simpler four-isolate community with two isolates from each of the early and late stages. We found that, although the short-term dynamics of the four-isolate communities in F. vesiculosus leachate was dependent on initial isolate ratios, in the long term, the four isolates stably coexist in F. vesiculosus leachate, albeit with some strains at low abundance. We therefore explored the potential for nonredundant substrate use by genomic content analysis and RNA expression patterns. This analysis revealed that the four isolates mainly differed in peripheral metabolic pathways, such as the ability to degrade pyrimidine, leucine, and tyrosine, as well as aromatic substrates. These results highlight the importance of fine-scale differences in metabolic strategies for supporting the frequently observed coexistence of large numbers of rare organisms in natural microbiomes.

GDPF: a data resource for the distribution of prokaryotic protein families across the global biosphere

Microorganisms encode most of the functions of life on Earth. However, conventional research has primarily focused on specific environments such as humans, soil and oceans, leaving the distribution of functional families throughout the global biosphere poorly comprehended. Here, we present the database of the global distribution of prokaryotic protein families (GDPF, http://bioinfo.qd.sdu.edu.cn/GDPF/), a data resource on the distribution of functional families across the global biosphere. GDPF provides global distribution information for 36 334 protein families, 19 734 superfamilies and 12 089 KEGG (Kyoto Encyclopedia of Genes and Genomes) orthologs from multiple source databases, covering typical environments such as soil, oceans, animals, plants and sediments. Users can browse, search and download the distribution data of each entry in 10 000 global microbial communities, as well as conduct comparative analysis of distribution disparities among multiple entries across various environments. The GDPF data resource contributes to uncovering the geographical distribution patterns, key influencing factors and macroecological principles of microbial functions at a global level, thereby promoting research in Earth ecology and human health.

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.

From microbiome composition to functional engineering, one step at a time

SUMMARYCommunities of microorganisms (microbiota) are present in all habitats on Earth and are relevant for agriculture, health, and climate. Deciphering the mechanisms that determine microbiota dynamics and functioning within the context of their respective environments or hosts (the microbiomes) is crucially important. However, the sheer taxonomic, metabolic, functional, and spatial complexity of most microbiomes poses substantial challenges to advancing our knowledge of these mechanisms. While nucleic acid sequencing technologies can chart microbiota composition with high precision, we mostly lack information about the functional roles and interactions of each strain present in a given microbiome. This limits our ability to predict microbiome function in natural habitats and, in the case of dysfunction or dysbiosis, to redirect microbiomes onto stable paths. Here, we will discuss a systematic approach (dubbed the N+1/N-1 concept) to enable step-by-step dissection of microbiome assembly and functioning, as well as intervention procedures to introduce or eliminate one particular microbial strain at a time. The N+1/N-1 concept is informed by natural invasion events and selects culturable, genetically accessible microbes with well-annotated genomes to chart their proliferation or decline within defined synthetic and/or complex natural microbiota. This approach enables harnessing classical microbiological and diversity approaches, as well as omics tools and mathematical modeling to decipher the mechanisms underlying N+1/N-1 microbiota outcomes. Application of this concept further provides stepping stones and benchmarks for microbiome structure and function analyses and more complex microbiome intervention strategies.

Strain-level characterization of health-associated bacterial consortia that colonize the human gut during infancy

Background

The human gut microbiome develops rapidly during infancy, a key window of development coinciding with maturation of the adaptive immune system. However, little is known of the microbiome growth dynamics over the first few months of life and whether there are any generalizable patterns across human populations. We performed metagenomic sequencing on stool samples (n=94) from a cohort of infants (n=15) at monthly intervals in the first six months of life, augmenting our dataset with seven published studies for a total of 4,441 metagenomes from 1,162 infants.

Results

Strain-level de novo analysis was used to identify 592 of the most abundant organisms in the infant gut microbiome. Previously unrecognized consortia were identified which exhibited highly correlated abundances across samples and were composed of diverse species spanning multiple genera. Analysis of a cohort of infants with cystic fibrosis identified one such novel consortium of diverse Enterobacterales which was positively correlated with weight gain. While all studies showed an increased community stability during the first year of life, microbial dynamics varied widely in the first few months of life, both by study and by individual.

Conclusion

By augmenting published metagenomic datasets with data from a newly established cohort we were able to identify novel groups of organisms that are correlated with measures of robust human development. We hypothesize that the presence of these groups may impact human health in aggregate in ways that individual species may not in isolation.

Life history strategies of soil bacterial communities across global terrestrial biomes

The life history strategies of soil microbes determine their metabolic potential and their response to environmental changes. Yet these strategies remain poorly understood. Here we use shotgun metagenomes from terrestrial biomes to characterize overarching covariations of the genomic traits that capture dominant life history strategies in bacterial communities. The emerging patterns show a triangle of life history strategies shaped by two trait dimensions, supporting previous theoretical and isolate-based studies. The first dimension ranges from streamlined genomes with simple metabolisms to larger genomes and expanded metabolic capacities. As metabolic capacities expand, bacterial communities increasingly differentiate along a second dimension that reflects a trade-off between increasing capacities for environmental responsiveness or for nutrient recycling. Random forest analyses show that soil pH, C:N ratio and precipitation patterns together drive the dominant life history strategy of soil bacterial communities and their biogeographic distribution. Our findings provide a trait-based framework to compare life history strategies of soil bacteria.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses

BACKGROUND

A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes.

RESULTS

Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables.

CONCLUSIONS

Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Viral activation and ecological restructuring characterize a microbiome axis of spaceflight-associated immune activation

Maintenance of astronaut health during spaceflight will require monitoring and potentially modulating their microbiomes, which play a role in some space-derived health disorders. However, documenting the response of microbiota to spaceflight has been difficult thus far due to mission constraints that lead to limited sampling. Here, we executed a six-month longitudinal study centered on a three-day flight to quantify the high-resolution microbiome response to spaceflight. Via paired metagenomics and metatranscriptomics alongside single immune profiling, we resolved a microbiome "architecture" of spaceflight characterized by time-dependent and taxonomically divergent microbiome alterations across 750 samples and ten body sites. We observed pan-phyletic viral activation and signs of persistent changes that, in the oral microbiome, yielded plaque-associated pathobionts with strong associations to immune cell gene expression. Further, we found enrichments of microbial genes associated with antibiotic production, toxin-antitoxin systems, and stress response enriched universally across the body sites. We also used strain-level tracking to measure the potential propagation of microbial species from the crew members to each other and the environment, identifying microbes that were prone to seed the capsule surface and move between the crew. Finally, we identified associations between microbiome and host immune cell shifts, proposing both a microbiome axis of immune changes during flight as well as the sources of some of those changes. In summary, these datasets and methods reveal connections between crew immunology, the microbiome, and their likely drivers and lay the groundwork for future microbiome studies of spaceflight.

Unraveling the functional dark matter through global metagenomics

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics

Owing to the ineffectiveness of traditional culture techniques for the vast majority of microbial species, culture-independent analyses utilizing next-generation sequencing and bioinformatics have become essential for gaining insight into microbial ecology and function. This mini-review focuses on two essential methods for obtaining genetic information from uncultured prokaryotes, metagenomics and single-cell genomics. We analyzed the registration status of uncultured prokaryotic genome data from major public databases and assessed the advantages and limitations of both the methods. Metagenomics generates a significant quantity of sequence data and multiple prokaryotic genomes using straightforward experimental procedures. However, in ecosystems with high microbial diversity, such as soil, most genes are presented as brief, disconnected contigs, and lack association of highly conserved genes and mobile genetic elements with individual species genomes. Although technically more challenging, single-cell genomics offers valuable insights into complex ecosystems by providing strain-resolved genomes, addressing issues in metagenomics. Recent technological advancements, such as long-read sequencing, machine learning algorithms, and in silico protein structure prediction, in combination with vast genomic data, have the potential to overcome the current technical challenges and facilitate a deeper understanding of uncultured microbial ecosystems and microbial dark matter genes and proteins. In light of this, it is imperative that continued innovation in both methods and technologies take place to create high-quality reference genome databases that will support future microbial research and industrial applications.

Computational exploration of the global microbiome for antibiotic discovery

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine learning-based approach to predict prokaryotic antimicrobial peptides (AMPs) by leveraging a vast dataset of 63,410 metagenomes and 87,920 microbial genomes. This led to the creation of AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, the majority of which were previously unknown. We observed that AMP production varies by habitat, with animal-associated samples displaying the highest proportion of AMPs compared to other habitats. Furthermore, within different human-associated microbiota, strain-level differences were evident. To validate our predictions, we synthesized and experimentally tested 50 AMPs, demonstrating their efficacy against clinically relevant drug-resistant pathogens both in vitro and in vivo. These AMPs exhibited antibacterial activity by targeting the bacterial membrane. Additionally, AMPSphere provides valuable insights into the evolutionary origins of peptides. In conclusion, our approach identified AMP sequences within prokaryotic microbiomes, opening up new avenues for the discovery of antibiotics.

A comprehensive genomic catalog from global cold seeps

Cold seeps harbor abundant and diverse microbes with tremendous potential for biological applications and that have a significant influence on biogeochemical cycles. Although recent metagenomic studies have expanded our understanding of the community and function of seep microorganisms, knowledge of the diversity and genetic repertoire of global seep microbes is lacking. Here, we collected a compilation of 165 metagenomic datasets from 16 cold seep sites across the globe to construct a comprehensive gene and genome catalog. The non-redundant gene catalog comprised 147 million genes, and 36% of them could not be assigned to a function with the currently available databases. A total of 3,164 species-level representative metagenome-assembled genomes (MAGs) were obtained, most of which (94%) belonged to novel species. Of them, 81 ANME species were identified that cover all subclades except ANME-2d, and 23 syntrophic SRB species spanned the Seep-SRB1a, Seep-SRB1g, and Seep-SRB2 clades. The non-redundant gene and MAG catalog is a valuable resource that will aid in deepening our understanding of the functions of cold seep microbiomes.

Towards monitoring of antimicrobial resistance in the environment: For what reasons, how to implement it, and what are the data needs?

Antimicrobial resistance (AMR) is a global threat to human and animal health and well-being. To understand AMR dynamics, it is important to monitor resistant bacteria and resistance genes in all relevant settings. However, while monitoring of AMR has been implemented in clinical and veterinary settings, comprehensive monitoring of AMR in the environment is almost completely lacking. Yet, the environmental dimension of AMR is critical for understanding the dissemination routes and selection of resistant microorganisms, as well as the human health risks related to environmental AMR. Here, we outline important knowledge gaps that impede implementation of environmental AMR monitoring. These include lack of knowledge of the 'normal' background levels of environmental AMR, definition of high-risk environments for transmission, and a poor understanding of the concentrations of antibiotics and other chemical agents that promote resistance selection. Furthermore, there is a lack of methods to detect resistance genes that are not already circulating among pathogens. We conclude that these knowledge gaps need to be addressed before routine monitoring for AMR in the environment can be implemented on a large scale. Yet, AMR monitoring data bridging different sectors is needed in order to fill these knowledge gaps, which means that some level of national, regional and global AMR surveillance in the environment must happen even without all scientific questions answered. With the possibilities opened up by rapidly advancing technologies, it is time to fill these knowledge gaps. Doing so will allow for specific actions against environmental AMR development and spread to pathogens and thereby safeguard the health and wellbeing of humans and animals.

SemiBin2: self-supervised contrastive learning leads to better MAGs for short- and long-read sequencing

MOTIVATION

Metagenomic binning methods to reconstruct metagenome-assembled genomes (MAGs) from environmental samples have been widely used in large-scale metagenomic studies. The recently proposed semi-supervised binning method, SemiBin, achieved state-of-the-art binning results in several environments. However, this required annotating contigs, a computationally costly and potentially biased process.

RESULTS

We propose SemiBin2, which uses self-supervised learning to learn feature embeddings from the contigs. In simulated and real datasets, we show that self-supervised learning achieves better results than the semi-supervised learning used in SemiBin1 and that SemiBin2 outperforms other state-of-the-art binners. Compared to SemiBin1, SemiBin2 can reconstruct 8.3-21.5% more high-quality bins and requires only 25% of the running time and 11% of peak memory usage in real short-read sequencing samples. To extend SemiBin2 to long-read data, we also propose ensemble-based DBSCAN clustering algorithm, resulting in 13.1-26.3% more high-quality genomes than the second best binner for long-read data.

AVAILABILITY AND IMPLEMENTATION

SemiBin2 is available as open source software at https://github.com/BigDataBiology/SemiBin/ and the analysis scripts used in the study can be found at https://github.com/BigDataBiology/SemiBin2_benchmark.

Global analysis of the biosynthetic chemical space of marine prokaryotes

BACKGROUND

Marine prokaryotes are a rich source of novel bioactive secondary metabolites for drug discovery. Recent genome mining studies have revealed their great potential to bio-synthesize novel secondary metabolites. However, the exact biosynthetic chemical space encoded by the marine prokaryotes has yet to be systematically evaluated.

RESULTS

We first investigated the secondary metabolic potential of marine prokaryotes by analyzing the diversity and novelty of the biosynthetic gene clusters (BGCs) in 7541 prokaryotic genomes from cultivated and single cells, along with 26,363 newly assembled medium-to-high-quality genomes from marine environmental samples. To quantitatively evaluate the unexplored biosynthetic chemical space of marine prokaryotes, the clustering thresholds for constructing the biosynthetic gene cluster and molecular networks were optimized to reach a similar level of the chemical similarity between the gene cluster family (GCF)-encoded metabolites and molecular family (MF) scaffolds using the MIBiG database. The global genome mining analysis demonstrated that the predicted 70,011 BGCs were organized into 24,536 mostly new (99.5%) GCFs, while the reported marine prokaryotic natural products were only classified into 778 MFs at the optimized clustering thresholds. The number of MF scaffolds is only 3.2% of the number of GCF-encoded scaffolds, suggesting that at least 96.8% of the secondary metabolic potential in marine prokaryotes is untapped. The unexplored biosynthetic chemical space of marine prokaryotes was illustrated by the 88 potential novel antimicrobial peptides encoded by ribosomally synthesized and post-translationally modified peptide BGCs. Furthermore, a sea-water-derived Aquimarina strain was selected to illustrate the diverse biosynthetic chemical space through untargeted metabolomics and genomics approaches, which identified the potential biosynthetic pathways of a group of novel polyketides and two known compounds (didemnilactone B and macrolactin A 15-ketone).

CONCLUSIONS

The present bioinformatics and cheminformatics analyses highlight the promising potential to explore the biosynthetic chemical diversity of marine prokaryotes and provide valuable knowledge for the targeted discovery and biosynthesis of novel marine prokaryotic natural products. Video Abstract.

ResFinderFG v2.0: a database of antibiotic resistance genes obtained by functional metagenomics

Metagenomics can be used to monitor the spread of antibiotic resistance genes (ARGs). ARGs found in databases such as ResFinder and CARD primarily originate from culturable and pathogenic bacteria, while ARGs from non-culturable and non-pathogenic bacteria remain understudied. Functional metagenomics is based on phenotypic gene selection and can identify ARGs from non-culturable bacteria with a potentially low identity shared with known ARGs. In 2016, the ResFinderFG v1.0 database was created to collect ARGs from functional metagenomics studies. Here, we present the second version of the database, ResFinderFG v2.0, which is available on the Center of Genomic Epidemiology web server (https://cge.food.dtu.dk/services/ResFinderFG/). It comprises 3913 ARGs identified by functional metagenomics from 50 carefully curated datasets. We assessed its potential to detect ARGs in comparison to other popular databases in gut, soil and water (marine + freshwater) Global Microbial Gene Catalogues (https://gmgc.embl.de). ResFinderFG v2.0 allowed for the detection of ARGs that were not detected using other databases. These included ARGs conferring resistance to beta-lactams, cycline, phenicol, glycopeptide/cycloserine and trimethoprim/sulfonamide. Thus, ResFinderFG v2.0 can be used to identify ARGs differing from those found in conventional databases and therefore improve the description of resistomes.

Quantifying Shared and Unique Gene Content across 17 Microbial Ecosystems

Measuring microbial diversity is traditionally based on microbe taxonomy. Here, in contrast, we aimed to quantify heterogeneity in microbial gene content across 14,183 metagenomic samples spanning 17 ecologies, including 6 human associated, 7 nonhuman host associated, and 4 in other nonhuman host environments. In total, we identified 117,629,181 nonredundant genes. The vast majority of genes (66%) occurred in only one sample (i.e., "singletons"). In contrast, we found 1,864 sequences present in every metagenome, but not necessarily every bacterial genome. Additionally, we report data sets of other ecology-associated genes (e.g., abundant in only gut ecosystems) and simultaneously demonstrated that prior microbiome gene catalogs are both incomplete and inaccurately cluster microbial genetic life (e.g., at gene sequence identities that are too restrictive). We provide our results and the sets of environmentally differentiating genes described above at http://www.microbial-genes.bio. IMPORTANCE The amount of shared genetic elements has not been quantified between the human microbiome and other host- and non-host-associated microbiomes. Here, we made a gene catalog of 17 different microbial ecosystems and compared them. We show that most species shared between environment and human gut microbiomes are pathogens and that prior gene catalogs described as "nearly complete" are far from it. Additionally, over two-thirds of all genes only appear in a single sample, and only 1,864 genes (0.001%) are found in all types of metagenomes. These results highlight the large diversity between metagenomes and reveal a new, rare class of genes, those found in every type of metagenome, but not every microbial genome.

Microbiome-based enrichment pattern mining has enabled a deeper understanding of the biome-species-function relationship

Microbes live in diverse habitats (i.e. biomes), yet their species and genes were biome-specific, forming enrichment patterns. These enrichment patterns have mirrored the biome-species-function relationship, which is shaped by ecological and evolutionary principles. However, a grand picture of these enrichment patterns, as well as the roles of external and internal factors in driving these enrichment patterns, remain largely unexamined. In this work, we have examined the enrichment patterns based on 1705 microbiome samples from four representative biomes (Engineered, Gut, Freshwater, and Soil). Moreover, an "enrichment sphere" model was constructed to elucidate the regulatory principles behind these patterns. The driving factors for this model were revealed based on two case studies: (1) The copper-resistance genes were enriched in Soil biomes, owing to the copper contamination and horizontal gene transfer. (2) The flagellum-related genes were enriched in the Freshwater biome, due to high fluidity and vertical gene accumulation. Furthermore, this enrichment sphere model has valuable applications, such as in biome identification for metagenome samples, and in guiding 3D structure modeling of proteins. In summary, the enrichment sphere model aims towards creating a bluebook of the biome-species-function relationships and be applied in many fields.

Linking prokaryotic genome size variation to metabolic potential and environment

While theories and models have appeared to explain genome size as a result of evolutionary processes, little work has shown that genome sizes carry ecological signatures. Our work delves into the ecological implications of microbial genome size variation in benthic and pelagic habitats across environmental gradients of the brackish Baltic Sea. While depth is significantly associated with genome size in benthic and pelagic brackish metagenomes, salinity is only correlated to genome size in benthic metagenomes. Overall, we confirm that prokaryotic genome sizes in Baltic sediments (3.47 Mbp) are significantly bigger than in the water column (2.96 Mbp). While benthic genomes have a higher number of functions than pelagic genomes, the smallest genomes coded for a higher number of module steps per Mbp for most of the functions irrespective of their environment. Some examples of this functions are amino acid metabolism and central carbohydrate metabolism. However, we observed that nitrogen metabolism was almost absent in pelagic genomes and was mostly present in benthic genomes. Finally, we also show that Bacteria inhabiting Baltic sediments and water column not only differ in taxonomy, but also in their metabolic potential, such as the Wood-Ljungdahl pathway or the presence of different hydrogenases. Our work shows how microbial genome size is linked to abiotic factors in the environment, metabolic potential and taxonomic identity of Bacteria and Archaea within aquatic ecosystems.

Evolutionary ecology of microbial populations inhabiting deep sea sediments associated with cold seeps

Deep sea cold seep sediments host abundant and diverse microbial populations that significantly influence biogeochemical cycles. While numerous studies have revealed their community structure and functional capabilities, little is known about genetic heterogeneity within species. Here, we examine intraspecies diversity patterns of 39 abundant species identified in sediment layers down to 430 cm below the sea floor across six cold seep sites. These populations are grouped as aerobic methane-oxidizing bacteria, anaerobic methanotrophic archaea and sulfate-reducing bacteria. Different evolutionary trajectories are observed at the genomic level among these physiologically and phylogenetically diverse populations, with generally low rates of homologous recombination and strong purifying selection. Functional genes related to methane (pmoA and mcrA) and sulfate (dsrA) metabolisms are under strong purifying selection in most species investigated. These genes differ in evolutionary trajectories across phylogenetic clades but are functionally conserved across sites. Intrapopulation diversification of genomes and their mcrA and dsrA genes is depth-dependent and subject to different selection pressure throughout the sediment column redox zones at different sites. These results highlight the interplay between ecological processes and the evolution of key bacteria and archaea in deep sea cold seep extreme environments, shedding light on microbial adaptation in the subseafloor biosphere.

The Epidermal Microbiome Within an Aggregation of Leopard Sharks (Triakis semifasciata) Has Taxonomic Flexibility with Gene Functional Stability Across Three Time-points

The epidermis of Chondrichthyan fishes consists of dermal denticles with production of minimal but protein-rich mucus that collectively, influence the attachment and biofilm development of microbes, facilitating a unique epidermal microbiome. Here, we use metagenomics to provide the taxonomic and functional characterization of the epidermal microbiome of the Triakis semifasciata (leopard shark) at three time-points collected across 4 years to identify links between microbial groups and host metabolism. Our aims include (1) describing the variation of microbiome taxa over time and identifying recurrent microbiome members (present across all time-points); (2) investigating the relationship between the recurrent and flexible taxa (those which are not found consistently across time-points); (3) describing the functional compositions of the microbiome which may suggest links with the host metabolism; and (4) identifying whether metabolic processes are shared across microbial genera or are unique to specific taxa. Microbial members of the microbiome showed high similarity between all individuals (Bray-Curtis similarity index = 82.7, where 0 = no overlap, 100 = total overlap) with the relative abundance of those members varying across sampling time-points, suggesting flexibility of taxa in the microbiome. One hundred and eighty-eight genera were identified as recurrent, including Pseudomonas, Erythrobacter, Alcanivorax, Marinobacter, and Sphingopxis being consistently abundant across time-points, while Limnobacter and Xyella exhibited switching patterns with high relative abundance in 2013, Sphingobium and Sphingomona in 2015, and Altermonas, Leeuwenhoekiella, Gramella, and Maribacter in 2017. Of the 188 genera identified as recurrent, the top 19 relatively abundant genera formed three recurrent groups. The microbiome also displayed high functional similarity between individuals (Bray-Curtis similarity index = 97.6) with gene function composition remaining consistent across all time-points. These results show that while the presence of microbial genera exhibits consistency across time-points, their abundances do fluctuate. Microbial functions however remain stable across time-points; thus, we suggest the leopard shark microbiomes exhibit functional redundancy. We show coexistence of microbes hosted in elasmobranch microbiomes that encode genes involved in utilizing nitrogen, but not fixing nitrogen, degrading urea, and resistant to heavy metal.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Stingray epidermal microbiomes are species-specific with local adaptations

Marine host-associated microbiomes are affected by a combination of species-specific (e.g., host ancestry, genotype) and habitat-specific features (e.g., environmental physiochemistry and microbial biogeography). The stingray epidermis provides a gradient of characteristics from high dermal denticles coverage with low mucus to reduce dermal denticles and high levels of mucus. Here we investigate the effects of host phylogeny and habitat by comparing the epidermal microbiomes of Myliobatis californica (bat rays) with a mucus rich epidermis, and Urobatis halleri (round rays) with a mucus reduced epidermis from two locations, Los Angeles and San Diego, California (a 150 km distance). We found that host microbiomes are species-specific and distinct from the water column, however composition of M. californica microbiomes showed more variability between individuals compared to U. halleri. The variability in the microbiome of M. californica caused the microbial taxa to be similar across locations, while U. halleri microbiomes were distinct across locations. Despite taxonomic differences, Shannon diversity is the same across the two locations in U. halleri microbiomes suggesting the taxonomic composition are locally adapted, but diversity is maintained by the host. Myliobatis californica and U. halleri microbiomes maintain functional similarity across Los Angeles and San Diego and each ray showed several unique functional genes. Myliobatis californica has a greater relative abundance of RNA Polymerase III-like genes in the microbiome than U. halleri, suggesting specific adaptations to a heavy mucus environment. Construction of Metagenome Assembled Genomes (MAGs) identified novel microbial species within Rhodobacteraceae, Moraxellaceae, Caulobacteraceae, Alcanivoracaceae and Gammaproteobacteria. All MAGs had a high abundance of active RNA processing genes, heavy metal, and antibiotic resistant genes, suggesting the stingray mucus supports high microbial growth rates, which may drive high levels of competition within the microbiomes increasing the antimicrobial properties of the microbes.

Bio-Organic Fertilizer Promotes Pear Yield by Shaping the Rhizosphere Microbiome Composition and Functions

Bio-organic fertilizers (BOF) containing both organic amendments and beneficial microorganisms have been consistently shown to improve soils fertility and yield. However, the exact mechanisms which link amendments and yields remain disputed, and the complexity of bio-organic fertilizers may work in parallel in several ways. BOF may directly improve yield by replenishing soil nutrients or introducing beneficial microbial genes or indirectly by altering the soil microbiome to enrich native beneficial microorganisms. In this work, we aim to disentangle the relative contributions of direct and indirect effects on pear yield. We treated pear trees with either chemical fertilizer or organic fertilizer with/without the plant-beneficial bacterium Bacillus velezensis SQR9. We then assessed, in detail, soil physicochemical and biological properties (metagenome sequencing) as well as pear yield. We then evaluated the relative importance of direct and indirect effects of soil amendments on pear yield. Both organic treatments increased plant yield by up to 20%, with the addition of bacteria tripling the increase driven by organic fertilizer alone. This increase could be linked to alterations in soil physicochemical properties, bacterial community function, and metabolism. Supplementation of organic fertilizer SQR9 increased rhizosphere microbiome richness and functional diversity. Fertilizer-sensitive microbes and functions responded as whole guilds. Pear yield was most positively associated with the Mitsuaria- and Actinoplanes-dominated ecological clusters and with gene clusters involved in ion transport and secondary metabolite biosynthesis. Together, these results suggested that bio-organic fertilizers mainly act indirectly on plant yield by creating soil chemical properties which promote a plant-beneficial microbiome. IMPORTANCE Bio-organic fertilization is a widely used, eco-friendly, sustainable approach to increasing plant productivity in the agriculture and fruit industries. However, it remains unclear whether the promotion of fruit productivity is related to specific changes in microbial inoculants, the resident microbiome, and/or the physicochemical properties of rhizosphere soils. We found that bio-organic fertilizers alter soil chemical properties, thus manipulating specific microbial taxa and functions within the rhizosphere microbiome of pear plants to promote yield. Our work unveils the ecological mechanisms which underlie the beneficial impacts of bio-organic fertilizers on yield promotion in fruit orchards, which may help in the design of more efficient biofertilizers to promote sustainable fruit production.