Experiment materials and reagents preparation
The basal feed used in our study was purchased from Hunan Slaccas Jingda Laboratory Animal Company (Hunan, China). The high-sugar and high-fat feed (80% basal feed mixed with 12% lard and 8% honey) was purchased from Jiangsu Synergetic Pharmaceutical Biological Engineering Co., LTD with license number Su Feeding Certificate (2014) 01008. Hongxing Erguotou (56 Degrees) was produced by Beijing Hongxing Co., LTD.
Animals
Ten one-week-old specific pathogen free (SPF) male Kunming mice, weighing 18−20 g, were purchased from Hunan Slaccas Jingda Laboratory Animal Company (Hunan, China) with license number SCXK (Xiang) 2016-0002. Experimental animals were kept in the Experimental Animal Center of Hunan University of Chinese Medicine (Experimental Unit Use License Number: SYXK (Xiang) 2019-0009) under the controlled conditions (temperature 23-25°C, humidity 50-70%). All procedures involving animals were performed according to protocols approved by the Institutional Animal Care and Use Committee of Hunan University of Chinese Medicine. The trial procedures for the care and use of experimental animals were carried out in accordance with the European Community guidelines (directive 2010/63/EU).
Animal experimental process
After a 3-days acclimatization period, 10 mice were randomly distributed into the control (ccm) group and model (cmm) group, with 5 mice in each group. We have improved the modeling method of diarrhea with IDHS construction in reference literature[25]: mice in the cmm group were fed with high sugar and high fat feed continuously 11 days, whereas the ccm mice were fed with basal feed. Starting from day 12, mice in the cmm group were placed in an artificial climate chamber at 31.5 °C-32.5 °C and 95 % relative humidity for 8 h/d. At the same time, the cmm group were gavaged 0.35 mL of Hongxing Erguotou at 9 am and 8 pm every day, and gavaged 0.35 mL of ice water at 3 pm, lasted 4 days. The ccm group was administered with equal amount of distilled water.
Model evaluation criteria
We referred to the "Consensus Opinions of Diarrhea Traditional Chinese Medicine Diagnosis and Treatment Experts (2017)" on the main points of the differentiation of the diarrhea with IDHS, and formulated the general behavioral diagnostic criteria of the diarrhea with IDHS model[3]: ①increased defecation frequency; ②yellowish-brown loose stools; ③foul-smelling stools; ④perianal soiling.
Collection of mice intestinal mucosa
Mice in each group were sacrificed by cervical dislocation, immediately placed on an ultra-clean bench and whole sections of small intestine were taken. The small intestine was dissected, the contents were washed from the intestinal wall with saline and the water was blotted dry with sterile filter paper. We scraped the intestinal mucosa with a coverslip, weighed and labelled it in the sterile eppendorf (EP) tubes, added 2 times the weight of the intestine in saline and stored it at -80°C in the refrigerator.
16S rRNA gene high-throughput sequencing
Total metagenomic DNA of intestinal mucosal bacteria were extracted with CTAB and SDS. Then, the quantity and quality of the extracted DNA were measured using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and 3.0 Flurometer Qubit (Life Technologies, CA, USA) and agarose gel electrophoresis method. The full length of bacterial 16S rDNA gene sequence was amplified using the extracted DNA as template. The primers used for amplification were 27F (5'-AGRGTTYGATYMTGGCTCAG-3') and 1492R (5'-RGYTACCTTGTTACGACTT-3'). PCR amplification system was prepared as follows (25 L): 5 μL of KAPA HiFi Buffer (5 ×), 0.75μL (1 U/ μL) of KAPA HiFi Hot Start DNA Polymerase, 0.75 μL (10 mM) of dNTPs, 0.75 μL (10 uM) each of the forward and reverse primers, 2 μL of DNA template, and 15 μL of ddH2O. PCR amplification conditions were as follows: initial denaturation at 95 °C for 5 sec, followed by denaturation at 95 °C for 30 sec, annealing at 57 °C for 30 sec, extension at 72 °C for 60 sec, a total of 25 cycles, and a final extension at 72 °C. PCR amplicons were purified with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN) and quantied using the PicoGreen dsDNA AssayKit (Invitrogen, Carlsbad, CA, USA). After the individual quantification step, amplicons were pooled in equal amounts. The pooled sample was then used to generate a library by using SMRTbell Template Prep Kit 1.0-SPv3, and the PCR products were then sequenced by PacBio Platform using DNA/Polymerase Binding Kit 3.0 (PacBio). The extraction, amplification, database construction and sequencing of the sample DNA were all completed by Wuhan Frasergen Bioinformatics Co., LTD.
Bioinformatics and statistical analysis
Bioinformatics analysis
Based on the SMRT (Single Molecule Real-Time) single molecular real-time sequencing technology and the PacBio Sequel sequencing platform, the raw data were stripped of the adapter sequences to obtain validated inserts. The software SMRT Link v8.0 was used to pre-process and filter the raw sequencing output data to obtain Circular Consensus Aequences (CCS), i.e. raw reads. Finally, clean reads were obtained for subsequent analysis by primer removal and length filtering (1300-1600bp) performed on raw read (the corresponding raw data has been uploaded to National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/) database and the accession number is PRJNA673123). Qiime software (v.1.8.0, http://www.qiime.org/) was used to conduct OTU clustering with 97% identity for Clean CCS sequences of all samples. The longest sequence in each OTU was selected as the representative sequence of the OTU, and the sequence with the similarity over 97% with the representative sequence of the OTU was selected to generate the original OTU table. Qiime was used to flatten the original OTU table according to the sample with the lowest data volume, so that all the samples had the same data volume after flattening, and the final OTU table was obtained. Community diversity was reflected by a variety of different indexes, which measured the diversity with different emphasis. Chao 1 and Observed species index focused on estimating the abundance of bacteria, while Simpson and Shannon indexes were usually used to describe the diversity of the community. In this study, Chao 1, Observed species, Simpson and Shannon indexes were calculated by MOTHUR (version v.1.30.1, http://www.mother.org/). Beta diversity was used to investigate the structural variation of microbial communities across samples using UniFrac distance metrics and used R software to analyze Principal coordinate analysis (PCoA). Linear discriminant analysis effect size (LEfSe) was used to screen key biomarkers based on linear discriminant analysis (LDA) effect size.
Statistical analysis
SPSS 24.00 software (IBM Corp., Armonk, NY, USA) was used for statistical analysis, and data were expressed as mean ± standard deviation. Data were tested for normality using Shapiro-Wilk test and for variance homogeneity by Levene test. Means between the two groups conformed to normality and chi-squared were compared using the independent samples t-test, otherwise the non-parametric Mann-Whitey test was applied. Then, p<0.05 was regarded as a significant difference, p<0.01 as the extremely significant difference.