Gait characterization in golden retriever muscular dystrophy dogs using linear discriminant analysis
- Bodvaël Fraysse†1Email author,
- Inès Barthélémy†2, 3, 4,
- El Mostafa Qannari5,
- Karl Rouger6,
- Chantal Thorin7,
- Stéphane Blot2, 3, 4,
- Caroline Le Guiner1,
- Yan Chérel^6 and
- Jean-Yves Hogrel8
© The Author(s). 2017
Received: 29 July 2016
Accepted: 20 March 2017
Published: 12 April 2017
Accelerometric analysis of gait abnormalities in golden retriever muscular dystrophy (GRMD) dogs is of limited sensitivity, and produces highly complex data. The use of discriminant analysis may enable simpler and more sensitive evaluation of treatment benefits in this important preclinical model.
Accelerometry was performed twice monthly between the ages of 2 and 12 months on 8 healthy and 20 GRMD dogs. Seven accelerometric parameters were analysed using linear discriminant analysis (LDA). Manipulation of the dependent and independent variables produced three distinct models. The ability of each model to detect gait alterations and their pattern change with age was tested using a leave-one-out cross-validation approach.
Selecting genotype (healthy or GRMD) as the dependent variable resulted in a model (Model 1) allowing a good discrimination between the gait phenotype of GRMD and healthy dogs. However, this model was not sufficiently representative of the disease progression. In Model 2, age in months was added as a supplementary dependent variable (GRMD_2 to GRMD_12 and Healthy_2 to Healthy_9.5), resulting in a high overall misclassification rate (83.2%). To improve accuracy, a third model (Model 3) was created in which age was also included as an explanatory variable. This resulted in an overall misclassification rate lower than 12%. Model 3 was evaluated using blinded data pertaining to 81 healthy and GRMD dogs. In all but one case, the model correctly matched gait phenotype to the actual genotype. Finally, we used Model 3 to reanalyse data from a previous study regarding the effects of immunosuppressive treatments on muscular dystrophy in GRMD dogs. Our model identified significant effect of immunosuppressive treatments on gait quality, corroborating the original findings, with the added advantages of direct statistical analysis with greater sensitivity and more comprehensible data representation.
Gait analysis using LDA allows for improved analysis of accelerometry data by applying a decision-making analysis approach to the evaluation of preclinical treatment benefits in GRMD dogs.
KeywordsMuscular dystrophy GRMD Treatment evaluation Gait assessment Discriminant analysis Animal model Accelerometry
Duchenne muscular dystrophy (DMD) is an X-linked disorder caused by various mutations in the gene encoding for dystrophin, resulting in the absence of the functional protein in muscle fibres . DMD patients display progressive muscle weakness leading to the loss of independent mobility in young adolescents and respiratory and heart failure in young adults. While gene, cell, and pharmacological therapies have all been investigated [2, 3], there is currently no curative therapy for DMD. However, several issues such as efficacy of studied drugs, locoregional or systemic medication pathways or dosage, remain to be further explored. The preclinical DMD model of choice is the dystrophin-deficient golden retriever muscular dystrophy (GRMD) dog, which closely mimics many aspects of the human disease [4, 5]. Scale is a key factor that influences the translation of data from animal models to humans. Thus, when studying mechanical impacts, molecular diffusion and/or cell migration data acquired in GRMD dogs is much more relevant to human DMD than that obtained in smaller animal models such as the mdx mutant mouse . The downside of the GRMD model is that it is more expensive to purchase and house, and the use of "man's best friend" for research purposes entails additional political and ethical considerations. Accordingly, limited numbers of GRMD dogs are generally used in studies. At first glance, this appears to be a major impediment to the preclinical evaluation of therapies, particularly given that GRMD dogs exhibit considerable inter-individual phenotypic variability [6–8]. However, there are numerous similarities between canine and human diseases, and considerable inter-individual variability is also observed among human patients . Moreover, given the rarity of muscular dystrophies, and for obvious ethical reasons, clinical studies in man generally involve limited numbers of patients and pose similar challenges to those performed in GRMD dogs. Tools that allow the establishment of better readouts of disease progression and treatment response are essential to overcome limitations imposed by small sample sizes and wide inter-individual variability in studies using the GRMD model, and to better predict the pathogenesis of muscular dystrophy and treatment efficacy in humans.
Studies by several groups have conducted gait analysis in GRMD dogs [10–14]. Using accelerometry analysis in these animals we have demonstrated less regular and less powerful acceleration, decreased stride length and frequency, and a redistribution of power from the cranio-caudal to the medio-lateral axis . Using the main gait variables, we developed a global gait index that consistently detected early changes in gait patterns in GRMD, as well as the progressive deterioration of gait quality. This index was based on the use of principal component analysis (PCA) and the computation of Euclidean distances at multiple time points with respect to an age-matched control group . Given the inherent complexity of the methodology, we believe that this approach is not best suited to the problem at hand, and in fact may hinder the evaluation of therapies in preclinical studies in GRMD dogs. To increase sensitivity and aid interpretation of the outcomes, we sought to design a simpler and more appropriate analytical method using linear discriminant analysis (LDA). Like PCA, LDA is an orthogonal transformation and data reduction technique, but unlike PCA, LDA seeks to minimize intra-group variance and maximize inter-group variance. Moreover, LDA yields a predictive model based on control group data. This is achieved by computation of group membership based on experimental data and assignment rules, which allow the prediction of group membership for future observations. Here, we describe a new method of 3D accelerometric gait analysis using LDA. We discuss the choice of dependent and independent variables and describe how to represent the results in a manner that can be understood by a broad range of users, including those unfamiliar with LDA. Finally, we demonstrate the validity of this method by reanalysing data from a previous study regarding immunosuppressive treatment in GRMD dogs.
All procedures were carried out in accordance with Guidelines for the Care and Use of Laboratory Animals, and approved by the common Ethical Committee of the National Veterinary School of Alfort, ANSES, and UPEC. Eight healthy golden retrievers and 20 GRMD dogs (all males) were included in the study. All animals came from the French GRMD colony. The healthy dogs were littermates of some of the GRMD dogs used. Six of the healthy dogs and 12 of the GRMD dogs had participated in a previous study . All dogs were housed in the same facilities, and were genotyped as previously described . Only the GRMD dogs that were still ambulatory after 9 months of age were included in this study.
Evaluation of gait quality using 3D accelerometry
As previously described , the 3-dimensional accelerometer recorder used was a Locometrix® gait analysis system, composed of three orthogonally positioned accelerometers, which can record accelerations along the dorso-ventral, cranio-caudal, and medio-lateral axes.
To assess the efficiency of our method, we analysed gait data acquired for four immunosuppressed dogs that had been published in a previous paper . These dogs had been treated with high doses of oral prednisolone (2 mg/kg/d) and cyclosporine A (initial dose of 20 mg/kg/d) between 2 and 9 months of age.
Dogs were carried from the kennel to a 45-meter-long testing corridor located close to the laboratory facilities, as previously described . The belt to which the accelerometric device was attached was fastened around the thorax of the dog, near the centre of gravity at rest. Each animal was tested twice per month from 2 months of age (when motor clinical signs appear in animals that survive the neonatal period) to 9.5 months of age, thus covering the period of growth and disease progression . To enhance the discriminatory power of the method, values obtained for GRMD dogs that survived to 12 months were also included in the analysis. Young puppies were familiarized with the corridor and the belt before the test. The height at withers (HW) was measured at the end of each test. All tests were performed by the same experimenter (IB). In each test, the dog was encouraged to walk or run at its preferred gait, and its speed calculated over a distance of five meters, as previously described .
For quadrupedal gait analysis, acceleration curves were analysed using the software provided by the manufacturer of the recording device (Equimetrix®, Centaure Metrix, Evry, France). A 10-second sequence of steady-state locomotion, which was easily identifiable in the dorso-ventral acceleration curves, was analysed. The following variables, which have been previously described in detail [10, 11], were computed: stride frequency (SF, /s), stride regularity (Reg, dimensionless), total power of accelerations (TP, W/kg), relative components of the total power along the three axes (%) (calculated by dividing cranio-caudal (CCP), dorso-ventral (DVP), or medio-lateral power (MLP) by total power (TP)), and stride length (calculated by dividing the speed by SF), which was normalized to height at withers (SL/HW) in order to circumvent the effect of limb length on this variable. For the sake of consistency, only observations with a Reg value >70 were considered for analysis. Gait testing consisted of two consecutive round trips in the corridor. However, in contrast to healthy dogs, it was very difficult, if not impossible, for some GRMD dogs to complete the second round trip.
The main objective of the present study was to assess the capacity of discriminant analysis (DA) to evaluate gait, detect functional alterations, and evaluate treatment benefits during the growth period in GRMD dogs. We used gait data obtained from healthy and GRMD dogs and performed DA using XLSTAT software (Addinsoft™). This Excel add-in extends the analytical functions of Excel and covers the key requirements for data analysis and statistics. DA is a commonly used multivariate data analysis method. The aim of this supervised method is to predict group memberships of a set of individuals based on multivariate data. In this scenario, the groups of individuals are assumed to be known a priori. DA reduces the dimensionality of the data at hand by computing synthetic variables, often called canonical variables or factors, the aim of which is to maximize inter-group variance while minimizing intra-group variance. DA yields new variables, which are linear combinations of the original variables. The maximum number of such variables is equal to the number of groups minus one and are usually noted F1 to F(n-1) where n is the number of groups. However, in practice, only the first few canonical variables are used for the purpose of discrimination, since the remaining canonical variables may be predominantly associated with noise present in the data. The graphical displays generated using the retained canonical variables are useful to depict the separation of groups. DA methods include linear DA (LDA) and quadratic DA (QDA) . LDA is a parametric method that assumes Gaussian distributions with the same variance-covariance matrix within the various groups. In practice however, this is often not the case. Nonetheless, even in the case of a slight deviation from this requirement, LDA performs reasonably well , and was thus the method selected for the present study.
Blind test for model validation
We previously demonstrated that PCA, although a non-supervised analysis strategy, can distinguish the gait phenotype of GRMD dogs from that of healthy controls with reasonable accuracy, and does so independently of age . We show in the present study that similar results could be obtained using LDA. However, when the age of the dogs were not take into account, the model built using LDA (Model 1) was found not to adequately represent the progression of the disease, which evolves over the course of the postnatal growth phase. On the other hand, addition of age, in months, as a supplementary dependent variable in LDA, the model obtained (Model 2) exhibited a very high rate of misclassification. This likely reflected the slow rate of disease progression over the short age intervals analysed.
Nonetheless, in Model2, the discriminant axis mainly associated with age and disease progression, F2, presented a low canonical correlation, indicating that the model was not sufficiently sensitive to accurately evaluate the impact of progressive dystrophy on gait in growing GRMD dogs. One means of improving the model was to increase the age interval (e.g. from 0.5 to 1 month), which would certainly reduce the misclassification error, but would also decrease the sensitivity of the model. The aim of this study was not to demonstrate altered gait in GRMD versus healthy dogs, but rather to develop a method to better identify gait alterations and their progression with age in GRMD dogs. To improve the descriptive properties of the model while keeping its predictive capabilities, we have artificially stretched the model (Model 3) in the direction of the progression of age, by introducing age (in days) as an explanatory (or dependent) variable in the LDA analysis, in addition to including age (in months) as a dependent variable. The dogs included in the study were not born on the same day. Additionally, sometimes some dogs were not able to walk, thus their acquisition session was delayed. By contrast, acquisition sessions were always done on the same day of the week, every 15 days. Thus, the age of the animals in days in the same age group in months could vary from about 7 days. The intrinsic variability of the age in days spread the model in the factorial plan helping to discriminate the different groups of age, although it could be assumed that this variability does not significantly impact the measurements. In Model 3, the two first canonical discriminant factors, F1 and F2, accounted for 99.7% of variance. F1 was predominantly correlated with age whereas F2 was discriminant for gait phenotype. Model 3 presented a low overall misclassification error and high canonical correlations of 0.99 and of 0.95 were calculated for F1 and F2, respectively.
In order to investigate further the suitability of Model 3 as a tool for gait analysis in preclinical studies with GRMD dogs, we used it to predict the genotypes of a cohort of 81 healthy and GRMD dogs starting from their gait accelerometry characteristics. For each dog of the cohort, a number was randomly attributed by an external scientist blinded to the data. Using this approach, when “the data were unblinded”, the predicted gait phenotype matched the actual phenotype and genotype in all cases but one. These findings strongly support the robustness and accuracy of the method.
Some members of our group tested, in a previous study, the effects of immunosuppressive treatments (oral administration of cyclosporine A and corticosteroids) on muscular dystrophy and overall health in GRMD dogs . Indeed, immunomodulatory treatments have been employed in several studies assessing the effectiveness of gene, cell, or pharmacological therapies in dog models of DMD to supress the immune response to the viral vector, donor cells, and/or the transgene product [18–20]. The obtained results vary considerably depending on the variables measured. Although Barthélémy and colleagues  reported a more severe disease progression in terms of isometric force and histology, they also found an improvement in gait classification using a PCA-based gait index calculated from the seven accelerometric variables used in the present study. This longitudinal analysis was complex, as PCA, unlike DA, is an unsupervised approach. The authors performed PCA for each age category and calculated the corresponding Euclidean distance for each GRMD dog to the centroid of age-matched healthy dogs, and, thereafter, plotted the evolution of this distance with age. Using Model 3, we found that immunosuppression in GRMD dogs had a significant beneficial effect on gait, although this was both partial and temporary. This finding was in line with those of Barthélémy and coworkers . Furthermore, LDA revealed that significant effects of treatment on gait observed at 4 months of age were preceded by signs of improvement that were evident as early as 3 months, supporting the sensitivity of our model.
Using LDA we have generated a highly sensitive model of gait alterations due to muscular dystrophy in the GRMD dog. Our model shows a high degree of discriminatory accuracy, distinguishing the gait phenotype of GRMD dogs from that of healthy dogs as early as 2.5 months of age, and thus overcoming some of the difficulties in analysing a progressive disease that occurs during the growth phase of postnatal development. Moreover, we designed a new means of representing the outcomes of our analysis that allows for easier interpretation of the results. This is a key strength of our study: because preclinical results are the base upon which phase 1/2 clinical trials in human patients are prepared and designed, preclinical outcomes should be completely understandable to all those involved in the design and the testing of potential treatments from preclinical through to clinical phases. In our experience, graphical displays depicting only the two variables (F1 and F2) are very likely to be unclear, and sometimes misleading, to practitioners unfamiliar with DA or PCA.
We previously demonstrated that accelerometry combined with PCA constitutes a reliable follow-up tool for gait analysis in preclinical therapeutic trials using GRMD dogs. However, this approach has some limitations, including limited sensitivity, and high complexity of the data generated. The model presented here greatly improves upon this method by employing LDA.
Duchenne muscular dystrophy
Golden retriever muscular dystrophy
- HW :
Height at the Withers
Linear discriminant analysis
Principal component analysis
Quadratic discriminant analysis
- Reg :
- SF :
- TP :
We thank the Centre d’Elevage du Domaine des Souches to have bred the GRMD dogs for this study, and Dr Pablo Aguilar, Xavier Cauchois and the whole team of the ENVA neurobiology research unit, for their very good daily care to the dogs.
This work was supported by the Association Française contre les Myopathies (AFM), and by the Ministère de l’Agriculture, de l’Agroalimentaire et de la Forêt. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and material
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.
BF designed the experiments, analysed the data and drafted the manuscript. IB designed and performed the experiments, analysed the data and drafted the manuscript. EMQ analysed the data and drafted the manuscript. KR designed the experiments and drafted the manuscript. CT analysed the data and drafted the manuscript. SB, CLG, YC and JYH designed the experiments, analysed the data and drafted the manuscript. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
The present study does not involve human participants, human data or human tissue. Concerning animal data, all procedures were carried out in accordance with Guidelines for the Care and Use of Laboratory Animals, and approved by the common Ethical Committee of the National Veterinary School of Alfort, ANSES, and UPEC.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Koenig M, Hoffman EP, Bertelson CJ, Monaco AP, Feener C, Kunkel LM. Complete cloning of the Duchenne muscular dystrophy (DMD) cDNA and preliminary genomic organization of the DMD gene in normal and affected individuals. Cell. 1987;50(3):509–17.View ArticlePubMedGoogle Scholar
- Wein N, Alfano L, Flanigan KM. Genetics and emerging treatments for Duchenne and Becker muscular dystrophy. Pediatr Clin North Am. 2015;62(3):723–42.View ArticlePubMedGoogle Scholar
- Sienkiewicz D, Kulak W, Okurowska-Zawada B, Paszko-Patej G, Kawnik K. Duchenne muscular dystrophy: current cell therapies. Ther Adv Neurol Disord. 2015;8(4):166–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Valentine BA, Winand NJ, Pradhan D, Moise NS, de Lahunta A, Kornegay JN, Cooper BJ. Canine X-linked muscular dystrophy as an animal model of Duchenne muscular dystrophy: a review. Am J Med Genet. 1992;42(3):352–6.View ArticlePubMedGoogle Scholar
- McGreevy JW, Hakim CH, McIntosh MA, Duan D. Animal models of Duchenne muscular dystrophy: from basic mechanisms to gene therapy. Dis Model Mech. 2015;8(3):195–213.View ArticlePubMedPubMed CentralGoogle Scholar
- Kornegay JN, Bogan JR, Bogan DJ, Childers MK, Li J, Nghiem P, Detwiler DA, Larsen CA, Grange RW, Bhavaraju-Sanka RK, et al. Canine models of Duchenne muscular dystrophy and their use in therapeutic strategies. Mamm Genome. 2012;23(1–2):85–108.View ArticlePubMedPubMed CentralGoogle Scholar
- Willmann R, Possekel S, Dubach-Powell J, Meier T, Ruegg MA. Mammalian animal models for Duchenne muscular dystrophy. Neuromuscul Disord. 2009;19(4):241–9.View ArticlePubMedGoogle Scholar
- Vieira NM, Guo LT, Estrela E, Kunkel LM, Zatz M, Shelton GD. Muscular dystrophy in a family of Labrador Retrievers with no muscle dystrophin and a mild phenotype. Neuromuscul Disord. 2015;25(5):363–70.View ArticlePubMedGoogle Scholar
- Zatz M, Pavanello RC, Lazar M, Yamamoto GL, Lourenco NC, Cerqueira A, Nogueira L, Vainzof M. Milder course in Duchenne patients with nonsense mutations and no muscle dystrophin. Neuromuscul Disord. 2014;24(11):986–9.View ArticlePubMedGoogle Scholar
- Barthelemy I, Barrey E, Aguilar P, Uriarte A, Le Chevoir M, Thibaud JL, Voit T, Blot S, Hogrel JY. Longitudinal ambulatory measurements of gait abnormality in dystrophin-deficient dogs. BMC Musculoskelet Disord. 2011;12:75.View ArticlePubMedPubMed CentralGoogle Scholar
- Barthelemy I, Barrey E, Thibaud JL, Uriarte A, Voit T, Blot S, Hogrel JY. Gait analysis using accelerometry in dystrophin-deficient dogs. Neuromuscul Disord. 2009;19(11):788–96.View ArticlePubMedGoogle Scholar
- Barthelemy I, Uriarte A, Drougard C, Unterfinger Y, Thibaud JL, Blot S. Effects of an immunosuppressive treatment in the GRMD dog model of Duchenne muscular dystrophy. PLoS One. 2012;7(11), e48478.View ArticlePubMedPubMed CentralGoogle Scholar
- Marsh AP, Eggebeen JD, Kornegay JN, Markert CD, Childers MK. Kinematics of gait in golden retriever muscular dystrophy. Neuromuscul Disord. 2010;20(1):16–20.View ArticlePubMedGoogle Scholar
- Shin JH, Greer B, Hakim CH, Zhou Z, Chung YC, Duan Y, He Z, Duan D. Quantitative phenotyping of Duchenne muscular dystrophy dogs by comprehensive gait analysis and overnight activity monitoring. PLoS One. 2013;8(3), e59875.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartlett RJ, Winand NJ, Secore SL, Singer JT, Fletcher S, Wilton S, Bogan DJ, Metcalf-Bogan JR, Bartlett WT, Howell JM, et al. Mutation segregation and rapid carrier detection of X-linked muscular dystrophy in dogs. Am J Vet Res. 1996;57(5):650–4.PubMedGoogle Scholar
- Stewart S, Ivy MA, Anslyn EV. The use of principal component analysis and discriminant analysis in differential sensing routines. Chem Soc Rev. 2014;43(1):70–84.View ArticlePubMedGoogle Scholar
- Nathan R, Spiegel O, Fortmann-Roe S, Harel R, Wikelski M, Getz WM. Using tri-axial acceleration data to identify behavioral modes of free-ranging animals: general concepts and tools illustrated for griffon vultures. J Exp Biol. 2012;215(Pt 6):986–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Dell'Agnola C, Wang Z, Storb R, Tapscott SJ, Kuhr CS, Hauschka SD, Lee RS, Sale GE, Zellmer E, Gisburne S, et al. Hematopoietic stem cell transplantation does not restore dystrophin expression in Duchenne muscular dystrophy dogs. Blood. 2004;104(13):4311–8.View ArticlePubMedGoogle Scholar
- Sampaolesi M, Blot S, D'Antona G, Granger N, Tonlorenzi R, Innocenzi A, Mognol P, Thibaud JL, Galvez BG, Barthelemy I, et al. Mesoangioblast stem cells ameliorate muscle function in dystrophic dogs. Nature. 2006;444(7119):574–9.View ArticlePubMedGoogle Scholar
- Rouger K, Larcher T, Dubreil L, Deschamps JY, Le Guiner C, Jouvion G, Delorme B, Lieubeau B, Carlus M, Fornasari B, et al. Systemic delivery of allogenic muscle stem cells induces long-term muscle repair and clinical efficacy in duchenne muscular dystrophy dogs. Am J Pathol. 2011;179(5):2501–18.View ArticlePubMedPubMed CentralGoogle Scholar