Skip to main content

Emergence of nosocomial associated opportunistic pathogens in the gut microbiome after antibiotic treatment

Abstract

Introduction

According to the Centers for Disease Control’s 2015 Hospital Acquired Infection Hospital Prevalence Survey, 1 in 31 hospital patients was infected with at least one nosocomial pathogen while being treated for unrelated issues. Many studies associate antibiotic administration with nosocomial infection occurrence. However, to our knowledge, there is little to no direct evidence of antibiotic administration selecting for nosocomial opportunistic pathogens.

Aim

This study aims to confirm gut microbiota shifts in an animal model of antibiotic treatment to determine whether antibiotic use favors pathogenic bacteria.

Methodology

We utilized next-generation sequencing and in-house metagenomic assembly and taxonomic assignment pipelines on the fecal microbiota of a urinary tract infection mouse model with and without antibiotic treatment.

Results

Antibiotic therapy decreased the number of detectable species of bacteria by at least 20-fold. Furthermore, the gut microbiota of antibiotic treated mice had a significant increase of opportunistic pathogens that have been implicated in nosocomial infections, like Acinetobacter calcoaceticus/baumannii complex, Chlamydia abortus, Bacteroides fragilis, and Bacteroides thetaiotaomicron. Moreover, antibiotic treatment selected for antibiotic resistant gene enriched subpopulations for many of these opportunistic pathogens.

Conclusions

Oral antibiotic therapy may select for common opportunistic pathogens responsible for nosocomial infections. In this study opportunistic pathogens present after antibiotic therapy harbored more antibiotic resistant genes than populations of opportunistic pathogens before treatment. Our results demonstrate the effects of antibiotic therapy on induced dysbiosis and expansion of opportunistic pathogen populations and antibiotic resistant subpopulations of those pathogens. Follow-up studies with larger samples sizes and potentially controlled clinical investigations should be performed to confirm our findings.

Introduction

Patients who acquire a nosocomial infection have approximately a 1 in 10 chance of mortality during their hospitalization [1]. In 2009, hospital acquired infections (HAIs) were estimated to cost about 45 billion dollars in direct costs [2]. One of the most common primary nosocomial infections is urinary tract infection (UTI). Recent research suggests oral antibiotic treatment of a primary nosocomial infection, like a UTI, is indirectly linked to secondary nosocomial antibiotic resistant infections [3, 4]. Secondary nosocomial infections are often multi-drug resistant opportunistic pathogens, such as, the Acinetobacter calcoaceticus/baumannii complex or Clostridioides difficile [5,6,7]. The human gut and intestinal microbiome are heavily populated with diverse microbial organisms. It is well known that antibiotic treatment reduces the diversity and number of bacteria in the digestive microbiome. Furthermore, it is well established that normal microbiota helps prevent the growth of pathogenic bacteria [8, 9]. We acknowledge the great amount of work showcasing antimicrobial therapy induced Clostridiodes difficile infection [10]. However, to our knowledge there is very little other direct experimental evidence that antibiotic treatment specifically selects for opportunistic pathogens associated with HAI in the gut microbiome.

The introduction of high-throughput sequencing technology has led to powerful advances in microbial community ecology composition and function assessments. Despite the extensive advancements in high-throughput sequencing there is a large portion of the human gut microbiome composition still uncharacterized. The two major sequencing methods used to assess microbial communities of human systems are amplicon sequencing and whole metagenome sequencing (WMS). While amplicon sequencing may provide enough information to assess composition, there are biases present, from amplification and primer selection, that cannot be ignored [11]. Additionally, WMS is the preferred method for strain assessment and functional profiling of diverse microbiomes. Traditionally, taxonomy profiling of WMS data relies on mapping each individual read to reference databases. However, a common issue with such a method is the large portion of unmappable reads from uncharacterized species and subspecies present in the microbiome. Most reference databases are built from cultivated microbes. Unfortunately, many species present in microbial ecosystems are not able to be cultivated and therefore reference genome databases only represent a portion of species present. Even culture-independent genomic approaches have a large presence of unexplored microbial populations [12, 13].

Recently we applied a traditional read-mapping based method to the metagenome analysis of antibiotic treated mice and revealed a decrease of gut microbiota diversity and enrichment of antibiotic resistance genes (ARGs) in the gut [14]. However, as with many other studies, the read-mapping method left a large portion of the reads uncharacterized. In this paper we implemented an inhouse pipeline that incorporates metagenome assembly, gene prediction, and lowest common ancestor taxonomic assignment (Fig. 1). Such a method not only greatly increased the portion of characterizable reads, but also expanded upon the direct evidence that antibiotic treatment for a common healthcare-associated infection selects for a gut antibiotic resistant reservoir of opportunistic pathogens associated with HAIs, in addition to the well-established digestive microbiome dysbiosis.

Fig. 1
figure1

Bioinformatics pipeline. While traditional metagenome analysis methods align reads with sequences in reference databases and perform taxonomy assignment directly, the method we used assembles reads into longer scaffolds before mapping to reference databases, and relies on lowest common ancestor (LCA) analysis to perform taxonomy assignment because scaffolds can potentially harbor multiple predicted genes. The longer length of scaffolds also makes it possible to classify scaffolds into one of the known, divergent, or novel categories based on their overall sequence similarity with reference databases. Publicly available tools used are annotated on top of boxes where used in the pipeline. Boxes in the center refer to publicly available databases used in the pipeline. NCBI: National Center for Biotechnology Information

Materials and methods

Animal studies

Experiments were conducted on BALB/c female mice (Taconic Biosciences, Derwood, MD) 8 – 10 weeks old. The same strain, sex and age mice were obtained from the same vendor and the same location in the vendor facility. Mice were randomly allocated to control or treatment groups (Ampicillin, Ciprofloxacin, or Fosfomycin treated groups), each with 4 animals. Not all animals yielded sufficient amount of sequencing reads (> ~ 10 million per sample) and in the end 3 animals per group were used for subsequent sequencing analysis. To recapitulate the most common nosocomial infection, a urinary tract infection model previously described [15,16,17] was used. All mice, except two of the three control mice, were inoculated with CFT073 uropathogenic E.coli (ATCC, Manassas, VA) to create the UTI model. Our previous study showed no discernable difference in the gut microbiome of naïve and urinary tract infection model mice [14]. Treatment groups include antibiotics Ampicillin, Ciprofloxacin, and Fosfomycin. Ampicillin treatment was with ampicillin trihydrate (Sigma, St. Louis, MO) 200 mg/kg in 0.1 M HCl twice daily with an 8-h interval for three consecutive days. Ciprofloxacin treatment was with a 5% oral suspension 50 mg/kg (Bayer HeathCare, Whippany, NJ) twice daily with an 8-h interval for three consecutive days. Fosfomycin treatment was with Monurol® Fosfomycin tromethamine 1000 mg/kg (Forest Pharmaceutical Inc., St Louis, MO) once daily for three consecutive days. Fecal samples to be sequenced were collected from the distal ileum and proximal colon after the third day of treatment. All animal studies were done in accordance with an approved protocol from the Institutional Animal Care and Use Committee of the White Oak Federal Research Center.

DNA extraction and sequencing

Genomic DNA extraction and sequencing was completed as previously described [14]. Briefly, collected fecal samples had DNA extracted using QIAamp DNA stool mini kit (Qiagen, Germantown, MD) with slight modifications to the protocol. Libraries were prepared with the Nextera DNA Library Prep Kit (Illumina, San Diego, CA) and sequenced on the NextSeq 500 Sequencing system (Illumina, San Diego, CA).

Host reads removal and metagenomic assembly

Sample outputs were processed into raw read fastq files with bcl2fastq2 (version 2.18.0.12, Illumina Inc.). Host reads were removed by first aligning reads to an indexed mouse genome (generated by downloading all sequences from the National Center for Biotechnology Information (NCBI) databases with the taxonomy ID of 10,090) using BWA-MEM (version 0.7.12), then unmapped reads were extracted using SAMtools (version 0.1.18) with parameters “-h -f 4”, and the output was converted to fastq with Picard-Tools (version 2.1.1) SamToFastq function. Assembly was performed with SPAdes (version 3.12.0) using the recommended parameters “-k 21,33,55,77” with the –meta flag to initiate MetaSPAdes mode for paired end samples (Control, Ampicillin and Ciprofloxacin cohorts) [18] and -s flag for single end samples (Fosfomycin) [19]. All scaffolds with greater than 500 bps were retained for downstream analysis. To further select for microorganism scaffolds another round of host read removal was performed.

Gene prediction and scaffold divergence classification

The MetaGeneAnnotator program [20] was used to predict genes within retained scaffolds for each sample. The resulting text file output was converted to BED format using an inhouse python script (https://github.com/FDA/metagenome). Gene sequences were extracted from the retained scaffold using BEDTools [21] and the BED file of results, generating a gene prediction results fasta file. We characterized each scaffold using methods described by Kowarsky et al. [22]. Briefly, scaffolds were characterized as either novel, divergent, or known. Novel scaffolds were all those with BLASTn alignments that spanned less than 20% of bases and had an average gene identity below 60%. Scaffolds were assigned to the known category if their average gene identity was greater than 80%. Divergent scaffolds were all those that neither fit into novel or known. To determine gene identity of each scaffold each samples’ predicted genes results fasta file was aligned with BLASTx to the NCBI nonredundant protein database (https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz). To determine alignment coverage scaffolds greater than 500bps in length were aligned with BLASTn to NCBI nucleotide database (https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nt.gz). BLASTn alignments to determine coverage were completed with the following parameters “-task megablast -evalue 1e-6 -max_target_seqs 1 -best_hit_score_edge 0.05 -best_hit_overhang 0.05 -window_size 0 -percent_identity 90”. BLASTx alignments to determine percent identity were completed with the following parameters “-max_target_seqs 1 -evalue 1e-6”. All BLAST alignments were completed using a local installation of the BLAST command line application BLAST + (version 2.3.0) [23].

Taxonomy, opportunistic pathogen, and antibiotic resistance profiling

Taxonomic assignment for scaffolds was performed in three steps. First, the top 10 hits from each predicted gene in a BLASTx alignment to the NCBI RefSeq database sorted by bitscore were retained and stored in an in-house developed SQLite database through the R package RSQLite (https://cran.r-project.org/package=RSQLite). The top scoring hits for a predicted gene with at least 60% positive scoring matches were retained for lowest common ancestor analysis. Second, the full length of each scaffold was mapped to the NCBI NT database through BLASTn, and any scaffold with over 90% percent identity with any mouse sequences were removed. Third, the taxon id of the lowest common ancestor of all the retained top BLASTx hits was assigned to the scaffold. To determine the abundance of opportunistic pathogens we counted how many reads of each sample aligned to the genome of opportunistic pathogens using Bowtie2 (version 2.1.0). For plotting purpose, normalization was performed by multiplying mapped read counts by 1 million and then dividing by the total number of reads from the respective sample (reads per million). To determine the profile of antibiotic resistance genes for specific species present we aligned predicted genes using BLASTx to the CARD database (August, 2019) [24]. Counts were tallied for every read that was aligned to the CARD database and was associated with a scaffold from an opportunistic pathogen.

Statistical analysis

ARG counts and counts of reads mapping to opportunistic pathogens were normalized using edgeR’s default normalization method (TMM normalization) and edgeR’s analysis functions, glmFit, and glmLRT were used to estimate the dispersion and test for differential abundance [25].

Results

Scaffold assembly increased percent of reads mapped

To characterize the microbiome, filtered reads were mapped using BLASTn to the NCBI nucleotide database. To increase mapping rates, we performed de novo assembly of reads into scaffolds. A total of 313,519 scaffolds greater than 500 bps were created for the Control cohort. Ampicillin, Ciprofloxacin, and Fosfomycin-treated groups had 12,438, 25,171, and 5,331 scaffolds greater than 500 bps generated, respectively. Assembly in conjunction with de novo gene prediction (see Methods) increased the average percent of reads mapped in the Control cohort by ~ 77%, and ~ 67% in Fosfomycin cohort (Fig. 2). Little change occurred in mapping rates for the Ampicillin and Ciprofloxacin Cohorts.

Fig. 2
figure2

Comparison of Mapped Reads before and after assembly. Pre-assembly mapped reads (Grey): percentage of reads that can be mapped to the NCBI nucleotide database. Post-assembly mapped reads (Gold): percentage of read that can either be mapped to the NCBI nucleotide database or be assembled into a gene that can be mapped to the NCBI nucleotide or protein databases. Control, Ampicillin, and Ciprofloxacin cohorts had scaffolds assembled using pair-end reads. Fosfomycin cohort contigs were assembled using single end reads. See Methods for details

A large portion of uncharacterized scaffolds are not present in antibiotic treated mice gut microbiomes

Through a characterization method described in Kowarsky et al.[22], we classified the assembled scaffolds into three categories: novel (BLASTn < 20% of nts and gene identity < 60%), known (BLASTn > 80% of nts), and divergent (neither novel or known; see Methods). We determined that 129,272 of the 313,519 Control scaffolds were collected from novel scaffolds likely to originate from genomes not fully characterized. Ampicillin-treated mice had 2,877 of their 12,438 scaffolds classified as novel. Interestingly, the Ciprofloxacin group only had 274 of its 25,171 scaffolds classified as novel. Lastly, Fosfomycin-treated group had 192 of its 5,331 scaffolds classified as novel. Because different scaffolds had different lengths, we further analyzed the portion of base pairs assigned to each category to normalize the effect of scaffold length. Control had 42% of scaffold base pairs assigned to the novel characterization (Fig. 3). Ampicillin, Ciprofloxacin, and Fosfomycin had 28%, 2%, and 1% of scaffold base pairs assigned to novel, respectively. In addition, Control had ~ 40% of their predicted genes come from novel scaffolds (Additional file 1: Table 1). Ampicillin, Ciprofloxacin, and Fosfomycin treatment groups had 24%, 2%, and 2% of their predicted genes come from novel scaffolds, respectively (Table 1). Markedly, antibiotic treated mice had a large reduction in scaffolds, predicted genes, and reads coming from novels scaffolds (Fig. 3, Additional file 1: Table 1). Further break down of the characterization of scaffolds for each cohort can be found in Additional file 1: Table 1.

Fig. 3
figure3

Composition of different categories of genomic scaffolds in each cohort. For each cohort, the assembled scaffolds are classified into one of the three categories: known, novel, and divergent. The classification is based on the degree of similarity between the predicted genes on scaffolds and known genes in the NCBI database. Novel scaffolds have BLASTn alignments that spanned less than 20% of bases and had an average gene identity below 60%. Known scaffolds had an average gene identity greater than 80%. Divergent scaffolds were all those that neither fit into novel or known (in Gene Prediction and Scaffold Divergence Classification section of Methods). The percentage of base pairs from scaffolds in each category relative to the total number of base pairs from all scaffolds is shown. Of note samples from the Fosfomycin cohort went through single-end sequencing while all other cohorts went through paired-end sequencing (see Methods and Discussion). The three categories of scaffolds (Novel, Known, and Divergent) were separately analyzed in downstream analysis

Table 1 Top species by the number of contigs for the control and antibiotic treatment cohorts

Antibiotic treatment selects for opportunistic pathogens

Higher level analysis was performed through taxonomic profiling based on a lowest common ancestor approach briefly described in methods. Lowest common ancestor assignment revealed that most control cohort scaffolds only had enough information to be assigned to the phylum level, with Firmicutes being the top phylum. Some control scaffolds were able to be identified at the species level and the top 10 species were from the Firmicutes phylum and Clostridiales order with the initiation of divergent branching occurring at the family level (Table 1). A total of 1,878 unique bacterial species with at least 5 scaffolds each were found in the control cohort. Conversely, Ampicillin, Ciprofloxacin, and Fosfomycin-treated groups only had 44, 108, and 49 unique bacterial species assigned with at least 5 scaffolds, respectively. Unlike the control cohort, the top 10 species present in each antibiotic treatment cohort had opportunistic pathogens present. Most notably, Ampicillin treated animals presented with Acinetobacter calcoaceticus/baumannii complex and Chlamydia abortus species of opportunistic pathogens (Table 1, highlighted in red). The top 10 species for Ciprofloxacin treated animals included opportunistic pathogens Chlamydia abortus, Clostridioides difficile, and Staphylococcus aureus (Table 1, highlighted in red). Fosfomycin’s top 10 species included Bacteroides fragilis and Bacteroides thetaiotaomicron, common gut microbes known for pathogenesis outside of the gut system (Table 1, highlighted in red). Importantly, some scaffolds we identified are likely from uncharacterized genomes (or epichromosomal elements) of specific pathogen species, as they harbor genes that are similar to, but distinct from, known genes of pathogen origin (Fig. 4).

Fig. 4
figure4

Gene track visual of two scaffolds that are likely from uncharacterized pathogen genomes (or epichromosomal elements). Labelled percentages are percent identity of query (predicted genes on the scaffolds) to subject (target proteins in the database). Transparent red boxes represent query coverage to subject. Top: an Acinetobacter calcoaceticus/baumannii scaffold (node 218) from the Ampicillin-treated cohort with a length of 1671 nucleotides (nts), with 3 predicted genes (green track) mapped to their targets. Bottom: a Clostridioides difficile scaffold (node 14,106) from the control cohort with a length of 2406 nts with 4 predicted genes (green track)

Quantitative analysis of read counts reveals opportunistic pathogen selection following antibiotic treatment

Employing a quantitative analysis (see Methods) based on reads mapped to predicted genes harbored by opportunistic pathogens, we determined that Chlamydia abortus and Acinetobacter calcoaceticus/baumannii complex had roughly 40-fold and 7-fold more normalized read counts in the Ampicillin group than in the Control group, respectively (Fig. 5). Chlamydia trachomatis, which was undetectable in Control group, registered a modest increase in normalized read counts in Ampicillin and Ciprofloxacin treated mice. The Fosfomycin cohort had > 30-fold more normalized reads mapped to both Bacteroides fragilis and Bacteroides thetaiotaomicron than the Control cohort. According to our analysis the opportunistic pathogen with the highest relative abundance (read counts relative to/normalized by the total microbiome reads from each sample) following antibiotic treatment was Bacteroides fragilis, which had on average 17,678 reads per million in the Fosfomycin-treated samples. This is a 46-fold increase in the relative abundance from the control group, where this opportunistic pathogen had on average 376 reads per million (Fig. 5).

Fig. 5
figure5

Opportunistic pathogen and antibiotic resistant gene abundance. The height of each bar in each panel indicates the normalized read counts (reads per million) for the opportunistic pathogen in a specific antibiotic-treated cohort. The blue and orange portion of the bar indicates the normalized read counts associated with ARGs (antibiotic resistant genes) and non-ARGs, respectively. Reads associated with ARGs (blue portion) were obtained by mapping the predicted genes to CARD (Comprehensive Antibiotic Resistance Database). *: Benjamini-Hochberg (BH) adjusted p value < 0.05

Higher number of antibiotic resistant genes in antibiotic selected opportunistic pathogen populations

Some top opportunistic pathogen species selected by antibiotics (Table 1) do not show an increase in relative abundance after antibiotic treatment compared to control (Fig. 5). To examine other changes of these opportunistic pathogen populations in antibiotic treated cohorts, we assessed the Antibiotic Resistant Gene (ARG) profile. We determined that the Acinetobacter calcoaceticus/baumannii subpopulations selected for in Ampicillin treated mice harbored a greater number of antibiotic resistant genes (328 ARGs) than the entire Acinetobacter calcoaceticus/baumannii population in the Control cohort (13 ARGs). The difference in antibiotic resistant gene diversity is even more stark when comparing these numbers relative to the total number of predicted genes in their respective sample. For example, the average number of relative Acinetobacter calcoaceticus/baumannii antibiotic resistant gene counts was 520 (normalized as per 10 K gene) in Ampicillin treated mice and 0.17 in the Control cohort (Benjamini-Hochberg (BH) adj. p value = 0.00075) (Fig. 6). Additionally, antibiotic resistant gene counts in Chlamydia abortus (BH adj. p value = 0.0001) and Chlamydia trachomatis (BH adj. p value = 0.00075) exhibited a statistically significant increase in Ampicillin cohort compared to Controls. We further examined the diversity profile of each treatment group and found that Ciprofloxacin treated mice had well over 500-fold more antibiotic resistant genes associated with Campylobacter coli (BH adj. p value = 0.0037) than Control (57 ARGs vs 0.09 ARGs, Fig. 6). Similar to Ampicillin cohort, we observed a statistically significant uptick in the number of antibiotic resistant genes associated with both species of Chlamydia in Ciprofloxacin treated mice. Lastly, Fosfomycin treated mice had an average of 10.8 and 6.2 antibiotic resistant genes in Bacteroides fragilis and Bacteroides thetaiotaomicron, respectively. For both species of Bacteroides, the number of relative antibiotic resistant gene counts in Fosfomycin treated mice is > fivefold more than that in the control.

Fig. 6
figure6

Opportunistic pathogen’s antibiotic resistance diversity profile. Y axis: number of genes associated with antibiotic resistance (normalized as per 10 K predicted genes). X axis: the control cohort and the three antibiotic-treated cohorts. *: Benjamini-Hochberg (BH) adjusted p value < 0.05

Discussion

Urinary tract infection (UTI) is the leading cause of hospital acquired infections worldwide [26]. Previous studies suggest UTIs make up roughly 40% of all nosocomial infections in the USA [27, 28]. However, recent guidelines for prevention of catheter-associated UTIs and implementation of evidence-based interventions have resulted in significant decreases in nosocomial UTIs [29]. Despite the improvements in best practices and nosocomial UTI rates, the Centers for Disease Control (CDC) reports that UTIs are still among the most common type of hospital-associated infection in the USA (https://www.cdc.gov/hai/ca_uti/uti.html). Oral antibiotic administration to treat nosocomial infections, like hospital acquired UTIs, have been implicated in gastrointestinal microbiota dysbiosis, and secondary infections in human and mouse studies [30,31,32,33,34]. Furthermore, these treatments often select for an antibiotic resistant population, likely from an antibiotic resistant (AR) reservoir [3].

Our results strongly support the previously described hypothesis with direct evidence of gut microbiota dysbiosis, opportunistic pathogen selection, and a higher percentage of selected opportunistic pathogens harboring antibiotic resistant genes occurring after oral antibiotic administration. Here, the reduction of our total sequencing reads, characterized scaffolds, and number of unique species present define the state of dysbiosis after antibiotic administration. Furthermore, the shift from Firmicutes to Proteobacteria, that we present in our results, is a known state of dysbiosis after antibiotic treatment [35] and during disease [36,37,38]. Importantly, our scaffold assembly and lowest common ancestor assignment approach allowed us to identify which top species were selected for after antibiotic treatment. The assembly of short reads into longer scaffolds before aligning to reference sequences decreased the number of ambiguous alignments and increased the number of reads that can be mapped to reference databases. The lowest common ancestor approach increased taxonomic resolution for longer scaffolds with multiple predicted genes. These methods can be applied to other metagenomic studies to potentially gather more information from the complex metagenome datasets compared to traditional read-mapping based methods. Many studies have correlated nosocomial infection of Acinetobacter calcoaceticus/baumannii complex after antibiotic administration for a primary nosocomial infection [4, 39]. We provide direct evidence for an increase of relative abundance of Acinetobacter calcoaceticus/baumannii, Chlamydia abortus, Chlamydia trachomatis, Bacteroides fragilis, and Bacteroides thetaiotaomicron after various antibiotic treatment. For other pathogens studied, for example Campylobacter coli in the Ciprofloxacin cohort, we found that while the relative abundance is not increased, the numbers of ARGs are, suggesting more ARG-enriched subpopulations were selected by the treatment. Interestingly, our previous study suggested that ARGs enriched in cohorts treated by one class of antibiotics are often resistant to other classes of antibiotics [14]. We also observed that some opportunistic pathogens have more ARG-associated reads in control compared to antibody-treated groups (Fig. 5), probably because the relative abundance of these pathogens was reduced so much after antibiotic treatment that any reads coming from them, including ARG-associated reads, became almost undetectable.

We acknowledge several limitations of our studies. First, the number of animals in each cohort is small, and some animals were not included in the analysis due to insufficient number of sequencing reads. It is unknown whether the excluded animals would have any impact on the final results. Second, paired-end sequencing was completed on all samples except the entire Fosfomycin cohort, which was single-end sequenced. This potentially had several effects on our Fosfomycin results. For example, the number of scaffolds generated that passed our quality control measures were far lower in Fosfomycin than other cohorts. This was likely due to the method of assembly differences necessary to build the scaffolds, as reported in our methods. Additionally, the assembly process may have contributed to the differences seen in the characterization of scaffolds, where a larger portion of control and Fosfomycin scaffolds were characterized as novel. This large portion of novel scaffolds is expected for the control cohort as there is substantial evidence for a large unknown population of microbiota present in unperturbed gut microbiome [40]. However, the assembly pipeline we used, metaSPAdes, is optimized to use pair-end reads and there is a possibility the assembly process of metagenomic single-end reads into scaffolds won’t maintain the high confidence found in the pair-end assembly protocol [18]. The decision to sequence Fosfomycin cohorts single-ended was made in the beginning of this project as it was not clear if paired-end sequencing would generate enough reads for a good coverage of the gut microbiome. As subsequent experiments revealed paired-end sequencing not only generated enough reads but also provided higher quality reads for assembly, we switched to paired-end sequencing for all other samples. However, we feel it is still important to report the findings from the single-end-sequenced Fosfomycin cohorts as the observed pattens (selection for opportunistic pathogens and enrichment of ARGs) are consistent with other cohorts with slightly different (paired-end sequencing) techniques. A third limitation in this pipeline is the difficulty in differentiating divergent contigs from plasmids and free phages which may have been present. This limitation may be mitigated with further downstream analysis to include a plasmid and phage identifier tool. A final limitation is the use of a mouse model as there are many differences at the cellular level as well as the gross anatomy level from a mouse to a human. However, overall the gut microbiota of mice and humans is mostly comprised of Bacteroidetes and Firmicutes [41]. Furthermore, the amount of knowledge on mouse gastroenterology, genetics, and microbiome composition provide ample advantage on any other model for microbiota research [42]. Therefore, we believe the findings from our mouse UTI model potentially revealed a general mechanism where the use of single oral antibiotic treatment directly induces the expansion of opportunistic pathogens and ARG-enriched subpopulations in the gut microbiome, which may subsequently lead to secondary nosocomial infections. While more work, especially controlled clinical studies, are warranted to confirm this mechanism in humans, this raised a possibility that we may need to consider combination antibiotic treatments to suppress the expansion of opportunistic pathogens during an antimicrobial therapy. New studies are being conducted to investigate this possibility.

Availability of data and materials

All of the raw sequencing reads were deposited in Sequence Read Archive (sra@ncbi.nlm.nih.gov) with an accession number SRP152866.

References

  1. 1.

    Magill SS, et al. Changes in prevalence of health care-associated infections in US hospitals. N Engl J Med. 2018;379(18):1732–44.

    Article  Google Scholar 

  2. 2.

    Scott, R.D. The direct medical costs of healthcare-associated infections in US hospitals and the benefits of prevention. 2009; Available from: https://www.cdc.gov/hai/pdfs/hai/scott_costpaper.pdf.

  3. 3.

    Magill SS, et al. Prevalence of antimicrobial use in US acute care hospitals, May-September 2011. J Am Med Assoc. 2014;312(14):1438–46.

    CAS  Article  Google Scholar 

  4. 4.

    Gulen TA, et al. Clinical importance and cost of bacteremia caused by nosocomial multi drug resistant Acinetobacter baumannii. Int J Infect Dis. 2015;38:32–5.

    Article  Google Scholar 

  5. 5.

    Peng Z, et al. Update on antimicrobial resistance in Clostridium difficile: resistance mechanisms and antimicrobial susceptibility testing. J Clin Microbiol. 2017;55(7):1998–2008.

    CAS  Article  Google Scholar 

  6. 6.

    Perez F, et al. Global challenge of multidrug-resistant Acinetobacter baumannii. Antimicrob Agents Chemother. 2007;51(10):3471–84.

    CAS  Article  Google Scholar 

  7. 7.

    García-Garmendia J-L, et al. Risk factors for Acinetobacter baumannii nosocomial bacteremia in critically ill patients: a cohort study. Clin Infect Dis. 2001;33(7):939–46.

    Article  Google Scholar 

  8. 8.

    Lupp C, et al. Host-mediated inflammation disrupts the intestinal microbiota and promotes the overgrowth of Enterobacteriaceae. Cell Host. 2007;2(2):119–29.

    CAS  Article  Google Scholar 

  9. 9.

    Mazmanian SK, Round JL, Kasper DL. A microbial symbiosis factor prevents intestinal inflammatory disease. Nature. 2008;453(7195):620.

    CAS  Article  Google Scholar 

  10. 10.

    Freeman J, et al. Comparison of the efficacy of ramoplanin and vancomycin in both in vitro and in vivo models of clindamycin-induced Clostridium difficile infection. J Antimicrob Chemother. 2005;56(4):717–25.

    CAS  Article  Google Scholar 

  11. 11.

    Claesson MJ, et al. Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res. 2010;38(22):e200–e200.

    Article  Google Scholar 

  12. 12.

    Nielsen HB, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat Biotechnol. 2014;32(8):822.

    CAS  Article  Google Scholar 

  13. 13.

    Rinke C, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013;499(7459):431.

    CAS  Article  Google Scholar 

  14. 14.

    Xu, L.S., A.; Raplee, I.; Chockalingam, A.; Stewart, S.; Walker, L.; Sacks, L.; Patel, V.; Li, Z.; Rouse, R., The effect of antibiotics on the gut microbiome: A metagenomics analysis of microbial shift and gut antibiotic resistance in antibiotic treated mice. BMC Genom, 2020.

  15. 15.

    Kucheria R, et al. Urinary tract infections: new insights into a common problem. Postgrad Med J. 2005;81(952):83.

    CAS  Article  Google Scholar 

  16. 16.

    Thai KH, Thathireddy A, Hsieh MH. Transurethral induction of mouse urinary tract infection. J Vis Exp. 2010;42:e2070.

    Google Scholar 

  17. 17.

    Chockalingam, A., et al., Evaluation of Immunocompetent Urinary Tract Infected Balb/C Mouse Model For the Study of Antibiotic Resistance Development Using Escherichia Coli CFT073 Infection. Antibiotics (Basel), 2019. 8(4).

  18. 18.

    Nurk S, et al. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34.

    CAS  Article  Google Scholar 

  19. 19.

    Bankevich A, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.

    CAS  Article  Google Scholar 

  20. 20.

    Noguchi H, Taniguchi T, Itoh T. MetaGeneAnnotator: detecting species-specific patterns of ribosomal binding site for precise gene prediction in anonymous prokaryotic and phage genomes. DNA Res. 2008;15(6):387–96.

    CAS  Article  Google Scholar 

  21. 21.

    Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

    CAS  Article  Google Scholar 

  22. 22.

    Kowarsky M, et al. Numerous uncharacterized and highly divergent microbes which colonize humans are revealed by circulating cell-free DNA. Proc Natl Acad Sci. 2017;114(36):9623–8.

    CAS  Article  Google Scholar 

  23. 23.

    Camacho C, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10(1):421.

    Article  Google Scholar 

  24. 24.

    Jia, B., et al., CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Research, 2016: p. gkw1004.

  25. 25.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    CAS  Article  Google Scholar 

  26. 26.

    Tambyah PA. Catheter-associated urinary tract infections: diagnosis and prophylaxis. Int J Antimicrob Agents. 2004;24:44–8.

    CAS  Article  Google Scholar 

  27. 27.

    Haley RW, et al. The nationwide nosocomial infection rate: a new need for vital statistics. Am J Epidemiol. 1985;121(2):159–67.

    CAS  Article  Google Scholar 

  28. 28.

    Haley RW, et al. Nosocomial infections in US hospitals, 1975–1976: estimated frequency by selected characteristics of patients. Am J Med. 1981;70(4):947–59.

    CAS  Article  Google Scholar 

  29. 29.

    Saint S, et al. A program to prevent catheter-associated urinary tract infection in acute care. N Engl J Med. 2016;374(22):2111–9.

    CAS  Article  Google Scholar 

  30. 30.

    Dethlefsen L, et al. The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 2008;6(11):e280.

    Article  Google Scholar 

  31. 31.

    Robinson CJ, Young VB. Antibiotic administration alters the community structure of the gastrointestinal microbiota. Gut Microbes. 2010;1(4):279–84.

    Article  Google Scholar 

  32. 32.

    Bohnhoff M, Drake BL, Miller CP. Effect of streptomycin on susceptibility of intestinal tract to experimental Salmonella infection. Proc Soc Exp Biol. 1954;86(1):132–7.

    CAS  Article  Google Scholar 

  33. 33.

    Barbut F, Petit J-C. Epidemiology of clostridium difficile-associated infections. Clin Microbiol. 2001;7(8):405–10.

    CAS  Article  Google Scholar 

  34. 34.

    Fekety R, Shah AB. Diagnosis and treatment of Clostridium difficile colitis. J Am Med Assoc. 1993;269(1):71–5.

    CAS  Article  Google Scholar 

  35. 35.

    Isaac S, et al. Short-and long-term effects of oral vancomycin on the human intestinal microbiota. J Antimicrob Chemother. 2016;72(1):128–36.

    Article  Google Scholar 

  36. 36.

    Zhang C, et al. Structural resilience of the gut microbiota in adult mice under high-fat dietary perturbations. ISME J. 2012;6(10):1848.

    CAS  Article  Google Scholar 

  37. 37.

    Wang J, et al. Modulation of gut microbiota during probiotic-mediated attenuation of metabolic syndrome in high fat diet-fed mice. ISME J. 2015;9(1):1.

    Article  Google Scholar 

  38. 38.

    Everard A, et al. Microbiome of prebiotic-treated mice reveals novel targets involved in host response during obesity. ISME J. 2014;8(10):2116.

    CAS  Article  Google Scholar 

  39. 39.

    Villers D, et al. Nosocomial Acinetobacter baumannii infections: microbiological and clinical epidemiology. Ann Intern Med. 1998;129(3):182–9.

    CAS  Article  Google Scholar 

  40. 40.

    Pasolli E, et al. Extensive unexplored human microbiome diversity revealed by over 150,000 genomes from metagenomes spanning age, geography, and lifestyle. Cell Host. 2019;176(3):649–62.

    CAS  Google Scholar 

  41. 41.

    Eckburg PB, et al. Diversity of the human intestinal microbial flora. Science. 2005;308(5728):1635–8.

    Article  Google Scholar 

  42. 42.

    Nguyen TLA, et al. How informative is the mouse for human gut microbiota research? Disease Mod Mech. 2015;8(1):1–16.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Leonard Sacks, Heather Stone, and Office of Medical Policy at US Food and Drug Administration (FDA) for their support and encouragement in pursuing this work. The authors are also grateful to the FDA High Performance Computing team, especially Mike Mikailov, for their scientific input and help with setting up the high performance computing environment for this work. This work is supported by the Research Participation Program at the Center for Drug Evaluation and Research, administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the US Department of Energy and the US Food and Drug Administration.

Funding

All funding for this project came from internal FDA grants or appropriated annual operating budgets for the Division of Applied Regulatory Science. These funding sources had no input on the design, execution, analysis, interpretation, or presentation of this work.

Author information

Affiliations

Authors

Corresponding author

Correspondence to Zhihua Li.

Ethics declarations

Disclaimer

This report is not an official US Food and Drug Administration guidance or policy statement. No official support or endorsement by the US Food and Drug Administration is intended or should be inferred.

Transparency declarations

None to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table 1

. Descriptive statistics of data processed by the analysis pipeline. For each sample in each cohort, the number of basepairs (BPs) in assembled contigs, the number of reads mapped to scaffolds, the number of predicted genes on scaffolds, and number of scaffolds classified into each of the three categories (novel, known, divergent) are given.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Raplee, I., Walker, L., Xu, L. et al. Emergence of nosocomial associated opportunistic pathogens in the gut microbiome after antibiotic treatment. Antimicrob Resist Infect Control 10, 36 (2021). https://doi.org/10.1186/s13756-021-00903-0

Download citation

Keywords

  • Antibiotic resistance
  • Metagenomics
  • Nosocomial
  • Microbiota
  • Antibiotic therapy