ACLAME: a CLAssification of Mobile genetic Elements, update 2010

Water volume influences antibiotic resistomes and microbiomes during fish corpse decomposition

Corpse decomposition may cause serious pollution (e.g., releasing antibiotic resistance genes) to the water environment, thereby threatening public health. However, whether antibiotic resistance genes (ARGs) and microbiomes are affected by different water volumes during carcass decomposition remains unknown. Here, we investigated the effects of large/small water volumes on microbial communities and ARGs during fish cadaver decomposition by 16S rRNA high-throughput sequencing and high-throughput quantitative PCR. The results showed that the large water volume almost eliminated the effects of corpse decomposition on pH, total organic carbon (TOC), and total nitrogen (TN). When the water volume enlarged by 62.5 fold, the relative abundances of some ARGs resisting tetracycline and sulfonamide during carcass decomposition decreased by 217 fold on average, while there was also a mean 5267 fold increase of vancomycin resistance genes. Compared with the control group, the enriched types of ARGs varied between the large and small volume. Water volume, mobile genetic elements, and carcass decomposition were the most important factors affecting ARG profiles. Many opportunistic pathogens (like Bacteroides and Comamonas) were enriched in the corpse group. Bacteroides and Comamonas may be potential hosts of ARGs, indicating the potential for the spread of ARGs to humans by water pathogenic bacteria. This research highlights that the "dilution effect" can contribute to eliminating this adverse effect during corpse decomposition to a certain extent. It may provide references for environmental governance and public health.

Cross-Sectional Study on the Gut Microbiome of Parkinson's Disease Patients in Central China

Gastrointestinal dysfunction plays an important role in the occurrence and development of Parkinson's disease (PD). This study investigates the composition of the gut microbiome using shotgun metagenomic sequencing in PD patients in central China. Fecal samples from 39 PD patients (PD group) and the corresponding 39 healthy spouses of the patients (SP) were collected for shotgun metagenomics sequencing. Results showed a significantly altered microbial composition in the PD patients. Bilophila wadsworthia enrichment was found in the gut microbiome of PD patients, which has not been reported in previous studies. The random forest (RF) model, which identifies differences in microbiomes, reliably discriminated patients with PD from controls; the area under the receiver operating characteristic curve was 0.803. Further analysis of the microbiome and clinical symptoms showed that Klebsiella and Parasutterella were positively correlated with the duration and severity of PD, whereas hydrogen-generating Prevotella was negatively correlated with disease severity. The Cluster of Orthologous Groups of protein database, the KEGG Orthology database, and the carbohydrate-active enzymes of gene-category analysis showed that branched-chain amino acid-related proteins were significantly increased, and GH43 was significantly reduced in the PD group. Functional analysis of the metagenome confirmed differences in microbiome metabolism in the PD group related to short-chain fatty acid precursor metabolism.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

PHROG: families of prokaryotic virus proteins clustered using remote homology

Viruses are abundant, diverse and ancestral biological entities. Their diversity is high, both in terms of the number of different protein families encountered and in the sequence heterogeneity of each protein family. The recent increase in sequenced viral genomes constitutes a great opportunity to gain new insights into this diversity and consequently urges the development of annotation resources to help functional and comparative analysis. Here, we introduce PHROG (Prokaryotic Virus Remote Homologous Groups), a library of viral protein families generated using a new clustering approach based on remote homology detection by HMM profile-profile comparisons. Considering 17 473 reference (pro)viruses of prokaryotes, 868 340 of the total 938 864 proteins were grouped into 38 880 clusters that proved to be a 2-fold deeper clustering than using a classical strategy based on BLAST-like similarity searches, and yet to remain homogeneous. Manual inspection of similarities to various reference sequence databases led to the annotation of 5108 clusters (containing 50.6 % of the total protein dataset) with 705 different annotation terms, included in 9 functional categories, specifically designed for viruses. Hopefully, PHROG will be a useful tool to better annotate future prokaryotic viral sequences thus helping the scientific community to better understand the evolution and ecology of these entities.

Examining horizontal gene transfer in microbial communities

Bacteria acquire novel DNA through horizontal gene transfer (HGT), a process that enables an organism to rapidly adapt to changing environmental conditions, provides a competitive edge and potentially alters its relationship with its host. Although the HGT process is routinely exploited in laboratories, there is a surprising disconnect between what we know from laboratory experiments and what we know from natural environments, such as the human gut microbiome. Owing to a suite of newly available computational algorithms and experimental approaches, we have a broader understanding of the genes that are being transferred and are starting to understand the ecology of HGT in natural microbial communities. This Review focuses on these new technologies, the questions they can address and their limitations. As these methods are applied more broadly, we are beginning to recognize the full extent of HGT possible within a microbiome and the punctuated dynamics of HGT, specifically in response to external stimuli. Furthermore, we are better characterizing the complex selective pressures on mobile genetic elements and the mechanisms by which they interact with the bacterial host genome.

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

digIS: towards detecting distant and putative novel insertion sequence elements in prokaryotic genomes

Background

The insertion sequence elements (IS elements) represent the smallest and the most abundant mobile elements in prokaryotic genomes. It has been shown that they play a significant role in genome organization and evolution. To better understand their function in the host genome, it is desirable to have an effective detection and annotation tool. This need becomes even more crucial when considering rapid-growing genomic and metagenomic data. The existing tools for IS elements detection and annotation are usually based on comparing sequence similarity with a database of known IS families. Thus, they have limited ability to discover distant and putative novel IS elements.

Results

In this paper, we present digIS, a software tool based on profile hidden Markov models assembled from catalytic domains of transposases. It shows a very good performance in detecting known IS elements when tested on datasets with manually curated annotation. The main contribution of digIS is in its ability to detect distant and putative novel IS elements while maintaining a moderate level of false positives. In this category it outperforms existing tools, especially when tested on large datasets of archaeal and bacterial genomes.

Conclusion

We provide digIS, a software tool using a novel approach based on manually curated profile hidden Markov models, which is able to detect distant and putative novel IS elements. Although digIS can find known IS elements as well, we expect it to be used primarily by scientists interested in finding novel IS elements. The tool is available at https://github.com/janka2012/digIS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04177-6.

Infection with Bacteroides Phage BV01 Alters the Host Transcriptome and Bile Acid Metabolism in a Common Human Gut Microbe

Gut-associated phages are hypothesized to alter the abundance and activity of their bacterial hosts, contributing to human health and disease. Although temperate phages constitute a significant fraction of the gut virome, the effects of lysogenic infection are underexplored. We report that the temperate phage, Bacteroides phage BV01, broadly alters its host's transcriptome, the prominent human gut symbiont Bacteroides vulgatus. This alteration occurs through phage-induced repression of a tryptophan-rich sensory protein (TspO) and represses bile acid deconjugation. Because microbially modified bile acids are important signals for the mammalian host, this is a mechanism by which a phage may influence mammalian phenotypes. Furthermore, BV01 and its relatives in the proposed phage family Salyersviridae are ubiquitous in human gut metagenomes, infecting a broad range of Bacteroides hosts. These results demonstrate the complexity of phage-bacteria-mammal relationships and emphasize a need to better understand the role of temperate phages in the gut microbiome.

Sequence analysis of tyrosine recombinases allows annotation of mobile genetic elements in prokaryotic genomes

Mobile genetic elements (MGEs) sequester and mobilize antibiotic resistance genes across bacterial genomes. Efficient and reliable identification of such elements is necessary to follow resistance spreading. However, automated tools for MGE identification are missing. Tyrosine recombinase (YR) proteins drive MGE mobilization and could provide markers for MGE detection, but they constitute a diverse family also involved in housekeeping functions. Here, we conducted a comprehensive survey of YRs from bacterial, archaeal, and phage genomes and developed a sequence-based classification system that dissects the characteristics of MGE-borne YRs. We revealed that MGE-related YRs evolved from non-mobile YRs by acquisition of a regulatory arm-binding domain that is essential for their mobility function. Based on these results, we further identified numerous unknown MGEs. This work provides a resource for comparative analysis and functional annotation of YRs and aids the development of computational tools for MGE annotation. Additionally, we reveal how YRs adapted to drive gene transfer across species and provide a tool to better characterize antibiotic resistance dissemination.

Novel Siphoviridae Bacteriophages Infecting Bacteroides uniformis Contain Diversity Generating Retroelement

Intestinal phages are abundant and important components of gut microbiota, yet the isolated and characterized representatives that infect abundant gut bacteria are sparse. Here we describe the isolation of human intestinal phages infecting Bacteroidesuniformis. Bacteroides is one of the most common bacterial groups in the global human gut microbiota; however, to date not many Bacteroides specific phages are known. Phages isolated in this study belong to a novel viral genus, Bacuni, within the Siphoviridae family. Their genomes encode diversity-generating retroelements (DGR), which were shown in other bacteriophages to promote phage adaptation to rapidly changing environmental conditions and to broaden their host range. Three isolated phages showed 99.83% genome identity but one of them infected a distinct B. uniformis strain. The tropism of Bacuni phages appeared to be dependent on the interplay of DGR mediated sequence variations of gene encoding putative phage fimbrial tip proteins and mutations in host genes coding for outer-membrane proteins. We found prophages with up to 85% amino acid similarity over two-thirds of the Bacuni phage genome in the B. acidifaciens and Prevotella sp. genomes. Despite the abundance of Bacteroides within the human microbiome, we found Bacuni phages only in a limited subset of published gut metagenomes.

Expanding Diversity of Firmicutes Single-Strand Annealing Proteins: A Putative Role of Bacteriophage-Host Arms Race

Bacteriophage-encoded single strand annealing proteins (SSAPs) are recombinases which can substitute the classical, bacterial RecA and manage the DNA metabolism at different steps of phage propagation. SSAPs have been shown to efficiently promote recombination between short and rather divergent DNA sequences and were exploited for in vivo genetic engineering mainly in Gram-negative bacteria. In opposition to the conserved and almost universal bacterial RecA protein, SSAPs display great sequence diversity. The importance for SSAPs in phage biology and phage-bacteria evolution is underlined by their role as key players in events of horizontal gene transfer (HGT). All of the above provoke a constant interest for the identification and study of new phage recombinase proteins in vivo, in vitro as well as in silico. Despite this, a huge body of putative ssap genes escapes conventional classification, as they are not properly annotated. In this work, we performed a wide-scale identification, classification and analysis of SSAPs encoded by the Firmicutes bacteria and their phages. By using sequence similarity network and gene context analyses, we created a new high quality dataset of phage-related SSAPs, substantially increasing the number of annotated SSAPs. We classified the identified SSAPs into seven distinct families, namely RecA, Gp2.5, RecT/Redβ, Erf, Rad52/22, Sak3, and Sak4, organized into three superfamilies. Analysis of the relationships between the revealed protein clusters led us to recognize Sak3-like proteins as a new distinct SSAP family. Our analysis showed an irregular phylogenetic distribution of ssap genes among different bacterial phyla and specific phages, which can be explained by the high rates of ssap HGT. We propose that the evolution of phage recombinases could be tightly linked to the dissemination of bacterial phage-resistance mechanisms (e.g., abortive infection and CRISPR/Cas systems) targeting ssap genes and be a part of the constant phage-bacteria arms race.

Next generation sequencing approaches to evaluate water and wastewater quality

The emergence of next generation sequencing (NGS) is revolutionizing the potential to address complex microbiological challenges in the water industry. NGS technologies can provide holistic insight into microbial communities and their functional capacities in water and wastewater systems, thus eliminating the need to develop a new assay for each target organism or gene. However, several barriers have hampered wide-scale adoption of NGS by the water industry, including cost, need for specialized expertise and equipment, challenges with data analysis and interpretation, lack of standardized methods, and the rapid pace of development of new technologies. In this critical review, we provide an overview of the current state of the science of NGS technologies as they apply to water, wastewater, and recycled water. In addition, a systematic literature review was conducted in which we identified over 600 peer-reviewed journal articles on this topic and summarized their contributions to six key areas relevant to the water and wastewater fields: taxonomic classification and pathogen detection, functional and catabolic gene characterization, antimicrobial resistance (AMR) profiling, bacterial toxicity characterization, Cyanobacteria and harmful algal bloom identification, and virus characterization. For each application, we have presented key trends, noteworthy advancements, and proposed future directions. Finally, key needs to advance NGS technologies for broader application in water and wastewater fields are assessed.

Novel Soil-Derived Beta-Lactam, Chloramphenicol, Fosfomycin and Trimethoprim Resistance Genes Revealed by Functional Metagenomics

Antibiotic resistance genes (ARGs) in soil are considered to represent one of the largest environmental resistomes on our planet. As these genes can potentially be disseminated among microorganisms via horizontal gene transfer (HGT) and in some cases are acquired by clinical pathogens, knowledge about their diversity, mobility and encoded resistance spectra gained increasing public attention. This knowledge offers opportunities with respect to improved risk prediction and development of strategies to tackle antibiotic resistance, and might help to direct the design of novel antibiotics, before further resistances reach hospital settings or the animal sector. Here, metagenomic libraries, which comprise genes of cultivated microorganisms, but, importantly, also those carried by the uncultured microbial majority, were screened for novel ARGs from forest and grassland soils. We detected three new beta-lactam, a so far unknown chloramphenicol, a novel fosfomycin, as well as three previously undiscovered trimethoprim resistance genes. These ARGs were derived from phylogenetically diverse soil bacteria and predicted to encode antibiotic inactivation, antibiotic efflux, or alternative variants of target enzymes. Moreover, deduced gene products show a minimum identity of ~21% to reference database entries and confer high-level resistance. This highlights the vast potential of functional metagenomics for the discovery of novel ARGs from soil ecosystems.

Diversification of methanogens into hyperalkaline serpentinizing environments through adaptations to minimize oxidant limitation

Metagenome assembled genomes (MAGs) and single amplified genomes (SAGs) affiliated with two distinct Methanobacterium lineages were recovered from subsurface fracture waters of the Samail Ophiolite, Sultanate of Oman. Lineage Type I was abundant in waters with circumneutral pH, whereas lineage Type II was abundant in hydrogen rich, hyperalkaline waters. Type I encoded proteins to couple hydrogen oxidation to CO₂ reduction, typical of hydrogenotrophic methanogens. Surprisingly, Type II, which branched from the Type I lineage, lacked homologs of two key oxidative [NiFe]-hydrogenases. These functions were presumably replaced by formate dehydrogenases that oxidize formate to yield reductant and cytoplasmic CO₂ via a pathway that was unique among characterized Methanobacteria, allowing cells to overcome CO₂/oxidant limitation in high pH waters. This prediction was supported by microcosm-based radiotracer experiments that showed significant biological methane generation from formate, but not bicarbonate, in waters where the Type II lineage was detected in highest relative abundance. Phylogenetic analyses and variability in gene content suggested that recent and ongoing diversification of the Type II lineage was enabled by gene transfer, loss, and transposition. These data indicate that selection imposed by CO₂/oxidant availability drove recent methanogen diversification into hyperalkaline waters that are heavily impacted by serpentinization.

Combined effects of composting and antibiotic administration on cattle manure-borne antibiotic resistance genes

BACKGROUND

Research is needed to delineate the relative and combined effects of different antibiotic administration and manure management practices in either amplifying or attenuating the potential for antibiotic resistance to spread. Here, we carried out a comprehensive parallel examination of the effects of small-scale (> 55 °C × 3 days) static and turned composting of manures from dairy and beef cattle collected during standard antibiotic administration (cephapirin/pirlimycin or sulfamethazine/chlortetracycline/tylosin, respectively), versus from untreated cattle, on "resistomes" (total antibiotic resistance genes (ARGs) determined via shotgun metagenomic sequencing), bacterial microbiota, and indicator ARGs enumerated via quantitative polymerase chain reaction. To gain insight into the role of the thermophilic phase, compost was also externally heated to > 55 °C × 15 days.

RESULTS

Progression of composting with time and succession of the corresponding bacterial microbiota was the overarching driver of the resistome composition (ANOSIM; R = 0.424, p = 0.001, respectively) in all composts at the small-scale. Reduction in relative abundance (16S rRNA gene normalized) of total ARGs in finished compost (day 42) versus day 0 was noted across all conditions (ANOSIM; R = 0.728, p = 0.001), except when externally heated. Sul1, intI1, beta-lactam ARGs, and plasmid-associated genes increased in all finished composts as compared with the initial condition. External heating more effectively reduced certain clinically relevant ARGs (bla_OXA, bla_CARB), fecal coliforms, and resistome risk scores, which take into account putative pathogen annotations. When manure was collected during antibiotic administration, taxonomic composition of the compost was distinct according to nonmetric multidimensional analysis and tet(W) decayed faster in the dairy manure with antibiotic condition and slower in the beef manure with antibiotic condition.

CONCLUSIONS

This comprehensive, integrated study revealed that composting had a dominant effect on corresponding resistome composition, while little difference was noted as a function of collecting manure during antibiotic administration. Reduction in total ARGs, tet(W), and resistome risk suggested that composting reduced some potential for antibiotic resistance to spread, but the increase and persistence of other indicators of antibiotic resistance were concerning. Results indicate that composting guidelines intended for pathogen reduction do not necessarily provide a comprehensive barrier to ARGs or their mobility prior to land application and additional mitigation measures should be considered. Video Abstract.

Viruses in fermented foods: are they good or bad? Two sides of the same coin

The emergence of Coronavirus disease 2019 as a global pandemic has increased popular concerns about diseases caused by viruses. Fermented foods containing high loads of viable fungi and bacteria are potential sources for virus contamination. The most common include viruses that infect bacteria (bacteriophage) and yeasts reported in fermented milks, sausages, vegetables, wine, sourdough, and cocoa beans. Recent molecular studies have also associated fermented foods as vehicles for pathogenic human viruses. Human noroviruses, rotavirus, and hepatitis virus have been identified in different fermented foods through multiple routes. No severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) virus or close members were found in fermented foods to date. However, the occurrence/persistence of other pathogenic viruses reveals a potential vulnerability of fermented foods to SARS-CoV-2 contamination. On the other side of the coin, some bacteriophages are being suggested for improving the fermentation process and food safety, as well as owing potential probiotic properties in modern fermented foods. This review will address the diversity and characteristics of viruses associated with fermented foods and what has been changed after a short introduction to the most common next-generation sequencing platforms. Also, the risk of SARS-CoV-2 transmission via fermented foods and preventive measures will be discussed.

Metagenomic analysis of microbiota and antibiotic resistome in household activated carbon drinking water purifiers

Existing drinking water treatment systems have limited ability to control emerging contaminants such as antibiotic resistance genes (ARGs). Household activated carbon water purifiers (HWPs) are convenient measures to assure drinking water quality. However, ARGs distribution in HWPs has not been reported. Here, ARGs, mobile genetic elements (MGEs) and bacteria communities were profiled in tap water (TW), filter water (FW) and activated carbon (AC) biofilm from six kinds of HWPs after 80 days operation, using metagenomics. Results showed that the bacteria community diversities in FW and AC were higher than those in TW. A total of 88, 116 and 80 ARG subtypes were detected in TW, AC and FW, respectively. The AC structure was an important factor influencing the bacterial communities and ARG profiles in FW. The network analysis revealed the co-occurrence patterns between ARGs and bacteria. SourceTracker analyses showed AC biofilms were important contributors of microbes (29-79%) and ARGs (17-53%) in FW. Moreover, MGEs e.g. pBBta01, pMKMS02 and pMFLV01 plasmids, and ISMysp3 had significant co-occurrence patterns with ARGs in the AC biofilms. This study helps to understand the actual purification effect of HWPs and provides a theoretical reference for the management and control of ARGs pollution in domestic drinking water.

Critical evaluation of short, long, and hybrid assembly for contextual analysis of antibiotic resistance genes in complex environmental metagenomes

In the fight to limit the global spread of antibiotic resistance, the assembly of environmental metagenomes has the potential to provide rich contextual information (e.g., taxonomic hosts, carriage on mobile genetic elements) about antibiotic resistance genes (ARG) in the environment. However, computational challenges associated with assembly can impact the accuracy of downstream analyses. This work critically evaluates the impact of assembly leveraging short reads, nanopore MinION long-reads, and a combination of the two (hybrid) on ARG contextualization for ten environmental metagenomes using seven prominent assemblers (IDBA-UD, MEGAHIT, Canu, Flye, Opera-MS, metaSpades and HybridSpades). While short-read and hybrid assemblies produced similar patterns of ARG contextualization, raw or assembled long nanopore reads produced distinct patterns. Based on an in-silico spike-in experiment using real and simulated reads, we show that low to intermediate coverage species are more likely to be incorporated into chimeric contigs across all assemblers and sequencing technologies, while more abundant species produce assemblies with a greater frequency of inversions and insertion/deletions (indels). In sum, our analyses support hybrid assembly as a valuable technique for boosting the reliability and accuracy of assembly-based analyses of ARGs and neighboring genes at environmentally-relevant coverages, provided that sufficient short-read sequencing depth is achieved.

Free-floating extracellular DNA: Systematic profiling of mobile genetic elements and antibiotic resistance from wastewater

The free-floating extracellular DNA (exDNA) fraction of microbial ecosystems harbors antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs). Natural transformation of these xenogenetic elements can generate microbial cells resistant to one or more antibiotics. Isolating and obtaining a high yield of exDNA is challenging due to its low concentration in wastewater environments. Profiling exDNA is crucial to unravel the ecology of free-floating ARGs and MGEs and their contribution to horizontal genetransfer. We developed a method using chromatography to isolate and enrich exDNA without causing cell lysis from complex wastewater matrices like influent (9 µg exDNA out of 1 L), activated sludge (5.6 µg out of 1 L), and treated effluent (4.3 µg out of 1 L). ARGs and MGEs were metagenomically profiled for both the exDNA and intracellular DNA (iDNA) of activated sludge, and quantified by qPCR in effluent water. qPCR revealed that ARGs and MGEs are more abundant in the iDNA fraction while still significant on exDNA (100-1000 gene copies mL^-1) in effluent water. The metagenome highlighted that exDNA is mainly composed of MGEs (65%). According to their relatively low abundance in the resistome of exDNA, ARGs uptake by natural transformation is likely not the main transfer mechanism. Although ARGs are not highly abundant in exDNA, the prevalence of MGEs in the exDNA fraction can indirectly promote antibiotic resistance development. The combination of this method with functional metagenomics can help to elucidate the transfer and development of resistances in microbial communities. A systematic profiling of the different DNA fractions will foster microbial risk assessments across water systems, supporting water authorities to delineate measures to safeguard environmental and public health.

Acquisition of pcnB [poly(A) polymerase I] genes via horizontal transfer from the β, γ-Proteobacteria

Poly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3'-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ-Proteobacteria. Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε-Proteobacteria and from four other bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes and Aquificae). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ-Proteobacteria.

Environmental risk characterization and ecological process determination of bacterial antibiotic resistome in lake sediments

The increasing prevalence of antibiotic resistance genes (ARGs) in aquatic environments has attracted considerable concerns due to their potential threat to public health. For reducing environmental risk of ARGs, it is crucial to identify the pathogenic resistant bacteria, determine the driving forces governing the ARG community and apportion their sources, which is yet remained to explore. In this study, we developed a framework integrating high-throughput sequencing (HTS) analyses, null-model-based methods and machine-learning classification tool for understanding the environmental resistome risk and the ecological processes that control the ARG profile in aquatic sediments, and applied to two urban lakes (Lake Tai and Lake Baiyang) in China. The HTS-based metagenomic analyses revealed abundant and diverse resistome, mobilome and virulome in the two lakes, including some emerging ARGs such as mcr and carbapenemases types. Relatively, the diversities for ARGs, mobile genetic elements (MGEs) and virulence factor genes in Lake Baiyang were significantly higher than those in Lake Tai (p < 0.05). The metagenomic assembly and binning approaches tracked a number of potential pathogenic antibiotic resistant bacteria and found the co-occurrence of ARGs, MGEs and human bacterial pathogens in ~50% of the sediment samples, indicating a substantial resistome risk in the lakes. Comparison of multiple-site beta-diversity dissimilarity indexes suggested the ARG diversity was mainly explained by the spatial turnover rather than nestedness and exhibited significant distance-decay pattern. The results of using a novel null-model-based stochasticity ratio showed the stochastic processes made a higher contribution than the deterministic processes on the ARG profile in the environment, especially for Lake Baiyang (>65%). This was confirmed by the determination analyses of various ecological processes on ARG community by utilizing the null-model-based statistical framework for quantifying community assembly. That is, homogenizing dispersal (40%) dominated in Lake Baiyang, followed by homogeneous selection (32%) and ecological drift (15%), while ecological drift (33%) and homogenizing dispersal (31%) were the dominators in Lake Baiyang. SourceTracker analysis showed human sewage-associated sources were the largest contributor (~62%) of ARGs in the environment. The findings shed light on the dissemination risk and driver dynamics of antimicrobial resistance in the aquatic environment, which may help to make effective management strategies for controlling pollution of ARGs.

Will a Non-antibiotic Metalloid Enhance the Spread of Antibiotic Resistance Genes: The Selenate Story

The rapid emergence of antibiotic resistance genes (ARGs) has become an increasingly serious threat to public health. Previous studies illustrate the antibiotic-like effect of many substances. However, whether and how commonly used or existing non-antibiotic metalloids (e.g., selenate) would enhance ARG spread remains poorly known. Here, we tracked the long-term operation of a bioreactor continuously fed with selenate for more than 1000 days. Metagenomic sequencing identified 191 different ARGs, of which the total abundance increased significantly after the amendment of selenate. Network analyses showed that ARGs resisting multiple drugs had very similar co-occurrence patterns, implying a potentially larger health risk. Host classification not only indicated multidrug-resistant species but also distinguished the mechanism of ARG enrichment for vertical transfer and horizontal gene transfer. Genome reconstruction of an ARG host suggested that selenate and its bioreduction product selenite could stimulate the overproduction of intracellular reactive oxygen species, which was confirmed by the direct measurement. Bacterial membrane permeability, type IV pilus formation, and DNA repair and recombination were also enhanced, together facilitating the horizontal acquirement of ARGs. Overall, this study for the first time highlights the ARG emergence and dissemination induced by a non-antibiotic metalloid and identifies ARG as a factor to consider in selenate bioremediation.

mMGE: a database for human metagenomic extrachromosomal mobile genetic elements

Extrachromosomal mobile genetic elements (eMGEs), including phages and plasmids, that can move across different microbes, play important roles in genome evolution and shaping the structure of microbial communities. However, we still know very little about eMGEs, especially their abundances, distributions and putative functions in microbiomes. Thus, a comprehensive description of eMGEs is of great utility. Here we present mMGE, a comprehensive catalog of 517 251 non-redundant eMGEs, including 92 492 plasmids and 424 759 phages, derived from diverse body sites of 66 425 human metagenomic samples. About half the eMGEs could be further grouped into 70 074 clusters using relaxed criteria (referred as to eMGE clusters below). We provide extensive annotations of the identified eMGEs including sequence characteristics, taxonomy affiliation, gene contents and their prokaryotic hosts. We also calculate the prevalence, both within and across samples for each eMGE and eMGE cluster, enabling users to see putative associations of eMGEs with human phenotypes or their distribution preferences. All eMGE records can be browsed or queried in multiple ways, such as eMGE clusters, metagenomic samples and associated hosts. The mMGE is equipped with a user-friendly interface and a BLAST server, facilitating easy access/queries to all its contents easily. mMGE is freely available for academic use at: https://mgedb.comp-sysbio.org.

Genome-Resolved Metagenomics and Antibiotic Resistance Genes Analysis in Reclaimed Water Distribution Systems

Water reuse is increasingly pursued to alleviate global water scarcity. However, the wastewater treatment process does not achieve full removal of biological contaminants from wastewater, hence microorganisms and their genetic elements can be disseminated into the reclaimed water distribution systems (RWDS). In this study, reclaimed water samples are investigated via metagenomics to assess their bacterial diversity, metagenome-assembled genomes (MAGs) and antibiotic resistance genes (ARGs) at both point of entry (POE) and point of use (POU) in 3 RWDS. The number of shared bacterial orders identified by metagenome was higher at the POE than POU among the three sites, indicating that specific conditions in RWDS can cause further differentiation in the microbial communities at the end of the distribution system. Two bacterial orders, namely Rhizobiales and Sphingomonadales, had high replication rates in two of the examined RWDS (i.e., site A and B), and were present in higher relative abundance in POU than at POE. In addition, MAG and ARG relative abundance exhibited a strong correlation (R2 = 0.58) in POU, indicating that bacteria present in POU may have a high incidence of ARG. Specifically, resistance genes associated with efflux pump mechanisms (e.g., adeF and qacH) increased in its relative abundance from POU to POE at two of the RWDS (i.e., site A and B). When correlated with the water quality data that suggests a significantly lower dissolved organic carbon (DOC) concentration at site D than the other two RWDS, the metagenomic data suggest that low DOC is needed to maintain the biological stability of reclaimed water along the distribution network.

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Integrated into their bacterial hosts' genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches-bioinformatics, analytical modelling and computational simulation-to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.

Demonstrating an Integrated Antibiotic Resistance Gene Surveillance Approach in Puerto Rican Watersheds Post-Hurricane Maria

Comprehensive surveillance approaches are needed to assess sources, clinical relevance, and mobility of antibiotic resistance genes (ARGs) in watersheds. Here, we examined metrics derived from shotgun metagenomic sequencing and relationship to human fecal markers (HFMs; crAssphage, enterococci) and anthropogenic antibiotic resistance markers (AARMs; intI1, sul1) in three distinct Puerto Rican watersheds as a function of adjacent land use and wastewater treatment plant (WWTP) input 6 months after Hurricane Maria, a category V storm. Relative abundance and diversity of total ARGs increased markedly downstream of WWTP inputs, with ARGs unique to WWTP and WWTP-impacted river samples predominantly belonging to the aminoglycoside and β-lactam resistance classes. WWTP and other anthropogenic inputs were similarly associated with elevated resistome risk scores and mobility incidence (M%). Contig analysis indicated a wide variety of mobile β-lactam ARGs associated with pathogens downstream of WWTP discharge that were consistent with regional clinical concern, e.g., Klebsiella pneumoniae contigs containing KPC-2 within an ISKpn6-like transposase. HFMs and AARMs correlated strongly with the absolute abundance of total ARGs, but AARMs better predicted the majority of ARGs in general (85.4 versus <2%) and β-lactam ARGs in particular. This study reveals sensitive, quantitative, mobile, clinically relevant, and comprehensive targets for antibiotic resistance surveillance in watersheds.

Localized effect of treated wastewater effluent on the resistome of an urban watershed

Background

Wastewater treatment is an essential tool for maintaining water quality in urban environments. While the treatment of wastewater can remove most bacterial cells, some will inevitably survive treatment to be released into natural environments. Previous studies have investigated antibiotic resistance within wastewater treatment plants, but few studies have explored how a river’s complete set of antibiotic resistance genes (the “resistome") is affected by the release of treated effluent into surface waters.

Results

Here we used high-throughput, deep metagenomic sequencing to investigate the effect of treated wastewater effluent on the resistome of an urban river and the downstream distribution of effluent-associated antibiotic resistance genes and mobile genetic elements. Treated effluent release was found to be associated with increased abundance and diversity of antibiotic resistance genes and mobile genetic elements. The impact of wastewater discharge on the river’s resistome diminished with increasing distance from effluent discharge points. The resistome at river locations that were not immediately downstream from any wastewater discharge points was dominated by a single integron carrying genes associated with resistance to sulfonamides and quaternary ammonium compounds.

Conclusions

Our study documents variations in the resistome of an urban watershed from headwaters to a major confluence in an urban center. Greater abundances and diversity of antibiotic resistance genes are associated with human fecal contamination in river surface water, but the fecal contamination effect seems to be localized, with little measurable effect in downstream waters. The diverse composition of antibiotic resistance genes throughout the watershed suggests the influence of multiple environmental and biological factors.

The Gut Virome Database Reveals Age-Dependent Patterns of Virome Diversity in the Human Gut

The gut microbiome profoundly affects human health and disease, and their infecting viruses are likely as important, but often missed because of reference database limitations. Here, we (1) built a human Gut Virome Database (GVD) from 2,697 viral particle or microbial metagenomes from 1,986 individuals representing 16 countries, (2) assess its effectiveness, and (3) report a meta-analysis that reveals age-dependent patterns across healthy Westerners. The GVD contains 33,242 unique viral populations (approximately species-level taxa) and improves average viral detection rates over viral RefSeq and IMG/VR nearly 182-fold and 2.6-fold, respectively. GVD meta-analyses show highly personalized viromes, reveal that inter-study variability from technical artifacts is larger than any "disease" effect at the population level, and document how viral diversity changes from human infancy into senescence. Together, this compact foundational resource, these standardization guidelines, and these meta-analysis findings provide a systematic toolkit to help maximize our understanding of viral roles in health and disease.

Homogeneous selection drives antibiotic resistome in two adjacent sub-watersheds, China

Rivers are a significant reservoir of antibiotic resistance genes (ARGs), yet the biogeographic pattern of riverine ARGs and its underlying driving forces remain poorly understood. Here, we used metagenomic approach to investigate the spatio-temporal variation of ARGs in two adjacent sub-watersheds viz. North River (NR) and West River (WR), China. The results demonstrated that Bacitracin (22.8 % of the total ARGs), multidrug (20.7 %), sulfonamide (15.2 %) and tetracycline (10.9 %) were the dominant ARG types. SourceTracker analysis indicated that sewage treatment plants as the main source of ARGs, while animal feces mainly contributed in spreading the ARGs in the upstream of NR. Random forest and network analyses confirmed that NR was under the influence of fecal pollution. PCoA analysis demonstrated that the composition of ARGs changed along with the anthropogenic gradients, while the Raup-Crick null model showed that homogenizing selection mediated by class 1 integron intI1 resulted in stable ARG communities at whole watershed scale. Structural equation models revealed that microbial community, grassland and several non-antibiotic micropollutants may also play certain roles in influencing the distribution of ARGs. Overall, the observed deterministic formation of ARGs in riverine systems calls effective management strategies to mitigate the risks of antibiotic resistance on public health.

Gene Loss and Acquisition in Lineages of Pseudomonas aeruginosa Evolving in Cystic Fibrosis Patient Airways

Bacterial airway infections, predominantly caused by P. aeruginosa, are a major cause of mortality and morbidity of CF patients. While short insertions and deletions as well as point mutations occurring during infection are well studied, there is a lack of understanding of how gene loss and acquisition play roles in bacterial adaptation to the human airways. Here, we investigated P. aeruginosa within-host evolution with regard to gene loss and acquisition. We show that during long-term infection P. aeruginosa genomes tend to lose genes, in particular, genes related to virulence. This adaptive strategy allows reduction of the genome size and evasion of the host’s immune response. This knowledge is crucial to understand the basic mutational steps that, on the timescale of years, diversify lineages and adds to the identification of bacterial genetic determinants that have implications for CF disease.

Genome analyses have documented that there are differences in gene repertoire between evolutionary distant lineages of the same bacterial species; however, less is known about microevolutionary dynamics of gene loss and acquisition within bacterial lineages as they evolve over years. Here, we analyzed the genomes of 45 Pseudomonas aeruginosa lineages evolving in the lungs of cystic fibrosis (CF) patients to identify genes that are lost or acquired during the first years of infection. On average, lineage genome content changed by 88 genes (range, 0 to 473). Genes were more often lost than acquired, and prophage genes were more variable than bacterial genes. We identified convergent loss or acquisition of the same genes across lineages, suggesting selection for loss and acquisition of certain genes in the host environment. We found that a notable proportion of such genes are associated with virulence; a trait previously shown to be important for adaptation. Furthermore, we also compared the genomes across lineages to show that the within-lineage variable genes (i.e., genes that had been lost or acquired during the infection) often belonged to genomic content not shared across all lineages. In sum, our analysis adds to the knowledge on the pace and drivers of gene loss and acquisition in bacteria evolving over years in a human host environment and provides a basis to further understand how gene loss and acquisition play roles in lineage differentiation and host adaptation.

Structure of the putative long tail fiber receptor-binding tip of a novel temperate bacteriophage from the Antarctic bacterium Bizionia argentinensis JUB59

Tailed bacteriophages are one of the most widespread biological entities on Earth. Their singular structures, such as spikes or fibers are of special interest given their potential use in a wide range of biotechnological applications. In particular, the long fibers present at the termini of the T4 phage tail have been studied in detail and are important for host recognition and adsorption. Although significant progress has been made in elucidating structural mechanisms of model phages, the high-resolution structural description of the vast population of marine phages is still unexplored. In this context, we present here the crystal structure of C24, a putative receptor-binding tip-like protein from Bizionia argentinensis JUB59, a psychrotolerant bacterium isolated from the marine surface waters of Potter Cove, Antarctica. The structure resembles the receptor-binding tip from the bacteriophage T4 long tail fiber yet showing marked differences in its domain organization, size, sequence identity and metal binding nature. We confirmed the viral origin of C24 by induction experiments using mitomycin C. Our results reveal the presence of a novel uncharacterized prophage in the genome of B. argentinensis JUB59, whose morphology is compatible with the order Caudovirales and that carries the nucleotide sequence of C24 in its genome. This work provides valuable information to expand our current knowledge on the viral machinery prevalent in the oceans.

Characterization of putative circular plasmids in sponge-associated bacterial communities using a selective multiply-primed rolling circle amplification

Plasmid transfers among bacterial populations can directly influence the ecological adaptation of these populations and their interactions with host species and environment. In this study, we developed a selective multiply-primed rolling circle amplification (smRCA) approach to enrich and characterize circular plasmid DNA from sponge microbial symbionts via high-throughput sequencing (HTS). DNA (plasmid and total community DNA) obtained from sponge (Cinachyrella sp.) samples and a bacterial symbiont (Vibrio sp. CyArs1) isolated from the same sponge species (carrying unknown plasmids) were used to develop and validate our methodology. The smRCA was performed during 16 hr with 141 plasmid-specific primers covering all known circular plasmid groups. The amplified products were purified and subjected to a reamplification with random hexamer primers (2 hr) and then sequenced using Illumina MiSeq. The developed method resulted in the successful amplification and characterization of the sponge plasmidome and allowed us to detect plasmids associated with the bacterial symbiont Vibrio sp. CyArs1 in the sponge host. In addition to this, a large number of small (<2 kbp) and cryptic plasmids were also amplified in sponge samples. Functional analysis identified proteins involved in the control of plasmid partitioning, maintenance and replication. However, most plasmids contained unknown genes, which could potentially serve as a resource of unknown genetic information and novel replication systems. Overall, our results indicate that the smRCA-HTS approach developed here was able to selectively enrich and characterize plasmids from bacterial isolates and sponge host microbial communities, including plasmids larger than 20 kbp.

Isolation and Characterization of Nitrate-Reducing Bacteria as Potential Probiotics for Oral and Systemic Health

Recent evidence indicates that the reduction of salivary nitrate by oral bacteria can contribute to prevent oral diseases, as well as increase systemic nitric oxide levels that can improve conditions such as hypertension and diabetes. The objective of the current manuscript was to isolate nitrate-reducing bacteria from the oral cavity of healthy donors and test their in vitro probiotic potential to increase the nitrate-reduction capacity (NRC) of oral communities. Sixty-two isolates were obtained from five different donors of which 53 were confirmed to be nitrate-reducers. Ten isolates were selected based on high NRC as well as high growth rates and low acidogenicity, all being Rothia species. The genomes of these ten isolates confirmed the presence of nitrate- and nitrite reductase genes, as well as lactate utilization genes, and the absence of antimicrobial resistance, mobile genetic elements and virulence genes. The pH at which most nitrate was reduced differed between strains. However, acidic pH 6 always stimulated the reduction of nitrite compared to neutral pH 7 or slightly alkaline pH 7.5 (p < 0.01). We tested the effect of six out of 10 isolates on in vitro oral biofilm development in the presence or absence of 6.5 mM nitrate. The integration of the isolates into in vitro communities was confirmed by Illumina sequencing. The NRC of the bacterial communities increased when adding the isolates compared to controls without isolates (p < 0.05). When adding nitrate (prebiotic treatment) or isolates in combination with nitrate (symbiotic treatment), a smaller decrease in pH derived from sugar metabolism was observed (p < 0.05), which for some symbiotic combinations appeared to be due to lactate consumption. Interestingly, there was a strong correlation between the NRC of oral communities and ammonia production even in the absence of nitrate (R = 0.814, p < 0.01), which indicates that bacteria involved in these processes are related. As observed in our study, individuals differ in their NRC. Thus, some may have direct benefits from nitrate as a prebiotic as their microbiota naturally reduces significant amounts, while others may benefit more from a symbiotic combination (nitrate + nitrate-reducing probiotic). Future clinical studies should test the effects of these treatments on oral and systemic health.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Widespread transfer of mobile antibiotic resistance genes within individual gut microbiomes revealed through bacterial Hi-C

The gut microbiome harbors a ‘silent reservoir’ of antibiotic resistance (AR) genes that is thought to contribute to the emergence of multidrug-resistant pathogens through horizontal gene transfer (HGT). To counteract the spread of AR, it is paramount to know which organisms harbor mobile AR genes and which organisms engage in HGT. Despite methods that characterize the overall abundance of AR genes in the gut, technological limitations of short-read sequencing have precluded linking bacterial taxa to specific mobile genetic elements (MGEs) encoding AR genes. Here, we apply Hi-C, a high-throughput, culture-independent method, to surveil the bacterial carriage of MGEs. We compare two healthy individuals with seven neutropenic patients undergoing hematopoietic stem cell transplantation, who receive multiple courses of antibiotics, and are acutely vulnerable to the threat of multidrug-resistant infections. We find distinct networks of HGT across individuals, though AR and mobile genes are associated with more diverse taxa within the neutropenic patients than the healthy subjects. Our data further suggest that HGT occurs frequently over a several-week period in both cohorts. Whereas most efforts to understand the spread of AR genes have focused on pathogenic species, our findings shed light on the role of the human gut microbiome in this process.

Linking antibiotic resistance (AR) in the gut microbiome with their bacterial hosts remains challenging. Here, the authors apply bacterial Hi-C to map mobile genetic elements in metagenomes, and illustrate that genes are present in more diverse taxa in neutropenic patients than healthy subjects.

Horizontal transmission and recombination maintain forever young bacterial symbiont genomes

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation.

Environmental resistome risks of wastewaters and aquatic environments deciphered by shotgun metagenomic assembly

In this paper, we deciphered the core resistome disseminating from hospital wastewater to the aquatic environment by characterising the resistome, plasmidome, mobilome and virulome using metagenomic analysis. This study also elucidated different environmental resistome risks using shotgun-metagenomic assembly. The results showed that clinically relevant taxa were found in assessed matrices (Salmonella spp., Acinetobacter spp, Escherichia-Shigella spp., Pseudomonas spp., Staphylococcus spp. and Vibrio spp.). For the plasmidome, we found 249 core plasmidome sequences that were shared among all assessed matrices. The core mobilome of 2424 mobile genetic elements shared among all assessed matrices was found. Regarding the virulome, we found 148 core virulence factors shared among all assessed samples, and the core virulome content was consistently shared across the most abundant bacterial genera. Although influent of wastewater showed considerable higher relative bacterial abundance (P = 0.008), hospital wastewater showed significant higher environmental resistome risk scores against all other assessed matrices, with an average of 46.34% (P = 0.001). These results suggest hospital wastewater, effluent and sewage sludge should be subjected to stringent mitigating measures to minimise such dissemination.

Mobilization of Antibiotic Resistance: Are Current Approaches for Colocalizing Resistomes and Mobilomes Useful?

Antimicrobial resistance (AMR) poses a global human and animal health threat, and predicting AMR persistence and transmission remains an intractable challenge. Shotgun metagenomic sequencing can help overcome this by enabling characterization of AMR genes within all bacterial taxa, most of which are uncultivatable in laboratory settings. Shotgun sequencing, therefore, provides a more comprehensive glance at AMR "potential" within samples, i.e., the "resistome." However, the risk inherent within a given resistome is predicated on the genomic context of various AMR genes, including their presence within mobile genetic elements (MGEs). Therefore, resistome risk stratification can be advanced if AMR profiles are considered in light of the flanking mobilizable genomic milieu (e.g., plasmids, integrative conjugative elements (ICEs), phages, and other MGEs). Because such mediators of horizontal gene transfer (HGT) are involved in uptake by pathogens, investigators are increasingly interested in characterizing that resistome fraction in genomic proximity to HGT mediators, i.e., the "mobilome"; we term this "colocalization." We explored the utility of common colocalization approaches using alignment- and assembly-based techniques, on clinical (human) and agricultural (cattle) fecal metagenomes, obtained from antimicrobial use trials. Ordination revealed that tulathromycin-treated cattle experienced a shift in ICE and plasmid composition versus untreated animals, though the resistome was unaffected during the monitoring period. Contrarily, the human resistome and mobilome composition both shifted shortly after antimicrobial administration, though this rebounded to pre-treatment status. Bayesian networks revealed statistical AMR-MGE co-occurrence in 19 and 2% of edges from the cattle and human networks, respectively, suggesting a putatively greater mobility potential of AMR in cattle feces. Conversely, using Mobility Index (MI) and overlap analysis, abundance of de novo-assembled contigs supporting resistomes flanked by MGE increased shortly post-exposure within human metagenomes, though > 40 days after peak dose such contigs were rare (∼2%). MI was not substantially altered by antimicrobial exposure across all cattle metagenomes, ranging 0.5-4.0%. We highlight that current alignment- and assembly-based methods estimating resistome mobility yield contradictory and incomplete results, likely constrained by approach-specific data inputs, and bioinformatic limitations. We discuss recent laboratory and computational advancements that may enhance resistome risk analysis in clinical, regulatory, and commercial applications.

Reconstruction of plasmids by shotgun sequencing from environmental DNA: which bioinformatic workflow?

Plasmids play important roles in microbial evolution and also in the spread of antibiotic resistance. Plasmid sequences are extensively studied from clinical isolates but rarely from the environment with a metagenomic approach focused on the plasmid fraction referred to as the plasmidome. A clear challenge in this context is to define a workflow for discriminating plasmids from chromosomal contaminants existing in the plasmidome. For this purpose, we benchmarked existing tools from assembly to detection of the plasmids by reference-free methods (cBar and PlasFlow) and database-guided approaches. Our simulations took into account short-reads alone or combined with moderate long-reads like those actually generated in environmental genomics experiments. This benchmark allowed us to select the best tools for limiting false-positives associated to plasmid prediction tools and a combination of reference-guided methods based on plasmid and bacterial databases.

PHAST, PHASTER and PHASTEST: Tools for finding prophage in bacterial genomes

PHAST (PHAge Search Tool) and its successor PHASTER (PHAge Search Tool - Enhanced Release) have become two of the most widely used web servers for identifying putative prophages in bacterial genomes. Here we review the main capabilities of these web resources, provide some practical guidance regarding their use and discuss possible future improvements. PHAST, which was first described in 2011, made its debut just as whole bacterial genome sequencing and was becoming inexpensive and relatively routine. PHAST quickly gained popularity among bacterial genome researchers because of its web accessibility, its ease of use along with its enhanced accuracy and rapid processing times. PHASTER, which appeared in 2016, provided a number of much-needed enhancements to the PHAST server, including greater processing speed (to cope with very large submission volumes), increased database sizes, a more modern user interface, improved graphical displays and support for metagenomic submissions. Continuing developments in the field, along with increased interest in automated phage and prophage finding, have already led to several improvements to the PHASTER server and will soon lead to the development of a successor to PHASTER (to be called PHASTEST).

Beyond Just Bacteria: Functional Biomes in the Gut Ecosystem Including Virome, Mycobiome, Archaeome and Helminths

Gut microbiota refers to a complex network of microbes, which exerts a marked influence on the host’s health. It is composed of bacteria, fungi, viruses, and helminths. Bacteria, or collectively, the bacteriome, comprises a significant proportion of the well-characterized microbiome. However, the other communities referred to as ‘dark matter’ of microbiomes such as viruses (virome), fungi (mycobiome), archaea (archaeome), and helminths have not been completely elucidated. Development of new and improved metagenomics methods has allowed the identification of complete genomes from the genetic material in the human gut, opening new perspectives on the understanding of the gut microbiome composition, their importance, and potential clinical applications. Here, we review the recent evidence on the viruses, fungi, archaea, and helminths found in the mammalian gut, detailing their interactions with the resident bacterial microbiota and the host, to explore the potential impact of the microbiome on host’s health. The role of fecal virome transplantations, pre-, pro-, and syn-biotic interventions in modulating the microbiome and their related concerns are also discussed.

Comparative Genomic Analysis Confirms Five Genetic Populations of the Select Agent, Rathayibacter toxicus

Rathayibacter toxicus is a Gram-positive, nematode-vectored bacterium that infects several grass species in the family Poaceae. Unique in its genus, R. toxicus has the smallest genome, possesses a complete CRISPR-Cas system, a vancomycin-resistance cassette, produces tunicamycin, a corynetoxin responsible for livestock deaths in Australia, and is designated a Select Agent in the United States. In-depth, genome-wide analyses performed in this study support the previously designated five genetic populations, with a core genome comprising approximately 80% of the genome for all populations. Results varied as a function of the type of analysis and when using different bioinformatics tools for the same analysis; e.g., some programs failed to identify specific genomic regions that were actually present. The software variance highlights the need to verify bioinformatics results by additional methods; e.g., PCR, mapping genes to genomes, use of multiple algorithms). These analyses suggest the following relationships among populations: RT-IV ↔ RT-I ↔ RT-II ↔ RT-III ↔ RT-V, with RT-IV and RT-V being the most unrelated. This is the most comprehensive analysis of R. toxicus that included populations RT-I and RT-V. Future studies require underrepresented populations and more recent isolates from varied hosts and geographic locations.

Longitudinal survey of microbiome associated with particulate matter in a megacity

BACKGROUND

While the physical and chemical properties of airborne particulate matter (PM) have been extensively studied, their associated microbiome remains largely unexplored. Here, we performed a longitudinal metagenomic survey of 106 samples of airborne PM_2.5 and PM₁₀ in Beijing over a period of 6 months in 2012 and 2013, including those from several historically severe smog events.

RESULTS

We observed that the microbiome composition and functional potential were conserved between PM_2.5 and PM₁₀, although considerable temporal variations existed. Among the airborne microorganisms, Propionibacterium acnes, Escherichia coli, Acinetobacter lwoffii, Lactobacillus amylovorus, and Lactobacillus reuteri dominated, along with several viral species. We further identified an extensive repertoire of genes involved in antibiotic resistance and detoxification, including transporters, transpeptidases, and thioredoxins. Sample stratification based on Air Quality Index (AQI) demonstrated that many microbial species, including those associated with human, dog, and mouse feces, exhibit AQI-dependent incidence dynamics. The phylogenetic and functional diversity of air microbiome is comparable to those of soil and water environments, as its composition likely derives from a wide variety of sources.

CONCLUSIONS

Airborne particulate matter accommodates rich and dynamic microbial communities, including a range of microbial elements that are associated with potential health consequences.

Microbial and Viral Communities and Their Antibiotic Resistance Genes Throughout a Hospital Wastewater Treatment System

Antibiotic resistance poses a serious threat to global public health, and antibiotic resistance determinants can enter natural aquatic systems through discharge of wastewater effluents. Hospital wastewater in particular is expected to contain high abundances of antibiotic resistance genes (ARGs) compared to municipal wastewater because it contains human enteric bacteria that may include antibiotic-resistant organisms originating from hospital patients, and can also have high concentrations of antibiotics and antimicrobials relative to municipal wastewater. Viruses also play an important role in wastewater treatment systems since they can influence the bacterial community composition through killing bacteria, facilitating transduction of genetic material between organisms, and modifying the chromosomal content of bacteria as prophages. However, little is known about the fate and connections between ARGs, viruses, and their associated bacteria in hospital wastewater systems. To address this knowledge gap, we characterized the composition and persistence of ARGs, dsDNA viruses, and bacteria from influent to effluent in a pilot-scale hospital wastewater treatment system in Israel using shotgun metagenomics. Results showed that ARGs, including genes conferring resistance to antibiotics of high clinical relevance, were detected in all sampling locations throughout the pilot-scale system, with only 16% overall depletion of ARGs per genome equivalent between influent and effluent. The most common classes of ARGs detected throughout the system conferred resistance to aminoglycoside, cephalosporin, macrolide, penam, and tetracycline antibiotics. A greater proportion of total ARGs were associated with plasmid-associated genes in effluent compared to in influent. No strong associations between viral sequences and ARGs were identified in viral metagenomes from the system, suggesting that phage may not be a significant vector for ARG transfer in this system. The majority of viruses in the pilot-scale system belonged to the families Myoviridae, Podoviridae, and Siphoviridae. Gammaproteobacteria was the dominant class of bacteria harboring ARGs and the most common putative viral host in all samples, followed by Bacilli and Betaproteobacteria. In the total bacterial community, the dominant class was Betaproteobacteria for each sample. Overall, we found that a variety of different types of ARGs and viruses were persistent throughout this hospital wastewater treatment system, which can be released to the environment through effluent discharge.

Plasmid Reconstruction from Next-Gen Data: A Detailed Protocol for the Use of PLACNETw for the Reconstruction of Plasmids from WGS Datasets

Mobile Genetic Elements (MGE) play essential roles in adaptive bacterial evolution, facilitating genetic exchange for extrachromosomal DNA, especially antibiotic resistance genes and virulence factors. For this reason, high-throughput next-generation sequencing of bacteria is of great relevance, especially for clinical pathogenic bacteria. Accurate identification of MGE from whole-genome sequencing (WGS) datasets is one of the major challenges, still hindered by methodological limitations and high sequencing costs.This chapter encompasses the protocol used for plasmid reconstruction by applying the PLACNETw methodology, from raw reads to assembled plasmids and chromosome. PLACNETw is a graphical user-friendly interface to visualize and reconstruct MGE from short-read WGS datasets. No bioinformatic background or sophisticated computational resources are required and high precision and sensitivity are achieved.

Comparative Genomic Analysis of Rhodococcus equi: An Insight into Genomic Diversity and Genome Evolution

Rhodococcus equi, a member of the Rhodococcus genus, is a gram-positive pathogenic bacterium. Rhodococcus possesses an open pan-genome that constitutes the basis of its high genomic diversity and allows for adaptation to specific niche conditions and the changing host environments. Our analysis further showed that the core genome of R. equi contributes to the pathogenicity and niche adaptation of R. equi. Comparative genomic analysis revealed that the genomes of R. equi shared identical collinearity relationship, and heterogeneity was mainly acquired by means of genomic islands and prophages. Moreover, genomic islands in R. equi were always involved in virulence, resistance, or niche adaptation and possibly working with prophages to cause the majority of genome expansion. These findings provide an insight into the genomic diversity, evolution, and structural variation of R. equi and a valuable resource for functional genomic studies.

Comprehensive analysis of chromosomal mobile genetic elements in the gut microbiome reveals phylum-level niche-adaptive gene pools

Mobile genetic elements (MGEs) drive extensive horizontal transfer in the gut microbiome. This transfer could benefit human health by conferring new metabolic capabilities to commensal microbes, or it could threaten human health by spreading antibiotic resistance genes to pathogens. Despite their biological importance and medical relevance, MGEs from the gut microbiome have not been systematically characterized. Here, we present a comprehensive analysis of chromosomal MGEs in the gut microbiome using a method that enables the identification of the mobilizable unit of MGEs. We curated a database of 5,219 putative MGEs encompassing seven MGE classes called ImmeDB. We observed that many MGEs carry genes that could confer an adaptive advantage to the gut environment including gene families involved in antibiotic resistance, bile salt detoxification, mucus degradation, capsular polysaccharide biosynthesis, polysaccharide utilization, and sporulation. We find that antibiotic resistance genes are more likely to be spread by conjugation via integrative conjugative elements or integrative mobilizable elements than transduction via prophages. Horizontal transfer of MGEs is extensive within phyla but rare across phyla, supporting phylum level niche-adaptive gene pools in the gut microbiome. ImmeDB will be a valuable resource for future studies on the gut microbiome and MGE communities.

Quantifying the forces that maintain prophages in bacterial genomes

Genome sequencing has revealed that prophages, viral sequences integrated in a bacterial chromosome, are abundant, accounting for as much as 20% of the bacterial genome. These sequences can confer fitness benefits to the bacterial host, but may also instigate cell death through induction. Several recent investigations have revealed that the distribution of prophage lengths is bimodal, with a clear distinction between small and large prophages. Here we develop a mathematical model of the evolutionary forces affecting the prophage size distribution, and fit this model to three recent data sets. This approach offers quantitative estimates for the relative rates of lysogeny, induction, mutational degradation and selection acting on a wide class of prophage sequences. The model predicts that large prophages are predominantly maintained by the introduction of new prophage sequences through lysogeny, whereas shorter prophages can be enriched when they no longer encode the genes necessary for induction, but still offer selective benefits to their hosts.

Comprehensive Analysis Reveals the Evolution and Pathogenicity of Aeromonas, Viewed from Both Single Isolated Species and Microbial Communities

The genus Aeromonas is a common gastrointestinal pathogen associated with human and animal infections. Due to the high level of cross-species similarity, their evolutionary dynamics and genetic diversity are still fragmented. Hereby, we investigated the pan-genomes of 29 Aeromonas species, as well as Aeromonas species in microbial communities, to clarify their evolutionary dynamics and genetic diversity, with special focus on virulence factors and horizontal gene transfer events. Our study revealed an open pan-genome of Aeromonas containing 10,144 gene families. These Aeromonas species exhibited different functional constraints, with the single-copy core genes and most accessory genes experiencing purifying selection. The significant congruence between core genome and pan-genome trees revealed that core genes mainly affected evolutionary divergences of Aeromonas species. Gene gains and losses revealed a high level of genome plasticity, exhibited by hundreds of gene expansions and contractions, horizontally transferred genes, and mobile genetic elements. The selective constraints shaped virulence gene pools of these Aeromonas strains, where genes encoding hemolysin were ubiquitous. Of these strains, Aeromonas aquatica MX16A seemed to be more resistant, as it harbored most resistance genes. Finally, the virulence factors of Aeromonas in microbial communities were quite dynamic in response to environment changes. For example, the virulence diversity of Aeromonas in microbial communities could reach levels that match some of the most virulent Aeromonas species (such as A. hydrophila) in penetrated-air and modified-air packaging. Our work shed some light onto genetic diversity, evolutionary history, and functional features of Aeromonas, which could facilitate the detection and prevention of infections.IMPORTANCE Aeromonas has long been known as a gastrointestinal pathogen, yet it has many species whose evolutionary dynamics and genetic diversity had been unclear until now. We have conducted pan-genome analysis for 29 Aeromonas species and revealed a high level of genome plasticity exhibited by hundreds of gene expansions and contractions, horizontally transferred genes, and mobile genetic elements. These species also contained many virulence factors both identified from single isolated species and microbial community. This pan-genome study could elevate the level for detection and prevention of Aeromonas infections.

Comparative Study of Gut Microbiota in Wild and Captive Giant Pandas (Ailuropoda melanoleuca)

Captive breeding has been used as an effective approach to protecting endangered animals but its effect on the gut microbiome and the conservation status of these species is largely unknown. The giant panda is a flagship species for the conservation of wildlife. With integrated efforts including captive breeding, this species has been recently upgraded from "endangered" to "vulnerable" (IUCN 2016). Since a large proportion (21.8%) of their global population is still captive, it is critical to understand how captivity changes the gut microbiome of these pandas and how such alterations to the microbiome might affect their future fitness and potential impact on the ecosystem after release into the wild. Here, we use 16S rRNA (ribosomal RNA) marker gene sequencing and shotgun metagenomics sequencing to demonstrate that the fecal microbiomes differ substantially between wild and captive giant pandas. Fecal microbiome diversity was significantly lower in captive pandas, as was the diversity of functional genes. Additionally, captive pandas have reduced functional potential for cellulose degradation but enriched metabolic pathways for starch metabolism, indicating that they may not adapt to a wild diet after being released into the wild since a major component of their diet in the wild will be bamboo. Most significantly, we observed a significantly higher level of amylase activity but a lower level of cellulase activity in captive giant panda feces than those of wild giant pandas, shown by an in vitro experimental assay. Furthermore, antibiotic resistance genes and virulence factors, as well as heavy metal tolerance genes were enriched in the microbiomes of captive pandas, which raises a great concern of spreading these genes to other wild animals and ecosystems when they are released into a wild environment. Our results clearly show that captivity has altered the giant panda microbiome, which could have unintended negative consequences on their adaptability and the ecosystem during the reintroduction of giant pandas into the wild.