- National Food Institute, Technical University of Denmark, Kongens Lyngby, Denmark
The intestinal gut microbiota is essential for maintaining host health. Concerns have been raised about the possible connection between antibiotic use, causing microbiota disturbances, and the increase in allergic and autoimmune diseases observed during the last decades. To elucidate the putative connection between antibiotic use and immune regulation, we have assessed the effects of the antibiotic amoxicillin on immune regulation, protein uptake, and bacterial community structure in a Brown Norway rat model. Daily intra-gastric administration of amoxicillin resulted in an immediate and dramatic shift in fecal microbiota, characterized by a reduction of within sample (α) diversity, reduced variation between animals (β diversity), increased relative abundance of Bacteroidetes and Gammaproteobacteria, with concurrent reduction of Firmicutes, compared to a water control group. In the small intestine, amoxicillin also affected microbiota composition significantly, but in a different way than observed in feces. The small intestine of control animals was vastly dominated by Lactobacillus, but this genus was much less abundant in the amoxicillin group. Instead, multiple different genera expanded after amoxicillin administration, with high variation between individual animals, thus the small intestinal α and β diversity were higher in the amoxicillin group compared to controls. After 1 week of daily amoxicillin administration, total fecal IgA level, relative abundance of small intestinal regulatory T cells and goblet cell numbers were higher in the amoxicillin group compared to controls. Several bacterial genera, including Escherichia/Shigella, Klebsiella (Gammaproteobacteria), and Bifidobacterium, for which the relative abundance was higher in the small intestine in the amoxicillin group than in controls, were positively correlated with the fraction of small intestinal regulatory T cells. Despite of epidemiologic studies showing an association between early life antibiotic consumption and later prevalence of inflammatory bowel diseases and food allergies, our findings surprisingly indicated that amoxicillin-induced perturbation of the gut microbiota promotes acute immune regulation. We speculate that the observed increase in relative abundance of small intestinal regulatory T cells is partly mediated by immunomodulatory lipopolysaccharides derived from outgrowth of Gammaproteobacteria.
Introduction
The gut associated lymphoid tissue is continuously exposed to a vast amount of bacterial and food derived antigens. Failure to develop tolerance toward these antigens may lead to inflammatory bowel diseases (IBD) or food allergies, both characterized by an adverse immune response directed against either commensal bacteria or dietary components. A healthy gut microbiota is important to avoid detrimental responses against harmless antigens (Stephen-Victor and Chatila, 2019), and epidemiological studies suggest that perturbation of the early life gut microbiota by pre- or post-natal exposure to antibiotics increases the risk of food allergies (Metsälä et al., 2013; Hirsch et al., 2017; Mitre et al., 2018) and autoimmune disorders including IBD (Hviid et al., 2011; Shaw et al., 2011; Kronman et al., 2012).
The gut microbiota promotes immune tolerance by stimulating the generation of regulatory T cells (Tregs) (Atarashi et al., 2011, 2013; Furusawa et al., 2013) and by reinforcing the intestinal barrier function, e.g., through regulation of tight junction permeability (Peng et al., 2009). Additionally, intestinal bacteria influence the secretion of mucus by goblet cells (Jung et al., 2015) and of immunoglobulin A (IgA) by plasma cells (Moreau et al., 1978; Shroff et al., 1995), both of which contribute to (I) antigen presentation and tolerogenic immune priming (Weltzin et al., 1989; McDole et al., 2012), (II) intestinal barrier function, and (III) balancing the microbial community (Martens et al., 2009; Donaldson et al., 2018). While multiple studies have assessed the influence of gut microbiota perturbation on immune regulation using broad-spectrum antibiotic cocktails (Ochoa-Repáraz et al., 2009; Hill et al., 2010; Nakamura et al., 2016) and germ-free (GF) animals (Ishikawa et al., 2008; Cahenzli et al., 2013), only few studies have investigated how administration of a single, clinically relevant antibiotic affects immune regulation (Tulstrup et al., 2015; Benakis et al., 2016).
Brown Norway (BN) rats constitute an well-established animal model in food allergy research (Abril-Gil et al., 2015; Jensen et al., 2019), as they mount an immune response resembling that of atopic human individuals (Knippels et al., 1999). The used BN rats originate from our in-house breeding colony, and in contrast to most animals from commercial vendors they have a natural microbiota and may thereby better resemble the human immune response (Rosshart et al., 2019). Amoxicillin is a bactericidal β-lactam antibiotic that affects both Gram-positive and Gram-negative bacteria by inhibiting an essential enzyme for bacterial cell wall synthesis (Bodey and Nance, 1972). It is the most widely used penicillin in Europe (European Centre for Disease Prevention and Control, 2018), and is frequently prescribed for treatment of pediatric infections (Dolk et al., 2018), and the impact of microbiota alterations caused by this antibiotic on early development of the immune system is therefore highly relevant to address. The gut microbiota of BN rats has not previously been analyzed by NGS based methods. In order to elucidate the interplay between microbiota composition and immune regulation, we characterized the temporal and spatial effects of 7 days of amoxicillin administration in BN rats by analyzing the gut microbiota, host immune regulation and intestinal permeability to cow’s milk protein. Interestingly, the observed rapid amoxicillin-driven perturbation of the gut microbiota was associated with an acute immune regulatory response.
Materials and Methods
Animals
BN rats from the in-house breeding colony (National Food Institute, Technical University of Denmark, Denmark) were housed in Makrolon cages and kept at a 12 h light:dark cycle, at a temperature of 22 ± 1°C and a relative humidity of 55 ± 5%. Rats were observed twice daily and clinical signs recorded. The rats were fed a milk-free diet for ≥10 generation that was produced in-house and based on rice flour, potato protein and fish meal as protein sources, as previously described (Bøgh et al., 2009), with the exception of maize flakes being substituted with rice flour. Diet and water were given ad libitum. Ethical approval was given by the Danish Animal Experiments Inspectorate and the authorization number given 2015-15-0201-00553-C1. The experiment was overseen by the National Food Institute’s in-house Animal Welfare Committee for animal care and use.
Animal Experiment
BN rats of both sex were allocated into two groups of 12 rats. Starting from 4–5 weeks of age, rats were daily gavaged with 30 mg amoxicillin in 0.5 mL milli Q water or water alone from days 0 to 6. Fecal samples were collected 3 days before first gavage (day −3) and daily throughout the study and stored at −80°C for microbiota analysis. From days −3 and 7 one additional fecal sample was collected and stored on ice for total IgA quantification in fecal water. Intestinal uptake of β-lactoglobulin (BLG) from cow’s milk was assessed by dosing animals with 100 mg of whey protein concentrate (kindly provided by Arla Foods Ingredients) in 1 mL milli Q water by oral gavage 15 min prior to sacrifice.
Dissection
Animals were euthanized by exsanguination using carbon dioxide inhalation as anesthesia. Whole blood was collected in sodium-heparin coated tubes (BD Bioscience, Franklin Lakes, NJ, United States) that were mixed for 10 min at room temperature (RT) and then stored on ice until processing immediately after the dissection, and in non-coated tubes stored 1 h at RT for preparation of serum. The intestine was excised and divided into multiple fractions for different analyses: mesenteric lymph nodes (mLNs) were stored in phosphate buffered saline (PBS, 137 mmol/L NaCl, 2.7 mmol/L KCl, 10 mmol/L Na2HPO4, and 1.8 mmol/L KH2PO4, pH 7.2) on ice, and processed for flow cytometry immediately after the dissection. From the first 20 cm of the small intestine (7 cm distal from the stomach), content was snap frozen in liquid nitrogen and stored at −80°C for microbiota analysis. The tissue was rinsed with 0.9% (w/v) NaCl before Peyer’s patches (PPs) and epithelium was removed, and the tissue was snap frozen in liquid nitrogen and stored at −80°C until analysis of permeability. PPs from the rest of the small intestine were excised, washed in 0.9% (w/v) NaCl and stored in PBS on ice until processing immediately after the dissection. The next 10 cm of the small intestine (without PPs) was rinsed with 0.9% (w/v) NaCl and stored in wash buffer [PBS with 2% (v/v) heat-inactivated fetal calf serum (FCS) and 15 mM HEPES] on ice until processing immediately after the dissection. The next 1 cm piece was rinsed with 0.9% (w/v) NaCl and stored in RNAlater (Invitrogen, Carlsbad, CA, United States) at −20°C until analysis of gene expression. Finally, 1 cm pieces of the small intestine and of colon were fixed overnight in 4% (w/v) paraformaldehyde for histology. Cecum was weighted, and cecum content was snap frozen in liquid nitrogen and stored at −80°C for microbiota analysis.
Bacterial DNA Extraction
DNA was extracted from up to 200 mg of feces and cecum content and up to 250 mg of small intestine content by DNeasy PowerLyzer PowerSoil Kit (Qiagen, Hilden, Germany) according to the manufacture’s protocol. Mechanical lysis of bacteria was conducted at 30 cycles/s twice for 5 min on bead beater MM300 (Retsch, VWR, Haan, Germany). DNA concentrations were measured by the Qubit ds DNA BR kit (Life Technologies, Carlsbad, CA, United States) and DNA was stored at −20°C.
Bacterial Load
Bacterial load was estimated by quantification of 16S rRNA gene copies by qPCR as previously described (Tulstrup et al., 2015). Briefly, the V3-region of the 16S rRNA gene was amplified in triplicate for each sample, using universal primers HDA1 and HDA2 (Walter et al., 2000). Each qPCR reactions contained 5.5 μL LightCycler1 480 II SYBR Green I Master (Roche Applied Science, Penzberg, Germany), 0.2 μM of each primer and 0.2 μL of diluted template DNA in a total reaction volume of 11 μL. Reaction conditions were: Pre-incubation at 95°C for 5 min followed by 45 cycles of 95°C for 10 s, 60°C for 15 s and 72°C for 45 s. Lastly, a melting curve was generated (95°C for 5 s, 68°C for 1 min and increasing the temperature to 98°C with a rate of 0.11°C/s with continuous fluorescence detection). The qPCR was run in 384-well format on a LightCycler® 480 II (Roche Applied Science) and analyzed using the LightCycler® 480 software. Tenfold dilutions of a linearized (SphI-digested) plasmid standard, construction by cloning the 199bp V3-region of the 16S rRNA gene of Escherichia coli (ATCC 25922) into the pCR14Blunt-TOPO vector (Invitrogen), was used for quantification of 16S rRNA genes.
Amplicon Sequencing of the 16S rRNA Gene
The bacterial community composition was analyzed by sequencing of the V3-region of the 16S rRNA gene in the extracted DNA. Amplification of the V3-region and subsequent sequencing was performed using the Ion Torrent PGM platform (Life Technologies) as previously described (Christensen et al., 2014). In short, the V3-region of the 16S rRNA gene was amplified using a universal forward primer (PBU 5′-A- adapter-TCAG-barcode-CCTACGGGAGGCAGCAG-3′) with a unique 10–12 bp barcode for each sample (Ion Xpress barcode as suggested by the supplier, Life Technologies) and a universal reverse primer (PBR 5′-trP1-adapter-ATTACCGCGGCTGCTGG-3′). The PCR reactions were conducted with 4 μL HF-buffer, 0.4 μL dNTP (10 mM of each base), 1 μM forward primer, 1 μM reverse primer, 5 ng template DNA in 1 μL, and 0.2 μL Phusion High-Fidelity DNA polymerase (Thermo Fisher Scientific, Waltham, MA, United States) in a total reaction volume of 20 μL. Reaction conditions were as follows: Initial 98°C for 30 s followed by 24 (feces and cecum content) or 30 (small intestine content) cycles of 98°C for 15 s and 72°C for 30 s, finally, 72°C for 5 min before cooling to 4°C. PCR products were purified by HighPrepTM PCR Clean-up System (Magbio, Gaithersburg, MD, United States) according to the manufacture’s protocol, and DNA concentrations were determined with Qubit HS assay. Finally, a library was constructed by mixing an equal amount of PCR products from each sample. Sequencing was performed on a 318-chip for Ion Torrent sequencing using the Ion OneTouchTM 200 Template Kit v2 DL.
Sequence Data Handling
Sequence data was obtained in FASTQ format and further processed using CLC bio genomic workbench (Qiagen) in order to de-multiplex samples and remove sequencing primers. Further quality trimming using default settings (remove low quality nucleotides pbase-calling error = 0.05, trim ambiguous nucleotides = 2) and filtering only reads with a final length between 125–180 bp were exported in FASTQ format. Further quality trimming was performed in Divisive Amplicon Denoising Algorithm 2 (DADA2, Version 1.10.1) with default settings as described elsewhere (Callahan et al., 2016). Finally, an amplicon sequencing variant (ASV) table was constructed which contains the counts of each sequence variant in each sample. All sequence reads were taxonomically classified using the Ribosomal Database Project Multi-classifier tool (Wang et al., 2007).
The ASV table was imported into the “Quantitative Insights Into Microbial Ecology” (Qiime) 2 pipeline (Version 2019.1), and α and β diversity metrics were calculated by the function “diversity core-metrics-phylogenetic” based on a rooted phylogenetic tree. For analysis of fecal samples only, samples were rarefied to 20,000 reads to eliminate bias from uneven sampling depth, and for analyses that also included small intestine, samples were rarefied to 10,000 reads.
Preparation of Serum and Fecal Water
To prepare serum, blood samples were allowed to coagulate for 1 h at RT and subsequently overnight at 4°C. The following day the coagulated blood was removed and samples were centrifuged at 1,800 g for 20 min at 4°C. The supernatants were transferred to clean tubes and stored at −20°C until analysis.
To prepare fecal water, fecal pellets were mixed with 10 μL cold PBS with 0.05% (w/v) NaN3 (Sigma, St. Louis, MO, United States) per mg sample in a bead beater MM300 (Retsch) for 15 min at 30 cycles/s followed by centrifugation at 16,000 g at 4°C for 10 min. The supernatants were transferred to clean tubes and stored at −20°C until analysis.
Sandwich ELISA for Detection of Total IgA
For detection of total IgA antibodies sandwich ELISA was performed. Plates (96 wells MaxiSorp, NUNC, Roskilde, Denmark) were coated with 100 μL/well of mouse anti-rat IgA (MCA191, Bio-Rad, Oxford, United Kingdom) diluted 1:2,000 in carbonate buffer (15 mm Na2CO3, 35 mm NaHCO3, pH 9.6) and incubated overnight at 4°C. To block unspecific binding 200 μL/well of 3% (w/v) egg white protein (E0500, Sigma) in PBS with 0.01% (w/v) Tween 20 (PBS-T) was added and plates were incubated at 37°C for 1 h. Next, 50 μL/well of twofold serial dilutions of serum or fecal water samples and positive and negative control samples diluted in PBS-T were added, and plates were incubated for 1 h at RT. For detection, plates were incubated with 100 μL/well of horseradish peroxidase-conjugated goat anti-rat IgA (STAR111P, Bio-Rad, United Kingdom) diluted 1:10,000 in PBS-T. Between each step, plates were washed five times in PBS-T. After the last wash, plates were additionally washed twice in running tap water. The reaction was visualized by adding 100 μL/well of 3,3′,5,5-tetramethylbenzidine (TMB)-one substrate (Kem-En-Tec, Taastrup, Denmark) for 12 min and stopped with 100 μL/well of 0.2 M H2SO4. Absorbance was measured at 450 nm with a reference wavelength of 630 nm using a microtiter reader (Gen5, BioTek Instruments, Winooski, VT, United States). Antibody titers were expressed as the log2 titer values and defined as the interpolated dilution of the given serum sample leading to the mean absorbance for the negative control +3 standard deviation.
Single Cell Suspensions of Small Intestinal and Lymphoid Tissues
Small intestine lamina propria (LP) samples stored on ice were opened longitudinally and cut into 5 mm pieces. The pieces were washed five times in wash buffer and tissue was digested in RPMI medium (Sigma) with 10% (v/v) FCS (Sigma), 15 mM HEPES (Sigma), 100 units/mL penicillin (Sigma), 100 μg/mL streptomycin (Sigma), 250 μg/mL collagenase (Sigma), and 1 mg/mL dispase II (Sigma) for 45 min on a shaker at 37°C. Afterward, the tissue digest was filtered through a 70 μm cell strainer, which was regularly rinsed with FACS buffer [PBS with 2% (v/v) FCS and 0.05% (v/v) NaN3]. PPs and mLNs were minced through 70 μm cell strainers, which were regularly rinsed with FACS buffer. The cell suspensions were centrifuged at 400 g (LP and mLNs) or 600 g (PP), 4°C for 10 min. The supernatants were discarded, and the cells re-suspended in 1 mL FACS buffer, and cells were counted in NucleoCassettes by use of a NucleoCounter (Chemometec, Allerød, Denmark).
Staining of Lymphocytes for Flow Cytometry
Approximate 106 cells/well were plated in 96-well plates. Non-specific binding was prevented by incubation 5 min on ice with 50 μL blocking solution with 10% (v/v) rat serum and 5 μg/mL anti-CD32 (D34-485, BD Biosciences) in FACS buffer. Cells were surface stained by addition of 50 μL/well antibody cocktail containing 2 μg/mL of each of the following: anti-B220-PE (HIS24, BD Biosciences), anti-CD3-PerCp-Cy5 (eBioG4.18, Thermo Fisher Scientific), anti-CD4-PE-Cy7 (OX-35, BD Biosciences), anti-CD45-APC-Cy7 (OX-1, Thermo Fisher Scientific), anti-CD25-BV421 (OX-39, BD Biosciences) in FACS buffer, and incubated 20 min on ice. Cells were washed with 100 μL FACS buffer and plates were centrifuged at 400 g, 4°C for 5 min. Supernatants were discarded and pellets resuspended in 200 μL FACS buffer.
Whole blood samples were diluted two times in FACS buffer and 100 μL were stained with 50 μL of the above described antibody cocktail and incubated for 20 min on ice. Afterward, 1 mL VersaLyse (Beckman Coulter, Brea, CA, United States) was added to lyse the red blood cells. The samples were incubated in dark at RT for 10 min before 2 mL FACS buffer was added and samples were centrifuge at 350 g, 4°C for 6 min. Supernatants were discarded and pellets resuspended in 200 μL FACS buffer and transferred to 96-well plates.
For identification of the transcription factors Forkhead box P3 (FoxP3) and Helios, cells were subsequently stained by using the BD Transcription Factor buffer set (BD Biosciences) according to the manufacture’s protocol. After fixation and permeabilization, the cells were blocked as described above. Cells were stained with 50 μL of 5 μg/mL anti-FoxP3-FITC (JFK-16s, Thermo Fisher Scientific) and 2 μg/mL anti-Helios- Alexa Fluor 647 (22F6, BD Biosciences) in TF Perm/Wash Buffer from the BD Transcription Factor buffer set, and incubated for 40 min on ice. After two washes in TF Perm/Wash Buffer, the cells were resuspended in FACS buffer. The stained cells were stored at 4°C until data was acquired. Just before accusation, 10 μL CountBright counting beads (Thermo Fisher Scientific) was added to each well. Data was acquired on BD FACSCanto II system (BD Bioscience) and analyzed by FlowJo (Version 10.4.2, BD Bioscience).
Intestinal Protein Uptake
Total proteins were extracted from small intestine tissue samples by mixing with 10 μL tris-lysis buffer [150 mM NaCl, 20 mM Tris, 1 mM EGTA, 1% (v/v) Triton X-100, 1 mM EDTA] with 2% (v/v) Halt protease inhibitor cocktail (Thermo Fisher Scientific) per mg tissue. One stainless steel bead (Qiagen) were added to each sample, and samples were homogenized in buffer by TissueLyser II (Qiagen) at 30 cycles/s for 2 min. Samples were incubated on ice for 20 min and mixed by vortexing every 5 min. Samples were centrifuged at 15,000 g for 20 min at 4°C and the supernatants were frozen at −80°C until analysis. Concentrations of the cow’s milk protein BLG were quantified in tissue extracts and serum samples by a commercial bovine BLG ELISA kit (Bethyl Laboratories, Montgomery, AL, United States) in 96 wells MaxiSorp plates (NUNC) according to the manufacture’s protocol with the exception that plates were coated overnight. The concentrations were determined in triplicates and calculated from the duplicated standard curves generated for each plate.
Histology
After overnight fixation in 4% formalin (CellPath, Newtown, United Kingdom), sections of small intestine and colon were dehydrated in a graded ethanol series from 77 to 99% ethanol (VWR Chemicals, Radnor, PS, United States). Xylene (VWR Chemicals) was used as clearing agent to replace the ethanol before the tissues were embedded in paraffin (Hounisen, Skanderborg, Denmark). Histological sections of 2 μm were stained with Hematoxylin (Ampliqon, Odense, Denmark) and Eosin (Merck, Darmstadt, Germany) to identify eosinophils, and Periodic acid–Schiff (PAS) (Merck) to identify goblet cells. The slides were examined using a Leica DMR upright microscope (Leica Microsystems GmbH, Wetzlar, Germany) and the software ImagePro Plus 7.0 (Media Cybernetics, Rockville, MD, United States) was used for images and measurements. Villus length was measured from the villus tip to the crypt-villus junction and in colon crypt depth was measured, with three villi or crypts measured/counted per animal. Analyses of histological sections were performed blinded.
Tissue RNA Extraction, cDNA Synthesis, and RT-qPCR
Approximately 20–40 mg small intestine tissue, stored in RNAlater, was homogenized using TissueLyser II (Qiagen) and total RNA was extracted with RNeasy Mini Kit (Qiagen) with on-column DNase digestion (RNase Free DNase Kit, Qiagen) according to manufacturer’s protocol. RNA quality and quantity were assessed using NanoDrop. cDNA was synthesized from 500 ng RNA with Omniscript RT Kit (Qiagen) in accordance to manufacture’s protocol and in addition of Random Primer Mix (BioNordika, Herlev, Denmark) and Anti-RNAse (Ambion, Life Technologies, Carlsbad, CA, United States). Quantitative RT-PCR (RT-qPCR) reactions were run in technical duplicates using a Quantstudio 7 Flex Real Time PCR system (Applied Biosystems, Thermo Fisher Scientific) in 11 μL reactions containing 3 μL diluted cDNA (1:20), TaqMan Fast Universal PCR Master Mix (2×) (Applied Biosystems), and TaqMan gene Expression Assay (Applied Biosystems). TaqMan gene assays used were Cdh1 (Cadherin-1 Rn00580109_m1), Ocln (Occludin Rn00580064_m1), Cldn2 (Claudin 2 Rn02063575_s1), Tjp1 (Tight junction protein 1 (ZO-1) Rn02116071_s1), TSLP (Thymic stromal lymphopoietin Rn01761072_m1), Il33 (Interleukin 33 Rn01759835_m1), and Muc2 (Mucin 2 Rn01498206_m1). Reaction conditions were: An initial cycle at 95°C for 20 s followed by 45 two-step thermal cycles at 95°C for 1 s and at 60°C for 20 s. The relative gene expression was calculated using the 2–ΔCT method using B2m (Beta-2-microglobulin Rn00560865_m1) and Sdha (Succinate dehydrogenase complex Rn00590475_m1) as normalization genes. Data was acquired with Quantstudio 7 Flex software (Applied Biosystems, Foster City, CA, United States).
Statistics
Differences between β diversity of groups were assessed by applying ANalysis Of SImilarities ANOSIM (Clarke, 1993) to weighted and unweighted UniFrac distances, and within group dispersion was analyzed with PERMDISP (Anderson et al., 2006) both in Qiime 2. Differential abundant genera between the amoxicillin and control group were determined by ANalysis of COMposition of Microbiomes (ANCOM) (Mandal et al., 2015) in R. This was only applied to small intestine samples, because more than 25% of genera varied between the two groups in cecum and feces. Graphs and additionally statistically analyses were generated in Prism version 8.1.1 (GraphPad, San Diego, CA, United States). Difference between initial and later time points within the same group were analyzed by one-way ANOVA followed by Dunnett’s multiple comparisons test. Difference between the two groups at the same time point were analyzed by t-test (when data passed D’Agostino-Pearson normality test) or non-parametric Mann–Whitney test (when not normally distributed). Non-parametric Spearman correlations were calculated between all pairs of relative abundance of small intestine genera (present in minimum 10 animals) and selected host response parameters. When indicated, p-values were false discovery rate (FDR) corrected by a two-stage sharpened method (Benjamini et al., 2006).
Results
Amoxicillin Transiently Reduces Fecal Bacterial Load and Persistently Decreases Fecal Microbial Diversity
BN rats were gavaged daily with either amoxicillin or water for 1 week, and the temporal effects on fecal bacterial load and composition were followed. The 1st day after onset, a statistically significant reduction of total bacterial load was observed in the amoxicillin group as compared to both the baseline level (day −3) and to the level in the control group (Figure 1A). The effect of amoxicillin on bacterial load was transient since no differences were observed compared to the initial level nor to the control group from day 2 and onward. In order to assess the overall effect of amoxicillin on gut microbiota, the within sample (α) diversity was estimated by the Faith’s phylogenetic diversity index. A statistically significant reduction in diversity was observed in the amoxicillin group as compared to the initial level from day 1, as well as between the amoxicillin and control group from day 2 (Figure 1B).
Figure 1. Temporal effects of amoxicillin on fecal bacterial load and composition. (A,B) Total bacterial load (A) and Faith’s phylogenetic diversity index (B) of control (CTR, black) and amoxicillin (AMX, red) groups. (C) Mean relative abundance of bacterial classes. (D,E) Principal coordinate analysis (PCoA) plots of unweighted UniFrac distances in the AMX (D) and CTR (E) groups colored according to sample day. (F) Heat-map of pairwise comparisons between days within the control group (upper right corner) and within the amoxicillin group (lower left corner) identified by Analysis of similarities (ANOSIM) between unweighted UniFrac distances. Intensity of the color indicates the value of the R statistics, where intense colors indicate large differences. Statistically significant differences between CTR and AMX on the same day analyzed by multiple t-test and FDR-correction (A,B) or between days analyzed by ANOSIM (F) are indicated as *p < 0.05, **p ≤ 0.01, ***p ≤ 0.001, and statistically significant differences relative to the initial level in the same group analyzed by one-way ANOVA followed by Dunnett’s multiple comparisons test (A,B) are indicated as ##p ≤ 0.01, ###p ≤ 0.001.
Amoxicillin Rapidly Disturbs Fecal Microbiota and Promotes Gammaproteobacteria Outgrowth
From day 1, the mean relative abundance of Firmicutes and Actinobacteria were lower in the amoxicillin group than in the control group, while Proteobacteria, especially Gammaproteobacteria were higher (Figure 1C and Supplementary Figures S1C, S2). In the amoxicillin group, the fecal microbiota continued changing over the course of the study, with the relative abundance of Bacteroidetes gradually expanding over time. This continuous shift in composition was also observed by principal coordinate analysis of unweighted UniFrac distances measuring the overall diversity between samples, while the microbiota compositions remained relatively stable in the control group (Figures 1D,E and Supplementary Figures S1A,B,D,E). These observations were confirmed by pairwise comparisons of unweighted UniFrac distances between different days within the same group, which, however, also revealed a minor temporal effect in the control group (Figure 1F).
Amoxicillin Affects Small Intestine and Cecal/Fecal Microbiota Compositions Differently
After 7 days of amoxicillin administration, animals were euthanized and the effects on intestinal microbiota in different intestinal compartments were analyzed. Amoxicillin was found to have a profound physiological effect on the ceca, which were dramatically enlarged (Figure 2A), indicating changes in microbiota composition and activity. Indeed, similar to the fecal microbiota, the Faith’s phylogenetic diversity of the cecal microbiota was much lower in the amoxicillin group compared to the control group on day 7 (Figure 2B). In contrast, the small intestine diversity was higher in the amoxicillin group compared to the control group.
Figure 2. Effects of amoxicillin on intestinal microbiota. (A,B) Cecum weight (A) and Faith’s phylogenetic diversity (B) of control (CTR, black) and amoxicillin (AMX, red) groups. (C) Relative abundance of bacterial classes in individual animals. (D,E) Principal coordinate analysis (PCoA) plot of weighted (D) and unweighted (E) and UniFrac distances between samples from small intestine (SI) (squares), cecum (triangles), and fecal (circles) content day 7 colored according to CTR and AMX treatment. (F) Relative abundance of bacterial genera in individual animals indicated by the cell color ranging from highly abundant genera in red and less abundant genera in blue. Genera for which the abundance was significantly different between CTR and AMX groups in the small intestine are highlighted with arrows marking whether the genera were higher or lower after amoxicillin administration. Statistically significant differences between groups analyzed by non-parametric Mann–Whitney test (A) or t-tests (B) are indicated as *p ≤ 0.05, ***p ≤ 0.001.
The distribution of bacterial classes was found to vary between the small intestinal, cecal and fecal compartments within the control group at day 7 (Figure 2C). The small intestine was dominated almost exclusively by Bacilli (lactobacilli), whereas three different classes dominated the cecal (Bacteroidia, Clostridia, and Deltaproteobacteria) and fecal (Bacteroidia, Clostridia, and Bacilli) microbiota. Amoxicillin administration was found to dramatically, and differently alter the microbiota in the different intestinal compartments. In the small intestine, the relative abundance of Bacilli was reduced while Gammaproteobacteria, Clostridia, and Actinobacteria were increased in all animals, but to a varying degree, resulting in a high variation between individual animals in the amoxicillin group. In contrast, the cecal and fecal microbiota composition were found to be uniform among the individual animals in the amoxicillin group, and the distribution of bacterial classes remarkably similar between the cecal and fecal microbiota in the amoxicillin group. In feces and cecum, the relative abundances of Clostridia and Bacilli were lower in the amoxicillin group than in the control group, while the relative abundance of Bacteroidia and Gammaproteobacteria were increased.
In line with the dramatic effect of amoxicillin on bacterial class distribution, the overall bacterial community structure (β diversity) was significantly affected by amoxicillin in both the small intestine, cecum and feces as determined by the ANOSIM method (Figures 2D,E). Furthermore, analysis of sample dispersion by PERMDISP of weighted UniFrac distances confirmed that the individual animals’ small intestinal samples were significantly more dispersed in the amoxicillin group than seen in the controls, while in contrast, cecal and fecal samples were significantly less dispersed in the amoxicillin group than in controls (Figure 2D). In summary, amoxicillin reduced within sample (α) and between samples (β) diversity with regard to the overall bacterial composition (weighted UniFrac distances) in fecal and cecal samples, but increased α and β diversity in small intestine samples.
Opposite, analysis of unweighted UniFrac distances based on presence/absence status of ASVs and not relative abundances revealed that the cecal and fecal samples from individual animals were significantly more dispersed in the amoxicillin group than in controls (Figure 2E). This suggests that while amoxicillin caused the amount of dominating species to become more similar within the group, the presence/absence status of rare species became more dissimilar within the group.
Statistical analysis of compositional changes (ANCOM method) of the microbiota caused by amoxicillin was performed only for the small intestinal compartment due to the very pronounced effects in the cecum and feces. The abundance of Lactobacillus was significantly lower in the amoxicillin group, as were Veillonella and Rothia (Figure 2F and Supplementary Table S1). A significant outgrowth of Escherichia/Shigella and Klebsiella (Gammaproteobacteria) in response to amoxicillin was also observed in the small intestine similar to the cecal and fecal compartments. Furthermore, the relative abundances of the Clostridial genera Romboutsia and Clostridium sensu stricto were significantly higher within the small intestine in the amoxicillin group than in the control group.
Amoxicillin Promotes Host Mucosal Immune Regulation
The effects of amoxicillin on humoral and mucosal immune regulation were analyzed after 7 days of administration. This revealed that total fecal IgA levels were higher in the amoxicillin group than in the control group, while total serum IgA levels were lower (Figure 3A). Various lymphocyte populations were analyzed in blood, mLNs, small intestine PPs and LP (Supplementary Figures S3A–C). Of these, the relative amount of CD25+FoxP3+ Tregs were significantly higher in small intestine LP (both CD4+ and CD4– Tregs) and PPs (only CD4– Tregs) in the amoxicillin group compared to the control group, while amoxicillin had no effect on Tregs in blood and mLNs (Figures 3B,C). Small intestinal uptake of the cow’s milk protein BLG was assessed by quantifying BLG in gut tissue and serum 15 min after oral gavage with whey protein. This revealed a tendency for higher uptake in the amoxicillin group (p = 0.183 and p = 0.139 for LP and serum, respectively; Figure 3D). Finally, mucus-producing goblet cells were quantified in histological slides (representative pictures in Supplementary Figure S4). Goblet cell numbers in both the small intestine and colon were elevated in the amoxicillin group compared to controls (Figure 3E). Relative expression of various immune related host genes were assessed in small intestine LP (Supplementary Figure S3D), and a tendency (p = 0.131) for higher expression of the Muc2 gene encoding mucin 2 was observed in the amoxicillin group (Figure 3F). A strong positive correlation between small intestine goblet cells and Muc2 gene expression supported these findings (Figure 3G).
Figure 3. Host responses to amoxicillin. (A) Total IgA levels in feces and serum. (B,C) Relative amounts of CD25+ FoxP3+ T regulatory cells of CD4+ T helper cells (B) and CD4– T cells (C) in blood, and small intestine mesenteric lymph nodes (mLN), Peyer’s patches (PPs) and lamina propria (LP). Cells in LP are plotted according to the right y-axis while cell from the remaining tissues are plotted according to the left y-axis. (D) β-lactoglobulin (BLG) concentration in small intestine LP (left y-axis) and serum (right y-axis). (E) Average goblet cells per villi or crypt in histological slides of small intestine (left y-axis) and colon (right y-axis). (F) Expression of the mucin 2 (Muc2) gene relative to two housekeeping genes. (G) Correlation between number small intestine goblet cells and relative expression of Muc2. Each symbol represents one animal and horizontal lines indicate mean (A,E,F) or median (B–D). Statistically significant differences between control (CTR, black) and amoxicillin (AMX, red) groups analyzed by t-tests (A,E) or non-parametric Mann–Whitney test (B,C) are indicated as *p ≤ 0.05, ***p ≤ 0.001.
Specific Bacterial Genera Are Linked to Immune Regulatory Responses
To investigate the interplay between the bacterial composition of the small intestine and the host immune response, pairwise correlations between relative abundance of genera and the host immune parameters shown in Figure 3 and Supplementary Figure S3 were performed (Figures 4A,B). This revealed negative correlations between three bacterial genera with reduced relative abundance in the amoxicillin group, namely Lactobacillus, Rothia, and Veillonella, and the fraction of both CD25+FoxP3+ Tregs in the small intestine LP (both CD4+ and CD4– Tregs) and PPs (only CD4– Tregs) as well as with the fraction of B cells of all PPs lymphocytes. Further, a consortium of Escherichia/Shigella, Klebsiella, Bifidobacterium, and four Firmicutes genera for which the relative abundance was higher in the amoxicillin group compared to controls were positively correlated with the same lymphocyte populations. These observations were supported by the additional finding of a significant correlation between LP CD4+ Tregs and Escherichia/Shigella in the control group only (p = 0.035) and for Turicibacter in the amoxicillin group only (p = 0.014). In addition, a negative correlation was found between the relative abundance of Propionibacterium and the expression of the Cldn2 gene encoding the tight junction protein Claudin 2 in the small intestine tissue, and a positive correlation between the relative abundance of Weissella and the fraction of Helios– induced Tregs in circulation.
Figure 4. Correlations between small intestine bacterial genera and host immune responses. (A) Spearman’s Rank correlation matrix between relative abundance of small intestinal bacterial genera present in minimum 10 animals, and all host response parameters indicated on Figure 3 and Supplementary Figure S3 including small intestinal (SI) gene expression and lymphocyte populations in lamina propria (LP), Peyer’s patches (PP) and blood. Only those parameters for which correlations are significant after applying FDR-correction are included. Spearman’s Rank correlation coefficients are indicated by the cell color. Statistically significant differences are indicated as *q < 0.1, **q < 0.05. (B) Visualization of selected correlations between relative abundance of small intestine genera and host responses in control (CTR, black) and amoxicillin (AMX, red) administered group. Each symbol represents one animal.
Discussion
To investigate the putative relationship between microbiota composition and immune regulation, BN rats were orally administered with the antibiotic amoxicillin or water as control for 7 days, and the gut microbiota, host immune regulation and intestinal permeability to cow’s milk protein were assed. Despite of epidemiologic studies showing an association between early life antibiotic consumption and prevalence of food allergies (Metsälä et al., 2013; Hirsch et al., 2017; Mitre et al., 2018) and autoimmune disorders such as IBD (Hviid et al., 2011; Shaw et al., 2011; Kronman et al., 2012), reflecting a failure to develop appropriate tolerance, the current study surprisingly indicated that a 7 days amoxicillin-induced perturbation of the gut microbiota promoted multiple parameters associated with immune regulation. These included higher relative abundance of small intestine Tregs, total fecal IgA and goblet cell numbers in the amoxicillin group compared to controls.
The observed increase in goblet cell numbers is in agreement with a previous study showing that administration the board-spectrum antibiotic ceftriaxone to Wistar rats causes increased goblet cell numbers and size (Holota et al., 2019). It is well established that bacterial colonization affects mucin gene expression, and in agreement with the present study, a negative correlation between abundance of Lactobacillus spp. and Muc2 gene expression in newborn mice has previously been reported (Bergström et al., 2012). Goblet cells promote oral tolerance by delivering luminal antigens to tolerogenic CD103+ dendritic cells in the small intestine to generate antigen specific Tregs (McDole et al., 2012; Shan et al., 2013) and by producing mucus, which contributes to intestinal barrier function and provides attachment sites and nutrients for commensal gut bacteria (Martens et al., 2009).
In the current study, the relative abundance of CD25+FoxP3+ Tregs in the small intestine LP (both CD4– and CD4+) and PPs (only CD4– Tregs) were found to be higher in the amoxicillin group than in controls. While conventional CD4+ Tregs have been extensively studied, the function of CD8+CD4– Tregs remains less characterized, but are thought to play an important role in preventing food allergy (Yamada et al., 2009). In line with the results from the current study, oral administration of broad-spectrum antibiotic cocktails is often reported to increase the relative abundance of CD4+ Treg cells in small intestine (Ivanov et al., 2008; Ochoa-Repáraz et al., 2009; Ichinohe et al., 2011; Janelsins et al., 2014; Benakis et al., 2016; Nakamura et al., 2016; Strzępa et al., 2017), and to reduce relative abundance of Th1 (Ochoa-Repáraz et al., 2009; Hill et al., 2010) and/or Th17 effector cells (Ivanov et al., 2008; Ochoa-Repáraz et al., 2009; Hill et al., 2010; Janelsins et al., 2014; Benakis et al., 2016). This shift in the balance between Tregs and effector T cells has been shown to reduce severity of multiple autoimmune (Ochoa-Repáraz et al., 2009; Benakis et al., 2016; Nakamura et al., 2016) and allergic diseases (Strzępa et al., 2017). On the other hand, the antibiotic-induced Treg promotion has been shown to prevent the development of proper immune responses to oral vaccination (Hall et al., 2008) and enteric infections (Sekirov and Finlay, 2009; Thackray et al., 2018).
Tregs are central stimulators of mucosal IgA production by plasma cells (Feng et al., 2011), which is reflected in the current study, since the higher relative abundance of small intestinal Tregs in the amoxicillin group than in controls, was accompanied by higher fecal IgA levels. Mucosal IgA is known to be important for the immunological defense against pathogenic bacteria and toxins, while the role of IgA in immune tolerance is less understood (Tordesillas and Berin, 2018). Still, results from human (Kukkonen et al., 2010) and animal (Frossard et al., 2004) studies indicate that IgA might play an important role for tolerance development, and alterations in IgA coating patterns of intestinal bacteria have been associated with multiple disorders including atopic diseases (Dzidic et al., 2017) and IBD (Van Der Waaij et al., 2004), although the results are ambiguous. Opposite to what was observed for fecal IgA, total serum IgA levels were significantly lower in the amoxicillin group than in controls. This is in line with previous studies showing that fecal and serum IgA levels are not affected in the same way by oral exposure to allergens (Frossard et al., 2004) and high-fat diet (Luck et al., 2019) in mice. Amoxicillin-induced perturbation of the microbiota may affect the transport of IgA from circulation to the intestinal lumen and/or tissue specific differences in regulation of IgA secretion.
To investigate the interplay between microbiota composition and immune regulation, the fecal microbiota was analyzed throughout the study, and the small intestinal and cecal microbiota at termination. Amoxicillin was found to dramatically, and differently alter the microbiota in the different intestinal compartments. The small intestine was vastly dominated by Lactobacillus in the control group thus this compartment had lower α diversity compared to feces and cecum. This genus was much less abundant in the amoxicillin group, where multiple different genera were more abundant with high variation between individual animals causing a higher α and β diversity compared to controls. This was in contrast to the fecal and cecal microbiota, for which both α and β diversity were significantly lower in the amoxicillin group compared to controls. These differences in response probably result from differences in bacterial load, diversity and composition in the different compartments (Hayashi et al., 2005; Wang et al., 2005; Seekatz et al., 2019). Differences in the transit time, which is much faster in the small intestine compared to cecum and colon (Wang et al., 2015), imply a shorter contact time between bacteria and antibiotics, which may also have affected the response. The current study highlights that changes in fecal bacterial population do not provide a good proxy for small intestinal responses to antibiotics.
Common for all analyzed intestinal compartments was the reduction in relative abundance of Lactobacillus and increase in relative abundance of Gammaproteobacteria. This was also seen in feces after 7 days treatment of healthy adults with amoxicillin in combination with β-lactamase inhibitor clavulanate (Kabbani et al., 2017; MacPherson et al., 2018). The elevated level of facultative anaerobes within the Gammaproteobacteria may in part be explained by depletion of butyrate producing clostridia, since reduction in butyrate has been reported to promote expansion of oxygen tolerating species due to increased epithelial oxygenation (Rivera-Chávez et al., 2016). In line with this explanation, previous studies report reduced butyrate levels and elevated cecal pH after amoxicillin administration in Wistar rats (Tulstrup et al., 2015, 2018).
Multiple studies have revealed that bacterial stimulation is important for mucosal T cell responses and in particular important for balancing the ratio between Tregs and effector T (Lee and Kim, 2017). In the current study, a consortium consisting of Escherichia/Shigella, Klebsiella (Gammaproteobacteria), Bifidobacterium, and four genera belonging to Firmicutes, for which the relative abundance in the small intestine were higher in the amoxicillin group than in controls, positively correlated with the fraction of small intestinal Tregs. This is in accordance with previous studies reporting positive effects of administration with various probiotic Bifidobacterium spp. on the Treg/T effector cell balance in humans (Konieczna et al., 2012) and mice, which thereby ameliorate allergic and autoimmune diseases (Lyons et al., 2010; Zuo et al., 2014; Liu et al., 2018; Kim et al., 2019). Furthermore, a growing body of evidence suggests an important immunomodulatory role of hexa-acylated lipopolysaccharide (LPS), a major constituent of the outer membrane of Gammaproteobacteria, which exhibit superior Toll-like receptor 4 (TLR4) stimulatory capacity compared to the less acylated LPS present from other Gram-negative bacteria (Larsen et al., 2012, 2015; Li et al., 2013; Brix et al., 2015). We speculate that bacterial derived signals such as hexa-acylated LPS from Gammaproteobacteria could directly promote the expansion of the intestinal Tregs, which in turn stimulate mucosal IgA production by local plasma cells (Figure 5). Indeed, in vitro studies have shown that LPS can activate murine (Caramalho et al., 2011) and human (Lewkowicz et al., 2016) Tregs directly via TLR4 activation. Furthermore, animal studies have indicated a link between high dietary LPS levels and higher numbers of CD4+ Tregs in mLN, as well as higher CD8+ Tregs in mLN and PPs (Hrncir et al., 2008). However, further studies are need to address this hypothesis, such as applying the same experimental treatment to LPS-non-responsive animals, to see if Treg and IgA induction fail to appear. One LPS-non-responsive animal model is TLR4–/–B6129PF mice, which are in fact prone to develop food allergy (Bashir et al., 2004).
Figure 5. Microbiota host interactions in the small intestine. The relative abundance of Lactobacillus was dramatically lower in the small intestine of amoxicillin administered rats, while multiple different genera were more abundant. Among these were a consortium including Escherichia/Shigella, Klebsiella (Gammaproteobacteria), and Bifidobacterium, which strongly correlated with the relative abundance of regulatory T cells (Tregs) in the small intestine. We speculate that bacterial derived signals such as hexa-acylated lipopolysaccharide (LPS) from Gammaproteobacteria, could directly promote the expansion of the intestinal Tregs, which in turn stimulate mucosal IgA production by local plasma cells and possibly goblet cell activation and mucus secretion. The figure is created with BioRender.com.
Conclusion
This study showed that microbiota perturbation by amoxicillin promotes acute intestinal immune regulation. The observed immune regulation may represent an acute response to overt inflammatory signals derived from the microbiota, e.g., released from dead bacteria and/or new or expanded bacterial species. The long-term immunological and clinical consequences of this effect remain to be investigated to determine whether this condition affects development of oral tolerance to luminal antigens and thus prevents inflammatory and allergic diseases.
Data Availability Statement
The 16S rRNA gene sequence data are deposited in the NCBI Sequence Read Archive with the accession number PRJNA599292.
Ethics Statement
Ethical approval was given by the Danish Animal Experiments Inspectorate and the authorization number given 2015-15-0201-00553-C1. The experiment was overseen by the National Food Institute’s in-house Animal Welfare Committee for animal care and use.
Author CONTRIBUTIONS
KB and TL conceived of the study. KG analyzed serum antibody levels, protein uptake, and histology. MB and KG analyzed microbiota data. JL analyzed the lymphocyte populations and participated in interpretation of the results. A-SB analyzed the gene expression. KG, MB, TL, and KB participated in the design of the study and the interpretation of the results. KG drafted the manuscript. All authors made substantial intellectual contributions, revised the manuscript critically, and approved the final version of the manuscript.
Funding
This work was supported by the Danish Dairy Research Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Elise Navntoft, Kenneth Rene Worm, Maja Danielsen, Olav Dahlgaard, and Anne Ørngreen for their great assistance with the animal experiment, and Bodil Madsen, Juliane Margrethe Gregersen, Sarah Grundt Simonsen, and Heidi Broksø Letting for their great assistance with the laboratory work. Furthermore, the authors wish to thank the students Louise Kristensen, Line Quist Stauersbøl, Sofie Emilie Hornslet, Helena Utzon, and Louise Almer for their excellent assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00496/full#supplementary-material
References
Abril-Gil, M., Garcia-Just, A., Pérez-Cano, F. J., Franch, À., and Castell, M. (2015). Development and characterization of an effective food allergy model in Brown Norway rats. PLoS One 10:e0125314. doi: 10.1371/journal.pone.0125314
Anderson, M. J., Ellingsen, K. E., and McArdle, B. H. (2006). Multivariate dispersion as a measure of beta diversity. Ecol. Lett. 9, 683–693. doi: 10.1111/j.1461-0248.2006.00926.x
Atarashi, K., Tanoue, T., Oshima, K., Suda, W., Nagano, Y., Nishikawa, H., et al. (2013). Treg induction by a rationally selected mixture of Clostridia strains from the human microbiota. Nature 500, 232–236. doi: 10.1038/nature12331
Atarashi, K., Tanoue, T., Shima, T., Imaoka, A., Kuwahara, T., Momose, Y., et al. (2011). Induction of colonic regulatory T cells by indigenous Clostridium species. Science 331, 337–341. doi: 10.1126/science.1198469
Bashir, M. E. H., Louie, S., Shi, H. N., and Nagler-Anderson, C. (2004). Toll-like receptor 4 signaling by intestinal microbes influences susceptibility to food allergy. J. Immunol. 172, 6978–6987. doi: 10.4049/jimmunol.172.11.6978
Benakis, C., Brea, D., Caballero, S., Faraco, G., Moore, J., Murphy, M., et al. (2016). Commensal microbiota affects ischemic stroke outcome by regulating intestinal γδ T cells. Nat. Med. 22, 516–523. doi: 10.1038/nm.4068
Benjamini, Y., Krieger, A. M., and Yekutieli, D. (2006). Adaptive linear step-up procedures that control the false discovery rate. Biometrika 93, 491–507. doi: 10.1093/biomet/93.3.491
Bergström, A., Kristensen, M. B., Bahl, M. I., Metzdorff, S. B., Fink, L. N., Frøkiær, H., et al. (2012). Nature of bacterial colonization influences transcription of mucin genes in mice during the first week of life. BMC Res. Notes 5:402. doi: 10.1186/1756-0500-5-402
Bodey, G. P., and Nance, J. (1972). Amoxicillin: in vitro and pharmacological studies. Antimicrob. Agents Chemother. 1, 358–362. doi: 10.1128/AAC.1.4.358
Bøgh, K. L., Kroghsbo, S., Dahl, L., Rigby, N. M., Barkholt, V., Mills, E. N. C., et al. (2009). Digested Ara h 1 has sensitizing capacity in Brown Norway rats. Clin. Exp. Allergy 39, 1611–1621. doi: 10.1111/j.1365-2222.2009.03333.x
Brix, S., Eriksen, C., Larsen, J. M., and Bisgaard, H. (2015). Metagenomic heterogeneity explains dual immune effects of endotoxins. J. Allergy Clin. Immunol. 135, 277–280. doi: 10.1016/j.jaci.2014.09.036
Cahenzli, J., Köller, Y., Wyss, M., Geuking, M. B., and McCoy, K. D. (2013). Intestinal microbial diversity during early-life colonization shapes long-term IgE levels. Cell Host Microbe 14, 559–570. doi: 10.1016/j.chom.2013.10.004
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
Caramalho, I., Rodrigues-Duarte, L., Perez, A., Zelenay, S., Penha-Gonçalves, C., and Demengeot, J. (2011). Regulatory T cells contribute to diabetes protection in lipopolysaccharide-treated non-obese diabetic mice. Scand. J. Immunol. 74, 585–595. doi: 10.1111/j.1365-3083.2011.02627.x
Christensen, E. G., Licht, T. R., Leser, T. D., and Bahl, M. I. (2014). Dietary Xylo-oligosaccharide stimulates intestinal bifidobacteria and lactobacilli but has limited effect on intestinal integrity in rats. BMC Res. Notes 7:660. doi: 10.1186/1756-0500-7-660
Clarke, K. R. (1993). Non−parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 18, 117–143. doi: 10.1111/j.1442-9993.1993.tb00438.x
Dolk, F. C. K., Pouwels, K. B., Smith, D. R. M., Robotham, J. V., and Smieszek, T. (2018). Antibiotics in primary care in England: which antibiotics are prescribed and for which conditions? J. Antimicrob. Chemother. 73(Suppl. 2), ii2–ii10. doi: 10.1093/jac/dkx504
Donaldson, G. P., Ladinsky, M. S., Yu, K. B., Sanders, J. G., Yoo, B. B., Chou, W. C., et al. (2018). Gut microbiota utilize immunoglobulin a for mucosal colonization. Science 360, 795–800. doi: 10.1126/science.aaq0926
Dzidic, M., Abrahamsson, T. R., Artacho, A., Björkstén, B., Collado, M. C., Mira, A., et al. (2017). Aberrant IgA responses to the gut microbiota during infancy precede asthma and allergy development. J. Allergy Clin. Immunol. 139, 1017–1025.e14. doi: 10.1016/j.jaci.2016.06.047
European Centre for Disease Prevention and Control, (2018). Antimicrobial Consumption Annual Epidemiological Report for 2017. Stockholm: European Centre for Disease Prevention and Control.
Feng, T., Elson, C. O., and Cong, Y. (2011). Treg cell-IgA axis in maintenance of host immune homeostasis with microbiota. Int. Immunopharmacol. 11, 589–592. doi: 10.1016/j.intimp.2010.11.016
Frossard, C. P., Hauser, C., and Eigenmann, P. A. (2004). Antigen-specific secretory IgA antibodies in the gut are decreased in a mouse model of food allergy. J. Allergy Clin. Immunol. 114, 377–382. doi: 10.1016/j.jaci.2004.03.040
Furusawa, Y., Obata, Y., Fukuda, S., Endo, T. A., Nakato, G., Takahashi, D., et al. (2013). Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature 504, 446–450. doi: 10.1038/nature12721
Hall, J. A., Bouladoux, N., Sun, C. M., Wohlfert, E. A., Blank, R. B., Zhu, Q., et al. (2008). Commensal DNA limits regulatory T cell conversion and is a natural adjuvant of intestinal immune responses. Immunity 29, 637–649. doi: 10.1016/j.immuni.2008.08.009
Hayashi, H., Takahashi, R., Nishi, T., Sakamoto, M., and Benno, Y. (2005). Molecular analysis of jejunal, ileal, caecal and recto-sigmoidal human colonic microbiota using 16S rRNA gene libraries and terminal restriction fragment length polymorphism. J. Med. Microbiol. 54(Pt 11), 1093–1101. doi: 10.1099/jmm.0.45935-0
Hill, D. A., Hoffmann, C., Abt, M. C., Du, Y., Kobuley, D., Kirn, T. J., et al. (2010). Metagenomic analyses reveal antibiotic-induced temporal and spatial changes in intestinal microbiota with associated alterations in immune cell homeostasis. Mucosal Immunol. 3, 148–158. doi: 10.1038/mi.2009.132
Hirsch, A. G., Pollak, J., Glass, T. A., Poulsen, M. N., Bailey-Davis, L., Mowery, J., et al. (2017). Early-life antibiotic use and subsequent diagnosis of food allergy and allergic diseases. Clin. Exp. Allergy. 47, 236–244. doi: 10.1111/cea.12807
Holota, Y., Dovbynchuk, T., Kaji, I., Vareniuk, I., Dzyubenko, N., Chervinska, T., et al. (2019). The long-term consequences of antibiotic therapy: role of colonic short-chain fatty acids (SCFA) system and intestinal barrier integrity. PLoS One 14:e0220642. doi: 10.1371/journal.pone.0220642
Hrncir, T., Stepankova, R., Kozakova, H., Hudcovic, T., and Tlaskalova-Hogenova, H. (2008). Gut microbiota and lipopolysaccharide content of the diet influence development of regulatory T cells: studies in germ-free mice. BMC Immunol. 9:65. doi: 10.1186/1471-2172-9-65
Hviid, A., Svanström, H., and Frisch, M. (2011). Antibiotic use and inflammatory bowel diseases in childhood. Gut 60, 49–54. doi: 10.1136/gut.2010.219683
Ichinohe, T., Pang, I. K., Kumamoto, Y., Peaper, D. R., Ho, J. H., Murray, T. S., et al. (2011). Microbiota regulates immune defense against respiratory tract influenza a virus infection. Proc. Natl. Acad. Sci. U.S.A. 108, 5354–5359. doi: 10.1073/pnas.1019378108
Ishikawa, H., Tanaka, K., Maeda, Y., Aiba, Y., Hata, A., Tsuji, N. M., et al. (2008). Effect of intestinal microbiota on the induction of regulatory CD25 + CD4+ T cells. Clin. Exp. Immunol. 153, 127–135. doi: 10.1111/j.1365-2249.2008.03668.x
Ivanov, I. I., Frutos, R., de, L., Manel, N., Yoshinaga, K., Rifkin, D. B., et al. (2008). Specific microbiota direct the differentiation of IL-17-producing T-helper cells in the mucosa of the small intestine. Cell Host Microbe 4, 337–349. doi: 10.1016/j.chom.2008.09.009
Janelsins, B. M., Lu, M., and Datta, S. K. (2014). Altered inactivation of commensal LPS due to acyloxyacyl hydrolase deficiency in colonic dendritic cells impairs mucosal Th17 immunity. Proc. Natl. Acad. Sci. U.S.A. 111, 373–378. doi: 10.1073/pnas.1311987111
Jensen, L. H., Larsen, J. M., Madsen, C. B., Laursen, R. R., Jacobsen, L. N., and Bøgh, K. L. (2019). Preclinical Brown Norway rat models for the assessment of infant formulas in the prevention and treatment of cow’s milk allergy. Int. Arch. Allergy Immunol. 178, 307–314. doi: 10.1159/000495801
Jung, T. H., Park, J. H., Jeon, W. M., and Han, K. S. (2015). Butyrate modulates bacterial adherence on LS174T human colorectal cells by stimulating mucin secretion and MAPK signaling pathway. Nutr. Res. Pract. 9, 343–349. doi: 10.4162/nrp.2015.9.4.343
Kabbani, T. A., Pallav, K., Dowd, S. E., Villafuerte-Galvez, J., Vanga, R. R., Castillo, N. E., et al. (2017). Prospective randomized controlled study on the effects of Saccharomyces boulardii CNCM I-745 and amoxicillin-clavulanate or the combination on the gut microbiota of healthy volunteers. Gut Microbes 8, 17–32. doi: 10.1080/19490976.2016.1267890
Kim, W. G., Kang, G. D., Kim, H. I, Han, M. J., and Kim, D. H. (2019). Bifidobacterium longum IM55 and Lactobacillus plantarum IM76 alleviate allergic rhinitis in mice by restoring Th2/Treg imbalance and gut microbiota disturbance. Benef. Microbes 10, 55–67. doi: 10.3920/BM2017.0146
Knippels, L. M. J., Penninks, A. H., Van Meeteren, M., and Houben, G. F. (1999). Humoral and cellular immune responses in different rat strains on oral exposure to ovalbumin. Food Chem. Toxicol. 37, 881–888. doi: 10.1016/S0278-6915(99)00064-2
Konieczna, P., Groeger, D., Ziegler, M., Frei, R., Ferstl, R., Shanahan, F., et al. (2012). Bifidobacterium infantis 35624 administration induces Foxp3 T regulatory cells in human peripheral blood: potential role for myeloid and plasmacytoid dendritic cells. Gut 61, 354–366. doi: 10.1136/gutjnl-2011-300936
Kronman, M. P., Zaoutis, T. E., Haynes, K., Feng, R., and Coffin, S. E. (2012). Antibiotic exposure and IBD development among children: a population-based cohort study. Pediatrics 130, e794–e803. doi: 10.1542/peds.2011-3886
Kukkonen, K., Kuitunen, M., Haahtela, T., Korpela, R., Poussa, T., and Savilahti, E. (2010). High intestinal IgA associates with reduced risk of IgE-associated allergic diseases. Pediatr. Allergy Immunol. 21(1 Pt 1), 67–73. doi: 10.1111/j.1399-3038.2009.00907.x
Larsen, J. M., Musavian, H. S., Butt, T. M., Ingvorsen, C., Thysen, A. H., and Brix, S. (2015). Chronic obstructive pulmonary disease and asthma-associated Proteobacteria, but not commensal Prevotella spp., promote Toll-like receptor 2-independent lung inflammation and pathology. Immunology 144, 333–342. doi: 10.1111/imm.12376
Larsen, J. M., Steen-Jensen, D. B., Laursen, J. M., Søndergaard, J. N., Musavian, H. S., Butt, T. M., et al. (2012). Divergent pro-inflammatory profile of human dendritic cells in response to commensal and pathogenic bacteria associated with the airway microbiota. PLoS One 7:e31976. doi: 10.1371/journal.pone.0031976
Lee, N., and Kim, W. U. (2017). Microbiota in T-cell homeostasis and inflammatory diseases. Exp. Mol. Med. 49:e340. doi: 10.1038/emm.2017.36
Lewkowicz, N., Mycko, M. P., Przygodzka, P., Ćwiklińska, H., Cichalewska, M., Matysiak, M., et al. (2016). Induction of human IL-10-producing neutrophils by LPS-stimulated Treg cells and IL-10. Mucosal Immunol. 9, 364–378. doi: 10.1038/mi.2015.66
Li, Y., Wang, Z., Chen, J., Ernst, R. K., and Wang, X. (2013). Influence of lipid a acylation pattern on membrane permeability and innate immune stimulation. Mar. Drugs 11, 3197–3208. doi: 10.3390/md11093197
Liu, Q., Jing, W., and Wang, W. (2018). Bifidobacterium lactis ameliorates the risk of food allergy in Chinese children by affecting relative percentage of Treg and Th17 cells. Can. J. Infect. Dis. Med. Microbiol. 2018:4561038. doi: 10.1155/2018/4561038
Luck, H., Khan, S., Kim, J. H., Copeland, J. K., Revelo, X. S., Tsai, S., et al. (2019). Gut-associated IgA+ immune cells regulate obesity-related insulin resistance. Nat. Commun. 10:3650. doi: 10.1038/s41467-019-11370-y
Lyons, A., O’Mahony, D., O’Brien, F., MacSharry, J., Sheil, B., Ceddia, M., et al. (2010). Bacterial strain-specific induction of Foxp3+ T regulatory cells is protective in murine allergy models. Clin. Exp. Allergy 40, 811–819. doi: 10.1111/j.1365-2222.2009.03437.x
MacPherson, C. W., Mathieu, O., Tremblay, J., Champagne, J., Nantel, A., Girard, S. A., et al. (2018). Gut bacterial microbiota and its resistome rapidly recover to basal state levels after short-term amoxicillin-clavulanic acid treatment in healthy adults. Sci. Rep. 8:11192. doi: 10.1038/s41598-018-29229-5
Mandal, S., Van Treuren, W., White, R. A., Eggesbø, M., Knight, R., and Peddada, S. D. (2015). Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb. Ecol. Health Dis. 26:27663. doi: 10.3402/mehd.v26.27663
Martens, E. C., Roth, R., Heuser, J. E., and Gordon, J. I. (2009). Coordinate regulation of glycan degradation and polysaccharide capsule biosynthesis by a prominent human gut symbiont. J. Biol. Chem. 284, 18445–18457. doi: 10.1074/jbc.M109.008094
McDole, J. R., Wheeler, L. W., McDonald, K. G., Wang, B., Konjufca, V., Knoop, K. A., et al. (2012). Goblet cells deliver luminal antigen to CD103 + dendritic cells in the small intestine. Nature 483, 345–349. doi: 10.1038/nature10863
Metsälä, J., Lundqvist, A., Virta, L. J., Kaila, M., Gissler, M., and Virtanen, S. M. (2013). Mother’s and offspring’s use of antibiotics and infant allergy to cow’s milk. Epidemiology 24, 303–309. doi: 10.1097/EDE.0b013e31827f520f
Mitre, E., Susi, A., Kropp, L. E., Schwartz, D. J., Gorman, G. H., and Nylund, C. M. (2018). Association between use of acid-suppressive medications and antibiotics during infancy and allergic diseases in early childhood. JAMA Pediatr. 172:e180315. doi: 10.1001/jamapediatrics.2018.0315
Moreau, M. C., Ducluzeau, R., Guy-Grand, D., and Muller, M. C. (1978). Increase in the population of duodenal immunoglobulin A plasmocytes in axenic mice associated with different living or dead bacterial strains of intestinal origin. Infect. Immun. 21, 532–539. doi: 10.1128/iai.21.2.532-539.1978
Nakamura, Y. K., Metea, C., Karstens, L., Asquith, M., Gruner, H., Moscibrocki, C., et al. (2016). Gut microbial alterations associated with protection from autoimmune uveitis. Invest. Ophthalmol. Vis. Sci. 57, 3747–3758. doi: 10.1167/iovs.16-19733
Ochoa-Repáraz, J., Mielcarz, D. W., Ditrio, L. E., Burroughs, A. R., Foureau, D. M., Haque-Begum, S., et al. (2009). Role of Gut commensal microflora in the development of experimental autoimmune encephalomyelitis. J. Immunol. 183, 6041–6050. doi: 10.4049/jimmunol.0900747
Peng, L., Li, Z.-R., Green, R. S., Holzman, I. R., and Lin, J. (2009). Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. J. Nutr. 139, 1619–1625. doi: 10.3945/jn.109.104638
Rivera-Chávez, F., Zhang, L. F., Faber, F., Lopez, C. A., Byndloss, M. X., Olsan, E. E., et al. (2016). Depletion of butyrate-producing clostridia from the Gut microbiota drives an aerobic luminal expansion of Salmonella. Cell Host Microbe 19, 443–454. doi: 10.1016/j.chom.2016.03.004
Rosshart, S. P., Herz, J., Vassallo, B. G., Hunter, A., Wall, M. K., Badger, J. H., et al. (2019). Laboratory mice born to wild mice have natural microbiota and model human immune responses. Science 365:eaaw4361. doi: 10.1126/science.aaw4361
Seekatz, A. M., Schnizlein, M. K., Koenigsknecht, M. J., Baker, J. R., Hasler, W. L., Bleske, B. E., et al. (2019). Spatial and temporal analysis of the stomach and small-intestinal microbiota in fasted healthy humans. mSphere 4:e00126-19. doi: 10.1128/msphere.00126-19
Sekirov, I., and Finlay, B. B. (2009). The role of the intestinal microbiota in enteric infection. J. Physiol. 587(Pt 17), 4159–4167. doi: 10.1113/jphysiol.2009.172742
Shan, M., Gentile, M., Yeiser, J. R., Walland, A. C., Bornstein, V. U., Chen, K., et al. (2013). Mucus enhances gut homeostasis and oral tolerance by delivering immunoregulatory signals. Science 342, 447–453. doi: 10.1126/science.1237910
Shaw, S. Y., Blanchard, J. F., and Bernstein, C. N. (2011). Association between the use of antibiotics and new diagnoses of Crohn’s disease and ulcerative colitis. Am. J. Gastroenterol. 106, 2133–2142. doi: 10.1038/ajg.2011.304
Shroff, K. E., Meslin, K., and Cebra, J. J. (1995). Commensal enteric bacteria engender a self-limiting humoral mucosal immune response while permanently colonizing the gut. Infect. Immun. 63, 3904–3913. doi: 10.1128/iai.63.10.3904-3913.1995
Stephen-Victor, E., and Chatila, T. A. (2019). Regulation of oral immune tolerance by the microbiome in food allergy. Curr. Opin. Immunol. 60, 141–147. doi: 10.1016/j.coi.2019.06.001
Strzępa, A., Majewska-Szczepanik, M., Lobo, F. M., Wen, L., and Szczepanik, M. (2017). Broad spectrum antibiotic enrofloxacin modulates contact sensitivity through gut microbiota in a murine model. J. Allergy Clin. Immunol. 140, 121–133.e33. doi: 10.1016/j.jaci.2016.11.052
Thackray, L. B., Handley, S. A., Gorman, M. J., Poddar, S., Bagadia, P., Briseño, C. G., et al. (2018). Oral antibiotic treatment of mice exacerbates the disease severity of multiple Flavivirus infections. Cell Rep. 22, 3440–3453.e6. doi: 10.1016/j.celrep.2018.03.001
Tordesillas, L., and Berin, M. C. (2018). Mechanisms of oral tolerance. Clin. Rev. Allergy Immunol. 62, 1523–1529. doi: 10.1007/s12016-018-8680-5
Tulstrup, M. V. L., Christensen, E. G., Carvalho, V., Linninge, C., Ahrné, S., Højberg, O., et al. (2015). Antibiotic treatment affects intestinal permeability and gut microbial composition in wistar rats dependent on antibiotic class. PLoS One 10:e0144854. doi: 10.1371/journal.pone.0144854
Tulstrup, M. V.-L., Roager, H. M., Thaarup, I. C., Frandsen, H. L., Frøkiær, H., Licht, T. R., et al. (2018). Antibiotic treatment of rat dams affects bacterial colonization and causes decreased weight gain in pups. Commun. Biol. 1:145. doi: 10.1038/s42003-018-0140-5
Van Der Waaij, L. A., Kroese, F. G. M., Visser, A., Nelis, G. F., Westerveld, B. D., Jansen, P. L. M., et al. (2004). Immunoglobulin coating of faecal bacteria in inflammatory. Eur. J. Gastroenterol. Hepatol. 16, 669–674. doi: 10.1097/01.meg.0000108346.41221.19
Walter, J., Tannock, G. W., Tilsala-Timisjarvi, A., Rodtong, S., Loach, D. M., Munro, K., et al. (2000). Detection and identification of gastrointestinal Lactobacillus species by using denaturing gradient gel electrophoresis and species-specific PCR primers. Appl. Environ. Microbiol. 66, 297–303. doi: 10.1128/AEM.66.1.297-303.2000
Wang, M., Ahrné, S., Jeppsson, B., and Molin, G. (2005). Comparison of bacterial diversity along the human intestinal tract by direct cloning and sequencing of 16S rRNA genes. FEMS Microbiol. Ecol. 54, 219–231. doi: 10.1016/j.femsec.2005.03.012
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, Y. T., Mohammed, S. D., Farmer, A. D., Wang, D., Zarate, N., Hobson, A. R., et al. (2015). Regional gastrointestinal transit and pH studied in 215 healthy volunteers using the wireless motility capsule: influence of age, gender, study country and testing protocol. Aliment. Pharmacol. Ther. 42, 761–772. doi: 10.1111/apt.13329
Weltzin, R., Lucia-Jandris, P., Michetti, P., Fields, B. N., Kraehenbuhl, J. P., and Neutra, M. R. (1989). Binding and transepithelial transport of immunoglobulins by intestinal M cells: demonstration using monoclonal IgA antibodies against enteric viral proteins. J. Cell Biol. 108, 1673–1685. doi: 10.1083/jcb.108.5.1673
Yamada, A., Ohshima, Y., Yasutomi, M., Ogura, K., Tokuriki, S., Naiki, H., et al. (2009). Antigen-primed splenic CD8+ T cells impede the development of oral antigen-induced allergic diarrhea. J. Allergy Clin. Immunol. 123, 889–894. doi: 10.1016/j.jaci.2008.12.1115
Keywords: antibiotics, microbiome, host response, Treg, regulatory T cells, IgA, goblet cells, mucus
Citation: Graversen KB, Bahl MI, Larsen JM, Ballegaard A-SR, Licht TR and Bøgh KL (2020) Short-Term Amoxicillin-Induced Perturbation of the Gut Microbiota Promotes Acute Intestinal Immune Regulation in Brown Norway Rats. Front. Microbiol. 11:496. doi: 10.3389/fmicb.2020.00496
Received: 24 January 2020; Accepted: 06 March 2020;
Published: 26 March 2020.
Edited by:
Markus M. Heimesaat, Charité – Universitätsmedizin Berlin, GermanyReviewed by:
Bong-Soo Kim, Hallym University, South KoreaJohn MacSharry, University College Cork, Ireland
Copyright © 2020 Graversen, Bahl, Larsen, Ballegaard, Licht and Bøgh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katrine Lindholm Bøgh, kalb@food.dtu.dk