Skip to main content

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

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.

Peer Review reports

Introduction

Tendinopathy is commonly seen after degenerative changes in the tendons, often characterized by pain and limited movement, and its pathogenesis is still unclear [1]. In recent years, relevant studies have shown that diabetes is closely related to tendon injury, and diabetes increases the risk of tendon injury. Longo et al. [2] retrospectively analyzed 194 patients who underwent arthroscopic surgery from 2007 to 2008, and found that patients with diabetes were more likely to have rotator cuff injury, and hyperglycemia was a risk factor for rotator cuff injury (P = 0.007). Spoendlin et al. [3] divided 7,895 patients with achilles tendon injury or biceps tendon injury from 1995 to 2013 into diabetic group and non-diabetic group, and found through comparative analysis that the degree of control of blood sugar level in female patients was closely related to tendon injury. Park et al. [4] selected 634 subjects diagnosed with rotator cuff tears and found that diabetes was one of the risk factors for rotator cuff tears. Other studies have found that the likelihood of tendinopathy in patients with diabetes is 4 times that of nondiabetic, and the likelihood of tendon tear or rupture is 5 times that of nondiabetic [5].

Previous studies have confirmed that diabetes promotes tendon degeneration and promotes the occurrence of tendinopathy. Batista et al. [6] found that 62 of the 70 patient with diabetes had tendon fiber disorder and 17 had tendon calcification through color ultrasound. Another study [7] conducted shoulder joint ultrasound examination on 48 elderly patient with diabetes without shoulder pain or dysfunction and 32 control patients. The results showed that the mean thickness of supraspinatus muscle tendon in patient with diabetes was 6.2 mm, which was higher than 5.2 mm in control group (P < 0.05). At the same time, the results of ultrasound examination showed that 42.7% of patient with diabetes with supraspinatus tendon degenerative changes, much higher than the control group. All the above experiments showed that diabetes was closely related to the occurrence of tendinopathy. Diabetes damages the normal structure of tendons from many aspects such as tendon collagen, blood vessels and nerves, and then changes the biomechanical properties of tendons. However, many of the specific influencing mechanisms are still unclear and need more in-depth research.

A comprehensive understanding of the variations occurring in diseased tendons is essential to prevent the onset and progression of tendinopathy and to relieve symptoms. Due to continuous technological innovation and continuous cost reduction, transcriptomic sequencing has been increasingly used as a powerful alternative to bulk RNA sequencing over the past few years. Transcriptomics is based on Illumina sequencing platform to study all mRNA transcribed by a specific tissue or cell ina certain period, which is the basis of gene function and structure research, and plays an important role in understanding the development of organisms and the occurrence of diseases. To date, transcriptome changes in injured tendons have been reported in several studies using bulk RNA sequencing [8,9,10]. In addition, single-cell and spatial RNA sequencing methods were used to further understand the cellular heterogeneity and molecular mechanism of tendinopathy progression, and the study results provided new insights into the control of tendinopathy [11]. Single-cell and spatial transcriptomics have shown that the pathogenesis of human tendinopathy is closely related to dysregulation of immune homeostasis [12]. Multi-omics techniques have been widely used to explore the pathogenesis of tendinopathy, but until now, transcriptomic techniques have not been applied to the study of diabetic tendinopathy.

In previous studies, we successfully established an animal model of rotator cuff tendinopathy in diabetic rats by comparing the histological and biomechanical changes of supraspinatus tendon at 2,4,8 and 12 weeks after diabetes induction [13]. More and more studies have explored the changes of gene transcription in diabetes-related complications through transcriptomics technology, and most of them have used rat models as experimental objects [14,15,16,17], so the rat model can be used as experimental objects close to the human model. In this study, we collected tendons of the supraspinatus tendon from T2DM and normal rats and then performed transcriptome sequencing. By constructing the of the supraspinatus muscle tendon of diabetic rats, we found the dynamic change of the “gene profile” in the occur of rotator cuff tendinopathy with the progression of diabetes, and found the key genes and potential pathogenic pathways that may be closely related to the pathogenic mechanism. These findings will help us to more thoroughly understand the pathogenesis of diabetes leading to supraspinatus tendinopathy and provide potential targets for the treatment of diabetic tendinopathy.

Methods

Experimental animals

SPF-grade healthy male Sprague-Dawley rat, aged 8 weeks and weighing 220  250 g, was provided by Animal Research Center of Qingdao University. Animal License Number: SCXK (Zhejiang) 2019-0001. This experiment was approved by the Animal Ethics Committee of Qingdao University (Approval No. 20220505SD8020221210126). The study is reported in accordance with ARRIVE guidelines.

Animal grouping

A total of 40 Sprague-Dawley rats were randomly divided into normal (NG, n = 10) and type 2 diabetes mellitus 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. High-fat diet(#MD12033; MediScience Diets Co., Ltd., Yangzhou, China)was fed for 2 weeks. The rats were housed in clean conditions with temperatures maintained at 20–24 °C. All rats ate 3 regular meals with 12 h of light per day, and turned off the light and stopped providing food at night. After 2 weeks, rats were deprived of water for 12 h and intraperitoneally injected with Streptozotocinn (STZ; Aladdin Co., Ltd., Beijing, China). After 72 h, tail venous blood was taken to detect fasting blood glucose of rats [18]. If the blood glucose level was ≥ 16.7 mmol/l for 3 consecutive days, the diabetic rat model was considered to be successfully established [19]. The normal group was fed with ordinary diet for 2 weeks, and after 2 weeks, the same amount of sodium citrate buffer ( BIOISCO Co., Ltd., Nanjing, China).

Obtain sample

At the 4th, 8th and 12th week after modeling, the bilateral supraspinatus muscle of diabetic rats was completely dissected, the excess muscle tissue around the tendon was separated, and the complete supraspinatus tendon was obtained, rinsed with phosphate-buffered saline (PBS), and the tissue samples were stored in a -80°refrigerator after drying. The tendons of supraspinatus muscle of normal rats were taken by the same method. In order to meet the requirements of transcriptomics for sample quality, supraspinatus tendon tissues of two animals (four shoulders) were combined to construct a transcriptomic sample.

RNA extraction and bioinformatics analysis of RNA-seq

Total RNA was isolated and purified using Trizol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). After quality control, libraries were constructed as previously described [20]. mRNA-seq were performed on the Hiseq4000 platform (Illumina, CA, USA) with 150 strategy. The method of RNA-seq bioinformatics analysis was consistent with previous studies [21]. The differentially expressed mRNAs were selected with fold change>2 or fold change<0.5 and with parametric F-test comparing nested linear models (adjusted p-value<0.05) by R package edgeR (https://bioconductor.org/packages/release/bioc/html/edgeR.html). Sample processing, computer sequencing and data visualization analysis were all assisted by Lc-Bio Technologies (Hangzhou, China).

Results

Three comparison groups of DEGs

20 supraspinatus tendon samples were extracted from normal and diabetic rats in this study to reveal a list of genes associated with the occur of diabetic tendinopathy (Fig. 1). In the process of data analysis, we found that two samples had little correlation with other samples in the group, and we eliminated these two samples. Therefore, the T2DM-4w and T2DM-12w groups are 4 samples. In order to determine the effect of T2DM on the transcriptome of supraspinatus tendon in rats and its dynamic change process, the “gene profile” of supraspinatus tendon was compared 4, 8 and 12w after T2DM induction and in the control group. 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. The upregulated and downregulated DEGs identified in the 3 comparison groups were ranked in order of change multiples from largest to smallest, and the top 10 DEGs induced by T2DM were listed (Table 1). Cluster analysis was performed on DEGs of the three comparison groups. The cluster heat map (Fig. 2) and volcano map (Fig. 3A-C) clearly showed the distribution of DEGs at different time points after T2DM induction.

Fig. 1
figure 1

Study design and transcriptomics analysis of the supraspinatus tendon in diabetic rats. Overview of the cohort (including 5 normal samples, 5 type 2 diabetes 4-week samples, 5 type 2 diabetes 8-week samples, and 5 type 2 diabetes 12-week samples) and study design (including RNA extraction, library construction, bioinformatics analysis of RNA-seq, and screening of DEGs)

Table 1 Three comparison groups up-regulated and down-regulated DEGs (Top10)
Fig. 2
figure 2

Heatmap of annotated genes with increasing and decreasing trend, including T2DM-4w group vs. NG (A), T2DM-8w group vs. NG (B), and T2DM-12w group vs. NG (C). Each column represents a sample, and each gene is visualized in a row red indicates a high abundance, and blue indicates a relatively low abundance of genes

Fig. 3
figure 3

Volcano plot of DEGs between three comparison groups. Each dot represents one gene red dots represent the significantly upregulated genes, and blue dots represent the significantly downregulated genes gray dots represent no significant DEGs. (A) T2DM-4w group vs. NG: (B) T2DM-8w group vs. NG: (C) T2DM-12w group vs. NG. According to the principle of Gene ontology (GO) analysis, the biological function of genes is divided into three main parts, namely biological process (BP), cellular component (CC) and molecular function (MF). (D) T2DM-4w group vs. NG: (E) T2DM-8w group vs. NG: (F) T2DM-12w group vs. NG

GO enrichment analysis

In order to study the potential biological functions and pathways of tendinopathy caused by diabetes, GO enrichment analysis of these DEGs was further conducted (Fig. 3D-F). According to the principle of GO analysis, the biological function of genes is divided into three main parts, namely biological process, cellular component and molecular function. In the comparison of T2DM-4w and NG, and the biological processes, cellular component and molecular function enrichment of DEGs are biological process, membrane and protein binding, there are 78, 198, 133 DEGs, respectively. In comparison of T2DM-8w and NG, and the biological processes, cellular component and molecular function enriched DEGs were the signal transduction, membrane and protein binding, there are 45, 176 and 120 DEGs, respectively. In the comparison of T2DM-12w and NG, and the biological processes, cellular component and molecular function enriched DEGs were the signal transduction, membrane and protein binding, there were 31, 129 and 91 DEGs, respectively. In summary, compared with NG, with the progression of diabetes, DEGs in T2DM group was mainly enriched in signal transduction, membrane and protein binding.

KEGG annotation classification and pathway enrichment

KEGG database was used for functional annotation classification and pathway enrichment of DEGs (Fig. 4). In the comparison of T2DM-4w and NG, KEGG pathway enrichment analysis results showed that the first three affected pathways included ECM-receptor interaction, cAMP signaling pathway, and PI3K-Akt signaling pathway showed higher enrichment and number of differential genes. In the comparison of T2DM-8w and NG, KEGG pathway enrichment analysis results showed that the first three affected pathways included ECM-receptor interaction, Fc epsilon RI signaling pathway, and Rap1 signaling pathway had higher enrichment degree and number of differential genes. In the comparison of T2DM-12w and NG, KEGG pathway enrichment analysis results showed that the first three affected pathways included B cell receptor signaling pathway, cAMP signaling pathway, and Chemokine signaling pathway had higher enrichment degree and number of differential genes.

Fig. 4
figure 4

The KEGG functional enrichment analysis in three comparison groups. (A) T2DM-4w group vs. NG: (B) T2DM-8w group vs. NG: (C) T2DM-12w group vs. NG

Transcriptomics alterations associated with the occur of diabetic tendinopathy

The result of the Venn diagram shows 103 key genes of sustained changes in the supraspinatus tendon following induction of diabetes, which are the first identified biomarkers of the supraspinatus tendinopathy as it occurs through the course of diabetes (Fig. 5A). The cluster heat map (Fig. 5B) clearly showed the distribution of DEGs at different time points after T2DM induction. The GO analysis results (Fig. 5C) 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 biological processes, cellular component and molecular function enriched DEGs were the signal transduction, membrane and protein binding, there were 10, 34 and 31 DEGs, respectively. 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 Fig. 5D; Table 2.

Fig. 5
figure 5

(A) The Venn diagram compares the number of different genes among all three comparison groups, including the paired comparison groups of T2DM-4w group vs. NG, T2DM-8w group vs. NG, and T2DM-12w group vs. NG. (B) Clustering heatmap of 103 DEGs. (C) Gene ontology (GO) enrichment analysis. (D) KEGG pathway analysis

Table 2 Pathways of major enrichment of 103 critical genes (p < 0.05)

Discussion

To date, transcriptome changes in injured tendons have been reported in several studies using bulk RNA sequencing [8,9,10]. In addition, single-cell and spatial RNA sequencing methods were used to further understand the cellular heterogeneity and molecular mechanism of tendinopathy progression, and the study results provided new insights into the control of tendinopathy [11]. Single-cell and spatial transcriptomics have shown that the pathogenesis of human tendinopathy is closely related to dysregulation of immune homeostasis [12]. At present, more and more studies are trying to discover the potential pathogenic small molecules of diabetes mellitus and its related complications through transcriptomics. Transcriptomics has a wide range of applications, including diabetic nephropathy, atherosclerosis, and retinopathy [22,23,24], which has resulted in the accumulation of large amount of biological data. However, the alterations in supraspinatus muscle tendon genes and the signaling pathways involved in their role in diabetic rats are still unclear.

This is the first study to establish a list of biomolecules that may be involved in the dynamic changes of the supraspinatus tendon of the rotator cuff with the progression of diabetes using transcriptomic techniques. In addition, we identified 103 key genes that may be closely related to the development of rotator cuff tendinopathy, and KEGG enrichment to 17 pathways, including the cAMP signaling pathway and Calcium signaling pathway. These results may provide a new target and theoretical basis for exploring the pathogenesis of diabetic tendinopathy.

The study contains a wealth of biological information that other researchers can use to explore the molecular changes associated with diabetes and tendinopathy. At the same time, KEGG analysis identified novel pathways that are differentially regulated in diabetic tendinopathy, including the cAMP signaling pathway and Calcium signaling pathway. Previous studies have shown that Cyclic adenosine monophosphate (cAMP) is an important intracellular second messenger, which plays various physiological roles through the activation of protein kinase A(PKA), and is a key regulatory factor in regulating metabolism and maintaining glucose homeostasis [25]. cAMP is known to play an important regulatory role in a variety of important physiological and pathological processes, including endothelial angiogenesis [26], osteoclast function [27], and reduction of oxidative stress [28]. cAMP signaling pathway also plays an important role in various complications caused by diabetes. cAMP signaling pathway is one of the important pathways in anti-diabetic nephropathy with renal fibrosis. Deb et al. [29] confirmed that activation of cAMP signaling pathway is an important mechanism to induce the transformation of renal tubular epithelial cells into stromal cells, thus promoting the progression of renal fibrosis. The activation of cAMP signaling pathway is closely related to increased glucose uptake in skeletal muscle. This physiological process is independent of insulin regulation, and the regulation of cAMP signaling pathway in skeletal muscle can be used as an alternative pathway for glucose processing and blood sugar control [30].

Of course, cAMP has a stimulating effect on insulin secretion and plays an active coordinating role in glucose-stimulated insulin secretion [31]. Hellman and Grill found that increased cAMP concentration plays an important role in glucose-stimulated insulin secretion in beta cells [32, 33]. Intracellular Ca2+ elevation is considered to be the ultimate activation mechanism of insulin secretion by beta cells. Later studies confirmed that glucose stimulation increases the level of cAMP, which is related to Ca2+ concentration [34, 35]. Our study also found that calcium signaling pathway may be associated with the development of diabetic tendinopathy, but whether there is a correlation between the two pathways requires further investigation of the mechanism. In addition, protecting islets by targeting cAMP signaling pathway is also an important strategy to promote glucose homeostasis. Because of the key role of cAMP signaling in regulating physiological and pathological processes, many drugs targeting cAMP are used in various diseases [36,37,38,39]. However, the correlation between cAMP signaling pathway and tendinopathy has not been reported. From the perspective of interfering with cAMP signaling pathway, it may be a new idea to actively develop prevention and treatment drugs to inhibit diabetic tendinopathy.

In addition, we found that the mechanism of tendinopathy induced by T2DM may be related to the PI3K-Akt signaling pathway and AMPK signaling pathway. Previous literature has also confirmed that the PI3K-Akt signaling pathway is closely related to impaired wound healing related to diabetes, and the up-regulation of the PI3K/Akt signaling pathway may be a promising target for accelerating healing [40]. AMPK and PI3K/Akt signaling pathways may be involved in diabetes-induced cognitive decline, and DL-3-n-butylphthalide ameliorates diabetes-associated cognitive decline by enhancing PI3K/Akt signaling [41].AMPK may be a fundamental cellular pathway for metabolic disorders and diabetes [42]. AMPK-mediated inflammation and oxidative stress are associated with decreased angiogenesis of endothelial progenitor cells associated with diabetes [43]. Meanwhile, AMPK and PI3K/Akt signaling pathways are closely related to mediating autophagy [44]. Diabetes is considered to be a mitochondrial disease at least to a certain extent, and because mitochondria play an important role in maintaining environmental homeostasis within rotator cuff tendons, mitochondrial damage may be associated with the occurrence of diabetic tendinopathy. However, whether autophagy mediated by AMPK and PI3K/Akt signaling pathways is involved in diabetes-induced rotator cuff tendinopathy is unclear. In the future, more work is needed to validate the role of AMPK and PI3K/Akt signaling pathways in diabetic tendinopathy.

In conclusion, the “gene profile” of supraspinatum tendon in rats was changed at different time points after T2DM induction, and the 103 DEGs identified may be potential molecular markers of diabetic tendinopathy. The cAMP signaling pathway and Calcium signaling pathway may be a potential pathogenic pathway for the occur of diabetic tendinopathy. Our work provides complete and detailed data support for exploring the relationship between diabetes and tendinopathy, which is one of the biggest innovations and contributions of this study. Our dataset provides a wealth of differentially expressed genes, each of which provides a basis for future research. The data could be useful to future researchers, giving researchers access to a huge subset of genes. Finally, these data provide raw data for identifying potential biological targets for future targeted therapies.

However, our study still has certain limitations. First of all, these 103 DEGs may play different roles in different stages of the pathogenesis of diabetic tendinopathy, so it is time-sensitive to some extent. However, our study did not carry out more in-depth verification and biological mechanism discussion based on these target genes. Second, due to the inevitable heterogeneity of experimental animals and the limited size of the sample, future studies need to replicate the experiment in parallel controls with larger sample cohorts. Finally, false positive results of some gene expression changes cannot be ruled out. After all, the transcriptomic sequencing method has certain limitations, and more accurate methods are needed to verify the results. At the same time, combining various platforms and multi-omics in further studies is crucial.

Conclusion

In this study, we established for the first time a biomolecule list of dynamic changes in rat supraspinatus tendon with the progression of diabetes using transcriptomic techniques. 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.

Data availability

The datasets generated and/or analyzed during the current study are not publicly available due the confidentiality of the participants’data, but are available from the corresponding author on reasonable requests.

Abbreviations

T2DM:

Type 2 diabetes mellitus

NG:

Normal group

DEGs:

Differentially expressed genes

GO:

gene ontology

KEGG:

Kyoto encyclopedia of genes and gnomes

STZ:

Streptozotocinn

PBS:

Phosphate-buffered saline

References

  1. Millar NL, Silbernagel KG, Thorborg K, et al. Tendinopathy Nat Rev Dis Primers. 2021;7(1):1.

    Article  PubMed  Google Scholar 

  2. Longo UG, Franceschi F, Ruzzini L, et al. Higher fasting plasma glucose levels within the normoglycaemic range and rotator cuff tears. Br J Sports Med. 2009;43:284–7.

    Article  CAS  PubMed  Google Scholar 

  3. Spoendlin J, Meier C, Jick SS, et al. Achilles or biceps tendon rupture in women and men with type 2 diabetes: a populationbased case-control study. Diabetes Complicat. 2016;30(5):903–9.

    Article  Google Scholar 

  4. Park HB, Gwark JY, Im JH, et al. Factors associated with atraumatic posterosuperior rotator cuff tears. J Bone Joint Surg Am. 2018;100(16):1397–405.

    Article  PubMed  Google Scholar 

  5. Nichols AEC, Oh I, Loiselle AE. Effects of type II diabetes mellitus on tendon homeostasis and healing. J Orthop Res. 2020;38:13–22.

    Article  PubMed  Google Scholar 

  6. Batista F, Nery C, Pinzur M, et al. Achilles tendinopathy in diabetes mellitus. Foot Ankle Int. 2008;29(5):498–501.

    Article  PubMed  Google Scholar 

  7. Abate M, Schiavone C, Salini V. Sonographic evaluation of the shoulder in asymptomatic elderly subjects with diabetes[J]. BMC Musculoskelet Disord. 2010;11:278.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tashjian RZ, Lock I, Granger EK, et al. Gene expression in torn rotator cuff tendons determined by RNA sequencing. Orthop J Sports Med. 2020;8(6):2325967120927480.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Cho N, Lee SG, Kim JO, et al. Identification of differentially expressed genes associated with extracellular matrix degradation and inflammatory regulation in calcific tendinopathy using RNA sequencing. Calcif Tissue Int. 2020;107(5):489–98.

    Article  CAS  PubMed  Google Scholar 

  10. Jelinsky SA, Rodeo SA, Li J, et al. Regulation of gene expression in human tendinopathy. BMC Musculoskelet Disord. 2011;12:86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Fu W, Yang R, Li J. Single-cell and spatial transcriptomics reveal changes in cell heterogeneity during progression of human tendinopathy. BMC Biol. 2023;21(1):132.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Akbar M, MacDonald L, Crowe LAN, et al. Single cell and spatial transcriptomics in human tendon disease indicate dysregulated immune homeostasis. Ann Rheum Dis. 2021;80(11):1494–7.

    Article  CAS  PubMed  Google Scholar 

  13. Xu K, Zhang L, Ren Z, et al. Evaluating the role of type 2 diabetes mellitus in rotator cuff tendinopathy: development and analysis of a novel rat model. Front Endocrinol (Lausanne). 2022;13:1042878.

    Article  PubMed  Google Scholar 

  14. Liu Y, Zhang X, Yang L, et al. Proteomics and transcriptomics explore the effect of mixture of herbal extract on diabetic wound healing process. Phytomedicine. 2023;116:154892.

    Article  CAS  PubMed  Google Scholar 

  15. Guo NL, Poh TY, Pirela S, et al. Integrated Transcriptomics, Metabolomics, and Lipidomics Profiling in Rat Lung, blood, and serum for Assessment of Laser Printer-emitted nanoparticle inhalation exposure-Induced Disease risks. Int J Mol Sci. 2019;20(24):6348.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Mao Y, Sun J, Wang Z, et al. Combining transcriptomic analysis and network pharmacology to explore the mechanism by which Shaofu Zhuyu decoction improves diabetes mellitus erectile dysfunction. Phytomedicine. 2023;119:155006.

    Article  CAS  PubMed  Google Scholar 

  17. Dumas ME, Domange C, Calderari S, et al. Topological analysis of metabolic networks integrating co-segregating transcriptomes and metabolomes in type 2 diabetic rat congenic series. Genome Med. 2016;8(1):101.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Xu K, Zhang L, Wang T, et al. Untargeted metabolomics reveals dynamic changes in metabolic profiles of rat supraspinatus tendon at three different time points after diabetes induction. Front Endocrinol (Lausanne). 2023;14:1292103.

    Article  PubMed  Google Scholar 

  19. Srinivasan K, Viswanad B, Asrat L, et al. Combination of high-fat diet-fed and low-dose streptozotocin-treated rat: a model for type 2 diabetes and pharmacological screening. Pharmacol Res. 2005;52(4):313–20.

    Article  CAS  PubMed  Google Scholar 

  20. Hong W, Li S, Cai Y, et al. The Target MicroRNAs and potential underlying mechanisms of Yiqi-Bushen-Tiaozhi recipe against-non-alcoholic steatohepatitis. Front Pharmacol. 2020;11:529553.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhou ZY, Zhao WR, Xiao Y, et al. Antiangiogenesis effect of timosaponin AIII on HUVECs in vitro and zebrafish embryos in vivo. Acta Pharmacol Sin. 2020;41(2):260–9.

    Article  CAS  PubMed  Google Scholar 

  22. Yu T, Xu B, Bao M, et al. Identification of potential biomarkers and pathways associated with carotid atherosclerotic plaques in type 2 diabetes mellitus: a transcriptomics study. Front Endocrinol (Lausanne). 2022;13:981100.

    Article  PubMed  Google Scholar 

  23. Wang R, Jian Q, Hu G, et al. Integrated Metabolomics and Transcriptomics reveal metabolic patterns in retina of STZ-Induced Diabetic Retinopathy Mouse Model. Metabolites. 2022;12(12):1245.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Rangasamy S, Monickaraj F, Legendre C, et al. Transcriptomics analysis of pericytes from retinas of diabetic animals reveals novel genes and molecular pathways relevant to blood-retinal barrier alterations in diabetic retinopathy. Exp Eye Res. 2020;195:108043.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang Y, Liu Q, Kang SG, et al. Dietary bioactive ingredients modulating the cAMP signaling in Diabetes Treatment. Nutrients. 2021;13(9):3038.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kim JJ, Kim DH, Lee JY, et al. cAMP/EPAC signaling enables ETV2 to induce endothelial cells with high angiogenesis potential. Mol Ther. 2020;28(2):466–78.

    Article  CAS  PubMed  Google Scholar 

  27. Jeevaratnam K, Salvage SC, Li M, et al. Regulatory actions of 3’,5’-cyclic adenosine monophosphate on osteoclast function: possible roles of Epac-mediated signaling. Ann N Y Acad Sci. 2018;1433(1):18–28.

    Article  CAS  PubMed  Google Scholar 

  28. Stokman G, Qin Y, Booij TH, et al. Epac-Rap signaling reduces oxidative stress in the tubular epithelium. J Am Soc Nephrol. 2014;25(7):1474–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Deb DK, Bao R, Li YC. Critical role of the cAMP-PKA pathway in hyperglycemia-induced epigenetic activation of fibrogenic program in the kidney. FASEB J. 2017;31(5):2065–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sato M, Dehvari N, Oberg AI, et al. Improving type 2 diabetes through a distinct adrenergic signaling pathway involving mTORC2 that mediates glucose uptake in skeletal muscle. Diabetes. 2014;63(12):4115–29.

    Article  CAS  PubMed  Google Scholar 

  31. Hagren OI, Tengholm A. Glucose and insulin synergistically activate phosphatidylinositol 3-kinase to trigger oscillations of phosphatidylinositol 3,4,5-trisphosphate in beta-cells. J Biol Chem. 2006;281(51):39121–7.

    Article  CAS  PubMed  Google Scholar 

  32. Grill V, Cerasi E. Activation by glucose of adenyl cyclase in pancreatic islets of the rat. FEBS Lett. 1973;33(3):311–4.

    Article  CAS  PubMed  Google Scholar 

  33. Hellman B, Idahl LA, Lernmark A, et al. The pancreatic beta-cell recognition of insulin secretagogues: does cyclic AMP mediate the effect of glucose? Proc Natl Acad Sci U S A. 1974;71(9):3405–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Nakagawa Y, Nagasawa M, Medina J, et al. Glucose evokes Rapid Ca2 + and cyclic AMP signals by activating the cell-surface glucose-sensing receptor in pancreatic β-Cells. PLoS ONE. 2015;10(12):e0144053.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tian G, Sol ER, Xu Y, et al. Impaired cAMP generation contributes to defective glucose-stimulated insulin secretion after long-term exposure to palmitate. Diabetes. 2015;64(3):904–15.

    Article  CAS  PubMed  Google Scholar 

  36. Pierre S, Eschenhagen T, Geisslinger G, et al. Capturing adenylyl cyclases as potential drug targets. Nat Rev Drug Discov. 2009;8(4):321–35.

    Article  CAS  PubMed  Google Scholar 

  37. Tong T, Ryu SE, Min Y, et al. Olfactory receptor 10J5 responding to α-cedrene regulates hepatic steatosis via the cAMP-PKA pathway. Sci Rep. 2017;7(1):9471.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Tong T, Park J, Moon C, et al. Regulation of Adipogenesis and Thermogenesis through mouse olfactory receptor 23 stimulated by α-Cedrene in 3T3-L1 cells. Nutrients. 2018;10(11):1781.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Tong T, Shen Y, Lee HW, et al. Adenylyl cyclase 3 haploinsufficiency confers susceptibility to diet-induced obesity and insulin resistance in mice. Sci Rep. 2016;6:34179.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Qu KS, Ru Y, Yang D, et al. Fu-Huang ointment ameliorates impaired wound healing associated with diabetes through PI3K-AKT signalling pathway activation. Comput Biol Med. 2023;155:106660.

    Article  CAS  PubMed  Google Scholar 

  41. Wang BN, Wu CB, Chen ZM, et al. DL-3-n-butylphthalide ameliorates diabetes-associated cognitive decline by enhancing PI3K/Akt signaling and suppressing oxidative stress. Acta Pharmacol Sin. 2021;42(3):347–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Maiese K. Cornerstone Cellular pathways for metabolic disorders and Diabetes Mellitus: non-coding RNAs, wnt signaling, and AMPK. Cells. 2023;12(22):2595.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Luo L, Dong B, Zhang J, et al. Dapagliflozin restores diabetes-associated decline in vasculogenic capacity of endothelial progenitor cells via activating AMPK-mediated inhibition of inflammation and oxidative stress. Biochem Biophys Res Commun. 2023;671:205–14.

    Article  CAS  PubMed  Google Scholar 

  44. Zheng X, Li W, Xu H, et al. Sinomenine ester derivative inhibits glioblastoma by inducing mitochondria-dependent apoptosis and autophagy by PI3K/AKT/mTOR and AMPK/mTOR pathway. Acta Pharm Sin B. 2021;11(11):3465–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

KSX and YZZ conceived the study design. TRW supervised the data collection. KSX and LZ drafted the manuscript. TBY and XZ contributed to the revision. YZZ is responsible for this article. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Xia Zhao or Yingze Zhang.

Ethics declarations

Ethics approval and consent to participate

The animal study was reviewed and approved by Ethics Committee for Laboratory Animal Welfare of Qingdao University(No. 20220505SD8020221210126).

Consent for publication

Not applicable.

Competing interests

The authors declare that there is no competing of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, K., Zhang, L., Wang, T. et al. Transcriptomics reveals dynamic changes in the “gene profiles” of rat supraspinatus tendon at three different time points after diabetes induction. BMC Med Genomics 17, 122 (2024). https://doi.org/10.1186/s12920-024-01899-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-01899-3

Keywords