Impact Factor 4.225 | CiteScore 5.0
More on impact ›

ORIGINAL RESEARCH article

Front. Pharmacol., 14 January 2021 | https://doi.org/10.3389/fphar.2020.574607

Gut Microbiota Combined With Metabolomics Reveals the Repeated Dose Oral Toxicity of β-Cyclodextrin in Mice

www.frontiersin.orgShuangyu Lv1, www.frontiersin.orgXiaomei Zhang1, www.frontiersin.orgYu Feng1, www.frontiersin.orgQiying Jiang1, www.frontiersin.orgChenguang Niu2, www.frontiersin.orgYanjie Yang1* and www.frontiersin.orgXinchun Wang2*
  • 1Institute of Molecular Medicine, School of Basic Medical Sciences, Henan University, Kaifeng, China
  • 2Key Laboratory of Clinical Resources Translation, The First Affiliated Hospital of Henan University, Kaifeng, China

Βeta-cyclodextrin (β-CD) with a hydrophobic cavity enables the formation of inclusion complexes with organic molecules. The formation of host–guest complexes makes the application of β-CD popular in many fields, but their interaction with organisms is poorly understood. In the present study, the effect of β-CD on gut microbiota (16S rRNA gene sequencing), serum metabolites (gas chromatography–mass spectrometry platform), and their correlation (Pearson correlation analysis) was investigated after 14 days repeated oral exposure in mice. β-CD did not significantly affect the α-diversity indexes, including Richness, Chao1, Shannon and Simpson indexes, but disturbed the structure of the gut bacteria according to the result of principal component analysis (PCA). After taxonomic assignment, 1 in 27 phyla, 2 in 48 classes, 3 in 107 orders, 6 in 192 families, and 8 in 332 genera were significantly different between control and β-CD treated groups. The serum metabolites were significantly changed after β-CD treatment according to the result of unsupervized PCA and supervised partial least squares-discriminant analysis (PLS-DA). A total of 112 differential metabolites (89 downregulated and 23 upregulated) were identified based on the VIP >1 from orthogonal PLS-DA and p <0.05 from Student’s t-test. The metabolic pathways, including ABC transporters, pyrimidine metabolism, purine metabolism, glucagon signaling pathway, insulin signaling pathway, and glycolysis/gluconeogenesis, were enriched by KEGG pathway analysis. Our study provides a general observation of gut microbiota, serum metabolites and their correlation after exposure to β-CD in mice, which will be helpful for future research and application of β-CD.

Introduction

Βeta-cyclodextrin (β-CD) is a cyclic oligosaccharide with a truncated cone shape. β-CD has a hydrophobic cavity and hydrophilic outer edge (Wu J.-W. et al., 2020; Zhong et al., 2020). The hydrophobicity of β-CD cavities enables formation of inclusion complexes with organic molecules by noncovalent intermolecular interactions, such as electrostatic, hydrophobic and Van der Waals interactions (Bouhadiba et al., 2020). The formation of host–guest complexes is important for their application (Bouhadiba, et al., 2020), including pharmaceuticals (Mura 2020), attachment and degradation of contaminants (Wu, et al., 2020; Zhou et al., 2020), and nanosponge preparation (Pawar et al., 2019). β-CD can reduce the agglomeration of nanomaterials and increase their specific surface area (Zhou, et al., 2020). The application and inevitable release of β-CD increase the opportunity for exposure, which may be harmful to humans and other organisms. However, the interaction between β-CD and organisms is poorly understood.

After oral exposure, β-CD passes through the digestive tract, where it may interact with the gut microbiota in the luminal environment. The gut microbiota is composed of different types of bacteria, viruses and fungi, and comprises 10-fold more cells than the human body (Yang et al., 2020; Zhang L. et al., 2020). The composition of gut microbiota plays crucial roles in the maintenance of host health (Chen et al., 2019; Yue et al., 2019). Dysbiosis of the gut microbiome can affect energy metabolism, nutritional digestion and absorption, immune status, and occurrence of many diseases (e.g. obesity, diabetes, arteriosclerosis, enteritis and hepatitis) (Chen et al., 2019; Vojinovic et al., 2019; Yue et al., 2019; Luo et al., 2020). The composition and diversity of the gut microbial community is highly sensitive to exogenous stressors (Chen et al., 2019), but the influence of β-CD on the gut microbiota is little known.

Symbiotic bacteria can regulate host metabolism, and circulating metabolite levels often act as an intermediary between the gut microbiota and host biology (Luo et al., 2020; Wu W. et al., 2020). Therefore, we speculate that oral exposure to β-CD may affect the gut microbiota and then alter serum metabolism. The metabolome as the final downstream product of the genome, transcriptome, and proteome, provides global and system-wide metabolic changes (Lindeque et al., 2018; Yang et al., 2019). Gas chromatography–mass spectrometry (GC–MS) with higher sensitivity and resolution is a robust platform for analysis of small molecules, and has been widely applied in metabolomics studies (Duan et al., 2020). We selected the GC–MS platform for detecting serum metabolites in this study. The combination of metabolome and microbiome is a promising approach to evaluate the relationship between host metabolism and gut microbiota (Luo et al., 2020).

The aim of this study was to explore whether oral administration of β-CD affected the composition and diversity of the gut microbiota in mice, and then indirectly caused a series of changes in circulating metabolites. We integrated gut microbiota, serum metabolites, and their correlation to explore the possible effects of β-CD on host health.

Materials and Methods

Animals and Treatment

Male Kunming mice of 22 ± 2 g (age 6–8 weeks) were supplied by the Animal Center of Henan Province (Zhengzhou, China). The mice were housed (4 or 5/cage) in a controlled environment of 22 ± 1 °C, 50–60% relative humidity, and 12 h light–dark cycle, with free access to food and water. The study and protocol were approved by the Committee of Medical Ethics and Welfare for Experimental Animals of Henan University School of Medicine (HUSOM2020-257).

Mice were randomly divided into the control group (normal saline) and β-CD treated (#C8510; Solarbio Science and Technology Co. Ltd., Beijing, China) group, with 10 animals per group. Mice were gavaged consecutively for 14 days with a total volume of 100 μL (17.5 g/L β-CD). The dose was selected based on previous studies of Jung et al. research on microbial communities and the aqueous solubility of 1.85 g/100 ml (Yu et al., 2008; Jung et al., 2015). At 2 h after the last gavage, fresh feces were collected, snap-frozen in liquid nitrogen, and stored at −80 °C for subsequent analysis. The body weight was measured, and blood was collected before killing by cervical dislocation.

Coefficients of Organs

After the mice were sacrificed, heart, liver, spleen, lung, and kidney tissues were carefully removed and weighed. The coefficients of organs were calculated as the ratio of organs weight (mg) to body weight (g).

Gut Microbiota Analysis

Gut microbiota were analyzed by 16S rRNA gene sequencing. Details on DNA extraction, 16S rRNA gene sequencing, sequence processing and analysis are shown below.

DNA Extraction and 16S rRNA Gene Sequencing Genomic DNA was extracted using a QIAamp 96 PowerFecal QIAcube HT kit (Hilden, Germany). DNA quantity and quality were determined with a NanoDrop 2000 UV-Vis Spectrophotometer (Thermo Scientific, Wilmington, DE, United States) and agarose gel electrophoresis. DNA samples were diluted to 1 ng/μL and used as template for PCR amplification of bacterial 16S rRNA genes. Universal primers 343F (5′-TACGGRAGGCAGCAG-3′) and 798R (5′-AGG​GTA​TCT​AAT​CCT-3′) for V3–V4 variable regions were used. Amplicon quality was detected by gel electrophoresis. After purification with Agencourt AMPure XP beads, another round of PCR was done. The final purified amplicon was quantified by Life Technologies Qubit dsDNA assay kit (Carlsbad, CA, United States). Equal amounts of purified amplicon were pooled, and sequencing was carried out using an Illumina MiSeq platform (San Diego, CA, United States).

Sequence Processing and Analysis Raw sequences data were preprocessed using Trimmomatic software (Bolger et al., 2014) to cut off ambiguous bases (N) and low quality sequences with average quality scores below 20. FLASH software was used to assemble the paired-end reads following the parameters of 10 bp of minimal overlapping, 200 bp of maximum overlapping, and 20% maximum mismatch rate (Reyon et al., 2012). Further denoizing was performed using QIIME software (version 1.8.0), and we abandoned the sequences with ambiguous, homologous sequences, below 200 bp or chimeras (Caporaso et al., 2010). Removed primer sequences and clean reads were clustered to generate operational taxonomic units (OTUs) based on 97% similarity cutoff by Vsearch software (version 2.4.2) (Edgar et al., 2011). The representative read (maximum richness) of each OTU was selected using the QIIME package, and then blasted and annotated against Silva database (version 123) using RDP classifier to obtain taxonomic information (Wang et al., 2007). The confidence threshold was set to 70%. The community compositions of each sample were counted at the level of phylum, class, order, family, and genus. The α-diversity was calculated for each sample, including Richness, Chao1, Simpson, and Shannon indexes. Principal component analysis (PCA) was used to reveal the β diversity. Analysis of variance (ANOVA) was used to compare the species differences between the control and β-CD treated groups.

Metabolomics Research

Serum metabolites were analyzed by GC–MS, and detailed information of sample preparation, data preprocessing, and GC–MS and multivariate analyses are shown below.

Sample Preparation Serum was obtained by centrifugation at 3,000 rpm for 10 min, followed by further centrifugation (Eppendorf 5810R, Germany) 12,000 rpm for 10 min at 4 °C. The supernatant was collected and stored at −80 °C. Thawed serum (80 μL) and 10 μL 0.3 mg/ml methanol solution of 2-chloro-l-phenylalanine (internal standard) were mixed and vortexed for 10 s. We added 240 μL cold mixture of methanol and acetonitrile (2/1, v/v) and vortexed for 1 min, then ultrasonicated in ice water for 5 min, and stored at −20 °C for 10 min. Samples were centrifuged at 12,000 rpm for 10 min at 4 °C, and 150 μL supernatant was collected and dried under vacuum. Subsequently, we added 80 μL 15 mg/ml methoxyamine hydrochloride in pyridine, vortexed vigorously for 2 min, and incubated at 37 °C for 90 min. We added 80 μL N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA, with 1% trimethylchlorosilane), 20 μL n-hexane and another 11 fatty acid methyl ester (C8/C9/C10/C12/C14/C16, 0.8 mg/ml; C18/C20/C22/C24/C26, 0.4 mg/ml) internal standard in chloroform, and derivatized at 70 °C for 1 h before GC–MS analysis. Quality control (QC) samples were prepared by mixing aliquots of all the samples and analyzing with the same procedure.

GC–MS Analysis The derivatized samples were analyzed using an Agilent 7890B GC system coupled to an Agilent 5977A MSD system (Agilent Technologies Inc., Santa Clara, CA, United States). We used a DB-5MS fused-silica capillary column (30 m × 0.25 mm × 0.25 μm; Agilent J and W Scientific, Folsom, CA, United States). Helium (>99.999%) was used as the carrier gas at a constant flow rate of 1.0 ml/min through the column. The temperature of the injector was maintained at 260 °C, and injection volume was 1 μL by splitless mode. The initial oven temperature was 60 °C, increased to 125 °C at a rate of 8 °C/min, to 210 °C at 5 °C/min, to 270 °C at 10 °C/min, to 305 °C at 20 °C/min and finally held at 305 °C for 5 min. The temperatures of the MS quadrupole and electron impact ion source were set to 150 and 230 °C, respectively. The collision energy was 70 eV. MS data were acquired in a full-scan mode (m/z 50–500), The time of solvent delay was set to 5 min. QC samples were detected at regular intervals (every 10 samples) throughout the analysis.

Data Preprocessing and Multivariate Statistical Analysis Raw data (.D format) were converted to.abf format by Analysis base File Converter software. MS-DIAL software was used to analyze the.abf data. Metabolites were annotated through Untarget database of GC–MS from Luming Biotechnology Co. Ltd. (Shanghai, China). All peak areas were processed by normalization of multi-interior label (RSD <0.3) according to the retention time partition period. Data were imported into R ropls package, where unsupervized PCA, supervised partial least squares-discriminant analysis (PLS-DA) and orthogonal PLS-DA (OPLS-DA) were performed. The differential metabolites were selected based on variable influence on projection (VIP) values of OPLS-DA model larger than 1.0, as well as p values of two-tailed Student’s t-test <0.05. Enrichment analysis of differential metabolites was performed by KEGG (http://www.kegg.jp/), and p < 0.05 was considered significantly enriched.

Integration Analysis Between the Modified Gut Microbiota and Metabolites

All the differential phyla and genera according to ANOVA (p < 0.05), and significantly changed metabolites (VIP >1 and p < 0.05) were selected. The Pearson correlation coefficient between bacteria and metabolites was calculated; *p < 0.05 and **p < 0.01 according to Pearson correlation.

Results

Coefficients of Organs

No abnormal behavior including eating, drinking, or activity, was found in the β-CD treated group. The body weights of the β-CD treated group were comparable to those of the control group (22.79 ± 0.76 g to 31.99 ± 0.72 g) (Supplementary Figure S1). The coefficients of lung (control, 6.82 ± 0.21 mg/g), liver (control, 52.93 ± 1.67 mg/g), spleen (control, 2.63 ± 0.23 mg/g), heart (control, 5.42 ± 0.23 mg/g), and kidney (control, 13.78 ± 0.34 mg/g) are shown in Supplementary Table S1. None of them was affected by oral administration of β-CD for 14 days.

Effects of Oral Exposure to β-CD on Gut Microbiota

Fresh fecal samples were collected after oral administration of β-CD for 14 days 16S rRNA gene sequencing was performed to construct the microbial community profiles. There were 17,004–33,139 valid tags in the fecal samples (Supplementary Table S2). A total of 3,031 OTUs were generated based on 97% similarity, including 1,064–1,463 OTUs per sample. The α-diversity indexes, including Richness, Chao1, Shannon, and Simpson indexes, were calculated (Figures 1A–D). None of the above indexes was disturbed by repeated oral administration of β-CD. β-CD affected the structure of the gut microbiota according to the result of PCA analysis (Figure 1E).

FIGURE 1
www.frontiersin.org

FIGURE 1. Alpha- and beta-diversity of the gut microbiota in mice after oral exposure to β-CD for 14 days. (A) Richness, (B) Chao1, (C) Shannon, (D) Simpson, (E) Principal component analysis (PCA). All data are presented as mean ± SEM (n = 10). p value vs. control according to Kruskal–Wallis. No significant difference was found between control and β-CDtreated groups.

According to the taxonomic assignment, the top 10 relative abundances of bacteria at the levels of phylum, class, order, family, and genus are shown in Figures 2A–E. At the phylum level, ∼90% of the sequences were within the top five phyla of Bacteroidetes, Firmicutes, Proteobacteria, Actinobacteria and Acidobacteria. The abundance of Proteobacteria was significantly increased in the β-CD treated group (p < 0.05, Figure 3A). The classification of bacteria tended to be dispersed in other levels, and 2 in 48 classes (uncultured Candidatus Saccharibacteria bacterium, JG37 AG 4), 3 in 107 orders (uncultured Candidatus Saccharibacteria bacterium, Subgroup 13, NKB5), 6 in 192 families (Clostridiales vadinBB60 group, uncultured Candidatus Saccharibacteria bacterium, Sandaracinaceae, Oligoflexaceae, Cystobacteraceae, Acidimicrobiaceae), and 8 in 332 genera (Inquilinus, Erysipelatoclostridium, Lactococcus, uncultured Candidatus Saccharibacteria bacterium, Paenibacillus, Ohtaekwangia, Ferruginibacter, Niastella) differed significantly between the control and β-CD treated group (Figures 3B–E). Taxonomic distributions and significant differences of bacteria at the species level were shown in Supplementary Figure S2.

FIGURE 2
www.frontiersin.org

FIGURE 2. Taxonomic distributions of bacteria (Top 10) at the (A) Phylum, (B) Class, (C) Order, (D) Family, and (E) Genus.

FIGURE 3
www.frontiersin.org

FIGURE 3. Scatter plots of significantly altered bacterial taxa, including (A) Phylum, (B) Class, (C) Order, (D) Family, (E) Genus by repeated oral exposure to β-CD for 14 days *p < 0.05 and **p < 0.01 according to ANOVA (n = 10).

Effects of Oral Exposure to β-CD on Serum Metabolites

Total ion current (TIC) chromatographs of control and β-CD treated samples are shown in Supplementary Figure S3. The reproducibility of methods was favorable according to the QC samples in PCA score plots (Supplementary Figure S4). The major spectrum can be assigned to specific metabolites by matching with the Untarget database of GC–MS from Lumingbio. Unsupervized PCA and supervised PLS-DA were selected to analyze the metabolic differences between the control and β-CD treated groups. There was an obvious separation tendency in PCA score plots (Figure 4A), and better separation was obtained in the PLS-DA model (Figure 4B), suggesting a significant difference after oral administration of β-CD for 14 days. Permutation test (200 times) was used to validate the reliability of the PLS-DA model (Figure 4D). OPLS-DA was performed and shown in Figure 4C.

FIGURE 4
www.frontiersin.org

FIGURE 4. Metabolic profiling between control and β-CD treated group. (A) score plots of the PCA model, (B) score plots of the PLS-DA model, (C) scores plot of the OPLS-DA model, (D) plot of the permutation test (200 times) of the PLS-DA model.

Combined with VIP > 1 from OPLS-DA and p < 0.05 from two-tailed Student’s t-test, 112 differential metabolites (89 downregulated and 23 upregulated) were identified (Table 1). To show the differences in expression of metabolites, cluster analysis of all the changed metabolites was conducted (Figure 5A). The abscissa represents the sample names, and the ordinate represents the differential metabolites. The color from green to red indicates expression abundance of metabolites from low to high. KEGG pathway analyses were performed on 112 differential metabolites to find out the metabolic pathway regulated by β-CD. The top 10 significantly enriched pathways (p < 0.05) were ATP-binding cassette (ABC) transporters, pyrimidine metabolism, glucagon signaling pathway, nonalcoholic fatty liver disease (NAFLD), glycolysis/gluconeogenesis, purine metabolism, central carbon metabolism in cancer, insulin signaling pathway, galactose metabolism, and FoxO signaling pathway (Figure 5B).

TABLE 1
www.frontiersin.org

Table 1. Significantly changed metabolites found in GC/MS-based metabolomic profiling.

FIGURE 5
www.frontiersin.org

FIGURE 5. Heat map (A) and scatterplot of significantly enriched KEGG pathways (B) of differential metabolites induced by repeated oral administration of β-CD for 14 days. In the heat map, the x axis represents the sample name and the y axis represents the differential metabolite. Green color indicates a low abundance of metabolites and red color indicates a high abundance. In the scatterplot, the size of the dots represents the metabolites number.

Association Analysis of Microbial Diversity With Serum Metabolites

Pearson correlation analysis was performed to detect correlations between gut microbiota and metabolite profiles modified by β-CD treatment. The phylum Proteobacteria was significantly negatively related to 3-pyridinol, nicotianamine, xanthosine and 2-deoxypentonic acid. At the genus level, Ohtaekwangia (mainly positive), Paenibacillus (mainly negative), uncultured Candidatus Saccharibacteria bacterium (mainly negative), and Erysipelatoclostridium (mainly positive), were closely related to many serum differential metabolites. However, the changes in Ferruginibacter were only significantly positively related to changes in putrescine (p < 0.05), and Lactococcus was significantly related to γ-aminobutyric acid (p < 0.05), O-phosphoethanolamine (p < 0.01), and bisphenol A (p < 0.01). The association analysis between the top 20 metabolites and microbial genera is shown in Figure 6, and the results at the levels of phylum, class, order, and family are shown in Supplementary Figure S5.

FIGURE 6
www.frontiersin.org

FIGURE 6. Pearson correlation analysis between modified metabolites and microbial genera by repeated oral administration of β-CD for 14 days *p < 0.05 and **p < 0.01 according to Pearson correlation.

Discussion

β-CD with a hydrophobicity cavity is popular in many fields, but the interaction between β-CD and organisms is little understood. In this study, the effects of oral exposure to β-CD on serum metabolites and gut microbiota were detected. A total of 1 phylum, 2 classes, 3 orders, 6 families and 8 genera were significantly different between the control and β-CDtreated groups. A total of 112 serum differential metabolites (89 downregulated and 23 upregulated) were detected after 14 days repeated oral administration of β-CD. Body weight as an important indicator for health was recorded, and no difference was found between the control and β-CD groups.

Proteobacteria contains several opportunistic pathogens and is suggested to exert proinflammatory activity both locally and systemically (Li et al., 2018; Ceccarani et al., 2020). The increased abundance of Proteobacteria induced by β-CD was significantly negatively related to serum differential metabolites of 3-pyridinol, nicotianamine, xanthosine and 2-deoxypentonic acid. Nicotianamine, a nonprotein amino acid, could inhibit the angiotensin-I-converting enzyme in the renin–angiotensin system and regulate blood pressure (Hayashi and Kimoto, 2007). Xanthosine is the intermediate in the metabolism of purine, and plays an important role in gene synthesis and metabolic regulation (Zhang et al., 2019). The effect of β-CD exposure on Clostridiales vadinBB60 group was highly correlated with urinary metabolites, and could also be changed by chronic cadmium exposure (He et al., 2020). The Cystobacteraceae family is well known for chitinase activity, and was upregulated by β-CD (Sharma and Subramanian, 2017). Acidimicrobiaceae is frequently encountered in acidic, metal-laden environments where they characteristically oxidize ferrous iron (Pinto et al., 2016). Fan et al. (2018) reported that Acidimicrobiaceae is capable of iron metabolism in the cyanobacteria-dominated open water zone. We speculated that the increased abundance of Acidimicrobiaceae induced by β-CD may affect iron metabolism in animals. Paenibacillus spp. upregulated by β-CD have been recognized as spoilage microorganisms and harmful to health, and they can survive conventional pasteurization regimens and grow during refrigerated storage (Nakano 2020). Paenibacillus was negatively related to many serum differential metabolites according to Pearson correlation analysis in our study. Erysipelatoclostridium is part of the normal gut microbiota but could become an opportunistic pathogen (Khan and Chousalkar, 2020). Lactococci are Gram-positive cocci and the Lactococcus genus includes five species, in which Lactococcus lactis is commonly used as a probiotic in food production and medical applications (Casalta and Montel, 2008; Rodrigues et al., 2016). β-CD-induced downregulation of Erysipelatoclostridium and upregulation of Lactococcus could be beneficial to the health of the host. Moreover, Lactococcus is significantly related to γ-aminobutyric acid, which is the principal inhibitory neurotransmitter and regulates neuronal excitability (Yi et al., 2017). Ferruginibacter upregulated by β-CD plays a key role in degrading organic components and denitrification (Liu et al., 2018). Ferruginibacter was only significantly positively related to the changed metabolites of putrescine, which is a naturally occurring polycation and plays an important role in cell growth, differentiation, and gene expression (Kwak et al., 2003). Uncultured Candidatus Saccharibacteria bacterium contributes to ammonium nitrogen (NH4+-N) removal (Ouyang et al., 2019), and was upregulated and negatively related to many changed serum metabolites in our study.

Pathways including ABC transporters, pyrimidine metabolism, purine metabolism, glucagon signaling, insulin signaling, glycolysis/gluconeogenesis and NAFLD were enriched by analysis of differential metabolites. ABC transporters mediate the ATP-dependent cellular export of a plethora of toxic metabolites and xenobiotic substances, and play key roles in multidrug resistance (Salustiano et al., 2020; Xiao et al., 2020). Dysregulation of the ABC transporters pathway indicates that animals might be under toxic stress after exposure to β-CD. It has also been reported that β-CD upregulates ABCA1 expression in vascular smooth muscle cells (Zhang C.-J. et al., 2020). Pyrimidines and purines are important components of DNA replication and RNA synthesis (Liu et al., 2019). Purines play a key role in neurotransmission and neuromodulation. Signaling molecules such as ATP and adenosine all belong to the purinergic system (Wang et al., 2019). Aberrations in purine metabolism by β-CD can affect neurological function and cell signal transduction. Abnormality of pyrimidine metabolism induced by β-CD may increase the risk of Alzheimer’s disease, immunodeficiency, growth retardation and aging (Wan et al., 2019). Glucagon promotes glycogenolysis and gluconeogenesis to elevate blood glucose levels in response to fasting, while insulin has the opposite effect (Jia et al., 2020; Liu et al., 2020). The glucagon signaling pathway, glycolysis/gluconeogenesis and insulin signaling pathway, as well as serum glucose level, were disturbed in this study. Therefore, we inferred that β-CD, a cyclic oligosaccharide consisting of seven glucose units (Alonso et al., 2018), could affect glucose metabolism after repeated oral administration. A previous study reported that intermediates of the glycolysis/gluconeogenesis pathway were related to insulin resistance, prediabetes, and diabetes (Guasch-Ferre et al., 2020). Diabetes and insulin resistance are both risk factors related to NAFLD (Komazaki et al., 2017). Path analysis showed that β-CD exposure also increased the risk of NAFLD in our study.

Conclusion

In this study, the effects of oral exposure to β-CD on serum metabolites and gut microbiota were investigated in mice. Body weight or coefficients of organs were not changed after 14 days repeated oral administration. β-CD did not significantly affect the α-diversity indexes, but disturbed the structure of the gut microbiota according to PCA. The regulated genera of Ohtaekwangia, Paenibacillus, uncultured Candidatus Saccharibacteria bacterium, and Erysipelatoclostridium were closely related to many serum differential metabolites according to Pearson correlation analysis. The metabolic pathway of ABC transporters, pyrimidine metabolism, purine metabolism, glucagon signaling pathway, and insulin signaling pathway were enriched by KEGG pathway analysis. Our study revealed the effects of β-CD on gut microbiota, serum metabolites and their correlation, which should be considered before their application.

Data Availability Statement

The 16s sequencing data has been deposited into BioProject (accession: PRJNA675411).

Ethics Statement

The animal study was reviewed and approved by The Committee of Medical Ethics and Welfare for Experimental Animals, Henan University School of Medicine, China.

Author Contributions

YY and XW developed the idea and designed the research. SL and XZ performed the experiments. SL, YF, QJ, and CN analyzed the data. SL wrote the draft of the manuscript. YY and XW contributed to revise the writing. All authors read and approved the submitted version.

Funding

This work was supported by research funds from the National Natural Science Foundation of China (Grant No. 81971280, and No. 81600974), the Key Science and Technology Program of Henan Province in China (Grant No. 192102310151, No. 202102310213, and No. 192102310080), the Program for Young Key Teacher of Henan Province (Grant No. 2020GGJS037), and the Youth Talent Promotion Plan of Henan Association for Science and Technology (Grant No. 2020HYTP054).

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.

Supplementary Material

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

References

Alonso, L., Fox, P. F., Calvo, M. V., and Fontecha, J. (2018). Effect of beta cyclodextrin on the reduction of cholesterol in Ewe’s milk manchego cheese. Molecules 23, 1789. doi:10.3390/molecules23071789

CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouhadiba, A., Rahali, S., Belhocine, Y., Allal, H., Nouar, L., and Rahim, M. (2020). Structural and energetic investigation on the host/guest inclusion process of benzyl isothiocyanate into β-cyclodextrin using dispersion-corrected DFT calculations. Carbohydr. Res. 491, 107980. doi:10.1016/j.carres.2020.107980

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi:10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Casalta, E., and Montel, M. C. (2008). Safety assessment of dairy microorganisms: the lactococcus genus. Int. J. Food Microbiol. 126, 271–273. doi:10.1016/j.ijfoodmicro.2007.08.013

Ceccarani, C., Bassanini, G., Montanari, C., Casiraghi, M. C., Ottaviano, E., Morace, G., et al. (2020). Proteobacteria overgrowth and butyrate-producing taxa depletion in the gut microbiota of glycogen storage disease type 1 patients. Metabolites 10, 133. doi:10.3390/metabo10040133

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Han, S., Zhou, D., Zhou, S., and Jia, G. (2019). Effects of oral exposure to titanium dioxide nanoparticles on gut microbiota and gut-associated metabolism in vivo. Nanoscale 11, 22398–22412. doi:10.1039/c9nr07580a

CrossRef Full Text | Google Scholar

Dong, X., Kong, L., Jing, Z., Wang, S., Lai, Y., Xie, M., et al. (2020). Facile synthesis of superparamagnetic beta-CD-MnFe2O4 as a peroxymonosulfate activator for efficient removal of 2,4- dichlorophenol: structure, performance, and mechanism. J. Hazard Mater. 394, 122528. doi:10.1016/j.jhazmat.2020.122528

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, L., Ma, A., Meng, X., Shen, G.-a., and Qi, X. (2020). QPMASS: a parallel peak alignment and quantification software for the analysis of large-scale gas chromatography-mass spectrometry (GC-MS)-based metabolomics datasets. J. Chromatogr., A 1620, 460999. doi:10.1016/j.chroma.2020.460999

CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi:10.1093/bioinformatics/btr381

CrossRef Full Text | Google Scholar

Fan, X., Ding, S., Gong, M., Chen, M., Gao, S., Jin, Z., et al. (2018). Different influences of bacterial communities on Fe (III) reduction and phosphorus availability in sediments of the cyanobacteria- and macrophyte-dominated zones. Front. Microbiol. 9, 2636. doi:10.3389/fmicb.2018.02636

PubMed Abstract | CrossRef Full Text | Google Scholar

Guasch-Ferré, M., Santos, J. L., Martínez-González, M. A., Clish, C. B., Razquin, C., Wang, D., et al. (2020). Glycolysis/gluconeogenesis- and tricarboxylic acid cycle-related metabolites, mediterranean diet, and type 2 diabetes. Am. J. Clin. Nutr. 111, 835–844. doi:10.1093/ajcn/ngaa016

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, P. (2020). Advantages of the combined use of cyclodextrins and nanocarriers in drug delivery: a review. Int. J. Pharm. 579, 119181. doi:10.1016/j.ijpharm.2020.119181

CrossRef Full Text | Google Scholar

Hayashi, A., and Kimoto, K. (2007). Nicotianamine preferentially inhibits angiotensin I‐converting enzyme. J. Nutr. Sci. Vitaminol. 53, 331–336. doi:10.3177/jnsv.53.331

He, X., Qi, Z., Hou, H., Gao, J., and Zhang, X. X. (2020). Effects of chronic cadmium exposure at food limitation-relevant levels on energy metabolism in mice. J. Hazard Mater. 388, 121791. doi:10.1016/j.jhazmat.2019.121791

CrossRef Full Text | Google Scholar

Izumi, M. K., Kensler, T. W., and Casero, R. A. (2003). Induction of phase 2 enzymes by serum oxidized polyamines through activation of Nrf2: effect of the polyamine metabolite acrolein. Biochem. Biophys. Res. Commun. 305, 662–670. doi:10.1016/s0006-291x(03)00834-9

CrossRef Full Text | Google Scholar

Ji, R., Zhao, S., Kong, N., Zhang, J., Loganathan, D., Mérette, S., et al. (2017). Quantitation of γ-aminobutyric acid in equine plasma by hydrophilic interaction liquid chromatography with tandem mass spectrometry. J. Sep. Sci. 40, 3239–3247. doi:10.1002/jssc.201700245

CrossRef Full Text | Google Scholar

Jia, L., Jiang, Y., Li, X., and Chen, Z. (2020). Purβ promotes hepatic glucose production by increasing Adcy6 transcription. Mol. Metab. 31, 85–97. doi:10.1016/j.molmet.2019.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, K., Kim, W., Park, G. W., Seo, C., Chang, H. N., and Kim, Y. C. (2015). Optimization of volatile fatty acids and hydrogen production from Saccharina japonica: acidogenesis and molecular analysis of the resulting microbial communities. Appl. Microbiol. Biotechnol. 99, 3327–3337. doi:10.1007/s00253-015-6419-2

CrossRef Full Text | Google Scholar

Khan, S., and Chousalkar, K. K. (2020). Salmonella typhimurium infection disrupts but continuous feeding of bacillus based probiotic restores gut microbiota in infected hens. J. Anim. Sci. Biotechnol. 11, 29. doi:10.1186/s40104-020-0433-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Komazaki, R., Katagiri, S., Takahashi, H., Maekawa, S., Shiba, T., Takeuchi, Y., et al. (2017). Periodontal pathogenic bacteria, aggregatibacter actinomycetemcomitans affect non-alcoholic fatty liver disease by altering gut microbiota and glucose metabolism. Sci. Rep. 7, 13950. doi:10.1038/s41598-017-14260-9

CrossRef Full Text | Google Scholar

Li, J., Lei, R., Li, X., Xiong, F., Zhang, Q., Zhou, Y., et al. (2018). The antihyperlipidemic effects of fullerenol nanoparticles via adjusting the gut microbiota in vivo. Part. Fibre Toxicol. 15, 5. doi:10.1186/s12989-018-0241-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Gu, C., Huang, W., Sheng, X., Du, J., and Li, Y. (2019). Targeted UPLC-MS/MS high-throughput metabolomics approach to assess the purine and pyrimidine metabolism. J. Chromatogr. B. 1113, 98–106. doi:10.1016/j.jchromb.2019.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Zhang, P., Li, H., Tian, Y., Wang, S., Song, Y., et al. (2018). Denitrification of landfill leachate under different hydraulic retention time in a two-stage anoxic/oxic combined membrane bioreactor process: performances and bacterial community. Bioresour. Technol. 250, 110–116. doi:10.1016/j.biortech.2017.11.026

CrossRef Full Text | Google Scholar

Morrissey, Y., Li, J., Sun, Y., Liang, Q., Peng, X., Liu, Y., et al. (2008). Solubility of beta-cyclodextrin in different mixed solvents. Petrol. Sci. 5, 263–268. doi:10.1007/s12182-008-0044-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakano, M. (2020). Development of a multiplex real-time PCR assay for the identification and quantification of group-specific Bacillus spp. and the genus Paenibacillus. Int. J. Food Microbiol. 323, 108573. doi:10.1016/j.ijfoodmicro.2020.108573

CrossRef Full Text | Google Scholar

Ouyang, E., Liu, Y., Ouyang, J., and Wang, X. (2019). Effects of different wastewater characteristics and treatment techniques on the bacterial community structure in three pharmaceutical wastewater treatment systems. Environ. Technol. 40, 329–341. doi:10.1080/09593330.2017.1393010

CrossRef Full Text | Google Scholar

Patterson, Y., Yang, J., Wang, J., Chen, L., Huang, H., Xiong, Y., et al. (2019). Quantification of serum purine metabolites for distinguishing patients with hepatitis B from hepatocellular carcinoma. Bioanalysis 11, 1003–1013. doi:10.4155/bio-2018-0319

PubMed Abstract | CrossRef Full Text | Google Scholar

Pawar, S., Shende, P., and Trotta, F. (2019). Diversity of β-cyclodextrin-based nanosponges for transformation of actives. Int. J. Pharm. 565, 333–350. doi:10.1016/j.ijpharm.2019.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, A. J., Sharp, J. O., Yoder, M. J., and Almstrand, R. (2016). Draft genome sequences of two novel acidimicrobiaceae members from an acid mine drainage biofilm metagenome. Genome Announc. 4, e01563. doi:10.1128/genomeA.01563-15

CrossRef Full Text | Google Scholar

Qin, L., Rimal, B., Nichols, R. G., Tian, Y., Smith, P. B., Hatzakis, E., et al. (2020). Perfluorooctane sulfonate alters gut microbiota-host metabolic homeostasis in mice. Toxicology 431, 152365. doi:10.1016/j.tox.2020.152365

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyon, D., Tsai, S. Q., Khayter, C., Foden, J. A., Sander, J. D., and Joung, J. K. (2012). FLASH assembly of TALENs for high-throughput genome editing. Nat. Biotechnol. 30, 460–465. doi:10.1038/nbt.2170

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, M. X., Lima, S. F., Higgins, C. H., Canniatti-Brazaca, S. G., and Bicalho, R. C. (2016). The lactococcus genus as a potential emerging mastitis pathogen group: a report on an outbreak investigation. J. Dairy Sci. 99, 9864–9874. doi:10.3168/jds.2016-11143

PubMed Abstract | CrossRef Full Text | Google Scholar

Salustiano, E. J., da Costa, K. M., Freire-de-Lima, L., Mendonça-Previato, L., and Previato, J. O. (2020). Inhibition of glycosphingolipid biosynthesis reverts multidrug resistance by differentially modulating ABC transporters in chronic myeloid leukemias. J. Biol. Chem. 295, 6457–6471. doi:10.1074/jbc.RA120.013090

CrossRef Full Text | Google Scholar

Sharma, G., and Subramanian, S. (2017). Unravelling the complete genome of archangium gephyra DSM 2261T and evolutionary insights into myxobacterial chitinases. Genome Biol. Evol. 9, 1304–1311. doi:10.1093/gbe/evx066

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, Y., Yang, L., Zhang, Y., Liu, X., Wu, Z., Gilbert, R. G., et al. (2020). Dendrobium officinale polysaccharide ameliorates diabetic hepatic glucose metabolism via glucagon-mediated signaling pathways and modifying liver-glycogen structure. J. Ethnopharmacol. 248, 112308. doi:10.1016/j.jep.2019.112308

CrossRef Full Text | Google Scholar

van Duijn, Q. L., Meng, X., Fu, X., Chen, B., Yang, J., Yang, H., et al. (2019). Intermediate metabolites of the pyrimidine metabolism pathway extend the lifespan of C. elegans through regulating reproductive signals. Aging (Albany NY) 11, 3993–4010. doi:10.18632/aging.102033

PubMed Abstract | CrossRef Full Text | Google Scholar

Vojinovic, D., Radjabzadeh, D., Kurilshikov, A., Amin, N., Wijmenga, C., Franke, L., et al. (2019). Relationship between gut microbiota and circulating metabolites in population-based cohorts. Nat. Commun. 10, 5813–5817. doi:10.1038/s41467-019-13721-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Chen, K., Li, J., Fang, Z., Pang, H., Yin, Y., et al. (2020). Gut microbiota combined with metabolomics reveals the metabolic profile of the normal aging process and the anti-aging effect of FuFang Zhenshu TiaoZhi(FTZ) in mice. Biomed. Pharmacother. 121, 109550. doi:10.1016/j.biopha.2019.109550

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi:10.1128/aem.00062-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J. W., Wu, C. R., Zhou, C. S., Dong, L. L., Liu, B. F., Xing, D. F., et al. (2020). Fate and removal of antibiotic resistance genes in heavy metals and dye co-contaminated wastewater treatment system amended with β-cyclodextrin functionalized biochar. Sci. Total Environ. 723, 137991. doi:10.1016/j.scitotenv.2020.137991

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W., Zhang, L., Xia, B., Tang, S., Liu, L., Xie, J., et al. (2020). Bioregional alterations in gut microbiome contribute to the plasma metabolomic changes in pigs fed with inulin. Microorganisms 8, 111. doi:10.3390/microorganisms8010111

Xiao, Q., Zhou, Y., and Lauschke, V. M. (2020). Ethnogeographic and inter-individual variability of human ABC transporters. Hum. Genet. 139, 623–646. doi:10.1007/s00439-020-02150-6

CrossRef Full Text | Google Scholar

Xie, Y., Li, W., Ran, L., Hou, R., Han, P., Lu, S., et al. (2020). Inclusion complexes of tea polyphenols with HP-β-cyclodextrin:Preparation, characterization, molecular docking, and antioxidant activity. J. Food Sci. 85, 1105–1113. doi:10.1111/1750-3841.15083

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, J. Z., Matthyser, A., Mason, S., Louw, R., and Taute, C. J. F. (2018). Metabolomics reveals the depletion of intracellular metabolites in HepG2 cells after treatment with gold nanoparticles. Nanotoxicology 12, 251–262. doi:10.1080/17435390.2018.1432779

CrossRef Full Text | Google Scholar

Yang, D., Zhao, D., Shah, S. Z. A., Wu, W., Lai, M., Zhang, X., et al. (2020). Implications of gut microbiota dysbiosis and metabolic changes in prion disease. Neurobiol. Dis. 135, 104704. doi:10.1016/j.nbd.2019.104704

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Lv, S., Wang, F., An, Y., Fang, N., Zhang, W., et al. (2019). Toxicity and serum metabolomics investigation of Mn-doped ZnS quantum dots in mice. Int. J. Nanomed. 14, 6297–6311. doi:10.2147/ijn.s212355

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, S., Zhao, D., Peng, C., Tan, C., Wang, Q., and Gong, J. (2019). Effects of theabrownin on serum metabolites and gut microbiome in rats with a high-sugar diet. Food Funct 10, 7063–7080. doi:10.1039/c9fo01334b

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C. J., Zhu, N., Liu, Z., Shi, Z., Long, J., Zu, X. Y., et al. (2020). Wnt5a/Ror2 pathway contributes to the regulation of cholesterol homeostasis and inflammatory response in atherosclerosis. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1865, 158547. doi:10.1016/j.bbalip.2019.158547

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, F., Liu, F., and Chen, W. (2019). Exposure to triclosan changes the expression of microRNA in male juvenile zebrafish (Danio rerio). Chemosphere 214, 651–658. doi:10.1016/j.chemosphere.2018.09.163

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: β-cyclodextrin, gut microbiota, metabolomics, oral exposure, mice

Citation: Lv S, Zhang X, Feng Y, Jiang Q, Niu C, Yang Y and Wang X (2021) Gut Microbiota Combined With Metabolomics Reveals the Repeated Dose Oral Toxicity of β-Cyclodextrin in Mice. Front. Pharmacol. 11:574607. doi: 10.3389/fphar.2020.574607

Received: 20 June 2020; Accepted: 27 November 2020;
Published: 14 January 2021.

Edited by:

Eleonore Fröhlich, Medical University of Graz, Austria

Reviewed by:

Yan Wang, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Sekena Hassanien Abdel-Aziem, National Research Center, Egypt

Copyright © 2021 Lv, Zhang, Feng, Jiang, Niu, Yang and Wang. 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: Yanjie Yang, yangyj@henu.edu.cn; Xinchun Wang, wxc619533@sina.com

These authors have contributed equally to this work