Abstract:Kuihan Tang, 1 Li Sun, 2 Long Chen, 2 Xiaobo Feng, 1 Jiarui Wu, 2 Hao Guo, 2 Yong Zheng 1 1 Department of Orthopedics, Beijing Jishuitan Hospital Guizhou Hospital, Guiyang, 550014, People's Republic of China; 2 Department of Orthopedics, Guizhou Provincial People's Hospital, Guiyang, 550000, People's Republic of China Correspondence: Yong Zheng, Email Background: Mitochondrial autophagy is closely related to the pathogenesis of osteoarthritis, In order to explore the role of mitochondrial autophagy related genes in knee osteoarthritis (KOA) and its molecular mechanism. Methods: KOA-related transcriptome data were extracted from the Gene Expression Omnibus (GEO) database. Differentially expressed mitochondrial autophagy gene (DEMGs) were screened in patients with KOA by differential expression analysis. The STRING website was used to construct a protein-protein interaction (PPI) network among DEMGs. Molecular complex detection (MCODE) method in Cytoscape software was performed to identify hub DEMGs. Support vector machine recursive feature elimination (SVM-RFE) method was used to construct the hub DEMG diagnosis model. Genes with diagnostic value were identified as biomarkers by plotting receiver operating characteristic (ROC) curves and Expression validation. CIBERSORT algorithm was used to calculate the proportion of 22 immune cells in each sample in the GSE114007 dataset. Finally, biomarker expression was verified by qPCR. Results: A total of 15 DEMGs were obtained and enrichment analyses showed that these DEMG strains were mainly enriched in the mitophagy-animal, shigellosis, autophagy-animal and FoxO signal pathways. The PPI network unveiled 13 DEMGs with interactions. In addition, 8 hub DEMGs (ULK1, CALCOCO2, MAP1LC3B, BNIP3L, GABARAPL1, BNIP3, FKBP8 and FOXO3) were obtained for KOA. And 5 model DEMGs (BNIP3L, BNIP3, MAP1LC3B, ULK1 and FOXO3) were screened. The ROC curves revealed that BNIP3 and FOXO3 has strong diagnostic value in these models of DEMG. Immune-infiltration and correlation analysis showed that BNIP3 and FOXO3 were significantly correlated with three different immune cells, including primary B cells, M0 macrophage and M2 macrophage. The cartilage tissue samples qPCR verification results show that FOXO3 and BNIP3 were all down-regulated in KOA ( p < 0.01), and the validation results are consistent with the above analysis. Conclusion: BNIP3 and FOXO3 have been identified as biomarkers for the diagnosis of KOA, which might supply a new insight for the pathogenesis and treatment of KOA. Keywords: Gene Expression Omnibus, diagnostic, bioinformatics analysis, immune infiltration, biomarkers Knee Osteoarthritis (KOA) is a chronic and progressive disease characterized by loss of articular cartilage, subchondral sclerosis, and synovial and periarticular structural abnormalities that cause pain, swelling and even joint deformity in the surrounding tissues of the knee. 1 Its pathogenic factors include abnormal biologic factors, genetic factors, cell aging and apoptosis, local inflammatory factors, free radicals and proteases, etc. Each of these risk factors may cause KOA, KOA is a degenerative, attritional disease, 2 its pathology mainly involves articular cartilage and synovium, the key factor of KOA cartilage degeneration is chondrocyte death, 3 which is caused by apoptosis and autophagy. 4,5 Under normal conditions, articular chondrocytes maintain a dynamic balance between extracellular matrix (ECM) synthesis and degradation, 6 which is broken in response to injury or inflammation and other stimuli, triggering chondrocyte stress and ECM degradation, it leads to abnormal accumulation of damaged proteins and dysfunction of organelles such as endoplasmic reticulum and mitochondria, 7,8 thus, cartilage degradation, chondrocyte apoptosis, subchondral bone dysfunction, and gradually develop into knee osteoarthritis. 9 According to the latest Global Burden of Disease study, about 250 million people worldwide suffer from osteoarthritis, and knee joint is the most common site of osteoarthritis. 2 Due to the lack of early diagnostic indicators, patients with KOA often miss the opportunity for optimal treatment, resulting in a poor prognosis. Mitochondria are membrane-coated organelles present in eukaryotic cells that regulate important cell functions and cell survival. Mitochondria have many important functions, including cell adenosine triphosphate (ATP) production, Ca2+ buffering, and being the main source of endogenous reactive oxygen species (ROS) under oxidative stress. 10 Mitochondrial autophagy is a major intracellular prote -Abstract Truncated-

Identification of genetic association between mitochondrial dysfunction and knee osteoarthritis through integrating multi-omics: a summary data-based Mendelian randomization study

Characterizing mitochondrial features in osteoarthritis through integrative multi-omics and machine learning analysis

Machine learning algorithm-based biomarker exploration and validation of mitochondria-related diagnostic genes in osteoarthritis

Causality of Genetically Determined Metabolites and Metabolic Pathways on Osteoarthritis: A Bi-Directional Two-Sample Mendelian Randomization Study

Proteome-Wide Mendelian Randomization and Colocalization Analysis Identify Therapeutic Targets for Knee and Hip Osteoarthritis

Evidence supported by Mendelian randomization: impact on inflammatory factors in knee osteoarthritis

Comprehensive Analysis of Key Genes, Signaling Pathways and miRNAs in Human Knee Osteoarthritis: Based on Bioinformatics

Bioinformatics Analysis and Experimental Validation of Mitochondrial Autophagy Genes in Knee Osteoarthritis

Identification and Verification of Hub Mitochondrial Dysfunction Genes in Osteoarthritis Based on Bioinformatics Analysis

Causal relationship between modifiable risk factors and knee osteoarthritis: a Mendelian randomization study

Causal relationship between controllable risk factors and knee osteoarthritis: a Mendelian randomized study

Gut microbiota metabolic pathways: Key players in knee osteoarthritis development

Identification of biomarkers and potential drug targets in osteoarthritis based on bioinformatics analysis and mendelian randomization

The role of mitochondrial autophagy in osteoarthritis

Causal link between metabolic related factors and osteoarthritis: a Mendelian randomization investigation

Mitochondrial dysfunction and onset of type 2 diabetes along with its complications: a multi-omics Mendelian randomization and colocalization study

Identification of key hub genes in knee osteoarthritis through integrated bioinformatics analysis

Mendelian randomization based on genome-wide association studies and expression quantitative trait loci, predicting gene targets for the complexity of osteoarthritis as well as the clinical prognosis of the condition

Mendelian randomization and transcriptome analysis identified immune-related biomarkers for osteoarthritis

Circulating cytokines levels and osteoarthritis: A Mendelian randomization study

The causal relationship between human brain morphometry and knee osteoarthritis: a two-sample Mendelian randomization study