Statistical analysis
Data were represented as mean ± standard deviation (SD) or count (%). All data were input into SPSS (version 25.0) to analyze. GraphPad prism (version 8.0) was used to draw diagrams. χ2 and Fisher’s exact test was used for categorical variables, and t-test or non-parametric Wilcoxon test was used for continuous variables where appropriate. P < 0.05 was considered to be statistically significant. And bioinformatics analysis for microbiome used R software (Bell Laboratories). Alpha and beta diversities were generated in the Quantitative Insights Into Microbial Ecology (QIIME) and calculated based on weighted or unweighted Unifrac distance matrices. We used the linear discriminant analysis (LDA) effect size (LEfSe) method to identify species that show statistically significant differential abundances between groups.