Genetic differences linked to movement, depression, and pain sensitivity after spinal cord injury in three mouse types

Young adult male mice of the A/J, BALB/c, and C57BL/6 inbred strains were purchased from Jackson Labs at 9 to 12 weeks of age. All animals were housed in the University of Maryland School of Medicine animal facility in a controlled environment with a 12 h (6 a.m. to 6 p.m.) light/dark cycle and ad libitum access to food and water at 4 to 5 mice per cage. Following isoflurane anesthesia, mice were subjected to a laminectomy at the T9 to T11 region, and the spinal column was placed onto bilateral steel clamps for stabilization. The mice received T10 spinal cord contusions with 35 (mild injury, m) or 60 (moderate/severe injury, mo/se) kDyne force using the Infinite Horizon Spinal Cord Impactor (Precision Systems and Instrumentation) [11,12]. Manual bladder expression was performed on SCI mice a minimum of three times per day in the first two weeks of injury until reflex bladder emptying was re-established. Sham mice underwent the same procedure except for contusion. All surgical procedures and behavioral experiments performed were under protocols approved by the University of Maryland School of Medicine Institutional Animal Care and Use Committee.

Open field (OF): This test was used to assess spontaneous motor activity and anxiety-like behavior under dim light [11,12]. At the start of the test, each mouse was individually placed in a corner facing the wall of the open-field arena (40 cm × 40 cm). The mice were allowed to explore the chamber freely for a total duration of 5 min (min). On Any-maze, the open field apparatus was divided into a grid of 10 cm × 10 cm squares. The outer 10 cm zone was defined as the outer zone, whereas the center 20 × 20 cm² region was defined as the inner zone. Parameters such as total distance travelled, average and max speed, and time spent in inside and outside zones were recorded with the Any-maze tracking program (Stoelting Co., Wood Dale, IL, USA).

Basso mouse scale (BMS) for locomotion: For assessment of hindlimb locomotion recovery, mice were observed for at least 4 min by two trained observers after being placed in a flat, enclosed box (62 cm × 42 cm). Both observers were blinded to the genetic background of the mice. Mice were tested for BMS scores on day 1 and 3, followed by weekly assessments for up to 8 weeks post injury. The BMS was used to score hindlimb locomotion on a scale of 0–9 based on joint movement, weight support, plantar stepping, and coordination of the hind limbs, with 0 being complete paralysis and 9 being normal function [1].

CatWalk XT automated gait analysis: Gait analysis with CatWalk XT (Noldus, http://www.noldus.com/↗, accessed on 26 April 2024; RRID: SCR_004074) was performed in a darkened room with the protocol as mentioned in previous publications [10,13]. The CatWalk apparatus features a green illuminated walkway and records multiple parameters with a fast-action camera placed directly under the glass walkway. Parameters recorded include regularity index, print position, stride length, weight distribution, and gait postures. A single researcher blinded to the group information was tasked with animal handling and data analysis. At least three valid runs or complete walkway crossings were obtained for each mouse. Runs with mid-way stops, turns, or scaling of sidewalls were considered invalid and excluded from analysis.

Tail suspension (TS) test: This test was used to assess depressive-like behavior on week 8 as described previously [11,12]. Briefly, each mouse was suspended at a height of 28 cm with 3 M adhesive tape placed approximately ~1 cm from its tail tip. A timer of 5 min was set for the test period, while the duration of time spent immobile (passive hanging, motionlessness) was recorded with a stopwatch.

Forced swim (FS) test: For assessment of depressive-like behavior, we performed FS test as previously described [11,12]. Mice were placed in a transparent glass cylinder (25 cm high and 20 cm diameter) filled with water (22 cm in depth) at room temperature (22–23 °C). The duration of time spent immobile, which is characterized as lack of vigorous activity and movement only to maintain its head above water, was recorded with a stopwatch throughout the 6 min test period.

Von Frey filament test: For determination of mechanical allodynia development, the von Frey filament test was used to detect hind paws withdrawal from a mechanical stimulus according to a method published previously [9,14]. Prior to testing, each mouse was placed in individual plexiglass cubicles on a wire mesh platform and underwent an acclimation period of 1 h. The von Frey filaments (MyNeuroLab, St. Louis, MO, USA), which have an incremental stiffness ranging from 0.04 to 2.0 g, were applied serially to the plantar surface of each hind paw in ascending order. A total of five von Frey filament applications were made perpendicularly to the hind paw in each trial, each lasting 2–3 s and with a 5 s interval between each application. A positive response is defined as brisk withdrawal of the paw (at least 3 times out of 5 applications) upon application of the filament. Whereas a negative response was defined as no paw withdrawal, and the filament was changed to the next one in ascending stiffness. Furthermore, the application of a filament was repeated when paw response was ambiguous. The threshold was defined as the filament with the lowest bending force that elicited at least 3 positive responses out of 5 trials.

Hot plate test: Thermal allodynia of the hindpaws was tested with an Incremental Hot/Cold Plate Analgesia Meter (PE34, IITC Life Science, Woodland Hills, CA, USA) as previously published [9,14]. In brief, mice were placed within a transparent plexiglass box (22 cm length, 12 cm width, and 20 cm height) on the contact probe of a computerized thermal stimulator. The plate temperature was set to increase from 30 to 50 °C with the incremental ramp rate set at 10 °C per minute. Testing was stopped when the mouse licked either one of its hind paws, and the threshold temperature was recorded. The test was conducted twice with an interval of 3 h between the first and second trials.

At the study endpoint, anesthetized mice were transcardially perfused with ice-cold normal saline. Mouse spinal cord tissue (~0.5 cm) centered on the injury site was dissected out and fast-frozen in dry ice. Total RNA was extracted using an RNeasy mini kit (Cat# 74104, Qiagen, Hilden, Germany) with the concentration measured and sent to the Institute for Genome Sciences (IGS) at the University of Maryland School of Medicine for RNA sequencing [9,10].

Library construction and sequencing: Sequencing libraries were prepared by the IGS with the Illumina mRNA-Seq Sample Prep Kit (Illumina, San Diego, CA, USA. #RS-122-2101), followed by sequencing on an Illumina HiSeq4000 Sequencer that used the 150 bp paired-end protocol. Raw data from the Sequencer were processed using Illumina’s Real-Time Analysis (RRID: SCR_014332) and CASAVA pipeline software version 1.8.2 (RRID:SCR_001802) by IGS, which includes image analysis, base calling, sequence quality scoring, and index demultiplexing. Data were then processed through both FastQC version 0.12.0 (RRID:SCR_005539) and in-house pipelines for sequence assessment and quality control. We assessed basecall quality, and reads with a median Phred-like quality score of lower than Q20 were truncated.

Sequencing read alignment and quantification: We used STAR aligner version 2.4.2a [15], with the default parameters, to map the short sequencing reads onto the Genome Reference Consortium Mouse Build 39 (GRCm39) [16]. We used mouse gene annotations from Ensembl Release 100 [17] for the mapping. RSEM version 1.2.22 [18] was used for the quantification of gene expression in Transcripts Per Million (TPM) at both gene and isoform levels. We used DESeq2 version 1.42.0 [19] for the differential expression analysis between the SCI and sham samples. R “stats” package version 4.3.2 was applied for the k-means clustering. We used clusterProfiler 4.10.0 [20] for the Gene Ontology analysis. The analysis of Ntrk2 (trkB) splice variants were performed as in Pattwell et al. [21] using the RSEM-produced isoform-level TPM values. The full length and truncated isoforms of Ntrk2 gene have Ensembl transcript ID ENSMUST00000079828 and ENSMUST00000109838, respectively. Non-major isoforms (median TPM < 5) are not included in the figure.

All quantitative data are presented as individual data points in column graphs with mean ± SEM. Statistical analysis was performed on Sigma-Plot, version 12 (Systat Software; RRID: SCR_003210) for BMS scores and Graphpad Prism Version 9.5.0 for Windows (Graphpad Software, RRID: SCR_002798) for all other tests. For normality assessment of data distribution, the Shapiro–Wilk test was used.

For multiple comparisons between groups, one-way ANOVA analyses were performed followed by Tukey’s or Newman–Keuls multiple comparisons post hoc test for parametric (normality and equal variance passed) data. BMS scores were analyzed with two-way ANOVA for repeated measurements, followed by Holm–Sidak’s post hoc test for multiple comparisons. Statistical analysis in each assay was detailed in figure legends and a p-value ≤ 0.05 was considered statistically significant.

Multivariate data analysis was used to gain a comprehensive understanding of all the behavior tests in the three strains after mSCI. The behavior omics data included BMS, CatWalk, TS, FS, von Frey filament, and hot plate. The partial least squares discriminant analysis (PLS-DA) was conducted to classify mice behaviors based on the group effect of strain or injury using self-coded R language based on the ropls and mixOmics (version 6.20.0) packages [22].

To assess baseline locomotion differences between the three mouse strains, the open field (OF) test was used in these naïve mice. Baseline motor function was significantly lower in both albino mice strains, A/J and BALB/c, compared to C57BL/6 mice (Figure 1). In the test period of 5 min, C57BL/6 showed the highest total distance travelled (p < 0.001, n = 25 mice/group, Figure 1A), while no differences were observed between A/J and BALB/c mice. Moreover, the motor function of albino mice is also reflected in movement speed, with both mean (p < 0.0001, Figure 1B) and max (p < 0.0001, Figure 1C) movement speeds in the OF test being significantly lower than mice from the C57BL/6 strain. In addition to gross motor function, the OF test has also been widely used to assess anxiety-like behavior in animals [23]. Thus, we analyzed the time each mouse spent in the inside and outside zones of the OF arena (Figure 1D,E). Compared to C57BL/6 mice, the two albino strains spent more time in the outside zone (p < 0.0001, Figure 1D), suggesting that the two strains may have inherent tendencies towards anxiety-like behavior. Taken together, our findings suggest that strain differences occur in naïve mice.

To assess how differing genetic backgrounds affect functional recovery following SCI, we subjected young adult mice to either a moderate/severe or mild injury at the T10 spinal cord. Mice were tested for BMS scores on day 1, day 3, and weekly thereafter for up to 8 weeks post-injury. BMS scores analysis showed that significant genotypes were observed between the three strains. Statistical analysis with two-way repeated measurements ANOVA showed significant main effects in strain (F (1,39) = 9.507, p < 0.001, Figure 2A). Starting at 2 weeks (w) following a moderate/severe contusion-based SCI, C57BL/6 mice showed better recovery than the A/J strain (p = 0.004), with an average score of 2.9, indicating that most C57BL/6 mice exhibited extensive ankle movement and plantar paw placement. The differences in hindlimb motor function became even more prominent in the weeks thereafter (p < 0.001, Figure 2A), as the A/J strain continued to exhibit an average score of less than 2, indicating that the majority of mice only showed slight ankle movement. The other albino strain we tested, BALB/c, also demonstrated faster endogenous recovery than A/J mice, showing an average score of 2.8 and significant differences in statistical analysis (p = 0.01).

Due to the lack of endogenous recovery in the A/J mouse strain, we reasoned that a moderate/severe SCI may have been too high in injury severity. In addition to affecting the results of BMS scoring, the lack of plantar paw placement and stepping in A/J mice was an impediment for further neurological assessment of fine motor function and pain sensitivity. To this end, we subjected a separate set of mice to mild SCI and evaluated their recovery progress through weekly BMS. Overall, significant main effects of strain (F (1,33) = 8.829, p < 0.001, Figure 2B) were also observed in this second set. Moreover, strain differences were significant starting as early as 1 d post-injury, with C57BL/6 showing higher scores than either BALB/c (p < 0.001) or A/J (p < 0.001) strains. This high level of hindlimb motor function recovery remained in the weeks thereafter, as the scores for C57BL/6 mice plateaued around 8 in week 4, indicating an almost full recovery with slight deficits in fine motor coordination. In contrast, mice of the BALB/c and A/J strains were only able to reach a maximum score of 7.75 and 6.79, respectively, suggesting high prevalence of hindlimb rotation and trunk instability. In general, the overall trends of motor function recovery remained the same in both moderate/severe and mild-SCI groups, with C57BL/6 showing the fastest progress and the highest extent of recovery.

Next, we used Catwalk automated gait analysis to detect potential motor differences beyond that recognized with BMS scores in mild-SCI groups where animals can place the plantar surface of the hindpaws on a flat surface. For evaluation of motor coordination, we used a regularity index, a parameter that tracks the order of paw placement in a step cycle (Figure 2C). A single step cycle is defined as each of the four paws being placed on the walking surface in sequence. Through analysis of step sequencing, the program then attributes each set of steps into either a normal stepping pattern or abnormal gait, with the result being a percentage of normal stepping out of all step cycles analyzed. Ideally, the regularity index of healthy C57BL/6 mice would be 100% [10,24]; though in practice, the results vary. Whereas, both C57BL/6 and BALB/c mice had a regularity index of 80 to 90% under normal conditions, A/J mice had an index score close to 50%, indicating baseline deficits in motor coordination. Following SCI, the regularity index of neither C57BL/6 nor A/J showed significant differences between sham and injury groups, while only BALB/c mice showed a significant injury effect (p < 0.05, Figure 2C). Print position is defined as the distance between a pair of hindpaw and forepaw of the same side. Ideally, healthy C57BL/6 mice should be able to place their hindpaw next to the location of the forepaw that has just been lifted from the walkway. After a mild contusion, the print position of A/J mice showed no SCI-mediated differences, which is likely due to the already prominent deficits observed in the sham group (p < 0.001 vs. BALB/c/Sham, p < 0.0001 vs. C57BL/6/Sham, Figure 2D). For C57BL/6 and BALB/c mice, however, print position was dramatically increased in injured mice (p < 0.05). Stride length, which is defined as the total distance between steps of the same paw, was significantly lower in A/J mice compared to BALB/c and C57BL/6 strains, regardless of injury. However, no injury effects were observable in any of the three strains for this parameter (Figure 2E). Hindlimb base-of-support is a parameter that measures the average width of the track (distance between RH and LH) made by the animal, in which the farther apart the feet are placed during locomotion, the less likely the animal is to fall and the larger the base-of-support [24]. Strain differences in base-of-support were observed between C57BL/6 and BALB/c in sham mice (p < 0.001, Figure 2F), with C57BL/6 showing a smaller track width. After injury, C57BL/6 mice showed a dramatic increase in base-of-support (p < 0.01, Figure 2F), suggesting that mice showed compensation for instable gait. As an extension for gait stability assessment, we examined weight distribution, which calculates how many feet are simultaneously on the ground during a complete stride cycle. Injury effects were prominent in BALB/c mice, which saw a significant decrease in time spent in a diagonal stance (p < 0.01, Figure 2G), but no differences were seen in the time spent on girdle or four-paw stances. For A/J mice, significant SCI-induced differences were seen in time spent on a girdle stance (p < 0.0001, Figure 2H) and a four-paw stance (p < 0.01, Figure 2I), suggesting a redistribution of weight to compensate for motor deficits following injury. Taken together, the recovery process following SCI is significantly affected by the genetic background of mice strains, with C57BL/6 showing faster recovery and less deficits in overall motor function and coordination.

When the BMS scores reached the plateau, a battery of neuropsychiatric behaviors was applied to those animals. To further explore the effects of genetic background on the prevalence of SCI-induced depression, we ran learned helplessness tests on the two sets of mice subjected to SCI. In the moderate/severe injured groups, mice from all three strains showed increased immobility time in the tail suspension (TS) test following injury (p < 0.01 for A/J and C57BL/6, p < 0.001 for BALB/c mice, Figure 3A), but no strain differences were observed in either the sham or SCI groups. Similarly, the immobility time was also increased in the forced swim (FS) test after SCI (p < 0.001 for A/J and C57BL/6, p < 0.01 for BALB/c mice, Figure 3B), with no strain differences detected. However, in mice subjected to mild SCI, the test results show dramatic differences depending on genetic background. The TS test showed higher immobility time in BALB/c mice following SCI (p < 0.05, Figure 3C), but not in A/J or C57BL/6 mice. However, strain differences were detectable when comparing the SCI mice of A/J to BALB/c (p < 0.01, Figure 3C) and C57BL/6 (p < 0.01). In the FS test, SCI-induced depression-like behavior was observed in both A/J (p < 0.0001) and BALB/c (p < 0.01) mice (Figure 3D). Like the TS test, dramatic strain differences were detected between the SCI mice of the A/J strain compared to either BALB/c (p < 0.01) or C57BL/6 (p < 0.05) mice. In contrast to moderate/severe SCI, C57BL/6 mice following mild SCI did not show significant depression-like behaviors. In the FS test, we also observed dramatic strain differences between the sham groups of each strain, suggesting varying degrees of susceptibility to depression. Taken together, our data suggest that genetic background and injury severity are the two crucial factors involved in the development of depression-like behaviors following SCI.

To test whether genetic background is involved in the development of post-injury allodynia evoked by mechanical and thermal stimuli, we conducted nocifensive behavioral testing in mild SCI groups, in which animals were capable of placing the plantar surface of the hindpaws on a flat surface. In the von Frey test for mechanical allodynia, the withdrawal threshold of all three mice strains were significantly lower at 5 weeks post-injury compared to their respective sham groups (p < 0.01 A/J and C57BL/6, p < 0.05 in BALB/c, Figure 4A). However, when comparing strains, no differences were found for the mechanical threshold. The thermal threshold obtained via an incremental hot plate test showed no injury-induced hyperalgesia in A/J and BALB/c mice, but a dramatic decrease in C57BL/6 mice at 6 weeks post-injury (p < 0.01, Figure 4B) was observed. Additionally, there is a significant strain difference in both the sham and SCI groups, with C57BL/6 mice having a significantly lower threshold than the two albino strains.

Next, pain sensitivity was further explored with relevant parameters obtained from Catwalk gait analysis, which reflects spontaneous pain-type behaviors. Maximum contact area was defined as the total area of the hindpaw at the point of maximum placement, which showed no injury effects in either A/J or BALB/c mice at 8 weeks post-injury (Figure 4C). However, strain differences between sham groups should be noted, with C57BL/6 mice showing the highest contact area and the A/J mice showing the lowest. Following SCI, C57BL/6 mice showed a dramatic decrease in max contact area, which suggests the development of spontaneous pain. In tandem with this parameter was the mean print area, which also saw a dramatic decrease in C57BL/6 mice after SCI (p < 0.01, Figure 4D), but no injury effects in either A/J or BALB/c mice were observed. Similarly, strain differences were also observed between both the sham and SCI groups, suggesting fundamental variations in weight-bearing and gait patterns based on genetic background.

To collectively delineate the total behavioral profiles of the animals from different strains and injured groups, all the behavioral data from the mild-SCI set of the mice were applied for a machine learning approach, partial least squares discriminant analysis (PLS-DA) [25]. Remarkably, we observed that mice from each strain are tightly clustered in the result (Figure 5). This indicates that the genotypes of the three strains are highly correlated with their behaviors of locomotion, depression, and neuropathic pain.

To examine the genomics between mouse strains following SCI, after completion of the behavioral tests at 8 weeks post-injury, bulk RNAseq from sham or injured spinal cord tissues was employed after moderate/severe SCI. Gene expression changes were analyzed with DEseq2 following normalization [19]. Principal component analysis (PCA) of all normalized gene counts revealed a clear separation of sham and SCI samples into individual groups across the first two principal components (PC, Figure S1). PC1 was the major principal component with a contribution of 67.59% of gene variations separating injury groups from sham groups, which may represent injury effects. PC2 contained 17.4% of gene variations and separated different strains of mice from sham or injured groups, which may represent strain effects. Differentially expressed genes (DEGs) with a Benjamini–Hochberg adjusted p-value of less than 0.05 in the Wald test were used for downstream analysis. DEGs between the SCI mice and sham group of each strain were presented in the heatmap for visualization after k-means clustering (Figure 6A), with each column representing the fold difference in log2 scale. Based on the heatmap results, we were able to identify 11 clusters of DEGs for further pathway enrichment analysis with the gene ontology (GO) biological process database. GO analysis of cluster 1, which consisted of DEGs downregulated in all the three strains by SCI, revealed genes involved with the regulation of membrane potential as the top enriched pathway (Figure 6B). Other pathways indicated for cluster 1 include learning or memory, dendrite development, and regulation of synapse activity or organization. These pathways suggest downregulation of genes contributing to the normal homeostasis of membrane potential and synaptic function after SCI. In contrast, the upregulated genes that comprised cluster 6 were mostly involved with the regulation of immune effector and innate immune response (Figure 6C), in which DEGs overlapped with more than 10% of the genes in this functional category. Of the other top enriched GO terms, leukocyte proliferation, regulation of response to biotic simulations, and leukocyte cell adhesion were in the top five. The GO terms for this cluster indicate the presence of persistent neuroinflammation in the injured spinal cord for all three strains. Consistently, GO analysis of cluster 5 showed highly upregulated immune genes in response to injury in all three mouse strains (Figure S2A). However, cluster 2 showed upregulated genes in the other strains, but not in A/J mice (Figure S2B). Of note, some genes were downregulated most prominently in A/J or BALB/c strains, which can be observed in cluster 11 (Figure 6D) and cluster 9 (Figure S2C). GO analysis of these clusters showed DEGs were highly enriched with genes regulating the cholesterol metabolic process and muscle organ development. Other functions implicated via GO analysis include 2nd alcohol metabolic process, sterol metabolic process, and potassium (K+) ion transmembrane transport and K+ ion transport. Cluster 7 represented downregulated genes in all three mouse strains after SCI (Figure S2D). The genes in this cluster are about 50% more downregulated in C57BL/6 mice.

Delving deeper into the specific genes of each cluster, we found significant upregulation of inflammatory genes (Figure 6E). DEGs of interest include Atf3 (activating transcription factor 3), a member of the mammalian activation transcription factor/cAMP responsive element-binding (CREB) family and a gene involved with cellular stress response and neuronal damage [26]. Moreover, Ccl2 (C-C Motif Chemokine Ligand 2) and its corresponding receptor Ccr2 (C-C Motif Chemokine Receptor 2) were upregulated, both of which are involved with the development of nociceptive behavior following SCI [27]. An interesting phenomenon is the downregulation of Il6 (Interleukin 6) in C57BL/6 mice, which studies found are crucial to the development of neuropathic pain following SCI [28,29], whereas it was upregulated in both A/J and BALB/c mice. In addition, the proinflammatory cytokines Il1b (Interleukin 1 beta) and Tnf (Tumor necrosis factor) were upregulated in all three strains. Furthermore, the genes Prkcd (Protein Kinase C Delta) and Ptgs2 (Prostaglandin-Endoperoxide Synthase 2), were upregulated more prominently in A/J mice compared to other strains. In order to explore the effects that differing genetic backgrounds may have on classic receptors involved with neuropathic pain, we examined the expression changes in TRP (transient receptor potential) receptors (Figure 6F). All three genes examined (Trpa1, Trpm1, and Trpv1) were upregulated more significantly in A/J mice following injury compared to BALB/c and C57BL/6 mice. Notably, Trpa1 (Transient Receptor Potential Cation Channel Subfamily A Member 1), a gene that is involved in familial episodic pain syndrome [30] and other forms of nociceptive behavior [31,32,33], was highly upregulated in A/J mice, but only minor changes were induced by SCI in the other strains. In contrast, all DEGs from the NMDA/GABA receptor family were found in clusters 1 and 4, showing significant downregulations after SCI (Figure 6G). All DEGs from the potassium channels were found in cluster 11, showing significant downregulations following SCI in all three mouse strains (Figure 6H).

Moreover, we found significant upregulation of cell cycle genes after SCI in all three mouse strains (Figure 7A), which includes Ccnd1 (cyclin D1), Ccnd2 (cyclin D2), and Cdk4. Strain differences were most prominent in Ccnd2, known to be involved in secondary injury damage following SCI [6]. In agreement with our previous studies [9,10], SCI significantly upregulated a truncated isoform of Ntrk2 (trkB.T1) but downregulated the full length of trkB (Figure 7B). Strain differences were most prominent in Ntrk2, known to be involved in SCI-PAIN [9,10]. Together, our data demonstrate differential gene expression in three mouse strains after chronic SCI which is associated with genetic background.