Whole Genome Expression Profiling of Semitendinosus Tendons from Children with Diplegic and Tetraplegic Cerebral Palsy

Cerebral palsy (CP) is the most common movement disorder in children with a prevalence rang-ing from 1.5 to 4 per 1000 live births. CP is caused by a non-progressive lesion of the developing brain, leading to progressive alterations of the musculoskeletal system, including spasticity, often leading to the development of fixed contractures, necessitating tendon lengthening surgery. Total RNA-sequencing analysis was performed on semitendinosus tendons from diplegic and tetraple-gic CP patients subjected to tendon lengthening surgery compared to control patients undergoing anterior cruciate ligament reconstructive surgery. Tetraplegic CP patients showed increased ex-pression of genes implicated in collagen synthesis and extracellular matrix (ECM) turnover, while only minor changes were observed in diplegic CP patients. In addition, tendons from tetra-plegic CP patients showed an enrichment for upregulated genes involved in vesicle-mediated transport and downregulated genes involved in cytokine and apoptotic signaling. Overall, our results indicate increased ECM turnover with increased net synthesis of collagen in tetraplegic CP patients without activation of inflammatory and apoptotic pathways similar to observations in athletes where ECM remodeling results in increased tendon stiffness and tensile strength. Never-theless, the resulting increased tendon stiffness is an important issue in clinical practice, where surgery is often required to restore joint mobility.

Keywords:

Subject: Biology and Life Sciences - Biochemistry and Molecular Biology

1. Introduction

Tendons are fibrous connective tissues that connect skeletal muscle to bone and transmit the mechanical force generated by muscle contraction to the bones, driving the movement of the skeletal system. Tendons are primarily composed of cells and extracellular matrix (ECM), principally comprised of collagen type I, which accounts for about 65-80% of the dry weight and is organized into parallel fiber bundles arranged into a three-dimensional network surrounded by a thin connective tissue sheath known as the epitenon [1]. Additionally, the ECM contains proteoglycans (e.g. aggrecan, versican, decorin, biglycan, fibromodulin, lumican, keratocan, osteoglycin, syndecans), glycoproteins (e.g. elastin, cartilage oligomeric matrix protein (COMP/TSP5), tenascin-C, lubricin, tenomodulin, fibronectin), glycosaminoglycans, and smaller quantities of other collagens with different functions. Aggrecan and versican are large aggregating proteoglycans, whose main function is to provide compression resistance to the tendon by increasing water content [2]. Decorin, biglycan, fibromodulin, lumican, and osteoglycin are small leucine-rich proteoglycans, which play a role in fibril assembly and tendon integrity [3,4,5], while syndecans are transmembrane heparan sulfate proteoglycans whose roles in tendons have remained elusive [6]. Thus, while collagen 1 fibers are responsible for the tensile strength of tendons, proteoglycans are regulating their viscoelastic properties. Tenoblasts and tenocytes constitute 90-95% of tendon cells and are fibroblast-like cells responsible for the production of collagen and other ECM components [1]. Tenocytes respond to mechanical load of tendons by altering the expression of ECM proteins, thereby modulating tendon structure, composition, and mechanical properties [7].

Cerebral palsy (CP) is a group of disorders of the development of movement and posture, causing activity limitations, which are attributed to non-progressive brain lesions that occurred in the developing fetus or infant [8]. CP is the most common chronic movement disorder in children with a prevalence of 1.5 to 4 per 1000 live births [9,10,11,12,13,14,15], and one of the most disabling and costly chronic disorders of children and adults. Although CP is the result of a non-progressive lesion of the developing brain, it leads to progressive alterations of the musculoskeletal system, which can manifest with different clinical presentations, such as spasticity, dyskinesia, dystonia, hypertonia, often leading to the development of fixed contractures, i.e. a permanent shortening and stiffening of the muscle-tendon unit, resulting in the loss of joint motility (reviewed in [16,17]). Patients are diagnosed in respect to the anatomical distribution of their deformity as monoplegic, hemiplegic, diplegic, or tetraplegic (reviewed in [16]). Involvement of only one limb is referred to as monoplegia, unilateral involvement of an upper and lower limb is referred to as hemiplegia, predominant lower limb involvement is referred to as diplegia, involvement of the lower limbs and one upper limb is referred to as triplegia, and involvement of all four limbs and the trunk is referred to as quadriplegia. Various treatment strategies, including physiotherapy, pharmacological interventions, neurectomy, and orthosis, are used to reduce spasticity, limit pain, and prevent contractures [16,17]. Nevertheless, most CP patients eventually develop fixes contractures, necessitating tendon lengthening surgery to restore the joint range of motion and relieve symptoms. In previous studies, our coworkers found that tendons of CP patients show ECM remodeling [18] and altered expression of selected genes related to collagen turnover [19,20]. However, the effect of CP on genome-wide gene expression has not been studied and it remains poorly understood how the exposure of tendons to the spasticity-induced increased mechanical loading and functional demands in CP patients affect tendon homeostasis at the molecular level. In the present study we performed total RNA-sequencing (RNA-Seq) analysis on semitendinosus tendons from diplegic and tetraplegic CP patients undergoing tendon-lengthening surgery compared to age-matched controls undergoing cruciate ligament (ACL) reconstruction.

2. Materials and Methods

2.1. Sample collection and preparation

A cohort of 30 children were included in the study: 12 children with diplegic cerebral palsy (10 males and 2 females; mean age: 13.9 years), 13 children with tetraplegic cerebral palsy (8 males and 5 females; mean age: 12.6 years), and 5 typically developing pediatric patients as control patients (4 males and 1 female; mean age: 17.9). The characteristics of the subjects are listed in Table 1. Semitendinosus tendon biopsies were obtained during surgery of cerebral palsy patients undergoing tendon lengthening or control patients undergoing ACL reconstruction. Biopsies were immediately frozen on dry ice and stored at –80°C for further processing.

2.2. RNA extraction

For extraction of tendon RNA, tendons were first powderized with a hammer on dry ice after which TRIzol Reagent (Thermo Fisher Scientific) was added and tendons were homogenized using a TissueLyser II (Qiagen). RNA was subsequently extracted using the RNA-ZOL Direct Clean-up Kit (Fisher Molecular Biology) following the instructions of the manufacturer. After DNase I treatment, RNA purity and concentration were evaluated using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific) and a Qubit Fluorometer (Thermo Fisher Scientific). RNA quality was analyzed on an Agilent Tape Station 4200 system and the RNA integration number (RIN) varied from 2.7 – 9.1.

2.3. RNA-Seq analysis

Indexed sequencing libraries were generated from 140 ng of tendon RNA using the KAPA RNA HyperPrep with RiboErase (Roche), which allows for library construction from partially degraded samples. Paired-end multiplexed sequencing of libraries to generate reads of 150 bp (PE 2x250) was performed on a NovaSeq 6000 instrument (Illumina Inc.). 22-40 million reads per sample were obtained. RNA-Seq samples were first demultiplexed and FASTQ files were created from BCL files using bcl2fastq v2.17.1.14 software (Illumina Inc.). Quality control and assessment were performed on FASTQ files using FastQC v0.11.9; http://www.bioinformatics.babraham.ac.uk/projects/fastqc). Paired-end reads of 150 bp were aligned to the Gencode Human reference genome (build GRCh38) using STAR v2.7.2b [21]. Annotation of genes was performed using the biomaRt package v2.54 [22] based on annotations from Ensembl release 110. The raw read counts were normalized with TMM implemented in the edgeR v3.40 package [23] in R/Bioconductor [24] and lesser expressed genes were filtered out with the filterByExpr (min.count=6) function. Differential expression analysis of read counts was performed using the voom, lmFit, and eBayes (robust=T) functions of the limma v3.46 package [25] in R. Significantly differentially expressed genes were chosen based on a false discovery rate (FDR) < 0.1 and |log2 fold change (FC)| ≥ 0.4. The Sva v3.38 package [26] was used to estimate artifacts (n.sv=1) and correct the CPM values. Hierarchical clustering of significantly modulated genes was performed using hclust and dist functions on sva corrected log2CPM in R. Clustering was performed with the ward.D2 method and Euclidean distance to generate a heatmap using pheatmap (https://cran.r-project.org/web/packages/pheatmap/index.html) and scaling the rows. Principal components analysis (PCA) was performed using the PCA function of the FactoMineR v2.7 [27] package in R. The Venn Diagrams were obtained using the R package VennDiagram (https://cran.r-project.org/web/packages/VennDiagram/index.html).

Raw data files for the RNA-Seq analysis have been deposited in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/) under accession code PRJNA1004310.

2.4. Gene ontology and pathway analysis

For gene enrichment analysis, the gene ontology (GO) Biological Process 2021 and Reactome 2022 databases within EnrichR [28] were used to find significant enriched GO terms and pathways with an adjusted p value ≤ 0.05.

3. Results

Total RNA-Seq analysis was performed on semitendinosus tendons from 12 diplegic, 13 tetraplegic and 5 control patients. Unfortunately, due to muscle contamination 5 control samples had to be excluded. Principal component analysis (PCA) (Figure 1A) and hierarchical clustering showed clustering of the control group distant from the CP groups (Figure 1A,B). 1170 genes (676 upregulated and 494 downregulated) were differentially expressed in tetraplegic vs. control patients, 112 in diplegic vs. control patients (61 upregulated and 51 downregulated), and 101 (70 upregulated and 31 downregulated) in tetraplegic vs. diplegic patients. As illustrated in the Venn diagram in Figure 1C, 101 genes (58 upregulated and 43 downregulated) were commonly altered in diplegic and tetraplegic patients, while 1068 and 11 genes were differentially expressed specifically in tetraplegic vs. control patients and diplegic vs. control patients, respectively. Thus, most genes that were differentially expressed in diplegic patients were altered also in tetraplegic patients. The full lists of differentially expressed (DE) genes are shown in Table S1-S3. For genes that were differentially expressed in tetraplegic vs. control patients, gene ontology analysis for biological processes as well as Reactome pathway analysis were performed (Figure 2 and Tables S4–S7). For upregulated genes, there was an enrichment for genes involved in collagen formation, ECM organization, and categories related to vesicle-mediated transport (Figure 2A,B, left), while for downregulated genes, genes related to cytokine and apoptotic signaling, including TNF-related apoptosis-inducing agent (TRAIL) signaling, were most significantly affected (Figure 2A,B, right). The differentially expressed genes in diplegic patients did not show significant enrichment for genes involved to any specific process. Differentially expressed genes in tetraplegic CP patients vs. control patients involved in ECM organization, vesicle-mediated transport, cytokine signaling, and apoptotic signaling are listed in Figure 3. Many genes encoding collagens or involved in collagen formation and organization were upregulated, including proteoglycans, glycoproteins, integrins, and growth factors (Figure 3A). Increased expression of several genes encoding metallopeptidases (MMPs), involved in ECM turnover, was also found, including MMP2, which is the most expressed MMP in tendons. Furthermore, many genes involved in membrane trafficking were upregulated (Figure 3B). Among the downregulated genes were many genes involved in cytokine signaling, including interferon signaling, interleukin signaling, and chemokines (Figure 3C) as well as a number of genes related to TRAIL signaling and apoptosis (Figure 3D).

4. Discussion

The present study is the first to report genome-wide expression profiling on tendons from CP patients. Consistent with our previous findings [19], tetraplegic CP patients showed increased expression of genes encoding collagens and ECM proteins. The most upregulated genes were the collagen 1-encoding genes COL1A1 and COL1A2, which were 15- and 8-fold upregulated, respectively. Furthermore, genes encoding both fibrillar (COL5A1/2, COL11A1/2) and non-fibrillar collagens (COL6A1/2/3 (beaded filament forming collagen), COL8A1/2 (hexagonal network forming collagen), COL10A1 (hexagonal network forming collagen), COL14A1 (fibril-associated collagen with interrupted triple-helix (FACIT)), COL15A1 (FACIT)) were highly upregulated as shown in Figure 3A [29].

Collagen is composed of a triple helix, consisting of three α chains with the repeating (Gly-X-Y)n amino acid sequence, where X and Y are frequently proline or hydroxyproline but can be any amino acid [29]. While some collagens contain three identical α chains, others contain two or three different chains referred to as α1, α2, and α3. In tendons, fibrillar collagens are bundled into fibrils, which form larger fibers. During collagen biosynthesis, collagen undergoes post-translational modification by hydroxylation of proline residues into hydroxyproline by prolyl 3-hydroxylases (P3H1, P3H2, P3H3, P3H4) and prolyl 4-hydroxylases (a tetramer composed of two α (P4H1, P4H2, P4H3) and one β (P4HB) subunit) as well as hydroxylation of lysine into hydroxylysine by lysyl hydroxylases (LH1, LH2, LH3), encoded by procollagen-lysine,2-oxoglutarate 5-dioxygenases (PLOD1, PLOD2, PLOD3) [30,31]. Hydroxylysine residues are often subject to glycosylation and P4HB is responsible for the formation of disulfide bands between the α chains. Peptidyl-prolyl cis-trans isomerases (FKBP10, FKBP14, PPIB) catalyze the cis-trans isomerization of peptide bonds, which is a rate limiting step for triple helix formation. After secretion, lysyl oxidases (LOX, LOXLs) oxidize hydroxylysine residues to induce cross-link formation [32], which provide the collagen fibrils with mechanical stability and tensile strength as well as contribute to the stiffness of the collagen fibril. In tendons of tetraplegic CP patients, genes involved in all steps of collagen formation and cross-linking were upregulated, including P3H3, P3H4, P4HA1, P4HA3, P4HB, PLOD1, FKBP10, FKBP14, PPIB, LOXL1, LOXL2, and LOXL3, while only PPIB was increased in diplegic CP patients. This suggests that collagen synthesis is strongly induced in tetraplegic CP patients.

In tetraplegic CP patients, genes encoding the small leucine-rich proteoglycans biglycan (BGN), keratocan (KERA), and osteoglycin (OGN) were highly upregulated by 2.4-, 7.4- and 4.7-fold, respectively. Biglycan was recently shown to be involved in the maintenance of tendon structure and mechanics in mature tendons [33] and osteoglycin has been implicated in collagen fibrillar organization and tendon mechanical function [5,34], while the role of keratocan has remained elusive [35]. Decorin, which constitutes about 80% of the total amount of small leucine-rich proteoglycans in tendons, was not affected. This may be explained by a recent study, which demonstrated that while decorin plays a more important role than biglycan (the second most abundant small leucine-rich proteoglycan in tendons) in the modulation of collagen fibril structure and viscoelastic mechanics during tendon development, biglycan is more important for the maintenance of tendon structure and mechanical properties during homeostasis in mature tendons [33]. While the expression of none of the genes encoding small leucine-rich proteoglycans was altered in diplegic CP patients, SDC2, encoding the transmembrane heparan sulfate proteoglycan syndecan-2, was upregulated in both diplegic and tetraplegic CP patients by 1.7- and 1.6-fold, respectively.

Several genes encoding glycoproteins, including tenascin-C (TNC), fibronectin 1 (FN1), tenomodulin (TNMD), secreted protein acidic and rich in cysteine (SPARC), and thrombospondin 4 (THBS4) were upregulated in tetraplegic CP patients. In particular, SPARC and THBS4 were strongly upregulated by 6.0- and 5.6-fold, respectively. SPARC has been shown to be essential for load-induced tendon tissue maturation and homeostasis, affecting ECM composition and tendon biomechanical properties [36], while tenomodulin plays a role in tenocyte proliferation and collagen fibril maturation [37]. Thrombospondin 4, encoded by THBS4, which was 2.9-fold upregulated, is involved in the organization of collagen fibrils [37], while tenascin C and fibronectin, whose genes were 3.3- and 2.2-fold upregulated, respectively, both contribute to the mechanical stability of the ECM [38,39]. Among the glycoproteins, only TNC was upregulated (2.2-fold) in diplegic CP patients. In tetraplegic CP patients, alterations in the expression of genes encoding different subunits of the large multidomain heterotrimeric glycoprotein laminin was also found. Laminin is located in the basement membrane, where it interacts with collagen type IV, integrins, and dystroglycans, and play important roles in cell adhesion, differentiation, and migration [40]. LAMA1 (encoding laminin α1) and LAMB2 (encoding laminin β2) were 3.3- and 1.5-fold upregulated, respectively, while LAMC2 (encoding laminin γ2) was 3.3-fold downregulated.

Altered mRNA expression of several integrins was also found in the tendons of tetraplegic CP patients, including ITGA2, ITGB1, and ITGBL1, which were 3.5-, 1.4,- and 2.2-fold upregulated, respectively as well as ITGA3, which was 1.8-fold downregulated. Integrins are heterodimeric transmembrane receptors, which play a major role in linking the ECM to the cytoskeleton [41]. It is believed that integrins can sense and transmit mechanical stimuli from the ECM to tenocytes, thereby triggering intracellular signaling pathways leading to adaptive regulation of gene expression. In particular, mechanical stretch was found to activate the AKT/mTOR pathway via β1 integrin, thereby regulating collagen expression [42]. Also, GJA1, encoding connexin 43, a gap junction protein with an important role in the communications between tenocytes allowing for the passage of passage of free metabolites and ions [43], was 2.7-fold upregulated in tetraplegic CP patients.

The expression of number of genes encoding growth factors were increased in tendons of tetraplegic CP patients, including integrin like growth factor 1 (IGF1), integrin like growth factor binding protein 3 (IGFBP3), fibroblast growth factor 12 (FGF12), transforming growth factors (TGFB1, TGFB2, TGFB3), and vascular endothelial growth factor B (VEGFB), while FGF11 was downregulated (see Figure 3A). IGF1 has been shown to be required for adult tendon growth in response to increased mechanical overload through stimulation of tenocyte proliferation and protein synthesis [44]. Similarly, TGFβ was shown to promote collagen synthesis and matrix remodeling during tendon healing, thereby enhancing mechanical strength [45]. On the other hand, excessive TGFβ activation as a result of mechanical overload or repetitive mechanical loading can lead to tendinopathy [46]. VEGF promotes angiogenesis and increases vascular permeability and is important for neovascularization during tendon healing [47]. Also, the MKX gene encoding Mohawk, an atypical homeobox transcription factor involved in postnatal tendon maturation [48,49] and critical for the tendon response to mechanical stimuli [50], was 2.6-fold upregulated in tetraplegic CP patients.

Among the differentially expressed genes in the tendons of tetraplegic CP patients were many genes encoding metalloproteases, including matrix metalloproteinases (MMP2, MMP10, MMP14, MMP16) and members of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) protein family (ADAMTS2, ADAMTS3, ADAMTS17, ADAMTSL2). Metalloproteases are important for the turnover of the ECM, which is a finely balanced, dynamic process of protein synthesis and degradation taking place at a low rate during homeostasis and at higher rates during conditions of inflammation, tissue damage, or increased mechanical load [51]. Metalloproteinases are mediating the proteolytic degradation of components of the ECM and are regulated by tissue inhibitors of metalloproteinases (TIMS). While ECM turnover allows for tendon adaptation in response to altered mechanical loading conditions, an imbalance in the activities of MMPs and TIMS can lead to pathological conditions [52]. MMP2, MMP14, and MMP16, whose mRNA levels were 3.5-, 2.0-, and 3.6-fold upregulated, respectively, are able to cleave fibrillar collagen I, while MMP10, whose gene was 7.4-fold downregulated, does not cleave collagen I [52], but degrades various components of the ECM, including proteoglycans, fibronectin, and collagen III-V [53] as well as activates other MMPs [54]. ADAMTS2 and ADAMTS3, which were 3.4- and 3.3-fold upregulated at the mRNA level, respectively, are involved in fibrillar collagen maturation though cleavage of the N-terminal propeptide of procollagens [55,56]. Also, PCOLCE, encoding procollagen C-endopeptidase enhancer, which enhances the activity of procollagen C proteinases cleaving the C-terminal propeptide of procollagens, was 2.6-fold upregulated [57]. In addition, mRNA encoding ADAMTS17, whose function in tendons has remained elusive, was 2.5-fold upregulated. On the others hand, transcript levels of ADAMTSL2, which belongs to the subfamily of ADAMTS-like proteins without catalytic activity and is involved in the modulation of microfibril formation, were 2.1-fold reduced [58]. Among the metalloproteases, only MMP2 was 2.1-fold upregulated in diplegic CP patients.

Consistent with increased protein turnover and remodeling in tetraplegic CP patients, Reactome analysis and gene ontology analysis for biological processes revealed an enrichment for upregulated genes involved in vesicle-mediated transport (Figure 2A,B, left). Among the downregulated genes, there was an enrichment for genes involved in interferon and cytokine signaling, including a number of chemokines (Figure 2A,B, right, Figure 3C). In particular, CCL5, CCL16, CXCL10, and CXCL11 were strongly downregulated by respectively 3.5-, 2.1-, 5.6-, and 4.4-fold, while CXCL16 was 1.7-fold upregulated. Thus, the higher mechanical load experienced by tendons of CP patients does not appear to be associated with inflammation. Furthermore, apoptosis does not seem to be activated as there was an enrichment for downregulated genes associated with TRAIL signaling and apoptotic processes (Figure 3D).

The results of our RNA-Seq analysis are consistent with increased collagen synthesis and turnover in tendons from tetraplegic CP patients as a consequence of the spasticity-induced chronic mechanical stimulation. Many of the changes in gene expression observed in CP patients are similar to those observed in the tendons of patients with tendinopathy, but there are important differences suggesting that the altered gene expression in tetraplegic CP patients may be adaptive rather than pathological as in tendinopathy. Like in tendinopathy patients the expression of many collagens was increased. However, while the mRNA expression of a number of collagens was 1.6-2.7-fold upregulated in tendinopathy patients, COL1A1, COL1A2, COL10A1, and COL11A1/2 were 5.3-15.2-fold upregulated in tetraplegic CP patients, suggesting strongly increased collagen synthesis. In regard to the expression of metalloprotease-encoding genes, increased expression of MMP2, MMP14, MMP16, ADAMTS2, ADAMTS3 as well as downregulation of MMP10 has also been found in tendinopathy [59,60]. Also, upregulation of genes encoding proteoglycans (BGN) glycoproteins (TNC, FN1, SPARC), integrins (ITGB1), and growth factors has been reported in tendinopathy patients [60,61,62]. On the other hand, we found no alterations in genes associated with Notch3, Wnt, and JAK/STAT signaling, which were reported to be increased in tendinopathy patients ([60]). In addition, although the role of inflammation in tendinopathy remains controversial [63,64,65,66], a number of genes encoding pro-inflammatory cytokines have been reported to be upregulated in tendinopathy, whereas many chemokines and interleukins were downregulated in tetraplegic CP patients [60,62] (see Figure 3C). Moreover, whereas increased expression of markers of apoptosis was found in tendinopathy [62], many apoptosis-related genes were downregulated in tetraplegic CP patients. In addition, an enrichment for upregulated genes associated with vesicle-mediated transport has not been reported in tendinopathy patients. Overall, our results indicate that tetraplegic CP patients show increased ECM turnover with an increased net synthesis of collagen to allow for reorganization of the ECM as an adaptive response to increased mechanical loads and functional demands as a result of spasticity. This is similar to observations after exercise training, where ECM remodeling has been found to result in reduced tendon stress and increased tendon stiffness and tensile strength [67,68,69,70,71,72,73], likely generating more resistant tendons to support higher loads. The absence of pro-inflammatory cytokines and markers of apoptosis in CP patients compared to tendinopathy patients could be due to the difference in stress stimuli, being sub maximal and continuous in the former and over maximal and discontinuous in the ladder. Although tendon remodeling in CP patients does not appear to be directly pathological, increased tendon stiffness and shortening in CP patients as a result of continuous mechanical stimulation is an important issue in clinical practice, where tendon extension surgery is often required to restore joint mobility. In clinal practice, various kinds of muscle relaxants and neuromuscular blocking agents are also often used to reduce spasticity and consequently decrease pathological stimulation of tendons and prevent the development of contractures [16,17]. In future studies it would be interesting to perform mechanical tension studies and histological analyses of patients with different severity of CP to determine whether they show increased tendon stiffness and cross-sectional area as observed in tendons from long distance runners [74,75] and after high resistance training [68,69,71,72,76]. The much fewer gene expression changes in diplegic compared to tetraplegic CP patients implies that tendons of patients with a less severe form of CP are not affected by the disease to an extent that induces major tendon remodeling, consistent with the fact that only intense exercise is associated with substantial tendon remodeling [68,74,76]. This is also consistent with our personal observations in the clinic, where increased tendon stiffness is generally more evident in tetraplegic CP patients compared to diplegic CP patients.

In conclusion, the present study provides new insights into the alterations in gene expression in CP patients as a consequence of the chronically increased mechanical loads due to muscle hyperactivity. A better understanding of the molecular alterations in the tendons of CP patients during the course of the disease may be beneficial for disease monitoring and lead to better treatment strategies to improve quality of life in CP patients.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Table S1: Differentially expressed genes in tetraplegic CP patients vs. control patients. Table S2: Differentially expressed genes in diplegic CP patients vs. control patients. Table S3: Differentially expressed genes in tetraplegic CP patients vs. diplegic CP patients. Table S4: Gene ontology analysis for biological processes for upregulated genes in tetraplegic CP patients vs. control patients. Table S5: Gene ontology analysis for biological processes for downregulated genes in tetraplegic CP patients vs. control patients. Table S6: Reactome analysis for upregulated genes in tetraplegic CP patients vs. control patients. Table S7: Reactome analysis for downregulated genes in tetraplegic CP patients vs. control patients.

Author Contributions

Conceptualization, N.M.P. and M.L.B.; Data curation, S.S.; Validation, M.L.B.; Formal Analysis, S.S., M.L.B., G.B.; Investigation, S.N., V.L, M.L.B.; Resources, N.M.P., M.L.B.; Writing – Original Draft Preparation, M.L.B.; Writing – Review & Editing, M.L.B., N.M.P., S.S., G.B., S.N.; Visualization, M.L.B.; Supervision, M.L.B.; Project administration, M.L.B.; Funding Acquisition, N.M.P., G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Ariel Foundation (www.fondazioneariel.it).

Institutional Review Board Statement

The study was conducted in according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Humanitas Research Hospital, Rozzano (Milan), Italy (#1102, April 11, 2013).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are openly available in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/) under accession code PRJNA1004310. The data presented in this study are available in Table S1-S7.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

Kannus, P. Structure of the tendon connective tissue. Scand J Med Sci Sports 2000, 10, 312–320. [Google Scholar] [CrossRef] [PubMed]
Vogel, K.G. What happens when tendons bend and twist? Proteoglycans. J Musculoskelet Neuronal Interact 2004, 4, 202–203. [Google Scholar] [PubMed]
Ezura, Y.; Chakravarti, S.; Oldberg, A.; Chervoneva, I.; Birk, D.E. Differential expression of lumican and fibromodulin regulate collagen fibrillogenesis in developing mouse tendons. J Cell Biol 2000, 151, 779–788. [Google Scholar] [CrossRef] [PubMed]
Dunkman, A.A.; Buckley, M.R.; Mienaltowski, M.J.; Adams, S.M.; Thomas, S.J.; Satchell, L.; Kumar, A.; Pathmanathan, L.; Beason, D.P.; Iozzo, R.V.; et al. Decorin expression is important for age-related changes in tendon structure and mechanical properties. Matrix Biol 2013, 32, 3–13. [Google Scholar] [CrossRef]
Boote, C.; Ma, Q.; Goh, K.L. Age-dependent mechanical properties of tail tendons in wild-type and mimecan gene-knockout mice - A preliminary study. J Mech Behav Biomed Mater 2023, 139, 105672. [Google Scholar] [CrossRef]
Gopal, S.; Arokiasamy, S.; Pataki, C.; Whiteford, J.R.; Couchman, J.R. Syndecan receptors: pericellular regulators in development and inflammatory disease. Open Biol 2021, 11, 200377. [Google Scholar] [CrossRef]
Wang, J.H. Mechanobiology of tendon. J Biomech 2006, 39, 1563–1582. [Google Scholar] [CrossRef]
Bax M, G.M. , Rosenbaum P, Leviton A, Paneth N, Dan B, Jacobsson B, Damiano D. Proposed definition and classification of cerebral palsy, April 2005. Developmental Medicine & Child Neurology 2005, 47, 571–571. [Google Scholar]
Paneth, N.; Hong, T.; Korzeniewski, S. The Descriptive Epidemiology of Cerebral Palsy. Clinics in Perinatology 2006, 33, 251–267. [Google Scholar] [CrossRef]
Yeargin-Allsopp, M.; Van Naarden Braun, K.; Doernberg, N.S.; Benedict, R.E.; Kirby, R.S.; Durkin, M.S. Prevalence of cerebral palsy in 8-year-old children in three areas of the United States in 2002: a multisite collaboration. Pediatrics 2008, 121, 547–554. [Google Scholar] [CrossRef]
Europe. , S.o.C.P.i. Prevalence and characteristics of children with cerebral palsy in Europe. Dev Med Child Neurol 2002, 44, 633–640. [Google Scholar]
Van Naarden Braun, K.; Doernberg, N.; Schieve, L.; Christensen, D.; Goodman, A.; Yeargin-Allsopp, M. Birth Prevalence of Cerebral Palsy: A Population-Based Study. Pediatrics 2016, 137, 1–9. [Google Scholar] [CrossRef] [PubMed]
El-Tallawy, H.N.; Farghaly, W.M.; Shehata, G.A.; Rageh, T.A.; Metwally, N.A.; Badry, R.; Sayed, M.A.; Abd El Hamed, M.; Abd-Elwarth, A.; Kandil, M.R. Cerebral palsy in Al-Quseir City, Egypt: prevalence, subtypes, and risk factors. Neuropsychiatr Dis Treat 2014, 10, 1267–1272. [Google Scholar] [PubMed]
Chang, M.J.; Ma, H.I.; Lu, T.H. Estimating the prevalence of cerebral palsy in Taiwan: A comparison of different case definitions. Res Dev Disabil 2015, 36C, 207–212. [Google Scholar] [CrossRef]
Christensen, D.; Van Naarden Braun, K.; Doernberg, N.S.; Maenner, M.J.; Arneson, C.L.; Durkin, M.S.; Benedict, R.E.; Kirby, R.S.; Wingate, M.S.; Fitzgerald, R.; et al. Prevalence of cerebral palsy, co-occurring autism spectrum disorders, and motor functioning - Autism and Developmental Disabilities Monitoring Network, USA, 2008. Dev Med Child Neurol 2014, 56, 59–65. [Google Scholar] [CrossRef] [PubMed]
Graham, H.K.; Rosenbaum, P.; Paneth, N.; Dan, B.; Lin, J.P.; Damiano, D.L.; Becher, J.G.; Gaebler-Spira, D.; Colver, A.; Reddihough, D.S.; et al. Cerebral palsy. Nature reviews. Disease primers 2016, 2, 15082. [Google Scholar] [CrossRef] [PubMed]
Handsfield, G.G.; Williams, S.; Khuu, S.; Lichtwark, G.; Stott, N.S. Muscle architecture, growth, and biological Remodelling in cerebral palsy: a narrative review. BMC Musculoskelet Disord 2022, 23, 233. [Google Scholar] [CrossRef] [PubMed]
Gagliano, N.; Menon, A.; Martinelli, C.; Pettinari, L.; Panou, A.; Milzani, A.; Dalle-Donne, I.; Portinaro, N.M. Tendon structure and extracellular matrix components are affected by spasticity in cerebral palsy patients. Muscles Ligaments Tendons J 2013, 3, 42–50. [Google Scholar] [CrossRef]
Gagliano, N.; Pelillo, F.; Chiriva-Internati, M.; Picciolini, O.; Costa, F.; Schutt, R.C., Jr.; Gioia, M.; Portinaro, N. Expression profiling of genes involved in collagen turnover in tendons from cerebral palsy patients. Connect Tissue Res 2009, 50, 203–208. [Google Scholar] [CrossRef]
Gagliano, N.; Pelillo, F.; Grizzi, F.; Picciolini, O.; Gioia, M.; Portinaro, N. Gene expression profile of extracellular matrix of tendons in cerebral palsy. Dev Med Child Neurol 2007, 49, 557–558. [Google Scholar] [CrossRef]
Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
Durinck, S.; Spellman, P.T.; Birney, E.; Huber, W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc 2009, 4, 1184–1191. [Google Scholar] [CrossRef]
Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
Gentleman, R.C.; Carey, V.J.; Bates, D.M.; Bolstad, B.; Dettling, M.; Dudoit, S.; Ellis, B.; Gautier, L.; Ge, Y.; Gentry, J.; et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5, R80. [Google Scholar] [CrossRef]
Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015, 43, e47. [Google Scholar] [CrossRef]
Leek, J.T.; Johnson, W.E.; Parker, H.S.; Jaffe, A.E.; Storey, J.D. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 2012, 28, 882–883. [Google Scholar] [CrossRef]
Le, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. J Stat Softw 2008, 25, 1–18. [Google Scholar] [CrossRef]
Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res 2016, 44, W90–97. [Google Scholar] [CrossRef] [PubMed]
Ricard-Blum, S. The collagen family. Cold Spring Harb Perspect Biol 2011, 3, a004978. [Google Scholar] [CrossRef] [PubMed]
Heino, J. The collagen family members as cell adhesion proteins. Bioessays 2007, 29, 1001–1010. [Google Scholar] [CrossRef] [PubMed]
Gjaltema, R.A.; Bank, R.A. Molecular insights into prolyl and lysyl hydroxylation of fibrillar collagens in health and disease. Crit Rev Biochem Mol Biol 2017, 52, 74–95. [Google Scholar] [CrossRef]
Ellingson, A.J.; Pancheri, N.M.; Schiele, N.R. Regulators of collagen crosslinking in developing and adult tendons. Eur Cell Mater 2022, 43, 130–152. [Google Scholar] [CrossRef] [PubMed]
Beach, Z.M.; Bonilla, K.A.; Dekhne, M.S.; Sun, M.; Adams, T.H.; Adams, S.M.; Weiss, S.N.; Rodriguez, A.B.; Shetye, S.S.; Birk, D.E.; et al. Biglycan has a major role in maintenance of mature tendon mechanics. J Orthop Res 2022, 40, 2546–2556. [Google Scholar] [CrossRef] [PubMed]
Tasheva, E.S.; Koester, A.; Paulsen, A.Q.; Garrett, A.S.; Boyle, D.L.; Davidson, H.J.; Song, M.; Fox, N.; Conrad, G.W. Mimecan/osteoglycin-deficient mice have collagen fibril abnormalities. Mol Vis 2002, 8, 407–415. [Google Scholar] [PubMed]
Rees, S.G.; Waggett, A.D.; Kerr, B.C.; Probert, J.; Gealy, E.C.; Dent, C.M.; Caterson, B.; Hughes, C.E. Immunolocalisation and expression of keratocan in tendon. Osteoarthritis Cartilage 2009, 17, 276–279. [Google Scholar] [CrossRef] [PubMed]
Wang, T.; Wagner, A.; Gehwolf, R.; Yan, W.; Passini, F.S.; Thien, C.; Weissenbacher, N.; Lin, Z.; Lehner, C.; Teng, H.; et al. Load-induced regulation of tendon homeostasis by SPARC, a genetic predisposition factor for tendon and ligament injuries. Sci Transl Med 2021, 13. [Google Scholar] [CrossRef]
Alexandrov, V.P.; Naimov, S.I. A Prospectus of Tenomodulin. Folia Med (Plovdiv) 2016, 58, 19–27. [Google Scholar] [CrossRef]
Wang, J.H.; Guo, Q.; Li, B. Tendon biomechanics and mechanobiology--a minireview of basic concepts and recent advancements. J Hand Ther 2012, 25, 133–140; quiz 141. [Google Scholar] [CrossRef]
Jones, F.S.; Jones, P.L. The tenascin family of ECM glycoproteins: structure, function, and regulation during embryonic development and tissue remodeling. Dev Dyn 2000, 218, 235–259. [Google Scholar] [CrossRef]
Halper, J.; Kjaer, M. Basic components of connective tissues and extracellular matrix: elastin, fibrillin, fibulins, fibrinogen, fibronectin, laminin, tenascins and thrombospondins. Adv Exp Med Biol 2014, 802, 31–47. [Google Scholar]
Docheva, D.; Popov, C.; Alberton, P.; Aszodi, A. Integrin signaling in skeletal development and function. Birth Defects Res C Embryo Today 2014, 102, 13–36. [Google Scholar] [CrossRef] [PubMed]
Mousavizadeh, R.; Hojabrpour, P.; Eltit, F.; McDonald, P.C.; Dedhar, S.; McCormack, R.G.; Duronio, V.; Jafarnejad, S.M.; Scott, A. beta1 integrin, ILK and mTOR regulate collagen synthesis in mechanically loaded tendon cells. Sci Rep 2020, 10, 12644. [Google Scholar] [CrossRef] [PubMed]
Maeda, E.; Ye, S.; Wang, W.; Bader, D.L.; Knight, M.M.; Lee, D.A. Gap junction permeability between tenocytes within tendon fascicles is suppressed by tensile loading. Biomech Model Mechanobiol 2012, 11, 439–447. [Google Scholar] [CrossRef] [PubMed]
Disser, N.P.; Sugg, K.B.; Talarek, J.R.; Sarver, D.C.; Rourke, B.J.; Mendias, C.L. Insulin-like growth factor 1 signaling in tenocytes is required for adult tendon growth. FASEB J 2019, 33, 12680–12695. [Google Scholar] [CrossRef] [PubMed]
Hou, Y.; Mao, Z.; Wei, X.; Lin, L.; Chen, L.; Wang, H.; Fu, X.; Zhang, J.; Yu, C. The roles of TGF-beta1 gene transfer on collagen formation during Achilles tendon healing. Biochem Biophys Res Commun 2009, 383, 235–239. [Google Scholar] [CrossRef] [PubMed]
Wang, X.; Liu, S.; Yu, T.; An, S.; Deng, R.; Tan, X.; Crane, J.; Zhang, W.; Pan, D.; Wan, M.; et al. Inhibition of Integrin alphavbeta6 Activation of TGF-beta Attenuates Tendinopathy. Adv Sci (Weinh) 2022, 9, e2104469. [Google Scholar] [CrossRef]
Liu, X.; Zhu, B.; Li, Y.; Liu, X.; Guo, S.; Wang, C.; Li, S.; Wang, D. The Role of Vascular Endothelial Growth Factor in Tendon Healing. Front Physiol 2021, 12, 766080. [Google Scholar] [CrossRef]
Liu, W.; Watson, S.S.; Lan, Y.; Keene, D.R.; Ovitt, C.E.; Liu, H.; Schweitzer, R.; Jiang, R. The atypical homeodomain transcription factor Mohawk controls tendon morphogenesis. Mol Cell Biol 2010, 30, 4797–4807. [Google Scholar] [CrossRef]
Ito, Y.; Toriuchi, N.; Yoshitaka, T.; Ueno-Kudoh, H.; Sato, T.; Yokoyama, S.; Nishida, K.; Akimoto, T.; Takahashi, M.; Miyaki, S.; et al. The Mohawk homeobox gene is a critical regulator of tendon differentiation. Proc Natl Acad Sci U S A 2010, 107, 10538–10542. [Google Scholar] [CrossRef]
Kayama, T.; Mori, M.; Ito, Y.; Matsushima, T.; Nakamichi, R.; Suzuki, H.; Ichinose, S.; Saito, M.; Marumo, K.; Asahara, H. Gtf2ird1-Dependent Mohawk Expression Regulates Mechanosensing Properties of the Tendon. Mol Cell Biol 2016, 36, 1297–1309. [Google Scholar] [CrossRef]
Riley, G. Tendinopathy--from basic science to treatment. Nat Clin Pract Rheumatol 2008, 4, 82–89. [Google Scholar] [CrossRef] [PubMed]
Sbardella, D.; Tundo, G.R.; Fasciglione, G.F.; Gioia, M.; Bisicchia, S.; Gasbarra, E.; Ippolito, E.; Tarantino, U.; Coletta, M.; Marini, S. Role of metalloproteinases in tendon pathophysiology. Mini Rev Med Chem 2014, 14, 978–987. [Google Scholar] [CrossRef] [PubMed]
Galloway, W.A.; Murphy, G.; Sandy, J.D.; Gavrilovic, J.; Cawston, T.E.; Reynolds, J.J. Purification and characterization of a rabbit bone metalloproteinase that degrades proteoglycan and other connective-tissue components. Biochem J 1983, 209, 741–752. [Google Scholar] [CrossRef] [PubMed]
Nicholson, R.; Murphy, G.; Breathnach, R. Human and rat malignant-tumor-associated mRNAs encode stromelysin-like metalloproteinases. Biochemistry 1989, 28, 5195–5203. [Google Scholar] [CrossRef] [PubMed]
Fernandes, R.J.; Hirohata, S.; Engle, J.M.; Colige, A.; Cohn, D.H.; Eyre, D.R.; Apte, S.S. Procollagen II amino propeptide processing by ADAMTS-3. Insights on dermatosparaxis. J Biol Chem 2001, 276, 31502–31509. [Google Scholar] [CrossRef]
Colige, A.; Ruggiero, F.; Vandenberghe, I.; Dubail, J.; Kesteloot, F.; Van Beeumen, J.; Beschin, A.; Brys, L.; Lapiere, C.M.; Nusgens, B. Domains and maturation processes that regulate the activity of ADAMTS-2, a metalloproteinase cleaving the aminopropeptide of fibrillar procollagens types I-III and V. J Biol Chem 2005, 280, 34397–34408. [Google Scholar] [CrossRef]
Vadon-Le Goff, S.; Kronenberg, D.; Bourhis, J.M.; Bijakowski, C.; Raynal, N.; Ruggiero, F.; Farndale, R.W.; Stocker, W.; Hulmes, D.J.; Moali, C. Procollagen C-proteinase enhancer stimulates procollagen processing by binding to the C-propeptide region only. J Biol Chem 2011, 286, 38932–38938. [Google Scholar] [CrossRef]
Hubmacher, D.; Wang, L.W.; Mecham, R.P.; Reinhardt, D.P.; Apte, S.S. Adamtsl2 deletion results in bronchial fibrillin microfibril accumulation and bronchial epithelial dysplasia--a novel mouse model providing insights into geleophysic dysplasia. Dis Model Mech 2015, 8, 487–499. [Google Scholar] [CrossRef]
Jones, G.C.; Corps, A.N.; Pennington, C.J.; Clark, I.M.; Edwards, D.R.; Bradley, M.M.; Hazleman, B.L.; Riley, G.P. Expression profiling of metalloproteinases and tissue inhibitors of metalloproteinases in normal and degenerate human achilles tendon. Arthritis Rheum 2006, 54, 832–842. [Google Scholar] [CrossRef]
Jelinsky, S.A.; Rodeo, S.A.; Li, J.; Gulotta, L.V.; Archambault, J.M.; Seeherman, H.J. Regulation of gene expression in human tendinopathy. BMC Musculoskelet Disord 2011, 12, 86. [Google Scholar] [CrossRef]
Ireland, D.; Harrall, R.; Curry, V.; Holloway, G.; Hackney, R.; Hazleman, B.; Riley, G. Multiple changes in gene expression in chronic human Achilles tendinopathy. Matrix Biol 2001, 20, 159–169. [Google Scholar] [CrossRef] [PubMed]
Dean, B.J.; Franklin, S.L.; Carr, A.J. A systematic review of the histological and molecular changes in rotator cuff disease. Bone Joint Res 2012, 1, 158–166. [Google Scholar] [CrossRef] [PubMed]
Arvind, V.; Huang, A.H. Reparative and Maladaptive Inflammation in Tendon Healing. Front Bioeng Biotechnol 2021, 9, 719047. [Google Scholar] [CrossRef]
Ellis, I.; Schnabel, L.V.; Berglund, A.K. Defining the Profile: Characterizing Cytokines in Tendon Injury to Improve Clinical Therapy. J Immunol Regen Med 2022, 16. [Google Scholar] [CrossRef]
Jomaa, G.; Kwan, C.K.; Fu, S.C.; Ling, S.K.; Chan, K.M.; Yung, P.S.; Rolf, C. A systematic review of inflammatory cells and markers in human tendinopathy. BMC Musculoskelet Disord 2020, 21, 78. [Google Scholar] [CrossRef]
Russo, V.; El Khatib, M.; Prencipe, G.; Citeroni, M.R.; Faydaver, M.; Mauro, A.; Berardinelli, P.; Cervero-Varona, A.; Haidar-Montes, A.A.; Turriani, M.; et al. Tendon Immune Regeneration: Insights on the Synergetic Role of Stem and Immune Cells during Tendon Regeneration. Cells 2022, 11. [Google Scholar] [CrossRef]
Langberg, H.; Rosendal, L.; Kjaer, M. Training-induced changes in peritendinous type I collagen turnover determined by microdialysis in humans. J Physiol 2001, 534, 297–302. [Google Scholar] [CrossRef]
Arampatzis, A.; Karamanidis, K.; Albracht, K. Adaptational responses of the human Achilles tendon by modulation of the applied cyclic strain magnitude. J Exp Biol 2007, 210, 2743–2753. [Google Scholar] [CrossRef]
Couppe, C.; Kongsgaard, M.; Aagaard, P.; Hansen, P.; Bojsen-Moller, J.; Kjaer, M.; Magnusson, S.P. Habitual loading results in tendon hypertrophy and increased stiffness of the human patellar tendon. J Appl Physiol (1985) 2008, 105, 805–810. [Google Scholar] [CrossRef]
Arampatzis, A.; Peper, A.; Bierbaum, S.; Albracht, K. Plasticity of human Achilles tendon mechanical and morphological properties in response to cyclic strain. J Biomech 2010, 43, 7–7. [Google Scholar] [CrossRef]
Seynnes, O.R.; Erskine, R.M.; Maganaris, C.N.; Longo, S.; Simoneau, E.M.; Grosset, J.F.; Narici, M.V. Training-induced changes in structural and mechanical properties of the patellar tendon are related to muscle hypertrophy but not to strength gains. J Appl Physiol (1985) 2009, 107, 523–530. [Google Scholar] [CrossRef]
Bohm, S.; Mersmann, F.; Tettke, M.; Kraft, M.; Arampatzis, A. Human Achilles tendon plasticity in response to cyclic strain: effect of rate and duration. J Exp Biol 2014, 217, 4010–4017. [Google Scholar] [CrossRef]
Bohm, S.; Mersmann, F.; Arampatzis, A. Human tendon adaptation in response to mechanical loading: a systematic review and meta-analysis of exercise intervention studies on healthy adults. Sports Med Open 2015, 1, 7. [Google Scholar] [CrossRef] [PubMed]
Rosager, S.; Aagaard, P.; Dyhre-Poulsen, P.; Neergaard, K.; Kjaer, M.; Magnusson, S.P. Load-displacement properties of the human triceps surae aponeurosis and tendon in runners and non-runners. Scand J Med Sci Sports 2002, 12, 90–98. [Google Scholar] [CrossRef] [PubMed]
Sponbeck, J.K.; Perkins, C.L.; Berg, M.J.; Rigby, J.H. Achilles Tendon Cross Sectional Area Changes Over a Division I NCAA Cross Country Season. Int J Exerc Sci 2017, 10, 1226–1234. [Google Scholar] [PubMed]
Kongsgaard, M.; Reitelseder, S.; Pedersen, T.G.; Holm, L.; Aagaard, P.; Kjaer, M.; Magnusson, S.P. Region specific patellar tendon hypertrophy in humans following resistance training. Acta Physiol (Oxf) 2007, 191, 111–121. [Google Scholar] [CrossRef]

Figure 1. RNA-Seq analysis on semitendinosus tendons from 12 diplegic CP (D), 12 tetraplegic CP (T), and 5 control (C) patients. A) Principal components analysis (PCA) based on transcriptional profiles with ellipses of confidence intervals set at 0.9 for each group. B) Heat map of unsupervised hierarchical clustering of 1223 differentially expressed protein-coding genes (FDR ≤ 0.1; |log2FC| ≥ 0.4). C) Venn diagram of the pairwise differential expression (DE) analysis (FDR ≤ 0.1; |log2FC| ≥ 0.4).

Figure 2. EnrichR analysis for upregulated and downregulated genes in semitendinosus tendons from tetraplegic CP (T) patients vs. control (C) patients (FDR ≤ 0.1; |log2FC| ≥ 0.4). A) Gene ontology (GO) analysis for biological processes and B) Reactome analysis (adjusted p value ≤ 0.05). For upregulated genes, the 20 most significantly enriched GO terms and pathways are shown.

Figure 3. Differentially expressed genes in semitendinosus tendons from tetraplegic CP (T) patients vs. control (C) patients involved in ECM matrix organization (A), membrane trafficking (B), cytokine signaling (C), and apoptotic processes (D) (FDR ≤ 0.1; |log2FC| ≥ 0.4).

Table 1. Patient characteristics

Patient	Gender	Age	Type
1	F	15.8	C
2	M	26.4	C
3	M	14.6	C
4	M	16.1	C
5	M	16.8	C
6	M	21.1	D
7	M	15.5	D
8	M	17.1	D
9	M	15.0	D
10	F	5.6	D
11	M	16.5	D
12	M	10.5	D
13	M	11.0	D
14	M	15.7	D
15	F	11.4	D
16	M	17.1	D
17	M	10.1	D
18	M	11.5	T
19	F	11.0	T
20	F	14.5	T
21	M	16.2	T
22	M	4.9	T
23	F	14.0	T
24	M	14.9	T
25	F	12.8	T
26	F	11.6	T
27	M	13.3	T
28	M	9.2	T
29	M	16.2	T
30	M	13.0	T

M, male; F, female; C, control, D, diplegic; T, tetraplegic

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer