Skip to main content

Transcription and proteome changes involved in re-innervation muscle following nerve crush in rats

Abstract

Severe peripheral nerve injury leads to the irreparable disruption of nerve fibers. This leads to disruption of synapses with the designated muscle, which consequently go through progressive atrophy and damage of muscle function. The molecular mechanism that underlies the re-innervation process has yet to be evaluated using proteomics or transcriptomics. In the present study, multi-dimensional data were therefore integrated with transcriptome and proteome profiles in order to investigate the mechanism of re-innervation in muscles. Two simulated nerve injury muscle models in the rat tibial nerve were compared: the nerve was either cut (denervated, DN group) or crushed but with the nerve sheath intact (re-innervated, RN group). The control group had a preserved and intact tibial nerve. At 4 weeks, the RN group showed better tibial nerve function and recovery of muscle atrophy compared to the DN group. As the high expression of Myh3, Postn, Col6a1 and Cfi, the RN group demonstrated superior re-innervation as well. Both differentially expressed genes (DEGs) and proteins (DEPs) were enriched in the peroxisome proliferator-activated receptors (PPARs) signaling pathway, as well as the energy metabolism. This study provides basic information regarding DEGs and DEPs during re-innervation-induced muscle atrophy. Furthermore, the crucial genes and proteins can be detected as possible treatment targets in the future.

Peer Review reports

Background

Skeletal muscle denervation is a common disease caused by trauma, inflammation, or disease. It is well-known that denervated muscle undergoes atrophic changes, involving disruption of sarcomere and myofibrilla, and the degree of atrophy increases with the duration of denervation [1]. In order to analyze the gene expression profiles induced by denervated muscle atrophy, microarray analysis has been used in the previous survey, as well as the RNA sequencing (RNA seq) [2]. Several proteomics studies of denervated skeletal muscle have also revealed biomarkers and pathways involved in denervation-induced muscle atrophy [3]. In spite of previous understanding of the mechanism for denervated nerves, the re-innervation response to crushed nerves has not been well characterized.

The effects of axotomy on motoneurons, nerves and muscles implicate various morphological, biochemical and physiological changes. After a nerve is transected, the damage axons will regenerate following their unique pathway, and eventually re-innervate the muscles. As previously described, nerve transection will redirect regenerated axons to inappropriate pathways [4]. In contrast to axotomy, the experimental technique of crushing maintains the stability of nerve sheaths and basement membranes. Despite the fact that advances have been made in clarifying the mechanisms of axon and muscle regeneration, the underlying mechanism for re-innervation has yet to be evaluated using a proteomic or transcriptomic approach.

Therefore, the current study aims to investigate possible functional mechanisms by identifying differentially expressed mRNA and proteins in re-innervated muscle over a 4-week period following injury. Furthermore, possible treatment targets might be identified by using bioinformatic analysis, including Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The 4-week study period covers the time when neurogenic muscular atrophy is still reversible and nerve re-innervation and myogenic regeneration can still occur [5, 6]. This research brings new insight into molecular mechanism of nerve re-innervation, while identifying specific targets and pathways for the novel treatment of nerve re-innervation.

Results

General observation

Daily observations of animals in each experimental group revealed normal behavior and no systemic or local inflammation or post-operative complications.

Tibial function index (TFI) analysis

TFI analysis revealed that the TFI values decreased after denervation and crush, and repair increasingly with time. The TFI values in the RN and DN group were significantly lower than in the control group at every time point, respectively (P < 0.001, Fig. 1a). The TFI values in the RN group were significantly higher than in the DN group at every time point (P < 0.001, Fig. 1a).

Fig. 1
figure 1

RN group achieved a higher TFI value and better protection for the gastrocnemius muscle. a The tibial functional index (TFI) was determined after analysis of the walking trajectory. The higher TFI values represent better functional recovery. #RN group versus control group, *DN group versus control group. Results statistically significant at P < 0.05 (*, #), P < 0.01 (**, ##) and P < 0.001 (***, ###). b The ratio of muscle on the surgical side compared to the healthy side. * Differences were statistically significant at P < 0.05. By 4 weeks, RN group revealed a well recovery of muscle compared with the DN (P < 0.01)

Muscle wet weight

The gastrocnemius muscle wet weight was scaled at 4 weeks (Fig. 1b). To deeply know the muscle mass, we compared the ratio of muscle samples between surgical and healthy side. The muscle mass of the control group is 100% retained. By 4 weeks, DN group revealed a poor return of muscle recording only 86.3% ± 7% (P < 0.05). In RN group, it was 96.2% ± 0.6% which was well recovered compared with the DN (P < 0.05).

Differentially expressed genes (DEGs) and proteins (DEPs)

Untargeted transcriptome and proteome analyses of muscles were performed to explore repair mechanisms in the DN and RN groups (Fig. 2a). A total of 567 genes were identified as differentially expressed genes (DEGs) at the transcriptome level, out of which 151 were altered in DN mice relative to control mice, while 441 DEGs were altered through re-innervating (Fig. 2a). Upon comparing the DN group with the RN group, 344 DEGs were found to be significantly up-regulated, while 97 DEGs were significantly down-regulated. Compared with the control group, 19 genes were down-regulated in DN group, while restored through re-innervating. In the same way, 6 genes were up-regulated in DN group and restored through re-innervating. These results were visualized in the hierarchically clustered heatmap (Fig. 2b) and volcano plot (Fig. 2c). Heatmap representation of hierarchical clustering of DEGs that were significantly altered by RN.

Fig. 2
figure 2

Signatures of DN and RN group at transcript and protein levels. a Venn-diagrams showing the number of differentially expressed genes (DEGs) that were down-regulated in the DN group (compared with the control group) and up-regulated in the RN group (compared with the DN group), and DEGs that were up-regulated in DN rats and down-regulated in RN rats. b Heatmap representation of hierarchical clustering of DEGs that were significantly altered in all three groups. c Volcano plot revealing the distribution of DEGs in the RN group compared to the DN group. d Venn-diagram revealing the number of differentially expressed proteins (DEPs) in DN group compared to control group, and DEPs in DN group compared to RN group. e Heatmap representation of hierarchical clustering of DEPs showing the proteins altered in all three groups. f Volcano plot showing distribution of the DEPs in RN group compared to the DN group

At the proteome level, a total of 385 proteins were differentially expressed. The Venn diagram in Fig. 2d shows the number of differentially expressed proteins (DEPs) in all 3 groups. Among the 223 proteins that were significantly affected by RN, 70 were up-regulated and 153 were down-regulated. Compared with control group, 10 proteins were down-regulated in DN group, and then restored through re-innervating, while 69 proteins were up-regulated and restored through re-innervating. The corresponding heatmap and volcano plot are shown in Fig. 2e, f.

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways for the Differentially expressed genes (DEGs) and proteins (DEPs)

Pathway enrichment was analyzed against KEGG database to obtain the potential functions of the DEGs and DEPs. As the results shown in Figs. 3a, b, c, d, DEGs and DEPs were all enriched in pathways involving in PPAR signaling pathway, while DEPs were significantly enriched in pathways related to energy metabolism. Among 126 DEPs distributed in the profile 5 of the short time-series expression miner (STEM) analysis (Fig. 3f), 55 DEPs were related to energy metabolism and suggested RN group a positive effect on energy metabolism comparing with DN group.

Fig. 3
figure 3

Network analysis of transcriptome and proteome. a, b, c, d Analysis of KEGG pathway enrichment in DEGs (a, b) and DEPs (c, d). Rich factor refers to the ratio of the number of DEGs (or DEPs) in the pathway entry to the total number of genes (or proteins) in that pathway entry. Pathways were considered significantly enriched when P < 0.05. e, f Short time-series expression miner (STEM) clustering of DEGs (Fig. e) and DEPs (Fig. f). Each square represents the expression trend of a gene (or protein) in the control, DN and RN groups. The black polygonal line in the square indicates the overall expression profile of all DEGs (or DEPs) in the corresponding set. The profile ID and the number of DEGs and DEPs in each profile are shown on the top of the square. Colored squares are significantly different (P < 0.05). DEPs involved in the energy metabolic pathway were mainly distributed in the expression patterns of profiles 5 in STEM analysis

Analysis combined transcriptome with proteome data

Furthermore, we combined the transcriptome and proteome data to identify genes that were altered in both omics’ layers. The Venn diagram in Fig. 4a shows the number of co-regulated DEGs/ DEPs, and diff-regulated DEGs/ DEPs between the DN group and control group. In the transcriptome level, 209 DEGs were up-or-down-regulated, while 1344 DEPs were up-or-down-regulated in the proteome level. Combined with transcriptome and proteome data, 53 were co-regulated while 61 were diff-regulated. Compared with RN group, 96 were co-regulated while 71 were diff-regulated (Fig. 4d). Of the 5 genes were found to be up-regulated in DN group (Fig. 4c, Q3), out of a total of 2534 genes co-regulated (Fig. 3b, Q3). In RN group, 4 genes differed significantly at both the protein and mRNA levels compared with DN group (Fig. 4f, Q3). These included Myh3, Postn, Col6a1, Thbs4, Cfi and so on, which are widely known to be related with muscle development, energy metabolism and neuronal differentiation.

Fig. 4
figure 4

Correlational analysis of transcriptome and proteome. a Venn-diagrams showing the number of co-regulated DEGs/ DEPs, and diff-regulated DEGs/ DEPs between the DN group and control group. b DEGs and DEPs that were altered in both omics layers between the DN group and control group. c Correlational analysis of DEGs and DEPs in the DN group relative to the control group. The omics datasets were merged based on the gene symbol. Red dots in the Q3 or Q7 quadrants represent co-regulated entities that were consistently up-regulated or down-regulated. Green dots in the Q2 and Q8 quadrants represent genes that were significantly altered in the transcriptome, while blue dots in the Q4 and Q6 quadrants represent significant genes in the proteome (Fig. c, Q3). d Venn-diagram revealing the number of co-regulated DEGs/ DEPs, and diff-regulated DEGs/ DEPs between the DN group and RN group. e DEGs and DEPs that were altered in both omics layers between the RN group and DN group. f Correlational analysis of DEGs and DEPs in the RN group relative to the DN group

Discussion

Several studies support the notion that longer de-innervated periods cause reduced re-innervation and ultimately reduced the recovery of muscle function [7]. Significantly, the re-innervation of peripheral nerves is known to influence the recovery of denervated muscle [8]. These observations highlight the need to re-innervating muscles during the early recovery and before unrecoverable neural impairment occurs. Hence, there is considerable interest in identifying the molecular mechanisms of repair using nerve injury models so as to avoid irreversible outcomes post-operatively. Gene expression studies performed with various muscle injury models have reported that certain genes may be critical for achieving a good clinical result [9, 10]. To date, however, the importance of nerve integrity in the repair and recovery of nerve injury is not known. In the present work, the tibial nerve was crushed near the sciatic nerve branch, conducting a re-innervated model and studying gastrocnemius muscle specifically. Moreover, we conducted transcriptomic and proteomic analyses of muscle following nerve crush and transection in order to gain insights into the repair mechanisms following these injuries.

In the current study, transcriptomic and proteomic analysis of muscles were evaluated at 4 weeks after the nerve injury. The different degrees of nerve injury applied to the DN and RN groups appears to result in the different expression profiles of genes and proteins, as well as the muscle repair process. We proved that RN group achieved better functional results concerning the peripheral nerve injury model of rats than DN group. The tibial nerve function of the RN group recovered better with a higher TFI value and better protection for the gastrocnemius muscle. Transection of nerve caused a greater degree of axonal destruction, which required greater recovery and re-innervation at the lesion site. In contrast with transection, nerve crush preserved the nerve sheath conduit, so as to allow the axonal regeneration and ultimately evoke a different genes profile. Apart from the Venn-diagrams, volcano plot in Fig. 2c was also demonstrated the different gene profiles between DN and RN groups at 4 weeks post-repair. Furthermore, pathway enrichment analysis based on the KEGG database was performed to evaluate the potential function of DEGs and DEPs. As shown in Fig. 3c, DEPs were significantly enriched in pathways related to energy metabolism, such as fatty acid metabolism, degradation and elongation. Interestingly, the DEPs involved in energy metabolic pathway were mainly distributed in the expression patterns of profiles 5 in STEM analysis (Fig. 3f), which suggest that RN group had a positive effect on energy metabolism comparing to DN group.

When the peripheral nerve is damaged, the distal axons of the injured site degenerates, then Schwann cells of peripheral nerve are able to proliferate and re-differentiate, thus providing an environment for axonal and nerve regeneration [11]. As Zhang et al. reported, Peroxisome proliferator-activated receptors-γ (PPAR-γ) expression in SCs increases during the inflammation process and this may promote their activation as immune cells [12]. In the current study, the PPAR pathway was found to be one of the major signaling pathways according to the multi-omics analyses. They regulate genes that are critical for cellular differentiation, lipid and metabolic homeostasis, and energy balance [13, 14]. PPAR-γ plays a key role in muscle cell function and myogenesis, which expressed in the skeletal muscle of humans [15]. As a nuclear receptor, it can activate cellular metabolism to promote cell growth and differentiation [16, 17]. Given these observations, drugs that activate PPAR-γ signaling pathway may therefore be effective in reducing and preventing muscle atrophy.

Cellular respiration is a series of metabolic reactions that convert biochemical energy to ATP and then release the waste products [18, 19]. Skeletal muscle cells efficiently converted energy through several metabolic pathways to maintain continuous muscle contraction. In denervated muscle, the process by which mitochondria produce ATP is damaged. Consistent with earlier reports [20], we found lower expression levels for almost all metabolic pathways after nerve denervation. In pathological conditions, the ability to regulate energy metabolism of skeletal muscle is disrupted. This may explain why many pathways associated with energy metabolism were found to be enriched in the current research. Furthermore, the levels of ATP increase rapamycin (mTOR) activity, and the mammalian target of mTOR pathway is the latest intracellular fuel-sensing mechanism to be implicated in the regulation of energy balance [21,22,23]. In organism, mTOR was a key regulator in maintaining skeletal muscle mass which regulates protein anabolism and catabolism [24]. We have previously proven that negative regulation of mTORC1 signaling could protect the denervated muscle atrophy and improved the functional outcome [25,26,27,28]. The essence of denervated muscle atrophy results from more significant muscle protein degradation than synthesis [29, 30]. After denervation, mTORC1 signaling plays more important roles in protein degradation than in protein synthesis and leads to muscle atrophy [24]. Recent studies have also demonstrated that mTORC1 activation increases both oxygen consumption and mitochondrial number [31], it might have a significant role in the coordination of cell growth and energy metabolism. Combined with previous studies, the changes in genes and proteins involved in energy metabolism we investigated in the present study, may be associated with the mTOR pathway, leading to muscle atrophy and thus reduced muscle function.

In the rat model used here, a significant up-regulation of thrombospondin 4 (Thbs4) in RN muscles was observed in the 4-week post-surgery period, which has important functions for neurogenesis in early post-natal and adult rodent brain [32]. Thbs4 could enhance the vesicular trafficking of membrane attachment complexes to stabilize the sarcolemma in myofibers [33]. On the other hand, Thbs4 homozygous knockout mice (Thbs4-KO) resulted in abnormal glial scar formation after injury, and significantly increased microvascular haemorrhage into the brain parenchyma [32]. Loss of Thbs4 expression led to dystrophic changes in the muscular dystrophic mouse models, while overexpression protected muscle and alleviated dystrophic disease.

Other nerve growth factors with increased expression observed in the RN group included Myh3, as well as Postn, Col6a1, and Cfi. Myosin heavy chain 3 (Myh3) gene is an important component of myosin, which is related to the traction and sliding of muscle [34]. It is highly induced during in vitro muscle cell differentiation, and its expression is markedly increased in muscular dystrophy, where it is thought to represent a marker of muscle regeneration [35]. As an important marker gene for muscle regeneration, Myh3 was also found high expression in the RN group in our study, and its expression correlated with the time when most muscle repair occurs [36]. In addition, Myh3 had also been reported to express during fetal development and muscle regeneration [37]. Periostin (Postn) is an extracellular matrix protein secreted by osteoblast and osteoblast-like cell line, which can promote proliferation and differentiation of osteoblast and its precursors on the periosteum [38]. It can not only promote the reconstruction of damaged blood vessels [39], but also promote the hyperplasia and fibrosis of inflammatory tissue [40]. Basing on above evidences, this may be related to the process of muscle repair and nerve regeneration in the RN group.

Collagen type VI alpha 1 chain (Col6α1) is widely distributed in the extracellular matrix of systemic tissues including skeletal muscle. Deficiency of Col6α1 has been shown to cause mitochondrial defects (25–30) and impaired muscle regeneration [41]. As Sakurai et al. showed [42], Col6α1 supplementation promoted muscle regeneration and maturation in Col6a1-KO mice. Thus, the high expression of Col6α1 in RN group might relate to the better muscle recovery in the present study. As for the complement factor I (Cfi), local complement expression by skeletal muscle in vivo may be implicated in some muscular inflammatory or pathological processes. As shown in the above studies, these genes are involved in regulating the muscle development or regeneration in vitro, however, the specific mechanism needs to be further studied and allow them to be the possible test candidates for therapeutic application.

It should be mentioned that there are some limitations in this study. The rats used in the experiments were female at about 8 weeks, which are between the puberty stage and breeding age. The phase of estrous cycle may have an influence on the secretion of various hormones that regulate muscle metabolism. After crushing nerve injury, temporary activity limitations may also affect food intake, which may in turn affect gene and muscle protein expression. However, the rats used in each group are the same age, which can ensure the consistency of each group and minimize the influence of hormonal profile between each group. Another drawback of our study was genes and related pathways were not investigated. The detection of candidate genes and proteins are not enough to support the whole survey for obtaining more comprehensive information. Perhaps, hormone secretion levels of rats and functional analyses of key genes/ proteins could be included in the future study to obtain more comprehensive information.

In spite of the mechanisms of muscle atrophy are still not fully understood, new approaches are now available to comprehensively observe changes in gene and protein expression during skeletal muscle atrophy caused by crush or transection. This has prompted investigations into novel therapeutic and biological solutions to enhance neural repair at an earlier phase. We identified 441 DEGs and 223 DEPs that were specific to the way of nerve injury. Furthermore, several possible candidate genes were identified in the re-innervated nerve, including Thbs4, Myh3, Postn, Col6a1 and Cfi. In particular, the PPAR signaling and energy metabolism pathways seems to play a vital role and warrant further study. This study provides basic information on the changes that occur in gene and protein expression during re-innervation-induced muscle atrophy and identified several key proteins as possible treatment targets. The results form the basis for further research aimed at the prevention and therapy of muscle atrophy.

Methods and materials

Animals

Eighteen adult female Sprague–Dawley rats (eight weeks, 200–250 g) were provided by the Experimental Animal Center of South China University of Technology and randomly divided into 3 groups (n = 6) to undergo nerve operations or sham operations. All experiments were carried out in SPF (specific pathogen free) animal house, and the animal feed was formulated based on the nutritional requirements of rats issued by the Ministry of Health.

Surgical procedures

We randomly divided the rats into three groups: the denervated (DN) group, the re-innervated (RN) group, and the control group. Standard aseptic technique was conducting to perform the surgery on the right posterior limb of each animal. Briefly, adult rats were anesthetized with the 30 mg/kg dose of sodium pentobarbital by intraperitoneal injection. All procedures were performed on animals in an unconscious stage on the same day. A longitudinal incision was made posterior to expose the sciatic nerve, and free the tibial and peroneal nerves totally. For relieving post-surgery pain, the rats were given buprenorphine (0.2 mg/kg subcutaneously) when performing the operation.

In the control group, the tibial nerve of rats was exposed clearly without being cut (Fig. 5a). In the DN group, the tibial nerve was transected at 6.5 mm proximal to the gastrocnemius muscle, thus simulating a denervated muscle. In order to prevent spontaneous re-innervation of the tibial nerve, we covered the proximal and distal nerve stumps with a silicone cap (Fig. 5b). In the RN group, the exposed tibial nerve at the same location as the transection site in the DN group was crushed three times for 10 s each using a pair of forceps, with an interval of 10 s. This model simulates a re-innervated muscle (Fig. 5c). The surgical wounds of all rats were thoroughly irrigated and sutured in layers. After the surgical incisions closed, animals were housed in controlled, maintaining a 12-h light/12-h dark cycles in the pathogen-free condition and allowed free access to water and food.

Fig. 5
figure 5

Experimental nerve injury models. The two simulated nerve injuries were the DN group for denervated nerve (nerve was transected), and the RN group for re-innervated nerve (nerve was crushed to damage the axons). In the control group, the nerve was preserved intact. All muscle lacerations were repaired with sutures

Analysis of Tibial Function Index (TFI)

TFI were calculated before and every week after surgery. After dipping the hindlimbs into black ink, the rats walked on the confined walking track (10 cm × 60 cm) covering with white paper. After the test, the footprint on the track was measured with the standard of clarity and completeness. We recorded the paired footprint length (PL, distance from heel to toe), intermediary toe spread (IT, distance from second to fourth toes), and toe spread (TS, distance from first to fifth toes) of each animal’s hind limb. The lengths of normal control left hind leg were defined as NPL, NIT, and NTS, while the experimental right legs were EPL, EIT, and ETS. The TFI analysis was performed following the formula of Bain-Mackinnon-Hunter [43]:

$$\mathrm{TFI }=104.4 \left(\left[\mathrm{ETS}-\mathrm{NTS}\right]/\mathrm{NTS}\right)- 37.2 \left(\left[\mathrm{EPL}-\mathrm{NPL}\right]/\mathrm{NPL}\right)+ 45.6 \left(\left[\mathrm{EIT}-\mathrm{NIT}\right]/\mathrm{NIT}\right) -8.8$$

Muscle weight

Four weeks after the surgery, the animals were sacrificed by the lethal dose of ether. Three operated muscle samples from each group were randomly selected to peel the gastrocnemius muscle and tibial nerve completely. The muscle wet weight was scaled analytical balance (BSA124S, sartorius). The muscle samples were snap-frozen in liquid nitrogen for more than 15 min, and stored at -80℃ until bioinformatic analysis 3 days later.

Transcriptomic analysis

RNA extraction, library construction and sequencing

Trizol reagent kit (Invitrogen, Carlsbad, CA, USA) was used to extract all RNA. Oligo(dT) beads were used to enrich eukaryotic mRNA, while Ribo-ZeroTM Magnetic Kit (Epicentre, Madison, WI, USA) was used to remove rRNA and enrich prokaryotic mRNA. The enriched mRNA was reverse transcribed into cDNA using primers. DNA polymerase I, RNase H, dNTP and buffer were used to synthesize second-strand cDNA, and purified using PCR extraction kit (Qiagen, Venlo, The Netherlands). The sequencing platform of Genedenovo Biotechnology Co., Ltd (Guangzhou, China) is responsible for constructing the library and sequencing it.

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis

Differential expression analysis of mRNA between two groups was performed using the R package DESeq2, and between two samples by edgeR. Statistically significant differentially expressed genes were detected if parameter of false discovery rate (FDR) set at < 0.05 and an absolute fold change ≥ 2. The KEGG analysis was performed using the KEGG pathway database (Release 94, https://www.genome.jp/kegg/pathway.html) to analyze the selected DEGs at the functional level [44,45,46].

Proteomic analysis

Protein digestion

The gastrocnemius muscle samples were put into the lysis buffer (1 mg/ml protease inhibitor, 7 M urea, 2% SDS), The protein concentration of the supernatant was measured using BCA Protein Assay Kit. 50 μg proteins were suspended in 50μL solution, 1uL 1 M dithiotreitol was added for reduction at 55 °C for 1 h, and 5uL 20 mM iodoacetamide was added for alkylation at 37 °C for 1 h. Then 300uL acetone was chilled and precipitated the sample overnight at -20℃. After two washes with cold acetone, the precipitate was resuspended in 50 mM ammonium bicarbonate. Finally, sequence-grade modified trypsin (Promega, Madison, WI) was used to digest the proteins at a substrate/enzyme ratio of 50:1 (w/w) at 37 °C for 16 h.

High PH Reverse Phase Separation

By using high pH separation, the peptide mixture was fractionated on an Ultimate 3000 UPLC system (Thermo Fisher Scientific, Waltham, MA, USA), then analyzed on EASY-nLC 1200 system (Thermo Fisher Scientific, MA, USA).

Data analysis

Spectronaut X (Biognosys AG, Switzerland) analyzed the raw Data of DIA with default parameters. DEPs with Absolute AVG log2 ratio > 0.58 and q < 0.05 will be filtered after Student’s T-test. KEGG database were annotated to obtain function of proteins, while GO pathways were detected within DEPs with q < 0.05.

Statistical analysis

Continuous data were presented as mean ± standard deviation or median with interquartile range. Two-tailed tests were used, and differences were considered significant if p < 0.05. For experiments on more than two groups within one categorical variable, ANOVA analysis was performed followed by Tukey´s post hoc test. SPSS 20.0 statistical software (Chicago, IL, USA) was used for the analysis. The data that support the findings of this study are available from the corresponding author upon reasonable request.

Availability of data and materials

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository [47] with the dataset identifier PXD033031.

Abbreviations

DN:

Denervated

RN:

Re-innervated

RNA seq:

RNA sequencing

SCs:

Schwann cells

PPARs:

Peroxisome proliferator-activated receptors

DEGs:

Differentially expressed genes

DEPs:

Differentially expressed proteins

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

Myh3 :

Myosin Heavy Chain 3

Postn :

Periostin gene

Col6a1 :

Collagen type VI alpha 1 chain

Cfi :

Complement factor I

mTOR:

Rapamycin

References

  1. Russo TL, Peviani SM, Durigan JL, Gigo-Benato D, Delfino GB, Salvini TF. Stretching and electrical stimulation reduce the accumulation of MyoD, myostatin and atrogin-1 in denervated rat skeletal muscle. J Muscle Res Cell Motil. 2010;31(1):45–57.

    CAS  PubMed  Article  Google Scholar 

  2. Weng J, Zhang P, Yin X, Jiang B. The whole transcriptome involved in denervated muscle atrophy following peripheral nerve injury. Front Mol Neurosci. 2018;11:69.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  3. Lang F, Aravamudhan S, Nolte H, Türk C, Hölper S, Müller S, Günther S, Blaauw B, Braun T, Krüger M. Dynamic changes in the mouse skeletal muscle proteome during denervation-induced atrophy. Dis Model Mech. 2017;10(7):881–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Zhou CJ, Kawabuchi M, Wang S, Liu WT, Hirata K. Age differences in morphological patterns of axonal sprouting and multiple innervation of neuromuscular junctions during muscle reinnervation following nerve crush injury. Ann Anat. 2002;184(5):461–72.

    PubMed  Article  Google Scholar 

  5. Pereira BP, Han HC, Yu Z, Tan BL, Ling Z, Thambyah A, Nathan SS. Myosin heavy chain isoform profiles remain altered at 7 months if the lacerated medial gastrocnemius is poorly reinnervated: a study in rabbits. J Orthop Res. 2010;28(6):732–8.

    CAS  PubMed  Article  Google Scholar 

  6. Shen Y, Zhang R, Xu L, Wan Q, Zhu J, Gu J, Huang Z, Ma W, Shen M, Ding F, et al. Microarray analysis of gene expression provides new insights into denervation-induced skeletal muscle atrophy. Front Physiol. 2019;10:1298.

    PubMed  PubMed Central  Article  Google Scholar 

  7. Aydin MA, Mackinnon SE, Gu XM, Kobayashi J, Kuzon WM Jr. Force deficits in skeletal muscle after delayed reinnervation. Plast Reconstr Surg. 2004;113(6):1712–8.

    PubMed  Article  Google Scholar 

  8. Carlson BM. The denervated muscle: 45 years later. Neurol Res. 2008;30(2):119–22.

    PubMed  Article  Google Scholar 

  9. Zeman RJ, Zhao J, Zhang Y, Zhao W, Wen X, Wu Y, Pan J, Bauman WA, Cardozo C. Differential skeletal muscle gene expression after upper or lower motor neuron transection. Pflugers Arch. 2009;458(3):525–35.

    CAS  PubMed  Article  Google Scholar 

  10. Warren GL, Summan M, Gao X, Chapman R, Hulderman T, Simeonova PP. Mechanisms of skeletal muscle injury and repair revealed by gene expression studies in mouse models. J Physiol. 2007;582(Pt 2):825–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. Fu SY, Gordon T. The cellular and molecular basis of peripheral nerve regeneration. Mol Neurobiol. 1997;14(1–2):67–116.

    CAS  PubMed  Article  Google Scholar 

  12. Zhang F, Liu F, Yan M, Ji H, Hu L, Li X, Qian J, He X, Zhang L, Shen A, et al. Peroxisome proliferator-activated receptor-gamma agonists suppress iNOS expression induced by LPS in rat primary Schwann cells. J Neuroimmunol. 2010;218(1–2):36–47.

    CAS  PubMed  Article  Google Scholar 

  13. Corton JC, Anderson SP, Stauber A. Central role of peroxisome proliferator-activated receptors in the actions of peroxisome proliferators. Annu Rev Pharmacol Toxicol. 2000;40:491–518.

    CAS  PubMed  Article  Google Scholar 

  14. Kersten S, Desvergne B, Wahli W. Roles of PPARs in health and disease. Nature. 2000;405(6785):421–4.

    CAS  PubMed  Article  Google Scholar 

  15. Kruszynska YT, Mukherjee R, Jow L, Dana S, Paterniti JR, Olefsky JM. Skeletal muscle peroxisome proliferator- activated receptor-gamma expression in obesity and non- insulin-dependent diabetes mellitus. J Clin Invest. 1998;101(3):543–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. Hihi AK, Michalik L, Wahli W. PPARs: transcriptional effectors of fatty acids and their derivatives. Cell Mol Life Sci. 2002;59(5):790–8.

    CAS  PubMed  Article  Google Scholar 

  17. Murphy GJ, Holder JC. PPAR-gamma agonists: therapeutic role in diabetes, inflammation and cancer. Trends Pharmacol Sci. 2000;21(12):469–74.

    CAS  PubMed  Article  Google Scholar 

  18. O’Leary MF, Hood DA. Effect of prior chronic contractile activity on mitochondrial function and apoptotic protein expression in denervated muscle. J Appl Physiol (1985). 2008;105(1):114–20.

    CAS  Article  Google Scholar 

  19. Batt J, Bain J, Goncalves J, Michalski B, Plant P, Fahnestock M, Woodgett J. Differential gene expression profiling of short and long term denervated muscle. FASEB J. 2006;20(1):115–7.

    CAS  PubMed  Article  Google Scholar 

  20. Adhihetty PJ, O’Leary MF, Chabi B, Wicks KL, Hood DA. Effect of denervation on mitochondrially mediated apoptosis in skeletal muscle. J Appl Physiol (1985). 2007;102(3):1143–51.

    CAS  Article  Google Scholar 

  21. Dennis PB, Jaeschke A, Saitoh M, Fowler B, Kozma SC, Thomas G. Mammalian TOR: a homeostatic ATP sensor. Science. 2001;294(5544):1102–5.

    CAS  PubMed  Article  Google Scholar 

  22. Cota D, Proulx K, Smith KA, Kozma SC, Thomas G, Woods SC, Seeley RJ. Hypothalamic mTOR signaling regulates food intake. Science. 2006;312(5775):927–30.

    CAS  PubMed  Article  Google Scholar 

  23. Ropelle ER, Pauli JR, Fernandes MF, Rocco SA, Marin RM, Morari J, Souza KK, Dias MM, Gomes-Marcondes MC, Gontijo JA, et al. A central role for neuronal AMP-activated protein kinase (AMPK) and mammalian target of rapamycin (mT OR) in high-protein diet-induced weight loss. Diabetes. 2008;57(3):594–605.

    CAS  PubMed  Article  Google Scholar 

  24. Yoon MS. mTOR as a Key Regulator in Maintaining Skeletal Muscle Mass. Front Physiol. 2017;8:788.

    PubMed  PubMed Central  Article  Google Scholar 

  25. Chen Y, Yuan W, Zeng X, Ma Y, Zheng Q, Lin B, Li Q. Combining reverse end-to-side neurorrhaphy with rapamycin treatment on chronically denervated muscle in rats. J Integr Neurosci. 2021;20(2):359–66.

    PubMed  Article  Google Scholar 

  26. Li QT, Zhang PX, Yin XF, Han N, Kou YH, Deng JX, Jiang BG. Functional recovery of denervated skeletal muscle with sensory or mixed nerve protection: a pilot stu dy. PLoS One. 2013;8(11):e79746.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. Li Q, Zhang P, Yin X, Han N, Kou Y, Jiang B. Early sensory protection in reverse end-to-side neurorrhaphy to improve the functional recovery of ch ronically denervated muscle in rat: a pilot study. J Neurosurg. 2014;121(2):415–22.

    CAS  PubMed  Article  Google Scholar 

  28. Li Q, Zhang P, Yin X, Jiang B. Early nerve protection with anterior interosseous nerve in modified end-to-side neurorrhaphy repairs high ulnar nerve injury: a hypothesis of a novel surgical technique. Artif Cells Nanomed Biotechnol. 2015;43(2):103–5.

    PubMed  Article  Google Scholar 

  29. Kamei Y, Miura S, Suzuki M, Kai Y, Mizukami J, Taniguchi T, Mochida K, Hata T, Matsuda J, Aburatani H, et al. Skeletal muscle FOXO1 (FKHR) transgenic mice have less skeletal muscle mass, down-regulated Type I (s low twitch/red muscle) fiber genes, and impaired glycemic control. J Biol Chem. 2004;279(39):41114–23.

    CAS  PubMed  Article  Google Scholar 

  30. Bonaldo P, Sandri M. Cellular and molecular mechanisms of muscle atrophy. Dis Model Mech. 2013;6(1):25–39.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. Finley LW, Haigis MC. The coordination of nuclear and mitochondrial communication during aging and calorie restriction. Ageing Res Rev. 2009;8(3):173–88.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. Benner EJ, Luciano D, Jo R, Abdi K, Paez-Gonzalez P, Sheng H, Warner DS, Liu C, Eroglu C, Kuo CT. Protective astrogenesis from the SVZ niche after injury is controlled by Notch modulator Thbs4. Nature. 2013;497(7449):369–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. Vanhoutte D, Schips TG, Kwong JQ, Davis J, Tjondrokoesoemo A, Brody MJ, Sargent MA, Kanisicak O, Yi H, Gao QQ, et al. Thrombospondin expression in myofibers stabilizes muscle membranes. eLife. 2016;5:e17589.

    PubMed  PubMed Central  Article  Google Scholar 

  34. Knight AE, Molloy JE. Muscle, myosin and single molecules. Essays Biochem. 2000;35:43–59.

    CAS  PubMed  Article  Google Scholar 

  35. Haslett JN, Sanoudou D, Kho AT, Bennett RR, Greenberg SA, Kohane IS, Beggs AH, Kunkel LM. Gene expression comparison of biopsies from Duchenne muscular dystrophy (DMD) and normal skeletal muscle. Proc Natl Acad Sci USA. 2002;99(23):15000–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. Chargé SB, Rudnicki MA. Cellular and molecular regulation of muscle regeneration. Physiol Rev. 2004;84(1):209–38.

    PubMed  Article  Google Scholar 

  37. Karsch-Mizrachi I, Travis M, Blau H, Leinwand LA. Expression and DNA sequence analysis of a human embryonic skeletal muscle myosin heavy chain gene. Nucleic Acids Res. 1989;17(15):6167–79.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. Litvin J, Selim AH, Montgomery MO, Lehmann K, Rico MC, Devlin H, Bednarik DP, Safadi FF. Expression and function of periostin-isoforms in bone. J Cell Biochem. 2004;92(5):1044–61.

    CAS  PubMed  Article  Google Scholar 

  39. Li P, Oparil S, Feng W, Chen YF. Hypoxia-responsive growth factors upregulate periostin and osteopontin expression via distinct signaling pathways in rat pulmonary arterial smooth muscle cells. J Appl Physiol (1985). 2004;97(4):1550–8 discussion 1549.

    CAS  Article  Google Scholar 

  40. Jia G, Erickson RW, Choy DF, Mosesova S, Wu LC, Solberg OD, Shikotra A, Carter R, Audusseau S, Hamid Q, et al. Periostin is a systemic biomarker of eosinophilic airway inflammation in asthmatic patients. J Allergy Clin Immunol. 2012;130(3):647-654.e610.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. Alexeev V, Arita M, Donahue A, Bonaldo P, Chu ML, Igoucheva O. Human adipose-derived stem cell transplantation as a potential therapy for collagen VI-related congenital muscular dystrophy. Stem Cell Res Ther. 2014;5(1):21.

    PubMed  PubMed Central  Article  Google Scholar 

  42. Takenaka-Ninagawa N, Kim J, Zhao M, Sato M, Jonouchi T, Goto M, Yoshioka CKB, Ikeda R, Harada A, Sato T, et al. Collagen-VI supplementation by cell transplantation improves muscle regeneration in Ullrich congenital muscular dystrophy model mice. Stem Cell Res Ther. 2021;12(1):446.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. Bain JR, Mackinnon SE, Hunter DA. Functional evaluation of complete sciatic, peroneal, and posterior tibial nerve lesions in the rat. Plast Reconstr Surg. 1989;83(1):129–38.

    CAS  PubMed  Article  Google Scholar 

  44. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545-d551.

    CAS  PubMed  Article  Google Scholar 

  47. Ma J, Chen T, Wu S, Yang C, Bai M, Shu K, Li K, Zhang G, Jin Z, He F, et al. iProX: an integrated proteome resource. Nucleic Acids Res. 2019;47(D1):D1211-d1217.

    PubMed  Article  Google Scholar 

Download references

Acknowledgements

No acknowledgments.

ARRIVE guidelines

Our study is reported in accordance with ARRIVE guidelines. The ARRIVE Essential 10 are following:

1. Study design

We divided all animals into 3 groups: denervated group, re-innervated group and control group. Two simulated nerve injury models in the rat tibial nerve were compared: the nerve was either cut (denervated, DN group) or crushed but with the nerve sheath intact (re-innervated, RN group). The control group had a preserved and intact tibial nerve.

2. Sample size

We provided 6 rats in each group, and a total of 18 rats in 3 groups. We determined the sample size in the previous pre-experiment. The degree of freedom (E) of variance analysis should be estimated, and the value range of E should be between 10-20.

E= sum of experimental animals in each group (3*6=18) - number of groups (3) = 15.

3. Inclusion and exclusion criteria

All rats underwent rigorous screening by the animal center to ensure their health. We selected the rats randomly to control and treatment group. Daily observations of animals in each experimental group revealed normal behavior and no systemic or local inflammation or post-operative complications, so all the included rats were analyzed.

4. Randomization

Randomization was used to allocate rats to control and treatment groups. The rats were randomly divided into 3 groups with 18 rats, which were divided into 3 groups by using random number table. First, the rats were organized into 1, 2, 3... 18, and then arbitrarily from a random number table from a row of a certain number to copy 18 numbers and corresponding. Divide each number by 3 and represent A, B, and C with the remainder 1, 2, and 3.

We do not control the confounders including the order of treatments and measurements, or animal/cage location and so as.

5. Blinding

Qingtian Li, Haotao Li, and Yijian Chen designed conducted the experiment, so they were aware of the group allocation. Bofu Lin, Shuai Wang and Zhantao Deng analyzed and interpreted the data. They were not aware of the group allocation.

6. Outcome measures

The analysis of Tibial Function Index, muscle weight measurement, and the transcriptome and proteome analysis were described in the “method” part.

7. Statistical methods

The statistical methods used for each analysis were described in the “method” part.

8. Experimental animals

Adult female Sprague-Dawley rats (eight weeks, 200–250g) were provided by the Experimental Animal Center of South China University of Technology.

9. Experimental procedures

The procedure details were described in the “Surgical Procedures” part.

10. Results

Summary/descriptive statistics for each experimental group were described in the “Result” part.

Funding

The work was supported by the Natural Science Foundation of Guangdong Province [2018A0303130206, 2021A1515011008]; National Natural Science Foundation of China [8207050645, 8217051081]; Guangdong Medical Science and Technology Research Fund [A2019150]; Traditional Chinese Medicine Bureau Foundation of Guangdong Provincial [20221008]. No benefits in any form have been received or will be received from a commercial party related directly or indirectly to the subject of this article.

Author information

Authors and Affiliations

Authors

Contributions

Qingtian Li, Haotao Li, and Yijian Chen designed and conducted the study, and approved the final version of the manuscript. Haotao Li and Qingtian Li wrote the first draft of the manuscript. Wanqiong Yuan and Qiujian Zheng modified the manuscript. Bofu Lin, Shuai Wang and Zhantao Deng analyzed and interpreted the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qingtian Li.

Ethics declarations

Ethics approval and consent to participate

Ethics approval and consent to participate All work was performed in accordance with the Guiding Opinions on the Treatment of Experimental Animals, Guangdong Province Laboratory Animal Management Regulations and the Laboratory Animal Committee (LAC) of the South China University of Technology Policy on the Humane Care and Use of Vertebrate Animals. Ethics approval was obtained from the Administration Committee of the Guangdong Provincial People’s Hospital, Guangdong Academy of Medical Sciences (Permit Number: 2018537A).

Consent for publication

The manuscript does not include details, images, or videos relating to an individual person.

Competing interests

The study was financially supported by institutional funds only. All authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

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

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, H., Yuan, W., Chen, Y. et al. Transcription and proteome changes involved in re-innervation muscle following nerve crush in rats. BMC Genomics 23, 666 (2022). https://doi.org/10.1186/s12864-022-08895-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08895-w

Keywords

  • Re-innervation
  • Denervation
  • Peripheral nerve injury
  • Skeletal muscle atrophy
  • Transcriptomics
  • Proteomics