1
|
Sweeney EM, Shinohara RT, Dewey BE, Schindler MK, Muschelli J, Reich DS, Crainiceanu CM, Eloyan A. Relating multi-sequence longitudinal intensity profiles and clinical covariates in incident multiple sclerosis lesions. Neuroimage Clin 2015; 10:1-17. [PMID: 26693397 PMCID: PMC4660378 DOI: 10.1016/j.nicl.2015.10.013] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2015] [Revised: 10/22/2015] [Accepted: 10/25/2015] [Indexed: 11/19/2022]
Abstract
The formation of multiple sclerosis (MS) lesions is a complex process involving inflammation, tissue damage, and tissue repair - all of which are visible on structural magnetic resonance imaging (MRI) and potentially modifiable by pharmacological therapy. In this paper, we introduce two statistical models for relating voxel-level, longitudinal, multi-sequence structural MRI intensities within MS lesions to clinical information and therapeutic interventions: (1) a principal component analysis (PCA) and regression model and (2) function-on-scalar regression models. To do so, we first characterize the post-lesion incidence repair process on longitudinal, multi-sequence structural MRI from 34 MS patients as voxel-level intensity profiles. For the PCA regression model, we perform PCA on the intensity profiles to develop a voxel-level biomarker for identifying slow and persistent, long-term intensity changes within lesion tissue voxels. The proposed biomarker's ability to identify such effects is validated by two experienced clinicians (a neuroradiologist and a neurologist). On a scale of 1 to 4, with 4 being the highest quality, the neuroradiologist gave the score on the first PC a median quality rating of 4 (95% CI: [4,4]), and the neurologist gave the score a median rating of 3 (95% CI: [3,3]). We then relate the biomarker to the clinical information in a mixed model framework. Treatment with disease-modifying therapies (p < 0.01), steroids (p < 0.01), and being closer to the boundary of abnormal signal intensity (p < 0.01) are all associated with return of a voxel to an intensity value closer to that of normal-appearing tissue. The function-on-scalar regression model allows for assessment of the post-incidence time points at which the covariates are associated with the profiles. In the function-on-scalar regression, both age and distance to the boundary were found to have a statistically significant association with the lesion intensities at some time point. The two models presented in this article show promise for understanding the mechanisms of tissue damage in MS and for evaluating the impact of treatments for the disease in clinical trials.
Collapse
Key Words
- Biomarker
- CI, confidence interval
- Expert rater trial
- FLAIR, fluid-attenuated inversion recovery
- Function-on-scalar regression
- Longitudinal lesion behavior
- Longitudinal study
- MRI, magnetic resonance imaging
- MS, multiple sclerosis
- Multi-sequence imaging
- Multiple sclerosis
- NAWM, normal-appearing white matter
- NINDS, National Institute of Neurological Disease and Stroke
- PC, principal component
- PCA, principal component analysis
- PD, proton density-weighted
- Principal component analysis and regression
- RRMS, relapsing remitting MS
- SPMS, secondary progressive MS
- Structural magnetic resonance imaging
- T, Tesla
- T1, T1-weighted
- T2, T2-weighted
- sd, standard deviation
Collapse
Affiliation(s)
- Elizabeth M. Sweeney
- Department of Biostatistics, The Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD 21205, United States
- Translational Neuroradiology Unit, Division of Neuroimmunology and Neurovirology, National Institute of Neurological Disease and Stroke, National Institute of Health, Bethesda, MD 20892, United States
- Corresponding author.
| | - Russell T. Shinohara
- Department of Biostatistics and Epidemiology, Center for Clinical Epidemiology and Biostatistics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, United States
| | - Blake E. Dewey
- Translational Neuroradiology Unit, Division of Neuroimmunology and Neurovirology, National Institute of Neurological Disease and Stroke, National Institute of Health, Bethesda, MD 20892, United States
| | - Matthew K. Schindler
- Translational Neuroradiology Unit, Division of Neuroimmunology and Neurovirology, National Institute of Neurological Disease and Stroke, National Institute of Health, Bethesda, MD 20892, United States
| | - John Muschelli
- Department of Biostatistics, The Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD 21205, United States
| | - Daniel S. Reich
- Department of Biostatistics, The Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD 21205, United States
- Translational Neuroradiology Unit, Division of Neuroimmunology and Neurovirology, National Institute of Neurological Disease and Stroke, National Institute of Health, Bethesda, MD 20892, United States
| | - Ciprian M. Crainiceanu
- Department of Biostatistics, The Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD 21205, United States
| | - Ani Eloyan
- Department of Biostatistics, Brown University School of Public Health, RI 02912, United States
| |
Collapse
|