Development of an in vitro Model of Human Gut Microbiota for Screening the Reciprocal Interactions With Antibiotics, Drugs, and Xenobiotics

Altering the gut microbiota can negatively affect human health. Efforts may be sustained to predict the intended or unintended effects of molecules not naturally produced or expected to be present within the organism on the gut microbiota. Here, culture-dependent and DNA-based approaches were combined to UHPLC-MS/MS analyses in order to investigate the reciprocal interactions between a constructed Human Gut Microbiota Model (HGMM) and molecules including antibiotics, drugs, and xenobiotics. Our HGMM was composed of strains from the five phyla commonly described in human gut microbiota and belonging to Firmicutes, Bacteroidetes, Proteobacteria, Fusobacteria, and Actinobacteria. Relevantly, the bacterial diversity was conserved in our constructed human gut model through subcultures. Uneven richness distribution was revealed and the sensitivity of the HGMM was mainly affected by antibiotic exposure rather than by drugs or xenobiotics. Interestingly, the constructed model and the individual cultured strains respond with the same sensitivity to the different molecules. UHPLC-MS/MS analyses revealed the disappearance of some native molecules in the supernatants of the HGMM as well as in those of the individual strains. These results suggest that biotransformation of molecules occurred in the presence of our gut microbiota model and the coupled approaches performed on the individual cultures may emphasize new bacterial strains active in these metabolic processes. From this study, the new HGMM appears as a simple, fast, stable, and inexpensive model for screening the reciprocal interactions between the intestinal microbiota and molecules of interest.

Introduction

The human gut microbiota is defined by a set of microorganisms in close interaction with the human gut (Vrancken et al., 2019) that must be adapted to its environment changes (Ejtahed et al., 2018; Gosalbes et al., 2019) by using different strategies allowing the species that compose it to maintain themselves and to communicate with each other (Sung et al., 2017; Selber-Hnatiw et al., 2020). It is composed of thousands of strains of which nearly 50 are common to all individuals and considered as the core microbiota of the human intestinal tract (Tap et al., 2009; Turnbaugh et al., 2009; Huse et al., 2012). There is a close relationship of mutualism and reciprocal control of species in a balanced microbiota, leading to an efficient adaptive capacity (Yaffe and Relman, 2020). A disturbed balance can lead to a disorganization of the species making up this microbiota (called dysbiosis) which can lead, in the case of disruption of the human microbiota, to many pathologies such as type II diabetes, obesity, or even autism (Kump et al., 2013; Le Bastard et al., 2018; Zuo et al., 2018). Various factors can be involved in this imbalance such as lifestyle, genetic factors, infections, medication, and especially antibiotics (Clarke et al., 2019). The effect of xenobiotics on the bacterial composition of microbiota has become a major interest in understanding the unintended effects they can induce but also in the study of specific strains that may be impacted. Furthermore, it is already known that bacteria in the intestinal microbiota can interact with xenobiotics and metabolize them to activate, reactivate, or render them inactive (Haiser et al., 2014; Enright et al., 2016; Wilson and Nicholson, 2017).

Our knowledge of the diversity and function of the gut microbiota is highly dependent on technological developments. Recent advances in molecular techniques (DNA-based approaches such as high throughput sequencing) have made it possible to visualize the diversity, the function, and the dynamic of the gut microbial communities in intimate association with the lining of the gut (Audebert et al., 2016; Lau et al., 2016; Frankel et al., 2017; Panek et al., 2018). However, given the complexity of the interactions and the number of microbial species found in the different microbiota, it appears necessary to develop simplified models allowing the culture of certain model microcosms (von Martels et al., 2017; Silverman et al., 2018; Villa et al., 2020) and the study, through these models, of the activity and the becoming of exogenous molecules, developed by the pharmaceutical or agri-food industry, in contact with the human microbiota (Maier et al., 2018; Zimmermann et al., 2019).

The combination of in vitro and in vivo investigations in prebiotic and probiotic research is beneficial. However, each approach has its advantage and disadvantages. Animal models (mice, rats, and piglets) allow researchers to study host-microbe interactions in highly controlled environments (Kostic et al., 2013). They also provide direct access to colonic contents, making them more suitable for studying gut microbiota metabolic activity than studies using feces (Pham and Mohajeri, 2018). Animal models are expensive, labor-intensive, physiologically unrepresentative of human metabolisms (Cohen, 1995), distinct from the human microbiome, and therefore it remains a challenge to translate the obtained findings from animal models to conditions in humans (Abboud et al., 2021; De Boeck et al., 2021).

To overcome these challenges, in vitro models have been developed as effective tools for studying the human gut microbiota in a controlled environment. In vitro fermentation models are cheaper (de Carvalho et al., 2021), more reproducible, do not require elaborate ethical approval, and can be conducted in a shorter time and enable the cultivation of human gut microbiota under simulated physiological conditions (Tsang et al., 2021). The in vitro models could be divided into two groups (Pham and Mohajeri, 2018): (i) in vitro batch fermentation, used mainly for the screening of the pre-and probiotics; and (ii) in vitro continuous fermentation models which mimic the human gut microbiota and are used to study the mechanism of action of the pre- and probiotics. The study of the microbiota is very complex to set up and calls for heavy and expensive equipment that requires culture systems (Marzorati et al., 2014; Shah et al., 2016; von Martels et al., 2017; von Martels, 2019). In addition, these studies generally focus exclusively on the gastrointestinal model, effectively excluding the whole microbiota (Iwasaki et al., 2019; Nguyen et al., 2020).

One of the advantages of developing in vitro models is being able to control very precisely and reproducibly the experiment conditions (Daina and Zoete, 2016) as well as the analysis of certain mechanisms regulating the balance of the microbiota (Chelakkot et al., 2018), the impact of the use of various molecules on the microbiota (Walsh et al., 2018; Tsitko et al., 2019; Vich Vila et al., 2020), and could significantly facilitate and accelerate the early screening of molecules of interest (Nguyen et al., 1999; Tsitko et al., 2019). In this context, we have studied the reciprocal impact of different molecules from pharmaceutical, agri-food, and industrial purposes on an in vitro Human Gut Microbiota Model (HGMM) we designed. The construction, the evolution, and the stability of the microbial communities of the HGMM were illustrated. The high-throughput sequencing analysis of the V5-V7 region of the 16S rRNA gene was combined to the UHPLC-MS/MS to reveal and evaluate the impact of the molecules including antibiotics, drugs, and xenobiotics on the HGMM. In addition, this study consisted of studying the effects of the tested molecules on the growth of single strains (individual cultures). Conversely, the effect of the single strains as well as the HGMM on the tested molecules, in term of the quantity of the native molecule in the supernatant, was investigated.

Materials and Methods

Bacterial Strains, Cultivation Medium, Growth Conditions, and Molecules

A total of 39 strains was selected from the representative panel of human gut core bacteria species previously described in the literature (Maier et al., 2018; Johnson et al., 2019; Zimmermann et al., 2019) to construct a Human Gut Microbiota Model in order to study its response toward different molecules including antibiotics, drugs, and xenobiotics. The bacterial strains (Table 1), tested alone separately or mixed in a constructed HGMM, were purchased from ATCC (American Type Culture Collection), from the DSMZ (Deutsche Sammlung von Mikroorganismen und Zellkulturen), or obtained from the laboratory culture collection of the University Hospital Center (Poitiers, France). The culture of the pure strains and the constructed HGMM was carried out in sterile Modified Gifu Anaerobic Medium Broth (mGAM Broth, HyServe), allowing the growth of all the selected strains. The mGAM medium was supplemented with resazurin 0.1% (w/v) as an anaerobiosis indicator and then was anaerobically distributed in appropriate containers (Hungate tubes or penicillin flasks) under oxygen-free conditions using N₂/CO₂/H₂ anaerobic mixture (90%/5%/5%, v/v). The final pH of the medium, after sterilization by autoclaving at 115°C for 15 min, was 7.3. This pH was in agreement with those reported in intestines, important host sites for microbial metabolism (Fallingborg, 1999; Koziolek et al., 2015). All the strains were reanimated and stored at –80°C in glycerol/mGAM (50% v/v). Cultures and growth maintenance procedures for the HGMM were performed in anaerobic conditions in Hungate tubes as previously described (El Houari et al., 2017) using mGAM broth. The temperature and the time of incubation were 37°C and 48 h, respectively.

TABLE 1

Table 1. List of the bacterial strains selected from the literature to construct the human gut microbiota model.

Antibiotics, drugs, and xenobiotics were purchased from Sigma Aldrich, resuspended in mGAM broth or DMSO (Sigma Aldrich), and stored at –20°C for further use. The molecules and their respective concentrations in which were tested are described in Table 2. The spectrum of action of the tested antibiotics is reported in the Supplementary File Data 1. The concentrations of the tested molecules were chosen based the range of the screening concentration (below 20 μM) in which were found in the terminal ileum and colon (Maier et al., 2018).

TABLE 2

Table 2. Antibiotics, drugs, xenobiotics tested on the HGMM and single strains.

Construction of the Human Gut Microbiota Model

Each bacterial strain was cultivated separately under strict anaerobic conditions in a final volume of 5 ml using mGAM medium. Inoculations (5%, v/v) were performed using exponentially growing cells. After incubation at 37°C for 48 h, the 39 strains (1 mL each) were pooled together in a 100 mL penicillin flask to construct a representative in vitro HGMM (Supplementary Figure 1A). To ensure the stability of the established HGMM, the consortium (mixture of 39 strains) was sub-cultured three times (defined as consortium subculture C1, C2, and C3) in the same growth and incubation conditions until obtaining a stable inoculum for further antibiotic, drug, and xenobiotic experiments. To ensure the stability of the HGMM after the third subculture C3 (used as inoculum of the experiments), an extra subculture C4 was carried out, incubated while testing the molecules, to ensure that there is no deviation of individual bacterial strain when pooled together in the constructed HGMM consortium. Therefore, C4 was considered as positive control (PC) for our experiments to observe the deviation the microbial communities due to the addition exposure to the tested xenobiotics, drugs, and antibiotics.

Sensitivity of the Human Gut Microbiota Model and the Individual Bacterial Strains to Antibiotics, Drugs, and Xenobiotics

For the HGMM, all tests were performed in duplicate using Hungate tubes (final volume of 5 ml) at 37°C and pH 7.3. Inoculation (5%, v/v) was made using exponentially growing cells on mGAM medium from the third subculture (named consortium C3) of the HGMM. The sensitivity of the gut microbiota model (Supplementary Figure 1B) was determined with respect to the antibiotics (tested alone or in consortium), drugs, and xenobiotics added to the Hungate tubes at final concentrations defined in Table 2. Growth was determined by monitoring changes in OD_595nm using a spectrophotometer (Camspec spectrophotometer, M107) compared with appropriate negative (mGAM medium containing the molecule of interest but no bacterial consortium) and positive controls (named PC; mGAM medium without molecule but inoculated with the HGMM consortium C3). After 48 h of incubation, the final OD_595nm was recorded, then the bacterial cell pellets were harvested from cultures by centrifugation (15 min, 4°C, 5000 rpm) for further DNA-based approaches analysis. The supernatant from each tested molecule was stored at –80°C for UHPLC-MS/MS analysis.

In addition, the same molecules were tested separately on each single strain used to construct the HGMM. Precultures in final volume of 5 mL of mGAM broth were prepared from pure strains grown on solid medium (using mGAM agar) and incubated for 48 h. After an overnight incubation at 37°C under anaerobic conditions, the bacterial suspensions were adjusted to 1.10⁷ CFU/mL and placed in 96-well plates. The different molecules were added at the concentration previously mentioned (Table 2). After 48 h of incubation, the OD_595nm were measured and then the cultures were centrifuged (30 min, 4°C, 2500 rpm). The supernatants were harvested and transferred into a new 96-well plate and stored at –20°C for analysis by UHPLC-MS/MS. The entire procedure was performed in triplicate under anaerobic conditions (90% N₂, 5% H₂, 5% CO₂) using an anaerobic workstation (Bactron 300, Blanc Labo SA, Switzerland).

Extraction of Genomic DNA and High Throughput Sequencing

DNA was extracted from pelleted cells using the QIAamp PowerFecal Pro DNA Kit with the QIAcube workstation according to Qiagen’s instructions. Samples were previously homogenized twice for 20 s at speed 5.5 m/s with a bead homogenizer (FastPrep instrument, MP Biomedicals). The integrity of the DNA samples was observed with 0.8% agarose gel electrophoresis. The amount of extracted DNA was quantified using a NanoDrop™ One spectrophotometer (Thermo Fisher Scientific™). DNA samples were diluted to a standard concentration of 5 ng/μL and stored at –20°C until use for molecular applications. Bacterial community structure was assessed by sequencing the hypervariable V5–V7 region of the 16S rRNA gene amplified using the primers 799F (5′-AACAGGATTAGATACCCTGG-3′) and 1193R (5′-ACGTCYTCCCCACCTTCC-3′) (Beckers et al., 2016). This region was selected based on a preliminary in silico multiple sequence alignment of the 16S ribosomal RNA (rRNA) sequences from the 39 strains constructing the HGMM in order to allow their discrimination at the species level (Supplementary File Data 2). All PCR reactions were prepared in triplicate for each sample in a total volume of 25 μL PCR mix with 0.3 μM of each forward and reverse primer, 1X PCR mix buffer with 200 mM of each dNTP, and 1U of Q5 high fidelity DNA polymerase (NEB), using a vapo.protect Mastercycler (Eppendorf). PCR cycling was performed at 98°C for 30 s, followed by 25 cycles at 98°C for 10 s, 56°C for 15 s, and 72°C for 10 s, and a final elongation at 72°C for 2 min. The amplified products were checked by 2% agarose gel electrophoresis. Triplicate PCR products for each sample were pooled, purified using the PCR clean-up Kit with the QIAcube workstation according to the manufacturer’s instructions (Qiagen), and quantified using a Qubit 2 fluorimeter (Invitrogen). Sequencing was carried out on the Illumina MiSeq sequencing platform with the V3 Illumina kit (2 × 300 bp paired-end reads) at the ICM Institute (Paris, France¹) according to standard protocols.

Raw datasets of sequencing can be found into the NCBI Sequence Read Archive database under the project accession number PRJNA733465.

Processing of Sequencing Data From Illumina MiSeq, Phylogenetic Relationship, and Statistical Analysis

Bioinformatics data treatments were performed using QIIME2 (version qiime2/2021.2) pipeline (Bolyen et al., 2019). All of the sequences were quality filtered, trimmed, merged, and denoised using the DADA2 algorithm (Callahan et al., 2016), then followed by removing of chimeric sequences. The amplicon sequence variants (ASVs) were grouped at an identity threshold of 97% using q2-vsearch. The ASVs were then taxonomically classified with 90% threshold based on the HGMM reference database (designed in this study) build using the full 16S rRNA gene sequences of the 39 strains constructing the HGMM withing their respective majority taxonomy seven levels (Supplementary File Data 3). From the taxonomic affiliation using the HGMM database, a representative sequence of each ASV affiliated was extracted. Finally, a normalization of the abundances obtained for each ASV was carried out by the DESeq method (Anders and Huber, 2012). Diversity and evenness indices were calculated using the bacterial ASVs data after normalization. The evolutionary history was inferred using the maximum likelihood method (Saitou and Nei, 1987). An alignment was then carried out between ASVs, and the alignment was used to build a phylogenetic tree with PhyML algorithm (Guindon et al., 2010).

The statistical calculations and multivariate analyses were conducted using R software packages (² R Core Team, 2020) with the following libraries: phyloseq, vegan, ape, heatmap2, gclus, ggplot2, GUniFrac, grid, and optparse. Permutational multivariate analysis of variance (PERMANOVA) was used to compare the mean values with respect to the samples at p < 0.05 using stats package. A principal component analysis (PCA) was performed to evaluate structure and the distribution of the samples based on the bacterial abundances in the samples.

Antibiotic, Drug, and Xenobiotic Quantification by UHPLC-MS/MS

UHPLC-MS/MS analysis was performed by a 1290 Infinity Binary LC system (Agilent Technologies, Waldbronn, Germany) coupled to a triple quadrupole or Q-TOF, Q-TRAP 5500 mass spectrometer with an ESI Turbo V ion source (SCIEX, Foster City, CA, United States) or a 6550 iFunnel accurate mass quadrupole time of flight tandem mass spectrometer with a dual Agilent Jet Stream source (Agilent Technologies, Santa-Clara, CA, United States), respectively.

Chromatographic separation was performed on C₁₈ column, Poroshell 120 EC-C18 (Agilent), or Atlantis Premier BEH C18 AX (Waters) for analysis using Q-TRAP or Q-TOF, respectively. The mobile phase consisted of two solutions including solvent A (0.1% formic acid in water) and solvent B (0.1% formic acid in acetonitrile), the column was thermo stated in an oven at 35°C, and the flow rate was set to 650 μl/min. The chromatographic gradient used for each of the compounds was specific; all details are shown in Supplementary Tables 1, 2.

For mass spectrometry analysis using triple quadrupole, data was acquired using electrospray ionization (ESI) in the positive mode, and the ion Spray Voltage was set at 5 500 V. Data was also acquired in negative mode, and the ion spray voltage was set at –4 500 V. The desolvation in source was accomplished using the following set parameters: Temperature (TEM) at 600°C, Ion Source Gas 1 (GS1) at 40 psi, Ion Source Gas 2 (GS2) at 60 PSI, and Curtain Gas (CUR) at 30 psi. The specific parameters of multiple reaction monitoring which permit quantifying and monitoring the compounds are given Supplementary Table 2. Raw data were processed in Sciex Analyst, and individual AUC (Area Under the Curve) for each analyte in each sample was determined using the MultiQuant software.

For mass spectrometry analysis using Q-TOF, data was acquired using electrospray ionization (ESI) in positive and negative modes. The Capillary Voltage (Vcap), Nozzle Voltage, Fragmentor, and Octopole RF peak were set at 4 000 V in positive and in negative mode, 500 V, 350 V, and 750 V, respectively. The desolvation in source was accomplished by using the following parameters: gas temperature at 230°C, drying gas flow at 15 L/min, nebulizer at 40 psi, sheath gas temperature at 280°C, and sheath gas flow at 11 L/min. The profile and centroid data were collected from 100 to 1700 m/z with an acquisition rate of three spectra per second. Raw data were processed in Agilent MassHunter, and individual AUC (Area Under the Curve) for each analyte in each sample was determined using the MassHunter Quantitative software.

Results

Construction of the Human Gut Microbiota Model

The 39 strains from the representative panel of human gut microbial species were pooled together to construct the artificial HGMM used in this study (Supplementary Figure 1). To ensure the reproducibility of the experiments based on a stable model in terms of bacterial density and diversity, the initial constructed consortium was subcultured three times before its use as inoculum. The evolution of the HGMM was analyzed by high-throughput sequencing during the three subcultures named C1, C2, and C3 (Figure 1). In addition, to ensure that there is no deviation of individual bacterial strain when pooled together in the constructed HGMM consortium, an extra subculture C4 was performed and used as positive control (PC) for the experiments. The duplicates for each subculture were regrouped together, underscoring the reproducibility of the analyses (Figure 1A). These findings indicate the stability of the constructed HGMM after the third subculture. The C3 subculture and the PC resulting to a next 48 h incubation of the HGMM C3 subculture without exposure to molecules (thus equivalent to a C4 subculture) were closely related and clustered together but were significantly discriminated from C1 and C2 along the first axis which explained 81.3% of the total variability (Figure 1A). The comparison of the relative abundance of populations at phylum level of each consortium showed a decrease in the abundance of Firmicutes (from 72.3% to 40.0%; mean value of two replicates) in favor of Fusobacteria (increasing from 0.1% to 18.6%) and Bacteroidetes (increase from 6.4% to 22.5%) between C1 and C3. Importantly, no significant difference between the percentages of the most abundant five phyla was observed between C3 and PC (Figure 1B), in agreement with the p-Value from the PERMANOVA test (p-Value: 0.001) performed for analysis of variance using Bray-Curtis distance matrices. This trend was confirmed at the genus level (Figure 1C, see details below). These results argue for the stability of the constructed HGMM after the C3 subculture used as inoculum for the experiments with antibiotics, drugs, and xenobiotics.

FIGURE 1

Figure 1. Characterization and evolution of the bacterial community composition of the HGMM. (A) Principal component analysis based on weighted UniFrac β-diversity metric, showing the repartition of the three HGMM subcultures C1, C2, C3 compared to the positive control (PC) resulting to a next 48 h incubation of the HGMM C3 subculture without exposure to molecules (thus equivalent to a C4 subculture). Total variance explained by the two axes was 99.6%. (B) Barplots showing the relative abundances of the HGMM bacterial community through the subcultures at panel (C) phylum level and at panel (C) genus level. Results from two independent replicates analyses (named.1 and.2) are reported for each subculture and for PC.

Analysis of the Bacterial Community Structure in the Human Gut Microbiota Model

MiSeq sequencing generated a total of 3,729,686 sequences from the 40 samples implemented for HGMM-antibiotic, drug, and xenobiotic interactions. After bioinformatic filters and homogenization steps, 3,290,572 sequences with a mean length of 277 bp were retained, corresponding to 88.2% of the effective sequences. Eighty ASVs were defined at an identity threshold of 97%. The lower values found for Shannon, Simpson, and reciprocal Simpson indexes than for Chao1 indicated an uneven distribution of the bacterial population in the samples (Table 3). Interestingly, all the diversity indexes were lower for the sample exposed to metronidazole compared to the other molecules.

TABLE 3

Table 3. Number of sequences of amplicon sequence variants, and diversity index for experiments with molecules and the positive control.

The comparison of each duplicate showed that there is no significant difference between the percentages of the relative abundances at genus level inside the same sample. These results were confirmed by the Principal Component Analysis (PCA; Figure 2) based on the weighted UniFrac β-diversity metric and showing that the duplicates for each sample were closely clustered. As a result, from now on, the duplicate values will be grouped and presented together (mean values of two replicates). Considering all the molecules, the samples exposed to antibiotics could be separated from all the others along the first axis explaining 38.3% of the variability (Figure 2A). Another PCA, excluding the antibiotics, revealed a distribution of the remaining samples along the first axis, explaining 48.6% of the variability (Figure 2B).

FIGURE 2

Figure 2. Repartition of the HGMM bacterial community after exposure to antibiotics, drugs and xenobiotics. (A) Principal component analysis based on weighted UniFrac β-diversity metric considering all the antibiotics, drugs and xenobiotics tested in the study. Total variance explained by the two axes was 64.4%. (B) The same analysis excluding the eight antibiotics. Total variance explained by the two axes was 60.1%. Results from two independent replicates analyses are reported for each molecule (written “1” and “2” after the name of the molecule). Antibiotic, drug and xenobiotic abbreviations were mentioned in Table 2.

At species level, 28 strains (out of the 39 strains used to construct the HGMM; Table 1) were found present and could be affiliated at the species level in the HGMM by MiSeq sequencing (Supplementary Figure 2), and therefore conserving a 71.8% of the diversity. The non-detected species in the final HGMM belong mainly (7 out of 11) to the Firmicutes, which is consistent with the decrease in abundance of this phylum. Interestingly, the only strain of the Fusobacteria phylum, F. nucleatum, represented only 0.36% of the final diversity and occupies 18.6% of the bacterial population of HGMM. Unfortunately, the only representative of the Verrucomicrobia phylum, A. muciniphila, was not detected in the final HGMM, and therefore its contribution in the response of the model consortium to antibiotics, drugs, and xenobiotics could not be evaluated based only on the MiSeq sequencing.

Sensitivity of the Human Gut Microbiota Model to Antibiotics, Drugs, and Xenobiotics

The 80 ASVs were taxonomically affiliated (Supplementary Figure 2) leading to the identification of the five major phyla (Figure 3A). Thus, at the phylum level, Firmicutes (ranged from 22.3% to 89.3%) was the most abundant phyla in all samples, followed by Bacteroidetes (from 5.2% to 30.4%), Proteobacteria (from 2.4% to 28.8%), Fusobacteria (from 0.1% to 12.4%), and finally Actinobacteria (from 0.1% to 11.8%). The relative abundance of the sequences unassigned at the phylum level was below 0.1%.

FIGURE 3

Figure 3. Relative abundances of the HGMM bacterial community in presence of antibiotics, drugs and xenobiotics. (A) Barplots showing the relative abundances of the HGMM bacterial community at the phylum level and (B) at the genus level following the 48 h of incubation in the presence of the molecules and compared to the HGMM control cultures without molecule (PC).

At the genus level, the relative abundance of the microbial communities changed mainly in response to antibiotic exposure in comparison to the respective positive control (HGMM subcultured on mGAM without any molecule) (Figure 3B). These results were in accordance with the PCA represented on Figure 2. The Proteobacteria strain (Escherichia genus) was sensitive to the cefpodoxime, erythromycin, trimethoprim, and STX. However, the Fusobacteria (Fusobacterium genus) was inhibited by the amoxicillin, cefpodoxime, erythromycin, sulfamethoxazole, and STX. The Bacteroidetes genera were affected in presence of metronidazole, trimethoprim, sulfamethoxazole, and STX. The Actinobacteria (Eggerthella and Bifidobacterium genus) were sensitive to both erythromycin and moxifloxacin, whereas the Firmicutes (Enterococcus) were inhibited mainly in presence of the moxifloxacin. When excluding the antibiotics from the analyses, a profile similar to the PC was observed, indicating no significant changes in the relative abundances of the samples exposed to drugs or to xenobiotics that included allergy, analgesics, antidepressant, cardiology-angiology, hepatology-gastroenterology, metabolism-nutrient, nonsteroidal anti-inflammatory, oncology-hematology, pesticides, plastics industry, and preservatives (Figure 3B).

As an alternative approach to identify the active strains in each sample, the ASVs sequences found out in all the samples by high throughput sequencing were identified and classified after taxonomic affiliation at species level. The sensitivity of the 28 strains found out in the samples was evaluated and represented in a heatmap (Figure 4) based on the relative abundances of each strain into all samples exposed to antibiotics, drugs, xenobiotics, and in the positive control. Based on their relative abundances, the detected strains could be grouped into four clusters in response to the molecule exposure. Interestingly, results from this clusterization were not related to taxonomic relationship or molecule category (Tables 1, 2). The Cluster III was the less abundant in terms of bacterial species, followed by Cluster I and Cluster IV. The Cluster II was the most abundant in the HGMM. The maximum of relative abundance was observed in the Cluster I with C. difficile (15.1%) for the HGMM exposed to erythromycin. Species from Cluster II were the most abundant and detected in all samples, with the exception for erythromycin and metronidazole exposures that completely inhibited B. vulgatus and F. nucleatum, respectively. Within Cluster II, E. faecalis was the most abundant species (relative abundance of 14.7%). When compared to the positive control, strains from the Cluster III were mainly inhibited by erythromycin, metronidazole, trimethoprim, STX, and moxifloxacin. The maximum relative abundance in the Cluster III was observed for B. fragilis (10.4%) in the presence of erythromycin. Interestingly, B. thetaiotaomicron formed the separated Cluster IV. The maximum of B. thetaiotaomicron relative abundance was detected following erythromycin exposure (12.4%). Overall, these results were in accordance with the findings from Figures 2, 3 showing that the human gut model composition was affected mainly by antibiotics compared to the other drugs and xenobiotics. The unassigned sequences might be due to the lack of the reference sequences in the database or to the quality of sequences in the respective samples C1, C2, fludioxonil, pyrimethanil, and erythromycin (relative abundance of 0.02%, 0.04%, 0.05%, 0.1%, and 0.1% of unassigned sequences, respectively).

FIGURE 4

Figure 4. Human gut microbiota model sensitivity at species level to antibiotics, drugs and xenobiotics. Heatmap showing the sensitivity of the bacteria from the consortium model in presence of molecules, following the 48 h incubation compared to the HGMM control cultures without molecule. The color code indicates the relative abundance of the HGMM strains ranging from white (low abundance) to blue (high abundance) in response to the exposure to the different molecules. The clusterization of strains (on the left) was based on their sensitivity to the tested molecules (not on a taxonomic relationship). The clusterization of the molecules (on the top) was based on their effect on the HGMM strains.

Effects of Human Gut Microbiota Model on Antibiotics, Drugs, and Xenobiotics

Antibiotics, drugs, and xenobiotics were quantified from the HGMM supernatant using an UHPLC-MS/MS approach after 48 h of exposure. These results were combined with the global bacterial growth of the HGMM in the presence of the tested molecules and monitored by measuring OD_595nm. In Figure 5, the data have been transformed into percentage of growth and percentage of absence of xenobiotics as compared to the control without bacteria (corresponding to 0% absence). Five clusters were observed. Cluster I consists of the molecules that did not impact the growth of HGMM, but which were transformed and were not found in the culture supernatants (approximately 100% absence). This was the case for four antibiotics including metronidazole, cefpodoxime, and sulfamethoxazole alone or when coupled with trimethoprim, as well as for one drug (hesperidin). A second cluster appears when the growth of HGMM was not impacted and only 25 to 50% of the amount of the molecules were detected. Cluster II includes four drugs (olsalazine, nisoldipine, nifedipine, and diclofenac sodium), one antibiotic (erythromycin), and one xenobiotic (dioctyl phthalate). Cluster III is comprised of the molecules whose quantities were decreased by approximately 25% and consists of one drug (fexofenadine hydrochloride) and five xenobiotics (difenoconazole, pyrimethanil, fludioxonil, butylparaben, and propylparaben). Cluster IV represents the molecules that did not inhibit the growth of HGMM and whose amounts have not been modified. This cluster is composed of three antibiotics (amoxicillin, trimethoprim alone or coupled with sulfamethoxazole), ten drugs (acetaminophen, warfarin, bisoprolol, aceclofenac, irinotecan, mercaptopurine, ascorbic acid, omeprazol, topotecan, and clomipramine hydrochloride), and six xenobiotics (di-isobutyl phthalate, di-n-butyl phthalate solution, boscalid, bisphenol A, and glyphosate). Cluster V includes only one molecule: the moxifloxacin antibiotic. This is the only molecule that inhibited the global bacterial growth of our HGMM and, moreover, without modification of its amount, according to UHPLC-MS/MS analyses.

FIGURE 5

Figure 5. Human gut microbiota model bacterial growth versus absence of the native antibiotics, drugs and xenobiotics. Biplot showing the percentage of bacterial growth of the HGMM versus the percentage of absence of the native molecules in the HGMM culture supernatant. The percentage of bacterial growth of the HGMM was calculated in the presence of antibiotics, drugs and xenobiotics compared to HGMM cultures without molecule. The percentage of absence of native antibiotics, drugs and xenobiotics in HGMM culture supernatant were calculated from the respective control cultures without HGMM bacterial community. STX (sulfamethoxazole) and STX (trimethoprim) indicate that only sulfamethoxazole or trimethoprim was analyzed, respectively.

Effects of Antibiotics, Drugs, and Xenobiotics on the Individual Growth of Each Bacterial Strain Constructing the Human Gut Microbiota Model

Each bacterial strain selected to construct the HGMM (Table 2) was exposed to the different molecules to access their individual susceptibility or resistance by an estimation of the bacterial growth parameter. The percentage of growth was calculated for each strain in the presence of the tested molecule with respect to the control without this molecule (Figure 6). Results indicated that effects of antibiotics, drugs, and xenobiotics on the growth of each bacterial strain constructing the HGMM were not related to taxonomic relationship or molecule category (Tables 1, 2).

FIGURE 6

Figure 6. Sensitivity of each strain constructing the HGMM to antibiotics, drugs and xenobiotics in individual cultures. The growth percentage of each strain (exception with C. leptum, R. intestinalis, and R. bromii) in the presence of a molecule was calculated from the control cultures without this molecule. Boscalid was not tested on individual strains. The color code indicates the growth percentage for each strain ranging from blue (no growth) to red (increased growth). The clusterization of strains (on the left) was based on their sensitivity to the tested molecules (not on a taxonomic relationship). The clusterization of the molecules (on the top) was based on their effect on the growth of the strain.

The molecules limiting mostly the growth of almost all the strains were antibiotics such as amoxicillin, metronidazole, STX, moxifloxacin, erythromycin, and cefpodoxime. Overall, drugs and xenobiotics seemed to not impact the bacterial growth, with a few notable exceptions. Remarkably, the growth of E. lenta was higher than controls (growth rate superior to 150%) in the presence at least of 20 molecules belonging to different categories (Table 2) from pesticides, plastic industry, and preservatives to analgesics, cardio-angiology or non-steroidal anti-inflammatory, as example. Two other strains belonging also to the Actinobacteria phylum, C. aerofaciens and B. adolescentis, increased their growth compared to their respective controls in presence of eight molecules belonging mostly to plastic industry and preservatives categories. The growth of B. wadsworthia (Proteobacteria phylum) was stimulated in the presence of the two molecules tested from the nonsteroidal anti-inflammatory category and of three out of four molecules tested from the cardio-angiology category, as an example. The growth of C. saccharolyticum (Firmicutes phylum) was also impacted by the presence of some drugs and xenobiotics, as shown in Figure 6. Some species-specific responses to antibiotics, drugs, and xenobiotics were also observed such as the growth of A. muciniphila only stimulated in presence of butylparaben, as example.

Effects of Individual Bacterial Strain Constructing the Human Gut Microbiota Model on Antibiotics, Drugs, and Xenobiotics

The absence percentage of each molecule in the bacterial supernatants compared to the control (medium containing the molecule but not bacteria) for each bacterial strain constructing the HGMM was analyzed (Figure 7). Results indicated that effects of each bacterial strain constructing the HGMM on the biotransformation of molecules were not closely related to taxonomic relationship or molecule category (Tables 1, 2). The cluster of the molecules disappearing mostly from the culture supernatants (approximately 100% absence) included sulfamethoxazole, metronidazole, omeprazol, nisoldipine, and also olsalazine. Interestingly, these molecules belong to clusters I and II previously discriminated in Figure 5 based on the effects of HGMM on molecules. A large panel of strains belonging to all the phyla (Table 1) seemed to carry out the total biotransformation of these molecules. Eight molecules belonging to different categories and including cefpodoxime, fexofenadine hydrochloride, butylparaben, difenoconazole, diclofenac sodium, clomipramine chlorhydrate, di-isobutyl phthalate, and di-n-butyl phthalate solution appeared also altered by a panel of strains.

FIGURE 7

Figure 7. Analysis of the native antibiotics, drugs or xenobiotics for each strain constructing the HGMM. Heatmap highlighting the absence of the native molecules in the culture supernatants of each strain separately (exception with C. leptum, R. intestinalis, and R. bromii) using UHPLC-MS/MS analysis. Percentages of the absence of the native antibiotics, drugs or xenobiotics in the culture supernatant for each strain was calculated from the respective control cultures without bacteria. The color code indicates the absence of the molecule in the culture supernatant ranging from blue (total presence) to red (total absence). The clusterization of strains (on the left) was based on the absence of the different native molecules in the culture supernatant (not on a taxonomic relationship). The clusterization of the molecules (on the top) was based on their absence in the supernatant of each strain. Method failed with amoxicillin, ascorbic acid, bisphenol A and irinotecan. Boscalid was not tested on individual strains. STX (SMT) and STX (TMP) indicate that only sulfamethoxazole (SMT) or trimethoprim (TMP) was analyzed, respectively.

In contrast, a cluster displayed 100% molecule detection in the culture supernatants after 48 h exposure, highlighting the absence of potential biotransformation by gut species. This included acetaminophen, bisoprolol, glyphosate, aceclofenac, and topotecan. As well, these molecules belong to cluster V previously identified in Figure 5. Remarkably, the strains mostly implied in the alteration of a large panel of molecules belong all to the Bacteroidetes phylum with P. merdae, B. vulgatus, O. splanchnicus, B. thetaiotaomicron, B. caccae, B. uniformis, and P. distasonis, as observed in Figure 7. Interestingly, these strains were not listed among the five strains with a growth higher than controls identified from Figure 6. Also, these strains did not display the highest relative abundance in HGMM (Figure 4), with the exception of B. vulgatus. Some species-specific responses to antibiotics, drugs, and xenobiotics were also observed (Figure 7) such as the total biotransformation of methylparaben by only A. muciniphila and B. obeum or of moxifloxacin by only C. perfringens and P. parasanguinis, as example. Remarkably, some species which are not detected in the final HGMM (Table 2), but able to grow in mGAM medium (Figure 6), appeared also to carry out significant biotransformation (>60%) of several molecules, such as C. aerofaciens, O. splanchnicus, P. copri, or B. obeum, as examples.

Reciprocal Interactions With Antibiotics, Drugs, and Xenobiotics of Human Gut Microbiota Model and Individual Species

The reciprocal interactions between antibiotics, drugs, and xenobiotics and our HGMM are summarized in Table 4. The molecules did not impact, except moxifloxacin, the global bacterial growth of our HGMM. In addition, the molecules did not impact, intended exception to antibiotics, the community structure of our HGMM. This also demonstrated the stability of our constructed model. However, our HGMM composed of 28 strains harbored an enormous metabolic potential, which can partially or totally transform molecules, whatever its category. This may strongly alter the metabolism responses of orally administered drugs or favor elimination of adverse xenobiotics.

TABLE 4

Table 4. Summary of the HGMM reciprocal interactions with the 36 tested antibiotics, drugs, and xenobiotics, and prediction of bacterial species potentially implied in its transformation.

Discussion

In this study we analyzed the effects of various molecules including antibiotics, drugs, and xenobiotics on both single-cultured strains selected from the human gut microbiota core and from an in vitro human gut microbiota model (HGMM) that we developed. Reciprocally, we also determined the effect of these bacterial strains, alone or in consortium (i.e., HGMM), on the quantity of the native molecule in the supernatant after 48 h of exposure. The investigation of the reciprocal effect that the tested molecules may have on the HGMM, or conversely the effect of the HGMM on the molecules, if incubated for a long period (more than 48 h), was not evaluated in this study. Using 48 h as time of incubation was chosen based on the time needed for the absorption of drugs, which usually takes up to 6 h after an oral administration (Tobío et al., 2000; Umigai et al., 2011). Therefore, 48 h of incubation can serve as high throughput initial investigations with early predictive responses, providing important clues to guide further studies of new candidate pharmaceutical, food, or environmental molecules. The concentrations of the tested molecules belong to the range of the screening concentration (below 20 μM) in which were found in the terminal ileum and colon, where most gut microbes reside (Donaldson et al., 2016). However, drug concentrations are only quantified in blood, and human-targeted drugs have order of magnitude amounts lower than in our screen (Maier et al., 2018).

The HGMM was constructed on mGAM medium using 39 species (Maier et al., 2018; Johnson et al., 2019; Zimmermann et al., 2019) belonging to the most common dominant phyla in human gut microbiota (Lozupone et al., 2012; Almeida et al., 2019). mGAM medium was chosen for the culture of individual strains and for HGMM construction since all the selected species were reported to grow robustly in this medium in a manner that is reflective of their gut abundance (Maier et al., 2018; Tramontano et al., 2018). In order to be representative of the gut microbiota of healthy individuals, selection of the strains was performed based on the approach followed by Maier et al. (2018) to select a set of ubiquitous gut bacterial detected at a relative abundance of ≥ 1% and prevalence of ≥ 50% in fecal samples of asymptomatic humans from three continents. All the 39 species are found in the gut of healthy individuals and are part of a larger strain resource panel for the healthy human gut microbiome (Tramontano et al., 2018).

The gut microbiota consists of hundreds of bacterial species showing stability over time (Faith et al., 2013; Hisada et al., 2015). This stability can be disturbed by some perturbations such as dietary or drugs administration that may affect the function and the diversity of the gut microbiota (Relman, 2012; Hemarajata and Versalovic, 2013). Using a stable microbiota is mandatory to analyze and predict the bacterial community behavior in an in vitro human gut model. In this study, the 16S rRNA gene high throughput sequencing results demonstrated the stability of our constructed model through subculturing, and no significant difference was observed after the third subculturing (Figure 1A). During the subculturing, a decrease mainly in the abundance of the Firmicutes and Proteobacteria (Figure 1B) in favor of Bacteroidetes and Fusobacteria phyla was observed. These changes in abundances could be explained by syntrophic and competition interactions between the 39 strains. Results from literature indicate that the gut microbial communities were shaped by syntropy as well as competition interactions involving nutrients, but also other interactions concerning quorum sensing, as an example, causing cell lysis. Thus, some interactions will play a role in shaping the dynamics and metabolic efficiency of community, whereas others will stabilize cooperative networks by introducing negative feedbacks (Coyte et al., 2015; Wasielewski et al., 2016; Venturelli et al., 2018; Coyte and Rakoff-Nahoum, 2019). These interactions reveal positive and negative responses allowing the dominance of some phyla at the expense of others (Hibbing et al., 2010; Foster and Bell, 2012). A recent study reported that these changes in the Firmicutes/Bacteroidetes ratio are widely accepted to have an important influence in maintaining normal intestinal homeostasis (Stojanov et al., 2020) as the increases in the abundance of specific Firmicutes or Bacteroidetes species could lead to obesity and bowel inflammation, respectively (Shen et al., 2018; Abenavoli et al., 2019). Our results were in accordance with the findings obtained in an in vitro fermentation study from fecal microbiota transplantation donors illustrating the decrease in Firmicutes/Bacteroidetes ratio from 1.46 to 0.98 (Fu et al., 2019).

At a global scale, the diversity analysis revealed the presence of 80 ASVs belonging to five phyla (Firmicutes, Bacteroidetes, Proteobacteria, Fusobacteria, and Actinobacteria) which were found to be present in the stable HGMM (Figure 1A and Supplementary Figure 2). These results were consistent with other studies reporting that these phyla were known to dominate the human gut microbiota (Iizumi et al., 2017; Senghor et al., 2018; Rinninella et al., 2019). At the species level, 28 strains out of the 39 used to construct the HGMM were found present (Figure 4 and Supplementary Figure 2) and could be affiliated at the species level. Observing 28 strains out of 39 could be explained by the short length of the ASVs sequences obtained using the universal primers 799F and 1193R which generate only 394 bp. After quality trimming, an average sequences length of 277 bp was retained, which is not enough to discriminate the taxonomic affiliation of the strains at species level. These factors may explain the non-detection of some species such as Akkermansia muciniphila closely associated with the protective mucous lining of the human intestine (Schneeberger et al., 2015). Therefore, using other approaches (e.g., PacBio or shotgun metagenomic sequencing) to obtain a full length 16S rRNA gene sequences of the HGMM could be of interest to give more resolution on the composition of the HGMM and thus may increase the number of the strains found out of the 39 strains used to construct the HGMM. However, this decrease could be also related to the positive and negative responses allowing the dominance of some strains to the detriment of others. However, even with the 28 strains out of 39 used to construct the HGMM, our model still conserves 71.8% of the initial diversity in the HGMM, which is higher than other studies reporting that the operational taxonomic units detected by 16S rRNA gene sequencing of cultivable fecal samples were only in a maximum of 40–50% (Goodman et al., 2011; Rettedal et al., 2014; Lau et al., 2016).

Considering all samples (Table 2), the analysis of distribution of microbial populations in the HGMM based on principal component analysis (Figure 2) and the PERMANOVA test for analysis of variance using Bray–Curtis distance (p-Value: 0.001) showed no significant difference between the percentages of abundances (whatever is the taxonomic affiliation level of the ASVs) of each sample duplicate. The global profile of the five phyla abundances was similar for the majority of the tested molecule categories compared to the positive control. However, based on the Hill numbers (Chao et al., 2014) shown in Table 3 as well as the beta diversity, the HGMM was significantly affected in the presence of the antibiotics compared to the other tested molecule categories (Figures 4–6). These results were in accordance with other studies reporting that antibiotics can affect the abundances of 30% of the bacteria in the gut community, leading to rapid and significant discrepancies in taxonomic richness, diversity, and evenness (Dethlefsen et al., 2008; Dethlefsen and Relman, 2011). Molecular omics works have shown that, in addition to the composition of taxa alteration, antibiotics also affect the gene expression, protein activity, and overall metabolism of the gut microbiota (Franzosa et al., 2015; Francino, 2016).

At species level (Figure 4), C. difficile (Cluster I) was the most abundant of all strains in the presence of erythromycin (relative abundance of 15.1%). Recently, antimicrobial resistance patterns were observed in C. difficile strains against erythromycin, metronidazole, and moxifloxacin (Tilkorn et al., 2020). However, erythromycin and metronidazole completely inhibit B. vulgatus and F. nucleatum, respectively. It was reported that Bacteroides spp. were not affected by metronidazole and were carrying nim genes for metronidazole resistance (Nagy et al., 2001; Teng et al., 2002; Gal and Brazier, 2004). B. wadsworthia and P. distasonis were completely inhibited in the presence of erythromycin and metronidazole. These results were in accordance with previous works showing the susceptibility of both genera Bilophila and Parabacteroides to erythromycin and metronidazole (Schumacher and Single, 1998; Song et al., 2005; Awadel-Kariem et al., 2010; Wang et al., 2015). In Cluster II, E. faecalis was the most abundant strain (relative abundance of 14.7%). Previous studies reported prevalence of acquired resistance of E. faecalis to many antibiotics and drugs (Nicas et al., 1989; Aarestrup et al., 2000; Butaye et al., 2001; Ono et al., 2005). From cluster III, it was observed that B. fragilis was dominant in presence of erythromycin. B. thetaiotaomicron formed a separated cluster IV and was abundant in presence of erythromycin (12.44%). Previously, many studies reported Bacteroides conjugative transposons carry both tetracycline and erythromycin resistance genes, which means the use of tetracycline selects for erythromycin resistant strains and vice versa, adding to the antibiotic resistance problem (Privitera et al., 1979; Whittle et al., 2002; Johnsen et al., 2017; Niestêpski et al., 2019). Overall, these results agreed with the findings showing that the human gut model was affected mainly by antibiotics compared to the other drugs and xenobiotics (Figures 2–5).

The single-culture strains were shown to be particularly sensitive to moxifloxacin and trimethoprim. Moxifloxacin is an antibiotic known to decrease bacterial diversity in the feces of treated individuals (Burdet et al., 2019) which is consistent with our previous study showing that using moxifloxacin in vitro significantly inhibited the growth of 12 strains commonly found in the intestinal microbiota out of the 30 tested (Ecale et al., 2021). Surprisingly, moxifloxacin was the only antibiotic that inhibited the growth of HGMM (Figure 5) suggesting that its action on the bacteria of the intestinal microbiota core is deleterious and that it must be used wisely to avoid significant dysbiosis (Singh et al., 2017).

We have also shown that trimethoprim affected the growth of single strains but also had an impact on the growth of HGMM (Figures 2, 3, 5). STX also affects the composition of the HGMM, with an increase in Actinobacteria and Firmicutes but a decrease in Proteobacteria, Fusobacteria, and Bacteroidetes. However, a smaller number of strains were affected, indicating that the trimethoprim/sulfamethoxazole combination had less impact on the composition of the microbiota than trimethoprim alone. The results obtained with our model are in accordance with previous study showing that the treatment of a patient with STX for 2 years did not affect the quantity of culturable anaerobic bacteria (Kofteridis et al., 2004).

In the presence of erythromycin, we showed that the phylum Fusobacteria disappeared. In our model, this phylum is composed of only one F. nucleatum strain. Surprisingly, our F. nucleatum strain was not inhibited in the presence of erythromycin when grown alone. This result may indicate that bacteria behave differently when cultured alone or with other strains, hence the importance of studying strains in single culture but also in consortium. Overall, antibiotics had increased the ratio of Firmicutes in our HGMM except for moxifloxacin (Figure 3A). These results are consistent with a study conducted in a mouse model which showed that treatment with antibiotics (penicillin, vancomycin, and tetracycline) increased the abundance of Firmicutes (Cho et al., 2012). Many studies mention a balance between Firmicutes and Bacteroidetes deregulated by the presence of antibiotics, which we confirmed with our model (Panda et al., 2014; Dudek-Wicher et al., 2018).

Our results show that metronidazole and sulfamethoxazole (alone or with trimethoprim) were absent from the culture supernatants of HGMM (Figure 5 and Table 4) and that these two antibiotics were weakly found in the culture supernatants of the single cultured strains. Recent studies reported that xenobiotics can be impacted by bacterial enzymes of the intestinal microbiota (Maier et al., 2018; Collins and Patterson, 2020). Our results indicated that C. perfringens was able to completely transform metronidazole. Similarly, metronidazole was totally absent from the HGMM supernatant. These results were in accordance with previous studies reporting that metronidazole is well known to be metabolized by strains of the intestinal microbiota and in particular by C. perfringens (Koch et al., 1979; Sousa et al., 2008; Koppel et al., 2017).

The absence of some antibiotics in HGMM culture supernatants, such as cefpodoxime (100% absent) and erythromycin (75% absent), suggest that strains of the intestinal microbiota have a close interaction with xenobiotics (Clarke et al., 2014; Koppel et al., 2017; Collins and Patterson, 2020). Overall, the results of single strain culture supernatant analyses showed that in addition to metronidazole and sulfamethoxazole, olsalazine, omeprazol, diclofenac, fexofenadine, nisoldipine, nifedipine, hesperidin, di-isobutyl phthalate, and di-n-butyl phthalate were very rarely found in bacterial supernatants (Figures 6, 7). Some of these molecules are known to be metabolized by the intestinal microbiota but not all bacteria are involved (Sousa et al., 2008; Klaassen and Cui, 2015; Lu et al., 2015; Enright et al., 2016; Jourova et al., 2016; Koppel et al., 2017; Currò, 2018; Clarke et al., 2019; Sun et al., 2019). For example, we showed here that Bacteroidetes as well as Firmicutes strains were able to completely transform the olsalazine, and we then put forward the hypothesis that these phyla could metabolize this molecule.

Omeprazole and clomipramine chlorhydrate were absent from the supernatant of several individual strains of the intestinal microbiota core in our study. Surprisingly, these molecules did not disappear in HGMM (0% absence) which is consistent with the hypothesis of Jourova et al., which stated that omeprazol is unlikely metabolized by bacteria of the intestinal microbiota (Jourova et al., 2016). Clomipramine was totally absent in O. splanchnicus. It was reported that desmethylclomipramine was metabolized to an active metabolite without indicating the precise mechanism (Benfield et al., 1980). We can then ask whether the bacterium O. splanchnicus is part of this mechanism.

We observed that di-isobutyl phthalate and di-n-butyl phthalate were absent supernatants of single strains whereas their quantities were not modified by HGMM. Similar results were observed for diclofenac sodium but was 65% absent in HGMM. Some studies have suggested that diclofenac sodium is metabolized by bacterial beta-glucuronidase enzymes from the gut microbiota causing acute diarrhea in patients treated with this drug (Boelsterli et al., 2013; Saitta et al., 2014; Oda et al., 2015). However, this enzyme is not present in all bacteria we tested, indicating that another mechanism may be involved in the disappearance we were observing.

Fexofenadine was absent from half of the culture supernatants of the single strains except for five bacteria (B. ovatus, B. fragilis, E. coli, B. obeum, and C. perfringens) and was absent from 30% of the HGMM supernatants. These results were consistent with a previous study reporting that fexofenadine interact with the intestinal microbiota and consequently decrease its initial concentration (Zou et al., 2020). Nisoldipine and nifedipine used against hypertension were absent from the culture supernatants of C. perfringens and D. formicigenerans and also by about 75% in the culture supernatant of HGMM. Some authors have hypothesized that nifedipine could be metabolized by bacteria of the intestinal microbiota thanks to dehydrogenase and/or nitroreductase enzymes (Choi et al., 2018). We in turn hypothesize that these bacteria could be, among others, D. formicigenerans and C. perfringens, possessing the two enzymes mentioned (Rafii and Cerniglia, 1993; Matsunaga et al., 2018).

Hesperidin, a flavonoid, was absent from the culture supernatants of five strains in single culture and absent from the HGMM one. It is known that hesperidin could be metabolized by the intestinal microbiota (Boonpawa et al., 2017). Thanks to our model, we can also hypothesize that hesperidin is metabolized by bacteria of the intestinal microbiota, particularly by P. merdae, R. gnavus, and S. parasanguinis. Hence, our data highlight the need to consider several approaches (combining molecular and culture-based strategies), single culture, and culture of a set of strains for studying of the interaction between the intestinal microbiota and xenobiotics.

Our HGMM model is considered as an in vitro batch fermentation model. These models are simplest and the most frequently used model to test the ability of specific strains or fecal microbial communities to metabolize different substrates (Venema and Van den Abbeele, 2013). For instance, recent studies of in vitro batch culture models (Day-Walsh et al., 2021; Pérez-Burillo et al., 2021) were easy to set up, useful for fermentation studies, and especially substrate digestion assessment. However, the batch models were mainly limited by the short-term fermentation studies (Pompei et al., 2008; Gumienna et al., 2011). Conversely, the continuous cultures (e.g., PolyFermS, TIM-2) or multistage continuous cultures (e.g., Reading model, SHIME) were characterized by their ability of mimicking the conditions found in vivo and performed in well controlled environmental parameters. However, their weaknesses are related to the absence of host functionality, and the experiments were time limited (Duncan et al., 2009; Maccaferri et al., 2010; Van den Abbeele et al., 2010; Payne et al., 2012). Our HGMM, as well as other in vitro gut fermentation models, remain irreplaceable tools for screening as well as studying the mechanisms of action of prebiotics and probiotics (Pham and Mohajeri, 2018). The HGMM stability was investigated for 48 h and showed that our model was stable after the third subculture. Other models such us the immobilized continuous culture in vitro models were suggested for long-term investigations and reported for their long-term stability of continuous fermentation system with immobilized fecal microbiota functionality (Le Blay et al., 2009; Zihler et al., 2010). However, their weakness is mainly related to the absence of host functionality (Payne et al., 2012). To mimic the conditions found in vivo, advanced artificial digestive systems were suggested. However, these models were limited by the absence of the immune and neuroendocrine response and experiments are limited to few days’ time (Blanquet-Diot et al., 2009; Kovatcheva-Datchary et al., 2009). In summary, the goals of the study determine which model system to choose. Technical aspects, advantages, and limitations of each model should be taken into account for the selection of the suitable one, as they are the decisive determinants of the capabilities of each model. In this study, our HGMM can serve as a high throughput model for screening and short-term investigations with early predictive responses, providing important clues to guide further studies of new candidate pharmaceutical, food, or environmental molecules.

Conclusion

Numerous studies have shown that the human gut microbiota can interact with and metabolize antibiotics, drugs, and xenobiotics. These studies were mainly based on metagenomic analyses of feces, complex experimental models, and animal models, despite technical difficulties and high costs. It would be critical to design a simplified model of the intestinal microbiota for the study of the reciprocal interactions with molecules including antibiotics, drugs, and xenobiotics. In this study, we developed a HGMM, an in vitro model from the most abundant 39 strains of the core human gut microbiota. Our findings aimed at predicting the intended or unintended effects of antibiotics, drugs, and xenobiotics on the intestinal microbiota and vice versa. After its implementation, we showed that the bacterial relative abundances of HGMM were affected by the presence of antibiotics alongside the individual cultured strains composing it. We also showed that some native molecules were absent from culture supernatants of HGMM and individual strains, thus possibly metabolized. These findings qualified the constructed HGMM as a simple, stable, inexpensive, and reproducible model that could be used to predict the effects that molecules may have on the gut microbiota in preliminary studies to develop new pharmaceutical, food, or environmental molecules.

Data Availability Statement

Raw datasets of sequencing can be found into the NCBI Sequence Read Archive database under the project accession number PRJNA733465.

Author Contributions

AE, FE, AM, M-HR, J-MB, and MR conceived and designed the study. AE, FE, AM, SC, and JL conducted the experiments. AE and FE contributed to analyzing the data and wrote and revised the manuscript. All authors discussed the results and commented on the manuscript and approved the final version of the manuscript.

Funding

This present research was supported by the Eurofins Discovery, an internal funding of the CNRS UMR 7267 laboratory and, partly granted by two 2015–2020 programs: State-Region Planning 800 Contracts (CPER) and European Regional Development Fund (FEDER).

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank the CHU of Poitiers for kindly providing the strain of B. thetaiotaomicron.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.828359/full#supplementary-material

Abbreviations

HGMM, Human Gut Microbiota Model; UHPLC-MS/MS, Ultrahigh performance liquid chromatography - tandem mass spectrometer; HuMiX, human microbial cross-talk model; HMI, host-microbiota interaction module; ATCC, American Type Culture Collection; DSMZ, Deutsche Sammlung von Mikroorganismen und Zellkulturen; mGAM Broth, modified gifu anaerobic medium broth; PCA, principal component analysis; DMSO, dimethyl sulfoxide; PC, positive control; OD, optical density; CFU, colony forming unit; ASV, amplicon sequence variant.

Footnotes

1. ^ https://icm-institute.org/fr/
2. ^ https://www.R-project.org/

References

Aarestrup, F. M., Agerso, Y., Gerner–Smidt, P., Madsen, M., and Jensen, L. B. (2000). Comparison of antimicrobial resistance phenotypes and resistance genes in Enterococcus faecalis and Enterococcus faecium from humans in the community, broilers, and pigs in Denmark. Diagn. Microbiol. Infect. Dis. 37, 127–137. doi: 10.1016/S0732-8893(00)00130-9

CrossRef Full Text | Google Scholar

Abboud, C., Duveau, A., Bouali-Benazzouz, R., Massé, K., Mattar, J., Brochoire, L., et al. (2021). Animal models of pain: diversity and benefits. J. Neurosci. Methods 348:108997. doi: 10.1016/j.jneumeth.2020.108997

PubMed Abstract | CrossRef Full Text | Google Scholar

Abenavoli, L., Scarpellini, E., Colica, C., Boccuto, L., Salehi, B., Sharifi-Rad, J., et al. (2019). Gut microbiota and obesity: a role for probiotics. Nutrients 11:2690. doi: 10.3390/nu11112690

PubMed Abstract | CrossRef Full Text | Google Scholar

Almeida, A., Mitchell, A. L., Boland, M., Forster, S. C., Gloor, G. B., Tarkowska, A., et al. (2019). A new genomic blueprint of the human gut microbiota. Nature 568, 499–504. doi: 10.1038/s41586-019-0965-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Anders, S., and Huber, W. (2012). Differential Expression of RNA-Seq Data at the Gene Level–The DESeq Package. Heidelberg: European Molecular Biology Laboratory (EMBL).

Google Scholar

Audebert, C., Even, G., Cian, A., Loywick, A., Merlin, S., Viscogliosi, E., et al. (2016). Colonization with the enteric protozoa Blastocystis is associated with increased diversity of human gut bacterial microbiota. Sci. Rep. 6, 1–11. doi: 10.1038/srep25255

PubMed Abstract | CrossRef Full Text | Google Scholar

Awadel-Kariem, F. M., Patel, P., Kapoor, J., Brazier, J. S., and Goldstein, E. J. C. (2010). First report of Parabacteroides goldsteinii bacteraemia in a patient with complicated intra-abdominal infection. Anaerobe 16, 223–225. doi: 10.1016/j.anaerobe.2010.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Beckers, B., De Beeck, M., Thijs, S., Truyens, S., Weyens, N., Boerjan, W., et al. (2016). Performance of 16s rDNA primer pairs in the study of rhizosphere and endosphere bacterial microbiomes in metabarcoding studies. Front. Microbiol. 7:650. doi: 10.3389/fmicb.2016.00650

PubMed Abstract | CrossRef Full Text | Google Scholar

Benfield, D. P., Harries, C. M., and Luscombe, D. K. (1980). Some pharmacological aspects of desmethylclomipramine. Postgrad. Med. J. 56(Suppl. 1), 13–18.

PubMed Abstract | Google Scholar

Blanquet-Diot, S., Soufi, M., Rambeau, M., Rock, E., and Alric, M. (2009). Digestive stability of xanthophylls exceeds that of carotenes as studied in a dynamic in vitro gastrointestinal system. J. Nutr. 139, 876–883. doi: 10.3945/jn.108.103655

PubMed Abstract | CrossRef Full Text | Google Scholar

Boelsterli, U. A., Redinbo, M. R., and Saitta, K. S. (2013). Multiple NSAID-induced hits injure the small intestine: underlying mechanisms and novel strategies. Toxicol. Sci. 131, 654–667. doi: 10.1093/toxsci/kfs310

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Boonpawa, R., Spenkelink, A., Punt, A., and Rietjens, I. M. C. M. (2017). Physiologically based kinetic modeling of hesperidin metabolism and its use to predict in vivo effective doses in humans. Mol. Nutr. Food Res. 61:1600894. doi: 10.1002/mnfr.201600894

PubMed Abstract | CrossRef Full Text | Google Scholar

Burdet, C., Nguyen, T. T., Duval, X., Ferreira, S., Andremont, A., Guedj, J., et al. (2019). Impact of antibiotic gut exposure on the temporal changes in microbiome diversity. Antimicrob. Agents Chemother. 63:e00820-19. doi: 10.1128/AAC.00820-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Butaye, P., Devriese, L. A., and Haesebrouck, F. (2001). Differences in antibiotic resistance patterns of Enterococcus faecalis and Enterococcus faecium strains isolated from farm and pet animals. Antimicrob. Agents Chemother. 45, 1374–1378. doi: 10.1128/AAC.45.5.1374-1378.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao, A., Gotelli, N. J., Hsieh, T. C., Sander, E. L., Ma, K. H., Colwell, R. K., et al. (2014). Rarefaction and extrapolation with hill numbers: a framework for sampling and estimation in species diversity studies. Ecol. Monogr. 84, 45–67. doi: 10.1890/13-0133.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chelakkot, C., Choi, Y., Kim, D.-K., Park, H. T., Ghim, J., Kwon, Y., et al. (2018). Akkermansia muciniphila-derived extracellular vesicles influence gut permeability through the regulation of tight junctions. Exp. Mol. Med. 50:e450. doi: 10.1038/emm.2017.282

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, I., Yamanishi, S., Cox, L., Methé, B. A., Zavadil, J., Li, K., et al. (2012). Antibiotics in early life alter the murine colonic microbiome and adiposity. Nature 488, 621–626. doi: 10.1038/nature11400

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, M. S., Yu, J. S., Yoo, H. H., and Kim, D.-H. (2018). The role of gut microbiota in the pharmacokinetics of antihypertensive drugs. Pharmacol. Res. 130, 164–171. doi: 10.1016/j.phrs.2018.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, G., Sandhu, K. V., Griffin, B. T., Dinan, T. G., Cryan, J. F., and Hyland, N. P. (2019). Gut reactions: breaking down xenobiotic–microbiome interactions. Pharmacol. Rev. 71, 198–224. doi: 10.1124/pr.118.015768

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, S. F., Murphy, E. F., O’Sullivan, O., Lucey, A. J., Humphreys, M., Hogan, A., et al. (2014). Exercise and associated dietary extremes impact on gut microbial diversity. Gut 63, 1913–1920. doi: 10.1136/gutjnl-2013-306541

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, S. M. (1995). Human relevance of animal carcinogenicity studies. Regul. Toxicol. Pharmacol. 21, 75–80. doi: 10.1006/rtph.1995.1012

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, S. L., and Patterson, A. D. (2020). The gut microbiome: an orchestrator of xenobiotic metabolism. Acta Pharm. Sin. B 10, 19–32. doi: 10.1016/j.apsb.2019.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Coyte, K. Z., and Rakoff-Nahoum, S. (2019). Understanding competition and cooperation within the mammalian gut microbiome. Curr. Biol. 29, R538–R544. doi: 10.1016/j.cub.2019.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Coyte, K. Z., Schluter, J., and Foster, K. R. (2015). The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666. doi: 10.1126/science.aad2602

PubMed Abstract | CrossRef Full Text | Google Scholar

Currò, D. (2018). The role of gut microbiota in the modulation of drug action: a focus on some clinically significant issues. Expert Rev. Clin. Pharmacol. 11, 171–183. doi: 10.1080/17512433.2018.1414598

PubMed Abstract | CrossRef Full Text | Google Scholar

Daina, A., and Zoete, V. (2016). A boiled−egg to predict gastrointestinal absorption and brain penetration of small molecules. ChemMedChem 11:1117. doi: 10.1002/cmdc.201600182

PubMed Abstract | CrossRef Full Text | Google Scholar

Day-Walsh, P., Shehata, E., Saha, S., Savva, G. M., Nemeckova, B., Speranza, J., et al. (2021). The use of an in-vitro batch fermentation (human colon) model for investigating mechanisms of TMA production from choline, l-carnitine and related precursors by the human gut microbiota. Eur. J. Nutr. 60, 3987–3999. doi: 10.1007/s00394-021-02572-6

PubMed Abstract | CrossRef Full Text | Google Scholar

De Boeck, I., Spacova, I., Vanderveken, O. M., and Lebeer, S. (2021). Lactic acid bacteria as probiotics for the nose? Microb. Biotechnol. 14, 859–869. doi: 10.1111/1751-7915.13759

PubMed Abstract | CrossRef Full Text | Google Scholar

de Carvalho, N., Oliveira, D. L., Saleh, M. A. D., Pintado, M. E., and Madureira, A. R. (2021). Importance of gastrointestinal in vitro models for the poultry industry and feed formulations. Anim. Feed Sci. Technol. 271:114730. doi: 10.1016/j.anifeedsci.2020.114730

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Dethlefsen, L., and Relman, D. A. (2011). Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc. Natl. Acad. Sci. 108, 4554–4561. doi: 10.1073/pnas.1000087107

PubMed Abstract | CrossRef Full Text | Google Scholar

Donaldson, G. P., Lee, S. M., and Mazmanian, S. K. (2016). Gut biogeography of the bacterial microbiota. Nat. Rev. Microbiol. 14, 20–32. doi: 10.1038/nrmicro3552

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudek-Wicher, R. K., Junka, A., and Bartoszewicz, M. (2018). The influence of antibiotics and dietary components on gut microbiota. Przegląd Gastroenterol. 13, 85–92. doi: 10.5114/pg.2018.76005

PubMed Abstract | CrossRef Full Text | Google Scholar

Duncan, S. H., Louis, P., Thomson, J. M., and Flint, H. J. (2009). The role of pH in determining the species composition of the human colonic microbiota. Environ. Microbiol. 11, 2112–2122. doi: 10.1111/j.1462-2920.2009.01931.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ecale, F., El Houari, A., Crapart, S., Laparre, J., Ramnath, M., Berjeaud, J.-M., et al. (2021). In vitro sensitivity of 30 anaerobic bacterial strains of the human intestinal core microbiota to antibiotics: culture and LC-MS/MS approaches. Anaerobe 67:102314. doi: 10.1016/j.anaerobe.2020.102314

PubMed Abstract | CrossRef Full Text | Google Scholar

Ejtahed, H.-S., Angoorani, P., Hasani-Ranjbar, S., Siadat, S.-D., Ghasemi, N., Larijani, B., et al. (2018). Adaptation of human gut microbiota to bariatric surgeries in morbidly obese patients: a systematic review. Microb. Pathog. 116, 13–21. doi: 10.1016/j.micpath.2017.12.074

PubMed Abstract | CrossRef Full Text | Google Scholar

El Houari, A., Ranchou-Peyruse, M., Ranchou-Peyruse, A., Dakdaki, A., Guignard, M., Idouhammou, L., et al. (2017). Desulfobulbus oligotrophicus sp. nov., a sulfate-reducing and propionate-oxidizing bacterium isolated from a municipal anaerobic sewage sludge digester. Int. J. Syst. Evol. Microbiol. 67, 275–281. doi: 10.1099/ijsem.0.001615

PubMed Abstract | CrossRef Full Text | Google Scholar

Enright, E. F., Gahan, C. G., Joyce, S. A., and Griffin, B. T. (2016). Focus: microbiome: the impact of the gut microbiota on drug metabolism and clinical outcome. Yale J. Biol. Med. 89:375.

Google Scholar

Faith, J. J., Guruge, J. L., Charbonneau, M., Subramanian, S., Seedorf, H., Goodman, A. L., et al. (2013). The long-term stability of the human gut microbiota. Science 341:1237439. doi: 10.1126/science.1237439

PubMed Abstract | CrossRef Full Text | Google Scholar

Fallingborg, J. (1999). Intraluminal pH of the human gastrointestinal tract. Dan. Med. Bull. 46, 183–196.

Google Scholar

Foster, K. R., and Bell, T. (2012). Competition, not cooperation, dominates interactions among culturable microbial species. Curr. Biol. 22, 1845–1850. doi: 10.1016/j.cub.2012.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Francino, M. P. (2016). Antibiotics and the human gut microbiome: dysbioses and accumulation of resistances. Front. Microbiol. 6:1543.

Google Scholar

Frankel, A. E., Coughlin, L. A., Kim, J., Froehlich, T. W., Xie, Y., Frenkel, E. P., et al. (2017). Metagenomic shotgun sequencing and unbiased metabolomic profiling identify specific human gut microbiota and metabolites associated with immune checkpoint therapy efficacy in melanoma patients. Neoplasia 19, 848–855. doi: 10.1016/j.neo.2017.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Franzosa, E. A., Hsu, T., Sirota-Madi, A., Shafquat, A., Abu-Ali, G., Morgan, X. C., et al. (2015). Sequencing and beyond: integrating molecular “omics” for microbial community profiling. Nat. Rev. Microbiol. 13, 360–372. doi: 10.1038/nrmicro3451

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y., Zhang, J., Chen, K., Xiao, C., Fan, L., Zhang, B., et al. (2019). An in vitro fermentation study on the effects of Dendrobium officinale polysaccharides on human intestinal microbiota from fecal microbiota transplantation donors. J. Funct. Foods 53, 44–53. doi: 10.1016/j.jff.2018.12.005

CrossRef Full Text | Google Scholar

Gal, M., and Brazier, J. S. (2004). Metronidazole resistance in Bacteroides spp. carrying nim genes and the selection of slow-growing metronidazole-resistant mutants. J. Antimicrob. Chemother. 54, 109–116. doi: 10.1093/jac/dkh296

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodman, A. L., Kallstrom, G., Faith, J. J., Reyes, A., Moore, A., Dantas, G., et al. (2011). Extensive personal human gut microbiota culture collections characterized and manipulated in gnotobiotic mice. Proc. Natl. Acad. Sci. 108, 6252–6257. doi: 10.1073/pnas.1102938108

PubMed Abstract | CrossRef Full Text | Google Scholar

Gosalbes, M. J., Compte, J., Moriano-Gutierrez, S., Vallès, Y., Jiménez-Hernández, N., Pons, X., et al. (2019). Metabolic adaptation in the human gut microbiota during pregnancy and the first year of life. EBioMedicine 39, 497–509. doi: 10.1016/j.ebiom.2018.10.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Guindon, S., Dufayard, J.-F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010

PubMed Abstract | CrossRef Full Text | Google Scholar

Gumienna, M., Lasik, M., and Czarnecki, Z. (2011). Bioconversion of grape and chokeberry wine polyphenols during simulated gastrointestinal in vitro digestion. Int. J. Food Sci. Nutr. 62, 226–233. doi: 10.3109/09637486.2010.532115

PubMed Abstract | CrossRef Full Text | Google Scholar

Haiser, H. J., Seim, K. L., Balskus, E. P., and Turnbaugh, P. J. (2014). Mechanistic insight into digoxin inactivation by Eggerthella lenta augments our understanding of its pharmacokinetics. Gut Microbes 5, 233–238. doi: 10.4161/gmic.27915

PubMed Abstract | CrossRef Full Text | Google Scholar

Hemarajata, P., and Versalovic, J. (2013). Effects of probiotics on gut microbiota: mechanisms of intestinal immunomodulation and neuromodulation. Ther. Adv. Gastroenterol. 6, 39–51. doi: 10.1177/1756283X12459294

PubMed Abstract | CrossRef Full Text | Google Scholar

Hibbing, M. E., Fuqua, C., Parsek, M. R., and Peterson, S. B. (2010). Bacterial competition: surviving and thriving in the microbial jungle. Nat. Rev. Microbiol. 8, 15–25. doi: 10.1038/nrmicro2259

PubMed Abstract | CrossRef Full Text | Google Scholar

Hisada, T., Endoh, K., and Kuriki, K. (2015). Inter- and intra-individual variations in seasonal and daily stabilities of the human gut microbiota in Japanese. Arch. Microbiol. 197, 919–934. doi: 10.1007/s00203-015-1125-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Huse, S. M., Ye, Y., Zhou, Y., and Fodor, A. A. (2012). A core human microbiome as viewed through 16S rRNA sequence clusters. PLoS One 7:e34242. doi: 10.1371/journal.pone.0034242

PubMed Abstract | CrossRef Full Text | Google Scholar

Iizumi, T., Battaglia, T., Ruiz, V., and Perez Perez, G. I. (2017). Gut microbiome and antibiotics. Arch. Med. Res. 48, 727–734. doi: 10.1016/j.arcmed.2017.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwasaki, M., Akiba, Y., and Kaunitz, J. D. (2019). Recent advances in vasoactive intestinal peptide physiology and pathophysiology: focus on the gastrointestinal system. F1000Research 8: 1629. doi: 10.12688/f1000research.18039.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnsen, B. O., Handal, N., Meisal, R., Bjørnholt, J. V., Gaustad, P., and Leegaard, T. M. (2017). erm gene distribution among Norwegian Bacteroides isolates and evaluation of phenotypic tests to detect inducible clindamycin resistance in Bacteroides species. Anaerobe 47, 226–232. doi: 10.1016/j.anaerobe.2017.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J. S., Spakowicz, D. J., Hong, B.-Y., Petersen, L. M., Demkowicz, P., Chen, L., et al. (2019). Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis. Nat. Commun. 10:5029. doi: 10.1038/s41467-019-13036-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Jourova, L., Anzenbacher, P., and Anzenbacherova, E. (2016). Human gut microbiota plays a role in the metabolism of drugs. Biomed. Pap. 160, 317–326. doi: 10.5507/bp.2016.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Klaassen, C. D., and Cui, J. Y. (2015). Review: mechanisms of how the intestinal microbiota alters the effects of drugs and bile acids. Drug Metab. Dispos. 43, 1505–1521. doi: 10.1124/dmd.115.065698

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, R. L., Chrystal, E. J. T., Beaulieu, B. B., and Goldman, P. (1979). Acetamide—a metabolite of metronidazole formed by the intestinal flora. Biochem. Pharmacol. 28, 3611–3615. doi: 10.1016/0006-2952(79)90407-6

CrossRef Full Text | Google Scholar

Kofteridis, D. P., Maraki, S., Mixaki, I., Mantadakis, E., and Samonis, G. (2004). Impact of prolonged treatment with trimethoprim-sulfamethoxazole on the human gut flora. Scand. J. Infect. Dis. 36, 771–772. doi: 10.1080/00365540410020956a

CrossRef Full Text | Google Scholar

Koppel, N., Maini Rekdal, V., and Balskus, E. P. (2017). Chemical transformation of xenobiotics by the human gut microbiota. Science 356:eaag2770. doi: 10.1126/science.aag2770

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostic, A. D., Howitt, M. R., and Garrett, W. S. (2013). Exploring host-microbiota interactions in animal models and humans. Genes Dev. 27, 701–718. doi: 10.1101/gad.212522.112

PubMed Abstract | CrossRef Full Text | Google Scholar

Kovatcheva-Datchary, P., Egert, M., Maathuis, A., Rajilić-Stojanović, M., De Graaf, A. A., Smidt, H., et al. (2009). Linking phylogenetic identities of bacteria to starch fermentation in an in vitro model of the large intestine by RNA-based stable isotope probing. Environ. Microbiol. 11, 914–926. doi: 10.1111/j.1462-2920.2008.01815.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Koziolek, M., Schneider, F., Grimm, M. Modeβ, C., Seekamp, A., Roustom, T., et al. (2015). Intragastric pH and pressure profiles after intake of the high-caloric, high-fat meal as used for food effect studies. J. Control. Release 220(Pt A), 71–78. doi: 10.1016/j.jconrel.2015.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kump, P. K., Gröchenig, H.-P., Lackner, S., Trajanoski, S., Reicht, G., Hoffmann, K. M., et al. (2013). Alteration of intestinal dysbiosis by fecal microbiota transplantation does not induce remission in patients with chronic active ulcerative colitis. Inflamm. Bowel Dis. 19, 2155–2165. doi: 10.1097/MIB.0b013e31829ea325

PubMed Abstract | CrossRef Full Text | Google Scholar

Lau, J. T., Whelan, F. J., Herath, I., Lee, C. H., Collins, S. M., Bercik, P., et al. (2016). Capturing the diversity of the human gut microbiota through culture-enriched molecular profiling. Genome Med. 8, 1–10. doi: 10.1186/s13073-016-0327-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Bastard, Q., Ward, T., Sidiropoulos, D., Hillmann, B. M., Chun, C. L., Sadowsky, M. J., et al. (2018). Fecal microbiota transplantation reverses antibiotic and chemotherapy-induced gut dysbiosis in mice. Sci. Rep. 8, 1–11. doi: 10.1038/s41598-018-24342-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Blay, G., Rytka, J., Zihler, A., and Lacroix, C. (2009). New in vitro colonic fermentation model for Salmonella infection in the child gut. FEMS Microbiol. Ecol. 67, 198–207. doi: 10.1111/j.1574-6941.2008.00625.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K., and Knight, R. (2012). Diversity, stability and resilience of the human gut microbiota. Nature 489, 220–230. doi: 10.1038/nature11550

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, K., Mahbub, R., and Fox, J. G. (2015). Xenobiotics: interaction with the Intestinal Microflora. ILAR J. 56, 218–227. doi: 10.1093/ilar/ilv018

PubMed Abstract | CrossRef Full Text | Google Scholar

Maccaferri, S., Vitali, B., Klinder, A., Kolida, S., Ndagijimana, M., Laghi, L., et al. (2010). Rifaximin modulates the colonic microbiota of patients with Crohn’s disease: an in vitro approach using a continuous culture colonic model system. J. Antimicrob. Chemother. 65, 2556–2565. doi: 10.1093/jac/dkq345

PubMed Abstract | CrossRef Full Text | Google Scholar

Maier, L., Pruteanu, M., Kuhn, M., Zeller, G., Telzerow, A., Anderson, E. E., et al. (2018). Extensive impact of non-antibiotic drugs on human gut bacteria. Nature 555, 623–628. doi: 10.1038/nature25979

PubMed Abstract | CrossRef Full Text | Google Scholar

Marzorati, M., Vanhoecke, B., De Ryck, T., Sadabad, M. S., Pinheiro, I., Possemiers, S., et al. (2014). The HMI™ module: a new tool to study the host-microbiota interaction in the human gastrointestinal tract in vitro. BMC Microbiol. 14:133. doi: 10.1186/1471-2180-14-133

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsunaga, N., Shimizu, H., Fujimoto, K., Watanabe, K., Yamasaki, T., Hatano, N., et al. (2018). Expression of glyceraldehyde-3-phosphate dehydrogenase on the surface of Clostridium perfringens cells. Anaerobe 51, 124–130. doi: 10.1016/j.anaerobe.2018.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagy, E., Sóki, J., Urbán, E., Szõke, I., Fodor, E., and Edwards, R. (2001). Occurrence of metronidazole and imipenem resistance among Bacteroides fragilis group clinical isolates in Hungary. Acta Biol. Hung. 52, 271–280. doi: 10.1556/ABiol.52.2001.2-3.11

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, R., Bae, S. D. W., Zhou, G., Read, S. A., Ahlenstiel, G., George, J., et al. (2020). Application of organoids in translational research of human diseases with a particular focus on gastrointestinal cancers. Biochim. Biophys. Acta Rev. Cancer 1873:188350. doi: 10.1016/j.bbcan.2020.188350

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, V. D., Nguyen, T. T., Pham, T. T., Packianather, M., and Le, C. H. (1999). Molecular screening and genetic diversity analysis of anticancer Azurin-encoding and Azurin-like genes in human gut microbiome deduced through cultivation-dependent and cultivation-independent studies. Int. Microbiol. 22, 437–449. doi: 10.1007/s10123-019-00070-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicas, T. I., Wu, C. Y., Hobbs, J. N., Preston, D. A., and Allen, N. E. (1989). Characterization of vancomycin resistance in Enterococcus faecium and Enterococcus faecalis. Antimicrob. Agents Chemother. 33, 1121–1124. doi: 10.1128/aac.33.7.1121

PubMed Abstract | CrossRef Full Text | Google Scholar

Niestêpski, S., Harnisz, M., Korzeniewska, E., Aguilera-Arreola, M. G., Contreras-Rodríguez, A., Filipkowska, Z., et al. (2019). The emergence of antimicrobial resistance in environmental strains of the Bacteroides fragilis group. Environ. Int. 124, 408–419. doi: 10.1016/j.envint.2018.12.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Oda, S., Fujiwara, R., Kutsuno, Y., Fukami, T., Itoh, T., Yokoi, T., et al. (2015). Targeted screen for human UDP-glucuronosyltransferases inhibitors and the evaluation of potential drug-drug interactions with zafirlukast. Drug Metab. Dispos. 43, 812–818. doi: 10.1124/dmd.114.062141

PubMed Abstract | CrossRef Full Text | Google Scholar

Ono, S., Muratani, T., and Matsumoto, T. (2005). Mechanisms of resistance to imipenem and ampicillin in Enterococcus faecalis. Antimicrob. Agents Chemother. 49, 2954–2958. doi: 10.1128/AAC.49.7.2954-2958.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Panda, S., Khader, I. E., Casellas, F., Vivancos, J. L., Cors, M. G., Santiago, A., et al. (2014). Short-term effect of antibiotics on human gut microbiota. PLoS One 9:e95476. doi: 10.1371/journal.pone.0095476

PubMed Abstract | CrossRef Full Text | Google Scholar

Panek, M., Paljetak, H. È, Barešić, A., Perić, M., Matijašić, M., Lojkić, I., et al. (2018). Methodology challenges in studying human gut microbiota–effects of collection, storage, DNA extraction and next generation sequencing technologies. Sci. Rep. 8:5143. doi: 10.1038/s41598-018-23296-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Payne, A. N., Zihler, A., Chassard, C., and Lacroix, C. (2012). Advances and perspectives in in vitro human gut fermentation modeling. Trends Biotechnol. 30, 17–25. doi: 10.1016/j.tibtech.2011.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Burillo, S., Molino, S., Navajas-Porras, B., Valverde-Moya, ÁJ., Hinojosa-Nogueira, D., López-Maldonado, A., et al. (2021). An in vitro batch fermentation protocol for studying the contribution of food to gut microbiota composition and functionality. Nat. Protoc. 16, 3186–3209. doi: 10.1038/s41596-021-00537-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pham, V. T., and Mohajeri, M. H. (2018). The application of in vitro human intestinal models on the screening and development of pre- and probiotics. Benef. Microbes 9, 725–742. doi: 10.3920/BM2017.0164

PubMed Abstract | CrossRef Full Text | Google Scholar

Pompei, A., Cordisco, L., Raimondi, S., Amaretti, A., Pagnoni, U. M., Matteuzzi, D., et al. (2008). In vitro comparison of the prebiotic effects of two inulin-type fructans. Anaerobe 14, 280–286. doi: 10.1016/j.anaerobe.2008.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Privitera, G., Dublanchet, A., and Sebald, M. (1979). Transfer of multiple antibiotic resistance between subspecies of Bacteroides fragilis. J. Infect. Dis. 139, 97–101. doi: 10.1093/infdis/139.1.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Rafii, F., and Cerniglia, C. E. (1993). Comparison of the azoreductase and nitroreductase from Clostridium perfringens. Appl. Environ. Microbiol. 59, 1731–1734. doi: 10.1128/aem.59.6.1731-1734.1993

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2020). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Relman, D. A. (2012). The human microbiome: ecosystem resilience and health. Nutr. Rev. 70, S2–S9. doi: 10.1111/j.1753-4887.2012.00489.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rettedal, E. A., Gumpert, H., and Sommer, M. O. A. (2014). Cultivation-based multiplex phenotyping of human gut microbiota allows targeted recovery of previously uncultured bacteria. Nat. Commun. 5:4714. doi: 10.1038/ncomms5714

PubMed Abstract | CrossRef Full Text | Google Scholar

Rinninella, E., Raoul, P., Cintoni, M., Franceschi, F., Miggiano, G. A. D., Gasbarrini, A., et al. (2019). What is the healthy gut microbiota composition? a changing ecosystem across age, environment, diet, and diseases. Microorganisms 7:14. doi: 10.3390/microorganisms7010014

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.

Google Scholar

Saitta, K. S., Zhang, C., Lee, K. K., Fujimoto, K., Redinbo, M. R., and Boelsterli, U. A. (2014). Bacterial β-glucuronidase inhibition protects mice against enteropathy induced by indomethacin, ketoprofen or diclofenac: mode of action and pharmacokinetics. Xenobiotica 44, 28–35. doi: 10.3109/00498254.2013.811314

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneeberger, M., Everard, A., Gómez-Valadés, A. G., Matamoros, S., Ramírez, S., Delzenne, N. M., et al. (2015). Akkermansia muciniphila inversely correlates with the onset of inflammation, altered adipose tissue metabolism and metabolic disorders during obesity in mice. Sci. Rep. 5:16643. doi: 10.1038/srep16643

PubMed Abstract | CrossRef Full Text | Google Scholar

Schumacher, U. K., and Single, B. (1998). Susceptibility of the anaerobic gram-negative non-sporulating rod, Bilophila wadsworthia to beta-lactams, beta-lactamase inhibitors, meropenem, metronidaz. Zentralblatt Für Bakteriol. 287, 421–425. doi: 10.1016/S0934-8840(98)80179-5

CrossRef Full Text | Google Scholar

Selber-Hnatiw, S., Sultana, T., Tse, W., Abdollahi, N., Abdullah, S., Al Rahbani, J., et al. (2020). Metabolic networks of the human gut microbiota. Microbiology 166, 96–119. doi: 10.1099/mic.0.000853

PubMed Abstract | CrossRef Full Text | Google Scholar

Senghor, B., Sokhna, C., Ruimy, R., and Lagier, J.-C. (2018). Gut microbiota diversity according to dietary habits and geographical provenance. Hum. Microbiome J. 8, 1–9. doi: 10.1016/j.humic.2018.01.001

CrossRef Full Text | Google Scholar

Shah, P., Fritz, J. V., Glaab, E., Desai, M. S., Greenhalgh, K., Frachet, A., et al. (2016). A microfluidics-based in vitro model of the gastrointestinal human–microbe interface. Nat. Commun. 7, 1–15. doi: 10.1038/ncomms11535

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Z.-H., Zhu, C.-X., Quan, Y.-S., Yang, Z.-Y., Wu, S., Luo, W.-W., et al. (2018). Relationship between intestinal microbiota and ulcerative colitis: mechanisms and clinical application of probiotics and fecal microbiota transplantation. World J. Gastroenterol. 24, 5–14. doi: 10.3748/wjg.v24.i1.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Silverman, J. D., Durand, H. K., Bloom, R. J., Mukherjee, S., and David, L. A. (2018). Dynamic linear models guide design and analysis of microbiota studies within artificial human guts. Microbiome 6, 1–20.

Google Scholar

Singh, S. B., Young, K., and Silver, L. L. (2017). What is an “ideal” antibiotic? Discovery challenges and path forward. Biochem. Pharmacol. 133, 63–73. doi: 10.1016/j.bcp.2017.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Y., Liu, C., Lee, J., Bolaòos, M., Vaisanen, M.-L., and Finegold, S. M. (2005). Bacteroides goldsteinii sp. nov.” Isolated from clinical specimens of human intestinal origin. J. Clin. Microbiol. 43, 4522–4527. doi: 10.1128/JCM.43.9.4522-4527.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sousa, T., Paterson, R., Moore, V., Carlsson, A., Abrahamsson, B., and Basit, A. W. (2008). The gastrointestinal microbiota as a site for the biotransformation of drugs. Int. J. Pharm. 363, 1–25. doi: 10.1016/j.ijpharm.2008.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Stojanov, S., Berlec, A., and Štrukelj, B. (2020). The influence of probiotics on the Firmicutes/Bacteroidetes ratio in the treatment of obesity and inflammatory bowel disease. Microorganisms 8:1715. doi: 10.3390/microorganisms8111715

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, C., Chen, L., and Shen, Z. (2019). Mechanisms of gastrointestinal microflora on drug metabolism in clinical practice. Saudi Pharm. J. 27, 1146–1156. doi: 10.1016/j.jsps.2019.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, J., Kim, S., Cabatbat, J. J. T., Jang, S., Jin, Y.-S., Jung, G. Y., et al. (2017). Global metabolic interaction network of the human gut microbiota for context-specific community-scale analysis. Nat. Commun. 8, 1–12. doi: 10.1038/ncomms15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Tap, J., Mondot, S., Levenez, F., Pelletier, E., Caron, C., Furet, J.-P., et al. (2009). Towards the human intestinal microbiota phylogenetic core. Environ. Microbiol. 11, 2574–2584. doi: 10.1111/j.1462-2920.2009.01982.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Teng, L.-J., Hsueh, P.-R., Tsai, J.-C., Liaw, S.-J., Ho, S.-W., and Luh, K.-T. (2002). High incidence of cefoxitin and clindamycin resistance among anaerobes in Taiwan. Antimicrob. Agents Chemother. 46, 2908–2913. doi: 10.1128/AAC.46.9.2908-2913.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Tilkorn, F. K. M. T., Frickmann, H., Simon, I. S., Schwanbeck, J., Horn, S., Zimmermann, O., et al. (2020). Antimicrobial resistance patterns in Clostridioides difficile strains isolated from neonates in Germany. Antibiotics 9:481. doi: 10.3390/antibiotics9080481

PubMed Abstract | CrossRef Full Text | Google Scholar

Tobío, M., Sánchez, A., Vila, A., Soriano, I., Evora, C., Vila-Jato, J. L., et al. (2000). The role of PEG on the stability in digestive fluids and in vivo fate of PEG-PLA nanoparticles following oral administration. Colloids Surf. B Biointerfaces 18, 315–323. doi: 10.1016/S0927-7765(99)00157-5

CrossRef Full Text | Google Scholar

Tramontano, M., Andrejev, S., Pruteanu, M., Klünemann, M., Kuhn, M., Galardini, M., et al. (2018). Nutritional preferences of human gut bacteria reveal their metabolic idiosyncrasies. Nat. Microbiol. 3, 514–522. doi: 10.1038/s41564-018-0123-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsang, J. O.-L., Zhou, J., Zhao, X., Li, C., Zou, Z., Yin, F., et al. (2021). Development of three-dimensional human intestinal organoids as a physiologically relevant model for characterizing the viral replication kinetics and antiviral susceptibility of enteroviruses. Biomedicines 9:88. doi: 10.3390/biomedicines9010088

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsitko, I., Wiik-Miettinen, F., Mattila, O., Rosa-Sibakov, N., Seppänen-Laakso, T., Maukonen, J., et al. (2019). A small in vitro fermentation model for screening the gut microbiota effects of different fiber preparations. Int. J. Mol. Sci. 20:1925. doi: 10.3390/ijms20081925

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457, 480–484. doi: 10.1038/nature07540

PubMed Abstract | CrossRef Full Text | Google Scholar

Umigai, N., Murakami, K., Ulit, M. V., Antonio, L. S., Shirotori, M., Morikawa, H., et al. (2011). The pharmacokinetic profile of crocetin in healthy adult human volunteers after a single oral administration. Phytomedicine 18, 575–578. doi: 10.1016/j.phymed.2010.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Van den Abbeele, P., Grootaert, C., Marzorati, M., Possemiers, S., Verstraete, W., Gérard, P., et al. (2010). Microbial community development in a dynamic gut model is reproducible, colon region specific, and selective for Bacteroidetes and Clostridium cluster IX. Appl. Environ. Microbiol. 76, 5237–5246. doi: 10.1128/AEM.00759-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Venema, K., and Van den Abbeele, P. (2013). Experimental models of the gut microbiome. Best Pract. Res. Clin. Gastroenterol. 27, 115–126. doi: 10.1016/j.bpg.2013.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Venturelli, O. S., Carr, A. C., Fisher, G., Hsu, R. H., Lau, R., Bowen, B. P., et al. (2018). Deciphering microbial interactions in synthetic human gut microbiome communities. Mol. Syst. Biol. 14:e8157. doi: 10.15252/msb.20178157

PubMed Abstract | CrossRef Full Text | Google Scholar

Vich Vila, A., Collij, V., Sanna, S., Sinha, T., Imhann, F., Bourgonje, A. R., et al. (2020). Impact of commonly used drugs on the composition and metabolic function of the gut microbiota. Nat. Commun. 11, 1–11. doi: 10.1038/s41467-019-14177-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Villa, M. M., Bloom, R. J., Silverman, J. D., Durand, H. K., Jiang, S., Wu, A., et al. (2020). Interindividual variation in dietary carbohydrate metabolism by gut Bacteria revealed with droplet microfluidic culture. Msystems 5:e0864-19. doi: 10.1128/mSystems.00864-19

PubMed Abstract | CrossRef Full Text | Google Scholar

von Martels, J. Z., Sadabad, M. S., Bourgonje, A. R., Blokzijl, T., Dijkstra, G., Faber, K. N., et al. (2017). The role of gut microbiota in health and disease: in vitro modeling of host-microbe interactions at the aerobe-anaerobe interphase of the human gut. Anaerobe 44, 3–12. doi: 10.1016/j.anaerobe.2017.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

von Martels, J. Z. H. (2019). Towards Novel Biomarkers And Rational Nutritional Interventions In Inflammatory Bowel Disease. Groningen: University of Groningen.

Google Scholar

Vrancken, G., Gregory, A. C., Huys, G. R., Faust, K., and Raes, J. (2019). Synthetic ecology of the human gut microbiota. Nat. Rev. Microbiol. 17, 754–763. doi: 10.1038/s41579-019-0264-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, J., Griffin, B. T., Clarke, G., and Hyland, N. P. (2018). Drug–gut microbiota interactions: implications for neuropharmacology. Br. J. Pharmacol. 175, 4415–4429. doi: 10.1111/bph.14366

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Lv, Y., Cai, J., Schwarz, S., Cui, L., Hu, Z., et al. (2015). A novel gene, optrA, that confers transferable resistance to oxazolidinones and phenicols and its presence in Enterococcus faecalis and Enterococcus faecium of human and animal origin. J. Antimicrob. Chemother. 70, 2182–2190. doi: 10.1093/jac/dkv116

PubMed Abstract | CrossRef Full Text | Google Scholar

Wasielewski, H., Alcock, J., and Aktipis, A. (2016). Resource conflict and cooperation between human host and gut microbiota: implications for nutrition and health. Ann. N. Y. Acad. Sci. 1372, 20–28. doi: 10.1111/nyas.13118

PubMed Abstract | CrossRef Full Text | Google Scholar

Whittle, G., Shoemaker, N. B., and Salyers, A. A. (2002). The role of Bacteroides conjugative transposons in the dissemination of antibiotic resistance genes. Cell. Mol. Life Sci. 59, 2044–2054. doi: 10.1007/s000180200004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, I. D., and Nicholson, J. K. (2017). Gut microbiome interactions with drug metabolism, efficacy, and toxicity. Transl. Res. 179, 204–222. doi: 10.1016/j.trsl.2016.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yaffe, E., and Relman, D. A. (2020). Tracking microbial evolution in the human gut using Hi-C reveals extensive horizontal gene transfer, persistence and adaptation. Nat. Microbiol. 5, 343–353. doi: 10.1038/s41564-019-0625-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zihler, A., Gagnon, M., Chassard, C., Hegland, A., Stevens, M. J., Braegger, C. P., et al. (2010). Unexpected consequences of administering bacteriocinogenic probiotic strains for Salmonella populations, revealed by an in vitro colonic model of the child gut. Microbiology 156, 3342–3353. doi: 10.1099/mic.0.042036-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmermann, M., Zimmermann-Kogadeeva, M., Wegmann, R., and Goodman, A. L. (2019). Mapping human microbiome drug metabolism by gut bacteria and their genes. Nature 570, 462–467. doi: 10.1038/s41586-019-1291-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, L., Spanogiannopoulos, P., Pieper, L. M., Chien, H.-C., Cai, W., Khuri, N., et al. (2020). Bacterial metabolism rescues the inhibition of intestinal drug absorption by food and drug additives. Proc. Natl. Acad. Sci. 117, 16009–16018. doi: 10.1073/pnas.1920483117

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, T., Wong, S. H., Cheung, C. P., Lam, K., Lui, R., Cheung, K., et al. (2018). Gut fungal dysbiosis correlates with reduced efficacy of fecal microbiota transplantation in Clostridium difficile infection. Nat. Commun. 9, 1–11. doi: 10.1038/s41467-018-06103-6

PubMed Abstract | CrossRef Full Text | Google Scholar