Analysis of differentially expressed genes in torn rotator cuff tendon tissues in diabetic patients through RNA-sequencing

Abstract

Background Rotator cuff tears (RCT) is a common musculoskeletal disorder in the shoulder which cause pain and functional disability. Diabetes mellitus (DM) is characterized by impaired ability of producing or responding to insulin and has been reported to act as a risk factor of the progression of rotator cuff tendinopathy and tear. Long non-coding RNAs (lncRNAs) are involved in the development of various diseases, but little is known about their potential roles involved in RCT of diabetic patients. Methods RNA-Sequencing (RNA-Seq) was used in this study to profile differentially expressed lncRNAs and mRNAs in RCT samples between 3 diabetic and 3 nondiabetic patients. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis were performed to annotate the function of the differentially expressed genes (DEGs). LncRNA-mRNA co-expression network and competing endogenous RNA (ceRNA) network were constructed to elucidate the potential molecular mechanisms of DM affecting RCT. Results In total, 505 lncRNAs and 388 mRNAs were detected to be differentially expressed in RCT samples between diabetic and nondiabetic patients. GO functional analysis indicated that related lncRNAs and mRNAs were involved in metabolic process, immune system process and others. KEGG pathway analysis indicated that related mRNAs were involved in ferroptosis, PI3K-Akt signaling pathway, Wnt signaling pathway, JAK-STAT signaling pathway and IL-17 signaling pathway and others. LncRNA-mRNA co-expression network was constructed, and ceRNA network showed the interaction of differentially expressed RNAs, comprising 5 lncRNAs, 2 mRNAs, and 142 miRNAs. TF regulation analysis revealed that STAT affected the progression of RCT by regulating the apoptosis pathway in diabetic patients. Conclusions We preliminarily dissected the differential expression profile of lncRNAs and mRNAs in torn rotator cuff tendon between diabetic and nondiabetic patients. And the bioinformatic analysis suggested some important RNAs and signaling pathways regarding inflammation and apoptosis were involved in diabetic RCT. Our findings offer a new perspective on the association between DM and progression of RCT.