Skip to main content
Download PDF

Screening for autophagy/hypoxia/ferroptosis/pyroptosis-related genes of tendon injury and repair in a rat model after celecoxib and lactoferrin treatment

Journal of Orthopaedic Surgery and Research volume 18, Article number: 383 (2023) Cite this article

Abstract

Background

Tendon injuries are among the most common musculoskeletal disorders. Celecoxib possesses an effective anti-inflammatory activity in the tendon injury treatment. Lactoferrin has a great potential for the tendon regeneration. However, the efficacy of celecoxib combined with lactoferrin in the treatment of tendon injury has not been reported. In this study, we aimed to investigate the effect of celecoxib and lactoferrin on tendon injury and repair, and screen for the crucial genes associated with the tendon injury and repair.

Methods

The rat tendon injury models were established and divided into four groups: normal control group (n = 10), tendon injury model group (n = 10), celecoxib treatment group (n = 10), and celecoxib + lactoferrin treatment group (n = 10). Then, RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs) and mRNAs (DEmRNAs) in celecoxib treatment group and celecoxib + lactoferrin treatment group. Next, autophagy/hypoxia/ferroptosis/pyroptosis-related DEmRNAs were further identified. Subsequently, functional enrichment, protein–protein interaction (PPI) network and transcriptional regulatory network construction for these genes were performed.

Results

The animal study demonstrated that combinational administration of celecoxib with lactoferrin rescued the harmful effects caused by celecoxib in the treatment of tendon injury. Compared to tendon injury model group, 945 DEmRNAs, 7 DEmiRNAs and 34 DElncRNAs were obtained in celecoxib treatment group, and 493 DEmRNAs, 8 DEmiRNAs and 21 DElncRNAs were obtained in celecoxib + lactoferrin treatment group, respectively. Subsequently, 376 celecoxib + lactoferrin treatment group-specific DEmRNAs were determined. Then, 25 DEmRNAs associated with autophagy/hypoxia/ferroptosis/pyroptosis were identified.

Conclusions

Several genes, such as, Ppp1r15a, Ddit4, Fos, Casp3, Tgfb3, Hspb1 and Hspa8, were identified to be associated with tendon injury and repair.

Introduction

Musculoskeletal disorders exert the significant detrimental effect on the life quality of patients [1]. Multiple pathological factors are associated with tendon injury, including trauma, aging, inflammation, chronic overuse, and genetic factors. The mechanisms of injury include some types of overload or overuse of the tendon that most likely lead to tendon degeneration, cell phenotype changes and hypervascularization [2]. In the clinical work, nonsteroidal anti-inflammatory drugs (NSAIDs) are used to treat tendon injury, such as, selective COX-2 inhibitor and celecoxib. Celecoxib possesses an effective analgesic and anti-inflammatory activity in the tendon injury treatment [3]. The side effects of celecoxib include inhibition of tendon cell proliferation and migration, which can adversely affect tendon healing [4]. Lactoferrin, an iron-binding glycoprotein, functions as an osteogenic growth factor and enhance the proliferation and differentiation of osteoblasts, and has a great potential for the tendon regeneration [5]. Moreover, when lactoferrin is used in combination with non-steroidal anti-inflammatory drugs, the inhibitory effects on tenocyte proliferation, viability, and collagen formation are rescued [6]. We have previously shown that lactoferrin is anabolic to human tenocytes in vitro and reverses potential inhibitory effects of NSAIDs on human tenocytes [7]. To date, the efficacy of celecoxib combined with lactoferrin in the treatment of tendon injury has not been reported.

Tendon fibroblasts play an important role in remodeling phase of wound healing [8]. In tendon injury, high glucose can repress the proliferation of tendon fibroblasts by inhibiting autophagy activation [9]. It has been indicated that hypoxia has a critical function in chondrogenesis, osteogenesis and angiogenesis, and plays an essential role in the tissue repair process [10]. Zhao et al. reported that hypoxia was essential for bone–tendon junction healing [11]. Yu et al. suggested that hypoxia enhanced tenocyte differentiation of adipose-derived mesenchymal stem cells [12]. Chen et al. demonstrated that hypoxia-induced mesenchymal stem cells exhibited stronger tenogenic differentiation capacities [13]. Ferroptosis, a new form of regulated cell death (driven by iron-dependent lipid peroxidation), is involved in a variety of diseases [14]. Iron is a key element that plays a crucial role in mammalian cells (such as osteoclast) [15]. Ni et al. demonstrated that ferroptosis was involved in osteoclastogenesis [16]. Pyroptosis (also known as cell inflammatory necrosis) is a programmed cell death mode closely associated with the inflammatory response and mediated by caspase-1 or caspase-11 [17, 18]. In alveolar bone, high glucose concentration may activate pyroptosis to inhibit the proliferation and differentiation of osteoblasts [19]. Under certain pathological conditions, pyroptosis may occur in osteoblasts, affect their proliferation and differentiation and consequently affect the development and morphological changes of bone tissue [20].

In this study, we performed animal studies to investigate the efficacy of celecoxib combined with lactoferrin in the treatment of tendon injury. Then, differentially expressed mRNAs (DEmRNAs), differentially expressed miRNAs (DEmiRNAs) and differentially expressed lncRNAs (DElncRNAs) in celecoxib treatment group and celecoxib + lactoferrin treatment group were identified. As we mentioned above, autophagy/hypoxia/ferroptosis/pyroptosis may be associated with tendon injury and repair. Hence, autophagy/hypoxia/ferroptosis/pyroptosis-related gene sets were downloaded for subsequent analysis. In this study, we aimed to investigate the effect of celecoxib and lactoferrin on tendon injury and repair, and screen for the crucial genes associated with the tendon injury and repair.

Materials and methods

Animal model

A total of 40 six-week-old Sprague–Dawley rats weighting 200–300 g were used. The rats were housed under a 12-h light/dark cycle in a pathogen-free area with free access to water and food. All animals were treated according to institutional guidelines for laboratory animal treatment and care. All experimental procedures were approved by the Animal Research Ethics Committee of our hospital.

Rats were randomly divided into four groups: normal control group (n = 10), tendon injury model group (n = 10), celecoxib treatment group (n = 10), and celecoxib + lactoferrin treatment group (n = 10). Lactoferrin was purchased from Wuhan Nuohui Pharmaceutical & Chemical Co., LTD. Following 1 week of feeding and adaptation, the rats were anesthetized. The midpoint of the achilles tendon was transected and sutured immediately. The tendon was not cut in the normal control group, and other operations were the same as those in the experimental group. Celecoxib (10 mg/kg) were injected daily at the site of injury in each rat in celecoxib treatment group. Celecoxib (10 mg/kg) and lactoferrin (2 g/kg) were injected daily at the site of injury in each rat in celecoxib and lactoferrin treatment group.

Hematoxylin and eosin (HE) staining, Masson staining and immunohistochemistry

Rats were sacrificed on the 1st, 14th, and 28th day after surgery. Three rats in each group were sacrificed on the 1st and 14th day, and four rats were sacrificed on the 28th day. Then, the tendon of the injured site was collected. For histological analysis, tendons were fixed in 4% paraformaldehyde for 24 h. Sagittal paraffin sections were prepared by embedding the samples in paraffin and cutting into 4 µm thick sections. Four sections per sample were stained with HE and observed under the OLYMPUS EX-51 microscope (Tokyo, Japan) at 100 × magnification. Masson’s trichrome stain (Solarbio, G1340) was performed according to kit directions. RECA-1 is a cell surface antigen expressed by rat endothelial cells. In this study, RECA-1 was selected for immunohistochemistry analysis on days 1, 14 and 28. Then, 5-μm-thick continuous sections were incubated with an anti-RECA-1 antibody (abcam, ab22492, 1:1000) produced in rabbit followed by goat anti-rabbit immunoglobulin antibody conjugated by horseradish peroxidase. Then, slides were visualized using diaminobenzidine (DAB) substrate.

RNA isolation and sequencing

We selected samples from injury model group, celecoxib treatment group and celecoxib + lactoferrin treatment group on day 14 for subsequent RNA sequencing. Total RNA was harvested using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) from the tendons. With Agilent 2100 and Nanodrop2000, the quality of RNA was assessed. The quality of the libraries was determined using an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System. Illumina Hiseq x-ten platform was used to perform RNA sequencing for mRNA and lncRNA. The raw sequencing data were submitted to sequencing quality control by FastQC. Reads with low quality were removed. BGISEQ-500 platform was used for RNA sequencing for miRNA. The Fastx-Toolkit was used to trim 5’ and 3’ segments of reads to remove bases (with mass < 20 and delete reads with N > 10%). The Rfam was used for annotation analysis on measured small RNA with BLAST v2.3.0. The mature miRNA and miRNA precursor sequences were downloaded from miRBase. The expression of miRNA was quantified with miRDeep2. HISAT2 was used to align the clean reads with the reference genome, Ensemble Rnor_6.0.

Differential expression analysis

By using DESeq2, the DEmRNAs and DElncRNAs were identified with p-value < 0.05. With DEGseq2, DEmiRNAs were identified with p-value < 0.05 as well. With R package “pheatmap,” hierarchical clustering analysis was performed. David 6.8 was applied to perform GO and KEGG enrichment analysis for DEmRNAs with p-value < 0.05. Particularly, celecoxib + lactoferrin treatment group-specific DEmRNAs were further obtained.

Identification of genes associated with autophagy/hypoxia/ferroptosis/pyroptosis

Autophagy-related genes were extracted from Human Autophagy Database (HADb, http://www.autophagy.lu/index.html) and the GOBP REGULATION OF AUTOPHAGY gene set in Molecular Signatures Database (MSigDB). Totally, 516 autophagy-related genes were included for subsequent analysis. Then, 200 hypoxia-related genes were retrieved from the MSigDB. A total of 267 ferroptosis-related genes were retrieved from the FerrDb dataset (http://www.zhounan.org/ferrdb/) and the previous literature [21]. In addition, 41 pyroptosis-related genes were retrieved from the previous literature [22,23,24]. Finally, genes associated with autophagy/hypoxia/ferroptosis/pyroptosis were obtained by overlapping celecoxib + lactoferrin treatment group-specific DEmRNAs with autophagy/hypoxia/ferroptosis/pyroptosis-related genes, respectively.

Functional annotation and protein–protein interaction (PPI) network construction

In order to explore the biological functions and the potential pathways of the genes associated with autophagy/hypoxia/ferroptosis/pyroptosis, David 6.8 was utilized to perform GO and KEGG enrichment analysis. A p-value < 0.05 was considered statistically significant. Online database STRING (https://string-db.org) was used to analyze the PPI networks.

Construction of transcriptional regulatory networks

Potential transcription factors (TFs) targeted to the genes associated with autophagy/hypoxia/ferroptosis/pyroptosis were identified using the transcriptional regulatory relationships unraveled database. The TF-gene regulatory networks were visualized by using Cytoscape software.

Results

Celecoxib combined with lactoferrin treatment decreased inflammatory cell and the degree of fiber structure disorder

In Fig. 1A, HE staining confirmed that the injured tendon had more inflammatory cell nuclei, obvious inflammation infiltration, and disorder of fiber structure in injury model group compared with the other three groups. The number of inflammatory cell nuclei decreased in the celecoxib treatment group and celecoxib + lactoferrin treatment group compared with the injury model group. The degree of fiber structure disorder also decreased in the celecoxib treatment group and celecoxib + lactoferrin group compared with the injury model group. In Fig. 1B and 1D, Masson’s trichrome staining showed that the injury model group had fewer collagen fibers and the collagen fiber structure was obviously disordered compared with the normal control group. Compared with the celecoxib treatment group, celecoxib + lactoferrin group formed new collagen fibers and increased the proportion of positive staining area. In Fig. 1C and 1F, compared with the injury model group, the distribution of RECA-1 protein in the celecoxib + lactoferrin group significantly increased. The distribution range of RECA-1 protein in the celecoxib treatment group and celecoxib + lactoferrin group increased, and the color became darker, indicating that the inflammation recovery period led to the proliferation of new capillaries and increased blood vessels. According to the scope and color depth, the celecoxib + lactoferrin treatment group has a wider range and darker color than the celecoxib treatment group, indicating that the inflammatory recovery is better in celecoxib + lactoferrin treatment group compared with the celecoxib treatment group.

Fig. 1
figure 1

HE staining (A), Masson staining (B) and immunohistochemistry (C) of tendons in normal control group, injury model group, celecoxib treatment group, celecoxib + lactoferrin treatment group on day 14. Scale bar = 50 μm. D Quantification analysis of Masson staining. E Quantification analysis of immunohistochemistry. * indicates p < 0.05, ** indicates p < 0.01

Full size image

Screening of DEmRNAs, DEmiRNAs and DElncRNAs

Compared to tendon injury model group, 945 (548 up-regulated and 397 down-regulated) DEmRNAs, 7 (5 up-regulated and 2 down-regulated) DEmiRNAs and 34 (20 up-regulated and 14 down-regulated) DElncRNAs were obtained in celecoxib treatment group, and 493 (219 up-regulated and 274 down-regulated) DEmRNAs, 8 (3 up-regulated and 5 down-regulated) DEmiRNAs and 21 (8 up-regulated and 13 down-regulated) DElncRNAs were obtained in celecoxib + lactoferrin treatment group, respectively. The heatmap of DEmRNAs, DEmiRNAs and DElncRNAs is shown in Figs. 2 and 3. GO enrichment analysis revealed that these biological processes such as blood vessel development, bone development, muscle cell differentiation, and muscle tissue development were dysregulated in celecoxib treatment group (Additional file 1: Figure S1A). KEGG pathway analysis highlighted that MAPK signaling pathway, ECM-receptor interaction, T cell receptor signaling pathway and GnRH signaling pathway were dysregulated in celecoxib treatment group (Additional file 1: Figure S1B). GO enrichment analysis revealed that these biological processes such as skeletal muscle cell differentiation, response to hypoxia, positive regulation of smooth muscle cell proliferation and negative regulation of cell proliferation were dysregulated in celecoxib + lactoferrin treatment group (Additional file 2: Figure S2A). KEGG pathway analysis highlighted that MAPK signaling pathway, HIF-1 signaling pathway, cAMP signaling pathway and PI3K-Akt signaling pathway were dysregulated in celecoxib + lactoferrin treatment group (Additional file 2: Figure S2B). Due to the small number of DEmiRNAs and DElncRNAs, our next research focus is mainly on DEmRNAs. Subsequently, 376 (169 up-regulated and 207 down-regulated) celecoxib + lactoferrin treatment group-specific DEmRNAs were determined (Additional file 3: Figure S3). Also, the ceRNA network was built, but no results were available.

Fig. 2
figure 2

The heatmap of DEmRNAs (A), DEmiRNAs (B) and DElncRNAs (C) between injury model group and celecoxib treatment group on day 14

Full size image
Fig. 3
figure 3

The heatmap of DEmRNAs (A), DEmiRNAs (B) and DElncRNAs (C) between injury model group and celecoxib + lactoferrin treatment group on day 14

Full size image

Identification of genes associated with autophagy/hypoxia/ferroptosis/pyroptosis

In order to further investigate what role the celecoxib + lactoferrin treatment group-specific DEmRNAs play in tendon injury repair, the autophagy/hypoxia/ferroptosis/pyroptosis-related gene sets were retrieved from the corresponding databases and related literature. In total, 25 genes associated with autophagy/hypoxia/ferroptosis/pyroptosis were obtained by overlapping celecoxib + lactoferrin treatment group-specific DEmRNAs with autophagy/hypoxia/ferroptosis/pyroptosis-related genes, respectively (Table 1). In addition, 6 genes (Hspb1, Fos, Gapdh, Ppp1r15a, Casp3, and Ddit4) were associated with at least two of autophagy, hypoxia, ferroptosis, and pyroptosis.

Table 1 The list of autophagy/hypoxia/ferroptosis/pyroptosis-related genes
Full size table

Functional annotation and PPI network construction

GO enrichment analysis revealed that these 25 genes were significantly enriched in negative regulation of apoptotic process, aging and response to hypoxia (Fig. 4A). KEGG pathway analysis highlighted that MAPK signaling pathway, HIF-1 signaling pathway and TNF signaling pathway were dysregulated (Fig. 4B). The PPI network included 25 proteins, which consisted of two subgroups: 17 proteins with strong connections with others and 8 separated proteins (Fig. 5). Among them, Hspa8 and Hspb1 had the highest interaction score, which was 0.951.

Fig. 4
figure 4

GO (A) and KEGG (B) enrichment analysis of the genes associated with autophagy/hypoxia/ferroptosis/pyroptosis BP: biological process; CC: cytological component; MF: molecular function

Full size image
Fig. 5
figure 5

PPI network each node represents a protein, while each edge represents one protein–protein association

Full size image

Construction of TF-gene regulatory network

A total of 12 genes targeted by TFs were identified, including Casp3, Cited2, Cryab, Ddit4, Eif4e, Fos, Hspb1, Itgb4, Noct, Slc6a6, Slc40a1, and Tgfb3 (Fig. 6). Among which, Fos, Casp3, Slc6a6, Ddit4 can be regulated by Trp53. It is noted that Fos can be regulated by itself.

Fig. 6
figure 6

TF regulatory networks Red and green circles represent up- and down-regulated DEmRNAs. Rhombuses represent transcription factors (TFs). Edges indicate TF-DEmRNAs interactions

Full size image

Discussion

In this study, the rat tendon injury models were established and divided into four groups: normal control, tendon injury model, celecoxib treatment and celecoxib + lactoferrin treatment groups. Based on the results of HE staining, Masson staining and immunohistochemistry, we speculated that combinational administration of celecoxib with lactoferrin could rescue the harmful effects caused by celecoxib in the treatment of tendon injury. Then, the RNA sequencing and bioinformatics analysis indicated that compared to tendon injury model group, 945 DEmRNAs, 7 DEmiRNAs and 34 DElncRNAs were obtained in celecoxib treatment group, and 493 DEmRNAs, 8 DEmiRNAs and 21 DElncRNAs were obtained in celecoxib + lactoferrin treatment group. Given the small number DEmiRNAs and DElncRNAs obtained, we focused on DEmRNAs in the subsequent analysis. Then, 376 celecoxib + lactoferrin treatment group-specific DEmRNAs were determined. Next, 25 DEmRNAs associated with autophagy/hypoxia/ferroptosis/pyroptosis were obtained. Finally, several genes, such as, Ppp1r15a, Ddit4, Fos, Casp3, Tgfb3, Hspb1 and Hspa8, were identified to be associated with tendon injury and repair.

Aberrant mammalian target of rapamycin (mTOR) complex 1 (mTORC1) signaling is associated with altered bone homeostasis [25]. Protein phosphatase 1 regulatory subunit 15A (PPP1R15A), also known as GADD34, is a positive regulator of osteoclastogenesis and suppresses mTORC1 activity at the later stages of osteoclastogenesis [26]. DNA-damage-inducible transcript 4 (DDIT4) is an inhibitor of mTOR signaling [27]. DDIT4 confers a protective effect on radiation-induced premature senescence in osteoblast cells [28]. The JUN class of transcription factors is composed of heterdimers of Fos-related factors and Jun proteins, which is an important transcription factor in osteoblastic differentiation [29]. The c-Fos and c-Jun genes are the commonly studied member of the cellular immediate-early genes [30]. FOS is reported to improve tendon healing by promoting tendon cell proliferation and differentiation and regulating the inflammatory response [31]. A previous study suggested that ultrasound increased bone morphogenetic protein-2 expression in osteoblasts via the PI3K, Akt, c-Fos/c-Jun, and AP-1 signaling pathways [32].

Heat shock proteins are classified based on their molecular weights and include small HSPs, HSP40, HSP60, HSP70, HSP90, and large HSPs (HSP110 and glucoseregulated protein 170, GRP170) [33]. Heat shock protein family B (small) member 1 (HSPB1), also known as HSP27, is a member of the small heat shock protein family, which participates in the regulation of multiple physiological and pathophysiological cell functions. Elevated HSPB1 was related to estrogen-induced resistance to osteoblast apoptosis [34]. Unphosphorylated HSPB1 suppresses fibroblast growth factor-2-stimulated vascular endothelial growth factor release in osteoblasts [35]. HSPB1 suppressed platelet-derived growth factor-BB induced cell migration of osteoblasts [36]. Phosphorylated HSPB1 was involved in the pathogenesis of osteoporosis [37]. Heat shock protein family A (Hsp70) member 8 (HSPA8) is a member of the heat shock protein 70 family. The association of MNSFβ with HSPA8 may promote RANKL-induced osteoclastogenesis [38].

Caspase-3, one of the major activated cysteine proteases, constitutes the caspase family. You et al. investigated the effect of CASP3 inhibition on osteoblast differentiation capacity in high glucose conditions in mouse osteoblastic cell line MC3T3-E1 [39]. Mogi suggested that caspase activity could be required for osteogenic differentiation of osteoblastic cell [40]. Slawomir and his colleagues reported that CASP3, a pro-inflammatory factor involved in the TNF-α transduction pathway, was expressed in the injured tendon tissues of patients with rotator cuff tendinopathy [41]. CASP3 is suggested to be necessary for extracellular matrix remodeling in rat patellar tendon [42]. The transforming growth factor-β (TGFβ) family is expressed in various cell types and plays a crucial role in cellular division, migration, adhesion, differentiation, and programmed death [43]. The vascular smooth muscle cells proliferation pathways are known to be associated with various mRNAs and non-coding RNAs [44,45,46]. A certain level of active TGF-beta suppressed MMP-2 expression to promote the contractile phenotype of vascular smooth muscle cells [47]. TGF-β3 effectively protected against flexor tendon injury via regulating adhesion formation through the JNK/c-Jun pathway [48].

Conclusions

In the present study, multiple bioinformatics analysis methods were applied to identify genes that play a key role in tendon injury and repair. Hspb1 was identified as an autophagy/ferroptosis-related gene, Fos and Ppp1r15a were identified as autophagy/hypoxia-related genes, Casp3 was identified as an autophagy/pyroptosis-related gene, Hspa8 was identified as an autophagy-related gene, Tgfb3 was identified as a hypoxia-related gene, and Ddit4 was identified as a hypoxia/ferroptosis-related gene. Hspa8 and Hspb1 had the highest interaction score in PPI network, which was 0.951. In addition, transcriptional regulatory network indicated that Fos can be regulated by itself. The above analysis results demonstrate the important role of these genes in tendon injury and repair from different perspectives. Also, the mechanism remained largely unknown. Our further research endeavors will concentrate on elucidating the fundamental mechanisms underlying the effect of celecoxib and lactoferrin on tendon injury and repair, including the intricate signaling pathways and protein expression profiles involved.

Availability of data and materials

The RNA-seq data have been deposited in the GSE201161 dataset (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE201161).

Abbreviations

DAB:

Diaminobenzidine

DDIT4:

DNA-damage-inducible transcript 4

DElncRNAs:

Differentially expressed lncRNAs

DEmiRNAs:

Differentially expressed miRNAs

DEmRNAs:

Differentially expressed mRNAs

DUSP1:

Dual-specificity phosphatase 1

HADb:

Human autophagy database

HSPA5:

Heat shock protein family A (Hsp70) member 5

HE:

Hematoxylin and eosin

HSPA8:

Heat shock protein family A (Hsp70) member 8

HSPB1:

Heat shock protein family B (small) member 1

MSigDB:

Molecular signatures database

mTOR:

Mammalian target of rapamycin

mTORC1:

Mammalian target of rapamycin complex 1

NSAIDs:

Clinical work, nonsteroidal anti-inflammatory drugs

PPI:

Protein–protein interaction

PPP1R15A:

Protein phosphatase 1 regulatory subunit 15A

TFs:

Transcription factors

TGFβ:

Transforming growth factor-β

References

  1. Liu Q, Zhu Y, Zhu W, Zhang G, Yang YP, Zhao C. The role of MicroRNAs in tendon injury, repair, and related tissue engineering. Biomaterials. 2021;277:121083. https://doi.org/10.1016/j.biomaterials.2021.121083.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Skovgaard D, Siersma VD, Klausen SB, Visnes H, Haukenes I, Bang CW, et al. Chronic hyperglycemia, hypercholesterolemia, and metabolic syndrome are associated with risk of tendon injury. Scand J Med Sci Sports. 2021;31(9):1822–31. https://doi.org/10.1111/sms.13984.

    Article  PubMed  Google Scholar 

  3. Wang J, Wang L, Gao Y, Zhang Z, Huang X, Han T, et al. Synergistic therapy of celecoxib-loaded magnetism-responsive hydrogel for tendon tissue injuries. Front Bioeng Biotechnol. 2020;8:592068. https://doi.org/10.3389/fbioe.2020.592068.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Tsai WC, Hsu CC, Chou SW, Chung CY, Chen J, Pang JH. Effects of celecoxib on migration, proliferation and collagen expression of tendon cells. Connect Tissue Res. 2007;48(1):46–51. https://doi.org/10.1080/03008200601071295.

    Article  CAS  PubMed  Google Scholar 

  5. Choi HJ, Choi S, Kim JG, Song MH, Shim KS, Lim YM, et al. Enhanced tendon restoration effects of anti-inflammatory, lactoferrin-immobilized, heparin-polymeric nanoparticles in an Achilles tendinitis rat model. Carbohydr Polym. 2020;241:116284. https://doi.org/10.1016/j.carbpol.2020.116284.

    Article  CAS  PubMed  Google Scholar 

  6. Antoshin AA, Shpichka AI, Huang G, Chen K, Lu P, Svistunov AA, et al. Lactoferrin as a regenerative agent: The old-new panacea? Pharmacol Res. 2021;167:105564. https://doi.org/10.1016/j.phrs.2021.105564.

    Article  CAS  PubMed  Google Scholar 

  7. Zhang Y, Wang X, Qiu Y, Cornish J, Carr AJ, Xia Z. Effect of indomethacin and lactoferrin on human tenocyte proliferation and collagen formation in vitro. Biochem Biophys Res Commun. 2014;454(2):301–7. https://doi.org/10.1016/j.bbrc.2014.10.061.

    Article  CAS  PubMed  Google Scholar 

  8. Thomopoulos S, Parks WC, Rifkin DB, Derwin KA. Mechanisms of tendon injury and repair. J Orthop Res Off Publ Orthop Res Soc. 2015;33(6):832–9. https://doi.org/10.1002/jor.22806.

    Article  Google Scholar 

  9. Song FC, Yuan JQ, Zhu MD, Li Q, Liu SH, Zhang L, et al. High glucose represses the proliferation of tendon fibroblasts by inhibiting autophagy activation in tendon injury. Biosci Rep. 2021;42(3):20210640. https://doi.org/10.1042/bsr20210640.

    Article  CAS  Google Scholar 

  10. Sun X, Wei Y. The role of hypoxia-inducible factor in osteogenesis and chondrogenesis. Cytotherapy. 2009;11(3):261–7. https://doi.org/10.1080/14653240902824765.

    Article  CAS  PubMed  Google Scholar 

  11. Zhao J, Zhang P, Qin L, Pan XH. Hypoxia is essential for bone-tendon junction healing: the molecular biological evidence. Int Orthop. 2011;35(6):925–8. https://doi.org/10.1007/s00264-010-1157-7.

    Article  PubMed  Google Scholar 

  12. Yu Y, Zhou Y, Cheng T, Lu X, Yu K, Zhou Y, et al. Hypoxia enhances tenocyte differentiation of adipose-derived mesenchymal stem cells by inducing hypoxia-inducible factor-1α in a co-culture system. Cell Prolif. 2016;49(2):173–84. https://doi.org/10.1111/cpr.12250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Chen G, Zhang W, Zhang K, Wang S, Gao Y, Gu J, et al. Hypoxia-induced mesenchymal stem cells exhibit stronger tenogenic differentiation capacities and promote patellar tendon repair in rabbits. Stem Cells Int. 2020;2020:8822609. https://doi.org/10.1155/2020/8822609.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–72. https://doi.org/10.1016/j.cell.2012.03.042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Ishii KA, Fumoto T, Iwai K, Takeshita S, Ito M, Shimohata N, et al. Coordination of PGC-1beta and iron uptake in mitochondrial biogenesis and osteoclast activation. Nat Med. 2009;15(3):259–66. https://doi.org/10.1038/nm.1910.

    Article  CAS  PubMed  Google Scholar 

  16. Ni S, Yuan Y, Qian Z, Zhong Z, Lv T, Kuang Y, et al. Hypoxia inhibits RANKL-induced ferritinophagy and protects osteoclasts from ferroptosis. Free Radical Biol Med. 2021;169:271–82. https://doi.org/10.1016/j.freeradbiomed.2021.04.027.

    Article  CAS  Google Scholar 

  17. Fink SL, Cookson BT. Apoptosis, pyroptosis, and necrosis: mechanistic description of dead and dying eukaryotic cells. Infect Immun. 2005;73(4):1907–16. https://doi.org/10.1128/iai.73.4.1907-1916.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Chen X, He WT, Hu L, Li J, Fang Y, Wang X, et al. Pyroptosis is driven by non-selective gasdermin-D pore and its morphology is different from MLKL channel-mediated necroptosis. Cell Res. 2016;26(9):1007–20. https://doi.org/10.1038/cr.2016.100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Yang L, Liu J, Shan Q, Geng G, Shao P. High glucose inhibits proliferation and differentiation of osteoblast in alveolar bone by inducing pyroptosis. Biochem Biophys Res Commun. 2020;522(2):471–8. https://doi.org/10.1016/j.bbrc.2019.11.080.

    Article  CAS  PubMed  Google Scholar 

  20. Zhang J, Wei K. Necrosulfonamide reverses pyroptosis-induced inhibition of proliferation and differentiation of osteoblasts through the NLRP3/caspase-1/GSDMD pathway. Exp Cell Res. 2021;405(2):112648. https://doi.org/10.1016/j.yexcr.2021.112648.

    Article  CAS  PubMed  Google Scholar 

  21. Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, et al. A novel ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma. Int J Biol Sci. 2020;16(13):2430–41. https://doi.org/10.7150/ijbs.45050.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ye Y, Dai Q, Qi H. A novel defined pyroptosis-related gene signature for predicting the prognosis of ovarian cancer. Cell Death Discov. 2021;7(1):71. https://doi.org/10.1038/s41420-021-00451-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Shao W, Yang Z, Fu Y, Zheng L, Liu F, Chai L, et al. The pyroptosis-related signature predicts prognosis and indicates immune microenvironment infiltration in gastric cancer. Front Cell Dev Biol. 2021;9:676485. https://doi.org/10.3389/fcell.2021.676485.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Ju A, Tang J, Chen S, Fu Y, Luo Y. Pyroptosis-related gene signatures can robustly diagnose skin cutaneous melanoma and predict the prognosis. Front Oncol. 2021;11:709077. https://doi.org/10.3389/fonc.2021.709077.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Chen J, Long F. mTOR signaling in skeletal development and disease. Bone Res. 2018;6:1. https://doi.org/10.1038/s41413-017-0004-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bae S, Oh B, Tsai J, Park PSU, Greenblatt MB, Giannopoulou EG, et al. The crosstalk between MYC and mTORC1 during osteoclastogenesis. Front Cell Dev Biol. 2022;10:920683. https://doi.org/10.3389/fcell.2022.920683.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Lisse TS, Liu T, Irmler M, Beckers J, Chen H, Adams JS, et al. Gene targeting by the vitamin D response element binding protein reveals a role for vitamin D in osteoblast mTOR signaling. FASEB J Offi Publ Fed Am Soc Exper Biol. 2011;25(3):937–47. https://doi.org/10.1096/fj.10-172577.

    Article  CAS  Google Scholar 

  28. Li XH, Ha CT, Fu D, Xiao M. REDD1 protects osteoblast cells from gamma radiation-induced premature senescence. PLoS ONE. 2012;7(5):e36604. https://doi.org/10.1371/journal.pone.0036604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Jensen ED, Gopalakrishnan R, Westendorf JJ. Regulation of gene expression in osteoblasts. BioFactors (Oxford, England). 2010;36(1):25–32. https://doi.org/10.1002/biof.72.

    Article  CAS  PubMed  Google Scholar 

  30. Yang CM, Chien CS, Yao CC, Hsiao LD, Huang YC, Wu CB. Mechanical strain induces collagenase-3 (MMP-13) expression in MC3T3-E1 osteoblastic cells. J Biol Chem. 2004;279(21):22158–65. https://doi.org/10.1074/jbc.M401343200.

    Article  CAS  PubMed  Google Scholar 

  31. Hammerman M, Blomgran P, Dansac A, Eliasson P, Aspenberg P. Different gene response to mechanical loading during early and late phases of rat Achilles tendon healing (Bethesda, Md: 1985). J Appl Physiol. 2017;123(4):800–15. https://doi.org/10.1152/japplphysiol.00323.2017.

    Article  CAS  PubMed  Google Scholar 

  32. Hou CH, Hou SM, Tang CH. Ultrasound increased BMP-2 expression via PI3K, Akt, c-Fos/c-Jun, and AP-1 pathways in cultured osteoblasts. J Cell Biochem. 2009;106(1):7–15. https://doi.org/10.1002/jcb.21934.

    Article  CAS  PubMed  Google Scholar 

  33. Wu J, Liu T, Rios Z, Mei Q, Lin X, Cao S. Heat Shock Proteins and Cancer. Trends Pharmacol Sci. 2017;38(3):226–56. https://doi.org/10.1016/j.tips.2016.11.009.

    Article  CAS  PubMed  Google Scholar 

  34. Cooper LF, Tiffee JC, Griffin JP, Hamano H, Guo Z. Estrogen-induced resistance to osteoblast apoptosis is associated with increased hsp27 expression. J Cell Physiol. 2000;185(3):401–7. https://doi.org/10.1002/1097-4652(200012)185:3%3c401::aid-jcp10%3e3.0.co;2-c.

    Article  CAS  PubMed  Google Scholar 

  35. Kondo A, Tokuda H, Matsushima-Nishiwaki R, Kato K, Kuroyanagi G, Mizutani J, et al. Unphosphorylated heat shock protein 27 suppresses fibroblast growth factor-2-stimulated vascular endothelial growth factor release in osteoblasts. Mol Med Rep. 2013;8(2):691–5. https://doi.org/10.3892/mmr.2013.1533.

    Article  PubMed  Google Scholar 

  36. Kainuma S, Tokuda H, Yamamoto N, Kuroyanagi G, Fujita K, Kawabata T, et al. Heat shock protein 27 (HSPB1) suppresses the PDGF-BB-induced migration of osteoblasts. Int J Mol Med. 2017;40(4):1057–66. https://doi.org/10.3892/ijmm.2017.3119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Daswani B, Gupta MK, Gavali S, Desai M, Sathe GJ, Patil A, et al. Monocyte proteomics reveals involvement of phosphorylated HSP27 in the pathogenesis of osteoporosis. Dis Markers. 2015;2015:196589. https://doi.org/10.1155/2015/196589.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Notsu K, Nakagawa M, Nakamura M. Ubiquitin-like protein MNSFβ noncovalently binds to molecular chaperone HSPA8 and regulates osteoclastogenesis. Mol Cell Biochem. 2016;421(1–2):149–56. https://doi.org/10.1007/s11010-016-2795-x.

    Article  CAS  PubMed  Google Scholar 

  39. You L, Gu W, Chen L, Pan L, Chen J, Peng Y. MiR-378 overexpression attenuates high glucose-suppressed osteogenic differentiation through targeting CASP3 and activating PI3K/Akt signaling pathway. Int J Clin Exp Pathol. 2014;7(10):7249–61.

    PubMed  PubMed Central  Google Scholar 

  40. Mogi M, Togari A. Activation of caspases is required for osteoblastic differentiation. J Biol Chem. 2003;278(48):47477–82. https://doi.org/10.1074/jbc.M307055200.

    Article  CAS  PubMed  Google Scholar 

  41. Struzik S, Czarkowska-Paczek B, Wyczalkowska-Tomasik A, Maldyk P, Paczek L. Selected clinical features fail to predict inflammatory gene expressions for TNF-α, TNFR1, NSMAF, Casp3 and IL-8 in tendons of patients with rotator cuff tendinopathy. Arch Immunol Ther Exp. 2021;69(1):6. https://doi.org/10.1007/s00005-021-00610-z.

    Article  CAS  Google Scholar 

  42. Patel SH, Sabbaghi A, Carroll CC. Streptozotocin-induced diabetes alters transcription of multiple genes necessary for extracellular matrix remodeling in rat patellar tendon. Connect Tissue Res. 2018;59(5):447–57. https://doi.org/10.1080/03008207.2018.1470168.

    Article  CAS  PubMed  Google Scholar 

  43. Ghasempour G, Mohammadi A, Zamani-Garmsiri F, Najafi M. miRNAs through β-ARR2/p-ERK1/2 pathway regulate the VSMC proliferation and migration. Life Sci. 2021;279:119703. https://doi.org/10.1016/j.lfs.2021.119703.

    Article  CAS  PubMed  Google Scholar 

  44. Ghasempour G, Mohammadi A, Zamani-Garmsiri F, Soleimani AA, Najafi M. Upregulation of TGF-β type II receptor in high glucose-induced vascular smooth muscle cells. Mol Biol Rep. 2022;49(4):2869–75. https://doi.org/10.1007/s11033-021-07100-7.

    Article  CAS  PubMed  Google Scholar 

  45. Ghasempour G, Mahabadi VP, Shabani M, Mohammadi A, Zamani-Garmsiri F, Amirfarhangi A, et al. miR-181b and miR-204 suppress the VSMC proliferation and migration by downregulation of HCK. Microvasc Res. 2021;136:104172. https://doi.org/10.1016/j.mvr.2021.104172.

    Article  CAS  PubMed  Google Scholar 

  46. Soleimani AA, Ghasmpour G, Mohammadi A, Gholizadeh M, Abkenar BR, Najafi M. Focal adhesion kinase-related pathways may be suppressed by metformin in vascular smooth muscle cells in high glucose conditions. Endocrinol Diabetes Metab. 2022;5(4):e351. https://doi.org/10.1002/edm2.351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Risinger GM Jr, Updike DL, Bullen EC, Tomasek JJ, Howard EW. TGF-beta suppresses the upregulation of MMP-2 by vascular smooth muscle cells in response to PDGF-BB. Am J Physiol Cell Physiol. 2010;298(1):C191-201. https://doi.org/10.1152/ajpcell.00417.2008.

    Article  CAS  PubMed  Google Scholar 

  48. Jiang K, Li Y, Xiang C, Xiong Y, Jia J. TGF-β3 regulates adhesion formation through the JNK/c-Jun pathway during flexor tendon healing. BMC Musculoskelet Disord. 2021;22(1):843. https://doi.org/10.1186/s12891-021-04691-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research received no external funding.

Author information

Author notes

  1. Yaonan Zhang and Lei Shi contributed equally to this work and should be considered co-first authors.

Authors and Affiliations

  1. Orthopaedic Department, Beijing Hospital; National Center of Gerontology; Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, People’s Republic of China

    Yaonan Zhang, Lei Shi, Fei Wang, Lin Wang, Nan Min, Liangyuan Wen & Qingyun Xue

Authors
  1. Yaonan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lei Shi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fei Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Lin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nan Min
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Liangyuan Wen
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Qingyun Xue
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YZ and LS contributed to the drafting of the article. LW and QX contributed to the concept and design of the study. FW, LW and NM contributed to the interpretation of data, and the critical revision of the study methods. YZ, LS, LW and QX contributed to the critical revision of the article for relevant intellectual content. All authors approved the final version of the article, including the authorship list.

Corresponding authors

Correspondence to Liangyuan Wen or Qingyun Xue.

Ethics declarations

Ethics approval and consent to participate

All animals were treated according to institutional guidelines for laboratory animal treatment and care. The study is reported in accordance with ARRIVE guidelines. All experimental procedures were approved by the Animal Research Ethics Committee of our hospital.

Consent for publication

Not applicable.

Competing interests

The authors indicated no potential conflicts of interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1. Figure S1

GO (A) and KEGG (B) enrichment analysis of DEmRNAs between injury model and celecoxib treatment group BP: biological process; CC: cytological component; MF: molecular function.

Additional file 2. Figure S2

GO (A) and KEGG (B) enrichment analysis of DEmRNAs between injury model group and celecoxib + lactoferrin treatment group BP: biological process; CC: cytological component; MF: molecular function.

Additional file 3. Figure S3

Identification of celecoxib + lactoferrin treatment group-specific DEmRNAs (A) Venn diagram of up-regulated DEmRNAs in injury model group vs celecoxib treatment group and injury model group vs celecoxib + lactoferrin treatment group. (B) Venn diagram of down-regulated DEmRNAs in injury model group vs celecoxib treatment group and injury model group vs celecoxib + lactoferrin treatment group.

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

Zhang, Y., Shi, L., Wang, F. et al. Screening for autophagy/hypoxia/ferroptosis/pyroptosis-related genes of tendon injury and repair in a rat model after celecoxib and lactoferrin treatment. J Orthop Surg Res 18, 383 (2023). https://doi.org/10.1186/s13018-023-03856-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-023-03856-9

Keywords

Journal of Orthopaedic Surgery and Research

ISSN: 1749-799X