Analysis of gut microbiome composition, function, and phenotype in patients with osteoarthritis

Su Liu; Guoqing Li; Yuanchao Zhu; Chang Xu; Qi Yang; Ao Xiong; Jian Weng; Fei Yu; Hui Zeng

doi:10.3389/fmicb.2022.980591

Analysis of gut microbiome composition, function, and phenotype in patients with osteoarthritis

Front Microbiol. 2022 Nov 25:13:980591. doi: 10.3389/fmicb.2022.980591. eCollection 2022.

Authors

Su Liu^{1

2}, Guoqing Li^{1

2}, Yuanchao Zhu^{1

2}, Chang Xu³, Qi Yang⁴, Ao Xiong^{1

2}, Jian Weng^{1

2}, Fei Yu^{1

2}, Hui Zeng^{1

2}

Affiliations

¹ Department of Bone and Joint Surgery, Peking University Shenzhen Hospital, Shenzhen, China.
² National and Local Joint Engineering Research Center of Orthopaedic Biomaterials, Peking University Shenzhen Hospital, Shenzhen, China.
³ Peking University Shenzhen Hospital Intelligent Hospital Research Academy, Shenzhen, China.
⁴ Department of Ultrasonography, Peking University Shenzhen Hospital, Shenzhen, China.

Abstract

Gut microbiome (GMB) disturbance can induce chronic low-grade inflammation, which is closely related to the occurrence and development of osteoarthritis (OA). However, the relationship between GMB and OA remains unclear. In this study, we collected stool samples from OA patients and healthy people, and performed Alpha diversity, Beta diversity, MetaStat, and LEfSe analysis by 16S rRNA sequencing to find out the species with significant difference between the two groups. Random forest analysis was performed to find out biomarkers that could distinguish between OA patients and healthy people. PICRUSt and Bugbase analysis were used to compare the difference in functions and phenotypes. Multivariate linear regression analysis (MaAsLin) was used to adjust for gender, age, and body mass index (BMI). The results showed that there was a significant difference in the overall composition of GMB between the two groups (p = 0.005). After adjusting for gender, age, and BMI, we found that p_Bacteroidota (Q = 0.039), c_Bacteroidia (Q = 0.039), and o_Bacteroidales (Q = 0.040) were enriched in the OA group, while s_Prevotella_copri (Q = 0.001) was enriched in the healthy control group. Prevotella could distinguish between OA patients and healthy people with a better diagnostic power (AUC = 77.5%, p < 0.001, 95% CI: 66.9-88.1%). The functions of DNA transcription, amino acid metabolism (including histidine, lysine, and isoleucine), ATP metabolism, and phospholipid metabolism significantly decreased, while glucose metabolism, protein acetylation, and aspartate kinase activity significantly increased in the OA group. In terms of phenotypes, we found that the relative abundance of aerobic (p = 0.003) and Gram-negative (p < 0.001) was higher in the OA group, while contains mobile elements (p = 0.001) and Gram-positive (p < 0.001) were higher in the healthy control group. Our study preliminarily demonstrated that there were differences in the composition, function, and phenotype of GMB in stool samples between OA patients and healthy people, which provided a novel perspective on further study in OA.

Keywords: 16S rRNA sequencing; gut microbiome; inflammation; osteoarthritis; stool.