Multivariate description of gait changes in a mouse model of peripheral nerve injury and trauma

Abstract

Objective

Animal models of nerve injury are important for studying nerve injury and repair, particularly for interventions that cannot be studied in humans. However, the vast majority of gait analysis in animals has been limited to univariate analysis even though gait data is highly multi-dimensional. As a result, little is known about how various spatiotemporal components of the gait relate to each other in the context of peripheral nerve injury and trauma. We hypothesize that a multivariate characterization of gait will reveal relationships among spatiotemporal components of gait with biological relevance to peripheral nerve injury and trauma. We further hypothesize that legitimate relationships among said components will allow for more accurate classification among distinct gait phenotypes than if attempted with univariate analysis alone.

Methods

DigiGait data was collected of mice across groups representing increasing degrees of damage to the neuromusculoskeletal sequence of gait; that is (a) healthy controls, (b) nerve damage only via total nerve transection + reconnection of the femoral and sciatic nerves, and (c) nerve, muscle, and bone damage via total hind-limb transplantation. Multivariate relationships among the 30+ spatiotemporal measures were evaluated using exploratory factor analysis and forward feature selection to identify the features and latent factors that best described gait phenotypes. The identified features were then used to train classifier models and compared to a model trained with features identified using only univariate analysis.

Results

10–15 features relevant to describing gait in the context of increasing degrees of traumatic peripheral nerve injury were identified. Factor analysis uncovered relationships among the identified features and enabled the extrapolation of a set of latent factors that further described the distinct gait phenotypes. The latent factors tied to biological differences among the groups (e.g. alterations to the anatomical configuration of the limb due to transplantation or aberrant fine motor function due to peripheral nerve injury). Models trained using the identified features generated values that could be used to distinguish among pathophysiological states with high statistical significance (p < .001) and accuracy (>80%) as compared to univariate analysis alone.

Conclusion

This is the first performance evaluation of a multivariate approach to gait analysis and the first demonstration of superior performance as compared to univariate gait analysis in animals. It is also the first study to use multivariate statistics to characterize and distinguish among different gradations of gait deficit in animals. This study contributes a comprehensive, multivariate characterization pipeline for application in the study of any pathologies in which gait is a quantitative translational outcome metric.

Introduction

In the study of human gait, methods of statistical evaluation include multivariate processes of feature extraction and selection [1]. Animal models are valuable for studying conditions and treatments that are difficult to study and experiment with in humans (e.g. traumatic nerve injuries, limb replacement, limb transplantation, tissue regeneration strategies) [2].

More than 20 million people in the United States are estimated to have some form of peripheral nerve damage [3]. However, available treatments for peripheral neuropathies are limited to correcting underlying causes, promoting healthy lifestyle habits, controlling immune responses + inflammation with medication, engaging in physical therapy, using rehabilitative devices, receiving surgery, and/or adopting behavioral strategies to cope [3]. There are no available treatments that specifically aim to promote nerve repair and regeneration. Rodent models are versatile for assessing novel approaches to promote nerve repair ranging from tissue engineering constructs to cutting-edge surgical interventions.

There are widely adopted methods for evaluating success when studying nerve injury and repair. These include methods for assessing cellular appearance (histopathologic examination) [4], nerve conduction (nerve conduction velocity) [5], nerve-muscle connection (electromyography) [6], sensory function (Hargreaves) [7], and motor function (Rotarod, DigitGait, CatWalk, swimming, climbing) [8, 9].

Ultimately, the lower extremity’s prime function is to enable locomotion. In human gait analysis, studies have investigated multivariate techniques to describe relationships among measurable features of gait [1]. These relationships help statistically characterize gait phenotypes and are used to select features for measurement, classification, and prediction [1]. Implementing a similar level of statistical rigor to animal gait studies represents a novel problem requiring further refinement and investigation.

In rodents, the study of locomotion via treadmill gait (DigiGait™) or free ambulating (CatWalk™) systems provides high levels of spatial and temporal detail, measuring over 30+ individual spatial or temporal gait parameters [10–14]. Analyses of data measured by these systems are largely limited to traditional, univariate, multiple-hypothesis testing or feature selection [15–18]. Usually this involves comparing the means of individual output parameters for statistical significance and applying some correction for multiple-hypothesis testing [19, 20]. Such univariate study enables conclusions that do not account for possible relationships among individual gait parameters. A machine learning model can also be trained using feature selection techniques that omit the examination of multivariate relationships between features [21]. Thus, discoveries in rodent models, to date, have been primarily limited to a portion of the interplay that makes up the animal’s overall gait.

These limited, univariate gait studies span a variety of rodent models of gait deficit including those of metabolic (diabetic neuropathy) [22], degenerative (ALS) [23], traumatic (nerve transection / limb transplantation) [24], congenital [25], and other etiologies. While these diseases may result in dysfunctional gait, it is unlikely that the spatiotemporal details of their dysfunction are the same. For locomotion, indices such as the Sciatic Function Index (SFI) have been developed via multivariate linear regression models to compare healthy and dysfunctional gait [26]. However, the SFI is a dimensionless number and does not reveal additional insight into exact relationships between specific spatial and temporal measures [27]. In short, what does it really mean if the SFI index is statistically different between two gait states? How do we describe more than just a “difference” between two states?

A comprehensive description of the rodent’s gait would require a methodology that can factor relationships among all 30+ measures and prioritize those that contribute to gait phenomena rather than noise. We hypothesize that out of all spatiotemporal measures of rodent gait, there exists a discrete subset sufficient to uniquely describe the gait of three distinct physiological states: (1) healthy animals, (2) those with nerve injury due to a complete nerve transection and re-attachment, and (3) those with nerve, muscle, and bone injury due to total limb transplantation. Furthermore, we hypothesize that multivariate characterization will reveal relationships that encode and describe the different gait phenotypes in a biologically consistent manner unlike univariate analysis. Finally, we hypothesize that the revealed relationships will enable the training of a classifier model that can discriminate between gait states with higher accuracy than models trained with features identified via univariate analysis.

Methods

Animals and experimental groups

43 male B6 mice, 8 to 10 weeks old (20g), were obtained from The Jackson Laboratory (Bar Harbor, ME, USA) and were group-housed. Male mice were used due to their larger size, thereby having easier vascular access.

There were three groups: 17 animals in the negative control group received no treatment. 12 animals in the experimental group received only nerve damage (neurorrhaphy) i.e. a complete nerve transection of the femoral and sciatic nerves with re-coaptation. 14 animals in the second experimental group received a total hind-limb transplant. Animals receiving a procedure received only one procedure on their left hind limb.

The motivation for assessing these three groups was to establish a gradient of increasing neuromuscular damage. In doing so, the goal was to characterize how spatio-temporal gait parameters were altered accordingly. Thereby, investigating how musculoskeletal damage, in addition to nerve injury, alters the statistical characterization compared to just the contribution of nerve injury alone. In other words, groups were chosen to represent increasing degrees of damage to the neuromusculoskeletal sequence of healthy gait.

Additionally, while clinical experience to date suggests it is unlikely that clinical lower limb transplantation will become common, the hind-limb transplantation model was still used in this study for four main reasons: (1) as part of an increasing gradient of damage on which to study the utility of multivariate characterization to identifying biologically relevant features and relationships, (2) to investigate the value of multivariate characterization as a novel and proposed new standard for gait analysis in a therapeutic intervention and model that is difficult, albeit, impossible to study in humans, (3) to establish and characterize this model and statistical pipeline as one of potential benefit to research in functional nerve recovery for upper extremity transplantation, and (4) to maximize the degree of damage and thereby spatiotemporal changes to the limb in the interest of evaluating this new method of gait characterization at its limits.

All experimental procedures were conducted according to an IACUC-approved protocol (protocol IS00003195) and were conducted in accordance with the National Research Council’s Guide for the Care and Use of Laboratory Animals. All surgery was performed with animals kept on a self-regulating heating pad and under isofluorane anesthesia. The depth of anesthesia was examined before surgery with a toe pinch response and subsequent checks at 15-minute intervals during the surgery. Post-operatively, the animals were kept in their cages, with a portion of the cage on a warm pad. The animals were checked on every hour for the first 3 post-op hours and then 4 times a day for two days for any signs of surgical complications. To minimize suffering, pain was managed via local anesthetics (bupivacaine or lidocaine) at the surgical site with a multi-modal analgesic regimen of Buprenorphine-ER and meloxicam to last an additional 3 days. At the conclusion of all experimental timepoints animals were sacrificed via Isofluorane euthanasia with confirmatory bilateral thoracotomy.

Mouse hind-limb transplant

The murine hind-limb transplantation was performed using the surgical technique modified from the one previously described by our lab [24] and Furtmuller et al [28]. The donor procedure took ~45 minutes per procedure and the recipient procedure took ~2 hours per procedure. The isogenic transplants required no immunosuppressive regimen. Briefly, for donor right hind limb harvest, after adequate anesthesia, the right hind limb was shaved and prepped with 70% ethanol plus povidone iodine, followed by a circumferential incision in the mid-thigh. The femoral neurovascular bundle was carefully dissected. The nerve, artery, and vein were separated and were then transected distally at the level of the superficial epigastric artery. The ventral and dorsal muscle groups were divided proximally to the mid-thigh, thereby separating the remainder of the hind limb from the animal’s body. The hind-limb graft was perfused with heparinized Ringer’s solution and stored in a Ringer’s solution-containing dish at 4°C. For the recipient limb implantation, the recipient native right hind limb amputation was performed similarly to the donor right hind limb harvest. The donor limb graft was anatomically aligned and the recipient’s femur bone was connected with the donor femur bone using an intramedullary stainless-steel rod. Muscle coaptation was accomplished using 6–0 Polysorb suture material. Vessel anastomoses were established using the cuff technique. Continuity of the sciatic and femoral nerves was established with 11–0 sutures. Then, the skin was closed using 6–0 monofilament sutures (Ethilon). After the surgery, the animal was then regularly monitored for at least 4 hours before returning to the housing facility.

Mouse model of sciatic and femoral nerve transection injury

The sciatic and femoral nerve transection model took ~30 min. to complete. The mice were anesthetized and prepped as described in the previous section. An incision of approximately 8 mm in length in the right posterolateral hind limb was made. The intermuscular septum was gently dissected. The sciatic nerve was exposed and separated from its sheath before being cut using sharp micro scissors. The two nerve stumps were carefully repositioned, and the epineurium was connected using an 11–0 suture. No retraction or torsion was observed, and the stumps fit well without any gaps between them. Similar procedures were performed on the femoral nerve. Subsequently, the skin was suture closed as described in the previous section. The animal was placed under a heating lamp and allowed to recover in its cage. Regular monitoring was conducted for at least 4 hours before the animal was returned to the housing facility.

Digigait procedure

After the respective procedures required for each experimental group, the animals were allowed to recover and reacclimate for two weeks before beginning data collection on the treadmill system. The mice were kept in an animal holding room and acclimated in the Behavioral Phenotyping Core before experiments. Gait was measured using the DigiGait™ Imaging System (Mouse Specifics Inc.) at three different speeds (10, 17, and 24 cm/s). Footprints from a 3–4 second video clip were analyzed for each speed using Digigait™ Analysis version 15. We used the software to quantify 30+ gait indices.

Statistical analysis

Dimensionality reduction.

All analysis was conducted in MATLAB. Exploratory analysis to reduce and identify the primary dimensions responsible for describing rodent gait was conducted via feature selection and factor analysis techniques. Feature selection was conducted using two contrasting approaches: (a) a traditional univariate analysis with Bonferroni correction and (b) a multivariate, forward sequential feature selection with cross-validation [29].

Feature selection.

In this study we implement a simple filter by applying univariate criteria separately on each feature with Bonferroni correction. This was done as a baseline to replicate typical (i.e. non-multivariate) analysis of gait performed across the field. This baseline is important for evaluating the hypothesis that multivariate characterization will identify feature sets and feature relationships that better describe gait states and will lead to more accurate classification. We also apply a multivariate, forward sequential feature selection in a wrapper fashion to find important features with the goal of minimizing misclassification error (MCE) of our learning algorithm. This was done to explore the hypothesis that multivariate characterization will reveal relationships that encode and describe the different gait phenotypes in a biologically consistent manner unlike univariate analysis; and that the revealed relationships will enable the training of a classifier model that can discriminate between gait states with higher accuracy than models trained on features identified via univariate analysis. Embedded selection was then applied during model training to further confirm the ideal feature set. During the feature selection procedure we applied 10-fold cross-validation to the training set. We ensured that the same animal was not included in both the training and test sets of each individual fold. While feature selection identified the measured parameters to use in the model it did not describe how the selected parameters may relate to each other.

Factor analysis.

To learn more about the exact relationships among features, factor analysis was conducted. Understanding the exact relationships among features is important to exploring the hypotheses that said relationships encode and describe the different gait phenotypes in a biologically consistent manner; and that the revealed relationships will enable the training of a classifier model that can discriminate between gait states with higher accuracy than models trained with features identified via univariate analysis. In the factor analysis model measured variables are dependent upon a smaller number of unobserved or latent factors [30]. The coefficients on these latent factors are called “loadings”. An added component to account for noise called “specific variance” is included. Average specific variance and factor loading as a function of number of latent factors was plotted to determine the number of factors that would be reasonable to assume for analysis. A loading cutoff greater than 0.6 or less than -0.6 was then used to identify the important features within each factor grouping.

Machine learning pipeline: Cross validation, model training, and accuracy calculation.

Definitions within study context. Machine learning: computational statistical algorithms that learn from data labeled as coming from healthy, nerve transected, or limb transplanted mice and using those learnings to understand the features and relationships that encode patterns specific to those respective groups to enable the training of models that can be used for various purposes like classification.

Training of classifier models: computational statistical model development via machine learning methods using data labeled from healthy, nerve transected, or limb transplanted mice. The classifier models are encoded by features and relationships identified by the machine learning methods. Once trained, new, unlabeled data from animals of uncertain gait quality can be input into the model. The model will then return numerical values representing the likelihood that the input data is characteristic of a healthy, nerve-transected, or limb-transplanted animal and its gait. The output values can be used to classify animals with known identities into respective groups. They can also be used to quantify the quality of the statistical characterization conducted on each gait phenotype.

Cross validation. Training and testing sets were defined by randomly picking 80% as training and 20% as testing data. 10 randomly selected training and testing sets were defined to perform 10-fold cross validation. We ensured that the same animal was not included in both the training and test sets within a respective fold.

Model training and performance evaluation. Various classifier model architectures were evaluated using the features identified from the feature selection and dimensional reduction techniques described above. Based on the feature selection results 16 features were chosen to be included in model training. The evaluated architectures included discriminant analyses, random forest, and support vector machine. Regardless of the model type, training was performed using the results of the multivariate feature selection process described above. Outputs of all models were either 0 = dysfunctional gait (surgery) or 1 = healthy gait (control). Average model accuracy, precision, recall, and F-score were measured at the conclusion of the cross-validated model training process. Models trained via multivariate feature selection were compared for their performance with a model trained via univariate feature selection to explore the hypothesis that a classifier model trained on features identified via multivariate analysis can discriminate between gait states with higher accuracy than models trained with features identified via univariate analysis.

Results

Surgical results

A total of 30 surgeries were performed in this study with a high success rate (87%). The experimental group receiving neurorraphy contained 14 animals, the group receiving total hind-limb transplant contained 12 animals, and the control group contained 17 animals.

Attention to the vascular anastomoses, neural repair, and bone coaptation was crucial to maximizing the viability of the limb transplant model. The isografted mice lived for the course of the study and none of them showed signs of rejection. The peripheral nerve transection model involved femoral and sciatic nerve transection and coaptation only. All mice in this group also lived for the course of the study. Both models used a nerve cuff to hold the coapted nerves in place. Operative design and representative results are shown in Fig 1.

[Figure omitted. See PDF.]

(A) The donor’s femoral neurovascular bundle was carefully dissected. The femoral nerve, artery, and vein were carefully separated and transected distally at the level of the superficial epigastric artery. (B) The donor limb graft was anatomically aligned such that the recipient’s femur bone was connected with the donor femur bone using an intramedullary stainless-steel rod, vessel anastomoses were done with the cuff shown in panel A, and the sciatic and femoral nerves were connected with 11–0 sutures. (C) Panel C shows a mouse with a completed hind-limb transplantation. After the procedure the mice were allowed to recover and acclimate for two weeks before proceeding with DigiGait evaluation. The peripheral nerve transection model involved total transection of the femoral and sciatic nerves in the same way it was done in the hind-limb transplantation model.

Feature selection

Univariate analysis.

Traditional gait studies conduct univariate analysis, which was done here as a basis for comparison. We first plotted the empirical cumulative distribution function (CDF) of the p-values. The CDF of the p-values visualizes the percentage of feature means, when compared for statistical difference via t-test that lie under the critical value. Applying Bonferroni correction for multiple hypothesis testing determined a critical p-value of .001. In comparing pathological gait from limb transplantation to healthy controls, 28% of features had statistical significance (Fig 2), whicih amounts to 9 parameters (Table 1). In comparing pathological gait from nerve transection injury alone, 44% of features had statistical significance (Fig 2), which amounts to 14 parameters (Table 1). Comparing the means of every feature 1:1 with the control and then comparing the set of statistically significant features between the two degrees of gait deficit revealed a high degree of shared features with statistical significance (Table 1). However, via univariate analysis we learn little about how the identified features relate to each other to describe the respective gait states.

[Figure omitted. See PDF.]

Cumulative distribution function showing (A) the percent of features that are under a significant p-value when comparing their means between healthy versus transplanted animals with t-tests. Boneferroni correction determines a cutoff of 0.05 / 32 = 0.001. 28% of features have p-values below the .001 threshold and (B) the percent of features that are under a significant p-value when comparing their means between healthy versus nerve transection injury animals with t-tests. Boneferroni correction determines a cutoff of 0.05 / 32 = 0.001. 44% of features have p-values below the .001 threshold. The features below threshold in (a) and (b) are identified in columns 1 and 2 of Table 1 are establish a baseline of features selected univariately for comparison to those identified via multivariate selection. The value in multivariate feature selection will be explored quantitatively in subsequent figures.

[Figure omitted. See PDF.]

Univariate filter feature selection.

Measuring and plotting misclassification error (MCE) as a function of an increasing number included features is a univariate approach to determining the ideal number of features to include in model training. To illustrate, we compute MCE for a discriminant analysis model between 2 to 32 features and plot MCE accordingly (Fig 3). Note, in this simple, univariate feature addition approach, using 30 features minimizes MCE. Ultimately, the smallest MCE achieved with this method was 0.17 indicating a potential model accuracy of 83% in discriminating healthy gait from pathological gait due to limb transplantation. However, this is a model trained on one holdout set with likely overfitting from using 30 of the 32 available features. It is likely that this is modeling noise and it is unclear how much generalizable gait phenomena is being modeled. Thus, multivariate (wrapper) methods for feature selection were explored.

[Figure omitted. See PDF.]

Top: (A) Misclassification error (MCE) as a function of number of features calculated using a pseudo-quadratic discriminant analysis, holdout set, and a simple filter method that does not account for any interactions. The lowest achievable MCE was around 0.16 when including 30 features. Bottom: Multivariate, forward feature selection with cross-validation to minimize misclassification error when classifying healthy gait from pathological gait was conducted to find the optimal combination of features that minimizes MCE. This approach takes interactions between features into account and uses a pseudoquadratic discriminant analysis. The lowest MCE was (B) ~0.12 using 10–14 features when classifying limb transplantation injury and (C) ~0.12 using 10–14 features when classifying just nerve transection injury. Notice how the MCE has a much more consistent and predictable pattern as features are added in comparison to the univariate, simple filter approach referenced in (A).

Multivariate wrapper feature selection with cross validation.

The univariate approach to filter feature selection described above omits relationships between features and is representative of how most rodent gait studies have evaluated their data. Moreover, choosing features purely based on their ranking in filter selection can result in the modeling of redundant information. For example, features 2 and 8 have a linear correlation coefficient of 1.0. This is because these two features, are by definition, related (%SwingStride = 1 - %StanceStride).

Sequential feature selection is a multivariate approach to selecting a subset of features that works by sequentially adding features in various relationships to the model and monitoring exactly which feature relationships lead to decreases in MCE until it is minimized. When classifying between healthy gait and pathological gait due to limb transplantation, our results show that 10 to 12 features are optimal for minimizing MCE (smallest MCE achieved = 0.12 (Fig 3). For classifying between healthy and pathological gait due to nerve transection injury, our results show that 16 features achieve the minimum possible cross-validation MCE of 0.10. This is a 5–7% improvement in MCE achieved using ~20 fewer features than the simple filter approach above. Note how MCE increases beyond the critical numbers of features indicating overfitting and supporting the following hypotheses: (a) there is a critical subset of features that more accurately describes healthy versus pathological gait and (b) the subset identified via multivariate characterization more generally describes gait when compared to the subset identified via univariate characterization as evidenced by superior model performance and less overfitting. The selected features are shown in Table 1.

However, when examining the selected parameters, it is unclear how they may relate to each other. Despite the compelling classification results from adding multivariate feature selection, the process does not reveal exactly how these features relate to each. A technique like factor analysis is needed to identify precisely which features relate to each other.

Factor analysis

A prior gait study proposed a multivariate rodent gait model to represent groupings of variables that related to each other and better describe gait. Lambert et al. proposed that the multi-dimensional DigiGait data in rats with olivocerebellar ataxia can be reduced to three uncorrelated groupings of variables, or common factors, termed rhythmicity, thrust, and contact [31]. Similarly, factor analysis allows us to understand how the directly measured features (i.e. the collected data) compose potential latent (or unmeasurable) factors that more directly describe the phenomena of interest. To determine the number of likely factors, the average specific variation and loading were plotted as a function of factor number (Fig 4).

[Figure omitted. See PDF.]

The red line represents average specific variation and the blue line represents average loading (top two panels). At around 6 factors is when average loading plateaus and when average specific variation is at its inflection point. Thus, we examined the first 6 factors for their constituent features (bottom three panels). Using a cutoff of +/- 0.6, loadings were examined to identify feature groupings that composed the 6 most relevant factors. Features with loadings of greater than 0.6 or less than -0.6 were marked in bold in the bottom panels.

In Fig 4, the red line represents the average specific variation, and the blue line represents the average loading. The results confirm that at 6 latent factors, the loading does not increase, while the specific variation continues to decrease, which indicates overfitting. Maximal loading is viewed using 6 factors, while reaching an inflection point in minimizing specific variation. Thus, 6 factors were used to continue the comparative analysis between the latent factors that characterize pathological gait due to transplantation versus nerve transection injury (Table 2). Fig 4‘s bottom-right two panels show the calculated feature loadings of greater than 0.6 or less than -0.6, which was identified as a significant feature in the respective factor grouping. Note that Stance is the only feature shared between two factors. Otherwise, the factors are each composed of unique features.

[Figure omitted. See PDF.]

Table 2 shows the latent factors and the features they comprise for pathological gait due to limb transplantation. Table 2 also shows the latent factors and the features they are composed of for pathological gait due to nerve transection injury. The bolded latent factors in each sub-table (Factor 4 and Factor 6 respectively) are composed of different features between the two groups. The other 5 latent factors are composed of the same (or highly similar) features. The distinction is that pathological gait due to limb transplantation is distinguished by the relationship between the absolute paw angle and midline distance. Pathological gait due to nerve transection injury is distinguished by the relationship between paw area at peak stance in sq. cm, paw area variability at peak stance in sq. cm, and MIN dA/dT. The high degree of overlap of feature groupings between the two pathologies suggests that much of the statistical description of gait in these two sources of traumatic nerve injury is the same with a significant point of difference being factors 4 and 6 respectively (i.e. the two bolded columns) (Table 2).

It can be suggested that animals with limb transplantation have more variation in the anatomy of their transplanted paw in relation to the midline of their bodies. In contrast, animals that received nerve transection injury and reattachment do not demonstrate this same variation. This is consistent with the underlying biology: animals with a total limb transplantation likely have an alteration in the anatomical configuration of their paw related to their midline while animals receiving only nerve transection injury might not. Additionally, the data suggests that animals with exclusively nerve transection injury have a discernible variation in how the contact of their paw relates to their ability to make fine motor movements (i.e. jerks or dA/dT) while animals receiving a total limb transplant are likely to have this masked by gross motor defects related to the extensive musculoskeletal injury that accompanies the nerve injury from the transplant.

Training a machine learning classifier to discriminate between gait phenotypes

Using the features identified via multivariate feature selection and confirmed by factor analysis we assessed four different model architectures (Table 3). Accuracy in distinguishing between healthy gait and pathological gait due to any form of peripheral nerve injury (Table 3A) ranged from 0.75–0.91, precision from 0.78–0.93, recall from 0.77–0.91, and F-Score from 0.77–0.92. The ensemble model (i.e. boosted classification trees), had the highest-performing metrics in distinguishing between healthy gait and pathological gait due to any form of peripheral nerve injury.

[Figure omitted. See PDF.]

Accuracy in distinguishing between types of peripheral nerve injury (i.e., hind-limb transplant vs. total nerve transection) ranged from 0.64–0.88, precision from 0.71–0.93, recall from 0.69–0.88, and F-score from 0.69–0.90. The ensemble model again was the highest-performing architecture in distinguishing between different types of peripheral nerve injury.

Values outputted by the model on the holdout set were graphed as a box and whisker plot. A strong statistical difference in output distribution between healthy versus injured animals as well as between pathological states was observed (p < .0001) (Fig 5). Sample confusion matrices from one iteration were calculated and included as well (S1 and S2 Figs in S1 File).

[Figure omitted. See PDF.]

The discriminant analysis model was used to generate values for animals with various phenotypes and compare (A) healthy gait to dysfunctional gait due to limb transplantation (B), healthy gait to dysfunctional gait due to nerve transection alone, and (C) dysfunctional gait due to nerve transection from dysfunctional gait due to limb transplantation. Independent, two-tailed t-tests were conducted with correction for multiple hypothesis testing to reveal high statistical significance among groups (p < 0.0001).

Discussion

The goal of this study was to use multivariate statistics to learn about the features, relationships, and factors that are most descriptive of gait in different pathological states. In doing so, we explored the hypothesis that a multivariate characterization of gait will reveal novel relationships between spatiotemporal components specifically in the context of peripheral nerve injury and trauma. Gait data from DigiGait’s system includes 30+ different spatiotemporal features or outputs. It is important to understand which features are most descriptive of specific pathological states. We called these features primary dimensions and defined them as those features most relevant to describing relationships within the gait data. We hypothesized that these relationships were deterministic of specific phenotypes, which contrasts with features that are not related to describing specific phenotypes and may conversely model noise.

Modeling noise can result in overfitting that detracts from statistical learning of the patterns that describe the phenotype of interest. Feature selection algorithms enable one to test various feature subsets from the original feature set and learn which subset is most descriptive of the underlying phenomena. Feature selection algorithms can be roughly grouped into three categories: filter, wrapper, and embedded methods [29]. Univariate filter methods select the feature subsets via multiple hypothesis testing without involving a learning algorithm [32]. Multivariate wrapper methods iteratively evaluate a variety of relationships between features and use the performance of a learning algorithm to evaluate each candidate feature subset [33]. Embedded methods determine feature importance as part of a model training process [34].

In Fig 3 the wrapper-based, multivariate feature selection resulted in the usage of 10 to 15 select features in the model and thereby decreased the MCE of the model by more than 20%. This result supports the hypothesis that accounting for relationships between components of the gait would allow for more accurate classification between distinct gait phenotypes than univariate studies alone (Table 2).

This accounting for multivariate relationships is largely absent from the field of animal gait analysis. Examining the features selected in a multivariate manner in Table 2 suggests that the two pathologies differ in their relationship to jerk (dA/dT). Nerve transection injury was better defined by minimum jerk (e.g. fine motor movements) while limb transplant injury was better defined by maximum jerk (e.g. gross motor movements). Clinically, the ability to identify characteristic measures of different degrees of injury could be relevant to evaluation, monitoring, and rehabilitation. It may inform the development of sensors capable of more precisely monitoring parameters like jerk. It may also help suggest quantitatively discernible differences in underlying (patho)physiology for further evaluation and focus.

Further exploring the potential importance of the single feature of jerk in distinguishing between these two gait phenotypes reveals that minimum vs. maximum jerk may be of particular importance (Table 2). While the two gradations of injury share many features as shown in Table 2, they are distinguished by the presence of minimum vs. maximum jerk. This difference suggests a hypothesis to investigate further: that nerve transection injury’s single source of variation to jerk is the nerve injury itself, resulting in a more discernible impact on fine motor movements (i.e. minimum dA/dT), whereas limb transplant injury, being a product of nerve + muscle + bone injury, contributes to a wider set variation and deficiencies in jerk, resulting in a more discernible impact to gross motor movements (i.e. maximum dA/dT). Translating this idea to a clinical setting, the neurologic physical exam includes components for evaluating fine motor and gross motor components of the gait [35]. However, the typical exam is at best semi-quantitative and does not measure jerk in a manner that would lend itself to statistical evaluation like the assessments performed in this study [35]. There exist quantitative methods for measuring changes in locomotive jerk in humans [36]. Perhaps there is value to doing so when evaluating the degree of peripheral nerve injury in settings of new presenting trauma or in monitoring the recovery of gait over time in response to trauma.

The ranking of the features in Table 2 may provide additional insight. It is listed in rank order of features selected for their proportion of contribution to characterizing each pathology. It is interesting to see the high degree of overlap in ranking between the nerve transection and limb transplant. Additionally, it is interesting to see how the features have some natural groupings. Measures of phase are primary, then secondarily measures of the paw, and tertiarily measures of symmetry. Also tertiary are measures of fine and gross motor movements respectively. While the measures of phase and the paw could be mutually characteristic of healthy gait and/or both etiologies of traumatic injury, we see that the tertiary measures are where significant difference between the two gait phenotypes may be resolved. The level of exploratory statistical description in the prior few paragraphs is largely absent from rodent studies across the field, including some of the latest that employ gait analysis [18, 25, 37]. Thus, we hypothesize that multivariate statistics has generalizable potential in reducing observational bias to current pre-clinical gait research protocols in the scientific community.

Towards this hypothesis, a recent study contributes evidence toward using gait as a quantitative translational outcome metric for therapeutic development in Angelman syndrome and other genetic neurodevelopmental syndromes [25]. This study contributes a comprehensive, multivariate characterization pipeline for application in the study of any pathologies in which gait is a quantitative translational outcome metric. Another study was comparable in its level of statistical description [31], coincidentally doing so in a genetic neurodevelopmental model of cerebellar ataxia [31]. Lambert et al. supported the value of multivariate characterization of gait by conducting feature extraction on highly dimensional gait data from a central nerve injury model of cerebellar ataxia. The authors hypothesized certain indirect, common factors that characterize gait and underlie the more directly measurable features. These common factors included thrust, rhythmicity, and contact area and were directly useful in discriminating between animals with a central lesion vs. those without [31]. In this study, we observed potential patterns in the selected features and hypothesized that pathological gait due to traumatic nerve injury could also be characterized by certain latent factors that are determined by and composed of the selected features.

The identified latent factors fell into similar categories as the ones identified by sequential feature selection. These were phasic parameters that represent either absolute or relative measures of the phases of each stride, symmetric parameters representing some measure of the gait symmetry, and regional parameters, which represent some measure of features specific to certain anatomical regions of the limb (i.e. the paw). The common factors and their feature compositions suggest that rodent gait is statistically characterized by a careful coordination between sub-components of the limb (e.g. specific anatomical parameters), the limb itself (e.g. phasic parameters), and the torso/remainder of the animal (e.g. symmetric parameters).

When contrasting the factors that are statistically meaningful to describing both types of traumatic nerve injury, we learn that the relationship between stride length and speed is important to both injury models (Table 2). Moreover, variation in the length of the stride is also important. Measures of the paw are also important factors in both models. However, the specific measures of the paw are consistent with biological differences in the respective pathologies. For example, limb transplantation injury was better characterized by the relationship of the paw to the midline, which may be due to loss of mechanical integrity and thereby positioning of the limb in relation to the area surrounding the transplant itself. In contrast, the nerve transection injury model was better characterized by fine motor movements of the paw (i.e. MIN dA/dT). This suggests that finer measures are more discernible when there is only nerve injury as compared to nerve + muscle + bone injury (Table 2). The results from factor analysis were consistent with the results from sequential feature selection and helped validate that the two gradations of peripheral nerve injury shared numerous important features. The results also displayed key differences related to fine vs. gross motor movements in each model respectively. Factor analysis allowed for a deeper understanding of the exact relationships among features in contrast to a simple rank list of features as provided by many feature selection techniques. This level of statistical description is currently lacking from gait studies in animal models. The many insights revealed support the hypothesis that multivariate gait analysis provides a more comprehensive statistical description of individual gait states than univariate analysis. In the clinic, employing such factor analysis techniques could similarly provide a window into better describing nuances of different patients’ gaits. This study offers a methodology that may be applied to existing DigiGait datasets to conduct similar exploratory multivariate analysis to reduce observational bias from traditional univariate methods. Evaluating the generalizability of this methodology on other etiologies of gait deficit is the subject of a subsequent study.

Lastly, we explored the hypothesis that characteristic feature relationships of distinct gait phenotypes could be used to train a classifier capable of distinguishing between said phenotypes. The trained classifier distinguishing individual animals for whether their gait was more likely healthy or pathological and discriminated between the two pathologies with high accuracy (Table 3). The model outputted a distribution of values with highly significant difference in their means (p < .0001) (Fig 5). The model was built by encoding the features and relationships uncovered to be most important in describing the gait pathologies respectively. Via the outputted values, we statistically distinguished between the two different traumatic pathologies with high statistical significance (p<0.0001). This is the first system to multivariately describe and distinguish between gradations of limb injury in a murine gait model. As a new methodology within rodent gait research protocols, one limitation of this study is its limited scope. On one hand, the multivariate statistics demonstrated in this study offer a convenient way for researchers to conceptualize meaningful groupings of features. On the other hand, clear biological tie-ins and strong statistical significance in classifier performance when applied to gait research protocols across the scientific community stands to be evaluated. We have identified datasets of DigiGait data available at our Behavioral Phenotyping Core from investigators studying central causes of gait disorder (e.g. stroke, transient ischemic attack) for us to apply the pipeline to and are pursuing said inquiry. Academic researchers can do the same by accessing the code available at our GitHub repository and citing this manuscript should any code be used towards future publications (https://github.com/luoyuanlab/gait). Industry professionals must reach out to the corresponding author to make use of the code.

Conclusion

There is no precedent for using multivariate statistics to characterize and distinguish between different etiologies of gait deficit. According to a PubMed search, hundreds of gait studies in animals have been published over the past two decades. The majority of them being univariate studies. This investigation paves the way for future studies of gait pathology across a variety of etiologies to conduct similar characterization. The statistical pipeline developed in this study represents an application of state-of-the-art methods to multivariately describe and quantify functional gait outcomes in a biologically consistent manner to a gradation of increasing injury. Development of this method may allow for more detailed characterization of functional outcomes in response to various therapeutic strategies of repairing peripheral nerve injury. There is an opportunity for method developers to contribute to time series characterization and classification. However, a barrier to time series study is the lack of an animal model of gait deficit that is confirmed to heal and recover a statistically significant amount of gait function to use as a positive control. More broadly, there is also an opportunity for methods developers to apply the pipeline developed here in the study of other causes of gait deficit (e.g. central, metabolic, vascular, degenerative, (auto)immune, genetic) and evaluate this method’s generalizability. The means to do so being offered in this study.

Supporting information

S1 File. Confusion Matrices of Varying Classifier Architectures (S1 Fig. Peripheral Injury vs. Control and S2 Fig. Nerve Transection vs. Limb Transplant).

All confusion matrices are from a single iteration of ten randomly selected training-testing splits. Performance metrics reported in the manuscript are the average of those ten.

https://doi.org/10.1371/journal.pone.0312415.s001

(DOCX)

Acknowledgments

The mouse hind-limb transplantation and the peripheral nerve injury surgical procedures were performed by the “Microsurgery & Preclinical Research Core" at Northwestern University Comprehensive Transplant Center. Gait analysis was done in the Behavioral Phenotyping Core at Northwestern University.

References

1. 1. Chau T. A review of analytical techniques for gait data. Part 1: fuzzy, statistical and fractal methods. Gait & posture. 2001;13(1):49–66. pmid:11166554

* View Article

* PubMed/NCBI

* Google Scholar

2. 2. Gorantla VS, Zor F, Nasir S, Breidenbach WC, Davis MR. Lower Extremity Transplantation: Concepts, Challenges, and Controversies. In: Tepe V, Peterson CM, editors. Full Stride: Advancing the State of the Art in Lower Extremity Gait Systems. New York, NY: Springer New York; 2017. p. 195–212.

3. 3. Peripheral Neuropathy. In: Stroke NIoNDa, editor. NIH.gov: National Institutes of Health; 2024.

4. 4. Geuna S, Raimondo S, Ronchi G, Di Scipio F, Tos P, Czaja K, et al. Chapter 3: Histology of the peripheral nerve and changes occurring during nerve regeneration. Int Rev Neurobiol. 2009;87:27–46. Epub 2009/08/18. pmid:19682632.

* View Article

* PubMed/NCBI

* Google Scholar

5. 5. Hodes R, Larrabee MG, German W. The human electromyogram in response to nerve stimulation and the conduction velocity of motor axons; studies on normal and on injured peripheral nerves. Arch Neurol Psychiatry. 1948;60(4):340–65. Epub 1948/10/01. pmid:18111214.

* View Article

* PubMed/NCBI

* Google Scholar

6. 6. Quan D, Bird SJ. Nerve conduction studies and electromyography in the evaluation of peripheral nerve injuries. Univ Pa Orthop J. 1999;12:45–51.

* View Article

* Google Scholar

7. 7. Chato-Astrain J, Philips C, Campos F, Durand-Herrera D, Garcia-Garcia OD, Roosens A, et al. Detergent-based decellularized peripheral nerve allografts: An in vivo preclinical study in the rat sciatic nerve injury model. J Tissue Eng Regen Med. 2020;14(6):789–806. Epub 2020/04/16. pmid:32293801.

* View Article

* PubMed/NCBI

* Google Scholar

8. 8. Navarro X. Functional evaluation of peripheral nerve regeneration and target reinnervation in animal models: a critical overview. Eur J Neurosci. 2016;43(3):271–86. Epub 2015/08/01. pmid:26228942.

* View Article

* PubMed/NCBI

* Google Scholar

9. 9. Lakes EH, Allen KD. Gait analysis methods for rodent models of arthritic disorders: reviews and recommendations. Osteoarthritis and cartilage. 2016;24(11):1837–49. pmid:26995111

* View Article

* PubMed/NCBI

* Google Scholar

10. 10. Datto JP, Shah AK, Bastidas JC, Arheart KL, Marcillo AE, Dietrich WD, et al. Use of the CatWalk gait device to assess differences in locomotion between genders in rats inherently and following spinal cord injury. Dataset Papers in Science. 2016;2016.

* View Article

* Google Scholar

11. 11. Tung TH, Mackinnon SE. Stem cell-based approaches to enhance nerve regeneration and improve functional outcomes in vascularized composite allotransplantation. Current opinion in organ transplantation. 2018;23(5):577–81. pmid:30138147

* View Article

* PubMed/NCBI

* Google Scholar

12. 12. Heinzel JC, Oberhauser V, Keibl C, Swiadek N, Langle G, Frick H, et al. Evaluation of Functional Recovery in Rats After Median Nerve Resection and Autograft Repair Using Computerized Gait Analysis. Front Neurosci. 2020;14:593545. Epub 2021/02/09. pmid:33551723; PubMed Central PMCID: PMC7859340.

* View Article

* PubMed/NCBI

* Google Scholar

13. 13. Zhang BG, Quigley AF, Myers DE, Wallace GG, Kapsa RM, Choong PF. Recent advances in nerve tissue engineering. The International Journal of Artificial Organs. 2014;37(4):277–91. pmid:24811182

* View Article

* PubMed/NCBI

* Google Scholar

14. 14. Timotius IK, Canneva F, Minakaki G, Moceri S, Plank AC, Casadei N, et al. Systematic data analysis and data mining in CatWalk gait analysis by heat mapping exemplified in rodent models for neurodegenerative diseases. J Neurosci Methods. 2019;326:108367. Epub 2019/07/28. pmid:31351096.

* View Article

* PubMed/NCBI

* Google Scholar

15. 15. Heinzel J, Langle G, Oberhauser V, Hausner T, Kolbenschlag J, Prahm C, et al. Use of the CatWalk gait analysis system to assess functional recovery in rodent models of peripheral nerve injury—a systematic review. J Neurosci Methods. 2020;345:108889. Epub 2020/08/07. pmid:32755615.

* View Article

* PubMed/NCBI

* Google Scholar

16. 16. Matias Junior I, Medeiros P, de Freita RL, Vicente-Cesar H, Ferreira Junior JR, Machado HR, et al. Effective Parameters for Gait Analysis in Experimental Models for Evaluating Peripheral Nerve Injuries in Rats. Neurospine. 2019;16(2):305–16. Epub 2019/01/18. pmid:30653907; PubMed Central PMCID: PMC6603843.

* View Article

* PubMed/NCBI

* Google Scholar

17. 17. Costa LM, Simoes MJ, Mauricio AC, Varejao AS. Chapter 7: Methods and protocols in peripheral nerve regeneration experimental research: part IV-kinematic gait analysis to quantify peripheral nerve regeneration in the rat. Int Rev Neurobiol. 2009;87:127–39. Epub 2009/08/18. pmid:19682636.

* View Article

* PubMed/NCBI

* Google Scholar

18. 18. Lu Y, Lin Z, Li M, Zhuang Y, Nie B, Lei J, et al. Three-phase Enriched Environment Improves Post-stroke Gait Dysfunction via Facilitating Neuronal Plasticity in the Bilateral Sensorimotor Cortex: A Multimodal MRI/PET Analysis in Rats. Neurosci Bull. 2023. Epub 2023/12/06. pmid:38055107.

* View Article

* PubMed/NCBI

* Google Scholar

19. 19. Dorman CW, Krug HE, Frizelle SP, Funkenbusch S, Mahowald ML. A comparison of DigiGait™ and TreadScan™ imaging systems: assessment of pain using gait analysis in murine monoarthritis. Journal of pain research. 2014;7:25.

* View Article

* Google Scholar

20. 20. Gabriel A, Marcus M, Honig W, Walenkamp G, Joosten E. The CatWalk method: a detailed analysis of behavioral changes after acute inflammatory pain in the rat. Journal of neuroscience methods. 2007;163(1):9–16. pmid:17382401

* View Article

* PubMed/NCBI

* Google Scholar

21. 21. Umansky D, Hagen KM, Chu TH, Pathiyil RK, Alzahrani S, Ousman SS, et al. Functional Gait Assessment Using Manual, Semi-Automated and Deep Learning Approaches Following Standardized Models of Peripheral Nerve Injury in Mice. Biomolecules. 2022;12(10). Epub 2022/10/28. pmid:36291564; PubMed Central PMCID: PMC9599622.

* View Article

* PubMed/NCBI

* Google Scholar

22. 22. Muller KA. Characterization & treatment of large sensory fiber peripheral neuropathy in diabetic mice: University of Kansas; 2008.

23. 23. Zu T, Guo S, Bardhi O, Ryskamp DA, Li J, Tusi SK, et al. Metformin inhibits RAN translation through PKR pathway and mitigates disease in C9orf72 ALS/FTD mice. Proceedings of the National Academy of Sciences. 2020;117(31):18591–9. pmid:32690681

* View Article

* PubMed/NCBI

* Google Scholar

24. 24. Zheng F, Tully A, Koss KM, Zhang X, Qiu L, Wang JJ, et al. Taking the Next Step: a Neural Coaptation Orthotopic Hind Limb Transplant Model to Maximize Functional Recovery in Rat. J Vis Exp. 2020;(162). Epub 2020/09/15. pmid:32925888.

* View Article

* PubMed/NCBI

* Google Scholar

25. 25. Petkova SP, Adhikari A, Berg EL, Fenton TA, Duis J, Silverman JL. Gait as a quantitative translational outcome measure in Angelman syndrome. Autism Res. 2022;15(5):821–33. Epub 2022/03/12. pmid:35274462; PubMed Central PMCID: PMC9311146.

* View Article

* PubMed/NCBI

* Google Scholar

26. 26. Ganguly A, McEwen C, Troy EL, Colburn RW, Caggiano AO, Schallert TJ, et al. Recovery of sensorimotor function following sciatic nerve injury across multiple rat strains. J Neurosci Methods. 2017;275:25–32. Epub 2016/12/17. pmid:27984099.

* View Article

* PubMed/NCBI

* Google Scholar

27. 27. Varejao AS, Meek MF, Ferreira AJ, Patricio JA, Cabrita AM. Functional evaluation of peripheral nerve regeneration in the rat: walking track analysis. J Neurosci Methods. 2001;108(1):1–9. Epub 2001/07/19. pmid:11459612.

* View Article

* PubMed/NCBI

* Google Scholar

28. 28. Furtmuller GJ, Oh B, Grahammer J, Lin CH, Sucher R, Fryer ML, et al. Orthotopic Hind Limb Transplantation in the Mouse. J Vis Exp. 2016;(108):53483. Epub 2016/03/12. pmid:26967527; PubMed Central PMCID: PMC4828154.

* View Article

* PubMed/NCBI

* Google Scholar

29. 29. Remeseiro B, Bolon-Canedo V. A review of feature selection methods in medical applications. Comput Biol Med. 2019;112:103375. Epub 2019/08/06. pmid:31382212.

* View Article

* PubMed/NCBI

* Google Scholar

30. 30. Rummel RJ. Applied factor analysis: Northwestern University Press; 1988.

* View Article

* Google Scholar

31. 31. Lambert CS, Philpot RM, Engberg ME, Johns BE, Kim SH, Wecker L. Gait analysis and the cumulative gait index (CGI): Translational tools to assess impairments exhibited by rats with olivocerebellar ataxia. Behav Brain Res. 2014;274:334–43. Epub 2014/08/15. pmid:25116252; PubMed Central PMCID: PMC4179979.

* View Article

* PubMed/NCBI

* Google Scholar

32. 32. Bommert A, Welchowski T, Schmid M, Rahnenfuhrer J. Benchmark of filter methods for feature selection in high-dimensional gene expression survival data. Brief Bioinform. 2022;23(1). Epub 2021/09/10. pmid:34498681; PubMed Central PMCID: PMC8769710.

* View Article

* PubMed/NCBI

* Google Scholar

33. 33. Kohavi R, John GH. Wrappers for feature subset selection. Artificial intelligence. 1997;97(1–2):273–324.

* View Article

* Google Scholar

34. 34. Pudjihartono N, Fadason T, Kempa-Liehr AW, O’Sullivan JM. A Review of Feature Selection Methods for Machine Learning-Based Disease Risk Prediction. Front Bioinform. 2022;2:927312. Epub 2022/10/29. pmid:36304293; PubMed Central PMCID: PMC9580915.

* View Article

* PubMed/NCBI

* Google Scholar

35. 35. Pirker W, Katzenschlager R. Gait disorders in adults and the elderly: A clinical guide. Wien Klin Wochenschr. 2017;129(3–4):81–95. Epub 2016/10/23. pmid:27770207; PubMed Central PMCID: PMC5318488.

* View Article

* PubMed/NCBI

* Google Scholar

36. 36. Cicirelli G, Impedovo D, Dentamaro V, Marani R, Pirlo G, D’Orazio TR. Human Gait Analysis in Neurodegenerative Diseases: A Review. IEEE J Biomed Health Inform. 2022;26(1):229–42. Epub 2021/06/29. pmid:34181559.

* View Article

* PubMed/NCBI

* Google Scholar

37. 37. Moore LK, Lee CS, Agha O, Liu M, Zhang H, Dang ABC, et al. A novel mouse model of hindlimb joint contracture with 3D-printed casts. J Orthop Res. 2022;40(12):2865–72. Epub 2022/03/11. pmid:35266583; PubMed Central PMCID: PMC10289010.

* View Article

* PubMed/NCBI

* Google Scholar

Citation: Naved BA, Han S, Koss KM, Kando MJ, Wang J-J, Weiss C, et al. (2025) Multivariate description of gait changes in a mouse model of peripheral nerve injury and trauma. PLoS ONE 20(1): e0312415. https://doi.org/10.1371/journal.pone.0312415

About the Authors:

Bilal A. Naved

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliations: Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, United States of America, Comprehensive Transplant Center, Department of Surgery, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

Shuling Han

Roles: Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Comprehensive Transplant Center, Department of Surgery, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

Kyle M. Koss

Roles: Conceptualization, Investigation, Project administration, Validation

Affiliations: Comprehensive Transplant Center, Department of Surgery, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America, Department of Surgery, College of Medicine, University of Arizona, Tucson, Arizona, United States of America

Mary J. Kando

Roles: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Behavioral Phenotyping Core, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

Jiao-Jing Wang

Roles: Conceptualization, Investigation

Affiliation: Comprehensive Transplant Center, Department of Surgery, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

ORICD: https://orcid.org/0000-0001-7865-0634

Craig Weiss

Roles: Conceptualization, Investigation, Methodology, Project administration, Resources, Writing – review & editing

Affiliation: Behavioral Phenotyping Core, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

ORICD: https://orcid.org/0000-0002-9861-8727

Maya G. Passman

Roles: Data curation, Formal analysis, Writing – original draft, Writing – review & editing

Affiliation: Barnard College, Columbia University, New York, NY United States of America

ORICD: https://orcid.org/0000-0002-1235-8427

Jason A. Wertheim

Roles: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

Affiliation: Department of Surgery, College of Medicine, University of Arizona, Tucson, Arizona, United States of America

Yuan Luo

Roles: Conceptualization, Investigation, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

E-mail: [email protected] (YL); [email protected] (ZJZ)

Affiliation: Division of Health and Biomedical Informatics, Department of Preventive Medicine, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

Zheng J. Zhang

Roles: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

E-mail: [email protected] (YL); [email protected] (ZJZ)

Affiliation: Comprehensive Transplant Center, Department of Surgery, Northwestern University Feinberg School of Medicine, Chicago, IL, United States of America

ORICD: https://orcid.org/0000-0001-6469-5778

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Chau T. A review of analytical techniques for gait data. Part 1: fuzzy, statistical and fractal methods. Gait & posture. 2001;13(1):49–66. pmid:11166554

2. Gorantla VS, Zor F, Nasir S, Breidenbach WC, Davis MR. Lower Extremity Transplantation: Concepts, Challenges, and Controversies. In: Tepe V, Peterson CM, editors. Full Stride: Advancing the State of the Art in Lower Extremity Gait Systems. New York, NY: Springer New York; 2017. p. 195–212.

3. Peripheral Neuropathy. In: Stroke NIoNDa, editor. NIH.gov: National Institutes of Health; 2024.

4. Geuna S, Raimondo S, Ronchi G, Di Scipio F, Tos P, Czaja K, et al. Chapter 3: Histology of the peripheral nerve and changes occurring during nerve regeneration. Int Rev Neurobiol. 2009;87:27–46. Epub 2009/08/18. pmid:19682632.

5. Hodes R, Larrabee MG, German W. The human electromyogram in response to nerve stimulation and the conduction velocity of motor axons; studies on normal and on injured peripheral nerves. Arch Neurol Psychiatry. 1948;60(4):340–65. Epub 1948/10/01. pmid:18111214.

6. Quan D, Bird SJ. Nerve conduction studies and electromyography in the evaluation of peripheral nerve injuries. Univ Pa Orthop J. 1999;12:45–51.

7. Chato-Astrain J, Philips C, Campos F, Durand-Herrera D, Garcia-Garcia OD, Roosens A, et al. Detergent-based decellularized peripheral nerve allografts: An in vivo preclinical study in the rat sciatic nerve injury model. J Tissue Eng Regen Med. 2020;14(6):789–806. Epub 2020/04/16. pmid:32293801.

8. Navarro X. Functional evaluation of peripheral nerve regeneration and target reinnervation in animal models: a critical overview. Eur J Neurosci. 2016;43(3):271–86. Epub 2015/08/01. pmid:26228942.

9. Lakes EH, Allen KD. Gait analysis methods for rodent models of arthritic disorders: reviews and recommendations. Osteoarthritis and cartilage. 2016;24(11):1837–49. pmid:26995111

10. Datto JP, Shah AK, Bastidas JC, Arheart KL, Marcillo AE, Dietrich WD, et al. Use of the CatWalk gait device to assess differences in locomotion between genders in rats inherently and following spinal cord injury. Dataset Papers in Science. 2016;2016.

11. Tung TH, Mackinnon SE. Stem cell-based approaches to enhance nerve regeneration and improve functional outcomes in vascularized composite allotransplantation. Current opinion in organ transplantation. 2018;23(5):577–81. pmid:30138147

12. Heinzel JC, Oberhauser V, Keibl C, Swiadek N, Langle G, Frick H, et al. Evaluation of Functional Recovery in Rats After Median Nerve Resection and Autograft Repair Using Computerized Gait Analysis. Front Neurosci. 2020;14:593545. Epub 2021/02/09. pmid:33551723; PubMed Central PMCID: PMC7859340.

13. Zhang BG, Quigley AF, Myers DE, Wallace GG, Kapsa RM, Choong PF. Recent advances in nerve tissue engineering. The International Journal of Artificial Organs. 2014;37(4):277–91. pmid:24811182

14. Timotius IK, Canneva F, Minakaki G, Moceri S, Plank AC, Casadei N, et al. Systematic data analysis and data mining in CatWalk gait analysis by heat mapping exemplified in rodent models for neurodegenerative diseases. J Neurosci Methods. 2019;326:108367. Epub 2019/07/28. pmid:31351096.

15. Heinzel J, Langle G, Oberhauser V, Hausner T, Kolbenschlag J, Prahm C, et al. Use of the CatWalk gait analysis system to assess functional recovery in rodent models of peripheral nerve injury—a systematic review. J Neurosci Methods. 2020;345:108889. Epub 2020/08/07. pmid:32755615.

16. Matias Junior I, Medeiros P, de Freita RL, Vicente-Cesar H, Ferreira Junior JR, Machado HR, et al. Effective Parameters for Gait Analysis in Experimental Models for Evaluating Peripheral Nerve Injuries in Rats. Neurospine. 2019;16(2):305–16. Epub 2019/01/18. pmid:30653907; PubMed Central PMCID: PMC6603843.

17. Costa LM, Simoes MJ, Mauricio AC, Varejao AS. Chapter 7: Methods and protocols in peripheral nerve regeneration experimental research: part IV-kinematic gait analysis to quantify peripheral nerve regeneration in the rat. Int Rev Neurobiol. 2009;87:127–39. Epub 2009/08/18. pmid:19682636.

18. Lu Y, Lin Z, Li M, Zhuang Y, Nie B, Lei J, et al. Three-phase Enriched Environment Improves Post-stroke Gait Dysfunction via Facilitating Neuronal Plasticity in the Bilateral Sensorimotor Cortex: A Multimodal MRI/PET Analysis in Rats. Neurosci Bull. 2023. Epub 2023/12/06. pmid:38055107.

19. Dorman CW, Krug HE, Frizelle SP, Funkenbusch S, Mahowald ML. A comparison of DigiGait™ and TreadScan™ imaging systems: assessment of pain using gait analysis in murine monoarthritis. Journal of pain research. 2014;7:25.

20. Gabriel A, Marcus M, Honig W, Walenkamp G, Joosten E. The CatWalk method: a detailed analysis of behavioral changes after acute inflammatory pain in the rat. Journal of neuroscience methods. 2007;163(1):9–16. pmid:17382401

21. Umansky D, Hagen KM, Chu TH, Pathiyil RK, Alzahrani S, Ousman SS, et al. Functional Gait Assessment Using Manual, Semi-Automated and Deep Learning Approaches Following Standardized Models of Peripheral Nerve Injury in Mice. Biomolecules. 2022;12(10). Epub 2022/10/28. pmid:36291564; PubMed Central PMCID: PMC9599622.

22. Muller KA. Characterization & treatment of large sensory fiber peripheral neuropathy in diabetic mice: University of Kansas; 2008.

23. Zu T, Guo S, Bardhi O, Ryskamp DA, Li J, Tusi SK, et al. Metformin inhibits RAN translation through PKR pathway and mitigates disease in C9orf72 ALS/FTD mice. Proceedings of the National Academy of Sciences. 2020;117(31):18591–9. pmid:32690681

24. Zheng F, Tully A, Koss KM, Zhang X, Qiu L, Wang JJ, et al. Taking the Next Step: a Neural Coaptation Orthotopic Hind Limb Transplant Model to Maximize Functional Recovery in Rat. J Vis Exp. 2020;(162). Epub 2020/09/15. pmid:32925888.

25. Petkova SP, Adhikari A, Berg EL, Fenton TA, Duis J, Silverman JL. Gait as a quantitative translational outcome measure in Angelman syndrome. Autism Res. 2022;15(5):821–33. Epub 2022/03/12. pmid:35274462; PubMed Central PMCID: PMC9311146.

26. Ganguly A, McEwen C, Troy EL, Colburn RW, Caggiano AO, Schallert TJ, et al. Recovery of sensorimotor function following sciatic nerve injury across multiple rat strains. J Neurosci Methods. 2017;275:25–32. Epub 2016/12/17. pmid:27984099.

27. Varejao AS, Meek MF, Ferreira AJ, Patricio JA, Cabrita AM. Functional evaluation of peripheral nerve regeneration in the rat: walking track analysis. J Neurosci Methods. 2001;108(1):1–9. Epub 2001/07/19. pmid:11459612.

28. Furtmuller GJ, Oh B, Grahammer J, Lin CH, Sucher R, Fryer ML, et al. Orthotopic Hind Limb Transplantation in the Mouse. J Vis Exp. 2016;(108):53483. Epub 2016/03/12. pmid:26967527; PubMed Central PMCID: PMC4828154.

29. Remeseiro B, Bolon-Canedo V. A review of feature selection methods in medical applications. Comput Biol Med. 2019;112:103375. Epub 2019/08/06. pmid:31382212.

30. Rummel RJ. Applied factor analysis: Northwestern University Press; 1988.

31. Lambert CS, Philpot RM, Engberg ME, Johns BE, Kim SH, Wecker L. Gait analysis and the cumulative gait index (CGI): Translational tools to assess impairments exhibited by rats with olivocerebellar ataxia. Behav Brain Res. 2014;274:334–43. Epub 2014/08/15. pmid:25116252; PubMed Central PMCID: PMC4179979.

32. Bommert A, Welchowski T, Schmid M, Rahnenfuhrer J. Benchmark of filter methods for feature selection in high-dimensional gene expression survival data. Brief Bioinform. 2022;23(1). Epub 2021/09/10. pmid:34498681; PubMed Central PMCID: PMC8769710.

33. Kohavi R, John GH. Wrappers for feature subset selection. Artificial intelligence. 1997;97(1–2):273–324.

34. Pudjihartono N, Fadason T, Kempa-Liehr AW, O’Sullivan JM. A Review of Feature Selection Methods for Machine Learning-Based Disease Risk Prediction. Front Bioinform. 2022;2:927312. Epub 2022/10/29. pmid:36304293; PubMed Central PMCID: PMC9580915.

35. Pirker W, Katzenschlager R. Gait disorders in adults and the elderly: A clinical guide. Wien Klin Wochenschr. 2017;129(3–4):81–95. Epub 2016/10/23. pmid:27770207; PubMed Central PMCID: PMC5318488.

36. Cicirelli G, Impedovo D, Dentamaro V, Marani R, Pirlo G, D’Orazio TR. Human Gait Analysis in Neurodegenerative Diseases: A Review. IEEE J Biomed Health Inform. 2022;26(1):229–42. Epub 2021/06/29. pmid:34181559.

37. Moore LK, Lee CS, Agha O, Liu M, Zhang H, Dang ABC, et al. A novel mouse model of hindlimb joint contracture with 3D-printed casts. J Orthop Res. 2022;40(12):2865–72. Epub 2022/03/11. pmid:35266583; PubMed Central PMCID: PMC10289010.