Transcriptomics reveals dynamic changes in the “gene profiles” of rat supraspinatus tendon at three different time points after diabetes induction

Abstract

Abstract Objective There is increasing evidence that type 2 diabetes mellitus (T2DM) is an independent risk factor for the occur of tendinopathy. Therefore, this study is the first to explore the dynamic changes of the “gene profile” of supraspinatus tendon in rats at different time points after T2DM induction through transcriptomics, providing potential molecular markers for exploring the pathogenesis of diabetic tendinopathy. Methods A total of 40 Sprague-Dawley rats were randomly divided into normal (NG, n = 10) and T2DM groups (T2DM, n = 30) and subdivided into three groups according to the duration of diabetes: T2DM-4w, T2DM-8w, and T2DM-12w groups; the duration was calculated from the time point of T2DM rat model establishment. The three comparison groups were set up in this study, T2DM-4w group vs. NG, T2DM-8w group vs. NG, and T2DM-12w group vs. NG. Differentially expressed genes (DEGs) in 3 comparison groups were screened. The intersection of the three comparison groups’ DEGs was defined as key genes that changed consistently in the supraspinatus tendon after diabetes induction. Cluster analysis, gene ontology (GO) functional annotation analysis and Kyoto encyclopedia of genes and genomes (KEGG) functional annotation and enrichment analysis were performed for DEGs. Results T2DM-4w group vs. NG, T2DM-8w group vs. NG, and T2DM-12w group vs. NG detected 519 (251 up-regulated and 268 down-regulated), 459 (342 up-regulated and 117 down-regulated) and 328 (255 up-regulated and 73 down-regulated) DEGs, respectively. 103 key genes of sustained changes in the supraspinatus tendon following induction of diabetes, which are the first identified biomarkers of the supraspinatus tendon as it progresses through the course of diabetes.The GO analysis results showed that the most significant enrichment in biological processes was calcium ion transmembrane import into cytosol (3 DEGs). The most significant enrichment in cellular component was extracellular matrix (9 DEGs). The most significant enrichment in molecular function was glutamate-gated calcium ion channel activity (3 DEGs). The results of KEGG pathway enrichment analysis showed that there were 17 major pathways (p < 0.05) that diabetes affected supratinusculus tendinopathy, including cAMP signaling pathway and Calcium signaling pathway. Conclusions Transcriptomics reveals dynamic changes in the“gene profiles”of rat supraspinatus tendon at three different time points after diabetes induction. The 103 DEGs identified in this study may provide potential molecular markers for exploring the pathogenesis of diabetic tendinopathy, and the 17 major pathways enriched in KEGG may provide new ideas for exploring the pathogenesis of diabetic tendinopathy.