MetaCHIP: community-level horizontal gene transfer identification through the combination of best-match and phylogenetic approaches

Co-occurrence and co-expression of antibiotic, biocide, and metal resistance genes with mobile genetic elements in microbial communities subjected to long-term antibiotic pressure: Novel insights from metagenomics and metatranscriptomics

The burgeoning of antibiotic resistance has emerged as a pressing global challenge. To gain a deeper understanding of the interactions between antibiotic resistance genes (ARGs), biocide and metal resistance genes (BRGs&MRGs), and mobile genetic elements (MGEs), this study utilized metagenomics and metatranscriptomics to investigate their co-occurrence and co-expression in two consortia subjected to long-term exposure to chloramphenicol and lincomycin. Long-term exposure to these antibiotics resulted in significant disparities in resistance profiles: ConsortiumCAP harbored 130 ARGs and 150 BRGs&MRGs, while ConsortiumLIN contained 57 ARGs and 32 BRGs&MRGs. Horizontal gene transfer (HGT) events were predicted at 125 and 300 instances in ConsortiumCAP and ConsortiumLIN, respectively, facilitating the emergence of multidrug-resistant bacteria, such as Caballeronia (10 ARGs, 2 BRGs&MRGs), Cupriavidus (2 ARGs, 10 BRGs&MRGs), and Bacillus (14 ARGs, 21 BRGs&MRGs). Chloramphenicol exposure significantly enriched genes linked to phenicol resistance (floR, capO) and co-expressed ARGs and BRGs&MRGs, while lincomycin exerted narrower effects on resistance genes. Additionally, both antibiotics modulated the expression of degradation genes and virulence factors, highlighting their role in altering bacterial substrate utilization and pathogenic traits. This study provides quantitative insights into the impact of antibiotics on microbial resistance profiles and functions at both DNA and RNA levels, highlighting the importance of reducing antibiotic pollution and limiting the spread of resistance genes in the environment.

Antibiotic resistance patterns and cross-family ARG transfer in families Burkholderiaceae and Sphingomonadaceae: A large-scale genome-wide analysis of over 10 K genomes

Members of Burkholderiaceae and Sphingomonadaceae play an active role in pollutant degradation, yet their antibiotic resistance risks are frequently overlooked. This study analyzed 9406 Burkholderiaceae and 2343 Sphingomonadaceae genomes to investigate the distribution, horizontal gene transfer (HGT), and co-occurrence patterns of antibiotic resistance genes (ARGs) and metal resistance genes (MRGs). ARGs were prevalent in Burkholderiaceae (93.2 % of genomes), dominated by bacitracin (89.0 %), multidrug (88.1 %), and beta-lactam (40.5 %) resistance, while Sphingomonadaceae exhibited lower ARG prevalence (11.6 %). Notably, Burkholderia and Caballeronia displayed high multidrug resistance (10.1 ARGs per genome) and frequent ARG-MRG co-occurrence (84.4 %). Strong ARG-MRG-MGE correlations were observed in Burkholderiaceae, suggesting MGEs play a key role in resistance dissemination. Additionally, ARGs correlated with metabolic genes, linking metabolic versatility to resistance. Genes like capO (chloramphenicol oxidase) and blaTEM-116 (beta-lactamase) were shared among distantly related genera, while mcr-5.1 (MCR phosphoethanolamine transferase) co-occurred with MRGs across Cupriavidus species, highlighting HGT and co-selection risks. ARG transfer between Burkholderiaceae, Sphingomonadaceae and clinical pathogens was frequent (114-1306 events/10,000 genome pairs), with sulfonamide resistance dominating (51.3 % of HGT). These findings highlight Burkholderiaceae and Sphingomonadaceae as critical reservoirs of resistance genes and emphasize the need for enhanced surveillance and mitigation strategies to curb the spread of multidrug resistance.

Resistance induction potency assessment of antibiotic production wastewater and associated resistome shaping mechanisms

Antibiotic production wastewater (APW) contains multiple substances known to select for and facilitate horizontal gene transfer (HGT) of antibiotic resistance genes (ARGs); however, whether these substances can induce the evolution of antibiotic resistance in real wastewater environments and the fate of such resistance induction potency during the treatment process are poorly understood, nor are its relationship with wastewater quality parameters and function in shaping the antibiotic resistome. In this study, the impacts of filter-sterilized APW and municipal wastewater on the resistance selection of Escherichia coli and the transfer dynamics of conjugative RP4 plasmid-borne ARGs across indigenous sludge communities were evaluated. The resistance development and transfer processes were accelerated in APW owing to enhanced growth inhibition, oxidative stress, and membrane permeability, with antibiotic concentrations much lower than their minimum inhibition concentrations. The effects were reduced simultaneously with the removal of COD and NH3N, but APW effluents still exhibited significant resistance induction potency with wastewater quality parameters meeting discharge standards. In contrast, municipal wastewater did not result in any detectable changes. Based on the metagenomic assembly and binning, stronger resistance induction potency in the antibiotic production wastewater treatment plant endowed indigenous sludge and effluent with greater accumulation, genetic mobility, and pathogenic accessibility of ARGs than in the municipal wastewater treatment plant. Antibiotic resistome assembly was determined primarily by deterministic processes, driven jointly by resistance induction potency, mobilome variance, and microbiome shifts. These results provide novel insights into the application of bioassays to comprehensively evaluate the antibiotic resistance induction effects of APW and their relationships with the resistome to manage risks during the treatment process.

Bacterium-Phage Interactions Enhance Biofilm Resilience during Membrane Filtration Biofouling under Oxidative and Hydraulic Stresses

Microbial interactions on membrane surfaces can facilitate biofilm formation and biofouling, which poses a significant challenge for pressure-driven membrane filtration systems. This multiomics study investigates the adaptive responses of bacterium-phage interactions under varying oxidative and hydraulic stress during membrane backwashing and their biological contributions to biofouling. Oxidative and hydraulic stress distinctly shaped bacteria and phage diversity and community composition. Under moderate oxidative backwashing (300 ppm of NaClO), diversity was maintained, with increased antioxidant enzyme activities, extracellular polymeric substance (EPS) production, and quorum sensing (QS) signaling, promoting bacterial resilience and biofilm formation. In contrast, excessive oxidative stress (600 ppm of NaClO) reduced bacteria and phage diversity, disrupted antioxidant responses, and increased microbial sensitivity. Hydraulic stress predominantly influenced viral diversity and co-occurrence network topology, favoring the expansion of broad host-range phages and lysogenic lifestyles under combined stresses. Phage-bacterium interaction analyses highlighted phages' adaptive preferences for hosts with high network centrality and broad ecological niches, which enhanced microbial interactions and resilience. Transcriptomic profiling demonstrated the early enrichment of genes associated with energy metabolism, ROS detoxification, and biofilm formation, followed by stabilization as biofilms matured. Phage-encoded auxiliary metabolic genes were involved in DNA repair, QS, and EPS biosynthesis, contributing to microbial adaptation through oxidative stress resistance and biofilm stabilization. Overall, these findings provide mechanistic insights into biofouling dynamics and highlight the need to optimize chlorine dosing to prevent suboptimal levels of microbial adaptation and biofouling.

Biogeography and ecological functions of underestimated CPR and DPANN in acid mine drainage sediments

Recent genomic surveys have uncovered candidate phyla radiation (CPR) bacteria and DPANN archaea as major microbial dark matter lineages in various anoxic habitats. Despite their extraordinary diversity, the biogeographic patterns and ecological implications of these ultra-small and putatively symbiotic microorganisms have remained elusive. Here, we performed metagenomic sequencing on 90 geochemically diverse acid mine drainage sediments sampled across southeast China and recovered 282 CPR and 189 DPANN nonredundant metagenome-assembled genomes, which collectively account for up to 28.6% and 31.2% of the indigenous prokaryotic communities, respectively. We found that, remarkably, geographic distance represents the primary factor driving the large-scale ecological distribution of both CPR and DPANN organisms, followed by pH and Fe. Although both groups might be capable of iron reduction through a flavin-based extracellular electron transfer mechanism, significant differences are found in their metabolic capabilities (with complex carbon degradation and chitin degradation being more prevalent in CPR whereas fermentation and acetate production being enriched in DPANN), indicating potential niche differentiation. Predicted hosts are mainly Acidobacteriota, Bacteroidota, and Proteobacteria for CPR and Thermoplasmatota for DPANN, and extensive, unbalanced metabolic exchanges between these symbionts and putative hosts are displayed. Together, our results provide initial insights into the complex interplays between the two lineages and their physicochemical environments and host populations at a large geographic scale.IMPORTANCECandidate phyla radiation (CPR) bacteria and DPANN archaea constitute a significant fraction of Earth's prokaryotic diversity. Despite their ubiquity and abundance, especially in anoxic habitats, we know little about the community patterns and ecological drivers of these ultra-small, putatively episymbiotic microorganisms across geographic ranges. This study is facilitated by a large collection of CPR and DPANN metagenome-assembled genomes recovered from the metagenomes of 90 sediments sampled from geochemically diverse acid mine drainage (AMD) environments across southeast China. Our comprehensive analyses have allowed first insights into the biogeographic patterns and functional differentiation of these major enigmatic prokaryotic groups in the AMD model system.

Viral diversity and auxiliary metabolic genes in rare earth element mine drainage in South China

In extreme environments, viruses play a crucial role in regulating the structure and metabolic activities of microbial communities, thereby impacting the overall biogeochemical cycles. Previous research found that rare earth element acid mine drainage (REE-AMD) harbors a wide array of microbial species. However, our understanding of the viruses that infect these microorganisms remains limited. In this study, we utilized metagenomic analysis to explore the viral diversity, interactions between viruses and their hosts, as well as the viruses encoded auxiliary metabolic genes (AMGs) within REE-AMD. The results demonstrated that viral communities showed increased diversity with REEs pollution. Furthermore, AMGs exhibited habitat and host specificity. Viruses in water samples contaminated with REEs tended to encode AMGs related to cellular metabolic processes and stress responses to protect their hosts. In contrast, viruses in sediment samples were more likely to encode AMGs associated with nutrient competition, thereby expanding the ecological niches of hosts and viruses. Viruses would carry more AMGs from the dominant prokaryotes. Additionally, under REEs stress, viruses encode a greater number of carbon- and sulfur-related AMGs, influencing the carbon and sulfur cycles of microorganisms in REE-AMD. Overall, our study provides a first systematic characterization of the viral community in REE-AMD, which is crucial for understanding the intricate interactions among viruses, their hosts, and the surrounding environment.

Metagenomic Insights in Antimicrobial Resistance Threats in Sludge from Aerobic and Anaerobic Membrane Bioreactors

Sludge is a biohazardous solid waste that is produced during wastewater treatment. It contains antibiotic resistance genes (ARGs) that pose significant antimicrobial resistance (AMR) threats. Herein, aerobic and anaerobic membrane bioreactors (AeMBRs and AnMBRs, respectively) were compared in terms of the volume of waste sludge generated by them, the presence of ARGs in the sludge, and the potential for horizontal gene transfer (HGT) events using metagenomics to determine which treatment process can better address AMR concerns associated with the generation of waste sludge. The estimated abundance of ARGs in the suspended sludge generated by the AnMBR per treated volume is, on average, 5-55 times lower than that of sludge generated by the AeMBR. Additionally, the ratio of potential HGT in the two independent runs was lower in the anaerobic sludge (0.6 and 0.9) compared with that in the aerobic sludge (2.4 and 1.6). The AnMBR sludge exhibited reduced HGT of ARGs involving potential opportunistic pathogens (0.09) compared with the AeMBR sludge (0.27). Conversely, the AeMBR sludge displayed higher diversity and more transfer events, encompassing genes that confer resistance to quinolones, rifamycin, multidrug, aminoglycosides, and tetracycline. A significant portion of these ARGs were transferred to Burkholderia sp. By contrast, the AnMBR showed a lower abundance of mobile genetic elements associated with conjugation and exhibited less favorable conditions for natural transformation. Our findings suggest that the risk of potential HGT to opportunistic pathogens is greater in the AeMBR sludge than in AnMBR sludge.

Current state and future prospects of Horizontal Gene Transfer detection

Artificial intelligence (AI) has been shown to be beneficial in a wide range of bioinformatics applications. Horizontal Gene Transfer (HGT) is a driving force of evolutionary changes in prokaryotes. It is widely recognized that it contributes to the emergence of antimicrobial resistance (AMR), which poses a particularly serious threat to public health. Many computational approaches have been developed to study and detect HGT. However, the application of AI in this field has not been investigated. In this work, we conducted a review to provide information on the current trend of existing computational approaches for detecting HGT and to decipher the use of AI in this field. Here, we show a growing interest in HGT detection, characterized by a surge in the number of computational approaches, including AI-based approaches, in recent years. We organize existing computational approaches into a hierarchical structure of computational groups based on their computational methods and show how each computational group evolved. We make recommendations and discuss the challenges of HGT detection in general and the adoption of AI in particular. Moreover, we provide future directions for the field of HGT detection.

Unveiling the role of aeration tanks in the emission and enrichment of airborne antibiotic resistance genes in a wastewater treatment plant

The aeration tanks in wastewater treatment plants (WWTPs) are important sources of airborne antibiotic resistance genes (ARGs) and antibiotic-resistant bacteria (ARB) due to bubble bursts at the air-liquid interface. This study employed an integrated metagenomic workflow, encompassing reference-based, assembly-based, and binning-based modules, to investigate resistomes in a WWTP in northern Taiwan. The role of aeration tanks in emitting airborne ARGs and their associated risks was assessed. The findings revealed a strong similarity between the ARG profiles in aeration tank wastewater and surrounding PM2.5, indicating atmospheric transmission of ARGs. Notably, the ARG level in PM2.5 (0.83 ± 0.11 ARGs/cell) was 59.6 % higher than in wastewater (0.52 ± 0.01 ARGs/cell). The assembly-based analysis showed that foam-forming bacteria such as Mycobacterium and Gordonia dominated ARGs-carrying contigs in PM2.5, suggesting that higher atomization capabilities of ARB contribute to airborne ARG prevalence. Furthermore, a significant proportion of stress response genes and increasing efflux pump resistance (122.6 %) in PM2.5 imply that mechanical forces during aerosolization and harsh atmospheric conditions select for airborne ARB capable of overcoming stress induced by dramatic environmental changes. Overall, the study indicates that ARG risk is intensified in PM2.5 due to their abundance, mobility, and pathogenicity. In conclusion, aeration tanks not only emit airborne ARGs but also cause an unexpected enrichment effect and exposure risk during aeration, highlighting the critical water-to-air transmission route of ARGs in WWTPs.

Bacterium-Phage Symbiosis Facilitates the Enrichment of Bacterial Pathogens and Antibiotic-Resistant Bacteria in the Plastisphere

The plastisphere, defined as the ecological niche for microbial colonization of plastic debris, has been recognized as a hotspot of pathogenic and antibiotic-resistant bacteria. However, the interactions between bacteria and phages facilitated by the plastisphere, as well as their impact on microbial risks to public health, remain unclear. Here, we analyzed public metagenomic data from 180 plastisphere and environmental samples, stemming from four different habitats and two plastic types (biodegradable and nonbiodegradable plastics) and obtained 611 nonredundant metagenome-assembled genomes (MAGs) and 4061 nonredundant phage contigs. The plastisphere phage community exhibited decreased diversity and virulent proportion compared to those found in environments. Indexes of phage-host interaction networks indicated significant associations of phages with pathogenic and antibiotic-resistant bacteria (ARB), particularly for biodegradable plastics. Known phage-encoded auxiliary metabolic genes (AMGs) were involved in nutrient metabolism, antibiotic production, quorum sensing, and biofilm formation in the plastisphere, which contributed to enhanced competition and survival of pathogens and ARB hosts. Phages also carried transcriptionally active virulence factor genes (VFGs) and antibiotic resistance genes (ARGs), and could mediate their horizontal transfer in microbial communities. Overall, these discoveries suggest that plastisphere phages form symbiotic relationships with their hosts, and that phages encoding AMGs and mediating horizontal gene transfer (HGT) could increase the source of pathogens and antibiotic resistance from the plastisphere.

Structure and metabolic function of spatiotemporal pit mud microbiome

BACKGROUND

Pit mud (PM) hosts diverse microbial communities, which serve as a medium to impart flavor and quality to Baijiu and exhibit long-term tolerance to ethanol and acids, resulting in a unique ecosystem. However, the ecology and metabolic functions of PM remain poorly understood, as many taxa in PM represent largely novel lineages. In this study, we used a combination of metagenomic analysis and chemical derivatization LC-MS analysis to provide a comprehensive overview of microbial community structure, metabolic function, phylogeny, horizontal gene transfer, and the relationship with carboxyl compounds in spatiotemporal PM samples.

RESULTS

Our findings revealed three distinct stages in the spatiotemporal changes of prokaryotic communities in PM: an initial phase dominated by Lactobacillus, a transitional phase, and a final state of equilibrium. Significant variations in α- and β-diversity were observed across different spatial and temporal PM samples. We identified 178 medium- and high-quality non-redundant metagenome-assembled genomes (MAGs), and constructed their phylogenetic tree, depicting their roles in the carbon, nitrogen, and sulfur cycles. The Wood-Ljungdahl pathway and reverse TCA cycle were identified as the main carbon fixation mechanisms, with both hydrogenotrophic and aceticlastic methanogens playing a major role in methane production, and methylotrophic pathway observed in older PM. Furthermore, we identified relationships between prokaryotes and 29 carboxyl metabolites, including medium- and long-chain fatty acids. Horizontal gene transfer (HGT) was widespread in PM, particularly among clostridia, Bacteroidota, Bacilli, and Euryarchaeota, and was shown to play critical roles in fermentation dynamics, carbon fixation, methane production, and nitrogen and sulfur metabolism.

CONCLUSION

Our study provides new insights into the evolution and function of spatiotemporal PM, as well as its interactions with carboxyl metabolites. Lactobacillus dominated in new PM, while methanogens and clostridia were predominant in older or deeper PM layers. The three distinct stages of prokaryotic community development in PM and HGT played critical roles in metabolic function of spatiotemporal PM. Furthermore, this study highlights the importance of α-diversity, β-diversity, methanogens, and Clostridium as useful indicators for assessing PM quality in the production of high-quality Baijiu.

GutMetaNet: an integrated database for exploring horizontal gene transfer and functional redundancy in the human gut microbiome

Metagenomic studies have revealed the critical roles of complex microbial interactions, including horizontal gene transfer (HGT) and functional redundancy (FR), in shaping the gut microbiome's functional capacity and resilience. However, the lack of comprehensive data integration and systematic analysis approaches has limited the in-depth exploration of HGT and FR dynamics across large-scale gut microbiome datasets. To address this gap, we present GutMetaNet (https://gutmetanet.deepomics.org/), a first-of-its-kind database integrating extensive human gut microbiome data with comprehensive HGT and FR analyses. GutMetaNet contains 21 567 human gut metagenome samples with whole-genome shotgun sequencing data related to various health conditions. Through systematic analysis, we have characterized the taxonomic profiles and FR profiles, and identified 14 636 HGT events using a shared reference genome database across the collected samples. These HGT events have been curated into 8049 clusters, which are annotated with categorized mobile genetic elements, including transposons, prophages, integrative mobilizable elements, genomic islands, integrative conjugative elements and group II introns. Additionally, GutMetaNet incorporates automated analyses and visualizations for the HGT events and FR, serving as an efficient platform for in-depth exploration of the interactions among gut microbiome taxa and their implications for human health.

Profiling lateral gene transfer events in the human microbiome using WAAFLE

Lateral gene transfer (LGT), also known as horizontal gene transfer, facilitates genomic diversification in microbial populations. While previous work has surveyed LGT in human-associated microbial isolate genomes, the landscape of LGT arising in personal microbiomes is not well understood, as there are no widely adopted methods to characterize LGT from complex communities. Here we developed, benchmarked and validated a computational algorithm (WAAFLE or Workflow to Annotate Assemblies and Find LGT Events) to profile LGT from assembled metagenomes. WAAFLE prioritizes specificity while maintaining high sensitivity for intergenus LGT. Applying WAAFLE to >2,000 human metagenomes from diverse body sites, we identified >100,000 high-confidence previously uncharacterized LGT (~2 per microbial genome-equivalent). These were enriched for mobile elements, as well as restriction-modification functions associated with the destruction of foreign DNA. LGT frequency was influenced by biogeography, phylogenetic similarity of involved pairs (for example, Fusobacterium periodonticum and F. nucleatum) and donor abundance. These forces manifest as networks in which hub taxa donate unequally with phylogenetic neighbours. Our findings suggest that human microbiome LGT may be more ubiquitous than previously described.

Anaerobic Degradation of Aromatic and Aliphatic Biodegradable Plastics: Potential Mechanisms and Pathways

Biodegradable plastics (BDPs) have been widely used as substitutes for traditional plastics, and their environmental fate is a subject of intense research interest. Compared with the aerobic degradation of BDPs, their biodegradability under anaerobic conditions in environmental engineering systems remains poorly understood. This study aimed to investigate the degradability of BDPs composed of poly(butylene adipate-co-terephthalate) (PBAT), poly(lactide acid) (PLA), and their blends, and explore the mechanism underlying their microbial degradation under conditions of anaerobic digestion (AD). The BDPs readily depolymerized under thermophilic conditions but were hydrolyzed at a slow rate under conditions of mesophilic AD. After 45 days of thermophilic AD, a decrease in the molecular weight and significant increase in the production of methane and carbon dioxide production were observed. Network and metagenomics analyses identified AD as reservoirs of plastic-degrading bacteria that produce multiple plastic-degrading enzymes. PETase was identified as the most abundant plastic-degrading enzyme. A potential pathway for the anaerobic biodegradation of BDPs was proposed herein. The polymers of high molecular weight were subjected to abiotic hydrolysis to form oligomers and monomers, enabling subsequent microbial hydrolysis and acetogenesis. Ultimately, complete degradation was achieved predominantly via the pathway involved in the conversion of acetic acid to methane. These findings provide novel insight into the mechanism underlying the anaerobic degradation of BDPs and the microbial resources crucial for the efficient degradation of BDPs.

Metagenomic insights into microalgae-bacterium-virus interactions and viral functions in phycosphere facing environmental fluctuations

Despite the ecological and biotechnological significance of microalgae-bacterium symbionts, the response of host-virus interactions to external environmental fluctuations and the role of viruses in phycosphere remain largely unexplored. Herein, we employed algal-bacterial granular sludge (ABGS) with varying light intensity and organic carbon loading to investigate the mechanisms of microalgae-bacterium-virus symbionts in response to environmental fluctuations. Metagenomics revealed that enhanced light intensity decreased the diversity of microalgae, so did the diversity of symbiotic bacteria and viruses. As carbon sources decreased, bacteria prompted horizontal gene transfer in phycosphere by 12.76 %-157.40 %, increased the proportion of oligotrophs as keystone species (0.00 % vs 14.29 %) as well as viruses using oligotrophs as hosts (18.52 % vs 25.00 %). Furthermore, virus-carried auxiliary metabolic genes (AMGs) and biosynthetic gene clusters (BGCs) encoding vitamin B12 synthesis (e.g., cobS), antioxidation (e.g., queC), and microbial aggregation (e.g., cysE). Additionally, phylogenetic and similarity analysis further revealed the evolutionary origin and potential horizontal transfer of the AMGs and BGCs, which could potentially enhance the adaptability of bacteria and eukaryotic microalgae. Overall, our research demonstrates that environmental fluctuations have cascading effects on the microalgae-bacteria-virus interactions, and emphasizes the important role of viruses in maintaining the stability of the phycosphere symbiotic community.

Genome-Centric Metagenomics Insights into the Plastisphere-Driven Natural Degradation Characteristics and Mechanism of Biodegradable Plastics in Aquatic Environments

Biodegradable plastics (BPs) are pervasively available as alternatives to traditional plastics, but their natural degradation characteristics and microbial-driven degradation mechanisms are poorly understood, especially in aquatic environments, the primary sink of plastic debris. Herein, the three-month dynamic degradation process of BPs (the copolymer of poly(butylene adipate-co-terephthalate) and polylactic acid (PLA) (PBAT/PLA) and single PLA) in a natural aquatic environment was investigated, with nonbiodegradable plastics polyvinyl chloride, polypropylene, and polystyrene as controls. PBAT/PLA showed the weight loss of 47.4% at 50 days and severe fragmentation within two months, but no significant decay for other plastics. The significant increase in the specific surface area and roughness and the weakening of hydrophobicity within the first month promoted microbial attachment to the PBAT/PLA surface. Then, a complete microbial succession occurred, including biofilm formation, maturation, and dispersion. Metagenomic analysis indicated that plastispheres selectively enriched degraders. Based on the functional genes involved in BPs degradation, a total of 16 high-quality metagenome-assembled genomes of degraders (mainly Burkholderiaceae) were recovered from the PBAT/PLA plastisphere. These microbes showed the greatest degrading potential at the biofilm maturation stage and executed the functions by PLA_depolymerase, polyesterase, hydrolase, and esterase. These findings will enhance understanding of BPs' environmental behavior and microbial roles on plastic degradation.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Comparative genomics-based insights into Xanthomonas indica, a non-pathogenic species of healthy rice microbiome with bioprotection function

Xanthomonas species are major pathogens of plants and have been studied extensively. There is increasing recognition of the importance of non-pathogenic species within the same genus. With this came the need to understand the genomic and functional diversity of non-pathogenic Xanthomonas (NPX) at the species and strain level. This study reports isolation and investigation into the genomic diversity and variation in NPX isolates, chiefly Xanthomonas indica, a newly discovered NPX species from rice. The study establishes the relationship of X. indica strains within clade I of Xanthomonads with another NPX species, X. sontii, also associated with rice seeds. Identification of highly diverse strains, open-pan genome, and systematic hyper-variation at the lipopolysaccharide biosynthetic locus when compared to pathogenic Xanthomonas indicates the acquisition of new functions for adaptation. Furthermore, comparative genomics studies established the absence of major virulence genes such as type III secretion system and effectors, which are present in the pathogens, and the presence of a known bacterial-killing type IV secretion system (X-T4SS). The diverse non-pathogenic strains of X. indica and X. sontii were found to protect rice from bacterial leaf blight pathogen, X. oryzae pv. oryzae (Xoo). The absence of phenotype of an X-T4SS mutant suggests redundancy in the genetic basis of the mechanisms involved in the bioprotection function, which may include multiple genetic loci, such as putative bacteriocin-encoding gene clusters and involvement of other factors such as nutrient and niche competition apart from induction of innate immunity through shared microbial-associated molecular patterns. The rice-NPX community and its pathogenic counterpart can be a promising model for understanding plant-microbe-microbiome interaction studies.IMPORTANCEThe Xanthomonas group of bacteria is known for its characteristic lifestyle as a phytopathogen. However, the discovery of non-pathogenic Xanthomonas (NPX) species is a major shift in understanding this group of bacteria. Multi-strain, in-depth genomic, evolutionary and functional studies on each of these NPX species are still lacking. This study on diverse non-pathogenic strains provides novel insights into genome diversity, dynamics, and evolutionary trends of NPX species from rice microbiome apart from its relationship with other relatives that form a sub-clade. Interestingly, we also uncovered that NPX species protect rice from pathogenic Xanthomonas species. The plant protection property shows their importance as a part of a healthy plant microbiome. Furthermore, finding an open pan-genome and large-scale variation at lipopolysaccharide biosynthetic locus indicates a significant role of the NPX community in host adaptation. The findings and high-quality genomic resources of NPX species and the strains will allow further systematic molecular and host-associated microbial community studies for plant health.

Decoupling of strain- and intrastrain-level interactions of microbiomes in a sponge holobiont

Holobionts are highly organized assemblages of eukaryotic hosts, cellular microbial symbionts, and viruses, whose interactions and evolution involve complex biological processes. It is largely unknown which specific determinants drive similarity or individuality in genetic diversity between holobionts. Here, we combine short- and long-read sequencing and DNA-proximity-linkage technologies to investigate intraspecific diversity of the microbiomes, including host-resolved viruses, in individuals of a model marine sponge. We find strong impacts of the sponge host and the cellular hosts of viruses on strain-level organization of the holobiont, whereas substantial overlap in nucleotide diversity between holobionts suggests frequent exchanges of microbial cells and viruses at intrastrain level in the local sponge population. Immune-evasive arms races likely restricted virus-host co-evolution at the intrastrain level, generated holobiont-specific genome variations, and linked virus-host genetics through recombination. Our work shows that a decoupling of strain- and intrastrain-level interactions is a key factor in the genetic diversification of holobionts.

Long-read sequencing reveals extensive gut phageome structural variations driven by genetic exchange with bacterial hosts

Genetic variations are instrumental for unraveling phage evolution and deciphering their functional implications. Here, we explore the underlying fine-scale genetic variations in the gut phageome, especially structural variations (SVs). By using virome-enriched long-read metagenomic sequencing across 91 individuals, we identified a total of 14,438 nonredundant phage SVs and revealed their prevalence within the human gut phageome. These SVs are mainly enriched in genes involved in recombination, DNA methylation, and antibiotic resistance. Notably, a substantial fraction of phage SV sequences share close homology with bacterial fragments, with most SVs enriched for horizontal gene transfer (HGT) mechanism. Further investigations showed that these SV sequences were genetic exchanged between specific phage-bacteria pairs, particularly between phages and their respective bacterial hosts. Temperate phages exhibit a higher frequency of genetic exchange with bacterial chromosomes and then virulent phages. Collectively, our findings provide insights into the genetic landscape of the human gut phageome.

Enhancing insights into diseases through horizontal gene transfer event detection from gut microbiome

Horizontal gene transfer (HGT) phenomena pervade the gut microbiome and significantly impact human health. Yet, no current method can accurately identify complete HGT events, including the transferred sequence and the associated deletion and insertion breakpoints from shotgun metagenomic data. Here, we develop LocalHGT, which facilitates the reliable and swift detection of complete HGT events from shotgun metagenomic data, delivering an accuracy of 99.4%-verified by Nanopore data-across 200 gut microbiome samples, and achieving an average F1 score of 0.99 on 100 simulated data. LocalHGT enables a systematic characterization of HGT events within the human gut microbiome across 2098 samples, revealing that multiple recipient genome sites can become targets of a transferred sequence, microhomology is enriched in HGT breakpoint junctions (P-value = 3.3e-58), and HGTs can function as host-specific fingerprints indicated by the significantly higher HGT similarity of intra-personal temporal samples than inter-personal samples (P-value = 4.3e-303). Crucially, HGTs showed potential contributions to colorectal cancer (CRC) and acute diarrhoea, as evidenced by the enrichment of the butyrate metabolism pathway (P-value = 3.8e-17) and the shigellosis pathway (P-value = 5.9e-13) in the respective associated HGTs. Furthermore, differential HGTs demonstrated promise as biomarkers for predicting various diseases. Integrating HGTs into a CRC prediction model achieved an AUC of 0.87.

Global survey of hgcA-carrying genomes in marine and freshwater sediments: Insights into mercury methylation processes

Mercury (Hg) methylation is a microbially mediated process that produces methylmercury (MeHg), a bioaccumulative neurotoxin. A highly conserved gene pair, hgcAB, is required for Hg methylation, which provides a basis for identifying Hg methylators and evaluating their genomic composition. In this study, we conducted a large-scale omics analysis in which 281 metagenomic freshwater and marine sediment samples from 46 geographic locations across the globe were queried. Specific objectives were to examine the prevalence of Hg methylators, to identify horizontal gene transfer (HGT) events involving hgcAB within Hg methylator communities, and to identify associations between hgcAB and microbial biochemical functions/genes. Hg methylators from the phyla Desulfobacterota and Bacteroidota were dominant in both freshwater and marine sediments while Firmicutes and methanogens belonging to Euryarchaeota were identified only in freshwater sediments. Novel Hg methylators were found in the Phycisphaerae and Planctomycetia classes within the phylum Planctomycetota, including potential hgcA-carrying anammox metagenome-assembled genomes (MAGs) from Candidatus Brocadiia. HGT of hgcA and hgcB were identified in both freshwater and marine methylator communities. Spearman's correlation analysis of methylator genomes suggested that in addition to sulfide, thiosulfate, sulfite, and ammonia may be important parameters for Hg methylation processes in sediments. Overall, our results indicated that the biochemical drivers of Hg methylation vary between marine and freshwater sites, lending insight into the influence of environmental perturbances, such as a changing climate, on Hg methylation processes.

Reference-free structural variant detection in microbiomes via long-read co-assembly graphs

MOTIVATION

The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining.

RESULTS

We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux.

AVAILABILITY AND IMPLEMENTATION

rhea is open source and available at: https://github.com/treangenlab/rhea.

Toward nitrogen recovery: Co-cultivation of microalgae and bacteria enhances the production of high-value nitrogen-rich cyanophycin

The algal-bacterial wastewater treatment process has been proven to be highly efficient in removing nutrients and recovering nitrogen (N). However, the recovery of the valuable N-rich biopolymer, cyanophycin, remains limited. This research explored the synthesis mechanism and recovery potential of cyanophycin within two algal-bacterial symbiotic reactors. The findings reveal that the synergy between algae and bacteria enhances the removal of N and phosphorus. The crude contents of cyanophycin in the algal-bacterial consortia reached 115 and 124 mg/g of mixed liquor suspended solids (MLSS), respectively, showing an increase of 11.7 %-20.4 % (p < 0.001) compared with conventional activated sludge. Among the 170 metagenome-assembled genomes (MAGs) analyzed, 50 were capable of synthesizing cyanophycin, indicating that cyanophycin producers are common in algal-bacterial systems. The compositions of cyanophycin producers in the two algal-bacterial reactors were affected by different lighting initiation time. The study identified two intracellular synthesis pathways for cyanophycin. Approximately 36 MAGs can synthesize cyanophycin de novo using ammonium and glucose, while the remaining 14 MAGs require exogenous arginine for production. Notably, several MAGs with high abundance are capable of assimilating both nitrate and ammonium into cyanophycin, demonstrating a robust N utilization capability. This research also marks the first identification of potential horizontal gene transfer of the cyanophycin synthase encoding gene (cphA) within the wastewater microbial community. This suggests that the spread of cphA could expand the population of cyanophycin producers. The study offers new insights into recycling the high-value N-rich biopolymer cyanophycin, contributing to the advancement of wastewater resource utilization.

New insights into functional divergence and adaptive evolution of uncultured bacteria in anammox community by complete genome-centric analysis

Anaerobic ammonium-oxidation (anammox) bacteria play a crucial role in global nitrogen cycling and wastewater nitrogen removal, but they share symbiotic relationships with various other microorganisms. Functional divergence and adaptive evolution of uncultured bacteria in anammox community remain underexplored. Although shotgun metagenomics based on short reads has been widely used in anammox research, metagenome-assembled genomes (MAGs) are often discontinuous and highly contaminated, which limits in-depth analyses of anammox communities. Here, for the first time, we performed Pacific Biosciences high-fidelity (HiFi) long-read sequencing on the anammox granule sludge sample from a lab-scale bioreactor, and obtained 30 accurate and complete metagenome-assembled genomes (cMAGs). These cMAGs were obtained by selecting high-quality circular contigs from initial assemblies of long reads generated by HiFi sequencing, eliminating the need for Illumina short reads, binning, and reassembly. One new anammox species affiliated with Candidatus Jettenia and three species affiliated with novel families were found in this anammox community. cMAG-centric analysis revealed functional divergence in general and nitrogen metabolism among the anammox community members, and they might adopt a cross-feeding strategy in organic matter, cofactors, and vitamins. Furthermore, we identified 63 mobile genetic elements (MGEs) and 50 putative horizontal gene transfer (HGT) events within these cMAGs. The results suggest that HGT events and MGEs related to phage and integration or excision, particularly transposons containing tnpA in anammox bacteria, might play important roles in the adaptive evolution of this anammox community. The cMAGs generated in the present study could be used to establish of a comprehensive database for anammox bacteria and associated microorganisms. These findings highlight the advantages of HiFi sequencing for the studies of complex mixed cultures and advance the understanding of anammox communities.

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.

Tire particles and its leachates: Impact on antibiotic resistance genes in coastal sediments

Tire particles (TPs), a significant group of microplastics, can be discharged into the coastal environments in various ways. However, our understanding of how TPs impact the antibiotic resistance and pathogenic risks of microorganisms in coastal sediments remains limited. In this study, we used metagenomics to investigate how TPs and their leachates could affect the prevalence of antibiotic resistance genes (ARGs), virulence factor genes (VFGs), and their potential risks to the living creatures such as soil invertebrates and microorganisms in the coastal sediments. We discovered that TP addition significantly increased the abundance and diversity of ARGs and VFGs in coastal sediments, with raw TPs displayed higher impacts than TP leachates and TPs after leaching on ARGs and VFGs. With increasing TP exposure concentrations, the co-occurrence frequency of ARGs and mobile genetic elements (MGEs) in the same contig also increased, suggesting that TPs could enhance the dispersal risk of ARGs. Our metagenome-based binning analysis further revealed that exposure to TPs increased the abundance of potentially pathogenic antibiotic-resistant bacteria (PARB). In addition, chemical additives of TP leachates (e.g., Zn and N-cyclohexylformamide) significantly affected the changes of ARGs in the pore water. In summary, our study provides novel insights into the adverse effects of TP pollutions on aggravating the dissemination and pathogenic risks of ARGs and PARB in the coastal environment.

Horizontal gene transfer in activated sludge enhances microbial antimicrobial resistance and virulence

Activated sludge (AS) plays a vital role in removing organic pollutants and nutrients from wastewater. However, the risks posed by horizontal gene transfer (HGT) between bacteria in AS are still unclear. Here, a total of 478 high-quality non-redundant metagenome-assembled genomes (MAGs) were obtained. >50 % and 5 % of MAGs were involved in at least one HGT and recent HGT, respectively. Most of the transfers (82.4 %) of antimicrobial resistance genes (ARGs) occurred among the classes of Alphaproteobacteria and Gammaproteobacteria. The bacteria involved in the transfers of virulence factor genes (VFGs) mainly include Alphaproteobacteria (42.3 %), Bacteroidia (19.2 %), and Gammaproteobacteria (11.5 %). Moreover, the number of ARGs and VFGs in the classes of Alphaproteobacteria and Gammaproteobacteria was higher than that in other bacteria (P < 0.001). Mobile genetic elements were important contributors to ARGs and VFGs in AS bacteria. These results have implications for the management of antimicrobial resistance and virulence in activated sludge microorganisms.

Eco-evolutionary dynamics of gut phageome in wild gibbons (Hoolock tianxing) with seasonal diet variations

It has been extensively studied that the gut microbiome provides animals flexibility to adapt to food variability. Yet, how gut phageome responds to diet variation of wild animals remains unexplored. Here, we analyze the eco-evolutionary dynamics of gut phageome in six wild gibbons (Hoolock tianxing) by collecting individually-resolved fresh fecal samples and parallel feeding behavior data for 15 consecutive months. Application of complementary viral and microbial metagenomics recovers 39,198 virulent and temperate phage genomes from the feces. Hierarchical cluster analyses show remarkable seasonal diet variations in gibbons. From high-fruit to high-leaf feeding period, the abundances of phage populations are seasonally fluctuated, especially driven by the increased abundance of virulent phages that kill the Lachnospiraceae hosts, and a decreased abundance of temperate phages that piggyback the Bacteroidaceae hosts. Functional profiling reveals an enrichment through horizontal gene transfers of toxin-antitoxin genes on temperate phage genomes in high-leaf season, potentially conferring benefits to their prokaryotic hosts. The phage-host ecological dynamics are driven by the coevolutionary processes which select for tail fiber and DNA primase genes on virulent and temperate phage genomes, respectively. Our results highlight complex phageome-microbiome interactions as a key feature of the gibbon gut microbial ecosystem responding to the seasonal diet.

The microbiome, resistome, and their co-evolution in sewage at a hospital for infectious diseases in Shanghai, China

The emergence of antibiotic-resistant bacteria (ARB) caused by the overuse of antibiotics severely threatens human health. Hospital sewage may be a key transmission hub for ARB. However, the complex link between the microbiome and resistomeresistance in hospital sewage remains unclear. In this study, metagenomic assembly and binning methods were used to investigate the microbial community, resistome, and association of antibiotic resistance genes (ARGs) with ARB in sewage from 10 representative sites (outpatient building, surgery building, internal medicine buildings [IMB1-4], staff dormitory, laboratory animal building, tuberculosis building [TBB], and hospital wastewater treatment plant) of a hospital in Shanghai from June 2021 to February 2022. A total of 252 ARG subtypes, belonging to 17 antibiotic classes, were identified. The relative abundance of KPC-2 was higher at IMBs and TBB than at other sites. Of the ARG-carrying contigs, 47.3%-62.6% were associated with mobile genetic elements, and the proportion of plasmid-associated ARGs was significantly higher than that of chromosome-associated ARGs. Although a similar microbiome composition was shared, certain bacteria were enriched at different sites. Potential pathogens Enterococcus B faecium and Klebsiella pneumoniae were primarily enriched in IMB2 and IMB4, respectively. The same ARGs were identified in diverse bacterial hosts (especially pathogenic bacteria), and accordingly, the latter possessed multiple ARGs. Furthermore, gene flow was frequently observed in the sewage of different buildings. The results provide crucial information on the characterization profiles of resistomes in hospital sewage in Shanghai.IMPORTANCEEnvironmental antibiotic resistance genes (ARGs) play a critical role in the emergence and spread of antimicrobial resistance, which poses a global health threat. Wastewater from healthcare facilities serves as a significant reservoir for ARGs. Here, we characterized the microbial community along with the resistome (comprising all antibiotic resistance genes) in wastewater from a specialized hospital for infectious diseases in Shanghai. Potential pathogenic bacteria (e.g., Escherichia coli, Pseudomonas aeruginosa, Klebsiella pneumoniae, Enterococcus B faecium) were frequently detected in hospital wastewater and carried multiple ARGs. A complex link between microbiome and resistome was observed in the wastewater of this hospital. The monitoring of ARGs and antibiotic-resistant bacteria (ARB) in hospital wastewater might be of great significance for preventing the spread of ARB.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

Comparative genomic analyses of Cutibacterium granulosum provide insights into genomic diversity

Cutibacterium granulosum, a commensal bacterium found on human skin, formerly known as Propionibacterium granulosum, rarely causes infections and is generally considered non-pathogenic. Recent research has revealed the transferability of the multidrug-resistant plasmid pTZC1 between C. granulosum and Cutibacterium acnes, the latter being an opportunistic pathogen in surgical site infections. However, there is a noticeable lack of research on the genome of C. granulosum, and the genetic landscape of this species remains largely uncharted. We investigated the genomic features and evolutionary structure of C. granulosum by analyzing a total of 30 Metagenome-Assembled Genomes (MAGs) and isolate genomes retrieved from public databases, as well as those generated in this study. A pan-genome of 6,077 genes was identified for C. granulosum. Remarkably, the 'cloud genes' constituted 62.38% of the pan-genome. Genes associated with mobilome: prophages, transposons [X], defense mechanisms [V] and replication, recombination and repair [L] were enriched in the cloud genome. Phylogenomic analysis revealed two distinct mono-clades, highlighting the genomic diversity of C. granulosum. The genomic diversity was further confirmed by the distribution of Average Nucleotide Identity (ANI) values. The functional profiles analysis of C. granulosum unveiled a wide range of potential Antibiotic Resistance Genes (ARGs) and virulence factors, suggesting its potential tolerance to various environmental challenges. Subtype I-E of the CRISPR-Cas system was the most abundant in these genomes, a feature also detected in C. acnes genomes. Given the widespread distribution of C. granulosum strains within skin microbiome, our findings make a substantial contribution to our broader understanding of the genetic diversity, which may open new avenues for investigating the mechanisms and treatment of conditions such as acne vulgaris.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Phylogenetic reconciliation: making the most of genomes to understand microbial ecology and evolution

In recent years, phylogenetic reconciliation has emerged as a promising approach for studying microbial ecology and evolution. The core idea is to model how gene trees evolve along a species tree and to explain differences between them via evolutionary events including gene duplications, transfers, and losses. Here, we describe how phylogenetic reconciliation provides a natural framework for studying genome evolution and highlight recent applications including ancestral gene content inference, the rooting of species trees, and the insights into metabolic evolution and ecological transitions they yield. Reconciliation analyses have elucidated the evolution of diverse microbial lineages, from Chlamydiae to Asgard archaea, shedding light on ecological adaptation, host-microbe interactions, and symbiotic relationships. However, there are many opportunities for broader application of the approach in microbiology. Continuing improvements to make reconciliation models more realistic and scalable, and integration of ecological metadata such as habitat, pH, temperature, and oxygen use offer enormous potential for understanding the rich tapestry of microbial life.

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

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

Functional prediction of proteins from the human gut archaeome

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

Proton pump inhibitors increase the risk of carbapenem-resistant Enterobacteriaceae colonization by facilitating the transfer of antibiotic resistance genes among bacteria in the gut microbiome

Carbapenem-resistant Enterobacteriaceae (CRE) pose a global health threat; however, there is still limited understanding of the risk factors and underlying mechanisms of CRE colonization in the gut microbiome. We conducted a matched case-control study involving 282 intensive care unit patients to analyze influencing covariates on CRE colonization. Subsequently, their effects on the gut microbiome were analyzed in a subset of 98 patients (47 CRE carriers and 51 non-CRE carriers) using whole metagenome sequences. The concomitant use of proton pump inhibitors (PPIs) and antibiotics was a significant risk factor for CRE colonization. The gut microbiome differed according to PPI administration, even within the CRE and non-CRE groups. Moreover, the transfer of mobile genetic elements (MGEs) harboring carbapenem resistance genes (CRGs) between bacteria was higher in the PPI-treated group than in the PPI-not-treated group among CRE carriers. The concomitant use of PPIs and antibiotics significantly alters the gut microbiome and increases the risk of CRE colonization by facilitating the transfer of CRGs among bacteria of the gut microbiome. Based on these findings, improved stewardship of PPIs as well as antibiotics can provide strategies to reduce the risk of CRE colonization, thereby potentially improving patient prognosis.

Assessing the Influence of HGT on the Evolution of Stress Responses in Microbial Communities from Shark Bay, Western Australia

Pustular microbial mats in Shark Bay, Western Australia, are modern analogs of microbial systems that colonized peritidal environments before the evolution of complex life. To understand how these microbial communities evolved to grow and metabolize in the presence of various environmental stresses, the horizontal gene transfer (HGT) detection tool, MetaCHIP, was used to identify the horizontal transfer of genes related to stress response in 83 metagenome-assembled genomes from a Shark Bay pustular mat. Subsequently, maximum-likelihood phylogenies were constructed using these genes and their most closely related homologs from other environments in order to determine the likelihood of these HGT events occurring within the pustular mat. Phylogenies of several stress-related genes-including those involved in response to osmotic stress, oxidative stress and arsenic toxicity-indicate a potentially long history of HGT events and are consistent with these transfers occurring outside of modern pustular mats. The phylogeny of a particular osmoprotectant transport gene reveals relatively recent adaptations and suggests interactions between Planctomycetota and Myxococcota within these pustular mats. Overall, HGT phylogenies support a potentially broad distribution in the relative timing of the HGT events of stress-related genes and demonstrate ongoing microbial adaptations and evolution in these pustular mat communities.

Unraveling the antibiotic resistome in backwater zones of large cascade reservoirs: Co-occurrence patterns, horizontal transfer directions and health risks

The widespread occurrence of antibiotic resistant genes (ARGs) throughout aquatic environments has raised global concerns for public health. However, the profiles and patterns of antibiotic resistome in backwater zone of cascade reservoirs, where water flow is slowed down, are still poorly understood. Here, we proposed a metagenomic analysis framework to comprehensively reveal the diversity, abundance, co-occurrence patterns and transfer direction of ARGs in cascade reservoirs system and evaluated their health risks through a procedure based on contigs. A total of 364 ARGs subtypes conferring resistance to different antibiotics classes were detected in our water samples, and the dominant ARGs (macB, bacA, vanRA, bcrA) were similar in different reservoirs. Meanwhile, the distribution of ARGs was influenced by the presence of biotic factors such as metal resistant genes (MRGs) and mobile genetic elements (MGEs), as well as abiotic factors such as dissolved oxygen (DO) and pH. Remarkably, ARGs (vanR, rosB, MexT) co-occurred with plasmids and virulence factor genes (VFGs), which can lead to the emergence and spread of highly virulent and antibiotic resistant bacteria in microbial communities. Overall, this study helps administrators to better understand the complex patterns of ARGs in backwater zones of large cascade reservoirs and provides a proper procedure for detecting the presence of high-risk of ARGs.

Pet cats may shape the antibiotic resistome of their owner's gut and living environment

BACKGROUND

Companion animals can contribute to the physical and mental health of people and often live in very close association with their owners. However, the antibiotic resistome carried by companion animals and the impact they have on their owners and living environment remain unclear. In this study, we compared the ARG profiles of cats, humans, and their living environments using metagenomic analysis to identify the core ARGs in the cat and human gut and explore the potential impact of cats on ARGs in the human gut through the environment.

RESULTS

Results showed that the abundance of ARGs in the cat gut was significantly higher than that in the human gut (P < 0.0001), with aminoglycoside and tetracycline resistance genes being the dominant ARGs in the cat gut. There was no significant difference in the abundance of total ARGs in the guts of cat owners and non-owners (P > 0.05). However, the abundance of aminoglycoside resistance genes including APH(2'')-IIa and AAC(6')-Im was significantly higher in cat owners than that in non-cat owners (P < 0.001). Also, ARG abundance was positively correlated with the frequency of cat activity in the living environment. Enterobacteriaceae was the dominant ARG host co-occurring in the cat gut, human gut, and living environment.

CONCLUSIONS

Our results show that cats may shape the living environment resistome and thus the composition of some ARGs in the human gut, highlighting the importance of companion animal environment health. Video Abstract.

Pangenome analysis provides insights into the genetic diversity, metabolic versatility, and evolution of the genus Flavobacterium

Members of the genus Flavobacterium are widely distributed and produce various polysaccharide-degrading enzymes. Many species in the genus have been isolated and characterized. However, few studies have focused on marine isolates or fish pathogens, and in-depth genomic analyses, particularly comparative analyses of isolates from different habitat types, are lacking. Here, we isolated 20 strains of the genus from various environments in South Korea and sequenced their full-length genomes. Combined with published sequence data, we examined genomic traits, evolution, environmental adaptation, and putative metabolic functions in total 187 genomes of isolated species in Flavobacterium categorized as marine, host-associated, and terrestrial including freshwater. A pangenome analysis revealed a correlation between genome size and coding or noncoding density. Flavobacterium spp. had high levels of diversity, allowing for novel gene repertories via recombination events. Defense-related genes only accounted for approximately 3% of predicted genes in all Flavobacterium genomes. While genes involved in metabolic pathways did not differ with respect to isolation source, there was substantial variation in genomic traits; in particular, the abundances of tRNAs and rRNAs were higher in the host-associdated group than in other groups. One genome in the host-associated group contained a Microviridae prophage closely related to an enterobacteria phage. The proteorhodopsin gene was only identified in four terrestrial strains isolated for this study. Furthermore, recombination events clearly influenced genomic diversity and may contribute to the response to environmental stress. These findings shed light on the high genetic variation in Flavobacterium and functional roles in diverse ecosystems as a result of their metabolic versatility. IMPORTANCE The genus Flavobacterium is a diverse group of bacteria that are found in a variety of environments. While most species of this genus are harmless and utilize organic substrates such as proteins and polysaccharides, some members may play a significant role in the cycling for organic substances within their environments. Nevertheless, little is known about the genomic dynamics and/or metabolic capacity of Flavobacterium. Here, we found that Flavobacterium species may have an open pangenome, containing a variety of diverse and novel gene repertoires. Intriguingly, we discovered that one genome (classified into host-associated group) contained a Microviridae prophage closely related to that of enterobacteria. Proteorhodopsin may be expressed under conditions of light or oxygen pressure in some strains isolated for this study. Our findings significantly contribute to the understanding of the members of the genus Flavobacterium diversity exploration and will provide a framework for the way for future ecological characterizations.

Profiling novel lateral gene transfer events in the human microbiome

Lateral gene transfer (LGT) is an important mechanism for genome diversification in microbial populations, including the human microbiome. While prior work has surveyed LGT events in human-associated microbial isolate genomes, the scope and dynamics of novel LGT events arising in personal microbiomes are not well understood, as there are no widely adopted computational methods to detect, quantify, and characterize LGT from complex microbial communities. We addressed this by developing, benchmarking, and experimentally validating a computational method (WAAFLE) to profile novel LGT events from assembled metagenomes. Applying WAAFLE to >2K human metagenomes from diverse body sites, we identified >100K putative high-confidence but previously uncharacterized LGT events (~2 per assembled microbial genome-equivalent). These events were enriched for mobile elements (as expected), as well as restriction-modification and transport functions typically associated with the destruction of foreign DNA. LGT frequency was quantifiably influenced by biogeography, the phylogenetic similarity of the involved taxa, and the ecological abundance of the donor taxon. These forces manifest as LGT networks in which hub species abundant in a community type donate unequally with their close phylogenetic neighbors. Our findings suggest that LGT may be a more ubiquitous process in the human microbiome than previously described. The open-source WAAFLE implementation, documentation, and data from this work are available at http://huttenhower.sph.harvard.edu/waafle.

Risk assessment and dissemination mechanism of antibiotic resistance genes in compost

In recent years, the excessive of antibiotics in livestock and poultry husbandry, stemming from extensive industry experience, has resulted in the accumulation of residual antibiotics and antibiotic resistance genes (ARGs) in livestock manure. Composting, as a crucial approach for the utilization of manure resources, has the potential to reduce the levels of antibiotics and ARGs in manure, although complete elimination is challenging. Previous studies have primarily focused on the diversity and abundance of ARGs in compost or have solely examined the correlation between ARGs and their carriers, potentially leading to a misjudgment of the actual risk associated with ARGs in compost. To address this gap, this study investigated the transfer potential of ARGs in compost and their co-occurrence with opportunistic pathogenic bacteria by extensively analyzing metagenomic sequencing data of compost worldwide. The results demonstrated that the potential risk of ARGs in compost was significantly lower than in manure, suggesting that composting effectively reduces the risk of ARGs. Further analysis showed that the microbes shifted their life history strategy in manure and compost due to antibiotic pressure and formed metabolic interactions dominated by antibiotic-resistant microbes, increasing ARG dissemination frequency. Therefore, husbandry practice without antibiotic addition was recommended to control ARG evolution, dissemination, and abatement both at the source and throughout processing.

Metagenomic analysis reveals the short-term influences on conjugation of blaNDM-1 and microbiome in hospital wastewater by silver nanoparticles at environmental-related concentration

Hospital wastewater contains large amounts of antibiotic-resistant bacteria and serves as an important reservoir for horizontal gene transfer (HGT). However, the response of the microbiome in hospital wastewater to silver remains unclear. In this study, the short-term impacts of silver on the microbiome in hospital wastewater were investigated by metagenome next-generation sequencing. The influence of silver on the conjugation of plasmid carrying bla_NDM-1 was further examined. Our results showed that in hospital wastewater, high abundances of antibiotic resistance genes (ARGs) were detected. The distribution tendencies of certain ARG types on chromosomes or plasmids were different. Clinically important ARGs were identified in phage-like contigs, indicating potential transmission via transduction. Pseudomonadales, Enterobacterales, and Bacteroidales were the major ARG hosts. Mobile genetic elements were mainly detected in plasmids and associated with various types of ARGs. The binning approach identified 29 bins that were assigned to three phyla. Various ARGs and virulence factors were identified in 14 and 11 bins, respectively. MetaCHIP identified 49 HGT events. The transferred genes were annotated as ARGs, mobile genetic elements, and functional genes, and they mainly originated from donors belonging to Bacteroides and Pseudomonadales. In addition, 20 nm AgNPs reduced microbial diversity and enhanced the relative abundance of Acinetobacter. The changes induced by 20 nm AgNPs included increases in the abundances of ARGs and genes involved lipid metabolism pathway. Conjugation experiments showed that Ag⁺ and 20 nm AgNPs caused 2.38-, 3.31-, 4.72-, and 4.57-fold and 1.46-, 1.61-, 3.86-, and 2.16-fold increases in conjugation frequencies of plasmid with bla_NDM-1 at 0.1, 1, 10, and 100 μg/L, respectively. Our findings provide insight into the response of the microbiome in hospital wastewater to silver, emphasize the adaptation capability of Acinetobacter inhabiting hospitals against adverse environments, and highlight the promotion of silver for antibiotic resistance.

Antimony efflux underpins phosphorus cycling and resistance of phosphate-solubilizing bacteria in mining soils

Microorganisms play crucial roles in phosphorus (P) turnover and P bioavailability increases in heavy metal-contaminated soils. However, microbially driven P-cycling processes and mechanisms of their resistance to heavy metal contaminants remain poorly understood. Here, we examined the possible survival strategies of P-cycling microorganisms in horizontal and vertical soil samples from the world's largest antimony (Sb) mining site, which is located in Xikuangshan, China. We found that total soil Sb and pH were the primary factors affecting bacterial community diversity, structure and P-cycling traits. Bacteria with the gcd gene, encoding an enzyme responsible for gluconic acid production, largely correlated with inorganic phosphate (Pi) solubilization and significantly enhanced soil P bioavailability. Among the 106 nearly complete bacterial metagenome-assembled genomes (MAGs) recovered, 60.4% carried the gcd gene. Pi transportation systems encoded by pit or pstSCAB were widely present in gcd-harboring bacteria, and 43.8% of the gcd-harboring bacteria also carried the acr3 gene encoding an Sb efflux pump. Phylogenetic and potential horizontal gene transfer (HGT) analyses of acr3 indicated that Sb efflux could be a dominant resistance mechanism, and two gcd-harboring MAGs appeared to acquire acr3 through HGT. The results indicated that Sb efflux could enhance P cycling and heavy metal resistance in Pi-solubilizing bacteria in mining soils. This study provides novel strategies for managing and remediating heavy metal-contaminated ecosystems.

Coal-source acid mine drainage reduced the soil multidrug-dominated antibiotic resistome but increased the heavy metal(loid) resistome and energy production-related metabolism

A recent global scale study found that mining-impacted environments have multi-antibiotic resistance gene (ARG)-dominated resistomes with an abundance similar to urban sewage but much higher than freshwater sediment. These findings raised concern that mining may increase the risk of ARG environmental proliferation. The current study assessed how typical multimetal(loid)-enriched coal-source acid mine drainage (AMD) contamination affects soil resistomes by comparing with background soils unaffected by AMD. Both contaminated and background soils have multidrug-dominated antibiotic resistomes attributed to the acidic environment. AMD-contaminated soils had a lower relative abundance of ARGs (47.45 ± 23.34 ×/Gb) than background soils (85.47 ± 19.71 ×/Gb) but held high-level heavy metal(loid) resistance genes (MRGs, 133.29 ± 29.36 ×/Gb) and transposase- and insertion sequence-dominated mobile genetic elements (MGEs, 188.51 ± 21.81 ×/Gb), which was 56.26 % and 412.12 % higher than background soils, respectively. Procrustes analysis showed that the microbial community and MGEs exerted more influence on driving heavy metal(loid) resistome variation than antibiotic resistome. The microbial community increased energy production-related metabolism to fulfill the increasing energy needs required by acid and heavy metal(loid) resistance. Horizontal gene transfer (HGT) events primarily exchanged energy- and information-related genes to adapt to the harsh AMD environment. These findings provide new insight into the risk of ARG proliferation in mining environments.

Assembly and analytical validation of a metagenomic reference catalog of human gut microbiota based on co-barcoding sequencing

Human gut microbiota is associated with human health and disease, and is known to have the second-largest genome in the human body. The microbiota genome is important for their functions and metabolites; however, accurate genomic access to the microbiota of the human gut is hindered due to the difficulty of cultivating and the shortcomings of sequencing technology. Therefore, we applied the stLFR library construction method to assemble the microbiota genomes and demonstrated that assembly property outperformed standard metagenome sequencing. Using the assembled genomes as references, SNP, INDEL, and HGT gene analyses were performed. The results demonstrated significant differences in the number of SNPs and INDELs among different individuals. The individual displayed a unique species variation spectrum, and the similarity of strains within individuals decreased over time. In addition, the coverage depth analysis of the stLFR method shows that a sequencing depth of 60X is sufficient for SNP calling. HGT analysis revealed that the genes involved in replication, recombination and repair, mobilome prophages, and transposons were the most transferred genes among different bacterial species in individuals. A preliminary framework for human gut microbiome studies was established using the stLFR library construction method.

Metagenome-based metabolic modelling predicts unique microbial interactions in deep-sea hydrothermal plume microbiomes

Deep-sea hydrothermal vents are abundant on the ocean floor and play important roles in ocean biogeochemistry. In vent ecosystems such as hydrothermal plumes, microorganisms rely on reduced chemicals and gases in hydrothermal fluids to fuel primary production and form diverse and complex microbial communities. However, microbial interactions that drive these complex microbiomes remain poorly understood. Here, we use microbiomes from the Guaymas Basin hydrothermal system in the Pacific Ocean to shed more light on the key species in these communities and their interactions. We built metabolic models from metagenomically assembled genomes (MAGs) and infer possible metabolic exchanges and horizontal gene transfer (HGT) events within the community. We highlight possible archaea-archaea and archaea-bacteria interactions and their contributions to the robustness of the community. Cellobiose, D-Mannose 1-phosphate, O2, CO2, and H2S were among the most exchanged metabolites. These interactions enhanced the metabolic capabilities of the community by exchange of metabolites that cannot be produced by any other community member. Archaea from the DPANN group stood out as key microbes, benefiting significantly as acceptors in the community. Overall, our study provides key insights into the microbial interactions that drive community structure and organisation in complex hydrothermal plume microbiomes.

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.