The Gut Resistome during Hematopoietic Stem Cell Transplantation in Children ============================================================================ * Sarah M. Heston * Rebecca R. Young * Kirsten Jenkins * Paul L. Martin * Andre Stokhuyzen * Doyle V. Ward * Shakti K. Bhattarai * Vanni Bucci * Mehreen Arshad * Nelson J. Chao * Patrick C. Seed * Matthew S. Kelly ## Abstract **Background** Children undergoing hematopoietic stem cell transplantation (HCT) are at high risk of acquiring antibiotic-resistant bacteria. Few prior studies examined antibiotic resistance genes (ARGs) within the gut metagenomes of children undergoing HCT. **Methods** We conducted a longitudinal study of children (age <18 years) undergoing HCT at a single institution. We performed shotgun metagenomic sequencing of fecal samples collected between days -30 and +100 relative to HCT. We evaluated the effects of aerobic (cefepime, vancomycin, fluoroquinolones, aminoglycosides, macrolides, and trimethoprim-sulfamethoxazole) and anaerobic (piperacillin-tazobactam, carbapenems, metronidazole, and clindamycin) antibiotic exposures on the diversity and composition of the gut microbiome and resistome. **Findings** Using metagenomic data from 693 fecal samples collected from 80 children, we identified 350 unique ARGs. The most frequent ARGs identified encode resistance to tetracycline (n=91), beta-lactams (n=80), and fluoroquinolones (n=76). Both aerobic and anaerobic antibiotic exposures were associated with a decrease in the number of bacterial species (aerobic, β=0.72, 95% CI: 0.66, 0.79; anaerobic, β=0.68, 95% CI: 0.61, 0.76) and the number of unique ARGs (aerobic, β=0.83, 95% CI: 0.76, 0.91; anaerobic, β=0.84, 95% CI: 0.76, 0.93) within the gut metagenome. However, only anaerobic antibiotics were associated with an increase in the number of newly acquired ARGs (29%, 95% CI: 10%, 52%) and the abundance of ARGs (95%, 95% CI: 59%, 138%) in the gut resistome. Specific antibiotic exposures were associated with distinct changes in the number and abundance of resistance genes for individual antibiotic classes. **Interpretation** The gut metagenome and resistome of children are highly dynamic throughout HCT, driven largely by antibiotic exposures. Compared to antibiotics without anaerobic activity, anaerobic antibiotics were associated with increased microbiome instability and expansion of the gut resistome. **Funding** Antibacterial Resistance Leadership Group, National Institutes of Health, Duke Children’s Health & Discovery Initiative, Children’s Miracle Network Hospitals ## Introduction Antibiotic resistance is one of the most serious global public health threats. In the United States alone, more than 2.8 million antibiotic-resistant infections occur each year, resulting in an estimated $4.6 billion in healthcare expenditures and more than 35,000 deaths.1,2 Despite the World Health Organization’s Global Action Plan to combat antimicrobial resistance, the number of infections and deaths from antimicrobial-resistant organisms continues to rise.3 Currently, 18 antimicrobial-resistant organisms are considered national threats by the United States Centers for Disease Control and Prevention, and there has been a 50% increase in the number of infections caused by extended-spectrum beta lactamase (ESBL)-producing Enterobacteriaceae alone since 2013.1 A number of factors contribute to growing antibiotic resistance, but use of antibiotics is the major driving force, thereby making hospitals a setting for the spread of antibiotic resistance among patients.1,4 Patients undergoing hematopoietic stem cell transplantation (HCT) are at high risk of acquiring antibiotic-resistant organisms due to their prolonged hospitalizations, impaired immunity, and frequent receipt of broad-spectrum antibiotics. Colonization by multidrug-resistant bacteria is observed in more than 50% of adult HCT recipients and has been associated with a higher incidence of non-relapse-related mortality.5 Moreover, HCT recipients are at high risk for antibiotic-resistant infections arising from the gut because of antibiotic-induced changes in the gut microbiome, impaired host immunity, and mucosal barrier injury from chemotherapy.6,7 Overgrowth of potential pathogens that colonize the gut increases the risk of an antibiotic-resistant bloodstream infection.7,8 In addition, antibiotics with an anaerobic spectrum of activity are increasingly recognized to disrupt the gut microbiome and are associated with poor outcomes after HCT including a higher risk of graft-versus-host disease (GVHD).9 Historically, studies of HCT recipients identified colonization by antibiotic-resistant bacteria using culture-or PCR-based methods. These methods are limited by focusing only on cultivable bacteria or on a small number of specific resistance genes. With the increasing accessibility of shotgun metagenomic sequencing, we are now able to simultaneously detect the full complement of antibiotic resistance genes (ARGs) within a sample, collectively referred to as the “resistome.” In this study, we used shotgun metagenomic sequencing of 693 fecal samples from 80 children and adolescents undergoing HCT. We performed longitudinal analyses to evaluate the impact of aerobic and anaerobic antibiotics on gut microbiome composition and the gut resistome during HCT. In addition, we evaluated associations between specific antibiotic exposures and the number and abundances of resistance genes to specific antibiotic classes to describe the effect of specific antibiotic exposures on the gut resistome. ## Methods ### Study Design and Population We conducted a prospective cohort study of children and adolescents (<18 years of age) who underwent their first HCT through the Duke University Pediatric Transplant and Cellular Therapy Program between October 2015 and February 2018. Subjects were prospectively enrolled during the pretransplantation evaluation, and fecal samples were collected from subjects on as many days as possible from enrollment through 100 days after HCT. Daily clinical data, including information on receipt of antibiotics, were collected from the electronic medical record. We obtained informed consent from participants’ legal guardians prior to enrollment, and the study protocol was approved by the Duke University Health System Institutional Review Board. ### Transplant Practices Most children received trimethoprim-sulfamethoxazole (TMP-SMX) for *Pneumocystis jirovecii* prophylaxis from the start of the preparatory chemotherapy conditioning regimen through two days before HCT. Thereafter, children received monthly inhaled or IV pentamidine starting 30 days after HCT. After HCT, broad-spectrum antibiotics were initiated at onset of fever or concern for infection, but routine antibacterial prophylaxis was not used. Cefepime was the first-line antibiotic for subjects experiencing fever and neutropenia; vancomycin or antibiotics with an anaerobic spectrum of activity (e.g., piperacillin-tazobactam) were also frequently used if there was concern for a specific infectious source. Antibiotic selection and duration were at the discretion of the clinical provider. ### Metagenomic Sequencing and Data Processing Approximately weekly fecal samples from each subject underwent shotgun metagenomic sequencing. DNA was extracted from fecal samples using PowerSoil Pro Kits (Qiagen, Germantown, MD). DNA sequencing libraries were constructed using Nextera XT DNA Library Prep Kits (Illumina, San Diego, CA) and sequenced on NextSeq500 or NovaSeq6000 instruments (Illumina) as 150-bp paired-end reads. We trimmed reads using Trimmomatic (version 0.39) and removed host decontamination using Bowtie2 (version 2.3.5).10,11 Host-decontaminated reads were profiled for bacterial species relative abundances using MetaPhlAn2.12 Using the ShortBRED pipeline (version 0.9.4), we aligned sequencing reads to the Comprehensive Antibiotic Resistance Database (CARD; version 313).13,14 Reads mapped to antibiotic resistance genes were normalized by dividing the proportion of sequencing reads mapped to a gene sequence by the length of the gene sequence to determine the reads per kilobase per million (RPKM). We pruned samples with less than 500,000 quality-filtered metagenomic read-pairs. ### Statistical Analyses We used MaAsLin2 to fit linear mixed effects models with subject as a random effect to evaluate associations between antibiotic exposures and clinical characteristics and the relative abundances of specific bacterial species.15 For these models, we considered species that had a minimum mean relative abundance of 1% and a sample prevalence of 5%; the false discovery rate was set to 0.1. We considered cefepime, vancomycin, fluoroquinolones, aminoglycosides, macrolides, and TMP-SMX to have aerobic spectra of activity. We considered piperacillin-tazobactam, carbapenems, metronidazole, and clindamycin to have spectra of activity that included anaerobes. A participant was considered exposed to an antibiotic if they received at least one dose of the antibiotic since the last sequenced sample. If a subject started antibiotics on the day of fecal sample collection, the antibiotic exposure was assigned to the next sequenced sample. We assessed for the effects of time during HCT and antibiotic exposures on gut microbiome composition and the gut resistome using mixed effects models with day relative to HCT, aerobic antibiotic exposure, and anaerobic antibiotic exposure as fixed effects and subject as a random effect. Outcome measures for these models included the number of species and ARGs, the number of new species and ARGs acquired since the previous sample, microbiome and resistome instability measured by the Jaccard distance, and the abundance of ARGs.16 Models were adjusted for subject age, sex, underlying diagnosis, preparatory conditioning regimen, type of HCT (autologous or allogeneic), and sample sequencing depth. Models for the number of species and ARGs and the number of newly acquired species and ARGs were fit using a negative binomial distribution with the *glmmADMB* R package (version 0.8.3.3). Models for the Jaccard distance were fit using linear regression with the *lme4* R package (version 1.1-27.1).17,18 We performed similar analyses to evaluate associations between exposures to specific antibiotics and measures of the resistome using adjusted mixed effects models. We correlated measures of the resistome with corresponding measures of microbiome composition using repeated measures correlations implemented using the *rmcorr* R package (version 0.4.4).19 Finally, we considered the effect of individual antibiotics on the number and abundance of ARGs that confer resistance to specific classes of antibiotics using unadjusted mixed effect negative binomial models with individual antibiotic exposure as a fixed effect and subject as a random effect. We applied the Benjamini-Hochberg correction for multiple comparisons on the unadjusted models.20 All analyses were conducted using R statistical software (version 4.0.2). ## Results ### Patient and sample characteristics Our cohort included 80 children with a median (interquartile range [IQR]) age of 5.1 (2.2, 12.8) years (**Table 1**). The most frequent indication for HCT was hematologic malignancy (42%), and most transplants were allogeneic (88%), with umbilical cord blood being the most common donor source (64%). Cefepime, vancomycin, and TMP-SMX were the most commonly received antibiotics (**Table 2**). With regard to antibiotics with an anaerobic spectrum of activity, 44 (55%) subjects received metronidazole, 26 (33%) received piperacillin-tazobactam, 15 (19%) received a carbapenem, and four (5%) received clindamycin. Analyses included metagenomic sequencing data from 693 fecal samples [median (IQR) of 8 (6, 12) samples per child]. Shotgun metagenomic sequencing generated a total of 5,624,857,872 high-quality, paired-end metagenomic reads with a median (IQR) sequencing depth of 7.45 (3.41, 12.02) million reads per sample. View this table: [Table 1.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T1) Table 1. Characteristics of the study population of 80 children and adolescents undergoing HCT View this table: [Table 2.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T2) Table 2. Prevalences of antibiotic exposures among the study population and fecal samples collected after antibiotic exposures ### Composition of the gut microbiome during HCT We identified a total of 597 bacterial species with a median (IQR) of 28 (15, 50) unique bacterial species per sample. The most highly abundant bacterial species were *Escherichia coli* (mean relative abundance [MRA]=0.067), *Enterococcus faecalis* (MRA=0.060), *Bacteroides vulgatus* (MRA=0.056), *Enterococcus faecium* (MRA=0.039), and *Bacteroides stercoris* (MRA=0.033; **Figure 1**). The composition of the gut microbiome underwent substantial changes during the course of HCT; specifically, the relative abundances of several anaerobic species declined and the relative abundances of potential pathogens, including *Klebsiella pneumoniae, Enterococcus faecalis*, and several of the viridans group streptococci, increased during the course of HCT (**Supplementary Table 1**). The number of acquired bacterial species increased by 1% (95% confidence interval [CI]: 0.7%, 1.1%) per day after HCT (**Supplementary Table 2**). There were no significant changes in the total number of species or the instability of the gut microbiome with increasing time since HCT. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/07/07/2022.07.07.22277185/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/F1) Figure 1. Relative abundances of highly abundant genera over time. Relative abundances of highly abundant genera fecal samples of children by day relative to hematopoietic stem cell transplantation. Exposure to aerobic antibiotics was associated with a 28% (95% CI: 21%, 34%) decrease in the number of species and a 0.08 (95% CI: 0.04, 0.12) increase in instability. Anaerobic antibiotic exposures were similarly associated with a 32% (95% CI: 24%, 39%) decrease in the number of species and a 0.09 (0.05, 0.13) increase in instability of the gut microbiome composition (**Supplementary Table 2**). However, while aerobic antibiotics had no effect on the number of acquired species, there was a 19% (95% CI: 2%, 39%) increase in the number of acquired species associated with anaerobic antibiotics (**Supplementary Table 2**). Taken together, both aerobic and anaerobic antibiotic exposures were associated with a decrease in the number of bacterial species in the gut microbiome and an increase in instability of gut microbiome composition, but only anaerobic antibiotic exposures were associated with an increase in the acquisition of new species. We next sought to evaluate associations between exposures to specific antibiotics and changes in the composition of the gut microbiome (**Figure 2**). Cefepime exposure was associated with decreases in the relative abundances of both anaerobic bacteria from the genera *Bifidobacterium, Blautia*, and *Clostridium*, as well as potential pathogens including *Escherichia coli, Klebsiella pneumoniae*, and viridans group streptococci. Piperacillin-tazobactam and metronidazole, both anaerobic antibiotics, were associated with decreases in the relative abundances of several anaerobic bacterial species. Piperacillin-tazobactam was also associated with an increase in the relative abundance of *Enterococcus faecium*, and metronidazole was associated with increases in the relative abundances of *E. faecalis* and *Enterococcus gallinarum*. Notably, fluoroquinolone exposure tended to be associated with increases in the relative abundances of several bacterial species from the genera *Bacteroides, Bifidobacterium, Blautia*, and *Clostridium* (**Supplementary Table 1**). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/07/07/2022.07.07.22277185/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/F2) Figure 2. Heatmap depicting associations between antibiotic exposures and other clinical factors and gut microbiome composition. Associations were estimated from linear mixed effects models fit using MaAsLin2. We included bacterial species with a minimum mean relative abundance of 1% and a sample prevalence of 5%; the false discovery rate was set to 0.1. “Other” antibiotics includes exposures to tetracyclines, clindamycin, aminoglycosides, and non-cefepime beta-lactams. Betas indicate effect sizes and asterisks denote significant changes. Pip-Tazo, piperacillin-tazobactam; TMP-SMX, trimethoprim-sulfamethoxazole; Non-malignant heme, Non-malignant hematologic disorder; Non-myelo/RIC, non-myeloablative and reduced-intensity conditioning. ### The gut resistome during HCT We detected a total of 350 unique ARGs with a median (IQR) of 28 (15, 56) ARGs per sample. The most common ARGs detected were ErmB, tetO, tetW, and dfrF, each of which were present in more than 60% of samples and can confer resistance to macrolides and lincosamides (e.g., clindamycin), tetracyclines, and trimethoprim, respectively. The most common ARG classes were tetracyclines, beta-lactams, and fluoroquinolones (**Table 3**). View this table: [Table 3.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T3) Table 3. Antibiotic resistance genes identified in fecal samples that may confer resistance to clinically relevant antibiotics The gut resistome of children was also highly dynamic during the course of HCT. There was no significant change in the relative abundance of ARGs within the gut. However, there was a 0.2% (95% CI: 0.1%, 0.3%) increase in the number of ARGs and a 0.8% (95% CI: 0.6%, 1.0%) increase in the number of new ARGs with each day after HCT. Interestingly, the instability of the resistome decreased with each day after HCT (β=-0.001; 95% CI: -0.002, -0.000; **Figure 3; Supplementary Table 2**). While each year of increasing age was associated with a 5% (95% CI: 1%, 8%) decrease in gut resistome abundance and a decrease in resistome instability (β=-0.01; 95% CI: - 0.01, -0.00), other subject and transplant factors were not associated with measures of the resistome (**Supplementary Table 2**). Thus, when age and other clinical factors are held constant, new ARGs are acquired during HCT and the overall stability of the resistome increases with time. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/07/07/2022.07.07.22277185/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/F3) Figure 3. Changes to the gut resistome and microbiome over time. A. Number of ARGs from the resistome and number of bacterial species from the microbiome in fecal samples by the day relative to HCT; correlation of the number of bacterial species and number of ARGs for each sample. B. Number of new ARGs from the resistome and bacterial species from the microbiome by the day relative to HCT; correlation of the number of new bacterial species and ARGs since the previous sample. C. Jaccard distance in ARGs from the resistome and bacterial species from the microbiome by day relative to HCT; correlation of the Jaccard distances of the microbiome and resistome. The Jaccard distance is a measure of dissimilarity between two samples, with values closer to one representing increased instability.16 Points represent individual fecal samples; the smoothed line represents the median value over time, and shaded area represents the 95% confidence interval. Resistome is represented in blue, and the microbiome in red. ARG, antibiotic resistance gene; HCT, hematopoietic stem cell transplantation; rrm, repeated measures correlation coefficient. Antibiotic exposures were a main driver of changes in the gut resistome in pediatric HCT recipients. Compared to no aerobic antibiotic exposure, aerobic antibiotics were associated with a 17% (95% CI: 9%, 24%) decrease in the number of ARGs but not with a change in the number of acquired ARGs, the abundance of ARGs, or the instability of ARGs (**Supplementary Table 2**; **Figure 4**). Conversely, anaerobic antibiotic exposure led to a 16% (95% CI: 7%, 24%) decrease in the number of ARGs, a 29% (95% CI: 10%, 52%) increase in the number of new ARGs, a 95% (95% CI: 59%, 138%) increase in the abundance of ARGs, and a 0.07 (95% CI: 0.03, 0.12) increase in the instability of ARGs compared to no anaerobic antibiotic exposure (**Figure 4**). When considered individually, there were antibiotic-specific effects on measures of the resistome (**Figure 4**; **Supplementary Table 3**). Thus, while both aerobic and anaerobic antibiotics led to a decrease in the number of ARGs, only anaerobic antibiotic exposures were associated with more acquired ARGs and higher ARG abundances. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/07/07/2022.07.07.22277185/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/F4) Figure 4. Effect of aerobic and anaerobic antibiotics on measures of the resistome. The effect of aerobic and anaerobic antibiotics for each measure of the resistome: the number of ARGs per sample, the number of new ARGs per sample, the abundance of ARGs, and the instability of the ARGs as measured by the Jaccard distance. The Jaccard distance is a measure of dissimilarity between two samples, with values closer to one representing increased instability.16 Points represent estimates, and error bars denote the 95% confidence interval. ARG, antibiotic resistance gene; TMP-SMX, trimethoprim-sulfamethoxazole; Pip-Tazo, piperacillin-tazobactam. We next evaluated for correlations between measures of gut microbiome composition and the gut resistome using repeated measures correlation. There were moderate correlations between the number of ARGs and bacterial species (rrm=0.52; 95% CI: 0.46, 0.58), the number of acquired ARGs and bacterial species (rrm=0.53; 95% CI: 0.47, 0.58), and the stability of the gut microbiome and resistome (rrm=0.61; 95% CI: 0.56, 0.66; **Figure 3**). Taken together, these data suggest that the richness of the gut microbiome composition if a major driver of the gut resistome, and losses of gut microbial richness – as occurs with antibiotic exposure – in turn reduces the richness of the gut resistome. ### Antibiotic Exposures Affect Same and Different ARG Classes Finally, we determined the effect of specific antibiotic exposures on the number and abundance of ARGs of the same antibiotic class and other ARG classes (**Figure 5; Supplementary Tables 4 and 5**). Cefepime was associated with decreases in the number and abundance of ARGs for multiple ARG classes, including ARGs that confer resistance to beta-lactams, carbapenems, and fluoroquinolones. Piperacillin-tazobactam and metronidazole, both anaerobic antibiotics, led to a decrease in the number of ARGs from many antibiotic classes but with increases in the relative abundances of some classes of ARGs, especially glycopeptide ARGs that confer resistance to vancomycin. Fluoroquinolone exposure tended to be associated with an increased number of ARGs and an increased relative abundance of ARGs for several antibiotic classes. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/07/07/2022.07.07.22277185/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/F5) Figure 5. Heatmaps of clinically relevant antibiotic exposures’ effects on ARGs of the same and different antibiotic classes as the exposure. A. Antibiotic exposure effect on the number of ARGs for the exposure antibiotic class and others. B. Antibiotic exposure effect on the ARG abundance for the exposure antibiotic class and others. Beta indicates the effect size. ARGs, antibiotic resistance genes, FEP, cefepime; TZP, piperacillin-tazobactam; CBP, carbapenem; VAN, vancomycin; MTZ, metronidazole; FLQ, fluoroquinolone; SXT, trimethoprim-sulfamethoxazole; MAC, macrolide; Other, antibiotic exposure from tetracyclines, non-cefepime beta-lactams, clindamycin, and aminoglycosides. Tetracycline includes tetracycline and glycylcycline resistance genes; Beta-lactam includes cephalosporin, cephamycin, penam, and monobactam resistance genes; Lincosamide confers resistance to clindamycin; Carbapenem includes carbapenem and penem resistance genes; Glycopeptide confers resistance to vancomycin; Diaminopyrimidine confers resistance to trimethoprim; Nitromidazole confers resistance to metronidazole. Asterisks denote significant change. Gray boxes indicate there were not enough data to determine an effect. ## Discussion Using shotgun metagenomics, we described the dynamic changes in the gut microbiome and resistome that occur among children undergoing HCT. We described the effects of antibiotics on both the composition of the gut microbiome and resistome, including the distinct effects of aerobic and anaerobic antibiotics on the relative abundances of bacterial species within the gut microbiome and on the number and abundance of specific ARG classes. Compared to previous culture- and PCR-based methodologies, use of shotgun metagenomics for the characterization of the gut resistome is a novel way to study this complex, dynamic reservoir of ARGs.21 Considering ARGs collectively is particularly important because commensal bacteria, including those from the phyla Bacteroidetes and Firmicutes, likely harbor most of the ARGs in the healthy human gut and serve as a potential source for the horizontal transfer of ARGs to common enteric pathogens.21,22 Colonization by commensal, anaerobic bacteria with resistance genes may also explain how ARGs are detected even in the absence of antibiotic pressure.23 For instance, in our cohort, we detected a preponderance of ARGs that confer resistance to tetracyclines, macrolides, fluoroquinolones, and aminoglycosides despite the study population being infrequently exposed to these antibiotic classes. Knowledge of gut colonization by any bacteria with ARGs is especially important in the HCT population, as these patients are at high-risk of antibiotic-resistant bloodstream infections caused by gut-derived bacteria. Surprisingly little is known about the dynamics of the gut resistome following HCT. A small cohort study of eight children undergoing HCT for high-risk acute leukemia found that the most commonly detected ARGs conferred resistance to tetracyclines, macrolides, beta-lactams, and aminoglycosides, as was seen in our cohort.24 They also concluded that the four children who developed acute GVHD during the study had a different resistome signature, with acquisition of several new ARGs after HCT, compared to children without GVHD.24 However, this study did not assess the effect of antibiotics on the gut resistome.24 Our findings suggest that antibiotics significantly impact both gut microbiome composition and the gut resistome throughout HCT. Due to the temporal similarities in the microbiome composition and resistome, the correlations between the microbiome and the resistome, and the similar loss of ARGs and species associated with antibiotic exposures, we concluded that changes to the resistome are largely impacted by changes in the microbiome composition. Specifically, antibiotic exposures lead to loss of bacterial species from the gut microbiome and this loss of species and their accompanying ARGs results in contraction of the gut resistome. While both aerobic and anaerobic antibiotics were associated with a loss of microbiome and resistome richness, we found that aerobic and anaerobic antibiotics affected gut microbiome composition and the gut resistome differently. These findings are notable because anaerobic bacteria in the gut are primarily responsible for colonization resistance – the ability to prevent colonization by exogenous bacterial pathogens – through a variety of microbial and microbe-host interaction mechanisms.25 Thus, loss of anaerobic bacteria after anaerobic antibiotic exposures likely reduces colonization resistance, leading to expansion of existing antibiotic-resistant bacteria or the acquisition of exogenous antibiotic-resistant bacteria.26 Loss of colonization resistance with anaerobic antibiotics likely explains our findings of increased ARG acquisition, ARG abundance, and resistome instability despite the overall loss of ARGs associated with anaerobic antibiotics. When looking at specific anaerobic antibiotics, piperacillin-tazobactam and metronidazole were associated with lower relative abundances of *Bacteroides, Bifidobacterium, Blautia*, and *Clostridium* species, higher relative abundances of *Enterococcus* species, and an increased abundance of glycopeptide ARGs commonly found in enterococci. While we cannot directly assign ARGs to specific bacteria, these associations support the conclusion that exposure to these antibiotics results in loss of colonization resistance from anaerobes and colonization or expansion of *Enterococcus*. Taken together, our findings provide further evidence that anaerobic antibiotics should be used judiciously in this high-risk population.9 Our data also suggest that exposure to an antibiotic of one class can lead to acquisition of collateral antibiotic resistance to other antibiotic classes. In our cohort, vancomycin, piperacillin-tazobactam, and metronidazole exposures were associated with increases in the relative abundances of several other ARG classes, including fluoroquinolone, macrolide, tetracycline, and aminoglycoside ARGs. This phenomenon has been demonstrated in healthy adults, in which ciprofloxacin led to an increase in class D beta-lactamases and macrolide ARGs, and in premature neonates, in whom meropenem led to higher relative abundances of ARGs that confer resistance to fluoroquinolones, macrolides, tetracyclines, and trimethoprim.27,28 This expansion of collateral antibiotic resistance has likely also been observed clinically in the HCT population. In a meta-analysis of adults undergoing HCT, there was increased odds of a carbapenem-resistant *Klebsiella pneumoniae* infection not only with carbapenem exposure, but also with aminoglycoside, glycopeptide, and quinolone exposures.29 Similarly, infection with multidrug-resistant *Pseudomonas aeruginosa* was associated with vancomycin exposure in a separate meta-analysis of HCT recipients.30 Our data suggest that antibiotic pressure selects for organisms that are resistant to the exposure agent and also allows for the acquisition of new bacteria and expansion of existing bacteria with collateral ARGs of different antibiotic classes. The effect of fluoroquinolone exposure on collateral antibiotic resistance is of particular interest, given that levofloxacin is the preferred agent for antibacterial prophylaxis after HCT.31 Fluoroquinolone use after HCT has been associated with a lower risk of fever and neutropenia and bloodstream infections.31 However, it has more recently been associated with increased risks of colonization and infection by carbapenem-resistant Enterobacteriaceae.32 In a murine model, ciprofloxacin exposure led to an increase in the abundance of beta-lactamase ARGs.33 We found increases in both the richness and abundance of several different ARG classes with fluoroquinolone exposure. However, none of these associations met statistical significance, likely an effect of the lack of statistical power with relatively few fluoroquinolone-exposed samples from our cohort. On the contrary, in a study of 49 children with acute lymphoblastic leukemia, levofloxacin prophylaxis was associated with an increase in the prevalence of fluoroquinolone resistance, while a rise in collateral antibiotic resistance from other ARG classes was not detected.34 As levofloxacin prophylaxis is becoming more common in pediatric centers, further research is needed to determine the effect of levofloxacin on the resistome. Our study has several notable limitations. First, we are likely underestimating the effect of antibiotics on acquisition of ARGs given that antibiotic-induced changes in the gut microbiome are known to last weeks after cessation of antibiotic exposure.35 Second, as for many shotgun metagenomic sequencing analyses, the sequencing depth of our samples likely affected our ability to detect all ARGs present within a sample. To account for this, we included the sequencing depth in all models for our primary analyses. Unfortunately, analyses of associations between individual antibiotic exposures and ARG classes were not adjusted for sequencing depth because the adjusted models did not converge. Similarly, these models were not adjusted for time relative to HCT. Thus, while TMP-SMX exposure appears to reduce the acquisition of ARGs, this may reflect the fact that TMP-SMX was primarily administered prior to HCT and before the administration of broad-spectrum antibiotics for febrile neutropenia. Finally, with short-read sequencing data, we are unable to assign ARGs to specific bacterial species. Therefore, we can only make associations between gut microbiome composition and the gut resistome and cannot prove a direct link. In conclusion, this study provides the largest, longitudinal evaluation of the pediatric gut resistome to date. Understanding the complex dynamics of the gut resistome during HCT and its response to antibiotic exposures can impact clinical decisions for antibiotic prophylaxis and treatment. This study provides further evidence for the importance of antimicrobial stewardship efforts within the HCT community. Antibiotic exposures in this high-risk population should be thoughtful in order to prevent further disruption of the gut microbiome and resistome and to avoid promoting antibiotic resistance. ## Data Availability Data not included in the manuscript will be available by request to the corresponding authors. All metagenomic data will be uploaded to the Sequence Read Archive. ## Funding Research reported in this publication was supported in part by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health under Award Number UM1AI104681. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. This work was additionally supported by the National Institute of Child Health and Human Development of the National Institutes of Health [T32HD094671 to S.M.H.], the National Human Genome Research Institute of the National Institutes of Health [T32 HG008955 to S.M.H.], the Translating Duke Health Ending Disease Where it Begins Working Group [to S.M.H.], the local efforts of the Duke Children’s Office of Development and its Children’s Miracle Network Hospitals fundraising corporate partnerships and programs [to S.M.H.], and a National Institutes of Health Career Development Award [K23-AI135090 to M.S.K.]. ## Data Sharing Data not included in the manuscript will be available by request to the corresponding authors. All metagenomic data will be uploaded to the Sequence Read Archive. ## Figure Legends View this table: [Supplementary Table 1.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T4) Supplementary Table 1. Significant associations of antibiotic exposures and clinical characteristics with abundances of bacterial species within the gut microbiome View this table: [Supplementary Table 2.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T5) Supplementary Table 2. Effect of clinical factors on measures of the resistome and microbiome. View this table: [Supplementary Table 3.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T6) Supplementary Table 3. Effect of antibiotics on the resistome. View this table: [Supplementary Table 4.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T7) Supplementary Table 4. Effect of antibiotic exposure on the number of antibiotic resistance genes by gene class. View this table: [Supplementary Table 5.](http://medrxiv.org/content/early/2022/07/07/2022.07.07.22277185/T8) Supplementary Table 5. Effect of antibiotic exposure on the abundance of antibiotic resistance genes by gene class. ## Acknowledgements We thank the Duke University School of Medicine for the use of the Microbiome Core Facility, which provided DNA extraction, and the Sequencing and Genomic Technologies Shared Resource, which provided library preparation and shotgun metagenomic sequencing. We especially thank all the children and families who participated in this study. * Received July 7, 2022. * Revision received July 7, 2022. * Accepted July 7, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. Antibiotic resistance threats in the United States, 2019. Atlanta, Georgia, 2019 DOI:10.15620/cdc:82532. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15620/cdc:82532&link_type=DOI) 2. Nelson RE, Hatfield KM, Wolford H, et al. National estimates of healthcare costs associated with multidrug-resistant bacterial infections among hospitalized patients in the United States. Clin Infect Dis 2021; 72: S17–26. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa1581&link_type=DOI) 3. World Health Organization. Global action plan on antimicrobial resistance. 2015. 4. Chatterjee A, Modarai M, Naylor NR, et al. Quantifying drivers of antibiotic resistance in humans: a systematic review. Lancet Infect Dis 2018; 18: e368–78. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-099(18)30296-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 5. Scheich S, Lindner S, Koenig R, et al. Clinical impact of colonization with multidrug-resistant organisms on outcome after allogeneic stem cell transplantation in patients with acute myeloid leukemia. Cancer 2018; 124: 286–96. 6. Girmenia C, Bertaina A, Piciocchi A, et al. Incidence, Risk Factors and Outcome of Pre-engraftment Gram-Negative Bacteremia After Allogeneic and Autologous Hematopoietic Stem Cell Transplantation: An Italian Prospective Multicenter Survey. Clin Infect Dis 2017; 65: 1884–96. 7. Kelly MS, Ward D V, Severyn CJ, et al. Gut Colonization Preceding Mucosal Barrier Injury Bloodstream Infection in Pediatric Hematopoietic Stem Cell Transplantation Recipients. Biol Blood Marrow Transplant 2019; 25: 2274–80. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbmt.2019.07.019&link_type=DOI) 8. Taur Y, Xavier JB, Lipuma L, et al. Intestinal domination and the risk of bacteremia in patients undergoing allogeneic hematopoietic stem cell transplantation. Clin Infect Dis 2012; 55: 905–14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cis580&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22718773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 9. Tanaka JS, Young RR, Heston SM, et al. Anaerobic Antibiotics and the Risk of Graft-versus-Host Disease after Allogeneic Hematopoietic Stem Cell Transplantation. Biol Blood Marrow Transplant 2020; 26: 2053–60. 10. Bolger AM, Lohse M, Usadel B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014; 30: 2114–20. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu170&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24695404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340049100004&link_type=ISI) 11. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods 2012; 9: 357–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1923&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22388286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302218500017&link_type=ISI) 12. Truong DT, Franzosa EA, Tickle TL, et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods 2015; 12: 902–3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.3589&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26418763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 13. Kaminski J, Gibson MK, Franzosa EA, Segata N, Dantas G, Huttenhower C. High-Specificity Targeted Functional Profiling in Microbial Communities with ShortBRED. PLoS Comput Biol 2015; 11: e1004557. 14. Alcock BP, Raphenya AR, Lau TTY, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res 2020; 48: D517–25. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz935&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31665441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 15. Mallick H, Rahnavard A, McIver LJ, et al. Multivariable association discovery in population-scale meta-omics studies. PLoS Comput Biol 2021; 17: e1009442. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 16. Jaccard P. Etude comparative de la distribution florale dans une portion des Alpe’s et du Jura. Bull la Soci’et’e Vaudoise des Sci Nat 1901; 37: 547–79. 17. Fournier DA, Skaug HJ, Ancheta J, et al. AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models. Optim Methods Softw 2012; 27: 233–49. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/10556788.2011.597854&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305292200005&link_type=ISI) 18. Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw 2015; 67. DOI:10.18637/jss.v067.i01. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v067.i01&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23757445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 19. Bakdash JZ, Marusich LR. Repeated Measures Correlation. Front Psychol 2017; 8: 456. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpsyg.2017.00456&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28439244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 20. Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc Ser B 1995; 57: 289–300. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730414553029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25326247&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 21. van Schaik W. The human gut resistome. Philos Trans R Soc Lond B Biol Sci 2015; 370: 20140087. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2014.0087&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25918444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 22. Sommer MOA, Dantas G, Church GM. Functional characterization of the antibiotic resistance reservoir in the human microflora. Science (80-) 2009; 325: 1128–31. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjUvNTk0NC8xMTI4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDcvMDcvMjAyMi4wNy4wNy4yMjI3NzE4NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 23. Fouhy F, Ogilvie LA, Jones B V, et al. Identification of aminoglycoside and β-lactam resistance genes from within an infant gut functional metagenomic library. PLoS One 2014; 9: e108016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0108016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25247417&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 24. D’Amico F, Soverini M, Zama D, et al. Gut resistome plasticity in pediatric patients undergoing hematopoietic stem cell transplantation. Sci Rep 2019; 9: 1–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-019-40055-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30862893&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 25. Buffie CG, Pamer EG. Microbiota-mediated colonization resistance against intestinal pathogens. Nat Rev Immunol 2013; 13: 790–801. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri3535&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24096337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom) 26. Boutrot M, Azougagh K, Guinard J, Boulain T. Antibiotics with activity against intestinal anaerobes and the hazard of acquired colonization with ceftriaxone-resistant Gram-negative pathogens in ICU patients: A propensity score-based analysis. J Antimicrob Chemother 2019; 74: 3095–103. 27. Willmann M, El-Hadidi M, Huson DH, et al. Antibiotic Selection Pressure Determination through Sequence-Based Metagenomics. Antimicrob Agents Chemother 2015; 59: 7335–45. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjEwOiI1OS8xMi83MzM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDcvMDcvMjAyMi4wNy4wNy4yMjI3NzE4NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 28. Gibson MK, Wang B, Ahmadi S, et al. Developmental dynamics of the preterm infant gut microbiota and antibiotic resistome. Nat Microbiol 2016; 1: 16024. 29. Liu P, Li X, Luo M, et al. Risk Factors for Carbapenem-Resistant Klebsiella pneumoniae Infection: A Meta-Analysis. Microb Drug Resist 2018; 24: 190–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/mdr.2017.0061&link_type=DOI) 30. Raman G, Avendano EE, Chan J, Merchant S, Puzniak L. Risk factors for hospitalized patients with resistant or multidrug-resistant Pseudomonas aeruginosa infections: a systematic review and meta-analysis. Antimicrob Resist Infect Control 2018; 7: 79. 31. Lehrnbecher T, Fisher BT, Phillips B, et al. Guideline for antibacterial prophylaxis administration in pediatric cancer and hematopoietic stem cell transplantation. Clin Infect Dis 2020; 71: 226–36. 32. Predic M, Delano JP, Tremblay E, Iovine N, Brown S, Prins C. Evaluation of patient risk factors for infection with carbapenem-resistant Enterobacteriaceae. Am J Infect Control 2020; 48: 1028–31. 33. Xu L, Surathu A, Raplee I, et al. The effect of antibiotics on the gut microbiome: a metagenomics analysis of microbial shift and gut antibiotic resistance in antibiotic treated mice. BMC Genomics 2020; 21: 263. 34. Margolis EB, Hakim H, Dallas RH, et al. Antibiotic prophylaxis and the gastrointestinal resistome in paediatric patients with acute lymphoblastic leukaemia: a cohort study with metagenomic sequencing analysis. The Lancet Microbe 2021; 2: e159–67. 35. Suez J, Zmora N, Zilberman-Schapira G, et al. Post-Antibiotic Gut Mucosal Microbiome Reconstitution Is Impaired by Probiotics and Improved by Autologous FMT. Cell 2018; 174: 1406–1423.e16. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2018.08.047&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30193113&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F07%2F2022.07.07.22277185.atom)