Functions predict horizontal gene transfer and the emergence of antibiotic resistance

Ecological and evolutionary mechanisms driving within-patient emergence of antimicrobial resistance

The ecological and evolutionary mechanisms of antimicrobial resistance (AMR) emergence within patients and how these vary across bacterial infections are poorly understood. Increasingly widespread use of pathogen genome sequencing in the clinic enables a deeper understanding of these processes. In this Review, we explore the clinical evidence to support four major mechanisms of within-patient AMR emergence in bacteria: spontaneous resistance mutations; in situ horizontal gene transfer of resistance genes; selection of pre-existing resistance; and immigration of resistant lineages. Within-patient AMR emergence occurs across a wide range of host niches and bacterial species, but the importance of each mechanism varies between bacterial species and infection sites within the body. We identify potential drivers of such differences and discuss how ecological and evolutionary analysis could be embedded within clinical trials of antimicrobials, which are powerful but underused tools for understanding why these mechanisms vary between pathogens, infections and individuals. Ultimately, improving understanding of how host niche, bacterial species and antibiotic mode of action combine to govern the ecological and evolutionary mechanism of AMR emergence in patients will enable more predictive and personalized diagnosis and antimicrobial therapies.

How micro-/nano-plastics influence the horizontal transfer of antibiotic resistance genes - A review

Plastic debris such as microplastics (MPs) and nanoplastics (NPTs), along with antibiotic resistance genes (ARGs), are pervasive in the environment and are recognized as significant global health and ecological concerns. Micro-/nano-plastics (MNPs) have been demonstrated to favor the spread of ARGs by enhancing the frequency of horizontal gene transfer (HGT) through various pathways. This paper comprehensively and systematically reviews the current study with focus on the influence of plastics on the HGT of ARGs. The critical role of MNPs in the HGT of ARGs has been well illustrated in sewage sludge, livestock farms, constructed wetlands and landfill leachate. A summary of the performed HGT assay and the underlying mechanism of plastic-mediated transfer of ARGs is presented in the paper. MNPs could facilitate or inhibit HGT of ARGs, and their effects depend on the type, size, and concentration. This review provides a comprehensive insight into the effects of MNPs on the HGT of ARGs, and offers suggestions for further study. Further research should attempt to develop a standard HGT assay and focus on investigating the impact of different plastics, including the oligomers they released, under real environmental conditions on the HGT of ARGs.

Horizontal plasmid transfer promotes antibiotic resistance in selected bacteria in Chinese frog farms

The emergence and dissemination of antibiotic resistance genes (ARGs) in the ecosystem are global public health concerns. One Health emphasizes the interconnectivity between different habitats and seeks to optimize animal, human, and environmental health. However, information on the dissemination of antibiotic resistance genes (ARGs) within complex microbiomes in natural habitats is scarce. We investigated the prevalence of antibiotic resistant bacteria (ARB) and the spread of ARGs in intensive bullfrog (Rana catesbeiana) farms in the Shantou area of China. Antibiotic susceptibilities of 361 strains, combined with microbiome analyses, revealed Escherichia coli, Edwardsiella tarda, Citrobacter and Klebsiella sp. as prevalent multidrug resistant bacteria on these farms. Whole genome sequencing of 95 ARB identified 250 large plasmids that harbored a wide range of ARGs. Plasmid sequences and sediment metagenomes revealed an abundance of tetA, sul1, and aph(3″)-Ib ARGs. Notably, antibiotic resistance (against 15 antibiotics) highly correlated with plasmid-borne rather than chromosome-borne ARGs. Based on sequence similarities, most plasmids (62%) fell into 32 distinct groups, indicating a potential for horizontal plasmid transfer (HPT) within the frog farm microbiome. HPT was confirmed in inter- and intra-species conjugation experiments. Furthermore, identical mobile ARGs, flanked by mobile genetic elements (MGEs), were found in different locations on the same plasmid, or on different plasmids residing in the same or different hosts. Our results suggest a synergy between MGEs and HPT to facilitate ARGs dissemination in frog farms. Mining public databases retrieved similar plasmids from different bacterial species found in other environmental niches globally. Our findings underscore the importance of HPT in mediating the spread of ARGs in frog farms and other microbiomes of the ecosystem.

Augmenting bacterial similarity measures using a graph-based genome representation

Relationships between bacterial taxa are traditionally defined using 16S rRNA nucleotide similarity or average nucleotide identity. Improvements in sequencing technology provide additional pairwise information on genome sequences, which may provide valuable information on genomic relationships. Mapping orthologous gene locations between genome pairs, known as synteny, is typically implemented in the discovery of new species and has not been systematically applied to bacterial genomes. Using a data set of 378 bacterial genomes, we developed and tested a new measure of synteny similarity between a pair of genomes, which was scaled onto 16S rRNA distance using covariance matrices. Based on the input gene functions used (i.e., core, antibiotic resistance, and virulence), we observed varying topological arrangements of bacterial relationship networks by applying (i) complete linkage hierarchical clustering and (ii) K-nearest neighbor graph structures to synteny-scaled 16S data. Our metric improved clustering quality comparatively to state-of-the-art average nucleotide identity metrics while preserving clustering assignments for the highest similarity relationships. Our findings indicate that syntenic relationships provide more granular and interpretable relationships for within-genera taxa compared to pairwise similarity measures, particularly in functional contexts.

IMPORTANCE

Given the prevalence and necessity of the 16S rRNA measure in bacterial identification and analysis, this additional analysis adds a functional and synteny-based layer to the identification of relatives and clustering of bacteria genomes. It is also of computational interest to model the bacterial genome as a graph structure, which presents new avenues of genomic analysis for bacteria and their closely related strains and species.

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.

A pervasive large conjugative plasmid mediates multispecies biofilm formation in the intestinal microbiota increasing resilience to perturbations

Although horizontal gene transfer is pervasive in the intestinal microbiota, we understand only superficially the roles of most exchanged genes and how the mobile repertoire affects community dynamics. Similarly, little is known about the mechanisms underlying the ability of a community to recover after a perturbation. Here, we identified and functionally characterized a large conjugative plasmid that is one of the most frequently transferred elements among Bacteroidales species and is ubiquitous in diverse human populations. This plasmid encodes both an extracellular polysaccharide and fimbriae, which promote the formation of multispecies biofilms in the mammalian gut. We use a hybridization-based approach to visualize biofilms in clarified whole colon tissue with unprecedented 3D spatial resolution. These biofilms increase bacterial survival to common stressors encountered in the gut, increasing strain resiliency, and providing a rationale for the plasmid's recent spread and high worldwide prevalence.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Microbiome single cell atlases generated with a commercial instrument

Single cell sequencing is useful for resolving complex systems into their composite cell types and computationally mining them for unique features that are masked in pooled sequencing. However, while commercial instruments have made single cell analysis widespread for mammalian cells, analogous tools for microbes are limited. Here, we present EASi-seq (Easily Accessible Single microbe sequencing). By adapting the single cell workflow of the commercial Mission Bio Tapestri instrument, this method allows for efficient sequencing of individual microbes' genomes. EASi-seq allows thousands of microbes to be sequenced per run and, as we show, can generate detailed atlases of human and environmental microbiomes. The ability to capture large shotgun genome datasets from thousands of single microbes provides new opportunities in discovering and analyzing species subpopulations. To facilitate this, we develop a companion bioinformatic pipeline that clusters microbes by similarity, improving whole genome assembly, strain identification, taxonomic classification, and gene annotation. In addition, we demonstrate integration of metagenomic contigs with the EASi-seq datasets to reduce capture bias and increase coverage. Overall, EASi-seq enables high quality single cell genomic data for microbiome samples using an accessible workflow that can be run on a commercially available platform.

The integrative and conjugative element ICECiPOL15 mediates horizontal transfer of β-lactam resistance gene in Chryseobacterium indoltheticum POL15

OBJECTIVES

The dissemination of antibiotic resistance genes (ARGs) from the environment, including agricultural sources, is of increasing concern. In this study, we examined the antibiotic resistance profile and genomic sequence of a strain of Chryseobacterium indoltheticum obtained from an agricultural location.

METHODS

The multidrug-resistant bacterial strain POL15 was isolated from the wastewater of a livestock farm in China. Whole-genome sequencing was performed followed by bioinformatics analyses to identify integrative and conjugative elements (ICEs) and ARGs. Mating assays were performed to analyse ICE transferability.

RESULTS

Whole-genome sequencing and annotation showed that the genome of POL15 encodes ARGs. Additionally, an ICE named ICECiPOL15, which carries a class C β-lactamase-encoding gene blaAQU, was identified in the POL15 genome. Genes encoding an integrase, an excisionase, a relaxase, a type IV coupling protein and conjugative transposon proteins involved in a type IV secretion system were also identified in ICECiPOL15. Sequence alignment revealed that ICECiPOL15 might have evolved from other Chryseobacterium species. The horizontal transferability of ICECiPOL15 was demonstrated by mating experiments between C. indoltheticum POL15 and Escherichia coli DL21.

CONCLUSIONS

This study represents the first characterization of a mobilizable antibiotic resistance ICE in a species of C. indoltheticum and provides evidence that C. indoltheticum strains could be important reservoirs and vehicles for ARGs on livestock farms.

Dissemination of antibiotic resistance genes is regulated by iron oxides: Insight into the influence on bacterial transformation

The transportation of antibiotic resistance genes (ARGs) in manure-soil-plant continuums poses risks to human health. Horizontal gene transfer, particularly for bacterial transformation, is an important way for ARG dissemination. As crucial components in soils, iron oxides impacted the fates of various abiotic and biotic contaminants due to their active properties. However, whether they can influence the transformation of ARGs is unknown, which waits to be figured out to boost the assessment and control of ARG spread risks. In this study, we have investigated the effects of goethite, hematite, and magnetite (0-250 mg/L, with sizes < 100 nm and > 100 nm) on the transfer of ampicillin resistance genes to Escherichia coli cells. At lower iron oxide concentrations, the transformation of ARGs was first facilitated (transformation frequency reached up to 3.38-fold higher), but the facilitating effects gradually weakened and eventually disappeared as concentrations further increased. Particle size and iron oxide type were not the universal determinants controlling the transformation. At lower concentrations, iron oxides interacted with proteins and phospholipids in E. coli envelope structures, and induced the overgeneration of intracellular reactive oxygen species. Consequently, they led to pore formation and permeability enhancement on the cell membrane, thus promoting the transformation. The facilitation was also associated with the carrier-like effect of iron oxides for antibiotic resistance plasmids. At higher concentrations, the weakened facilitations were attributed to the aggregation of iron oxides. In this study, we highlight the crucial roles of the concentrations (contents) of iron oxides on the dissemination of ARGs in soils; this study may serve as a reference for ARG pollution control in future agricultural production.

Zinc oxide nanoparticles enhance plasmid transfer among growth-promoting endophytes in Arabidopsis thaliana

Nanoparticles (NPs) hold great promise for agricultural applications, yet their potential impact on exogenous gene transfer within plant remains poorly understood. In this study, we utilized the non-conjugative plasmid pCAMBIA1300, harboring the bialaphos resistance (bar) gene expressed in plant and the kanamycin resistance (kanR) gene as selectable marker in bacteria. Our results revealed a significant increase in the transfer of plasmid (via carrier Escherichia coli DH5α), both intra- and inter-species within the endophyte, when Arabidopsis thaliana was exposed to environmentally relevant level of zinc oxide (ZnO) NPs at a concentration of 0.7 μg/mL throughout its lifespan. Intriguingly, the plasmid exhibited selective transfer to growth-promoting endophytes, such as Enterobacter, Serratia, and Achromobacter, with the presence of ZnO NPs expanding the pool of potential recipients. This result is due to the facilitation of an endophytic and mutualistic lifestyle of invasive E. coli DH5α and the enrichment of beneficial bacteria aided by ZnO NPs. The plant's descendant generations did not express the bar gene, and the endophytes carrying the exogenous plasmid did not transmit it to sub sequent generation. This research provides crucial insights for assessing the potential risks associated with gene contamination and ensuring the safe and sustainable use of NPs in agriculture.

Antibiotic resistance genes associated with size-segregated bioaerosols from wastewater treatment plants: A review

The antibiotic-resistant pollution in size-segregated bioaerosols from wastewater treatment plants (WWTPs) is of increasing concern due to its public health risks, but an elaborate review is still lacking. This work overviewed the profile, mobility, pathogenic hosts, source, and risks of antibiotic resistance genes (ARGs) in size-segregated bioaerosols from WWTPs. The dominant ARG type in size-segregated bioaerosols from WWTPs was multidrug resistance genes. Treatment units that equipped with mechanical facilities and aeration devices, such as grilles, grit chambers, biochemical reaction tanks, and sludge treatment units, were the primary sources of bioaerosol antibiotic resistome in WWTPs. Higher enrichment of antibiotic resistome in particulate matter with an aerodynamic diameter of <2.5 μm, was found along the upwind-downwind-WWTPs gradient. Only a small portion of ARGs in inhalable bioaerosols from WWTPs were flanked by mobile genetic elements. The pathogens with multiple drug resistance had been found in size-segregated bioaerosols from WWTPs. Different ARGs or antibiotic resistant bacteria have different aerosolization potential associated with bioaerosols from various treatment processes. The validation of pathogenic antibiotic resistance bacteria, deeper investigation of ARG mobility, emission mechanism of antibiotic resistome, and development of treatment technologies, should be systematically considered in future.

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.

Inter-plasmid transfer of antibiotic resistance genes accelerates antibiotic resistance in bacterial pathogens

Antimicrobial resistance is a major threat for public health. Plasmids play a critical role in the spread of antimicrobial resistance via horizontal gene transfer between bacterial species. However, it remains unclear how plasmids originally recruit and assemble various antibiotic resistance genes (ARGs). Here, we track ARG recruitment and assembly in clinically relevant plasmids by combining a systematic analysis of 2420 complete plasmid genomes and experimental validation. Results showed that ARG transfer across plasmids is prevalent, and 87% ARGs were observed to potentially transfer among various plasmids among 8229 plasmid-borne ARGs. Interestingly, recruitment and assembly of ARGs occur mostly among compatible plasmids within the same bacterial cell, with over 88% of ARG transfers occurring between compatible plasmids. Integron and insertion sequences drive the ongoing ARG acquisition by plasmids, especially in which IS26 facilitates 63.1% of ARG transfer events among plasmids. In vitro experiment validated the important role of IS26 involved in transferring gentamicin resistance gene aacC1 between compatible plasmids. Network analysis showed four beta-lactam genes (blaTEM-1, blaNDM-4, blaKPC-2, and blaSHV-1) shuffling among 1029 plasmids and 45 clinical pathogens, suggesting that clinically alarming ARGs transferred accelerate the propagation of antibiotic resistance in clinical pathogens. ARGs in plasmids are also able to transmit across clinical and environmental boundaries, in terms of the high-sequence similarities of plasmid-borne ARGs between clinical and environmental plasmids. This study demonstrated that inter-plasmid ARG transfer is a universal mechanism for plasmid to recruit various ARGs, thus advancing our understanding of the emergence of multidrug-resistant plasmids.

A highly selective and sensitive electrochemical sensor for tetracycline resistant genes detection based on the non-covalent interaction of graphene oxide and nucleobase

Antibiotic resistance genes (ARGs) are causing worldwide environmental problems, however, the traditional analytical methods and test equipment for them are time-consuming and expensive. The electrochemical sensor using the non-covalent bond between graphene oxide (GO) and single-stranded tet (ss-tet) was established for specific tetracycline resistance genes (tet, composed of ss-tet and complementary ss-tet (ss-tet') in water) detection, which preparation time was only 35 min and far less than most reported sensors based on covalent bond. As the result of the detection for tet, the developed sensor not only had the low detection limit of 50.0 pM (8.1 × 102 copies·mL-1), the short detection time within 42 min, but also had satisfactory stability, excellent reproducibility, and highly selectivity (RSD < 4.43 %). Besides, it also had acceptable accuracy comparing to the real-time quantitative polymerase chain reaction (RT-qPCR) and PCR array in tet detection. Noticeably, it also had been successfully applied to tetA detection in different water samples. In brief, the prepared non-covalent bond sensor is simple, rapid, and suitable for highly selective and sensitive detection of the ARGs in actual water.

Clinically relevant antibiotic resistance genes are linked to a limited set of taxa within gut microbiome worldwide

The acquisition of antimicrobial resistance (AR) genes has rendered important pathogens nearly or fully unresponsive to antibiotics. It has been suggested that pathogens acquire AR traits from the gut microbiota, which collectively serve as a global reservoir for AR genes conferring resistance to all classes of antibiotics. However, only a subset of AR genes confers resistance to clinically relevant antibiotics, and, although these AR gene profiles are well-characterized for common pathogens, less is known about their taxonomic associations and transfer potential within diverse members of the gut microbiota. We examined a collection of 14,850 human metagenomes and 1666 environmental metagenomes from 33 countries, in addition to nearly 600,000 isolate genomes, to gain insight into the global prevalence and taxonomic range of clinically relevant AR genes. We find that several of the most concerning AR genes, such as those encoding the cephalosporinase CTX-M and carbapenemases KPC, IMP, NDM, and VIM, remain taxonomically restricted to Proteobacteria. Even cfiA, the most common carbapenemase gene within the human gut microbiome, remains tightly restricted to Bacteroides, despite being found on a mobilizable plasmid. We confirmed these findings in gut microbiome samples from India, Honduras, Pakistan, and Vietnam, using a high-sensitivity single-cell fusion PCR approach. Focusing on a set of genes encoding carbapenemases and cephalosporinases, thus far restricted to Bacteroides species, we find that few mutations are required for efficacy in a different phylum, raising the question of why these genes have not spread more widely. Overall, these data suggest that globally prevalent, clinically relevant AR genes have not yet established themselves across diverse commensal gut microbiota.

Prevalence and transmission risk of colistin and multidrug resistance in long-distance coastal aquaculture

Due to the wide use of antibiotics, intensive aquaculture farms have been recognized as a significant reservoir of antibiotic resistomes. Although the prevalence of colistin resistance genes and multidrug-resistant bacteria (MDRB) has been documented, empirical evidence for the transmission of colistin and multidrug resistance between bacterial communities in aquaculture farms through horizontal gene transfer (HGT) is lacking. Here, we report the prevalence and transmission risk of colistin and multidrug resistance in 27 aquaculture water samples from 9 aquaculture zones from over 5000 km of subtropical coastlines in southern China. The colistin resistance gene mcr-1, mobile genetic element (MGE) intl1 and 13 typical antibiotic resistance genes (ARGs) were prevalent in all the aquaculture water samples. Most types of antibiotic (especially colistin) resistance are transmissible in bacterial communities based on evidence from laboratory conjugation and transformation experiments. Diverse MDRB were detected in most of the aquaculture water samples, and a strain with high-level colistin resistance, named Ralstonia pickettii MCR, was isolated. The risk of horizontal transfer of the colistin resistance of R. pickettii MCR through conjugation and transformation was low, but the colistin resistance could be steadily transmitted to offspring through vertical transfer. The findings have important implications for the future regulation of antibiotic use in aquaculture farms globally to address the growing threat posed by antibiotic resistance to human health.

Metagenomic and Machine Learning Meta-Analyses Characterize Airborne Resistome Features and Their Hosts in China Megacities

Urban ambient air contains a cocktail of antibiotic resistance genes (ARGs) emitted from various anthropogenic sites. However, what is largely unknown is whether the airborne ARGs exhibit site-specificity or their pathogenic hosts persistently exist in the air. Here, by retrieving 1.2 Tb metagenomic sequences (n = 136), we examined the airborne ARGs from hospitals, municipal wastewater treatment plants (WWTPs) and landfills, public transit centers, and urban sites located in seven of China's megacities. As validated by the multiple machine learning-based classification and optimization, ARGs' site-specificity was found to be the most apparent in hospital air, with featured resistances to clinical-used rifamycin and (glyco)peptides, whereas the more environmentally prevalent ARGs (e.g., resistance to sulfonamide and tetracycline) were identified being more specific to the nonclinical ambient air settings. Nearly all metagenome-assembled genomes (MAGs) that possessed the site-featured resistances were identified as pathogenic taxa, which occupied the upper-representative niches in all the neutrally distributed airborne microbial community (P < 0.01, m = 0.22-0.50, R2 = 0.41-0.86). These niche-favored putative resistant pathogens highlighted the enduring antibiotic resistance hazards in the studied urban air. These findings are critical, albeit the least appreciated until our study, to gauge the airborne dimension of resistomes' features and fates in urban atmospheric environments.

Antibiotics and their associations with antibiotic resistance genes and microbial communities in estuarine and coastal sediment of Quanzhou Bay, Southeast China

The antibiotic concentrations spanned from 11.2 to 173.8 ng/g, with quinolones and tetracyclines being observed to be prevalent. The amount of microbial biomass as determined by Phospholipid fatty acid (PLFA) ranged from 2.92 to 10.99 mg kg-1, with G- bacteria dominating. A total of 254 distinct ARGs and 10 MEGs were identified, with multidrug ARGs having the highest relative abundance (1.18 × 10-2 to 3.00 × 10-1 copies/16S rRNA gene copies), while vancomycin and sulfonamide resistance genes were the least abundant. Results from canonical-correlation analyses combined with redundancy analysis indicated that macrolides were significantly related to the shifts of microbial community structure in sediments, particularly in G+ bacteria that were more sensitive to antibiotic residues. It was observed that sulfonamide ARGs had a greater correlation with residual antibiotics than other ARGs. This study provided a field evidence that multiple residual antibiotics from coastal sites could cause fundamental shifts in microbial community and their associated ARGs.

Antibiotic resistome associated with inhalable bioaerosols from wastewater to atmosphere: Mobility, bacterial hosts, source contributions and resistome risk

Antibiotic resistome can be carried by the bioaerosols and propagate from wastewater treatment plants (WWTPs) to the atmosphere, but questions remain regarding their mobility, bacterial hosts, source, and resistome risk. Here, fine particulate matter (PM2.5) was collected within and around a large WWTP and analyzed by the metagenomic assembly and binning. PM2.5 was discovered with increasing enrichment of total antibiotic resistance genes (ARGs), potentially mobile ARGs, and antibiotic-resistant bacteria (ARB) along the WWTP-downwind-upwind gradient. Some ARGs were found to be flanked by certain mobile genetic elements and generally mediated by plasmids in WWTP-PM2.5. Totally, 198 metagenome assembled genomes assigning to seven phyla were identified as the ARB, and a contig-based analysis indicated that 32 pathogens were revealed harboring at least two ARGs. Despite disparate aerosolization potentials of ARGs or ARB at different WWTP units, high resistome risks were found, along with the dominant contribution of wastewater for airborne ARGs (44.79-62.82%) and ARB (35.03-40.10%). Among the detected WWTP matrices, the sludge dewatering room was characterized by the highest resistome risk associated with PM2.5. This study underscores the dispersion of ARGs and ARB from WWTPs to the atmosphere and provides a reference for managing risks of antibiotic resistance.

Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota

The microbiota inhabiting soil plays a significant role in essential life-supporting element cycles. Here, we investigated the occurrence of horizontal gene transfer (HGT) and established the HGT network of carbon metabolic genes in 764 soil-borne microbiota genomes. Our study sheds light on the crucial role of HGT components in microbiological diversification that could have far-reaching implications in understanding how these microbial communities adapt to changing environments, ultimately impacting agricultural practices. In the overall HGT network of carbon metabolic genes in soil-borne microbiota, a total of 6,770 nodes and 3,812 edges are present. Among these nodes, phyla Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes are predominant. Regarding specific classes, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacteroidia, Actinomycetia, Betaproteobacteria, and Clostridia are dominant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments of glycosyltransferase (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporter (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%) and butanoate metabolism (3.0%) are primarily identified. Glycosyltransferase involved in cell wall biosynthesis, glycosylation, and primary/secondary metabolism (with 363 HGT entries), ranks first overwhelmingly in the list of most frequently identified carbon metabolic HGT enzymes, followed by pimeloyl-ACP methyl ester carboxylesterase, alcohol dehydrogenase, and 3-oxoacyl-ACP reductase. Such HGT events mainly occur in the peripheral functions of the carbon metabolic pathway instead of the core section. The inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized, as well as their involvement in the metabolism and regulation processes of carbon organic, suggest a pervasive and substantial effect of HGT on the evolution of microbes.

Whole-cell modeling of E. coli colonies enables quantification of single-cell heterogeneity in antibiotic responses

Antibiotic resistance poses mounting risks to human health, as current antibiotics are losing efficacy against increasingly resistant pathogenic bacteria. Of particular concern is the emergence of multidrug-resistant strains, which has been rapid among Gram-negative bacteria such as Escherichia coli. A large body of work has established that antibiotic resistance mechanisms depend on phenotypic heterogeneity, which may be mediated by stochastic expression of antibiotic resistance genes. The link between such molecular-level expression and the population levels that result is complex and multi-scale. Therefore, to better understand antibiotic resistance, what is needed are new mechanistic models that reflect single-cell phenotypic dynamics together with population-level heterogeneity, as an integrated whole. In this work, we sought to bridge single-cell and population-scale modeling by building upon our previous experience in "whole-cell" modeling, an approach which integrates mathematical and mechanistic descriptions of biological processes to recapitulate the experimentally observed behaviors of entire cells. To extend whole-cell modeling to the "whole-colony" scale, we embedded multiple instances of a whole-cell E. coli model within a model of a dynamic spatial environment, allowing us to run large, parallelized simulations on the cloud that contained all the molecular detail of the previous whole-cell model and many interactive effects of a colony growing in a shared environment. The resulting simulations were used to explore the response of E. coli to two antibiotics with different mechanisms of action, tetracycline and ampicillin, enabling us to identify sub-generationally-expressed genes, such as the beta-lactamase ampC, which contributed greatly to dramatic cellular differences in steady-state periplasmic ampicillin and was a significant factor in determining cell survival.

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

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

Addressing antibiotic resistance: computational answers to a biological problem?

The increasing prevalence of infections caused by antibiotic-resistant bacteria is a global healthcare crisis. Understanding the spread of resistance is predicated on the surveillance of antibiotic resistance genes within an environment. Bioinformatics and artificial intelligence (AI) methods applied to metagenomic sequencing data offer the capacity to detect known and infer yet-unknown resistance mechanisms, and predict future outbreaks of antibiotic-resistant infections. Machine learning methods, in particular, could revive the waning antibiotic discovery pipeline by helping to predict the molecular structure and function of antibiotic resistance compounds, and optimising their interactions with target proteins. Consequently, AI has the capacity to play a central role in guiding antibiotic stewardship and future clinical decision-making around antibiotic resistance.

Uncovering Bacterial Hosts of Class 1 Integrons in an Urban Coastal Aquatic Environment with a Single-Cell Fusion-Polymerase Chain Reaction Technology

Horizontal gene transfer (HGT) is a key driver of bacterial evolution via transmission of genetic materials across taxa. Class 1 integrons are genetic elements that correlate strongly with anthropogenic pollution and contribute to the spread of antimicrobial resistance (AMR) genes via HGT. Despite their significance to human health, there is a shortage of robust, culture-free surveillance technologies for identifying uncultivated environmental taxa that harbor class 1 integrons. We developed a modified version of epicPCR (emulsion, paired isolation, and concatenation polymerase chain reaction (PCR)) that links class 1 integrons amplified from single bacterial cells to taxonomic markers from the same cells in emulsified aqueous droplets. Using this single-cell genomic approach and Nanopore sequencing, we successfully assigned class 1 integron gene cassette arrays containing mostly AMR genes to their hosts in coastal water samples that were affected by pollution. Our work presents the first application of epicPCR for targeting variable, multigene loci of interest. We also identified the Rhizobacter genus as novel hosts of class 1 integrons. These findings establish epicPCR as a powerful tool for linking taxa to class 1 integrons in environmental bacterial communities and offer the potential to direct mitigation efforts toward hotspots of class 1 integron-mediated dissemination of AMR.

The emergence of novel macrolide resistance island in Macrococcus caseolyticus and Staphylococcus aureus of food origin

Food-derived Staphylococcaceae species with severe antimicrobial resistance, especially Staphylococcus aureus, is a major threat to public health. Macrococcus caseolyticus (M. caseolyticus) is a member of the Staphylococcaceae family which plays a vital role in fermented products and disease causation in animals. In our previous study, several Staphylococcus aureus antibiotic-resistant island msr (SaRI_msr) were found in multidrug-resistant S. aureus. In this study, novel SaRI_msr, SaRI_msr-III emerged from S. aureus. Another novel SaRI_msr-like further emerged in M. caseolyticus from food. These isolates' prevalence and genetic environment were investigated and characterized to understand the distribution and transmission of these novel SaRI_msr strains. All SaRI_msr-positive S. aureus isolates exhibited a multidrug resistance (MDR) phenotype, within which a series of antimicrobial resistance genes (ARGs) and virulence factor genes (VFs) were identified. In addition, three SaRI_msr types, SaRI_msr-I (15.1 kb), SaRI_msr-II (16-17 kb), and SaRI_msr-III (18 kb) carrying mef(D)-msr(F), were identified in these isolates' chromosomes. SaRI_msr-(I-III) contains a site-specific integrase gene int and operon mef(D)-msr(F). SaRI_msr-III has an additional orf3-orf4-IS30 arrangement downstream of mef(D) and msr(F). Moreover, the SaRI_msr-like and macrolide-resistant transposon Tn6776 forming a novel mosaic structure coexisted in one M. caseolyticus isolate. Within this mosaic structure, the macrolide-resistant genes mef(D)-msr(F) were absent in SaRI_msr-like, whereas an operon, mef(F)-msr(G), was identified in Tn6776. The SaRI_msr-(I-III) and SaRI_msr-like structure were inserted into the rpsI gene encoding the 30S ribosomal protein S9 in the chromosome. Excision and cyclisation of SaRI_msr-III, SaRI_msr-like, operon mef(D)-msr(F), and orf3-orf4-IS30 arrangements were confirmed using two-step PCR. This study is the first to report MDR S. aureus harbouring novel SaRI_msr-III and M. caseolyticus containing novel mosaic structures isolated from retail foods. Similar SaRI_msr-type resistant islands' occurrence and propagation in Staphylococcaceae species require continuous monitoring and investigation.

Advances in linking single-cell bacterial stress response to population-level survival

Stress response mechanisms can allow bacteria to survive a myriad of challenges, including nutrient changes, antibiotic encounters, and antagonistic interactions with other microbes. Expression of these stress response pathways, in addition to other cell features such as growth rate and metabolic state, can be heterogeneous across cells and over time. Collectively, these single-cell-level phenotypes contribute to an overall population-level response to stress. These include diversifying actions, which can be used to enable bet-hedging, and coordinated actions, such as biofilm production, horizontal gene transfer, and cross-feeding. Here, we highlight recent results and emerging technologies focused on both single-cell and population-level responses to stressors, and we draw connections about the combined impact of these effects on survival of bacterial communities.

Graph representation learning in biomedicine and healthcare

Networks-or graphs-are universal descriptors of systems of interacting elements. In biomedicine and healthcare, they can represent, for example, molecular interactions, signalling pathways, disease co-morbidities or healthcare systems. In this Perspective, we posit that representation learning can realize principles of network medicine, discuss successes and current limitations of the use of representation learning on graphs in biomedicine and healthcare, and outline algorithmic strategies that leverage the topology of graphs to embed them into compact vectorial spaces. We argue that graph representation learning will keep pushing forward machine learning for biomedicine and healthcare applications, including the identification of genetic variants underlying complex traits, the disentanglement of single-cell behaviours and their effects on health, the assistance of patients in diagnosis and treatment, and the development of safe and effective medicines.

Co-opted genes of algal origin protect C. elegans against cyanogenic toxins

Amygdalin is a cyanogenic glycoside enriched in the tissues of many edible plants, including seeds of stone fruits such as cherry (Prunus avium), peach (Prunus persica), and apple (Malus domestica). These plants biosynthesize amygdalin in defense against herbivore animals, as amygdalin generates poisonous cyanide upon plant tissue destruction.1,2,3,4 Poisonous to many animals, amygdalin-derived cyanide is detoxified by potent enzymes commonly found in bacteria and plants but not most animals.5 Here we show that the nematode C. elegans can detoxify amygdalin by a genetic pathway comprising cysl-1, egl-9, hif-1, and cysl-2. A screen of a natural product library for hypoxia-independent regulators of HIF-1 identifies amygdalin as a potent activator of cysl-2, a HIF-1 transcriptional target that encodes a cyanide detoxification enzyme in C. elegans. As a cysl-2 paralog similarly essential for amygdalin resistance, cysl-1 encodes a protein homologous to cysteine biosynthetic enzymes in bacteria and plants but functionally co-opted in C. elegans. We identify exclusively HIF-activating egl-9 mutations in a cysl-1 suppressor screen and show that cysl-1 confers amygdalin resistance by regulating HIF-1-dependent cysl-2 transcription to protect against amygdalin toxicity. Phylogenetic analysis indicates that cysl-1 and cysl-2 were likely acquired from green algae through horizontal gene transfer (HGT) and functionally co-opted in protection against amygdalin. Since acquisition, these two genes evolved division of labor in a cellular circuit to detect and detoxify cyanide. Thus, algae-to-nematode HGT and subsequent gene function co-option events may facilitate host survival and adaptation to adverse environmental stresses and biogenic toxins.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Informing plasmid compatibility with bacterial hosts using protein-protein interaction data

The compatibility of plasmids with new host cells is significant given their role in spreading antimicrobial resistance (AMR) and virulence factor genes. Evaluating this using in vitro screening is laborious and can be informed by computational analyses of plasmid-host compatibility through rates of protein-protein interactions (PPIs) between plasmid and host cell proteins. We identified large excesses of such PPIs in eight important plasmids, including pOXA-48, using most known bacteria (n = 4363). 23 species had high rates of interactions with four blaOXA-48-positive plasmids. We also identified 48 species with high interaction rates with plasmids common in Escherichia coli. We found a strong association between one plasmid and the fimbrial adhesin operon pil, which could enhance host cell adhesion in aqueous environments. An excess rate of PPIs could be a sign of host-plasmid compatibility, which is important for AMR control given that plasmids like pOXA-48 move between species with ease.

Ferric chloride further simplified the horizontal gene transfer network of antibiotic resistance genes in anaerobic digestion

The role of ferric chloride (FC) on the reduction of antibiotic resistance genes (ARGs) in anaerobic digestion (AD) system was investigated from the perspective of vertical (VGT) and horizontal gene transfer (HGT) network through the high-throughput qPCR (HT-qPCR). Although FC showed limited impacts on methane production in AD of swine manure, the tetracycline and MLSB resistance genes were specifically reduced at the end, where tetQ of antibiotic target protection and ermF of antibiotic target alteration contributed the most to the reduction. Both VGT and HGT network were divided into three modules, and the complexity of HGT network was largely reduced along with AD, where the HGT connection was reduced from 683 (Module III) to 172 (Module I), and FC addition could further reduce the relative abundance of ARG hosts in Module I. The contribution of VGT and HGT to the changes of ARGs in AD was further deciphered, and although the VGT reflected by the changes of microbial community contributed the most to the dynamics of ARGs (68.0 %), the HGT contribution could further be reduced by the FC addition. This study provided a new perspective on the fate of ARGs response to the FC addition in the AD system.

Link Between Antibiotic Persistence and Antibiotic Resistance in Bacterial Pathogens

Both, antibiotic persistence and antibiotic resistance characterize phenotypes of survival in which a bacterial cell becomes insensitive to one (or even) more antibiotic(s). However, the molecular basis for these two antibiotic-tolerant phenotypes is fundamentally different. Whereas antibiotic resistance is genetically determined and hence represents a rather stable phenotype, antibiotic persistence marks a transient physiological state triggered by various stress-inducing conditions that switches back to the original antibiotic sensitive state once the environmental situation improves. The molecular basics of antibiotic resistance are in principle well understood. This is not the case for antibiotic persistence. Under all culture conditions, there is a stochastically formed, subpopulation of persister cells in bacterial populations, the size of which depends on the culture conditions. The proportion of persisters in a bacterial population increases under different stress conditions, including treatment with bactericidal antibiotics (BCAs). Various models have been proposed to explain the formation of persistence in bacteria. We recently hypothesized that all physiological culture conditions leading to persistence converge in the inability of the bacteria to re-initiate a new round of DNA replication caused by an insufficient level of the initiator complex ATP-DnaA and hence by the lack of formation of a functional orisome. Here, we extend this hypothesis by proposing that in this persistence state the bacteria become more susceptible to mutation-based antibiotic resistance provided they are equipped with error-prone DNA repair functions. This is - in our opinion - in particular the case when such bacterial populations are exposed to BCAs.

Combined read- and assembly-based metagenomics to reconstruct a Dehalococcoides mccartyi genome from PCB-contaminated sediments and evaluate functional differences among organohalide-respiring consortia in the presence of different halogenated contaminants

Microbial communities that support respiration of halogenated organic contaminants by Dehalococcoides sp. facilitate full-scale bioremediation of chlorinated ethenes and demonstrate the potential to aid in bioremediation of halogenated aromatics like polychlorinated biphenyls (PCBs). However, it remains unclear if Dehalococcoides-containing microbial community dynamics observed in sediment-free systems quantitatively resemble that of sediment environments. To evaluate that possibility we assembled, annotated, and analyzed a Dehalococcoides sp. metagenome-assembled genome (MAG) from PCB-contaminated sediments. Phylogenetic analysis of reductive dehalogenase gene (rdhA) sequences within the MAG revealed that pcbA1 and pcbA4/5-like rdhA were absent, while several candidate PCB dehalogenase genes and potentially novel rdhA sequences were identified. Using a compositional comparative metagenomics approach, we quantified Dehalococcoides-containing microbial community structure shifts in response to halogenated organics and the presence of sediments. Functional level analysis revealed significantly greater abundances of genes associated with cobamide remodeling and horizontal gene transfer in tetrachloroethene-fed cultures as compared to halogenated aromatic-exposed consortia with or without sediments, despite little evidence of statistically significant differences in microbial community taxonomic structure. Our findings support the use of a generalizable comparative metagenomics workflow to evaluate Dehalococcoides-containing consortia in sediments and sediment-free environments to eludicate functions and microbial interactions that facilitate bioremediation of halogenated organic contaminants.

The resistance within: Antibiotic disruption of the gut microbiome and resistome dynamics in infancy

Intestinal host-microbiota interactions during the first year of life are critical for infant development. Early-life antibiotic exposures disrupt stereotypical gut microbiota maturation and adversely affect childhood health. Furthermore, antibiotics increase the abundance of resistant bacteria and enrich the resistome-the compendium of antibiotic resistance genes-within the gut microbiota. Here, we discuss acute and persistent impacts of antibiotic exposure during infancy on pediatric health, the gut microbiome, and, particularly, the resistome. Reviewing our current understanding of antibiotic resistance acquisition and dissemination within and between microbiomes, we highlight open questions, which are imperative to resolve in the face of rising bacterial resistance.

Good microbes, bad genes? The dissemination of antimicrobial resistance in the human microbiome

A global rise in antimicrobial resistance among pathogenic bacteria has proved to be a major public health threat, with the rate of multidrug-resistant bacterial infections increasing over time. The gut microbiome has been studied as a reservoir of antibiotic resistance genes (ARGs) that can be transferred to bacterial pathogens via horizontal gene transfer (HGT) of conjugative plasmids and mobile genetic elements (the gut resistome). Advances in metagenomic sequencing have facilitated the identification of resistome modulators, including live microbial therapeutics such as probiotics and fecal microbiome transplantation that can either expand or reduce the abundances of ARG-carrying bacteria in the gut. While many different gut microbes encode for ARGs, they are not uniformly distributed across, or transmitted by, various members of the microbiome, and not all are of equal clinical relevance. Both experimental and theoretical approaches in microbial ecology have been applied to understand differing frequencies of ARG horizontal transfer between commensal microbes as well as between commensals and pathogens. In this commentary, we assess the evidence for the role of commensal gut microbes in encoding antimicrobial resistance genes, the degree to which they are shared both with other commensals and with pathogens, and the host and environmental factors that can impact resistome dynamics. We further discuss novel sequencing-based approaches for identifying ARGs and predicting future transfer events of clinically relevant ARGs from commensals to pathogens.