Long-read metagenomic exploration of extrachromosomal mobile genetic elements in the human gut

The Third-Generation Sequencing Challenge: Novel Insights for the Omic Sciences

The understanding of the human genome has been greatly improved by the advent of next-generation sequencing technologies (NGS). Despite the undeniable advantages responsible for their widespread diffusion, these methods have some constraints, mainly related to short read length and the need for PCR amplification. As a consequence, long-read sequencers, called third-generation sequencing (TGS), have been developed, promising to overcome NGS. Starting from the first prototype, TGS has progressively ameliorated its chemistries by improving both read length and base-calling accuracy, as well as simultaneously reducing the costs/base. Based on these premises, TGS is showing its potential in many fields, including the analysis of difficult-to-sequence genomic regions, structural variations detection, RNA expression profiling, DNA methylation study, and metagenomic analyses. Protocol standardization and the development of easy-to-use pipelines for data analysis will enhance TGS use, also opening the way for their routine applications in diagnostic contexts.

Current examining methods and mathematical models of horizontal transfer of antibiotic resistance genes in the environment

The increasing prevalence of antibiotic resistance genes (ARGs) in the environment has garnered significant attention due to their health risk to human beings. Horizontal gene transfer (HGT) is considered as an important way for ARG dissemination. There are four general routes of HGT, including conjugation, transformation, transduction and vesiduction. Selection of appropriate examining methods is crucial for comprehensively understanding characteristics and mechanisms of different HGT ways. Moreover, combined with the results obtained from different experimental methods, mathematical models could be established and serve as a powerful tool for predicting ARG transfer dynamics and frequencies. However, current reviews of HGT for ARG spread mainly focus on its influencing factors and mechanisms, overlooking the important roles of examining methods and models. This review, therefore, delineated four pathways of HGT, summarized the strengths and limitations of current examining methods, and provided a comprehensive summing-up of mathematical models pertaining to three main HGT ways of conjugation, transformation and transduction. Finally, deficiencies in current studies were discussed, and proposed the future perspectives to better understand and assess the risks of ARG dissemination through HGT.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Detecting horizontal gene transfer with metagenomics co-barcoding sequencing

Horizontal gene transfer (HGT) is the process through which genetic information is transferred between different genomes and that played a crucial role in bacterial evolution. HGT can enable bacteria to rapidly acquire antibiotic resistance and bacteria that have acquired resistance is spreading within the microbiome. Conventional methods of characterizing HGT patterns include short-read metagenomic sequencing (short-reads mNGS), long-read sequencing, and single-cell sequencing. These approaches present several limitations, such as short-read fragments, high amounts of input DNA, and sequencing costs, respectively. Here, we attempt to circumvent present limitations to detect HGT by developing a metagenomics co-barcode sequencing workflow (MECOS) and applying it to the human and mouse gut microbiomes. In addition to that, we have over 10-fold increased contig length compared to short-reads mNGS; we also obtained exceeding 30 million paired reads with co-barcode information. Applying the novel bioinformatic pipeline, we integrated this co-barcoding information and the context information from long reads, and observed over 50-fold HGT events after we corrected the potential wrong HGT events. Specifically, we detected approximately 3,000 HGT blocks in individual samples, encompassing ~6,000 genes and ~100 taxonomic groups, including loci conferring tetracycline resistance through ribosomal protection. MECOS provides a valuable tool for investigating HGT and advance our understanding on the evolution of natural microbial communities within hosts.IMPORTANCEIn this study, to better identify horizontal gene transfer (HGT) in individual samples, we introduce a new co-barcoding sequencing system called metagenomics co-barcoding sequencing (MECOS), which has three significant improvements: (i) long DNA fragment extraction, (ii) a special transposome insertion, (iii) hybridization of DNA to barcode beads, and (4) an integrated bioinformatic pipeline. Using our approach, we have over 10-fold increased contig length compared to short-reads mNGS, and observed over 50-fold HGT events after we corrected the potential wrong HGT events. Our results indicate the presence of approximately 3,000 HGT blocks, involving roughly 6,000 genes and 100 taxonomic groups in individual samples. Notably, these HGT events are predominantly enriched in genes that confer tetracycline resistance via ribosomal protection. MECOS is a useful tool for investigating HGT and the evolution of natural microbial communities within hosts, thereby advancing our understanding of microbial ecology and evolution.

Exploring the effects of antimicrobial treatment on the gut and oral microbiomes and resistomes from elderly long-term care facility residents via shotgun DNA sequencing

Monitoring antibiotic-resistant bacteria (ARB) and understanding the effects of antimicrobial drugs on the human microbiome and resistome are crucial for public health. However, no study has investigated the association between antimicrobial treatment and the microbiome-resistome relationship in long-term care facilities, where residents act as reservoirs of ARB but are not included in the national surveillance for ARB. We conducted shotgun metagenome sequencing of oral and stool samples from long-term care facility residents and explored the effects of antimicrobial treatment on the human microbiome and resistome using two types of comparisons: cross-sectional comparisons based on antimicrobial treatment history in the past 6 months and within-subject comparisons between stool samples before, during and 2-4 weeks after treatment using a single antimicrobial drug. Cross-sectional analysis revealed two characteristics in the group with a history of antimicrobial treatment: the archaeon Methanobrevibacter was the only taxon that significantly increased in abundance, and the total abundance of antimicrobial resistance genes (ARGs) was also significantly higher. Within-subject comparisons showed that taxonomic diversity did not decrease during treatment, suggesting that the effect of the prescription of a single antimicrobial drug in usual clinical treatment on the gut microbiota is likely to be smaller than previously thought, even among very elderly people. Additional analysis of the detection limit of ARGs revealed that they could not be detected when contig coverage was <2.0. This study is the first to report the effects of usual antimicrobial treatments on the microbiome and resistome of long-term care facility residents.

Insights into the circular: The cryptic plasmidome and its derived antibiotic resistome in the urban water systems

Plasmids have been a concern in the dissemination and evolution of antibiotic resistance in the environment. In this study, we investigated the total pool of plasmids (plasmidome) and its derived antibiotic resistance genes (ARGs) in different compartments of urban water systems (UWSs) in three European countries representing different antibiotic usage regimes. We applied a direct plasmidome approach using wet-lab methods to enrich circular DNA in the samples, followed by shotgun sequencing and in silico contig circularisation. We identified 9538 novel sequences in a total of 10,942 recovered circular plasmids. Of these, 66 were identified as conjugative, 1896 mobilisable and 8970 non-mobilisable plasmids. The UWSs' plasmidome was dominated by small plasmids (≤10 Kbp) representing a broad diversity of mobility (MOB) types and incompatibility (Inc) groups. A shared collection of plasmids from different countries was detected in all treatment compartments, and plasmids could be source-tracked in the UWSs. More than half of the ARGs-encoding plasmids carried mobility genes for mobilisation/conjugation. The richness and abundance of ARGs-encoding plasmids generally decreased with the flow, while we observed that non-mobilisable ARGs-harbouring plasmids maintained their abundance in the Spanish wastewater treatment plant. Overall, our work unravels that the UWS plasmidome is dominated by cryptic (i.e., non-mobilisable, non-typeable and previously unknown) plasmids. Considering that some of these plasmids carried ARGs, were prevalent across three countries and could persist throughout the UWSs compartments, these results should alarm and call for attention.

Pan-genome analysis of six Paracoccus type strain genomes reveal lifestyle traits

The genus Paracoccus capable of inhabiting a variety of different ecological niches both, marine and terrestrial, is globally distributed. In addition, Paracoccus is taxonomically, metabolically and regarding lifestyle highly diverse. Until now, little is known on how Paracoccus can adapt to such a range of different ecological niches and lifestyles. In the present study, the genus Paracoccus was phylogenomically analyzed (n = 160) and revisited, allowing species level classification of 16 so far unclassified Paracoccus sp. strains and detection of five misclassifications. Moreover, we performed pan-genome analysis of Paracoccus-type strains, isolated from a variety of ecological niches, including different soils, tidal flat sediment, host association such as the bluespotted cornetfish, Bugula plumosa, and the reef-building coral Stylophora pistillata to elucidate either i) the importance of lifestyle and adaptation potential, and ii) the role of the genomic equipment and niche adaptation potential. Six complete genomes were de novo hybrid assembled using a combination of short and long-read technologies. These Paracoccus genomes increase the number of completely closed high-quality genomes of type strains from 15 to 21. Pan-genome analysis revealed an open pan-genome composed of 13,819 genes with a minimal chromosomal core (8.84%) highlighting the genomic adaptation potential and the huge impact of extra-chromosomal elements. All genomes are shaped by the acquisition of various mobile genetic elements including genomic islands, prophages, transposases, and insertion sequences emphasizing their genomic plasticity. In terms of lifestyle, each mobile genetic elements should be evaluated separately with respect to the ecological context. Free-living genomes, in contrast to host-associated, tend to comprise (1) larger genomes, or the highest number of extra-chromosomal elements, (2) higher number of genomic islands and insertion sequence elements, and (3) a lower number of intact prophage regions. Regarding lifestyle adaptations, free-living genomes share genes linked to genetic exchange via T4SS, especially relevant for Paracoccus, known for their numerous extrachromosomal elements, enabling adaptation to dynamic environments. Conversely, host-associated genomes feature diverse genes involved in molecule transport, cell wall modification, attachment, stress protection, DNA repair, carbon, and nitrogen metabolism. Due to the vast number of adaptive genes, Paracoccus can quickly adapt to changing environmental conditions.

De Novo Synthesis of a Conjugative System from Human Gut Metagenomic Data for Targeted Delivery of Cas9 Antimicrobials

Metagenomic sequences represent an untapped source of genetic novelty, particularly for conjugative systems that could be used for plasmid-based delivery of Cas9-derived antimicrobial agents. However, unlocking the functional potential of conjugative systems purely from metagenomic sequences requires the identification of suitable candidate systems as starting scaffolds for de novo DNA synthesis. Here, we developed a bioinformatics approach that searches through the metagenomic "trash bin" for genes associated with conjugative systems present on contigs that are typically excluded from common metagenomic analysis pipelines. Using a human metagenomic gut data set representing 2805 taxonomically distinct units, we identified 1598 contigs containing conjugation genes with a differential distribution in human cohorts. We synthesized de novo an entire Citrobacter spp. conjugative system of 54 kb containing at least 47 genes and assembled it into a plasmid, pCitro. We found that pCitro conjugates from Escherichia coli to Citrobacter rodentium with a 30-fold higher frequency than to E. coli, and is compatible with Citrobacter resident plasmids. Mutations in the traV and traY conjugation components of pCitro inhibited conjugation. We showed that pCitro can be repurposed as an antimicrobial delivery agent by programming it with the TevCas9 nuclease and Citrobacter-specific sgRNAs to kill C. rodentium. Our study reveals a trove of uncharacterized conjugative systems in metagenomic data and describes an experimental framework to animate these large genetic systems as novel target-adapted delivery vectors for Cas9-based editing of bacterial genomes.

Long-read assembled metagenomic approaches improve our understanding on metabolic potentials of microbial community in mangrove sediments

BACKGROUND

Mangrove wetlands are coastal ecosystems with important ecological features and provide habitats for diverse microorganisms with key roles in nutrient and biogeochemical cycling. However, the overall metabolic potentials and ecological roles of microbial community in mangrove sediment are remained unanswered. In current study, the microbial and metabolic profiles of prokaryotic and fungal communities in mangrove sediments were investigated using metagenomic analysis based on PacBio single-molecule real time (SMRT) and Illumina sequencing techniques.

RESULTS

Comparing to Illumina short reads, the incorporation of PacBio long reads significantly contributed to more contiguous assemblies, yielded more than doubled high-quality metagenome-assembled genomes (MAGs), and improved the novelty of the MAGs. Further metabolic reconstruction for recovered MAGs showed that prokaryotes potentially played an essential role in carbon cycling in mangrove sediment, displaying versatile metabolic potential for degrading organic carbons, fermentation, autotrophy, and carbon fixation. Mangrove fungi also functioned as a player in carbon cycling, potentially involved in the degradation of various carbohydrate and peptide substrates. Notably, a new candidate bacterial phylum named as Candidatus Cosmopoliota with a ubiquitous distribution is proposed. Genomic analysis revealed that this new phylum is capable of utilizing various types of organic substrates, anaerobic fermentation, and carbon fixation with the Wood-Ljungdahl (WL) pathway and the reverse tricarboxylic acid (rTCA) cycle.

CONCLUSIONS

The study not only highlights the advantages of HiSeq-PacBio Hybrid assembly for a more complete profiling of environmental microbiomes but also expands our understanding of the microbial diversity and potential roles of distinct microbial groups in biogeochemical cycling in mangrove sediment. Video Abstract.

Noise reduction strategies in metagenomic chromosome confirmation capture to link antibiotic resistance genes to microbial hosts

The gut microbiota is a reservoir for antimicrobial resistance genes (ARGs). With current sequencing methods, it is difficult to assign ARGs to their microbial hosts, particularly if these ARGs are located on plasmids. Metagenomic chromosome conformation capture approaches (meta3C and Hi-C) have recently been developed to link bacterial genes to phylogenetic markers, thus potentially allowing the assignment of ARGs to their hosts on a microbiome-wide scale. Here, we generated a meta3C dataset of a human stool sample and used previously published meta3C and Hi-C datasets to investigate bacterial hosts of ARGs in the human gut microbiome. Sequence reads mapping to repetitive elements were found to cause problematic noise in, and may importantly skew interpretation of, meta3C and Hi-C data. We provide a strategy to improve the signal-to-noise ratio by discarding reads that map to insertion sequence elements and to the end of contigs. We also show the importance of using spike-in controls to quantify whether the cross-linking step in meta3C and Hi-C protocols has been successful. After filtering to remove artefactual links, 87 ARGs were assigned to their bacterial hosts across all datasets, including 27 ARGs in the meta3C dataset we generated. We show that commensal gut bacteria are an important reservoir for ARGs, with genes coding for aminoglycoside and tetracycline resistance being widespread in anaerobic commensals of the human gut.

Metagenomic analysis reveals gut plasmids as diagnosis markers for colorectal cancer

Background

Colorectal cancer (CRC) is linked to distinct gut microbiome patterns. The efficacy of gut bacteria as diagnostic biomarkers for CRC has been confirmed. Despite the potential to influence microbiome physiology and evolution, the set of plasmids in the gut microbiome remains understudied.

Methods

We investigated the essential features of gut plasmid using metagenomic data of 1,242 samples from eight distinct geographic cohorts. We identified 198 plasmid-related sequences that differed in abundance between CRC patients and controls and screened 21 markers for the CRC diagnosis model. We utilize these plasmid markers combined with bacteria to construct a random forest classifier model to diagnose CRC.

Results

The plasmid markers were able to distinguish between the CRC patients and controls [mean area under the receiver operating characteristic curve (AUC = 0.70)] and maintained accuracy in two independent cohorts. In comparison to the bacteria-only model, the performance of the composite panel created by combining plasmid and bacteria features was significantly improved in all training cohorts (mean AUC_composite = 0.804 and mean AUC_bacteria = 0.787) and maintained high accuracy in all independent cohorts (mean AUC_composite = 0.839 and mean AUC_bacteria = 0.821). In comparison to controls, we found that the bacteria-plasmid correlation strength was weaker in CRC patients. Additionally, the KEGG orthology (KO) genes in plasmids that are independent of bacteria or plasmids significantly correlated with CRC.

Conclusion

We identified plasmid features associated with CRC and showed how plasmid and bacterial markers could be combined to further enhance CRC diagnosis accuracy.

A survey on computational strategies for genome-resolved gut metagenomics

Recovering high-quality metagenome-assembled genomes (HQ-MAGs) is critical for exploring microbial compositions and microbe-phenotype associations. However, multiple sequencing platforms and computational tools for this purpose may confuse researchers and thus call for extensive evaluation. Here, we systematically evaluated a total of 40 combinations of popular computational tools and sequencing platforms (i.e. strategies), involving eight assemblers, eight metagenomic binners and four sequencing technologies, including short-, long-read and metaHiC sequencing. We identified the best tools for the individual tasks (e.g. the assembly and binning) and combinations (e.g. generating more HQ-MAGs) depending on the availability of the sequencing data. We found that the combination of the hybrid assemblies and metaHiC-based binning performed best, followed by the hybrid and long-read assemblies. More importantly, both long-read and metaHiC sequencings link more mobile elements and antibiotic resistance genes to bacterial hosts and improve the quality of public human gut reference genomes with 32% (34/105) HQ-MAGs that were either of better quality than those in the Unified Human Gastrointestinal Genome catalog version 2 or novel.

Genomic and physiological characterization of Novosphingobium terrae sp. nov., an alphaproteobacterium isolated from Cerrado soil containing a mega-sized chromid

A novel bacterial strain, designated GeG2T, was isolated from soils of the native Cerrado, a highly biodiverse savanna-like Brazilian biome. 16S rRNA gene analysis of GeG2T revealed high sequence identity (100%) to the alphaproteobacterium Novosphingobium rosa; however, comparisons with N. rosa DSM 7285T showed several distinctive features, prompting a full characterization of the new strain in terms of physiology, morphology, and, ultimately, its genome. GeG2T cells were Gram-stain-negative bacilli, facultatively anaerobic, motile, positive for catalase and oxidase activities, and starch hydrolysis. Strain GeG2T presented planktonic-sessile dimorphism and cell aggregates surrounded by extracellular matrix and nanometric spherical structures were observed, suggesting the production of exopolysaccharides (EPS) and outer membrane vesicles (OMVs). Despite high 16S rDNA identity, strain GeG2T showed 90.38% average nucleotide identity and 42.60% digital DNA-DNA hybridization identity with N. rosa, below species threshold. Whole-genome assembly revealed four circular replicons: a 4.1 Mb chromosome, a 2.7 Mb extrachromosomal megareplicon, and two plasmids (212.7 and 68.6 kb). The megareplicon contains a few core genes and plasmid-type replication/maintenance systems, consistent with its classification as a chromid. Genome annotation shows a vast repertoire of carbohydrate-active enzymes and genes involved in the degradation of aromatic compounds, highlighting the biotechnological potential of the new isolate. Chemotaxonomic features, including polar lipid and fatty acid profiles, as well as physiological, molecular, and whole-genome comparisons showed significant differences between strain GeG2T and N. rosa, indicating that it represents a novel species, for which the name Novosphingobium terrae is proposed. The type strain is GeG2T (= CBMAI 2313T = CBAS 753 T).

Managing the Transition to Widespread Metagenomic Monitoring: Policy Considerations for Future Biosurveillance

The technological possibilities and future public health importance of metagenomic sequencing have received extensive attention, but there has been little discussion about the policy and regulatory issues that need to be addressed if metagenomic sequencing is adopted as a key technology for biosurveillance. In this article, we introduce metagenomic monitoring as a possible path to eventually replacing current infectious disease monitoring models. Many key enablers are technological, whereas others are not. We therefore highlight key policy challenges and implementation questions that need to be addressed for "widespread metagenomic monitoring" to be possible. Policymakers must address pitfalls like fragmentation of the technological base, private capture of benefits, privacy concerns, the usefulness of the system during nonpandemic times, and how the future systems will enable better response. If these challenges are addressed, the technological and public health promise of metagenomic sequencing can be realized.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Target-enriched long-read sequencing (TELSeq) contextualizes antimicrobial resistance genes in metagenomes

BACKGROUND

Metagenomic data can be used to profile high-importance genes within microbiomes. However, current metagenomic workflows produce data that suffer from low sensitivity and an inability to accurately reconstruct partial or full genomes, particularly those in low abundance. These limitations preclude colocalization analysis, i.e., characterizing the genomic context of genes and functions within a metagenomic sample. Genomic context is especially crucial for functions associated with horizontal gene transfer (HGT) via mobile genetic elements (MGEs), for example antimicrobial resistance (AMR). To overcome this current limitation of metagenomics, we present a method for comprehensive and accurate reconstruction of antimicrobial resistance genes (ARGs) and MGEs from metagenomic DNA, termed target-enriched long-read sequencing (TELSeq).

RESULTS

Using technical replicates of diverse sample types, we compared TELSeq performance to that of non-enriched PacBio and short-read Illumina sequencing. TELSeq achieved much higher ARG recovery (>1,000-fold) and sensitivity than the other methods across diverse metagenomes, revealing an extensive resistome profile comprising many low-abundance ARGs, including some with public health importance. Using the long reads generated by TELSeq, we identified numerous MGEs and cargo genes flanking the low-abundance ARGs, indicating that these ARGs could be transferred across bacterial taxa via HGT.

CONCLUSIONS

TELSeq can provide a nuanced view of the genomic context of microbial resistomes and thus has wide-ranging applications in public, animal, and human health, as well as environmental surveillance and monitoring of AMR. Thus, this technique represents a fundamental advancement for microbiome research and application. Video abstract.

Antibiotic Resistance Genes in River Biofilms: A Metagenomic Approach toward the Identification of Sources and Candidate Hosts

Treated wastewater is a major pathway by which antibiotic resistance genes (ARG) enter aquatic ecosystems. However, knowledge gaps remain concerning the dissemination of specific ARG and their association with bacterial hosts. Here, we employed shotgun metagenomics to track ARG and taxonomic markers in river biofilms along a gradient of fecal pollution depicted by crAssphage signatures. We found strong evidence for an impact of wastewater effluents on both community composition and resistomes. In the light of such simultaneity, we employed a model comparison technique to identify ARG-host relationships from nonassembled metagenomic DNA. Hereby, a major cause of spurious associations otherwise encountered in correlation-based ARG-host analyses was suppressed. For several families of ARG, namely those conferring resistance to beta-lactams, particular bacterial orders were identified as candidate hosts. The found associations of blaFOX and cphA with Aeromonadales or blaPER with Chromatiales support the outcome of independent evolutionary analyses and thus confirm the potential of the methodology. For other ARG families including blaIMP or tet, clusters of bacterial orders were identified which potentially harbor a major proportion of host species. For yet other ARG, like, for example, ant or erm, no particular host candidates were identifiable, indicating their spread across various taxonomic groups.

3CAC: improving the classification of phages and plasmids in metagenomic assemblies using assembly graphs

MOTIVATION

Bacteriophages and plasmids usually coexist with their host bacteria in microbial communities and play important roles in microbial evolution. Accurately identifying sequence contigs as phages, plasmids and bacterial chromosomes in mixed metagenomic assemblies is critical for further unraveling their functions. Many classification tools have been developed for identifying either phages or plasmids in metagenomic assemblies. However, only two classifiers, PPR-Meta and viralVerify, were proposed to simultaneously identify phages and plasmids in mixed metagenomic assemblies. Due to the very high fraction of chromosome contigs in the assemblies, both tools achieve high precision in the classification of chromosomes but perform poorly in classifying phages and plasmids. Short contigs in these assemblies are often wrongly classified or classified as uncertain.

RESULTS

Here we present 3CAC, a new three-class classifier that improves the precision of phage and plasmid classification. 3CAC starts with an initial three-class classification generated by existing classifiers and improves the classification of short contigs and contigs with low confidence classification by using proximity in the assembly graph. Evaluation on simulated metagenomes and on real human gut microbiome samples showed that 3CAC outperformed PPR-Meta and viralVerify in both precision and recall, and increased F1-score by 10-60 percentage points.

AVAILABILITY AND IMPLEMENTATION

The 3CAC software is available on https://github.com/Shamir-Lab/3CAC.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Extensive gut virome variation and its associations with host and environmental factors in a population-level cohort

Indigenous bacteriophage communities (virome) in the human gut have a huge impact on the structure and function of gut bacterial communities (bacteriome), but virome variation at a population scale is not fully investigated yet. Here, we analyse the gut dsDNA virome in the Japanese 4D cohort of 4198 deeply phenotyped individuals. By assembling metagenomic reads, we discover thousands of high-quality phage genomes including previously uncharacterised phage clades with different bacterial hosts than known major ones. The distribution of host bacteria is a strong determinant for the distribution of phages in the gut, and virome diversity is highly correlated with anti-viral defence mechanisms of the bacteriome, such as CRISPR-Cas and restriction-modification systems. We identify 97 various intrinsic/extrinsic factors that significantly affect the virome structure, including age, sex, lifestyle, and diet, most of which showed consistent associations with both phages and their predicted bacterial hosts. Among the metadata categories, disease and medication have the strongest effects on the virome structure. Overall, these results present a basis to understand the symbiotic communities of bacteria and their viruses in the human gut, which will facilitate the medical and industrial applications of indigenous viruses.

Here, Nishijima et al. perform a large-scale analysis of the human gut virome in the Japanese 4D cohort of 4198 deeply phenotyped individuals, revealing thousands of bacteriophage genomes, virus-bacteria interactions, and describing associations with various host and environmental factors.

Metagenomic assembly reveals hosts and mobility of common antibiotic resistome in animal manure and commercial compost

BACKGROUND

Antibiotics and antibiotic resistance genes (ARGs) used in intensive animal farming threaten human health worldwide; however, the common resistome, ARG mobility, and ARG host composition in different animal manures and mixed manure composts remain unclear. In the present study, metagenomic assembly and cross-sample mapping were used to comprehensively decipher the common resistome and its potential mobility and hosts in animal manure and composts.

RESULTS

In total, 201 ARGs were shared among different animal (layer, broiler, swine, beef cow, and dairy cow) manures and accounted for 86-99% of total relative abundance of ARGs. Except for multidrug, sulfonamide, and trimethoprim resistance genes, the relative abundance of most ARGs in composts was significantly lower than that in animal manure. Procrustes analysis indicated that antibiotic residues positively correlated with ARG composition in manure but not in composts. More than 75% ARG subtypes were shared between plasmids and chromosomes in our samples. Transposases could play a pivotal role in mediating the transfer of ARGs between different phyla in animal manure and composting. Cross-sample mapping to contigs carrying ARGs showed that the hosts of common resistome in manure had preference on animal species, and the dominant genus of ARG host shifted from Enterococcus in manure to Pseudomonas in composts. The broad host range and linking with diverse mobile genetic elements (MGEs) were two key factors for ARGs, such as sul1 and aadA, which could survive during composting. The multidrug resistance genes represented the dominant ARGs in pathogenic antibiotic-resistant bacteria in manure but could be effectively controlled by composting.

CONCLUSIONS

Our experiments revealed the common resistome in animal manure, classified and relative quantified the ARG hosts, and assessed the mobility of ARGs. Composting can mitigate ARGs in animal manure by altering the bacterial hosts; however, persistent ARGs can escape from the removal because of diverse host range and MGEs. Our findings provide an overall background for source tracking, risk assessment, and control of livestock ARGs.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Epidemiology of Plasmids in Escherichia coli and Klebsiella pneumoniae with Acquired Extended Spectrum Beta-Lactamase Genes Isolated from Chronic Wounds in Ghana

Little information is available on the local epidemiology of mobile genetic elements such as plasmids harboring acquired beta-lactamase genes in Western African Ghana. In the present study, we screened for plasmids in three Escherichia coli and four Klebsiella pneumoniae isolates expressing extended spectrum beta-lactamases (ESBL) mediated by the bla_CTX-M-15 gene from chronically infected wounds of Ghanaian patients. Bacterial isolates were subjected to combined short-read and long-read sequencing to obtain the sequences of their respective plasmids. In the bla_CTX-M-15-gene-carrying plasmids of the four ESBL-positive K. pneumoniae isolates, IncFIB/IncFII (n = 3) and FIA (n = 1) sequences were detected, while in the bla_CTX-M-15-gene-carrying plasmids of the three ESBL-positive E. coli isolates, IncFIA/IncFIB (n = 2) and IncFIB (n = 1) sequences were found. The three IncFIB/IncFII sequence-containing plasmids were almost identical to a K. pneumoniae plasmid reported from France. They belonged to the clonal lineages ST17, ST36 and ST39 of K. pneumoniae, suggesting transversal spread of this obviously evolutionary successful plasmid in Ghana. Other resistance gene-encoding plasmids observed in the assessed Enterobacterales harbored IncFIA/IncR and IncFII sequences. International spread was confirmed by the high genetic similarity to resistance-mediating plasmids published from Asia, Australia, Europe and Northern America, including a bla_CTX-M-15-gene-carrying plasmid isolated from a wild bird in Germany. In conclusion, the study contributed to the scarcely available information on the epidemiology of third-generation cephalosporine resistance-mediating plasmids in Ghana. Furthermore, the global spread of resistance-mediating plasmids provided hints on the evolutionary success of individual resistance-harboring plasmids by transversal spread among K. pneumoniae lineages in Ghana.

Precise genotyping of circular mobile elements from metagenomic data uncovers human-associated plasmids with recent common ancestors

Mobile genetic elements with circular genomes play a key role in the evolution of microbial communities. Their circular genomes correspond to circular walks in metagenome graphs, and yet, assemblies derived from natural microbial communities produce graphs riddled with spurious cycles, complicating the accurate reconstruction of circular genomes. We present DomCycle, an algorithm that reconstructs likely circular genomes based on the identification of so-called 'dominant' graph cycles. In the implementation we leverage paired reads to bridge assembly gaps and scrutinize cycles through a nucleotide-level analysis, making the approach robust to misassembly artifacts. We validated the approach using simulated and real sequencing data. Application of DomCycle to 32 publicly available DNA shotgun sequence data sets from diverse natural environments led to the reconstruction of hundreds of circular mobile genomes. Clustering revealed 20 highly prevalent and cryptic plasmids that have clonal population structures with recent common ancestors. This method facilitates the study of microbial communities that evolve through horizontal gene transfer.

Novel canine high-quality metagenome-assembled genomes, prophages and host-associated plasmids provided by long-read metagenomics together with Hi-C proximity ligation

The human gut microbiome has been extensively studied, yet the canine gut microbiome is still largely unknown. The availability of high-quality genomes is essential in the fields of veterinary medicine and nutrition to unravel the biological role of key microbial members in the canine gut environment. Our aim was to evaluate nanopore long-read metagenomics and Hi-C (high-throughput chromosome conformation capture) proximity ligation to provide high-quality metagenome-assembled genomes (HQ MAGs) of the canine gut environment. By combining nanopore long-read metagenomics and Hi-C proximity ligation, we retrieved 27 HQ MAGs and 7 medium-quality MAGs of a faecal sample of a healthy dog. Canine MAGs (CanMAGs) improved genome contiguity of representatives from the animal and human MAG catalogues – short-read MAGs from public datasets – for the species they represented: they were more contiguous with complete ribosomal operons and at least 18 canonical tRNAs. Both canine-specific bacterial species and gut generalists inhabit the dog’s gastrointestinal environment. Most of them belonged to Firmicutes, followed by Bacteroidota and Proteobacteria. We also assembled one Actinobacteriota and one Fusobacteriota MAG. CanMAGs harboured antimicrobial-resistance genes (ARGs) and prophages and were linked to plasmids. ARGs conferring resistance to tetracycline were most predominant within CanMAGs, followed by lincosamide and macrolide ones. At the functional level, carbohydrate transport and metabolism was the most variable within the CanMAGs, and mobilome function was abundant in some MAGs. Specifically, we assigned the mobilome functions and the associated mobile genetic elements to the bacterial host. The CanMAGs harboured 50 bacteriophages, providing novel bacterial-host information for eight viral clusters, and Hi-C proximity ligation data linked the six potential plasmids to their bacterial host. Long-read metagenomics and Hi-C proximity ligation are likely to become a comprehensive approach to HQ MAG discovery and assignment of extra-chromosomal elements to their bacterial host. This will provide essential information for studying the canine gut microbiome in veterinary medicine and animal nutrition.

viralFlye: assembling viruses and identifying their hosts from long-read metagenomics data

Although the use of long-read sequencing improves the contiguity of assembled viral genomes compared to short-read methods, assembling complex viral communities remains an open problem. We describe the viralFlye tool for identification and analysis of metagenome-assembled viruses in long-read assemblies. We show it significantly improves viral assemblies and demonstrate that long-reads result in a much larger array of predicted virus-host associations as compared to short-read assemblies. We demonstrate that the identification of novel CRISPR arrays in bacterial genomes from a newly assembled metagenomic sample provides information for predicting novel hosts for novel viruses.

Small investments with big returns: environmental genomic bioprospecting of microbial life

Microorganisms and their natural products are major drivers of ecological processes and industrial applications. Microbial bioprospecting has been critical for the advancement in various fields such as pharmaceuticals, sustainable industries, food security and bioremediation. Next generation sequencing has been paramount in the exploration of diverse environmental microbiomes. It presents a culture-independent approach to investigating hitherto uncultured taxa, resulting in the creation of massive sequence databases, which are available in the public domain. Genome mining searches available (meta)genomic data for target biosynthetic genes, and combined with the large-scale public data, this in-silico bioprospecting method presents an efficient and extensive way to uncover microbial bioproducts. Bioinformatic tools have progressed to a stage where we can recover genomes from the environment; these metagenome-assembled genomes present a way to understand the metabolic capacity of microorganisms in a physiological and ecological context. Environmental sampling been extensive across various ecological settings, including microbiomes with unique physicochemical properties that could influence the discovery of novel functions and metabolic pathways. Although in-silico methods cannot completely substitute in-vitro studies, the contextual information it provides is invaluable for understanding the ecological and taxonomic distribution of microbial genotypes and to form effective strategies for future microbial bioprospecting efforts.

Diverse DNA modification in marine prokaryotic and viral communities

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Hybrid, ultra-deep metagenomic sequencing enables genomic and functional characterization of low-abundance species in the human gut microbiome

A large number of microbial genomes have already been identified from the human gut microbiome, but the understanding of the role of the low-abundance species at the individual level remains challenging, largely due to the relatively shallow sequencing depth used in most studies. To improve genome assembling performance, a HiSeq-PacBio hybrid, ultra-deep metagenomic sequencing approach was used to reconstruct metagenomic-assembled genomes (MAGs) from 12 fecal samples. Such approach combined third-generation sequencing with ultra-deep second-generation sequencing to improve the sequencing coverage of the low-abundance subpopulation in the gut microbiome. Our study generated a total of 44 megabase-scale scaffolds, achieving four single-scaffolds of complete (circularized, no gaps) MAGs (CMAGs) that were the first circular genomes of their species. Moreover, 475 high-quality MAGs were assembled across all samples. Among them, 234 MAGs were currently uncultured, including 24 MAGs that were not found in any public genome database. Additionally, 287 and 77 MAGs were classified as low-abundance (0.1-1%) and extra-low-abundance (<0.1%) gut species in each individual, respectively. Our results also revealed individual-specific genomic features in the MAG profiles, including microbial genome growth rate, selective pressure, and frequency of chromosomal mobile genetic elements. Finally, thousands of extrachromosomal mobile genetic elements were identified from the metagenomic data, including 5097 bacteriophages and 79 novel plasmid genomes. Overall, our strategy represents an important step toward comprehensive genomic and functional characterization of the human gut microbiome at an individual level.

FLUID study: study protocol for an open-label, single-centre pilot study to investigate the efFect of Lemborexant on sleep management in Japanese sUbjects aged 50 years and older with Insomnia Disorder

INTRODUCTION

Bidirectional associations have been reported between sleep disturbance and both cognitive impairment, including Alzheimer's disease and amyloid beta-peptide (Aβ) accumulation. These relationships can be explained by the glymphatic system, which acts as a garbage drainage system in the brain. As interstitial fluid dynamics are suggested to increase during sleep, clearance of Aβ can be influenced by sleep disturbance or deprivation. We hypothesised that using lemborexant, an orexin receptor antagonist, to improve sleep quality would also improve the function of the glymphatic system. We plan to examine the effect of lemborexant on sleep quality and the glymphatic system among patients with insomnia disorder.

METHODS AND ANALYSIS

This pilot study is designed as an open-label, single-arm, single-centre trial. Thirty patients aged 50 years and over with insomnia will be recruited. The participants will take lemborexant (5 mg) at bedtime for 12 weeks and undergo a home-based sleep study at baseline and weeks 4 and 12, as well as MRI examinations to evaluate the glymphatic system at baseline and week 12. The primary outcome will be changes in objective sleep parameters as evaluated using a sleep monitoring system. The secondary outcomes will be changes in subjective sleep parameters. The relationships between changes in sleep parameters and the glymphatic system will be evaluated using diffusion tensor image analysis along the perivascular space, which is called the ALPS-index. Sleep parameters and the ALPS-index will be analysed using a paired t-test or Pearson's correlation coefficient.

ETHICS AND DISSEMINATION

The study protocol was approved by Nagoya University Certified Review Board. The findings from this research will be published in peer-reviewed journals and be presented at local, national and international conferences.

TRIAL REGISTRATION NUMBER

jRCTs041210024.

Long-read metagenomics of multiple displacement amplified DNA of low-biomass human gut phageomes by SACRA pre-processing chimeric reads

The human gut bacteriophage community (phageome) plays an important role in the host's health and disease; however, the entire structure is poorly understood, partly owing to the generation of many incomplete genomes in conventional short-read metagenomics. Here, we show long-read metagenomics of amplified DNA of low-biomass phageomes with multiple displacement amplification (MDA), involving the development of a novel bioinformatics tool, split amplified chimeric read algorithm (SACRA), that efficiently pre-processed numerous chimeric reads generated through MDA. Using five samples, SACRA markedly reduced the average chimera ratio from 72% to 1.5% in PacBio reads with an average length of 1.8 kb. De novo assembly of chimera-less PacBio long reads reconstructed contigs of ≥5 kb with an average proportion of 27%, which was 1% in contigs from MiSeq short reads, thereby dramatically improving contig length and genome completeness. Comparison of PacBio and MiSeq contigs found MiSeq contig fragmentations frequently near local repeats and hypervariable regions in the phage genomes, and those caused by multiple homologous phage genomes coexisting in the community. We also developed a reference-independent method to assess the completeness of the linear phage genomes. Overall, we established a SACRA-coupled long-read metagenomics robust to highly diverse gut phageomes, identifying high-quality circular and linear phage genomes with adequate sequence quantity.

Monitoring antibiotic resistance genes in wastewater treatment: Current strategies and future challenges

Antimicrobial resistance (AMR) is a growing threat to human and animal health. Progress in molecular biology has revealed new and significant challenges for AMR mitigation given the immense diversity of antibiotic resistance genes (ARGs), the complexity of ARG transfer, and the broad range of omnipresent factors contributing to AMR. Municipal, hospital and abattoir wastewater are collected and treated in wastewater treatment plants (WWTPs), where the presence of diverse selection pressures together with a highly concentrated consortium of pathogenic/commensal microbes create favourable conditions for the transfer of ARGs and proliferation of antibiotic resistant bacteria (ARB). The rapid emergence of antibiotic resistant pathogens of clinical and veterinary significance over the past 80 years has re-defined the role of WWTPs as a focal point in the fight against AMR. By reviewing the occurrence of ARGs in wastewater and sludge and the current technologies used to quantify ARGs and identify ARB, this paper provides a research roadmap to address existing challenges in AMR control via wastewater treatment. Wastewater treatment is a double-edged sword that can act as either a pathway for AMR spread or as a barrier to reduce the environmental release of anthropogenic AMR. State of the art ARB identification technologies, such as metagenomic sequencing and fluorescence-activated cell sorting, have enriched ARG/ARB databases, unveiled keystone species in AMR networks, and improved the resolution of AMR dissemination models. Data and information provided in this review highlight significant knowledge gaps. These include inconsistencies in ARG reporting units, lack of ARG/ARB monitoring surrogates, lack of a standardised protocol for determining ARG removal via wastewater treatments, and the inability to support appropriate risk assessment. This is due to a lack of standard monitoring targets and agreed threshold values, and paucity of information on the ARG-pathogen host relationship and risk management. These research gaps need to be addressed and research findings need to be transformed into practical guidance for WWTP operators to enable effective progress towards mitigating the evolution and spread of AMR.

Identification of Microbiota Biomarkers With Orthologous Gene Annotation for Type 2 Diabetes

Type 2 diabetes (T2D) is a systematic chronic metabolic condition with abnormal sugar metabolism dysfunction, and its complications are the most harmful to human beings and may be life-threatening after long-term durations. Considering the high incidence and severity at late stage, researchers have been focusing on the identification of specific biomarkers and potential drug targets for T2D at the genomic, epigenomic, and transcriptomic levels. Microbes participate in the pathogenesis of multiple metabolic diseases including diabetes. However, the related studies are still non-systematic and lack the functional exploration on identified microbes. To fill this gap between gut microbiome and diabetes study, we first introduced eggNOG database and KEGG ORTHOLOGY (KO) database for orthologous (protein/gene) annotation of microbiota. Two datasets with these annotations were employed, which were analyzed by multiple machine-learning models for identifying significant microbiota biomarkers of T2D. The powerful feature selection method, Max-Relevance and Min-Redundancy (mRMR), was first applied to the datasets, resulting in a feature list for each dataset. Then, the list was fed into the incremental feature selection (IFS), incorporating support vector machine (SVM) as the classification algorithm, to extract essential annotations and build efficient classifiers. This study not only revealed potential pathological factors for diabetes at the microbiome level but also provided us new candidates for drug development against diabetes.

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.

Recovery of human gut microbiota genomes with third-generation sequencing

Human gut microbiota modulates normal physiological functions, such as maintenance of barrier homeostasis and modulation of metabolism, as well as various chronic diseases including type 2 diabetes and gastrointestinal cancer. Despite decades of research, the composition of the gut microbiota remains poorly understood. Here, we established an effective extraction method to obtain high quality gut microbiota genomes, and analyzed them with third-generation sequencing technology. We acquired a large quantity of data from each sample and assembled large numbers of reliable contigs. With this approach, we constructed tens of completed bacterial genomes in which there were several new bacteria species. We also identified a new conditional pathogen, Enterococcus tongjius, which is a member of Enterococci. This work provided a novel and reliable approach to recover gut microbiota genomes, facilitating the discovery of new bacteria species and furthering our understanding of the microbiome that underlies human health and diseases.

Detection of a New Resistance-Mediating Plasmid Chimera in a blaOXA-48-Positive Klebsiella pneumoniae Strain at a German University Hospital

Mobile genetic elements, such as plasmids, facilitate the spread of antibiotic resistance genes in Enterobacterales. In line with this, we investigated the plasmid-resistome of seven bla_OXA-48 gene-carrying Klebsiella pneumoniae isolates, which were isolated between 2013 and 2014 at the University Medical Center in Göttingen, Germany. All isolates were subjected to complete genome sequencing including the reconstruction of entire plasmid sequences. In addition, phenotypic resistance testing was conducted. The seven isolates comprised both disease-associated isolates and colonizers isolated from five patients. They fell into two clusters of three sequence type (ST)101 and two ST11 isolates, respectively; and ST15 and ST23 singletons. The seven isolates harbored various plasmids of the incompatibility (Inc) groups IncF, IncL/M, IncN, IncR, and a novel plasmid chimera. All bla_OXA-48 genes were encoded on the IncL/M plasmids. Of note, distinct phenotypical resistance patterns associated with different sets of resistance genes encoded by IncL/M and IncR plasmids were observed among isolates of the ST101 cluster in spite of high phylogenetic relatedness of the bacterial chromosomes, suggesting nosocomial transmission. This highlights the importance of plasmid uptake and plasmid recombination events for the fast generation of resistance variability after clonal transmission. In conclusion, this study contributes a piece in the puzzle of molecular epidemiology of resistance gene-carrying plasmids in K. pneumoniae in Germany.

Massive expansion of human gut bacteriophage diversity

Bacteriophages drive evolutionary change in bacterial communities by creating gene flow networks that fuel ecological adaptions. However, the extent of viral diversity and its prevalence in the human gut remains largely unknown. Here, we introduce the Gut Phage Database, a collection of ∼142,000 non-redundant viral genomes (>10 kb) obtained by mining a dataset of 28,060 globally distributed human gut metagenomes and 2,898 reference genomes of cultured gut bacteria. Host assignment revealed that viral diversity is highest in the Firmicutes phyla and that ∼36% of viral clusters (VCs) are not restricted to a single species, creating gene flow networks across phylogenetically distinct bacterial species. Epidemiological analysis uncovered 280 globally distributed VCs found in at least 5 continents and a highly prevalent phage clade with features reminiscent of p-crAssphage. This high-quality, large-scale catalog of phage genomes will improve future virome studies and enable ecological and evolutionary analysis of human gut bacteriophages.

Metagenomics workflow for hybrid assembly, differential coverage binning, metatranscriptomics and pathway analysis (MUFFIN)

Metagenomics has redefined many areas of microbiology. However, metagenome-assembled genomes (MAGs) are often fragmented, primarily when sequencing was performed with short reads. Recent long-read sequencing technologies promise to improve genome reconstruction. However, the integration of two different sequencing modalities makes downstream analyses complex. We, therefore, developed MUFFIN, a complete metagenomic workflow that uses short and long reads to produce high-quality bins and their annotations. The workflow is written by using Nextflow, a workflow orchestration software, to achieve high reproducibility and fast and straightforward use. This workflow also produces the taxonomic classification and KEGG pathways of the bins and can be further used for quantification and annotation by providing RNA-Seq data (optionally). We tested the workflow using twenty biogas reactor samples and assessed the capacity of MUFFIN to process and output relevant files needed to analyze the microbial community and their function. MUFFIN produces functional pathway predictions and, if provided de novo metatranscript annotations across the metagenomic sample and for each bin. MUFFIN is available on github under GNUv3 licence: https://github.com/RVanDamme/MUFFIN.

Rebamipide ameliorates indomethacin-induced small intestinal damage and proton pump inhibitor-induced exacerbation of this damage by modulation of small intestinal microbiota

Non-steroidal anti-inflammatory drugs (NSAIDs) induce small intestinal damage. It has been reported that rebamipide, a mucoprotective drug, exerts a protective effect against NSAID-induced small intestinal damage; however, the underlying mechanism remains unknown. In this study, we investigated the significance of the small intestinal microbiota in the protective effect of rebamipide against indomethacin-induced small intestinal damage in mice. A comprehensive analysis of the 16S rRNA gene sequencing revealed an alteration in the composition of the small intestinal microbiota at the species level, modulated by the administration of rebamipide and omeprazole. The transplantation of the small intestinal microbiota of the mice treated with rebamipide suppressed the indomethacin-induced small intestinal damage. Omeprazole, a proton pump inhibitor, exacerbated the indomethacin-induced small intestinal damage, which was accompanied by the alteration of the small intestinal microbiota. We found that the transplantation of the small intestinal microbiota of the rebamipide-treated mice ameliorated indomethacin-induced small intestinal damage and the omeprazole-induced exacerbation of the damage. These results suggest that rebamipide exerts a protective effect against NSAID-induced small intestinal damage via the modulation of the small intestinal microbiota, and that its ameliorating effect extends also to the exacerbation of NSAID-induced small intestinal damage by proton pump inhibitors.

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.

Long-read metagenomics using PromethION uncovers oral bacteriophages and their interaction with host bacteria

Bacteriophages (phages), or bacterial viruses, are very diverse and highly abundant worldwide, including as a part of the human microbiomes. Although a few metagenomic studies have focused on oral phages, they relied on short-read sequencing. Here, we conduct a long-read metagenomic study of human saliva using PromethION. Our analyses, which integrate both PromethION and HiSeq data of >30 Gb per sample with low human DNA contamination, identify hundreds of viral contigs; 0-43.8% and 12.5-56.3% of the confidently predicted phages and prophages, respectively, do not cluster with those reported previously. Our analyses demonstrate enhanced scaffolding, and the ability to place a prophage in its host genomic context and enable its taxonomic classification. Our analyses also identify a Streptococcus phage/prophage group and nine jumbo phages/prophages. 86% of the phage/prophage group and 67% of the jumbo phages/prophages contain remote homologs of antimicrobial resistance genes. Pan-genome analysis of the phages/prophages reveals remarkable diversity, identifying 0.3% and 86.4% of the genes as core and singletons, respectively. Furthermore, our study suggests that oral phages present in human saliva are under selective pressure to escape CRISPR immunity. Our study demonstrates the power of long-read metagenomics utilizing PromethION in uncovering bacteriophages and their interaction with host bacteria.

CheckV assesses the quality and completeness of metagenome-assembled viral genomes

Millions of new viral sequences have been identified from metagenomes, but the quality and completeness of these sequences vary considerably. Here we present CheckV, an automated pipeline for identifying closed viral genomes, estimating the completeness of genome fragments and removing flanking host regions from integrated proviruses. CheckV estimates completeness by comparing sequences with a large database of complete viral genomes, including 76,262 identified from a systematic search of publicly available metagenomes, metatranscriptomes and metaviromes. After validation on mock datasets and comparison to existing methods, we applied CheckV to large and diverse collections of metagenome-assembled viral sequences, including IMG/VR and the Global Ocean Virome. This revealed 44,652 high-quality viral genomes (that is, >90% complete), although the vast majority of sequences were small fragments, which highlights the challenge of assembling viral genomes from short-read metagenomes. Additionally, we found that removal of host contamination substantially improved the accurate identification of auxiliary metabolic genes and interpretation of viral-encoded functions.

metaFlye: scalable long-read metagenome assembly using repeat graphs

Long-read sequencing technologies have substantially improved the assemblies of many isolate bacterial genomes as compared to fragmented short-read assemblies. However, assembling complex metagenomic datasets remains difficult even for state-of-the-art long-read assemblers. Here we present metaFlye, which addresses important long-read metagenomic assembly challenges, such as uneven bacterial composition and intra-species heterogeneity. First, we benchmarked metaFlye using simulated and mock bacterial communities and show that it consistently produces assemblies with better completeness and contiguity than state-of-the-art long-read assemblers. Second, we performed long-read sequencing of the sheep microbiome and applied metaFlye to reconstruct 63 complete or nearly complete bacterial genomes within single contigs. Finally, we show that long-read assembly of human microbiomes enables the discovery of full-length biosynthetic gene clusters that encode biomedically important natural products.

Metagenome-assembled genome binning methods with short reads disproportionately fail for plasmids and genomic Islands

Metagenomic methods enable the simultaneous characterization of microbial communities without time-consuming and bias-inducing culturing. Metagenome-assembled genome (MAG) binning methods aim to reassemble individual genomes from this data. However, the recovery of mobile genetic elements (MGEs), such as plasmids and genomic islands (GIs), by binning has not been well characterized. Given the association of antimicrobial resistance (AMR) genes and virulence factor (VF) genes with MGEs, studying their transmission is a public-health priority. The variable copy number and sequence composition of MGEs makes them potentially problematic for MAG binning methods. To systematically investigate this issue, we simulated a low-complexity metagenome comprising 30 GI-rich and plasmid-containing bacterial genomes. MAGs were then recovered using 12 current prediction pipelines and evaluated. While 82-94 % of chromosomes could be correctly recovered and binned, only 38-44 % of GIs and 1-29 % of plasmid sequences were found. Strikingly, no plasmid-borne VF nor AMR genes were recovered, and only 0-45 % of AMR or VF genes within GIs. We conclude that short-read MAG approaches, without further optimization, are largely ineffective for the analysis of mobile genes, including those of public-health importance, such as AMR and VF genes. We propose that researchers should explore developing methods that optimize for this issue and consider also using unassembled short reads and/or long-read approaches to more fully characterize metagenomic data.

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.