Lei Liu1†
Qing Wang1,2†
Huai Lin1
Ranjit Das1
Siyi Wang1
Hongmei Qi1
Jing Yang1
Yingang Xue3
Daqing Mao4*
Yi Luo1*- 1College of Environmental Science and Engineering, Ministry of Education Key Laboratory of Pollution Processes and Environmental Criteria, Nankai University, Tianjin, China
- 2Hebei Key Laboratory of Air Pollution Cause and Impact (preparatory), College of Energy and Environmental Engineering, Hebei University of Engineering, Handan, China
- 3Key Laboratory of Environmental Protection of Water Environment Biological Monitoring of Jiangsu Province, Changzhou Environmental Monitoring Center, Changzhou, China
- 4School of Medicine, Nankai University, Tianjin, China
Antibiotics are frequently used to treat bacterial infections; however, they affect not only the target pathogen but also commensal gut bacteria. They may cause the dysbiosis of human intestinal microbiota and consequent metabolic alterations, as well as the spreading of antibiotic resistant bacteria and antibiotic resistance genes (ARGs). In vitro experiments by simulator of the human intestinal microbial ecosystem (SHIME) can clarify the direct effects of antibiotics on different regions of the human intestinal microbiota, allowing complex human microbiota to be stably maintained in the absence of host cells. However, there are very few articles added the antibiotics into this in vitro model to observe the effects of antibiotics on the human intestinal microbiota. To date, no studies have focused on the correlations between the bloomed pathogens caused by amoxicillin (AMX) exposure and increased functional pathway genes as well as ARGs. This study investigated the influence of 600 mg day–1 AMX on human intestinal microbiota using SHIME. The impact of AMX on the composition and function of the human intestinal microbiota was revealed by 16S rRNA gene sequencing and high-throughput quantitative PCR. The results suggested that: (i) AMX treatment has tremendous influence on the overall taxonomic composition of the gut microbiota by increasing the relative abundance of Klebsiella [linear discriminant analysis (LDA) score = 5.26] and Bacteroides uniformis (LDA score = 4.75), as well as taxonomic diversity (Simpson, P = 0.067, T-test; Shannon, P = 0.061, T-test), and decreasing the members of Parabacteroides (LDA score = 4.18), Bifidobacterium (LDA score = 4.06), and Phascolarctobacterium (LDA score = 3.95); (ii) AMX exposure significantly enhanced the functional pathway genes and beta-lactam resistance genes, and the bloomed pathogens were strongly correlated with the metabolic and immune system diseases gene numbers (R = 0.98, P < 0.001) or bl2_len and bl2be_shv2 abundance (R = 0.94, P < 0.001); (iii) the changes caused by AMX were “SHIME-compartment” different with more significant alteration in ascending colon, and the effects were permanent, which could not be restored after 2-week AMX discontinuance. Overall results demonstrated negative side-effects of AMX, which should be considered for AMX prescription.
Introduction
Human intestinal microbiota co-exists in symbiosis with human beings and comprises with about 150 times more genes than the human genome, which makes intestinal microbiota become “another” genome of human beings (Qin et al., 2010). Moreover, human intestinal microbiota has been demonstrated to provide numerous important functions for the human health, including fermentation of indigestible dietary polysaccharides, synthesis of essential amino acids, and vitamins, modulation of the immune function and protection from the pathogens, as well as metabolism of the xenobiotic drugs (Chow et al., 2010; Yatsunenko et al., 2012; Cabreiro et al., 2013). Both human and veterinary antibiotics were detected in the collective gut of the Chinese population through our previous research (Wang et al., 2020). The antibiotics in our gut would kill or prevent the growth of commensal beneficial bacteria, which would cause the dysbiosis of intestinal microbiota, resulting multiple human diseases (Blaser, 2016; Ianiro et al., 2016; Flandroy et al., 2018; Kho and Lal, 2018). Antibiotics may also promote the spreading of antibiotic resistant bacteria (ARB) and antibiotic resistance genes (ARGs) in the human gut (Stecher et al., 2013; Barraud et al., 2018). One of the most commonly prescribed antibiotic in clinical and residential applications is amoxicillin (AMX), an inexpensive oral penicillin-type, beta-lactam antibiotic that kills a broad spectrum of bacteria by interfering with the synthesis of bacterial cell wall peptidoglycan layer (Zapata and Quagliarello, 2015; Zhang et al., 2015). The influence of AMX on human intestinal microbiota by in vivo research has been well studied and recorded, including human (Ladirat et al., 2014a; Pallav et al., 2014; Vrieze et al., 2014; Zaura et al., 2015; Oh et al., 2016) and human microbiota associated-animal models (Barc et al., 2008; Collignon et al., 2008; Collignon et al., 2010). Oh et al. (2016) noticed that AMX increased the abundance of Klebsiella in human intestinal microbiota. As a typical Gram-negative bacterial pathogen, K. pneumoniae is ubiquitous in the environment and symbioses in the human gut (Rock et al., 2014). However, the bloom of K. pneumoniae in human gut is often related to unhealthy status (Pena et al., 1997; Wiener-Well et al., 2010; Gorrie et al., 2017). Furthermore, no studies have suggested the bloom of K. pneumoniae caused by AMX exposure contributed to the increase of functional pathway genes and beta-lactam resistance genes.
Previous reports have revealed that changes of intestinal microbiota by AMX exposure are region-dependent and more significant effects were observed in the proximal colon than in the distal colon (Ladirat et al., 2014b; Marzorati et al., 2017). In vivo experiments detected the fecal samples that generally stand for the distal intestinal microbiota, which could not reveal the impacts of antibiotics on different gut regions. Moreover, in vitro experiment using simulator of the human intestinal microbial ecosystem (SHIME) can elucidate the direct effects of medicines on different regions of the human intestinal microbiota, which allows the complex human microbiota to be stably maintained in the absence of host cells (Van de Wiele et al., 2015). The SHIME model is known to be a useful tool for the in vitro study, which has already been used to identify the influence of bacteria and compounds on the colon microbiota, including probiotics such as Clostridium cluster XIVa and Bifidobacterium longum (Van den Abbeele et al., 2013; Truchado et al., 2015), prebiotics and prebiotics like compounds such as inulin, polyphenols and orange juice (Kemperman et al., 2013; Duque et al., 2016; Selak et al., 2016), and other toxic compounds such as Chlorpyrifos and Arsenic (Reygner et al., 2016; Yu et al., 2016). However, to the best of our knowledge, there are very few research that applied antibiotics into the SHIME model (Van den Abbeele et al., 2012; Bussche et al., 2015; Marzorati et al., 2017; Ichim et al., 2018; El Hage et al., 2019; Liu et al., 2020). These articles focused on the benefit of the mucosal environment, high-fiber diets, probiotic, and propionate-producing consortium in human intestinal microbiota. The effects of antibiotics, including AMX and other antibiotics mixture, vancomycin, and clindamycin, were limited on microbiota composition and metabolite. Also, there was little information available to demonstrate the link between the bloomed pathogens and functional pathway genes or ARGs. Hence, there is a need to study the influence of AMX on human intestinal microbiota in the SHIME model that integrates the entire gastrointestinal tract and maintains microbiome stability over an extended timeframe (Van de Wiele et al., 2015).
Since the reasonable dosage of AMX for the adult human study is about 750 to 1,500 mg day–1, and only half volume of the adult gut can be simulated in the SHIME model, here the direct effects of 600 mg day–1 of AMX on the composition and function of the human fecal microbiota were followed by the previous studies (Pallav et al., 2014; Reijnders et al., 2016). Three reactors (representing the ascending, transverse, and descending colon, respectively) that inoculated with human intestinal microbiota were studied for the three groups: immediately before AMX administration during 0 to 21 days (a control group), AMX-exposure during 22 to 28 days (an AMX treated group), and after the AMX discontinuance during 29 to 42 days (a recovery group). The 16S rRNA gene sequencing and high-throughput quantitative PCR (HT-qPCR) results revealed that AMX exposure caused a tremendous impact on the overall taxonomic composition of the gut microbiota, and increased functional pathway genes as well as ARGs. The changes were “SHIME-compartment” different with more significant modulation in the ascending colon, which could not be restored after 2-week AMX discontinuance. Therefore, the results of our research demonstrated a severe impact and a negative side-effect of AMX related to health problems, which should be considered as a fundamental aspect of the cost-benefit equation for its prescription.
Materials and Methods
The SHIME Model Experimental Setup and Sampling
The SHIME set up was formed by five double-jacketed reactors designated as the stomach, small intestine, ascending colon, transverse colon, and descending colon, respectively (Supplementary Figure S1). The last three reactors were inoculated with a mixture of fecal microbiota from a healthy adult volunteer, who did not suffer from gastrointestinal diseases or take antibiotics in the previous 6 months according to previous classic studies (Yu et al., 2016; Wang et al., 2018). And the differences between individuals may be alleviated by same culture condition (Van den Abbeele et al., 2013). The study was approved by the Biomedical Ethics Committees of Nankai University. The participant has given written informed consent to understand the study purpose, procedures, risks, benefits, and rights. The details of the SHIME system and the startup process are summarized in the Supplementary Material.
As shown in Supplementary Figure S1, during the first 2 weeks (0 to 14 days) of the experiment, control nutritional medium was added into the reactors to stabilize the microbial community. After this period, the SHIME was sequentially exposed to nutritious medium (15 to 21 days), and nutritious medium + 600 mg day–1 AMX (22 to 28 days), each time system was maintained for 1-week. Then a nutritious medium was added and observed for another 2-weeks (29 to 42 days). Samples were collected and analyzed at six time points of 14, 21, 24, 28, 35, and 42 days from the ascending colon, transverse colon and descending colon, respectively. Each of the 18 samples is a mixture of three samples collected at specific time intervals in a day (Liu et al., 2020), and AMX was added after samples (C-A2, C-T2, and C-D2) had been collected in 22 days and discontinued after samples (AMX-A2, AMX-T2, and AMX-D2) had been collected in 29 days. Therefore, these samples were divided into three groups: before AMX administration during 0 to 21 days (control group), AMX exposure during 22 to 28 days (AMX treatment group) and after the AMX discontinuance during 29 to 42 days (recovery group), which were according to previous classic studies (Yu et al., 2016; Wang et al., 2018). Specifically, during each sampling time, three sterilized centrifuge tubes (50 ml) were used to collect the samples (10 ml) flow out from each colon vessel of the SHIME, respectively, and these samples were initially stored at −4°C (Van de Wiele et al., 2010; Yu et al., 2016). After all the samples were collected in each sampling day, three samples from the same vessel were mixed into one sterilized centrifuge tube (50 ml), which was operated in a super-clean bench. Then the samples were centrifuged at 10,400 g for 10 min, and the separated supernatants and pellets were stored at −80°C for further analyses.
16S rRNA Gene Sequencing and Analysis
Total DNA was extracted from the samples using the E.Z.N.A. stool DNA Kit (Omega, United States) according to the manufacturer’s protocols. The V3–V4 region of the bacterial 16S rRNA gene was amplified by polymerase chain reaction (PCR) using primers 341F 5′-CCTAYGGGRBGCASCAG-3′ and 806R 5′-GGACTACNNGGGTATCTAAT-3′ (95°C for 2 min, followed by 27 cycles at 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s, and a final extension at 72°C for 5 min). Amplicons were paired-end sequenced (PE250) on Illumina Hiseq2500 platform. In the end, a total of 651,502 tags were obtained. The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database under the accession number of SRR9330193–SRR9330210. The details of bacterial DNA extraction and PCR amplification of the 16S rRNA gene are described in the Supplementary Material.
Raw Illumina fastq files were de-multiplexed, quality-filtered, and analyzed using Quantitative Insights Into Microbial Ecology (QIIME) (Caporaso et al., 2010). The 16S rRNA gene sequences were classified taxonomically using the Ribosomal Database Project (RDP), and classifier 2.0.1 (Wang et al., 2007). The AMX exposure on alpha diversity was further measured by the taxon richness (Chao1 index), evenness (Simpson index), and diversity index (Shannon index) using all recommended samples. Besides, beta diversity of microbiota communities at baseline and after antibiotics were portrayed by nonmetric multidimensional scaling (NMDS) and principal coordinate analysis (PCoA) of weighted and unweighted UniFrac distances. Linear discriminant analysis effect size (LEfSe) analysis was performed to determine the bacterial taxa that significantly differed between the control and AMX exposure group using Galaxy application tool (Segata et al., 2011). Functional predictions of microbial community were performed to visualize the distribution of functional pathway genes in the three parts of the colon with different treatments using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) (Langille et al., 2013). The accuracy of PICRUSt for the detection of these more challenging functional groups was good (minimum accuracy = 0.82), suggesting that their inference of gene abundance across various types of functions was reliable, and PICRUSt predictions had high agreement with metagenome sample abundances across all body sites (Spearman r = 0.82, P < 0.001). PICRUSt has been successfully manipulated in many previous research for predicting microbial function of human intestinal microbiota (McHardy et al., 2013; Bunyavanich et al., 2016; Labus et al., 2017; Puri et al., 2018). The details of taxonomical classification, LEfSe analysis, and functional predictions are described in the Supplementary Material.
High-Throughput Quantitative PCR (HT-qPCR) and Analysis
High-throughput quantitative PCR reactions were performed to visualize the variation of ARGs during the treatment using Wafergen SmartChip Real-time PCR system, conducted by Anhui MicroAnaly Gene Technologies Co., Ltd. (Anhui, China). A total of 108 primer sets were used (Excel S1), including 102 primer sets to target the almost all major classes of ARGs found in the Chinese human gut microbiota (Hu et al., 2013), five mobile genetic elements (MGEs), and one 16S rRNA gene. The results were analyzed with SmartChip qPCR Software by excluding the wells with multiple melting peaks or amplification efficiency beyond the range (90–110%). Then data were screened with the conditions that a threshold cycle (CT) must be <31 and positive samples should have three replicates simultaneously. The details of HT-qPCR analysis are described in the Supplementary Material.
Data Analysis
All the results were expressed as mean values and standard deviations. The statistical analysis was performed with SPSS 17.0 software (SPSS Inc., Chicago, IL, United States). The T-test was conducted to compare the differences between the groups, and all the statistical tests were two-tailed. The statistical significance was set at three different levels (∗P < 0.05, ∗∗P < 0.01, and ∗∗∗P < 0.001). Spearman test, Mantel test, and Procrustes test for correlation analysis between the microbiota and the functional pathway genes or ARGs were performed in R with the vegan package. Correlations between the pairs of variables were considered to be significant at R > 0.6, and P values were <0.05. The Gephi (V 0.9.1) software was used to visualize the bipartite network graphs using the Force Atlas algorithm.
Results
AMX Exposure Increased Microbiota Diversity
The effects of 600 mg day–1 AMX treatment on the gut bacterial community were revealed by the 16S rRNA gene sequencing of fecal samples collected from three different vessels designated as ascending, transverse, and descending colon. The vessels were set up for three groups: before AMX administration during 0 to 21 days (control group), AMX exposure during 22 to 28 days (AMX treatment group) and after the AMX discontinuance during 29 to 42 days (recovery group). The alpha diversity of the fecal microbiota was assessed in each group. The taxon richness (Chao1 index), evenness (Simpson index), and diversity index (Shannon index) of the three groups are shown in Figure 1. As compared with the control group, a substantially rising in the evenness (Simpson, P = 0.067, T-test) and diversity (Shannon, P = 0.061, T-test) was observed in the AMX treatment groups; however, no change occurred in the microbiota richness (Chao1, P = 0.564, T-test). Moreover, data displayed that an increasing of microbiota evenness and diversity caused by the AMX treatments was continued after two weeks of AMX discontinuance (Simpson, P = 0.044, T-test; Shannon, P = 0.028, T-test).
Figure 1. Effects of AMX treatment on gut microbiota alpha diversity of the control (green), AMX (red), and recovery (blue) groups. The Chao1 index (A) was used to calculate the community richness, Simpson index (B) was used to calculate the community evenness, and Shannon index (C) was used to calculate the community diversity within each of the three groups (control, AMX, and recovery). Statistical significance between each of the three groups were analyzed using the T test at a significance level of 0.05.
In addition, the beta diversity of the microbial communities and weighted UniFrac distance was also affected by AMX treatment. As shown in Figure 2 and Supplementary Figure S2, the samples collected after the AMX treatment were differed from the control group, because the two groups clustered far away from each other in both UniFrac NMDS and PCoA analyses. The average weighted UniFrac distance between the AMX treatment group and control group (AMX vs control) was significantly higher (P < 0.001, T-test) than that within the control group (control vs control). Similarly, the beta diversity results showed that the gut microbial composition remained comparable after 2-week of AMX discontinuance because these samples were clustered together with the AMX treatment groups observed in both UniFrac NMDS and PCoA analyses. The average weighted UniFrac distance between the recovery group and the control group (recovery vs control) was also significantly higher (P < 0.001, T-test) than the control group (control vs control).
Figure 2. Effects of AMX treatment on gut microbiota beta diversity of the control (green), AMX (red), and recovery (blue) groups. Unweighted (A) and weighted (B) PCoA, and unweighted (C) and weighted (D) NMDS of UniFrac distances of samples in three groups (control, AMX, and recovery).
AMX Exposure Changed Microbiota Community Composition
The bacterial community compositions and their shifts at the phylum and genus levels after AMX exposure are shown in Figure 3 and Supplementary Figure S3. Based on the 16S rRNA gene analysis, the taxonomic assignment was related to most four dominant phyla such as Bacteroidetes, Firmicutes, Proteobacteria, and Synergistetes, which accounted for 96.5 to 97.9% of the total community (Figure 3A). However, a noticeable increase in abundance of Bacteroidetes (from 16.3 to 19.9%) and Synergistetes (from 8.6 to 13.2%), and decrease of Proteobacteria (from 59.6 to 50.7%) were seen after AMX treatment. The shifted phenomena (Bacteroidete from 13.3 to 23.2%; Proteobacteria from 68.8 to 57.8%) were more prominent in the ascending colon than the transverse and descending colons (Supplementary Figure S3A). A significant alteration in the communities was also observed at the phylum level (Bacteroidetes from 19.9 to 30.5%; Proteobacteria from 50.7 to 39.4%) after 2-week of AMX discontinued.
Figure 3. Composition shifts of gut microbial community at phylum (A) and genus level (B) between each of the three groups (control, AMX, and recovery).
At the genus level, the Bacteroides, Klebsiella, Megasphaera, and Pseudomonas were predominant (Figure 3B). The antibiotic-treated subjects were shown to be substantially overgrown by Bacteroides (from 9.8 to 15.0%), Klebsiella (from 1.5 to 18.1%) and Pyramidobacter (from 0.8 to 10.0%), and declined in the percentage of Cloacibacillus (from 6.2 to 2.5%) and Parabacteroides (from 4.5 to 1.7%). Similarly, the changes of Bacteroides (from 8.4 to 20.5%) and Klebsiella (from 1.6 to 32.9%) were more evident in the ascending colon (Supplementary Figure S3B). As shown in Figure 3B, during the recovery period, it observed that the gut microbiota was not fully recuperated. Only, Bacteroides was inclined by abundance from 15.0 to 23.4%; however, all others were less retrieved (Klebsiella from 18.1% to 10.9%; Pyramidobacte from 10.0 to 8.5%; Cloacibacillus from 2.5 to 4.3%; Parabacteroides from 1.7 to 3.7%).
The LEfSe showed that the AMX exposure caused an obvious decrease in several taxa, including the members of Parabacteroides [linear discriminant analysis (LDA) score = 4.18], Bifidobacterium (LDA score = 4.06), and Phascolarctobacterium (LDA score = 3.95) (Figure 4A). The variation was accompanied by significant increases in the relative abundances of Klebsiella (LDA score = 5.26), and Bacteroides uniformis (LDA score = 4.75). However, after 2-week of AMX discontinuance, the decrease of Bifidobacterium (LDA score = 3.88), and increase of Klebsiella (LDA score = 4.67) and Bacteroides uniformis (LDA score = 4.65) were still recognizable as compare with the control group (Figure 4B).
Figure 4. LDA score and cladogram of LEfSe comparison analysis between the control and AMX groups (A), and between the control and recovery groups (B). The red, green, or blue shading depicts bacterial taxa that were significantly higher in either the control, AMX or recovery groups, as indicated. Selection of discriminative taxa between the control and AMX groups or between the control and recovery groups were based on an LDA score cutoff of 3.0, and differences in the relative abundances of taxa were statistically determined based on a Mann–Whitney test at a significance level of 0.05.
AMX Exposure Increased Functional Pathway Genes
The hierarchy cluster heatmap analysis using metagenomic 16S rRNA gene sequencing was predicted by PICRUSt. Results showed that functional pathway genes, which included cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems were more abundant in the antibiotic exposure group than in the control group that represented by the ascending colon (Figure 5). For instance, as compared to the AMX-free sample C-A2, collected from the ascending colon before AMX administration, the gene numbers of human diseases related-pathways were 2.4–3.8 times enriched in the AMX exposed sample AMX-A2 that obtained after AMX treatment for 7 days. Similarly, the numbers of genes respect to functional pathways, for examples, membrane transport, transcription, xenobiotics biodegradation, metabolism, and excretory system were nearly 4.6–5.3 folds enriched in the AMX-A2 sample than C-A2. Moreover, these pathways were still maintained at a higher level (approximately 2.2–2.5 times enhanced than the control group) after the recovery phase.
Figure 5. Heatmap of human disease-related pathways in the three parts of colon within three groups and the network analysis revealing the co-occurrence patterns between microbial taxa abundance and functional pathway gene number. Heatmap colors reflect relative abundance from low (blue) to high (red). The nodes in network were colored according to types of functional pathway genes and microbial genera that increased (red) or decreased (green) after AMX exposure, and the edges were colored according to positive (red) or negative (blue) correlation. A connection represents strong and significant (P value < 0.05, R > 0.6) correlation. The size of each node is proportional to the number of connections, that is, degree.
To investigate whether the OTUs correlated with the functional pathway genes, the Mantel test and Procrustes analysis were used. The results showed that OTUs from the ascending, transverse, and descending colons were moderately correlated with the numbers of functional pathway genes (Mantel test, R = 0.743, P < 0.001). The Procrustes analysis indicated that OTUs obtained from the 16S rRNA gene data and the functional genes could be clustered according to the type of sample, which further exhibited a goodness-of-fit test (R = 0.669, P < 0.001, and 9,999 permutations) based on the Bray–Curtis dissimilarity metrics (Supplementary Figure S4A). As shown in Figure 5, the results demonstrated the co-occurrence patterns of significantly shifted microbial taxa and functional pathway genes during AMX-treatment. It can be seen that the abundances of most imminent bacteria after AMX exposure such as Bacteroides uniformis and Klebsiella were positively associated with the gene numbers of functional pathways; however, the suppressed bacteria like Phascolarctobacterium and Parabacteroides were negatively associated (Figure 5). Specifically, the correlation coefficients of Klebsiella with the gene numbers of the digestive system, immune system, and metabolic diseases were about 0.98 (P < 0.001) and the correlation coefficients of Phascolarctobacterium abundance with gene numbers of signaling molecules and interaction, transport and catabolism, endocrine system, cell growth, and death were about −0.77 (P < 0.01).
AMX Exposure Increased the Abundance of Beta-Lactam and Tetracycline Resistance Genes
The relative abundances of ARGs such as beta-lactam and tetracycline resistance genes were substantially higher in the AMX exposure group as compared to the control group, while the multidrug-resistant ARGs and transposase were lower in the AMX exposure group (Figure 6). Notably, the relative log abundance of bl2_len, bl2b_tem1, and bl2be_ctxm (aminoglycoside) were about 1.6–2.4 log units higher after AMX treatment (AMX-A2) than in control (C-A2). Similarly, the tetb, tetc, and tetr (tetracycline) were about 1.4–2.5 log units higher than the control group. Besides, the reductions of relative log abundance in the ARGs such as mdte and tolc (multidrug) were about 1.2 log units after AMX treatment (AMX-A2) than control (C-A2), and Tn22 (transposase) was 3.0 log units. Similarly, these ARGs were unable to return at the baseline level following the recovery procedure with about 0.6 to 3.2 log units change of relative abundance compared with the control group.
Figure 6. Heatmap of antibiotic resistance genes (ARGs) in the three parts of colon within three groups and the network analysis revealing the co-occurrence patterns between abundance of microbial taxa and that of ARG subtypes. Heatmap colors reflect relative abundance from low (blue) to high (red). The nodes in network were colored according to ARG types and microbial genera that increased (red) or decreased (green) after AMX exposure, and the edges were colored according to positive (red) or negative (blue) correlation. A connection represents strong and significant (P value < 0.05, R > 0.6) correlation. The size of each node is proportional to the number of connections, that is, degree.
To investigate whether the OTUs correlated with the resistome composition, Mantel test and Procrustes analysis were also performed to correlate the OTUs with the resistome using ascending, transverse and descending colons’ samples. Our results showed that OTUs were strongly correlated to the ARG profiles (Mantel test, R = 0.895, P < 0.001). The Procrustes analysis demonstrated that the bacterial OTUs and the ARGs in HT-qPCR data could be clustered by the type of sample and exhibited a goodness-of-fit-test (R = 0.941, P < 0.001, and 9,999 permutations) by Bray-Curtis dissimilarity metrics (Supplementary Figure S4B). The network analyses of co-occurrence patterns between the significantly changed microbial taxa after AMX treatment and ARG subtypes were shown in Figure 6. Interestingly, a very similar pattern of results was observed between the microbial taxa and functional pathway genes. The abundance of bacteria that significantly decreased after AMX treatment was positively associated with multidrug ARGs, while the abundance of increased bacteria was positively associated with beta-lactam resistance genes. For example, the moderate correlation coefficient of Bifidobacterium longum with acra abundance was 0.70 (P < 0.05), and for Phascolarctobacterium with acrb and mdte were 0.67 (P < 0.05). The strong correlations were found in the abundance of Klebsiella with the abundances of bl2_len, bl2b_tem1, and bl2be_shv2 (R = 0.9, P < 0.001), and in Bacteroides uniformis with bl2_len and bl2be_shv2 (R = 0.8, P < 0.01).
Discussion
AMX Treatment Has a Tremendous Influence on the Overall Taxonomic Composition and Biodiversity
The SHIME model was stably operated in this study because the predominant phyla of Bacteriodetes, Proteobacteria, Synergistetes, and Firmicutes in the gut microbiome was previously demonstrated Yu et al. (2016). Firmicutes and Bacteroidetes are usually dominate in the microbiota of a healthy subject, however in the control microbiota of SHIME, Proteobacteria is majoritarian. For in vivo studies, highest percentage of Proteobacteria had also been observed in fecal samples from healthy human and animals (Gao et al., 2015; Li et al., 2015). Our previous in vitro study also discovered this phenomenon (Liu et al., 2020).
It has been reported that in vivo experiments of AMX caused by tiny effects on human gut microbiota composition (Vrieze et al., 2014; Reijnders et al., 2016). However, a definite shift in the phylum level in our research may attribute to the absence of disturbances from neurohumoral regulation, the individual differences, dietary habits, and physiological status using in vitro SHIME (Karl et al., 2018). The decrease of Proteobacteria shown in our study may also attribute to AMX is more effective for sensitive Gram-negative bacteria that belonging to Proteobacteria by interfering with the synthesis of bacterial cell wall peptidoglycan layer (Zapata and Quagliarello, 2015). And Bacteroides uniformis belonging to Bacteroidetes may be resistant to AMX to make them survive, because of the strong correlations were found in their abundance with bl2_len and bl2be_shv2 (R = 0.8, P < 0.01). Some studies even showed the opposite phenomena like increasing of Proteobacteria and decreasing of Bacteroidetes, which might be caused by the combined effects with other antibiotics such as clarithromycin, fosfomycin, and metronidazole (Oh et al., 2016; Ishikawa et al., 2017). However, at the genus level, the blooming of intestinal microbiota such as Bacteroides uniformis and Klebsiella and decrease of Bifidobacterium occurred following AMX administration. These data were further supported by earlier reports (Ladirat et al., 2014a; Oh et al., 2016; Marzorati et al., 2017), while an increase of Parabacteroides was found in another study (Kabbani et al., 2017).
Amoxicillin treatment affected a vast number of sensitive intestinal bacterial species that directly engrossed their corresponding functions. As a result, the functionally redundant resistant intestinal bacterial species were initially present at low levels and became more abundant. Thus, to compensate this loss of function to maintain community function, and species with closer evolutionary relationship usually have more similar functions (Moya and Ferrer, 2016). For example, Van den Abbeele et al. (2012) have found that during antibiotic treatment, the abundance of antibiotic-resistant Pediococcus acidilactici was enlarged, while the amount of sensitive Lactobacillus mucosae was declined. Similarly, in this study, a reduction of sensitive Parabacteroides might be replaced by resistant Bacteroides species. Although the metabolic function of the intestinal microbiota could recover quickly, the unrestored composition of the intestinal microbiota was supported by increasing ARB, which may pose a significant therapeutic challenge and that needs more attention.
As demonstrated by several studies that administration of antibiotics is significantly associated with decrease in microbial community diversity and richness, our results that AMX exposure increased the microbial diversity seems counterintuitive (Dethlefsen et al., 2008). The study by Vrieze et al. (2014) discovered a trend toward increased taxonomic diversity after AMX exposure, which was in accordance with our study. However, in some other studies no significant effect on microbiome diversity was found, which were opposed to our study (Pallav et al., 2014; Zaura et al., 2015; Reijnders et al., 2016). For this phenomenon, the possible reason was that AMX resistant intestinal bacterial species could bloom and compensate for or even surpassed the loss of AMX sensitive species. Intestinal microbiota with higher biodiversity was usually more resistant to the perturbation and colonization by pathogens, which might help in maintaining the intestinal homeostasis (Li et al., 2016). However, the improved biodiversity by AMX exposure that discovered in our study is mainly caused by increased abundance of AMX resistant intestinal bacterial species, which included opportunistic pathogen of the human intestine. This phenomenon may not be considered healthy for people because the increasing of opportunistic pathogens may be related to some human diseases that may cause the transmission and diffusion of ARGs.
Pathogen Contribute to the Increased Human Disease Pathway Genes and Beta-Lactam Resistance Genes
Some research articles have not mentioned the significant alteration of metabolic functions of the gut microbiota after AMX treatment (Vrieze et al., 2014; Reijnders et al., 2016), while the disruption of the metabolic activity of microbiota (increased in succinate and monosaccharide and oligosaccharide levels) in the fecal samples were found elsewhere (Ladirat et al., 2014a). Moreover, this study found that Klebsiella was positively associated with the gene numbers of functional pathways including cancer, metabolic and immune diseases, and other human diseases, which have not been reported yet in other research articles. K. pneumoniae is a typical Gram-negative bacterial pathogen, which frequently colonizes in the human gut and processes to infection diseases (Pena et al., 1997; Wiener-Well et al., 2010; Gorrie et al., 2017). The increase of the relative abundance of intestinal Klebsiella genus has been reported to associate with diverse human diseases such as pneumonia, inflammation, Crohn’s disease, colitis, cystitis, liver abscess, and wound infections (Schneditz et al., 2014; Atarashi et al., 2017). All above findings suggested that Klebsiella might contribute to increase the pathway genes, including cancer, immune system diseases, infectious diseases, metabolic diseases, and neurodegenerative diseases.
It is also known from the literature that frequent use of antibiotics such as AMX may cause the transmission and proliferation of ARB and ARGs (Canton and Morosini, 2011; Zaura et al., 2015; Lenart-Boron et al., 2016). This study further demonstrated that the beta-lactamase resistance genes and the potential co-selected ARGs were increased after AMX exposure. The co-occurrence patterns between abundance of microbial taxa and that of ARG subtypes (Figure 6) showed that Bacteroides uniformis and Klebsiella might contributed to the increasing of ARGs including beta-lactam, tetracycline, and multidrug resistance. Therefore, treatment with AMX increases the number of resistance genes, even those not directly related to the antibiotic ingested. Bacteroides uniformis and Klebsiella are ubiquitous in the environment and are commensal in the human gut, notably the increased of their population, especially for Klebsiella is not usually considered as healthy. Moreover, resistance to beta-lactams via beta-lactamase production has mostly been described in Klebsiella spp. (Madhi et al., 2013; Marco et al., 2017). Bifidobacterium and Parabacteroides are most common or novel probiotics (Pinto-Sanchez et al., 2017; Wu et al., 2019), and the decreasing of these probiotics may cause dysbiosis of gut microbiota, which also leads to human health problem (Wischmeyer et al., 2016; George Kerry et al., 2018).
AMX Caused “SHIME-Compartment” Different and Permanent Alterations
In this study, Klebsiella and Bacteroides uniformis were found to be more significantly enriched in the ascending colon than other two colonic-regions, which was in keeping with the previous reports of the alteration in intestinal microbiota by exposure of antibiotics mixture and other toxic compounds like chlorpyrifos and arsenic (Reygner et al., 2016; Yu et al., 2016; Marzorati et al., 2017). These studies revealed the “SHIME-compartment” specific effects that could be due to the inconsistent biodegradation of multiple compounds, gut microbiome community, and pH in different colon regions (Reygner et al., 2016; Yu et al., 2016; Marzorati et al., 2017). Moreover, a large number of functional pathways related genes was also shown to be increased in the ascending colon after AMX exposure, which further clarified that the primary effect observed at the level of the microbiota could also be identified at the genomic and metabolic levels (Garcia-Villalba et al., 2017; Wang et al., 2018). The intestinal microbiota is a key “organ” for the individual’s health, and its response to AMX is started from the proximal colon and hence observed more obvious shifts there. This phenomenon suggested that more studies should focus on the proximal colon; however, it may be challenging to study by in vivo experiments that usually analysis the feces standing for the distal intestinal microbiota.
Furthermore, these AMX effects were still evident for at least 2 weeks after the AMX discontinuance, although the resilient tendency of microbial composition, functional pathway genes and ARGs were also observed. Pallav et al. (2014) revealed that AMX mediated opportunistic pathogen such as Escherichia/Shigella persisted up to 42 days after the interruption of antibiotic therapy. Similar phenomenon was also reported by Zwittink et al. (2018), whose work primarily discovered a 5 days treatment with combined of AMX and ceftazidime allowed Enterococcus to thrive and remain dominant up to 2 weeks, while the abundance of Bifidobacterium remained decreased till postnatal of 6 weeks after antibiotic treatment discontinuation. Therefore, our results demonstrated a serious negative side-effect of AMX, which could be persistent and different in the specific colon region and should be considered as an essential aspect of the risk assessment for AMX prescription.
PERSPECTIVES
Typical antibiotics have been detected in the collective gut of the Chinese population in our previous research, which provided a reference for this study of the effects of AMX on human gut microbiota (Wang et al., 2020). This study provides the originality of being able to control the effect of antibiotics in the different parts of the colon, evaluating in situ the alterations of the microbiota. Although SHIME has been successfully manipulated according to previous classic studies (Van den Abbeele et al., 2013; Yu et al., 2016; Wang et al., 2018), there were still three limitations of this study. First, AMX could not be absorbed by SHIME on the small intestine nor interact with the upper microbiota, which make the burden reaching the colon in the SHIME is much higher comparing to practical exposure to human colon. Considering that about 77 to 93% AMX is absorbed by the gastrointestinal tract and their interaction with the upper microbiota, the putative exposure dosage of AMX in colon is about 10% of the original oral dosage (Spyker et al., 1977; Arancibia et al., 1980; Huttner et al., 2019). This model is better to include the interaction with the human cell to get a closer result to in vivo, which is also a very good suggestion. Therefore, we would like to use 60 mg day–1 AMX (10% of this study) and include interaction with the human cell in our future work, with which the results would be close to in vivo. Second, this study lack of biological replicates and appropriate controls, which makes it difficult to interpret the results in relation to the microbial community variations and changes observed. Therefore, we would collect biological replicates and run a parallel experiment with blank medium at T0 and other timelines in our future study to negate any background noise and make the research scientifically reasonable. Another limitation is that the composition of the microbiota from the donor has not been analyzed before introducing it into the simulator. The composition of the initial microbiota from the donor may affect the stabilized microbiota in SHIME. Therefore, we would analyze the composition of the microbiota from the donor before introducing it into the simulator and compare the stabilized microbiota in SHIME that using multiple donors either by mixing their fecal samples or adding individually to explore the effects of initial microbiota composition in our future work, which would make this study more interesting. Although our work has these limitations, our findings would be valuable for directing future work. The findings in this study suggested several numbers of opportunities for additional study. One avenue is to expand the analysis to incorporate multi-omics approaches of the metagenome, metatranscriptome, and metabolome. Functional and metabolite analysis by multi-omics approaches would refine the results of predicting microbial function. Including interaction with the human cell in the fermentation vessels or manipulate in vivo study is also a good avenue to confirm the actual AMX’s effects on cancers, metabolic and immune diseases, and other human diseases that found in our in vitro experiments. On the other hand, the studies of virome and fungome may exert the substantial influences on the intestinal microbiome. Thus, an expanded analysis of other microorganisms will also be necessary. It is of interest to investigate the impacts of different kinds of antibiotics and a mixture of them on the gut microbiota. Moreover, research on the restoration effects of some prebiotics, probiotics, and synbiotics during and after antibiotic therapy is promising, which may help to discover several clinical strategies and restore side effects caused by antibiotics. As the increasing of opportunistic pathogens, functional genes and ARGs by antibiotics exposure may pose a significant therapeutic challenge, it should be more critical to take some efficient measures to reduce or even eliminate the effects caused by antibiotic treatment.
Conclusion
Exposure to AMX had significantly altered the overall taxonomic composition of the gut microbiota with increasing taxonomic richness, functional pathway genes, and beta-lactam resistance genes. The changes were “SHIME-compartment” different and more substantial effect was observed in the ascending colon. The shifted human gut microbiota could not be restored after 2 weeks’ of AMX discontinued. Importantly, most of the functional pathway genes and quantified beta-lactam resistance genes are positively associated with bacteria that increased after AMX exposure. Our results may open up new perspectives for the assessing of direct effects by antibiotics on the intestinal microbiota – a key “organ” in individual health. The results demonstrated the negative side-effects of AMX and should be considered for AMX prescription.
Data Availability Statement
The datasets generated for this study can be found in the NCBI SRA database, SRR9330193–SRR9330210.
Ethics Statement
The study was approved by the Biomedical Ethics Committees of Nankai University. The participant has given written informed consent to confirm they understand the study purpose, procedures, risks, benefits, and rights.
Author Contributions
LL carried out the laboratory work, analyzed the data, and wrote the manuscript. QW provided suggestions in the manuscript preparation and revised this manuscript. RD edited the language and improved the clarity of this manuscript. HL, SW, HQ, and JY carried out the laboratory work. YX provided funding support. DM and YL guided the laboratory work and revised this manuscript. All authors read and approved the final manuscript.
Funding
This study was sponsored by the Key projects of the National Natural Science Foundation of China (41831287), China National Funds for Distinguished Young Scientists (41525013), National Natural Science Foundation of China (31870351, 41703088, 31670509, and 21607016), Key projects of Research and Development of Hebei Province (19273707D), and Key projects of Tianjin Natural Science Foundation (19JCZDJC40800).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors express the sincerest thanks to Professor Bing Wu at Nanjing University for the guidance of constructing SHIME model, and Guangzhou Gene Denovo Co., Ltd. and Anhui MicroAnaly Gene Technologies Co., Ltd. for genome sequencing and analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01213/full#supplementary-material
References
Arancibia, A., Guttmann, J., Gonzalez, G., and Gonzalez, C. (1980). Absorption and disposition kinetics of amoxicillin in normal human subjects. Antimicrob. Agents Chemother. 17, 199–202. doi: 10.1128/aac.17.2.199
Atarashi, K., Suda, W., Luo, C., Kawaguchi, T., Motoo, I., Narushima, S., et al. (2017). Ectopic colonization of oral bacteria in the intestine drives TH1 cell induction and inflammation. Science 358, 359–365. doi: 10.1126/science.aan4526
Barc, M. C., Charrin-Sarnel, C., Rochet, V., Bourlioux, F., Sandre, C., Boureau, H., et al. (2008). Molecular analysis of the digestive microbiota in a gnotobiotic mouse model during antibiotic treatment: influence of Saccharomyces boulardii. Anaerobe 14, 229–233. doi: 10.1016/j.anaerobe.2008.04.003
Barraud, O., Peyre, M., Couve-Deacon, E., Chainier, D., Bahans, C., Guigonis, V., et al. (2018). Antibiotic resistance acquisition in the first week of life. Front. Microbiol. 9:1467. doi: 10.3389/fmicb.2018.01467
Blaser, M. J. (2016). Antibiotic use and its consequences for the normal microbiome. Science 352, 544–545. doi: 10.1126/science.aad9358
Bunyavanich, S., Shen, N., Grishin, A., Wood, R., Burks, W., Dawson, P., et al. (2016). Early-life gut microbiome composition and milk allergy resolution. J. Allergy Clin. Immunol. 138, 1122–1130. doi: 10.1016/j.jaci.2016.03.041
Bussche, J. V., Marzorati, M., Laukens, D., and Vanhaecke, L. (2015). Validated high resolution mass spectrometry-based approach for metabolomic fingerprinting of the human gut phenotype. Anal. Chem. 87, 10927–10934. doi: 10.1021/acs.analchem.5b02688
Cabreiro, F., Au, C., Leung, K. Y., Vergara-Irigaray, N., Cocheme, H. M., Noori, T., et al. (2013). Metformin retards aging in C. elegans by altering microbial folate and methionine metabolism. Cell 153, 228–239. doi: 10.1016/j.cell.2013.02.035
Canton, R., and Morosini, M. I. (2011). Emergence and spread of antibiotic resistance following exposure to antibiotics. FEMS Microbiol. Rev. 35, 977–991. doi: 10.1111/j.1574-6976.2011.00295.x
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chow, J., Lee, S. M., Shen, Y., Khosravi, A., and Mazmanian, S. K. (2010). Host-bacterial symbiosis in health and disease. Adv. Immunol. 107, 243–274. doi: 10.1016/S0065-2776(10)07001-X
Collignon, A., Bouttier, S., Lambert, S., Hoys, S., and Barc, M.-C. (2008). Gnotobiotic mouse model to study the effect of antibiotics on human faecal microbiota. Microb. Ecol. Health Dis. 20, 204–206. doi: 10.1080/08910600802408137
Collignon, A., Sandre, C., and Barc, M. C. (2010). Saccharomyces boulardii modulates dendritic cell properties and intestinal microbiota disruption after antibiotic treatment. Gastroenterol. Clin. Biol. 34, S71–S78. doi: 10.1016/S0399-8320(10)70024-5
Dethlefsen, L., Huse, S., Sogin, M. L., and Relman, D. A. (2008). The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 6:e280. doi: 10.1371/journal.pbio.0060280
Duque, A. L. R. F., Monteiro, M., Adorno, M. A. T., Sakamoto, I. K., and Sivieri, K. (2016). An exploratory study on the influence of orange juice on gut microbiota using a dynamic colonic model. Food Res. Int. 84, 160–169. doi: 10.1016/j.foodres.2016.03.028
El Hage, R., Hernandez-Sanabria, E., Calatayud Arroyo, M., Props, R., and Van de Wiele, T. (2019). Propionate-producing consortium restores antibiotic-induced dysbiosis in a dynamic in vitro model of the human intestinal microbial ecosystem. Front. Microbiol. 10:1206. doi: 10.3389/fmicb.2019.01206
Flandroy, L., Poutahidis, T., Berg, G., Clarke, G., Dao, M. C., Decaestecker, E., et al. (2018). The impact of human activities and lifestyles on the interlinked microbiota and health of humans and of ecosystems. Sci. Total Environ. 627, 1018–1038. doi: 10.1016/j.scitotenv.2018.01.288
Gao, Z., Guo, B., Gao, R., Zhu, Q., and Qin, H. (2015). Microbiota disbiosis is associated with colorectal cancer. Front. Microbiol. 6:20. doi: 10.3389/fmicb.2015.00020
Garcia-Villalba, R., Vissenaekens, H., Pitart, J., Romo-Vaquero, M., Espin, J. C., Grootaert, C., et al. (2017). Gastrointestinal simulation model TWIN-SHIME shows differences between human urolithin-metabotypes in gut microbiota composition, pomegranate polyphenol metabolism, and transport along the intestinal tract. J. Agric. Food Chem. 65, 5480–5493. doi: 10.1021/acs.jafc.7b02049
George Kerry, R., Patra, J. K., Gouda, S., Park, Y., Shin, H. S., and Das, G. (2018). Benefaction of probiotics for human health: a review. J. Food Drug Anal. 26, 927–939. doi: 10.1016/j.jfda.2018.01.002
Gorrie, C., Mirceta, M., Wick, R., Edwards, D., Thomson, N., Strugnell, R., et al. (2017). Gastrointestinal carriage is a major reservoir of Klebsiella pneumoniae infection in intensive care patients. Clin. Infect. Dis. 65, 208–215. doi: 10.1093/cid/cix270
Hu, Y., Yang, X., Qin, J., Lu, N., Cheng, G., Wu, N., et al. (2013). Metagenome-wide analysis of antibiotic resistance genes in a large cohort of human gut microbiota. Nat. Commun. 4, 2151–2249. doi: 10.1038/ncomms3151
Huttner, A., Bielicki, J., Clements, M. N., Frimodt-Møller, N., Muller, A. E., Paccaud, J. P., et al. (2019). Oral amoxicillin and amoxicillin–clavulanic acid: properties, indications and usage. Clin. Microbiol. Infect. doi: 10.1016/j.cmi.2019.11.028 [Epub ahead of print].
Ianiro, G., Tilg, H., and Gasbarrini, A. (2016). Antibiotics as deep modulators of gut microbiota: between good and evil. Gut 65, 1906–1915. doi: 10.1136/gutjnl-2016-312297
Ichim, T. E., Kesari, S., and Shafer, K. (2018). Protection from chemotherapy- and antibiotic-mediated dysbiosis of the gut microbiota by a probiotic with digestive enzymes supplement. Oncotarget 9, 30919–30935. doi: 10.18632/oncotarget.25778
Ishikawa, D., Sasaki, T., Osada, T., Kuwahara-Arai, K., Haga, K., Shibuya, T., et al. (2017). Changes in intestinal microbiota following combination therapy with fecal microbial transplantation and antibiotics for ulcerative colitis. Inflamm. Bowel Dis. 23, 116–125. doi: 10.1097/Mib.0000000000000975
Kabbani, T. A., Pallav, K., Dowd, S. E., Villafuerte-Galvez, J., Vanga, R. R., Castillo, N. E., et al. (2017). Prospective randomized controlled study on the effects of Saccharomyces boulardii CNCM I-745 and amoxicillin-clavulanate or the combination on the gut microbiota of healthy volunteers. Gut Microbes 8, 17–32. doi: 10.1080/19490976.2016.1267890
Karl, J. P., Hatch, A. M., Arcidiacono, S. M., Pearce, S. C., Pantoja-Feliciano, I. G., Doherty, L. A., et al. (2018). Effects of psychological, environmental and physical stressors on the gut microbiota. Front. Microbiol. 9:2013. doi: 10.3389/fmicb.2018.02013
Kemperman, R. A., Gross, G., Mondot, S., Possemiers, S., Marzorati, M., Van de Wiele, T., et al. (2013). Impact of polyphenols from black tea and red wine/grape juice on a gut model microbiome. Food Res. Int. 53, 659–669. doi: 10.1016/j.foodres.2013.01.034
Kho, Z. Y., and Lal, S. K. (2018). The human gut microbiome - a potential controller of wellness and disease. Front. Microbiol. 9:1835. doi: 10.3389/fmicb.2018.01835
Labus, J. S., Hollister, E. B., Jacobs, J., Kirbach, K., Oezguen, N., Gupta, A., et al. (2017). Differences in gut microbial composition correlate with regional brain volumes in irritable bowel syndrome. Microbiome 5:49. doi: 10.1186/s40168-017-0260-z
Ladirat, S. E., Schoterman, M. H. C., Rahaoui, H., Mars, M., Schuren, F. H. J., Gruppen, H., et al. (2014a). Exploring the effects of galacto-oligosaccharides on the gut microbiota of healthy adults receiving amoxicillin treatment. Br. J. Nutr. 112, 536–546. doi: 10.1017/S0007114514001135
Ladirat, S. E., Schuren, F. H. J., Schoterman, M. H. C., Nauta, A., Gruppen, H., and Schols, H. A. (2014b). Impact of galacto-oligosaccharides on the gut microbiota composition and metabolic activity upon antibiotic treatment during in vitro fermentation. FEMS Microbiol. Ecol. 87, 41–51. doi: 10.1111/1574-6941.12187
Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lenart-Boron, A., Wolny-Koladka, K., Stec, J., and Kasprowic, A. (2016). Phenotypic and molecular antibiotic resistance determination of airborne coagulase negative Staphylococcus spp. strains from healthcare facilities in Southern Poland. Microb. Drug Resist. 22, 515–522. doi: 10.1089/mdr.2015.0271
Li, D. T., Wang, P., Wang, P. P., Hu, X. S., and Chen, F. (2016). The gut microbiota: a treasure for human health. Biotechnol. Adv. 34, 1210–1224. doi: 10.1016/j.biotechadv.2016.08.003
Li, S., Zhang, C., Gu, Y., Chen, L., Ou, S., Wang, Y., et al. (2015). Lean rats gained more body weight than obese ones from a high-fibre diet. Br. J. Nutr. 114, 1188–1194. doi: 10.1017/S0007114515002858
Liu, L., Wang, Q., Wu, X., Qi, H., Das, R., Lin, H., et al. (2020). Vancomycin exposure caused opportunistic pathogens bloom in intestinal microbiome by simulator of the human intestinal microbial ecosystem (SHIME). Environ. Pollut. doi: 10.1016/j.envpol.2020.114399 [Epub ahead of print].
Madhi, F., Biscardi, S., Bingen, E., Jaby, O., Epaud, R., and Cohen, R. (2013). Combined relay therapy with oral cefixime and clavulanate for febrile urinary tract infection caused by extended-spectrum beta-lactamase-producing Escherichia coli. Pediatr. Infect. Dis. J. 32, 96–97. doi: 10.1097/INF.0b013e318271f369
Marco, R. H., Olmos, E. G., Breton-Martinez, J. R., Perez, L. G., Sanchez, B. C., Fujkova, J., et al. (2017). Community-acquired febrile urinary tract infection caused by extended-spectrum beta-lactamase-producing bacteria in hospitalised infants. Enferm. Infecc. Microbiol. Clin. 35, 287–292. doi: 10.1016/j.eimc.2016.01.012
Marzorati, M., Vilchez-Vargas, R., Bussche, J. V., Truchado, P., Jauregui, R., El Hage, R. A., et al. (2017). High-fiber and high-protein diets shape different gut microbial communities, which ecologically behave similarly under stress conditions, as shown in a gastrointestinal simulator. Mol. Nutr. Food Res. 61, 1600150. doi: 10.1002/mnfr.201600150
McHardy, I. H., Goudarzi, M., Tong, M., Ruegger, P. M., Schwager, E., Weger, J. R., et al. (2013). Integrative analysis of the microbiome and metabolome of the human intestinal mucosal surface reveals exquisite inter-relationships. Microbiome 1:17. doi: 10.1186/2049-2618-1-17
Moya, A., and Ferrer, M. (2016). Functional redundancy-induced stability of gut microbiota subjected to disturbance. Trends Microbiol. 24, 402–413. doi: 10.1016/j.tim.2016.02.002
Oh, B., Kim, B. S., Kim, J. W., Kim, J. S., Koh, S. J., Kim, B. G., et al. (2016). The effect of probiotics on gut microbiota during the Helicobacter pylori eradication: randomized controlled trial. Helicobacter 21, 165–174. doi: 10.1111/hel.12270
Pallav, K., Dowd, S. E., Villafuerte, J., Yang, X., Kabbani, T., Hansen, J., et al. (2014). Effects of polysaccharopeptide from Trametes versicolor and amoxicillin on the gut microbiome of healthy volunteers: a randomized clinical trial. Gut Microbes 5, 458–467. doi: 10.4161/gmic.29558
Pena, C., Pujol, M., Ricart, A., Ardanuy, C., Ayats, J., Linares, J., et al. (1997). Risk factors for faecal carriage of Klebsiella pneumoniae producing extended spectrum beta-lactamase (ESBL-KP) in the intensive care unit. J. Hosp. Infect. 35, 9–16. doi: 10.1016/s0195-6701(97)90163-8
Pinto-Sanchez, M. I., Hall, G. B., Ghajar, K., Nardelli, A., Bolino, C., Lau, J. T., et al. (2017). Probiotic Bifidobacterium longum NCC3001 reduces depression scores and alters brain activity: a pilot study in patients with irritable bowel syndrome. Gastroenterology 153, 448–459. doi: 10.1053/j.gastro.2017.05.003
Puri, P., Liangpunsakul, S., Christensen, J. E., Shah, V. H., Kamath, P. S., Gores, G. J., et al. (2018). The circulating microbiome signature and inferred functional metagenomics in alcoholic hepatitis. Hepatology 67, 1284–1302. doi: 10.1002/hep.29623
Qin, J. J., Li, R. Q., Raes, J., Arumugam, M., Burgdorf, K. S., Manichanh, C., et al. (2010). A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464, 59–65. doi: 10.1038/nature08821
Reijnders, D., Goossens, G. H., Hermes, G. D., Neis, E. P., van der Beek, C. M., Most, J., et al. (2016). Effects of gut microbiota manipulation by antibiotics on host metabolism in obese humans: a randomized double-blind placebo-controlled trial. Cell Metab. 24, 63–74. doi: 10.1016/j.cmet.2016.06.016
Reygner, J., Joly Condette, C., Bruneau, A., Delanaud, S., Rhazi, L., Depeint, F., et al. (2016). Changes in composition and function of human intestinal microbiota exposed to chlorpyrifos in oil as assessed by the SHIME® model. Int. J. Environ. Res. Public Health 13:1088. doi: 10.3390/ijerph13111088
Rock, C., Thom, K. A., Masnick, M., Johnson, J. K., Harris, A. D., and Morgan, D. J. (2014). Frequency of Klebsiella pneumoniae carbapenemase (KPC)-producing and non-KPC-producing Klebsiella species contamination of healthcare workers and the environment. Infect. Control Hosp. Epidemiol. 35, 426–429. doi: 10.1086/675598
Schneditz, G., Rentner, J., Roier, S., Pletz, J., Herzog, K. A., Bucker, R., et al. (2014). Enterotoxicity of a nonribosomal peptide causes antibiotic-associated colitis. Proc. Natl. Acad. Sci. U.S.A. 111, 13181–13186. doi: 10.1073/pnas.1403274111
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Selak, M., Riviere, A., Moens, F., Van den Abbeele, P., Geirnaert, A., Rogelj, I., et al. (2016). Inulin-type fructan fermentation by bifidobacteria depends on the strain rather than the species and region in the human intestine. Appl. Microbiol. Biotechnol. 100, 4097–4107. doi: 10.1007/s00253-016-7351-9
Spyker, D. A., Rugloski, R. J., Vann, R. L., and O’Brien, W. M. (1977). Pharmacokinetics of amoxicillin: dose dependence after intravenous, oral, and intramuscular administration. Antimicrob. Agents Chemother. 11, 132–141. doi: 10.1128/aac.11.1.132
Stecher, B., Maier, L., and Hardt, W. D. (2013). ‘Blooming’ in the gut: how dysbiosis might contribute to pathogen evolution. Nat. Rev. Microbiol. 11, 277–284. doi: 10.1038/nrmicro2989
Truchado, P., Van den Abbeele, P., Riviere, A., Possemiers, S., De Vuyst, L., and Van de Wiele, T. (2015). Bifidobacterium longum D2 enhances microbial degradation of long-chain arabinoxylans in an in vitro model of the proximal colon. Beneficial Microbes 6, 849–860. doi: 10.3920/Bm2015.0023
Van de Wiele, T., Gallawa, C. M., Kubachka, K. M., Creed, J. T., Basta, N., Dayton, E. A., et al. (2010). Arsenic metabolism by human gut microbiota upon in vitro digestion of contaminated soils. Environ. Health Perspect. 118, 1004–1009. doi: 10.1289/ehp.0901794
Van de Wiele, T., Van den Abbeele, P., Ossieur, W., Possemiers, S., and Marzorati, M. (2015). “The simulator of the human intestinal microbial ecosystem (SHIME®),” in The Impact of Food Bioactives on Health: In Vitro and Ex Vivo Models, eds K. Verhoeckx, P. Cotter, I. López-Expósito, C. Kleiveland, T. Lea, A. Mackie, et al. (Cham: Springer International Publishing), 305–317. doi: 10.1007/978-3-319-16104-4_27
Van den Abbeele, P., Belzer, C., Goossens, M., Kleerebezem, M., De Vos, W. M., Thas, O., et al. (2013). Butyrate-producing Clostridium cluster XIVa species specifically colonize mucins in an in vitro gut model. ISME J. 7, 949–961. doi: 10.1038/ismej.2012.158
Van den Abbeele, P., Roos, S., Eeckhaut, V., MacKenzie, D. A., Derde, M., Verstraete, W., et al. (2012). Incorporating a mucosal environment in a dynamic gut model results in a more representative colonization by Lactobacilli. Microb. Biotechnol. 5, 106–115. doi: 10.1111/j.1751-7915.2011.00308.x
Vrieze, A., Out, C., Fuentes, S., Jonker, L., Reuling, I., Kootte, R. S., et al. (2014). Impact of oral vancomycin on gut microbiota, bile acid metabolism, and insulin sensitivity. J. Hepatol. 60, 824–831. doi: 10.1016/j.jhep.2013.11.034
Wang, Q., Duan, Y. J., Wang, S. P., Wang, L. T., Hou, Z. L., Cui, Y. X., et al. (2020). Occurrence and distribution of clinical and veterinary antibiotics in the faeces of a Chinese population. J. Hazard. Mater. 383, 121129. doi: 10.1016/j.jhazmat.2019.121129
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/Aem.00062-07
Wang, Y., Rui, M., Nie, Y., and Lu, G. (2018). Influence of gastrointestinal tract on metabolism of bisphenol A as determined by in vitro simulated system. J. Hazard. Mater. 355, 111–118. doi: 10.1016/j.jhazmat.2018.05.011
Wiener-Well, Y., Rudensky, B., Yinnon, A. M., Kopuit, P., Schlesinger, Y., Broide, E., et al. (2010). Carriage rate of carbapenem-resistant Klebsiella pneumoniae in hospitalised patients during a national outbreak. J. Hosp. Infect. 74, 344–349. doi: 10.1016/j.jhin.2009.07.022
Wischmeyer, P. E., McDonald, D., and Knight, R. (2016). Role of the microbiome, probiotics, and ‘dysbiosis therapy’ in critical illness. Curr. Opin. Crit. Care 22, 347–353. doi: 10.1097/MCC.0000000000000321
Wu, T. R., Lin, C. S., Chang, C. J., Lin, T. L., Martel, J., Ko, Y. F., et al. (2019). Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis. Gut 68, 248–262. doi: 10.1136/gutjnl-2017-315458
Yatsunenko, T., Rey, F. E., Manary, M. J., Trehan, I., Dominguez-Bello, M. G., Contreras, M., et al. (2012). Human gut microbiome viewed across age and geography. Nature 486, 222–227. doi: 10.1038/nature11053
Yu, H., Wu, B., Zhang, X., Liu, S., Yu, J., Cheng, S., et al. (2016). Arsenic metabolism and toxicity influenced by ferric iron in simulated gastrointestinal tract and the roles of gut microbiota. Environ. Sci. Technol. 50, 7189–7197. doi: 10.1021/acs.est.6b01533
Zapata, H. J., and Quagliarello, V. J. (2015). The microbiota and microbiome in aging: potential implications in health and age-related diseases. J. Am. Geriatr. Soc. 63, 776–781. doi: 10.1111/jgs.13310
Zaura, E., Brandt, B. W., Teixeira de Mattos, M. J., Buijs, M. J., Caspers, M. P., Rashid, M. U., et al. (2015). Same exposure but two radically different responses to antibiotics: resilience of the salivary microbiome versus long-term microbial shifts in feces. mBio 6:e01693-15. doi: 10.1128/mBio.01693-15
Zhang, Q., Ying, G., Pan, C., Liu, Y., and Zhao, J. (2015). Comprehensive evaluation of antibiotics emission and fate in the river basins of China: source analysis, multimedia modeling, and linkage to bacterial resistance. Environ. Sci. Technol. 49, 6772–6782. doi: 10.1021/acs.est.5b00729
Zwittink, R. D., Renes, I. B., van Lingen, R. A., van Zoeren-Grobben, D., Konstanti, P., Norbruis, O. F., et al. (2018). Association between duration of intravenous antibiotic administration and early-life microbiota development in late-preterm infants. Eur. J. Clin. Microbiol. Infect. Dis. 37, 475–483. doi: 10.1007/s10096-018-3193-y
Keywords: amoxicillin, antibiotic resistance genes (ARGs), functional pathway genes, human intestinal microbiota, simulator of the human intestinal microbial ecosystem (SHIME)
Citation: Liu L, Wang Q, Lin H, Das R, Wang S, Qi H, Yang J, Xue Y, Mao D and Luo Y (2020) Amoxicillin Increased Functional Pathway Genes and Beta-Lactam Resistance Genes by Pathogens Bloomed in Intestinal Microbiota Using a Simulator of the Human Intestinal Microbial Ecosystem. Front. Microbiol. 11:1213. doi: 10.3389/fmicb.2020.01213
Received: 22 October 2019; Accepted: 12 May 2020;
Published: 04 June 2020.
Edited by:
Leonardo Neves de Andrade, University of São Paulo, BrazilReviewed by:
Rosa Del Campo, Ramón y Cajal Institute for Health Research, SpainPallavi Singh, Northern Illinois University, United States
Copyright © 2020 Liu, Wang, Lin, Das, Wang, Qi, Yang, Xue, Mao and Luo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daqing Mao, maodq@nankai.edu.cn
†These authors have contributed equally to this work