Complete hybrid genome assembly of clinical multidrug-resistant Bacteroides fragilis isolates enables comprehensive identification of antimicrobial-resistance genes and plasmids

Complete genome sequence of Peptostreptococcus porci isolated from porcine endocarditis in Japan

Peptostreptococcus porci is a recently described bacterium belonging to the Peptostreptococcaceae family, which was isolated in 2016 from pig intestine. Herein, we report the complete genome sequence of a clinical isolate of P. porci (GAI11004) obtained from porcine endocarditis in Japan. The genome contains a 2.4-Mb circular chromosome.

Pathogenic Bacteroides fragilis strains can emerge from gut-resident commensals

Bacteroides fragilis is a prominent member of the human gut microbiota, playing crucial roles in maintaining gut homeostasis and host health. Although it primarily functions as a beneficial commensal, B. fragilis can become pathogenic. To determine the genetic basis of its duality, we conducted a comparative genomic analysis of 813 B. fragilis strains, representing both commensal and pathogenic origins. Our findings reveal that pathogenic strains emerge across diverse phylogenetic lineages, due in part to rapid gene exchange and the adaptability of the accessory genome. We identified 16 phylogenetic groups, differentiated by genes associated with capsule composition, interspecies competition, and host interactions. A microbial genome-wide association study identified 44 genes linked to extra-intestinal survival and pathogenicity. These findings reveal how genomic diversity within commensal species can lead to the emergence of pathogenic traits, broadening our understanding of microbial evolution in the gut.

Dataset of 569 metagenome-assembled genomes from the caeca of multiple chicken breeds from commercial and backyard farming setups of Pakistan

This article focuses the recovery of prokaryotic organisms including bacteria and archaea from 9 different groups of chicken raised in different farm setups in Pakistan. The groups comprise of three different breeds (Broilers, White Layers, and Black Australorp) of chicken raised in different farming setups that include antibiotic-free control, commercial (open and controlled shed), and backyard farms. We have recovered 569 Metagenomics-Assembled Genomes (MAGs) with a completeness of ≥50 % and contamination of ≤10 %. For each MAG, functional annotations were obtained that include KEGG modules, carbohydrate active enzymes (CAZymes), peptidases, geochemical cycles, antibiotic resistance genes, stress genes, and virulence genes. Furthermore, two different sets of Single Copy Genes (SCGs) were used to construct the phylogenetic trees. Based on the reconstructed phylogeny, phylogenetic gain of each MAG is calculated to give an account of novelty.

Fifty years devoted to anaerobes: historical, lessons, and highlights

Renew interest and enthusiasm for anaerobes stem from both technological improvements (culture media, production of an adequate anaerobic atmosphere, identification methods) and greater awareness on the part of clinicians. Anaerobic infections were historically treated empirically, targeting the species known to be involved in each type of infection. Prevotella, fusobacteria, and Gram-positive cocci (GPAC) were considered responsible for infections above the diaphragm whereas for intra-abdominal infections, Bacteroides of the fragilis group (BFG), GPAC and clostridia were predominantly implicated. The antibiotic susceptibility of anaerobes was only taken into consideration by the clinician in the event of treatment failure or when faced with infections by multidrug-resistant bacteria (MDR). The evolution of antibiotic resistance together with clinical failures due to the absence of detection of hetero-resistant clones has resulted in a greater need for accessible antibiotic susceptibility testing (AST) and disc diffusion method. Improved isolation and identification of anaerobes, along with the availability of accessible and robust methods for performing AST, will ensure that treatment, whether empirical or guided by an antibiogram, will lead to better outcomes for anaerobic infections.

Roving methyltransferases generate a mosaic epigenetic landscape and influence evolution in Bacteroides fragilis group

Three types of DNA methyl modifications have been detected in bacterial genomes, and mechanistic studies have demonstrated roles for DNA methylation in physiological functions ranging from phage defense to transcriptional control of virulence and host-pathogen interactions. Despite the ubiquity of methyltransferases and the immense variety of possible methylation patterns, epigenomic diversity remains unexplored for most bacterial species. Members of the Bacteroides fragilis group (BFG) reside in the human gastrointestinal tract as key players in symbiotic communities but also can establish anaerobic infections that are increasingly multi-drug resistant. In this work, we utilize long-read sequencing technologies to perform pangenomic (n = 383) and panepigenomic (n = 268) analysis of clinical BFG isolates cultured from infections seen at the NIH Clinical Center over four decades. Our analysis reveals that single BFG species harbor hundreds of DNA methylation motifs, with most individual motif combinations occurring uniquely in single isolates, implying immense unsampled methylation diversity within BFG epigenomes. Mining of BFG genomes identified more than 6000 methyltransferase genes, approximately 1000 of which were associated with intact prophages. Network analysis revealed substantial gene flow among disparate phage genomes, implying a role for genetic exchange between BFG phages as one of the ultimate sources driving BFG epigenome diversity.

Rapid metagenomic sequencing for diagnosis and antimicrobial sensitivity prediction of canine bacterial infections

Antimicrobial resistance is a major threat to human and animal health. There is an urgent need to ensure that antimicrobials are used appropriately to limit the emergence and impact of resistance. In the human and veterinary healthcare setting, traditional culture and antimicrobial sensitivity testing typically requires 48-72 h to identify appropriate antibiotics for treatment. In the meantime, broad-spectrum antimicrobials are often used, which may be ineffective or impact non-target commensal bacteria. Here, we present a rapid, culture-free, diagnostics pipeline, involving metagenomic nanopore sequencing directly from clinical urine and skin samples of dogs. We have planned this pipeline to be versatile and easily implementable in a clinical setting, with the potential for future adaptation to different sample types and animals. Using our approach, we can identify the bacterial pathogen present within 5 h, in some cases detecting species which are difficult to culture. For urine samples, we can predict antibiotic sensitivity with up to 95 % accuracy. Skin swabs usually have lower bacterial abundance and higher host DNA, confounding antibiotic sensitivity prediction; an additional host depletion step will likely be required during the processing of these, and other types of samples with high levels of host cell contamination. In summary, our pipeline represents an important step towards the design of individually tailored veterinary treatment plans on the same day as presentation, facilitating the effective use of antibiotics and promoting better antimicrobial stewardship.

Comparative genomic assessment of members of genus Tenacibaculum: an exploratory study

Tenacibaculosis is an ulcerative skin disorder that affects finfish. It is caused by members of the genus Tenacibaculum, resulting in eccentric behavioural changes, including anorexia, lethargy, and abnormal swimming patterns that often result in mortality. Currently, species suspected of causing fish mortality include T. ovolyticum, T. gallaicum, T. discolor, T. finnmarkense, T. mesophilum, T. soleae, T. dicentrarchi, and T. maritimum. However, pathogenic members and the mechanisms involved in disease causation, progression, and transmission are limited due to the inadequate sequencing efforts in the past decade. In this study, we use a comparative genomics approach to investigate the characteristic features of 26 publicly available genomes of Tenacibaculum and report our observations. We propose the reclassification of "T. litoreum HSC 22" to the singaporense species and assignment of "T. sp. 4G03" to the species discolor (species with quotation marks have not been appropriately named). We also report the co-occurrence of several antimicrobial resistance/virulence genes and genes private to a few members. Finally, we mine several non-B DNA forming regions, operons, tandem repeats, high-confidence putative effector proteins, and sortase that might play a pivotal role in bacterial evolution, transcription, and pathogenesis.

Development and Clinical Application of a Multilocus Sequence Typing Scheme for Bacteroides fragilis Based on Whole-Genome Sequencing Data

Bacteroides fragilis is among the most abundant and pathogenic bacterial species in the gut microbiota and is associated with diarrheal disease in children, inflammatory bowel disease, and the development of colorectal cancer. It is increasingly resistant to potent antimicrobial agents such as carbapenems and metronidazole, making it among the most resistant anaerobic bacteria. These factors combined call for increased monitoring of B. fragilis and its population structure on a worldwide scale. Here, we present a possible solution through the development of a multilocus sequence typing scheme (MLST). The scheme is based on seven core gene fragments: groL (hsp60), rpoB, recA, dnaJ, rprX, prfA, and fusA. These gene fragments possess high discriminatory power while retaining concordance with whole core genome-based phylogenetic analysis. The scheme proved capable of differentiating B. fragilis isolates at the strain level. It offers a standardized method for molecular typing and can be applied to isolates from various sampling backgrounds, such as patient isolates, environmental samples, and strains obtained from food and animal sources. In total, 567 B. fragilis genomes were sequence typed and their isolate data collected. The MLST scheme clearly divided the B. fragilis population into two divisions based on the presence of the cfiA and cepA resistance genes. However, no other specific subpopulations within the analyzed genomes were found to be associated with any specific diseases or geographical location. With this MLST scheme, we hope to provide a powerful tool for the study and monitoring of B. fragilis on an international scale. IMPORTANCE Here, we present the first MLST scheme for Bacteroides fragilis, one of the most abundant pathogenic bacteria in the human gut microbiota. The scheme enables standard classification and monitoring of B. fragilis on a worldwide scale and groups the majority of current isolate data in one place. A more unified approach to the collection and analysis of B. fragilis data could provide crucial insights into how the pathogen operates and develops as a species. Close monitoring of B. fragilis is especially relevant, as it is increasingly resistant to potent antimicrobial agents and engages in horizontal gene transfer with other bacteria. Hopefully, this approach will guide new discoveries into how B. fragilis evolves and interacts with its human host. Additionally, the scheme could potentially be applied to other species of the genus Bacteroides.

Nanopore Electrochemistry for Pathogen Detection

Pathogen infections have seriously threatened human health, and there is an urgent demand for rapid and efficient pathogen identification to provide instructions in clinical diagnosis and therapeutic intervention. Recently, nanopore technology, a rapidly maturing technology which delivers ultrasensitive sensing and high throughput in real-time and at low cost, has achieved success in pathogen detection. Furthermore, the remarkable development of nanopore sequencing, for example, the MinION sequencer from Oxford Nanopore Technologies (ONT) as a competitive sequencing technology, has facilitated the rapid analysis of disease-related microbiomes at the whole-genome level and on a large scale. Here, we highlighted recent advances in nanopore approaches for pathogen detection at the single-molecule level. We also overviewed the applications of nanopore sequencing in pathogenic bacteria identification and diagnosis. In the end, we discussed the challenges and future developments of nanopore technology as promising tools for the management of infections, which may be helpful to aid understanding as well as decision-making.

Time for Some Group Therapy: Update on Identification, Antimicrobial Resistance, Taxonomy, and Clinical Significance of the Bacteroides fragilis Group

Bacteroides fragilis group (BFG) species are common members of the human microbiota that provide several benefits to healthy hosts, yet BFG are also the most common anaerobes isolated from human infections, including intra-abdominal infections, abscesses, and bloodstream infection. Compared to many other anaerobes associated with disease, members of the BFG are more likely to be resistant to commonly used antimicrobials, including penicillin (>90% resistant), carbapenems (2 to 20% resistant), and metronidazole (0.2 to 4% resistant). As a result, infection with BFG bacteria can be associated with poor clinical outcomes. Here, we discuss the role of BFG in human health and disease, proposed taxonomic reclassifications within the BFG, and updates in methods for species-level identification. The increasing availability of whole-genome sequencing (WGS) supports recent proposals that the BFG now span two families (Bacteroidaceae and "Tannerellaceae") and multiple genera (Bacteroides, Parabacteroides, and Phocaeicola) within the phylum Bacteroidota. While members of the BFG are often reported to "group" rather than "species" level in many clinical settings, new reports of species-specific trends in antimicrobial resistance profiles and improved resolution of identification tools support routine species-level reporting in clinical practice. Empirical therapy may not be adequate for treatment of serious infections with BFG, warranting susceptibility testing for serious infections. We summarize methods for antimicrobial susceptibility testing and resistance prediction for BFG, including broth microdilution, agar dilution, WGS, and matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). We examine global trends in BFG antimicrobial resistance and review genomics of BFG, revealing insights into rapid activation and dissemination of numerous antimicrobial resistance mechanisms.

Comparative Genomics of Bacteroides fragilis Group Isolates Reveals Species-Dependent Resistance Mechanisms and Validates Clinical Tools for Resistance Prediction

Bacteroides fragilis group (BFG) are the most frequently recovered anaerobic bacteria from human infections, and resistance to frontline antibiotics is emerging. In the absence of routine antimicrobial susceptibility testing (AST) for BFG in most clinical settings, we assessed the utility of clinical and modern genomics tools to determine BFG species-level identification and resistance patterns. A total of 174 BFG clinical isolates supplemented with 20 archived carbapenem-resistant B. fragilis sensu stricto (BFSS) isolates underwent antimicrobial susceptibility testing, MALDI-ToF mass-spectrometry, and whole-genome sequencing (WGS). Bruker BioTyper and VITEK-MS MALDI-ToF systems demonstrated accurate species-level identifications (91% and 90% agreement, respectively) compared to average nucleotide identity (ANI) analysis of WGS data. Distinct β-lactamase gene profiles were observed between BFSS and non-fragilis Bacteroides species, with significantly higher MICs to piperacillin-tazobactam in B. vulgatus and B. thetaiotaomicron relative to BFSS (P ≤ 0.034). We also uncovered phylogenetic diversity at the genomospecies level between division I and division II BFSS (ANI <0.95) and demonstrate that division II BFSS strains harbor an increased capacity to achieve carbapenem resistance through chromosomal activation of the CfiA carbapenemase. Finally, we report that CfiA detection by the Bruker BioTyper Subtyping Module accurately detected carbapenem resistance in BFSS with positive and negative percent agreement of 94%/90% and 95%/95% compared to ertapenem and meropenem susceptibility, respectively. These comparative analyses indicate that resistance mechanisms are distinct at both the phenotypic and genomic level across species within the BFG and that modern MALDI-ToF identification systems can be used for accurate species-level identification and resistance prediction of the BFG. IMPORTANCE Anaerobic infections present unique challenges in terms of detecting and identifying the etiologic agent and selecting the optimal antimicrobial therapy. Antimicrobial resistance is increasing in anaerobic pathogens, and it is critical to understand the prevalence and mechanisms of resistance to commonly prescribed antimicrobial therapies. This study uses comparative genomics to validate clinical tools for species-level identification and phenotypic resistance prediction in 194 isolates of Bacteroides fragilis group (BFG) bacteria, which represent the most commonly isolated organisms among anaerobic infections. We demonstrate species-specific patterns in antimicrobial resistance and validate new strategies for species-level organism identification and phenotypic resistance prediction in a routine clinical laboratory setting. These findings expand our understanding and management of anaerobic infections and justify further investigations into the molecular basis for species-specific resistance patterns observed within this study.

Paraphocaeicola brunensis gen. nov., sp. nov., Carrying Two Variants of nimB Resistance Gene from Bacteroides fragilis, and Caecibacteroides pullorum gen. nov., sp. nov., Two Novel Genera Isolated from Chicken Caeca

Three difficult-to-cultivate, strictly anaerobic strains, AN20^T, AN421^T, and AN502, were analyzed within a project studying possible probiotics for newly hatched chickens. Phylogenetic analyses showed that strains AN20^T, AN421^T, and AN502 formed two well-separated phylogenetic lineages in all phylogenetic and phylogenomic trees comprising members of the family Bacteroidaceae. Comparison to reference genomes of type species Bacteroides fragilis NCTC 9343^T, Phocaeicola abscessus CCUG 55929^T, and Capsularis zoogleoformans ATCC 33285^T showed low relatedness based on the calculated genome-to-genome distance and orthologous average nucleotide identity. Analysis of fatty acid profiles showed iso-C_15:0, anteiso-C_15:0, C_16:0, C_18:1ω9c, and iso-C_17:0 3OH as the major fatty acids for all three strains and additionally C_16:0 3OH for AN421^T and AN502. A specific combination of respiratory quinones different from related taxa was found in analyzed strains, MK-5 plus MK-11 in strain AN20^T and MK-5 plus MK-10 in strains AN421^T and AN502. Strains AN421^T and AN502 harbor complete CRISPR loci with CRISPR array, type II-C, accompanied by a set of cas genes (cas9, cas1, and cas2) in close proximity. Interestingly, strain AN20^T was found to harbor two copies of nimB gene with >95% similarity to nimB of B. fragilis, suggesting a horizontal gene transfer between these taxa. In summary, three isolates characterized in this study represent two novel species, which we proposed to be classified in two novel genera of the family Bacteroidaceae, for which the names Paraphocaeicola brunensis sp. nov. (AN20^T = CCM 9041^T = DSM 111154^T) and Caecibacteroides pullorum sp. nov. (AN421^T= CCM 9040^T = DSM 111155^T) are proposed.

IMPORTANCE This study represents follow-up research on three difficult-to-cultivate anaerobic isolates originally isolated within a project focused on strains that are able to stably colonize newly hatched chickens, thus representing possible probiotics. This project is exceptional in that it successfully isolates several miscellaneous strains that required modified and richly supplemented anaerobic media, as information on many gut-colonizing bacteria is based predominantly on metagenomic studies. Superior colonization of newly hatched chickens by Bacteroides spp., Phocaeicola spp., or related taxa can be considered of importance for development of future probiotics. Although different experiments can also be performed with provisionally characterized isolates, precise taxonomical definition is necessary for subsequent broad communication. The aim of this study is therefore to thoroughly characterize these isolates that represent novel genera and precisely determine their taxonomic position among related taxa to facilitate further research and communication involving these strains.

Complete Genome Sequences of Three Phocaeicola vulgatus Strains Isolated from a Healthy Japanese Individual

Phocaeicola vulgatus (formerly Bacteroides vulgatus) is a pathogenic anaerobic bacterium frequently involved in human infections. We present the complete genome sequences of three Phocaeicola vulgatus strains isolated from the same healthy person, determined by hybrid assembly using Nanopore long-read sequencing and DNBseq short-read sequencing.

Hybrid Assembly Provides Improved Resolution of Plasmids, Antimicrobial Resistance Genes, and Virulence Factors in Escherichia coli and Klebsiella pneumoniae Clinical Isolates

Emerging new sequencing technologies have provided researchers with a unique opportunity to study factors related to microbial pathogenicity, such as antimicrobial resistance (AMR) genes and virulence factors. However, the use of whole-genome sequence (WGS) data requires good knowledge of the bioinformatics involved, as well as the necessary techniques. In this study, a total of nine Escherichia coli and Klebsiella pneumoniae isolates from Norwegian clinical samples were sequenced using both MinION and Illumina platforms. Three out of nine samples were sequenced directly from blood culture, and one sample was sequenced from a mixed-blood culture. For genome assembly, several long-read, (Canu, Flye, Unicycler, and Miniasm), short-read (ABySS, Unicycler and SPAdes) and hybrid assemblers (Unicycler, hybridSPAdes, and MaSurCa) were tested. Assembled genomes from the best-performing assemblers (according to quality checks using QUAST and BUSCO) were subjected to downstream analyses. Flye and Unicycler assemblers performed best for the assembly of long and short reads, respectively. For hybrid assembly, Unicycler was the top-performing assembler and produced more circularized and complete genome assemblies. Hybrid assembled genomes performed substantially better in downstream analyses to predict putative plasmids, AMR genes and β-lactamase gene variants, compared to MinION and Illumina assemblies. Thus, hybrid assembly has the potential to reveal factors related to microbial pathogenicity in clinical and mixed samples.

MicroPIPE: validating an end-to-end workflow for high-quality complete bacterial genome construction

Background

Oxford Nanopore Technology (ONT) long-read sequencing has become a popular platform for microbial researchers due to the accessibility and affordability of its devices. However, easy and automated construction of high-quality bacterial genomes using nanopore reads remains challenging. Here we aimed to create a reproducible end-to-end bacterial genome assembly pipeline using ONT in combination with Illumina sequencing.

Results

We evaluated the performance of several popular tools used during genome reconstruction, including base-calling, filtering, assembly, and polishing. We also assessed overall genome accuracy using ONT both natively and with Illumina. All steps were validated using the high-quality complete reference genome for the Escherichia coli sequence type (ST)131 strain EC958. Software chosen at each stage were incorporated into our final pipeline, MicroPIPE.

Further validation of MicroPIPE was carried out using 11 additional ST131 E. coli isolates, which demonstrated that complete circularised chromosomes and plasmids could be achieved without manual intervention. Twelve publicly available Gram-negative and Gram-positive bacterial genomes (with available raw ONT data and matched complete genomes) were also assembled using MicroPIPE. We found that revised basecalling and updated assembly of the majority of these genomes resulted in improved accuracy compared to the current publicly available complete genomes.

Conclusions

MicroPIPE is built in modules using Singularity container images and the bioinformatics workflow manager Nextflow, allowing changes and adjustments to be made in response to future tool development. Overall, MicroPIPE provides an easy-access, end-to-end solution for attaining high-quality bacterial genomes. MicroPIPE is available at https://github.com/BeatsonLab-MicrobialGenomics/micropipe.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07767-z.

Genomic Background and Phylogeny of cfiA-Positive Bacteroides fragilis Strains Resistant to Meropenem-EDTA

Background: Bacteroides fragilis shows high antimicrobial resistance (AMR) rates and possesses numerous AMR mechanisms. Its carbapenem-resistant strains (metallo-β-lactamase cfiA-positive) appear as an emergent, evolving clade. Methods: This work examines the genomes, taxonomy, and phylogenetic relationships with respect to other B. fragilis genomes of two B. fragilis strains (CNM20180471 and CNM20200206) resistant to meropenem+EDTA and other antimicrobial agents. Results: Both strains possessed cfiA genes (cfiA14b and the new cfiA28), along with other AMR mechanisms. The presence of other efflux-pump genes, mexAB/mexJK/mexXY-oprM, acrEF/mdtEF-tolC, and especially cusR, which reduces the entry of carbapenem via the repression of porin OprD, may be related to meropenem–EDTA resistance. None of the detected insertion sequences were located upstream of cfiA. The genomes of these and other B. fragilis strains that clustered together in phylogenetic analyses did not meet the condition of >95% average nucleotide/amino acid identity, or >70% in silico genome-to-genome hybridization similarity, to be deemed members of the same species, although <1% difference in the genomic G+C content was seen with respect to the reference genome B. fragilis NCTC 9343T. Conclusions: Carbapenem-resistant strains may be considered a distinct clonal entity, and their surveillance is recommended given the ease with which they appear to acquire AMR.

Antimicrobial resistance in enteric bacteria: current state and next-generation solutions

Antimicrobial resistance is one of the largest threats to global health and imposes substantial burdens in terms of morbidity, mortality, and economic costs. The gut is a key conduit for the genesis and spread of antimicrobial resistance in enteric bacterial pathogens. Distinct bacterial species that cause enteric disease can exist as invasive enteropathogens that immediately evoke gastrointestinal distress, or pathobionts that can arise from established bacterial commensals to inflict dysbiosis and disease. Furthermore, various environmental reservoirs and stressors facilitate the evolution and transmission of resistance. In this review, we present a comprehensive discussion on circulating resistance profiles and gene mobilization strategies of the most problematic species of enteric bacterial pathogens. Importantly, we present emerging approaches toward surveillance of pathogens and their resistance elements as well as promising treatment strategies that can circumvent common resistance mechanisms.