Mice experiments
All animal experiment protocols were approved by the institutional Animal Ethics Committee of Experimental Animal Center of Nanjing Agricultural University (Approval number: PZ2019043, PZ2020024, PZ2020091, and NJAU.No20210425054). The animal experiments were carried out in Animal Center of Nanjing Agricultural University (SYXK<Jiangsu>2011-0037) in accordance with the Guidelines of the Ethical Committee of Experimental Animal Center of Nanjing Agricultural University and the National Guidelines for Experimental Animal Welfare. Mice were housed (two mice per cage) under specific pathogen free (SPF) conditions with free access to sterile food (irradiated) and sterile water in 12 h light/dark cycle. Five or six-week-old C57BL/6J male mice were purchased from Shanghai SLAC Laboratory Animal Co., Ltd. (Shanghai, China). The Tlr4-/- and WT C57BL/10J male mice with five weeks old were purchased from GemPharmatech Co., Ltd. (Nanjing, China). The ND D12450J (10% calories from fat) and HFD D12492 (60% calories from fat) were obtained from Research Diets, Inc. (New Brunswick, NJ, USA), whose compositions are shown in Supplementary Table 1. Before experiment, all animal underwent a 1-week acclimatization period fed with ND. The fluoride was administrated to mice by addition to drinking water, and the drinking water was refreshed twice one week. The food intake and body weight were recorded once a week during the animal experiments.
Animal protocol 1. The purpose of this first mice experiment was to evaluate the effect of fluoride on the obesity in HFD mice. 24 male C57BL/6J mice with 6 weeks of age were randomly divided into three groups including the ND, HFD, and HFD-F groups (HFD plus 50 mg/L of fluoride in drinking water) (n = 8 per group). The intervention continued for 10 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 2. The purpose of this mice experiment was to confirm the effect of fluoride on the obesity in HFD mice in a separate replicate experiment. Furthermore, the potential mechanisms were investigated based on the first mice experiment, which focused on the gut microbiota and intestinal barrier permeability. 24 male C57BL/6J mice with 5 weeks of age were randomly divided into three groups including ND, HFD, and HFD-F groups (n = 8 per group). The intervention continued for 10 weeks. The fecal samples were collected for 16S rRNA sequencing at the last week, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 3. The purpose of this mice experiment was to analyze the effect of fluoride on intestinal barrier permeability in vivo according to the previous work [67]. 16 male C57BL/6J mice with 5 weeks of age were randomly divided into two groups including HFD and HFD-F groups (n = 8 per group). The intervention continued for 10 weeks. At the end of the experiment, FITC-dextran solution (4 kDa (Sigma-Aldrich, St. Louis, MO) with dosage of 0.6 mg/g body weight) was administrated to mice by intragastric gavage after fasting for 4 h. Then, the mice were euthanized after 3 h to afford plasma for measurement of level of FITC-dextran by Fluorescence Spectrophotometer at an excitation wavelength of 485 nm and emission wavelength of 535 nm.
Animal protocol 4. The purpose of this mice experiment was to investigate whether the knockout of Tlr4 could block the exacerbation of obesity by fluoride in HFD mice. Sixteen WT and sixteen Tlr4-/- mice (C57BL/10J, 5 weeks old) were randomly divided into two groups, respectively, including HFD and HFD-F (HFD plus 50 mg/L of fluoride in drinking water) groups (n = 8 per group). The intervention continued for 8 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 5. The purpose of this mice experiment was to investigate whether fluoride-induced gut microbiota could be sufficient to induce some phenotypical changes of obesity by FMT. The mice fed HFD or HFD plus 50 mg/L of fluoride in drinking water described in Animal protocol 2 was chosen as donor. Another 16 male C57BL/6J mice with 5 weeks of age were randomly divided into two groups (n = 8). At first two weeks, all mice were treated by Abx (containing 10 g/L of vancomycin hydrochloride, 20 g/L of neomycin sulfate, 20 g/L of metronidazole and 20 g/L of ampicillin sodium salt) by gavage with dosage of 200 mL for each mouse once daily according to the previous work to obtain the antibiotic-treated mice [49, 68]. At third week, FMT was carried out by transplanting the fecal slurry from donor mice to the antibiotic-treated mice (n = 8) until the end of the experiment as described previously [2, 69]. Briefly, the fecal samples collected from donor mice (HFD and HFD-F groups) were immediately pooled and diluted in sterile saline containing 0.5 g/L cysteine hydrochloride at proportion of 100 mg/mL. Then, the fecal slurry was centrifuged at 1000 r/min for 1 min to obtain fecal supernatant. Each mouse received 200 mL of the fecal supernatant from HFD or HFD-F groups by gavage, respectively. The FMT continued for 8 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 6. The purpose of this mice experiment was to investigate whether fluoride could exacerbate the obesity in the absence of gut microbiota. Sixteen male C57BL/6J mice with 5 weeks of age were randomly divided into two groups including HFD fed with Abx (HFD-Abx) group, and HFD fed with Abx plus 50 mg/L of fluoride in drinking water (HFD-Abx-F) group (n = 8 per group). Abx was treated to mice by gavage once daily to clear gut microbiota according to the previous work with some modifications [68, 70]. The intervention continued for 10 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 7. The purpose of this mice experiment was to investigate whether fluoride could induce the obesity in ND mice. Sixteen male C57BL/6J mice with 5 weeks of age were randomly divided into two groups including ND group, and ND plus 50 mg/L of fluoride in drinking water (ND-F) group (n = 8 per group). The intervention continued for 10 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 8. The purpose of this mice experiment was to investigate whether E. ramosum could exacerbate the obesity in HFD mice (result obtained in animal experiment 2 and 5). Antibiotic-treated mice were obtained by treatment of Abx for first two weeks as described in Animal protocol 5. Then, antibiotic-treated mice were randomly divided into two groups including HFD fed with PBS and HFD fed with E. ramosum (HFD-ER) by intragastric gavage for two weeks. Each mouse received 0.2 mL of live E. ramosum (4 × 108 cfu/mL) daily. The intervention of E. ramosum continued for 2 weeks. After further 4 weeks, the mice were euthanized after fasting overnight at the end of the experiment.
Animal protocol 9. The purpose of last mice experiment was to present a potential strategy to prevent fluoride-induced exacerbation of obesity in HFD mice by TPS. Sixteen male C57BL/6J mice with 5 weeks of age were randomly divided into two groups including HFD plus 50 mg/L of fluoride in drinking water (HFD-F) group, and HFD plus 50 mg/L of fluoride in drinking water and 400 mg/kg/day of TPS (HFD-F-TPS) group by intragastric gavage (n = 8 per group). The mice in HFD-F group received same volume of water by intragastric gavage. The intervention continued for 8 weeks, and the mice were euthanized after fasting overnight at the end of the experiment.
Caco-2 cell monolayer permeability assay in vitro
The Caco-2 cell (HTB-37), obtained from American Type Culture Collection (ATCC, Manassas, VA, USA), was incubated in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% of fetal bovine serum, 1% of penicillin-streptomycin, 1% nonessential amino acids, and 1% L-glutamine at 37 ℃under 5% of CO2 atmosphere. The cytotoxity of fluoride on Caco-2 cell was investigated by MTT assay. Briefly, Caco-2 cell suspension was plated in a 96-well plate (200 mL/well), and cultured for 24 h. Then, the new medium containing various concentrations of fluoride (0, 6.25, 12.5, 25, 50, 100, 200 and 400 mg/mL) were added to replace the old medium, and further incubated for 24 h. After discarding the culture medium, MTT solution was added and incubated for a 4 h. Then, MTT solution was removed and 200 mL of DMSO was added. The cytotoxicity was calculated using the absorbance (Abs) measured at 570 nm. The effect of fluoride on the monolayer permeability was investigated according to the previous work [71, 72]. Caco-2 cell suspension was plated onto the polycarbonate membranes of Millicell inserts into a tween-well Transwell plate (0.4 mm pores; Corning, NY, USA) at a density of 1×105 cells per well. The cell medium in the apical and basolateral sides was renewed every two days at first 14 days and renewed daily at last 7 days. The transepithelial electrical resistance (TEER) was measured by a Millipore ERS-2 system (Millipore Crop., New Bedford, MA, USA). At 21th day, the wells with TEER values of more than 500 Ω cm2 were chosen for further experiments. After washing by PBS, 0.5 mL of cell medium (without fetal bovine serum) containing 0, 10 and 50 mg/L of fluoride were added into the apical side of well in Transwell plate, and 1.5 mL of cell medium without fetal bovine serum were added into the basolateral sides of well in Transwell plate. After incubation at 37 ℃ for 12 h, TEER values were measured again for evaluating the effect of fluoride on the monolayer permeability in vitro.
The lipid metabolism assay in HepG2 cells or primary murine hepatocytes model in vitro
HepG2 cells were purchased from National Collection of Authenticated Cell Cultures (Shanghai, China). HepG2 cells were incubated in DMEM supplemented with 10% of fetal bovine serum, and 1% of penicillin-streptomycin, at 37 ℃under 5% of CO2 atmosphere. The cytotoxity of fluoride on HepG2 cells was investigated by MTT assay as described for Caco-2 cell. The effect of fluoride on the lipid metabolism was evaluated in HepG2 cells in vitro as described previously with some modifications [73]. HepG2 cell suspension was plated in a 12-well plate with 1 mL/well, and cultured for 24 h. Then, the serum-free medium, containing 100 mM of palmitic acid, 200 mM of oleic acid and different dosages of fluoride (0, 10 and 50 mg/L, respectively), was added to replace the old medium, the cells were further incubated for 24 h to induce excess fat synthesis. After washed by PBS for three times, the cells were stained with Oil Red O solution to evaluate the level of intracellular fat droplets. Furthermore, the Oil Red O was extracted by 100% isopropyl alcohol, and the absorbance of the extracted isopropyl alcohol solution was measured at 490 nm to further quantify the intracellular lipid content. Moreover, the triglyceride (TG) and protein contents in HepG2 cells were measured by commercially available kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The primary murine hepatocytes were isolated from the liver of C57BL/6 mice using a collagenase perfusion method according to the reported work [74]. After isolation from mice, the cell suspension was inoculated in a 12-well plate with 1.0 mL/well. After cultured for 24 h, the effect of fluoride on the lipid metabolism in primary murine hepatocytes was investigated as mentioned above for HepG2 cells model.
Cultivation of Erysipelatoclostridium ramosum
E. ramosum CCUG35705 was purchased from Culture Collection University of Gothenburg (CCUG, Sweden). E. ramosum was cultured in brain heart infusion (BHI) medium at 37 ℃ at an anaerobic workstation (HYQX-Ⅲ, Shanghai Yuejin Medical Instrument Co., Ltd., Shanghai, China). After incubation for 12 h, the E. ramosum was centrifuged at 4000 g for 5 min, and resuspended in anaerobic PBS containing 4 × 108 cfu/mL of E. ramosum. Then, the E. ramosum would be gavaged to mice immediately.
Oral glucose tolerance test (OGTT)
One week before the end of animal experiment, the mice were fasted overnight for OGTT. After measurement of blood glucose, the mice were immediately given an oral gavage glucose at a dosage of 1.5 g/kg body weight. Then, the blood glucose was further measured after 15, 30, 60, 90 and 120 min. All blood samples were collected from the tip of the tail vein, and the blood glucose was measured using a glucose meter (Sinocare Inc., Changsha, Hunan, China).
Tissue sampling
The animals were anesthetized with carbon dioxide, and blood, liver, colon, perirenal fat, mesentery fat, and epididymal fat tissues were collected. The blood sample was collected in sterile anticoagulation tube (BD Biosciences, USA). Then, 50 mL of blood was immediately taken out and stored at -80 ℃ for quantification of Bacterial DNA. The residue blood was centrifuged at 3000 r/min for 15 min at 4 ℃ to afford plasma. After weighing and taking picture, a part of the colon and epididymal fat tissues were fixed in 4% paraformaldehyde solution for histological analysis. Other tissues were stored at -80 ℃ for further experiment.
Biochemical analyses
The contents of TG, TC and LDL-C in plasma and TG in liver were measured by commercially available kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The levels of TNF-a, IL-6 and IL-1b in plasma were determined using commercial ELISA kits (Neobioscience Technology Co, Ltd., Shenzhen, China). Plasma LPS was determined using commercial ELISA kit (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The insulin in plasma was determined using commercial ELISA kit (Mercodia AB, Uppsala, Sweden).
Quantification of bacterial DNA in blood
The quantification of bacterial DNA in blood were carried out according to previous work with some modifications [67]. The whole DNA in blood (50 mL) was extracted using a commercially available kit (TIANGEN BIOTECH (BEIJING) CO., LTD., Beijing, China). The concentrations of DNA were analyzed by Nanodrop 2000 (Thermo Fisher Scientific Inc., USA), and the bacterial DNA in blood was quantified by Femto™ Bacterial DNA Quantification Kit (ZYMO RESEARCH, USA).
RNA extraction and Quantitative real-time PCR
The total RNA of liver and colon was extracted by commercially available kits (YFXM0013, YIFEIXUE Biotechnology, Nanjing, China). After quantified by Nanodrop 2000 (Thermo Fisher Scientific Inc., USA), the RNA was reversed to cDNA by PrimeScript RT Master Mix (TaKaRa). The resulting cDNA was used for RT-qPCR using SYBR Green Master Mix (A25742, ABI, USA) and QuantStudio 6 Flex (Thermo Fisher Scientific Inc.). The expression of mRNA was normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by the method of 2-ΔΔCt. Primers were obtained from Sangon Biotech (Shanghai, China), which was presented in Supplementary Table 2.
Histology analysis and immunofluorescence assessment
The histological analyses of epididymal fat and colon were carried out according to the previous work [14, 38]. After fixed in 4% paraformaldehyde solution, the epididymal fat and colon tissues were embedded in paraffin, and stained with hematoxylin and eosin (H&E). The protein expressions of ZO-1, Occludin, Claudin-1 and MUC1 in the colon tissue were also evaluated by immunofluorescence assessment as previously described [74]. After deparaffinized and rehydrated, the colon samples were incubated twice in xylene for 15 min each, dehydrated twice in pure ethanol for 5 min each, dehydrated in gradient ethanol of 85% and 75% ethanol for 5 min each, and washed in distilled water, respectively. The samples were dipped in 3% of bull serum albumin (BSA) for 30 min to block non-specific binding. The BSA solution was discarded, the primary antibody, including Anti-ZO-1 (ab221547), Anti-Occludin (ab216327), Anti-Claudin-1 (ab211737) and Anti-MUC1 (ab45167) (Abcam, Cambridge, UK), diluted in PBS solution was added and incubated at 4 ℃ overnight. After washed by PBS for three times, the sample was further incubated in secondary antibody solution (Jackson ImmunoResearch, Philadelphia, USA) for 50 min in dark condition at room temperature. The sample was washed by PBS three times and incubated with 4',6-diamidino-2-phenylindole (DAPI, Beyotime Biotechnology, Shanghai, China) solution in dark condition at room temperature. After washed by PBS for three times, the sample was added into the spontaneous fluorescence quenching reagent, incubated for 5 min, and further washed in running tap water for 10 min. Then, the samples were detected and collected by Fluorescent Microscopy. The mean density of immunofluorescence assessment was quantified by ImageJ software (1.8.0_172).
Gut microbiota analysis by 16S rRNA amplicon sequencing
The extraction of fecal genomic DNA in fecal samples was carried out using E.Z.N.A. ®Stool DNA Kit (D4015, Omega, Inc., USA). The 16S rRNA gene sequencing was performed by LC-Bio Technology Co., Ltd (Hangzhou, China). Briefly, the variable V3-V4 region of 16S rDNA gene was amplified with the primers 341F CCTACGGGNGGCWGCAG, and 805R GACTACHVGGGTATCTAATCC by PCR using 32 amplification cycles. After detected using 2% of agarose gel electrophoresis, the target fragments were further purified by AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, USA) and quantified by Qubit (Invitrogen, USA). Then, the resulting PCR product was loaded onto NovaSeq PE250 platform using a paired-end sequencing protocol (2 × 250 bp) according to the manufacturer's recommendations. The paired-end reads were identified by unique barcode of each sample. After cutting off the barcode and primer sequence, paired-end reads were merged by FLASH (v1.2.8). According to the fqtrim (v0.94), the high-quality clean tags was obtained after the quality filtering of raw reads under specific filtering conditions. Vsearch software (v2.3.4) was used to filter chimeric sequences. Then, the feature table and feature sequence were obtained after dereplication by DADA2. The alpha diversity and beta diversity were analyzed using QIIME2. The sequence alignment of species was annotated by Blast based on databases of SILVA (Release 132, https://www.arb-silva.de/documentation/release-132/) and NT-16S. b-Diversity including PCA and PCoA was performed using the OmicStudio tools at https://www.omicstudio.cn/tool. LEfSe analysis was used to detect differentially abundant taxa in each group using the different parameters values of Kruskal-Wallis, Wilcoxon and linear discriminant analysis (LDA) score.
Quantification of E. ramosum in fecal samples
The quantification of E. ramosum in fecal samples was carried out according to the previous work [75]. Briefly, DNA in fecal samples was extracted by a commercially available kit (TIANGEN BIOTECH (BEIJING) CO., LTD.). The levels of DNA were quantified by Nanodrop 2000 (Thermo Fisher Scientific Inc., USA). Gene levels of E. ramosum were determined by RT-qPCR using SYBR Green Master Mix (A25742, ABI, USA) and QuantStudio 6 Flex (Thermo Fisher Scientific Inc.). 50 ng of bacterial DNA was used to perform the RT-qPCR for each sample. The primer of E. ramosum used in this work was as fellow: FW: ACAATGGATGGTGCAGAGGG, RV: TCAACTCTCTCGTGGTGTGACG [61]. Higher CT values suggest lower levels of E. ramosum.
Human cohort analyses
The data from the American Gut Project (AGP) was reanalyzed to gain insight into the relationship between the relative abundance of Erysipelotrichaceae and obesity according to the previous work [76, 77]. The subjects which contained fecal samples were chosen. The subjects which contained fecal samples were chosen. Subjects who were treated with antibiotic or had diabetes were excluded, and the samples with missing information such as BMI and sex were removed. Furthermore, the subjects whose BMI were more than 50 or less than 18.5 kg/m2 were also excluded. Finally, a total of 10,376 individuals aged from 20 to 99 years, including 6546 normal, 2869 overweight and 961 obesity subjects, were obtained in this work. The sequence read files of these subjects were obtained from European Bioinformatics Institute (EBI, PRJEB11419) [33]. The sequences were analyzed using pipeline as previously described [78] using Rstudio (1.4.1106), vsearch (v2.15.2) and usearch (v10.0.240). The ASV table was obtained by usearch, and the sequence alignment of species was annotated by vsearch according to the databases of silva_16s_v123. The subjects with less than 5000 clean reads were removed. After merging all annotation of ASVs to taxonomy tables, the relative abundance of Erysipelotrichaceae was obtained for each subject.
Statistical analysis.
The outliers in this work were removed by GraphPad Prism 9.3.1 based on the Grubbs test (San Diego, CA, USA). The results are expressed as means ± SEM. The normality of all data was checked by SPSS 22 software (IBM) according to Shapiro-Wilk test. If the data for multiple-group comparisons had normal distribution, statistical significance was carried out by one-way analysis of variance (ANOVA) with Tukey test, otherwise, the statistical significance was calculated using the Mann-Whitney test. Difference in two groups was calculated using the Mann-Whitney test or Kruskal-Wallis test. Statistical analysis was carried out by SPSS 22 software (IBM) and GraphPad Prism 9.3.1. Furthermore, adjusted p-values (q-value) were used to evaluate differences in analysis of gut microbiota based on false discovery rate (FDR) for multiple testing according to the Benjamini and Hochberg procedure [79]. A value of p or q < .05 was considered to be significant.