BackgroundMolecular components in blood, such as proteins, are used as biomarkers to detect or predict disease states, guide clinical interventions and aid in the development of therapies. While... Show moreBackgroundMolecular components in blood, such as proteins, are used as biomarkers to detect or predict disease states, guide clinical interventions and aid in the development of therapies. While multiplexing proteomics methods promote discovery of such biomarkers, their translation to clinical use is difficult due to the lack of substantial evidence regarding their reliability as quantifiable indicators of disease state or outcome. To overcome this challenge, a novel orthogonal strategy was developed and used to assess the reliability of biomarkers and analytically corroborate already identified serum biomarkers for Duchenne muscular dystrophy (DMD). DMD is a monogenic incurable disease characterized by progressive muscle damage that currently lacks reliable and specific disease monitoring tools.MethodsTwo technological platforms are used to detect and quantify the biomarkers in 72 longitudinally collected serum samples from DMD patients at 3 to 5 timepoints. Quantification of the biomarkers is achieved by detection of the same biomarker fragment either through interaction with validated antibodies in immuno-assays or through quantification of peptides by Parallel Reaction Monitoring Mass Spectrometry assay (PRM-MS).ResultsFive, out of ten biomarkers previously identified by affinity-based proteomics methods, were confirmed to be associated with DMD using the mass spectrometry-based method. Two biomarkers, carbonic anhydrase III and lactate dehydrogenase B, were quantified with two independent methods, sandwich immunoassays and PRM-MS, with Pearson correlations of 0.92 and 0.946 respectively. The median concentrations of CA3 and LDHB in DMD patients was elevated in comparison to those in healthy individuals by 35- and 3-fold, respectively. Levels of CA3 vary between 10.26 and 0.36 ng/ml in DMD patients whereas those of LDHB vary between 15.1 and 0.8 ng/ml.ConclusionsThese results demonstrate that orthogonal assays can be used to assess the analytical reliability of biomarker quantification assays, providing a means to facilitate the translation of biomarkers to clinical practice. This strategy also warrants the development of the most relevant biomarkers, markers that can be reliably quantified with different proteomics methods. Show less
Heezen, L.G.M.; Abdelaal, T.; Putten, M. van; Aartsma-Rus, A.; Mahfouz, A.; Spitali, P. 2023
Duchenne muscular dystrophy is caused by mutations in the DMD gene, leading to lack of dystrophin. Chronic muscle damage eventually leads to histological alterations in skeletal muscles. The... Show moreDuchenne muscular dystrophy is caused by mutations in the DMD gene, leading to lack of dystrophin. Chronic muscle damage eventually leads to histological alterations in skeletal muscles. The identification of genes and cell types driving tissue remodeling is a key step to developing effective therapies. Here we use spatial transcriptomics in two Duchenne muscular dystrophy mouse models differing in disease severity to identify gene expression signatures underlying skeletal muscle pathology and to directly link gene expression to muscle histology. We perform deconvolution analysis to identify cell types contributing to histological alterations. We show increased expression of specific genes in areas of muscle regeneration (Myl4, Sparc, Hspg2), fibrosis (Vim, Fn1, Thbs4) and calcification (Bgn, Ctsk, Spp1). These findings are confirmed by smFISH. Finally, we use differentiation dynamic analysis in the D2-mdx muscle to identify muscle fibers in the present state that are predicted to become affected in the future state. Show less
BackgroundBecker muscular dystrophy (BMD) is an X-linked disorder characterized by slow, progressive muscle damage and muscle weakness. Hallmarks include fibre-size variation and replacement of... Show moreBackgroundBecker muscular dystrophy (BMD) is an X-linked disorder characterized by slow, progressive muscle damage and muscle weakness. Hallmarks include fibre-size variation and replacement of skeletal muscle with fibrous and adipose tissues, after repeated cycles of regeneration. Muscle histology can detect these features, but the required biopsies are invasive, are difficult to repeat and capture only small muscle volumes. Diffusion-tensor magnetic resonance imaging (DT-MRI) is a potential non-invasive alternative that can calculate muscle fibre diameters when applied with the novel random permeable barrier model (RPBM). In this study, we assessed muscle fibre diameters using DT-MRI in BMD patients and healthy controls and compared these with histology. MethodsWe included 13 BMD patients and 9 age-matched controls, who underwent water-fat MRI and DT-MRI at multiple diffusion times, allowing RPBM parameter estimation in the lower leg muscles. Tibialis anterior muscle biopsies were taken from the contralateral leg in 6 BMD patients who underwent DT-MRI and from an additional 32 BMD patients and 15 healthy controls. Laminin and Sirius-red stainings were performed to evaluate muscle fibre morphology and fibrosis. Twelve ambulant patients from the MRI cohort underwent the North Star ambulatory assessment, and 6-min walk, rise-from-floor and 10-m run/walk functional tests. ResultsRPBM fibre diameter was significantly larger in BMD patients (P = 0.015): mean (SD) = 68.0 (25.3) mu m versus 59.4 (19.2) mu m in controls. Inter-muscle differences were also observed (P <= 0.002). Both inter- and intra-individual RPBM fibre diameter variability were similar between groups. Laminin staining agreed with the RPBM, showing larger median fibre diameters in patients than in controls: 72.5 (7.9) versus 63.2 (6.9) mu m, P = 0.006. However, despite showing similar inter-individual variation, patients showed more intra-individual fibre diameter variability than controls-mean variance (SD) = 34.2 (7.9) versus 21.4 (6.9) mu m, P < 0.001-and larger fibrosis areas: median (interquartile range) = 21.7 (5.6)% versus 14.9 (3.4)%, P < 0.001. Despite good overall agreement of RPBM and laminin fibre diameters, they were not associated in patients who underwent DT-MRI and muscle biopsy, perhaps due to lack of colocalization of DT-MRI with biopsy samples. ConclusionsDT-MRI RPBM metrics agree with histology and can quantify changes in muscle fibre size that are associated with regeneration without the need for biopsies. They therefore show promise as imaging biomarkers for muscular dystrophies. Show less
Signorelli, M.; Tsonka, R.; Aartsma-Rus, A.; Spitali, P. 2023
Duchenne muscular dystrophy (DMD) is caused by genetic mutations leading to lack of dystrophin in skeletal muscle. A better understanding of how objective biomarkers for DMD vary across subjects... Show moreDuchenne muscular dystrophy (DMD) is caused by genetic mutations leading to lack of dystrophin in skeletal muscle. A better understanding of how objective biomarkers for DMD vary across subjects and over time is needed to model disease progression and response to therapy more effectively, both in pre-clinical and clinical research. We present an in-depth characterization of disease progression in 3 murine models of DMD by multiomic analysis of longitudinal trajectories between 6 and 30 weeks of age. Integration of RNA-seq, mass spectrometry-based metabolomic and lipidomic data obtained in muscle and blood samples by Multi-Omics Factor Analysis (MOFA) led to the identification of 8 latent factors that explained 78.8% of the variance in the multiomic dataset. Latent factors could discriminate dystrophic and healthy mice, as well as different time-points. MOFA enabled to connect the gene expression signature in dystrophic muscles, characterized by pro-fibrotic and energy metabolism alterations, to inflammation and lipid signatures in blood. Our results show that omic observations in blood can be directly related to skeletal muscle pathology in dystrophic muscle. Show less
Velde, N.M. van de; Koeks, Z.; Signorelli, M.; Verwey, N.; Overzier, M.; Bakker, J.A.; ... ; Niks, E.H. 2023
Background and Objectives: The slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3... Show moreBackground and Objectives: The slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3 muscle-enriched biomarkers in serum of patients with BMD over 4-year time and studied associations with disease severity, disease progression, and dystrophin levels in BMD. Methods: We quantitatively measured creatine kinase (CK) using the International Federation of Clinical Chemistry reference method, creatine/creatinine(ratio) (Cr/Crn) using liquid chromatography-tandem mass spectrometry, and myostatin with ELISA in serum and assessed functional performance using the North Star Ambulatory Assessment (NSAA), 10-meter run velocity (TMRv), 6-Minute Walking Test (6MWT), and forced vital capacity in a 4-year prospective natural history study. Dystrophin levels were quantified in the tibialis anterior muscle using capillary Western immunoassay. The correlation between biomarkers, age, functional performance, mean annual change, and prediction of concurrent functional performance was analyzed using linear mixed models. Results: Thirty-four patients with 106 visits were included. Eight patients were nonambulant at baseline. Cr/Crn and myostatin were highly patient specific (intraclass correlation coefficient for both = 0.960). Cr/Crn was strongly negatively correlated, whereas myostatin was strongly positively correlated with the NSAA, TMRv, and 6MWT (Cr/Crn rho = -0.869 to -0.801 and myostatin rho = 0.792 to 0.842, all p < 0.001). CK showed a negative association with age (p = 0.0002) but was not associated with patients' performance. Cr/Crn and myostatin correlated moderately with the average annual change of the 6MWT (rho = -0.532 and 0.555, p = 0.02). Dystrophin levels did not correlate with the selected biomarkers nor with performance. Cr/Crn, myostatin, and age could explain up to 75% of the variance of concurrent functional performance of the NSAA, TMRv, and 6MWT. Discussion: Both Cr/Crn and myostatin could potentially serve as monitoring biomarkers in BMD, as higher Cr/Crn and lower myostatin were associated with lower motor performance and predictive of concurrent functional performance when combined with age. Future studies are needed to more precisely determine the context of use of these biomarkers. Show less
Velde, N.M. van de; Koeks, Z.; Signorelli, M.; Verwey, N.; Overzier, M.; Bakker, J.A.; ... ; Niks, E.H. 2023
Background and ObjectivesThe slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3... Show moreBackground and ObjectivesThe slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3 muscle-enriched biomarkers in serum of patients with BMD over 4-year time and studied associations with disease severity, disease progression, and dystrophin levels in BMD.MethodsWe quantitatively measured creatine kinase (CK) using the International Federation of Clinical Chemistry reference method, creatine/creatinineratio (Cr/Crn) using liquid chromatography–tandem mass spectrometry, and myostatin with ELISA in serum and assessed functional performance using the North Star Ambulatory Assessment (NSAA), 10-meter run velocity (TMRv), 6-Minute Walking Test (6MWT), and forced vital capacity in a 4-year prospective natural history study. Dystrophin levels were quantified in the tibialis anterior muscle using capillary Western immunoassay. The correlation between biomarkers, age, functional performance, mean annual change, and prediction of concurrent functional performance was analyzed using linear mixed models.ResultsThirty-four patients with 106 visits were included. Eight patients were nonambulant at baseline. Cr/Crn and myostatin were highly patient specific (intraclass correlation coefficient for both = 0.960). Cr/Crn was strongly negatively correlated, whereas myostatin was strongly positively correlated with the NSAA, TMRv, and 6MWT (Cr/Crn rho = −0.869 to −0.801 and myostatin rho = 0.792 to 0.842, all p < 0.001). CK showed a negative association with age (p = 0.0002) but was not associated with patients' performance. Cr/Crn and myostatin correlated moderately with the average annual change of the 6MWT (rho = −0.532 and 0.555, p = 0.02). Dystrophin levels did not correlate with the selected biomarkers nor with performance. Cr/Crn, myostatin, and age could explain up to 75% of the variance of concurrent functional performance of the NSAA, TMRv, and 6MWT.DiscussionBoth Cr/Crn and myostatin could potentially serve as monitoring biomarkers in BMD, as higher Cr/Crn and lower myostatin were associated with lower motor performance and predictive of concurrent functional performance when combined with age. Future studies are needed to more precisely determine the context of use of these biomarkers. Show less
Velde, N.M. van de; Koeks, Z.; Signorelli, M.; Verwey, N.; Overzier, M.; Bakker, J.A.; ... ; Niks, E.H. 2023
Background and Objectives The slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3... Show moreBackground and Objectives The slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3 muscle-enriched biomarkers in serum of patients with BMD over 4-year time and studied associations with disease severity, disease progression, and dystrophin levels in BMD.Methods We quantitatively measured creatine kinase (CK) using the International Federation of Clinical Chemistry reference method, creatine/creatinineratio (Cr/Crn) using liquid chromatography–tandem mass spectrometry, and myostatin with ELISA in serum and assessed functional performance using the North Star Ambulatory Assessment (NSAA), 10-meter run velocity (TMRv), 6-Minute Walking Test (6MWT), and forced vital capacity in a 4-year prospective natural history study. Dystrophin levels were quantified in the tibialis anterior muscle using capillary Western immunoassay. The correlation between biomarkers, age, functional performance, mean annual change, and prediction of concurrent functional performance was analyzed using linear mixed models.Results Thirty-four patients with 106 visits were included. Eight patients were nonambulant at baseline. Cr/Crn and myostatin were highly patient specific (intraclass correlation coefficient for both = 0.960). Cr/Crn was strongly negatively correlated, whereas myostatin was strongly positively correlated with the NSAA, TMRv, and 6MWT (Cr/Crn rho = −0.869 to −0.801 and myostatin rho = 0.792 to 0.842, all p < 0.001). CK showed a negative association with age (p = 0.0002) but was not associated with patients' performance. Cr/Crn and myostatin correlated moderately with the average annual change of the 6MWT (rho = −0.532 and 0.555, p = 0.02). Dystrophin levels did not correlate with the selected biomarkers nor with performance. Cr/Crn, myostatin, and age could explain up to 75% of the variance of concurrent functional performance of the NSAA, TMRv, and 6MWT.Discussion Both Cr/Crn and myostatin could potentially serve as monitoring biomarkers in BMD, as higher Cr/Crn and lower myostatin were associated with lower motor performance and predictive of concurrent functional performance when combined with age. Future studies are needed to more precisely determine the context of use of these biomarkers. Show less
Antisense-mediated exon skipping is one of the most promising therapeutic strategies for Duchenne muscular dystrophy (DMD), and some antisense oligonucleotide (ASO) drugs have already been approved... Show moreAntisense-mediated exon skipping is one of the most promising therapeutic strategies for Duchenne muscular dystrophy (DMD), and some antisense oligonucleotide (ASO) drugs have already been approved by the US FDA despite their low efficacy. The potential of this therapy is still limited by several challenges, including the reduced expression of the dystrophin transcript and the strong 5'-3' imbalance in mutated transcripts. We therefore hypothesize that increasing histone acetylation using histone deacetylase inhibitors (HDACi) could correct the transcript imbalance, offering more available pre-mRNA target and ultimately increasing dystrophin rescue. Here, we evaluated the impact of such a combined therapy on the Dmd transcript imbalance phenomenon and on dystrophin restoration levels in mdx mice. Analysis of the Dmd transcript levels at different exon-exon junctions revealed a tendency to correct the 5'-3' imbalance phenomenon following treatment with HDACi. Significantly higher levels of dystrophin restora-tion (up to 74% increase) were obtained with givinostat and valproic acid compared with mice treated with ASO alone. Additionally, we demonstrate an increase in H3K9 acetylation in human myocytes after treatment with valproic acid. These findings indicate that HDACi can improve the therapeutic potential of exon-skipping approaches, offering promising perspectives for the treatment of DMD. Show less
Velde, N.M. van de; Koeks, Z.; Signorelli, M.; Verwey, N.; Overzier, M.; Bakker, J.A.; ... ; Niks, E.H. 2022
Cardiac involvement is the main cause of death in Becker muscular dystrophy (BMD). Identification of left ventricular (LV) function is crucial, but standard echocardiographic measurements such as... Show moreCardiac involvement is the main cause of death in Becker muscular dystrophy (BMD). Identification of left ventricular (LV) function is crucial, but standard echocardiographic measurements such as LV ejection fraction (LVEF) might not be sensitive enough to detect early myocardial dysfunction. We explored the value of LV global longitudinal strain (GLS) as a more accurate echocardiographic parameter to detect and monitor LV dysfunction in BMD. Furthermore, we studied possible factors associated with LV dysfunction and progression. A total of 40 patients with BMD (age 39.0 +/- 13.2 years) and 21 matched controls were included. Clinical variables, pulmonary tests, serum biomarkers, and echocardiograms were collected at baseline and after 2 years. LV systolic function was assessed by LVEF and LV GLS; a significant progression in LV dysfunction was defined as an absolute LV GLS deterioration >= 15%. Responsiveness to cardiac disease progression was determined using standardized response means. Patients showed impaired LVEF and LV GLS compared with controls (p <0.001). Of interest, 31 patients (77.5%) showed impaired LV GLS (defined as greater than -18%), whereas only 24 patients (60%) had reduced LVEF. LV GLS and LVEF correlated with troponin I (r = 0.553 and -0.523) and N-terminal pro-b-type natriuretic peptide (r = 0.506 and -0.585), but not with skeletal muscle or pulmonary function. At follow-up (2.0 +/- 0.5 years, n = 29), LV GLS worsened significantly (-1.3 +/- 0.8%, p = 0.002, standardized response mean = 0.70, annually = 0.60%), whereas LVEF remained stable. No risk factors for LV dysfunction progression were identified. In BMD, LV GLS is frequently impaired and shows deterioration over time compared with LVEF. LV GLS could be used as a more sensitive parameter to identify and monitor LV dysfunction. (c) 2021 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license Show less
Background: In Becker muscular dystrophy evidence for neurocognitive and behavioral comorbidity is evolving. More insight into the extend of these problems is of great importance for early... Show moreBackground: In Becker muscular dystrophy evidence for neurocognitive and behavioral comorbidity is evolving. More insight into the extend of these problems is of great importance for early detection and remediation in clinical practice.Objective: In this study we aimed to describe the neurocognitive and behavioral features of a Dutch adult cohort of BMD patients, and to evaluate correlations to motor function outcomes.Methods: 28 adult BMD patients were included. Intelligence, executive functioning, verbal memory and reaction times were assessed cross-sectionally. Additionally, patients completed questionnaires on behavioral and emotional symptoms, psychosocial and executive functions. Results were compared to normative data and correlated with disease severity as measured by the 10-meter run/walk test and Performance of the Upper Limb version 1.2.Results: 15 patients (53.6%) had a high educational level despite frequent grade repeating (48.3%) during primary or secondary school. Neuropsychological testing revealed that intellectual abilities, verbal memory, processing speed and executive functioning were statistically significant below average, but still within normal range. Regarding outcomes of the behavioral questionnaires, no significant differences were reported compared to the norm population. No relevant correlations with disease severity were found.Conclusions: This cohort of adult BMD patients exhibits minor cognitive impairments and no significant behavioral problems. The lower outcomes on processing speed and verbal memory, combined with the relatively high prevalence of grade repeating during primary and secondary school, implies that these minor impairments played a role in childhood. However, the on average high educational levels suggests that they grow out of their cognitive impairments with ageing. Show less
Longitudinal and high-dimensional measurements have become increasingly common in biomedical research. However, methods to predict survival outcomes using covariates that are both longitudinal and... Show moreLongitudinal and high-dimensional measurements have become increasingly common in biomedical research. However, methods to predict survival outcomes using covariates that are both longitudinal and high-dimensional are currently missing. In this article, we propose penalized regression calibration (PRC), a method that can be employed to predict survival in such situations. PRC comprises three modeling steps: First, the trajectories described by the longitudinal predictors are flexibly modeled through the specification of multivariate mixed effects models. Second, subject-specific summaries of the longitudinal trajectories are derived from the fitted mixed models. Third, the time to event outcome is predicted using the subject-specific summaries as covariates in a penalized Cox model. To ensure a proper internal validation of the fitted PRC models, we furthermore develop a cluster bootstrap optimism correction procedure that allows to correct for the optimistic bias of apparent measures of predictiveness. PRC and the CBOCP are implemented in the R package pencal, available from CRAN. After studying the behavior of PRC via simulations, we conclude by illustrating an application of PRC to data from an observational study that involved patients affected by Duchenne muscular dystrophy, where the goal is predict time to loss of ambulation using longitudinal blood biomarkers. Show less
Signorelli, M.; Spitali, P.; Szigvarto, C.A.; Tsonoka, R. 2021
Longitudinal and high-dimensional measurements have become increasingly common in biomedical research. However, methods to predict survival outcomes using covariates that are both longitudinal and... Show moreLongitudinal and high-dimensional measurements have become increasingly common in biomedical research. However, methods to predict survival outcomes using covariates that are both longitudinal and high-dimensional are currently missing. In this article, we propose penalized regression calibration (PRC), a method that can be employed to predict survival in such situations. PRC comprises three modeling steps: First, the trajectories described by the longitudinal predictors are flexibly modeled through the specification of multivariate mixed effects models. Second, subject-specific summaries of the longitudinal trajectories are derived from the fitted mixed models. Third, the time to event outcome is predicted using the subject-specific summaries as covariates in a penalized Cox model. To ensure a proper internal validation of the fitted PRC models, we furthermore develop a cluster bootstrap optimism correction procedure that allows to correct for the optimistic bias of apparent measures of predictiveness. PRC and the CBOCP are implemented in the R package pencal, available from CRAN. After studying the behavior of PRC via simulations, we conclude by illustrating an application of PRC to data from an observational study that involved patients affected by Duchenne muscular dystrophy, where the goal is predict time to loss of ambulation using longitudinal blood biomarkers. Show less
Time-course RNAseq experiments, where tissues are repeatedly collected from the same subjects, e.g. humans or animals over time or under several different experimental conditions, are becoming more... Show moreTime-course RNAseq experiments, where tissues are repeatedly collected from the same subjects, e.g. humans or animals over time or under several different experimental conditions, are becoming more popular due to the reducing sequencing costs. Such designs offer the great potential to identify genes that change over time or progress differently in time across experimental groups. Modelling of the longitudinal gene expression in such time-course RNAseq data is complicated by the serial correlations, missing values due to subject dropout or sequencing errors, long follow up with potentially non-linear progression in time and low number of subjects. Negative Binomial mixed models can address all these issues. However, such models under the maximum likelihood (ML) approach are less popular for RNAseq data due to convergence issues (see, e.g. [1]). We argue in this paper that it is the use of an inaccurate numerical integration method in combination with the typically small sample sizes which causes such mixed models to fail for a great portion of tested genes. We show that when we use the accurate adaptive Gaussian quadrature approach to approximate the integrals over the random-effects terms, we can successfully estimate the model parameters with the maximum likelihood method. Moreover, we show that the boostrap method can be used to preserve the type I error rate in small sample settings. We evaluate empirically the small sample properties of the test statistics and compare with state-of-the-art approaches. The method is applied on a longitudinal mice experiment to study the dynamics in Duchenne Muscular Dystrophy. Contact:s.tsonaka@lumc.nlRoula Tsonaka is an assistant professor at the Medical Statistics, Department of Biomedical Data Sciences, Leiden University Medical Center. Her research focuses on statistical methods for longitudinal omics data. Pietro Spitali is an assistant professor at the Department of Human Genetics, Leiden University Medical Center. His research focuses on the identification of biomarkers for neuromuscular disorders. Show less
Tsonaka, R.; Seyer, A.; Aartsma-Rus, A.; Spitali, P. 2021
Duchenne muscular dystrophy (DMD) is a rare genetic disorder affecting paediatric patients. The disease course is characterized by loss of muscle mass, which is rapidly substituted by fibrotic and... Show moreDuchenne muscular dystrophy (DMD) is a rare genetic disorder affecting paediatric patients. The disease course is characterized by loss of muscle mass, which is rapidly substituted by fibrotic and adipose tissue. Clinical and preclinical models have clarified the processes leading to muscle damage and myofiber degeneration. Analysis of the fat component is however emerging as more evidence shows how muscle fat fraction is associated with patient performance and prognosis. In this article we aimed to study whether alterations exist in the composition of lipids in plasma samples obtained from mouse models. Analysis of plasma samples was performed in 4 mouse models of DMD and wild-type mice by LC-MS. Longitudinal samplings of individual mice covering an observational period of 7 months were obtained to cover the different phases of the disease. We report clear elevation of glycerolipids and glycerophospholipids families in dystrophic mice compared to healthy mice. Triacylglycerols were the strongest contributors to the signatures in mice. Annotation of individual lipids confirmed the elevation of lipids belonging to these families as strongest discriminants between healthy and dystrophic mice. A few sphingolipids (such as ganglioside GM2, sphingomyelin and ceramide), sterol lipids (such as cholesteryl oleate and cholesteryl arachidonate) and a fatty acyl (stearic acid) were also found to be affected in dystrophic mice. Analysis of serum and plasma samples show how several lipids are affected in dystrophic mice affected by muscular dystrophy. This study sets the basis to further investigations to understand how the lipid signature relates to the disease biology and muscle performance. Show less
We propose a top-down approach for pathway analysis of longitudinal metabolite data. We apply a score test based on a shared latent process mixed model which can identify pathways with... Show moreWe propose a top-down approach for pathway analysis of longitudinal metabolite data. We apply a score test based on a shared latent process mixed model which can identify pathways with differentially progressing metabolites. The strength of our approach is that it can handle unbalanced designs, deals with potential missing values in the longitudinal markers, and gives valid results even with small sample sizes. Contrary to bottom-up approaches, correlations between metabolites are explicitly modeled leveraging power gains. For large pathway sizes, a computationally efficient solution is proposed based on pseudo-likelihood methodology. We demonstrate the advantages of the proposed method in identification of differentially expressed pathways through simulation studies. Finally, longitudinal metabolite data from a mice experiment is analyzed to demonstrate our methodology. Show less