Abstract
Diarrhea constitutes one of the most common diseases affecting the survival of captive musk deer and is usually caused by an imbalance in intestinal microbiota. Currently, research regarding the structure and function of intestinal microbiota in diarrheic musk deer is lacking. Therefore, in the present study, high-throughput 16S-rRNA gene sequencing was used to analyze the intestinal microbiota in feces of healthy captive musk deer (HMD) (n = 8) and musk deer with mild (MMD) (n = 8), and severe (n = 5) (SMD) diarrhea to compare the difference in intestinal microbiota of musk deer under various physiological conditions. The results showed that the diversity of HMD fecal microbiota was significantly higher than that of the two diarrhea samples. β Diversity results indicated that there were extremely significant differences in bacterial communities between the HMD sample and the MMD and SMD samples. However, no significant difference was found between the two diarrhea samples. LefSe analysis showed that the degree of intestinal physiological dysfunction in musk deer was correlated with the types of major pathogens. The main pathogen in the MMD group is Escherichia-Shigella, whereas Fusobacterium is the main pathogen in the SMD group. PICRUSt functional profile prediction indicated that the intestinal microbiota disorder could also lead to changes in the abundance of genes in metabolic pathways of the immune system. Altogether, this study provides a theoretical basis for the exploration of treatments for diarrhea in captive musk deer, which is of considerable significance to the implementation of the musk deer release into the wild program.