Abstract
BACKGROUND: Heterotopic ossification is a frequent complication of soft tissue injuries, particularly in tendons. Although ossification in tendon tissue has been reported in a range of aging and disease models, the underlying biomarkers and mechanisms remain unknown. And the characterisation and sensitivity of previous diagnostic biomarkers for tendon ectopic ossification do not meet the demands of clinical use. The aim of this study was to characterise the effects of aging on ossification in the mouse Achilles tendon and to identify characteristic genes and therapeutic targets for tendon ossification in mice by using a machine learning approach. METHODS: We retrieved the transcriptome profile of GSE126118 from the Gene Expression Omnibus (GEO) database. Following background correction and normalization using the transcripts per million (TPM) method, differentially expressed genes (DEGs) were identified with the limma R package (p < 0.05, |log2FC| > 1). Subsequently, 468 senescence genes were downloaded from the Aging Atlas database, and senescence-associated DEGs (HO senescence genes) were identified. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and protein-protein interaction (PPI) network analyses were conducted on the identified DEGs. To further refine the HO aging signature, support vector machine (SVM) regression was employed. Additionally, we predicted transcription factors, miRNAs, and small molecule drugs potentially associated with the characterized genes. RESULTS: Three characterised genes were identified as biomarkers associated with ectopic ossification and aging in the mouse Achilles tendon, Atp5o, Mmp2 and Mmp13. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed significant enrichment in processes related to cartilage endochondral ossification, metalloendopeptidase activity, and mitochondrial proton transport ATP synthase complex. Additionally, HIF-1 and GnRH signaling pathways were prominently represented among the differentially expressed genes. CONCLUSION: Atp5o, Mmp2 and Mmp13 were identified as relevant signature genes for the effects of aging on Achilles tendon ossification in mice. Atp5o, Mmp2, and Mmp13 may influence tendon ossification by affecting mitochondrial function as well as extracellular matrix degradation to regulate senescence. This finding suggests a potential link between these processes, opening new avenues for research into diagnostic markers and therapeutic targets. These genes hold promise for the development of novel treatments for tendon ossification, a debilitating condition currently lacking effective therapeutic options.