Abstract
In many preclinical spinal cord injury (SCI) studies, assessment of locomotion recovery is key to understanding the effectiveness of the experimental intervention. In such rat SCI studies, the most basic locomotor recovery scoring system is a subjective observation of animals freely roaming in an open field, the Basso Beattie Bresnahan (BBB) score. In comparison, CatWalk is an automated gait analysis system, providing further parameter specifications. Although together the CatWalk parameters encompass gait, studies consistently report single parameters, which differ in significance from other behavioral assessments. Therefore, we believe no single parameter produced by the CatWalk can represent the fully-coordinated motion of gait. Typically, other locomotor assessments, such as the BBB score, combine several locomotor characteristics into a representative score. For this reason, we ranked the most distinctive CatWalk parameters between uninjured and SC injured rats. Subsequently, we combined nine of the topmost parameters into an SCI gait index score based on linear discriminant analysis (LDA). The resulting combination was applied to assess gait recovery in SCI experiments comprising of three thoracic contusions, a thoracic dorsal hemisection, and a cervical dorsal column lesion model. For thoracic lesions, our unbiased machine learning model revealed gait differences in lesion type and severity. In some instances, our LDA was found to be more sensitive in differentiating recovery than the BBB score alone. We believe the newly developed gait parameter combination presented here should be used in CatWalk gait recovery work with preclinical thoracic rat SCI models.
- CatWalk
- gait parameter
- linear discriminant analysis
- locomotion recovery
- preclinical development
- spinal cord injury
Significance Statement
As a quantitative locomotion analysis system, CatWalk provides an objective assessment of gait for rodents by computing numerous parameters. This gives an alternative to the current popular subjective locomotor assessment for rat spinal cord injury (SCI) models, i.e., Basso Beattie Bresnahan (BBB) score. As SCI affects multiple gait parameters, analyzing gait is challenging. Here, we developed a CatWalk gait parameter combination for sensitive and efficient gait recovery assessment in rat thoracic SCI models.
Introduction
Following a spinal cord injury (SCI), disruption to descending fine and gross motor control, as well as to ascending sensory input, leads to functional deficits, all the way up to paralysis. Standardized and consistent behavioral assessments are required to accurately monitor any spontaneous or induced recovery that follows.
Preclinical rat SCI studies used in developing treatments for SCI (Cheriyan et al., 2014; Silva et al., 2014; Kjell and Olson, 2016; Ahuja et al., 2017) commonly rely on a motor assessment based on a standardized scale involving the observation of various characteristics to generate a non-parametric score (Bhimani et al., 2017). For rat thoracic contusion SCI models, in particular, Basso et al. (1995) developed a standard non-parametric locomotion assessment named Basso Beattie Bresnahan (BBB) locomotion scale. This subjective locomotion assessment scores hindlimb (HL) movement, joint movement, weight supported plantar stepping, forelimb (FL)–HL coordination, toe clearance, trunk stability, and tail placement in an open field test setting. In the BBB locomotion assessment, the examiners’ experience attributes a significant aspect of the score’s reliability (Basso et al., 1996).
CatWalk (Noldus Information Technology), on the other hand, is a locomotion analysis system specifically assessing gait performance based on rodent footsteps and body silhouette. Introduced by Hamers, CatWalk was initially designed to assess FL–HL coordination (Hamers et al., 2001, 2006). Using this system, paw positions during walking can be recorded using a high-speed camera positioned beneath the glass walkway (Fig. 1). Based on the recorded video and paw positions, the CatWalk system computes a multitude of static and dynamic gait parameters, such as the base of support (BOS), stride length, swing speed, body speed, and regularity index (RI), to name a few. For example, the RI (Koopmans et al., 2005) is a gait parameter used to quantify FL–HL coordination based on the percentage of the defined normal step sequence patterns (NSSPs). The CatWalk apparatus provides further information than the traditional analysis of inked footprints collected on paper rolls, which can only deliver static gait parameters.
An example of CatWalk image captured from below an injured rat. The paws in contact with the glass walkway are indicated with the colored boxes (magenta: right hindpaw; yellow: left forepaw).
Since SCI affects several motor and sensory pathways, multiple gait parameters need to be analyzed simultaneously (Lakes and Allen, 2016). However, in the majority of published work, CatWalk parameters are evaluated individually (Garcia-Ovejero et al., 2014; Fagoe et al., 2016; Bhimani et al., 2017; Kappos et al., 2017; Slusarczyk et al., 2017; Bieler et al., 2018; Aceves et al., 2020; Heinzel et al., 2020). The large number of gait parameters as well as the variated gait effect of an injury often bring complications in the analysis of gait. Therefore, a combination of CatWalk gait parameters to faithfully characterize gait recovery in rat SCI models will be a valuable asset for future research studies.
Simple CatWalk parameter combinations have been used in a variety of disorders to assess gait, such as baseline parameter ratio (Datto et al., 2015, 2016; Hayakawa et al., 2015; Mondello et al., 2015), left-right-parameter ratio, left-right-parameter averaging, and subtraction (Liu et al., 2013; Parvathy and Masocha, 2013; Chen et al., 2014, 2017; Ishikawa et al., 2014; Muramatsu et al., 2014). In addition to these basic operations, researchers have used principal component analysis (PCA) on multiple CatWalk parameters to determine the most relevant parameters for locomotor phenotyping of mice (Zimprich et al., 2018). Also, linear regression of the Basso Mouse Scale (BMS) against each CatWalk parameter have been combined into a weighted score in thoracic mouse contusion SCI (Crowley et al., 2018), and classifiers developed based on high-dimensional CatWalk-derived data to classify gait pattern in Parkinson’s disease-relevant animal models and control groups (Fröhlich et al., 2018).
In contrast to this previous work, we propose a CatWalk parameter linear combination, which has been developed based on a linear discriminant analysis (LDA) of several gait parameters ranked according to the difference measured between uninjured and spinal cord injured rats. LDA was first introduced by Fisher (Fisher, 1936) as a method to determine a linear function that best discriminates two groups. Based on the most descriptive gait parameters identified, using LDA we were able to quantify gait recovery of different rat SCI models using a weighted linear combination of parameters. The designed LDA was then applied and validated to assess gait recovery in various SCI models, ranging from contusion, dorsal hemisection to dorsal column lesion. The gait recovery profile was reliably uncovered from most of these studies even when differences were not detectable by single CatWalk parameters or by the BBB score alone.
Material and Methods
Subjects and study protocols
In rats, the gait recovery assessment method was applied to three Th8/9 contusion SCI studies, one Th8/9 dorsal hemisection SCI, and one C4 bilateral dorsal column transection SCI, which are summarized in Table 1. The studies were conducted using two different versions of the CatWalk system, which have different camera frame rates (50 fps for CatWalk 7.1 and 100 fps for CatWalk XT).
Overview of the studies
Study 1
This study was performed as published previously (Sandner et al., 2018) to examine the effects of Epothilone D (Epo D) on the functional recovery and regenerative potential following contusion SCI. For the Th8/9 contusion SCI model, adult female Fischer 344 rats (Charles River Deutschland GmbH, Envigo, Janvier Labs) weighing 160–180 g were used (10 weeks of age). Experiments were conducted in accordance with the European Union Directive (2010/63/EU) and institutional guidelines. Rats were fed ad libitum. The rats were anesthetized using a cocktail of ketamine (62.5 mg/kg; Medistar), xylazine (3.175 mg/kg; Bayer), and acepromazine (0.625 mg/kg; Sanofi-Ceva) in 0.9% sterile saline solution. The rats received a contusion SCI at midthoracic Th9 level (representing spinal level Th11) using the Infinite Horizon (IH) Impactor SCI device (Precision Systems and Instrumentation) with an impact force of 150 kilodynes (kDyn) without any additional dwell time. Two days postoperatively, twice a day, buprenorphine (0.03 mg/kg; Reckitt Benckiser) was given subcutaneously. Ampicillin (167 mg/kg; Ratiopharm, twice a day) was given subcutaneously as long as manual bladder evacuation was necessary. Rats received intraperitoneal injections of Epo D [Abcam, catalog #ab143616; 1.5 mg/kg body weight (bw)] dissolved in dimethylsulfoxide (DMSO; 3 mg/ml) and diluted 1:1 with prewarmed saline before injections of Epo D or vehicle (1:1 mixture of DMSO and saline, control group) on days 1 and 15 postinjury by a blinded unbiased experimenter. An animal in each cohort was excluded because of inadequate force impact curves on spinal cord contusion. Based on the variability in spinal cord displacement, rats were divided into two groups. A displacement value of 1000 μm allowed for equal distribution of animals into two cohorts (moderate: <1000 μm and moderate-severe: >1000 μm). These rats were then assigned to receive either vehicle or Epo D treatment. Rats were killed eight weeks postinjury.
Behavioral testing
At 60 d postinjury (dpi), CatWalk (CatWalk XT), and BBB tests were performed. The number of rats and CatWalk runs collected in this study are listed in Table 2. The BBB scores ranged from 10 to 19 at 60 dpi.
Number of animals and CatWalk runs of Study 1
Study 2
This study was performed as published previously (Garcia-Ovejero et al., 2014) to examine the effects of progesterone on locomotor recovery and secondary damage following contusion SCI. Rats were handled in accordance with the guidelines published in the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals, the principles laid out in the Guidelines for the Use of Animals in Neuroscience Research published by the Society for Neuroscience, and European Union guidelines (Council Directive 86/609/EEC). Experimental procedures were approved by the Ethical Committee for Animal Welfare at the National Paraplegics Hospital (CEBA). For the Th8 contusion SCI model, young adult male Wistar rats (300–335 g, 12 weeks of age) were submitted to a moderate-severe contusion SCI. Briefly, rats were anesthetized with an intraperitoneal injection of sodium pentobarbital (45 mg/kg, Normon Veterinary Division) and Xilagesic (2% xylacine, 10 mg/kg, Calier). Once the absence of reflexes was confirmed, the rats were injected with a low dose of atropine (50 μg/kg bw; Brown Medical). After removing the Th8 vertebra, spinal cord contusion was performed with the IH impactor (Precision Systems and Instrumentation), applying a force of 200 kDyn without any additional dwell time. Postoperative care included a subcutaneous injection of Buprex (buprenorphine, 0.05 mg/kg; Schering Plough) and a prophylactic sub-cutaneous antibiotic injection 1 h after the lesion and on the following day (Baytril, Enrofloxacine, 1 mg/kg; Bayer). Injured rats received daily subcutaneous injections of either natural progesterone (16 mg/kg/d, Sigma-Aldrich, SCI+Prog. group) or vehicle (Castor oil, Sigma-Aldrich, SCI group) for 60 d until killing. The first injection was given to awake animals 1 h after injury.
Behavioral testing
Gait data were collected using CatWalk 7.1. The BBB scores (Basso et al., 1995) were additionally scored for coordination with the CatWalk RI (Koopmans et al., 2005), named RI-controlled BBB score. The rats with BBB scores of 9 or higher were further examined weekly by CatWalk 7.1 analyses. The number of rats and runs included in the data collection are listed in Table 3.
Number of animals and CatWalk runs of Study 2
Study 3
Th8 contusion SCI model
In this yet unpublished study, motor recovery after traumatic mid-thoracic SCI in three-month old female Fischer 344 rats (Charles River Deutschland GmbH; 160.1 ± 7.1 g) was compared with rats of 20–24 months of age (273.7 ± 23.7 g). Rats underwent a laminectomy at Th8 vertebral level followed by a 200 kDyn contusion generated by an IH Impactor (Precision Systems and Instrumentation) without additional dwell time. For surgical purposes, rats were under general inhalative anesthesia using a SomnoSuit unit (Kent Scientific) with 2% isoflurane/O2 gas mix. Analgesia was provided by subcutaneous injection of 0.05 mg/kg bw buprenorphine (Bupaq, Richter Pharma). Heart rate and oxygen saturation were monitored (SomnoSuit) and body temperature was maintained by a sensor-driven heating pad (Harvard Apparatus). Postoperative care included subcutaneous injection of 0.03 mg/kg bw buprenorphine twice daily for 2 d, 1–2 mg/kg bw meloxicam (Metacam; Richter Pharma) and 10 mg/kg bw Baytril (Bayer Pharma) daily for 7 d. The bladder was manually voided twice daily and once daily after reflexive bladder voiding was established (∼14 d).
Rats were housed in groups of 5 in a 12/12 h light/dark cycle and had food and water ad libitum. Experimental procedures were authorized by the ethical committee of the “Land Salzburg” (20910-TVG-79/17-2014) according to the European Directive 2010/63/EU on the protection of animals used for scientific purposes.
Behavioral testing
BBB tests were performed on days 1, 4, and 7 postinjury and weekly thereafter. The rats with BBB scores of 11 or higher were further examined weekly by CatWalk XT analyses. For the CatWalk XT analyses, rats were trained on the CatWalk device for at least two weeks (two times per week) before surgery and a baseline was created. Rats did not receive food rewards during the testing. The CatWalk data acquisition ended after 43 d for the young and 29 d for the aged group. The number of animals and runs included in the data collection are listed in Table 4.
Number of animals and CatWalk runs of Study 3
Study 4
This unpublished work was designed to study the effect of various undisclosed compounds on the recovery of locomotion following SCI. For the Th9 dorsal hemisection Scouten wire-knife transection SCI model, adult female Wistar rats (M&B Breeding) weighing 220–250 g at the time of operation were anesthetized with isoflurane (Forene, Abbott; 2–3% in O2 and N2O at a ratio of 1:2). Following laminectomy at Th8/9, the dura mater was opened at Th8/9 with a longitudinal cut and a dorsal hemisection injury was performed with a Scouten wire knife (Bilaney). After suture of the dura mater, a polyurethane rat intradural catheter (32 G, ReCathCo) was inserted at Th11 and epidurally guided to the lesion site. The catheter was filled with vehicle solution before it was inserted into the subarachnoid space in close proximity to the dura suture. Subsequently, it was connected to a prefilled osmotic minipump (Alzet pump model 2002), which was placed subcutaneously. Finally, the lesion area was covered with a piece of Nescofilm (Carl Roth), and overlying muscle and skin were sutured in layers. Immediately after surgery, animals received subcutaneous injections of 5 ml of physiological saline and 5 mg/kg carprofen (Rimadyl, Pfizer). Individual caging was provided until the animal had fully recovered from anesthesia. Postoperative care included manual bladder expression until normal bladder function returned, prophylactic antibiotic treatment (Baytril; Bayer Healthcare) for one week, and pain relief (Rimadyl; Pfizer) for 2 d postlesion. After a two-week infusion period, the osmotic minipump was removed under short anesthesia. During all surgical procedures, animals were placed on a heating pad to maintain body temperature. Institutional guidelines for animal safety and comfort were adhered to, and all surgical interventions and presurgical and postsurgical animal care were provided in compliance with the German Animal Protection law (State Office, Environmental and Consumer Protection of North Rhine-Westphalia, LANUV NRW, AZ 8.87-50.10.34.08.061). Experimental animals were housed in groups under standard conditions. Water and food were available ad libitum.
Behavioral testing
Four weeks before surgery, all animals were familiarized and pretrained in the behavioral tests. The overall HL function was assessed in an open field using the BBB score. Freely exploring rats were observed by two blinded examiners at 30 and 60 d postlesion and their HL movements were rated according to the BBB open field score. Differences in walking patterns were investigated using the CatWalk XT. The number of animals and runs involved in the CatWalk test are shown in Table 5. Food rewards were used throughout the entire study, and animals were trained to cross the horizontal glass runway without interruption.
Number of animals and CatWalk runs of Study 4
Exclusion criteria
Some animals showed signs of automutilation. In the case of severe automutilation, the respective animal was killed prematurely. If only minimal signs of automutilation were detected, the respective animal was not included in the behavioral tests at individual test time points.
Study 5
This study was performed as published previously (Bieler et al., 2018) and designed to study the motor deficits following a cervical cortical spinal tract SCI. For the C4 bilateral dorsal column tungsten wire-knife transection SCI model, this study was conducted on female Fischer 344 rats (169 ± 8 g bw). Experiments were performed in accordance with the Directive 2010/63/EU of the European Parliament and of the Council and were approved by the local animal health commission (Land Salzburg 20 901-TVG-65/8-2013). Rats were purchased from Charles River Laboratories and were housed in groups of five in standard conditions, i.e., a 12/12 h light/dark cycle, and food and water were provided ad libitum.
The lesion rat group underwent a bilateral transection of the dorsal column at the fourth cervical segment. Before surgery, rats were deeply anesthetized using an intramuscular injection of 46.5 mg/kg bw ketamine (Narketan 10%, Vétoquinol), 2.3 mg/kg bw xylazine (Rompun, Bayer Austria GmbH) and 0.46 mg/kg bw acepromazine (Vanastress, Vana GmbH). Rats were placed on a homeostatic heating pad with body temperature monitoring via a rectal sensor to prevent hypothermia. The dorsal spine of the rat was exposed and a laminectomy at C4 was performed to expose the spinal cord. Using a blunt tungsten wire-knife device (M122, David Kopf Instruments), the dorsal column was precisely transected bilaterally (2.5-mm width, 1.1-mm depth) as previously described (Weidner et al., 1999, 2001; Sandner et al., 2013; Bieler et al., 2018). To prevent infections after surgery, rats were treated with 10 mg/kg bw Enrofloxacin (Baytril, 2.5% injection solution, Bayer Austria GmbH) daily for five consecutive days starting peri-operatively. Analgesia was provided by daily subcutaneous administration of 1.25 mg/kg bw Meloxicam (Metacam, Boehringer Ingelheim Vetmedica GmbH) for five consecutive days after surgery. In the first 2 d after surgery, 0.02 mg/kg bw buprenorphine (Bupaq, Richter Pharma AG) was injected subcutaneously twice a day. Upon signs of dehydration, 1–2 ml of 0.9% NaCl solution was injected subcutaneously. In this SCI model, bladder function is sufficient and did not require manual voiding.
Behavioral testing
Differences in gait parameters were investigated using the CatWalk XT. Two weeks before injury, rats were familiarized with the device and testing conditions. Rats crossed the CatWalk voluntarily without the use of food rewards. Rats were tested 30 dpi. The number of animals and runs involved in this study are shown in Table 6.
Number of animals and CatWalk runs of Study 5
Establishment of the CatWalk parameter linear combination
The parameter linear combination was built on gait parameters extracted from the recorded videos by the CatWalk software. Only the gait parameters available in both the CatWalk 7.1. and CatWalk XT were considered for this study. Gait parameters from manually validated labeled runs were exported as run statistics. From 317 CatWalk XT parameters and 228 CatWalk 7.1. collected parameters, we included gait parameters related to mean values of paw statistics (60 parameters), step sequence (RI and sequence AB), BOS (two parameters), print position (two parameters), and number of paws supporting the walk ([Support_One + 2 × Support_Diagonal + 2 × Support_Girdle + 2 × Support_Lateral + 3 × Support_Three + 4 × Support_Four]/100). However, parameters directly related to green-pixel-intensity (ranged 0–255 in pixel-intensity units) were not used as they are often too sensitive to the experimental setting and walkway moisture. This intensity is closely related to paw pressure or weight support (Hamers et al., 2006). However, this exception does not include parameters that measure time (s) duration when max contact was realized. Gait parameters from the right and left paws were averaged except for body speed, which was averaged from all four paws, as performed in a recent PCA of CatWalk data from mice (Zimprich et al., 2018). This resulted in 22 parameters for consideration of ranked gait parameters by t test.
In order to develop a comprehensive measure which would efficiently examine the gait recovery differences between experimental and vehicle-treated rat groups, we had to first designate which gait parameters are the most descriptive of a spinal cord injured animal. To do this, we grouped uninjured rats from Studies 1 and 2 as controls. Our SCI group was pooled from vehicle moderate-severe SCI rats from Study 1 along with the vehicle SCI rats (60 d post-SCI) from Study 2. Next, we ranked gait parameters with a t test without assuming equal variances between these two groups for all individual gait parameters. To perform this t test, parameters from all runs of each rat were averaged. All gait parameters having p < 0.01 from the t test were used to build the linear combination of CatWalk parameters. Of the parameters that reached our defined threshold, the effect size examined by Cohen’s d was found to be >1.3, which is considered to be a very large effect size (Sullivan and Feinn, 2012).
The linear combination of CatWalk parameters pLDA can be expressed as Equation 1, where pi is the i-th CatWalk gait parameter, w1,i is the corresponding parameter-weight, and Np is the number of gait parameters.
(1)
The parameter-weights w1,i of the linear combination were computed based on a machine learning method, namely LDA, which aims to maximize the between-class scatter and minimize the within-class scatter. Here, the LDA used gait parameter values from all runs. The method starts with calculating the within-class scatter matrix Sw and between-class scatter matrix Sb as shown in Equations 2, 3, where
is the i-th run data of group j,
is the mean of group j, c = 2 is the number of group, Nj is the number of run data in class j, and
μ represents the mean of all classes.
(2)
(3)
To maximize the between-class scatter and minimize the within-class scatter, an eigen value problem as shown in Equation 4 was defined, where w is the eigen vector and
λ is the eigen value.
(4)
Here, we took the first eigen vector w1, which has the highest eigen value, and project our data to this eigen vector as given in Equation 5.
(5)
Analysis
Here, the pLDA parameter combinations were calculated for all run data from Studies 1 and 2 (total 107 run data, 57 run data from Study 1 and 50 run data from Study 2, total 25 animals, 4.28 ± 1.28 runs per animal). Thereafter, the differences between groups in each study were analyzed based on the resulting parameter combination pLDA from the average run data of each animal. For relative locomotor comparison, a statistical analysis of the thoracic SCI rats’ BBB score was also performed.
The statistical analyses between two groups were performed by t tests without assuming equal variances. The differences between three or more groups were investigated by ANOVA and multiple comparison tests with Bonferroni correction provided in MATLAB R2015a (8.5.0.). Groups were considered significantly different when the p values were below 0.05 (p < 0.05).
Results
The purpose of this study was to develop a combination parameter for CatWalk to be more predictive of gait recovery after SCI than any single CatWalk parameter considered alone or even more than the widespread BBB score (see Materials and Methods and below). Subsequently this parameter combination (built from Studies 1 and 2) was applied and validated in the assessment of gait recovery in the various experimental SCI paradigms shown here.
CatWalk parameter linear combination
The gait parameters differing between control and SCI rats from Studies 1 and 2 with a p < 0.01 from the t test without assuming equal variances (effect size, Cohen’s d > 1.3), as well as the values of the first eigen vector w1, are listed in Table 7. Although these SCI models directly affect HL function, the list of parameters contains several forepaw-related gait parameters. This confirms that FLs play an important and reinforced role after an SCI for the support and balance of bw during walking (Ghosh et al., 2010; Wilcox et al., 2017).
List of gait parameters pi (Np = 9) and their corresponding parameter-weight w1,i for the linear combination of parameters
The gait parameter combination,
, obtained by a machine-learning LDA-based computational method described above (see Establishment of the CatWalk parameter linear combination), aims to maximize the between-class scatter and minimize the within-class scatter of data collected from Studies 1 and 2. A machine learning method that generates a simple linear combination of parameters was chosen for ease of use. In applications of data classification, LDA performs optimally in normally distributed data (Liong and Foo, 2013). Therefore, normality tests were performed.
Normality tests using the Shapiro–Wilk W test (using R version 3.5.1) were done on the run data for these parameters (Studies 1 and 2 combined). From the nine parameters, the normality tests were run separately for uninjured and injured groups for the seven parameters. The normality test was not relevant for two parameters, RI (%) and sequence AB (%). All RI (%) values from the uninjured group were at their maximum value of 100% (equal to a normal distribution with a variance of zero). From this uninjured group, 80% of the sequence AB (%) values were at their maximum value of 100%. From the injured group, around 30% of the data shown at their maximum value of 100%. Therefore, the distributions resemble one-sided normal distributions. From the 14 normality tests, 10 parameters show no difference from a normal distribution with WShapiro–Wilk > 0.95, and four parameters show similarity with a normal distribution with WShapiro–Wilk > 0.84. These four parameters only show similarity with a normal distribution because of the combination of two studies involving characteristic differences of rat strains used.
Here, the LDA was used to generate a linear combination of gait parameters based on the resulting first eigen vector, pLDA. The LDA was not used to classify data nor to statistically analyze data. Therefore, we considered the similarity with the normal distribution (WShapiro–Wilk > 0.84) adequate for this objective. RI (%) and sequence AB (%) have been shown as important parameters in several studies (Hamers et al., 2006; Šedý et al., 2008; Datto et al., 2015; Crowley et al., 2018); therefore, we considered that their distributions would not hinder the efficacy of generating a linear combination of gait parameters. The gait parameters are interdependent on each other and some weights act to balance the influence of other gait parameters. The resulting parameter combination LDA was used to assess the progression of gait recovery in all following studies presented here. An example of the pLDA calculation can be found at https://github.com/mad-lab-fau/GaitParamCombSCI.
Specific definitions of the given parameters included within the pLDA are listed below:
Swing time (s) is calculated by the duration of non-contact with the walkway of a specific paw.
Stride length (cm) is calculated by the distance between the center points of two consecutive positions of the same paw.
Duty cycle (%) is the ratio of stand time to step cycle (duty cycle = stand time/step cycle). Stand time (s) is calculated by the duration of contact with the walkway of a specific paw. Step cycle (s) is calculated by the duration of two consecutive initial contacts of a specific paw (step cycle = stand time + swing time).
BOS (cm) is calculated by averaging the width on the y-axis between either forepaws or hindpaws (BOS forepaws = yRF – yLF; BOS hindpaws = yRH – yLH).
RI (%) is the ratio of the number of NSSP times four and the total number of paw placement (PP; Koopmans et al., 2005); RI = NSSP × 4/PP. There are six patterns that are considered as NSSP (Cheng et al., 1997): cruciate (CA: RF–LF–RH–LH; CB: LF–RF–LH–RH), alternate (AA: RF–RH–LF–LH; AB: LF–RH–RF–LH), and rotary (RA: RF–LF–LH–RH; RB: LF–RF–RH–LH). RF: right forepaw; LF: left forepaw; RH: right hindpaw; LH: left hindpaw.
Body speed (cm/s) is calculated from each paw by dividing the distance by the time of two consecutive initial contacts.
Sequence AB (%) is the percentage of a specific LF–RH–RF–LH paw sequence.
“Max contact at (%)” is the ratio of “max contact at (s)” to the stand time (s), where “max contact at (s)” is the time when maximum contact with the walkway is measured. This is calculated from each paw and related to the point in which the braking phase turns into the propulsion phase.
It should be noted that the ranked parameters used in our combination consist of spatial dimensions (stride length, BOS and speed), temporal dimensions (swing, duty cycle, speed and max contact at), interlimb dimensions (AB sequence) and variation of interlimb coordination (RI) that encompass gait. All of these parameters were previously reported to provide the most variance in CatWalk based on PCA in mouse phenotyping experiments (Zimprich et al., 2018).
The nine highest-ranked parameters between uninjured and injured groups cannot be represented in a single graph, however we show here coupled parameters that visualize the separation of injured and uninjured groups for Study 1 only. Figure 2 compares the two parameters (swing time and stride length) that were found to be significant in parameter ranking. Another two parameters (duty cycle and hindpaw BOS) are also shown in Figure 2, again depicting the separation of the uninjured and injured groups. Not all of the observed parameters depict differences. For example, two parameters (forepaw BOS and number of paw support) did not depict this separation, were not found to be significantly different and thus were not used (Fig. 2).
Scatter plot of several gait parameters obtained from the individual run data from Study 1. The injured (red) and uninjured (dark gray) animals are denoted by color here. The different shapes in each group represent the different individuals. a, Forepaw swing time and forepaw stride length. b, Forepaw duty cycle and hindpaw BOS. c, Forepaw BOS and number of paw support = [Support_One + 2 × Support_Diagonal + 2 × Support_Girdle + 2 × Support_Lateral + 3 × Support_Three + 4 × Support_Four]/100.
Testing of the CatWalk parameter linear combination in distinct contusion SCI severities (Study 1)
We tested the ability of the pLDA to distinguish differing lesion severities, as well as experimental treatment groups using the CatWalk data of Study 1 (Th8/9 moderate and moderate-severe contusion SCI). We observed that the uninjured rats had significantly higher values of pLDA (0.12) compared with the injured rats (60 dpi, pLDA 0.03–0.06), as shown in Figure 3A. Furthermore, pLDA distinguished gait performance of the moderate (0.06) and moderate-severe (0.03) injured rats significantly. Treatment used in this study did not improve gait according to pLDA, as shown in Figure 3C,E. These results were similar to those obtained with the BBB scores (Fig. 3B,D,F). The values of pLDA and nine CatWalk parameters of vehicle moderate and vehicle moderate-severe groups are given in Table 8, as well as their p values calculated from the t test without assuming equal variances. It is important to note that most of the single CatWalk parameters did not display significant differences (Table 8), while the combination of those parameters pLDA was able to detect a significant gait difference.
Mean and SD of pLDA and several CatWalk parameters from vehicle moderate SCI and vehicle moderate-severe SCI models in Study 1 at 60 dpi and their p values as shown in Figure 3A
(a, b) Combinations of gait parameters (pLDA) and BBB scores for each group in Study 1 without treatment (mean ± SEM). (c, d) Treated and untreated pLDA and BBB scores for moderate contusions. (e, f) Treated and untreated pLDA and BBB scores for moderate-severe contusions. UI: uninjured; Veh.Mod.: vehicle moderate SCI; Veh.Mod-Sev.: vehicle moderate-severe SCI; Exp.Mod.: experimental moderate SCI; Exp.Mod-Sev.: experimental moderate severe SCI; Exp.Mod-Sev.: experimental moderate severe SCI; *p < 0.05; **p < 0.01; ####pcompared with UI < 0.0001.
Examination of the CatWalk parameter linear combination across various studies with similar lesion type and severity (Studies 2 and 3)
Following validation that our newly developed pLDA is capable of distinguishing thoracic SCI severity, we re-examined a previous SCI study (Study 2; Garcia-Ovejero et al., 2014) using a similar lesion type (i.e., thoracic contusion) and reporting differences in single CatWalk parameters between experimental treatment groups (hindpaw duty cycle, hindpaw swing, hindpaw BOS, phase dispersions, and RI). In this study, the regularity-index controlled BBB score also detected differences between the treatment groups. In addition, we also re-examined unpublished data from a study (Study 3) in which no significant differences in recovery were detected neither according to the BBB score nor in half of the single CatWalk parameters selected for the pLDA. First, we examined vehicle SCI rats from Study 2 (male Wistar rats) and Study 3 (female Fischer 344 rats), after Th8/9 contusion injury with a displacement over 1000 mm, which showed similar gait performance according to our pLDA (Fig. 4A). Correspondingly, the BBB scores of the vehicle SCI rats in these two studies showed that they were equally impaired and did not recover within a period of 60 d (Study 2), or 43 d (Study 3), respectively. Importantly, it should be noted that the use of different CatWalk systems (i.e., CatWalk 7.1 and CatWalk XT) for these two studies did not influence the comparability of data acquired here.
(a, b) Combinations of gait parameters (pLDA) and BBB scores for vehicle and young groups in Studies 2 and 3, respectively (mean ± SEM). (c-d) The Two-way ANOVA in Study 2 shows significant effects on pLDA of both groups (p = 0.037) and time (p < 0.00001), without significant interaction. In the RI-controlled BBB score (Study 2), significant effects of both groups (p = 0.003) and time (p < 0.00001) were observed, as well as significant interaction (p = 0.049). (e-f) The Two-way ANOVA in Study 3 shows significant effect on pLDA of both group (p < 0.00001) and time (p < 0.00001), with significant interaction (p < 0.001). In the BBB score (Study 3), significant effect of time (p < 0.00001) was observed, but no significant effect of group and interaction; *p < 0.05, ***p < 0.001, ****p < 0.0001 Prog.: progesterone; Veh.: vehicle.
Examination of the pLDA and BBB scores from Study 2 (Fig. 4C,D) confirmed that treatment with natural progesterone improved gait performance, as previously reported (Garcia-Ovejero et al., 2014). Differences in gait performance between vehicle and treated rats were observed to be significant at 60 dpi with the pLDA and at 30 and 60 dpi with the RI-controlled BBB score. The values of pLDA and nine single CatWalk parameters for the SCI and SCI+PROG at 60 dpi are shown in Table 9, as well as their p values calculated from t tests without assuming equal variances.
Mean and SD of pLDA and several CatWalk parameters from the SCI and SCI+PROG models in Study 2 at 60 dpi and their p values as shown in Figure 4C
In Study 3, calculation of pLDA revealed a difference in gait recovery between rats undergoing contusion SCI at a young age (three months) versus at an old age (20–24 months; Fig. 4E). In contrast, no differences in the BBB score were detected between the two age groups (Fig. 4F). It should be noted that the young and old rats started at similar pLDA scores before injury, regardless of the differences in age and bodyweight at the given time. The values of pLDA and nine single CatWalk parameters of the young and old groups at 29 dpi are shown in Table 10, as well as their p values calculated from t tests without assuming equal variances.
Comparison of the CatWalk parameter linear combination between distinct SCI lesion models (Studies 2, 4, and 5)
Contusion injuries are the most common type of SCI experienced in patients, and thus this type of injury is commonly used in rodent SCI studies. However, several transection SCI models have been developed to examine the deficit brought on from injury to specific spinal tracts. Therefore, we challenged our new pLDA to distinguish various SCI lesion models, such as a bilateral dorsal column transection and a dorsal hemisection model, in comparison to the contusion SCI models described above. The results from the uninjured rats and vehicle SCI rats at 30 dpi acquired in Study 2 (Th8 contusion), Study 4 (Th8/9 dorsal hemisection), and Study 5 (C4 bilateral dorsal column wire knife lesion) are shown in Figure 5A.
a, Combinations of gait parameters (pLDA) of uninjured rats and vehicle SCI rats in Study 2 (S2: Th8 contusion), Study 4 (S4: Th8/9 dorsal hemisection), and Study 5 (S5: C4 bilateral dorsal column lesion). b, Combinations of gait parameters (pLDA) of vehicle SCI rats in Studies 2 and 4. c, BBB scores of vehicle SCI rats in Studies 2 and 4; mean ± SEM; ***p < 0.001, ****p < 0.0001, ###pcompared with UI < 0.001, ####pcompared with UI < 0.0001; UI: uninjured; Veh.: vehicle; Cont.: contusion; DH.: dorsal hemisection.
In the rat thoracic contusion and dorsal hemisection SCI studies (Studies 2 and 4), the pLDA was able to distinguish the gait differences between uninjured and injured rats. It could also discriminate locomotor differences between the two SCI lesion models at 30 dpi. The values of the pLDA and nine single CatWalk parameters of the vehicle SCI animals from Studies 2 and 4 at 30 dpi are shown in Table 11, as well as their p value calculated from an ANOVA with a multiple comparison test by Bonferroni correction. The values of the pLDA and nine single CatWalk parameters of the uninjured and injured rats (30 dpi) from Study 5 are shown in Table 12. The single parameters obtained with the CatWalk system depicted different gait properties following C4 bilateral dorsal column lesion in rats (Study 5) compared with the Th8/9 dorsal hemisection (Study 4) or contusion SCI rats (Study 2). Among others, the injured rats in Study 5 showed significantly longer forepaw stride lengths compared with the uninjured rats, while the injured rats in Studies 1–4 mostly showed smaller forepaw stride lengths compared with the uninjured rats. Consequently, our parameter combination pLDA is not adequate to reflect the progression of gait recovery following C4 bilateral dorsal column transection in rats.
Mean and SD of pLDA and several CatWalk parameters from the young and old SCI models in Study 3 at 29 dpi and their p values as shown in Figure 4E
Mean and SD of pLDA and several CatWalk parameters from Studies 2 and 4 at 30 dpi and their p values (ANOVA and multiple comparison tests by Bonferroni correction) shown in Figure 5A
Mean and SD of pLDA and several CatWalk parameters from Study 5 and their p values (t test without assuming equal variances) shown in Figure 5A
To compare the progression of gait recovery in rats following different SCI types performed at the same spinal level, we compared the data acquired in Th8 contusion SCI (Study 2) with Th8 dorsal hemisection SCI [Study 4; pLDA (Fig. 5B) and BBB score (Fig. 5C)]. The pLDA revealed significant differences of walking function between respective injury models, displaying its enhanced sensitivity in differentiation. This was not detectable through BBB scoring. Again, we observed that despite the use of male and female Wistar rats in Studies 2 and 4, respectively, their pLDA before surgery were nearly identical, regardless of bodyweight and sex differences. It should be noted here that, similar to the bilateral dorsal column lesion used in Study 5, the dorsal hemisection performed in Study 4 also led to longer hindpaw stride lengths.
Discussion
Sensitive tools to analyze gait recovery after SCI are necessary to better predict, from preclinical work, the potential effectiveness of treatments in clinical interventions. Gait represents a complex movement pattern, which can be divided into numerous movement components that can only be measured with a combination of single parameters which are dynamically interlinked. Current gait analyses predominantly rely on the BBB score, an observer-based open field locomotor scoring system which incorporates some static gait parameters from the HLs exclusively. Because of the non-linearity of this scoring system, ceiling effects at 8 and 13 points are often observed. A score of 8 has rhythmic movement of joints or plantar placement without weight support. A score of 13 has frequent to consistent weight supported plantar stepping with frequent FL–HL coordination, but the paw is rotated outwards. Moreover, it was found that experimental SCI targeting the dorsal region of the spinal cord has less influence on the BBB score (Schucht et al., 2002).
Here, we provide a novel unbiased combination of SCI-specific gait parameters acquired from the automated CatWalk gait analysis program, based on LDA (pLDA). It is noteworthy that our pLDA could be validated using parameters measured in the frame of five studies realized in various international laboratories, performed on different rat strains and from both sexes, by diverse experimenters using distinctive lesion types at different spinal levels. To the best of our knowledge, the robustness of this pLDA is unique for the analysis of CatWalk data in the field of preclinical SCI. The pLDA generated through the combination of nine SCI-related CatWalk static and dynamic gait parameters could reliably detect differences in gait performance, that researchers may have visually observed, although analysis of individual gait parameters or BBB score alone may not have disclosed differences to them. Although single CatWalk parameter analysis may be used to determine treatment effects following SCI, they lack the comprehensive analysis on the overall status of locomotion recovery as well as the comparison and standardization across studies currently used by the BBB score. Our pLDA indicates damage to the spinal cord involved in gross overground locomotor control by displaying a significant lowering of the overall value (Fig. 6). Thus, we believe this pLDA represents a significant improvement of our capacity and sensitivity to monitor the recovery of locomotion after SCI and the impact of therapeutic interventions.
Schematic representation of the pLDA scores for varying lesion types and severities.
Our study demonstrates that the comprehensive CatWalk-derived measure pLDA yields results, which are not affected by confounding factors such as sex, size and strain. This is not the case with individual CatWalk-derived single parameters, which have been shown to be influenced by sex, size and strain (Heglund et al., 1974; Taylor, 1978; Webb et al., 2003; Koopmans et al., 2007; Neckel et al., 2013; Jacobs et al., 2014; Machado et al., 2015; Neckel, 2015; Datto et al., 2016; Crowley et al., 2018). However, our goal with the pLDA approach was to add more power to our machine learning approach by neutralizing all these various sources of confounding variations. Therefore, data generated in studies performed in different laboratories on rats of different sex, strain, age, lesion type and lesion level were used to establish and validate the pLDA. In addition, it would be interesting to further examine the potential of our pLDA across species, e.g., in mice with similar thoracic SCI lesions.
It is intriguing that half of the nine parameters characterizing at best the uninjured and SCI rats were forepaw parameters. In addition to a simple functional adaption and weight transfer of HL to the FL (Wilcox et al., 2017), the neuroanatomical role played by propriospinal interneurons connecting the cervical and lumbar enlargements should also be considered (Ghosh et al., 2010; Flynn et al., 2011; Bareyre et al., 2004). These propriospinal connections mostly run along the ventrolateral tracts and are therefore more affected by a contusion injury than a dorsal hemisection. This could be appreciated through the drastic differences in the respective forepaw parameters following both lesions (Table 11, comparing Th8 Contusion, Study 2, vs Th8 dorsal hemisection, Study 4). In contrast, the BBB score is focused explicitly on HL alterations and may be why it is not as sensitive in distinguishing differences between these two SCI types (Fig. 5C).
In rodents, the reticulospinal and vestibulospinal tracts, located in the lateral to ventrolateral white matter, as well as the rubrospinal tracts, located in the dorsolateral white matter, are major elements of the locomotor control (Raineteau et al., 2002; Ballermann and Fouad, 2006; May et al., 2017; Asboth et al., 2018; Witts and Murray, 2019). Their preservation/disruption will be decisive in the locomotor parameters. White matter sparing was reported with progesterone treatment (Study 2) and forepaw parameters were directly affected. Further examination of direct sparing of ventrolateral descending tracts following progesterone treatment would be warranted. Moreover, plasticity of remaining spinal tracts has been reported to be crucial for locomotor recovery following SCI in rodents (Bareyre et al., 2004, 2005). Accordingly, our pLDA could detect the significantly different locomotor recovery following contusion in young and aged rats (Study 3), even in the presence of similar white matter sparing. It may be speculated that differences during recovery seen were age-related consequences in the capacity to undergo neuronal plasticity.
Our unbiased gait recovery analysis through the combination of SCI-specific CatWalk parameters relies heavily on the capacity of rats to retain or regain HL weight support (BBB score ≥10) required to transverse the runway. As a consequence, this pLDA method is not suited for the early phase of locomotor recovery during which weight support is not achieved, or in the case of very severe SCI. It should be stated, however, that this pLDA was specifically designed for thoracic lesions. When applied for the analysis of a bilateral dorsal column lesion at cervical level C4, the pLDA was unable to discern locomotor deficits compared with uninjured controls. We cannot rule out that other types of cervical lesions may be adequately characterized by our pLDA and this issue should be further examined.
pLDA represents a comprehensive linear behavioral measure reflecting the integrity/dysfunction of walking in SCI rats regardless of sex, strain and weight. Its strength is demonstrated in thoracic contusion SCI and defined partial transection models (dorsal hemisection). Moreover, certain conditions (progesterone treatment, young age) clearly improve gait function measure with pLDA. Although further work is necessary to fully understand the underlying mechanisms that lead to these differences. The fact that our pLDA integrates multidimensional parameters defining the quality of locomotion (spatial, temporal, interlimb and variation of interlimb coordination), as compared with the unidimensional and sequential BBB scoring system or single CatWalk parameters, explains its higher sensitivity. Therefore, we provide access to our gait recovery pLDA at https://github.com/mad-lab-fau/GaitParamCombSCI (simple downloadable Excel file) and hope that it will be used in the field of preclinical SCI for ease of comparisons in the recovery and quality of gait between various treatment groups.
Acknowledgments
Acknowledgements: We thank Dr. Ruediger Rupp and Dr. Christopher Sliwinski for discussions and reading of this manuscript that improved the overall quality of the work.
Footnotes
The authors declare no competing financial interests.
I.K.T. was supported by the DAAD Research Grant 57129429. B.E. was supported by the German Research Foundation (DFG), the Heisenberg Professorship Program Grant ES 434/8-1. Work in Study 1 was supported by the Olympia Morata Program at Heidelberg University Medical Faculty for B.S. and R.P. as well as the Interdisciplinary Neurobehavioral Core Grant AZ42-04HV.MED(15)/6/1 in Heidelberg. L.B. and S.C.-D. were supported by the Austrian Agency for International Cooperation in Education and Research (OeAD-GmbH) Grant HR02/2018 and Paracelsus Medical University Research Fund Grants PMU-FFF R-13/05/054-GRA, E-15/22/113-COS, and R-19/04/124-BIG.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.