In vivo MRI is sensitive to remyelination in a nonhuman primate model of multiple sclerosis

  1. Maxime Donadieu
  2. Nathanael J Lee
  3. María I Gaitán
  4. Seung-Kwon Ha
  5. Nicholas J Luciano
  6. Snehashis Roy
  7. Benjamin Ineichen
  8. Emily C Leibovitch
  9. Cecil C Yen
  10. Dzung L Pham
  11. Afonso C Silva
  12. Mac Johnson
  13. Steve Jacobson
  14. Pascal Sati  Is a corresponding author
  15. Daniel S Reich  Is a corresponding author
  1. Translational Neuroradiology Section, Laboratory of Functional and Molecular Imaging, National Institute of Neurological Disorders and Stroke, National Institutes of Health, United States
  2. Department of Neurology and Neurological Sciences, Stanford University, United States
  3. Department of Neurobiology, University of Pittsburgh, United States
  4. Section on Neural Function, National Institute of Mental Health, National Institutes of Health, United States
  5. Department of Neuroradiology, Clinical Neuroscience Center, University Hospital Zurich, University of Zurich, Switzerland
  6. Viral Immunology Section, Laboratory of Functional and Molecular Imaging, National Institute of Neurological Disorders and Stroke, National Institutes of Health, United States
  7. Cerebral Microcirculation Section, Laboratory of Functional and Molecular Imaging, National Institute of Neurological Disorders and Stroke, National Institutes of Health, United States
  8. Department of Radiology and Radiological Sciences, Uniformed Services University of the Health Sciences, United States
  9. Vertex Pharmaceuticals Incorporated, United States
  10. Neuroimaging Program, Department of Neurology, Cedars Sinai, United States
6 figures, 5 tables and 11 additional files

Figures

Example of an early active demyelinating EAE lesion with persistent hyperintensity on PDw MRI and active inflammatory demyelination on histopathology.

(A) In vivo PDw MRI acquired before EAE induction (baseline) and at various timepoints leading up to the terminal scan. Images were processed as described in Methods. Red arrows: focal white matter …

Example of a chronic, at least partially demyelinated EAE lesion core with long-lasting hyperintensity on serial PDw MRI and complete loss of myelin on histopathology.

(A) In vivo PDw MRI acquired before EAE induction (baseline) and at various timepoints leading up to the terminal scan. Images were processed as described in Methods. Red arrows: focal white matter …

Example of a nearly complete remyelinated EAE lesion with initial hyperintensity that returned to isointensity on serial PDw MRI.

(A) In vivo PDw MRI acquired before EAE induction (baseline) and at various timepoints leading up to the terminal scan. Images were processed as described in Methods. Red arrows: focal white matter …

Evolution of in vivo PDw MRI signal intensity of remyelinating EAE lesions shows a typical 4- to 9-week time course of remyelination.

Vertical axis: mean PDw signal intensity relative to gray matter. Horizontal axis: weeks post-immunization. Blue line corresponds to mean normalized PDw signal intensity (with standard deviation) of …

Histological quantification of myelin highlights different patterns across lesion categories and NAWM.

Cells were quantified in a region of interest of 1500 µm2 centered on each lesion core. Percentage of demyelination quantified for each lesion category and NAWM on PLP (A) and LFB (B). Number of …

Figure 6 with 2 supplements
In vivo (baseline) multicontrast marmoset brain magnetic resonance imaging (MRI).

The T1 subtraction image was obtained by voxelwise subtraction of pre-gadolinium from post-gadolinium T1-weighted (T1w) images. The magnetization transfer ratio (MTR) image was derived voxelwise as …

Figure 6—figure supplement 1
Workflow of the magnetic resonance imaging (MRI) postprocessing pipeline.

PD = proton density; MT = magnetization transfer; MTR = magnetization transfer ratio; N3=bias field correction.

Figure 6—figure supplement 2
Postprocessed proton-density weighted images used as input for the convolutional neural network-based lesion segmentation.

All hyperintense lesions that appeared on the follow-up magnetic resonance imaging (MRI) performed post-immunization (in this example, M#1, 22 weeks post-immunization) were captured by the automated …

Tables

Table 1
Classification of focal lesions by proton density-weighted magnetic resonance imaging.
AnimalPredicted early activePredicted chronic, at least partially demyelinatedPredicted remyelinatedTotal
M#11269
M#21168
M#31146
M#40123
M#53306
M#66208
Total12101840
Table 2
Classification of focal lesions by histopathology.
AnimalEarly activeChronic, at least partially demyelinatedRemyelinatedTotal
M#21348
M#31146
M#40213
M#53306
M#67108
Total1210931
Table 3
Confusion matrix for PDw MRI prediction of demyelinated vs. remyelinated lesions compared to histology.
n=31Predicted demyelination by MRIPredicted remyelination by MRI
Actual demyelination by histology19322
Actual remyelination by histology099
1912
Table 4
Demographics and experimental information for the six marmosets included in this study.

Immunizations used human white matter homogenate. Experiment duration corresponds to the time between immunization and terminal MRI.

AnimalSexAge (years)First immunizationSecond immunizationCorticosteroid treatmentExperiment duration (weeks)
M#1*Male3.7100 mg200 mgYes40
M#2*Male3.7100 mg200 mgNo54
M#3Male6.2200 mgYes61
M#4Male6.2200 mgNo60
M#5Female2.8200 mgYes18
M#6Female2.8200 mgNo16
  1. Denote the three different pairs of twin animals.

Table 5
Main parameters used for the different MRI contrasts acquired in vivo.
MRI contrastPDwT1wT2wT2*wMTR*
Sequence2D RARE2D MDEFT2D RARE2D MGE3D MGE
FOV (mm)32×2432×2432×2432×2438.4×38.4
Matrix214×160214×160214×160214×160256×256
Number of slices3636363636
Slice thickness (mm)11111
TR (ms)230012.58000215020
TE (ms)164.272185
TI (ms)N/A1200N/AN/AN/A
ETL1N/A4N/AN/A
Excitation pulse (shape, FA)Sinc3, 90°Sinc3, 12°Sinc3, 90°Sinc3, 70°Sinc3, 10°
Refocusing pulse (shape, FA)Sinc3, 180°N/ASinc3, 180°N/AN/A
Preparation pulse (type, shape, FA, offset)N/AExcitation/inversion, Sech, 90°/180°N/AN/AMT, Gauss, 1500 Hz
NEX11212
AT7 min 40 s6 min 56 s13 min 20 s7 min 10 s7 min 58 s
  1. FOV = field of view; TR = repetition tine; TE = echo time; TI = inversion time; ETL = echo train length (or RARE factor); FA = flip angle; NEX = number of repetitions; AT = acquisition time; RARE = rapid acquisition with relaxation and enhancement; MDEFT = modified drive equilibrium Fourier transform; MGE = multi-gradient echo.

  2. *

    Sequence was performed twice: with (MSAT) and without (M0) the MT pre-pulse.

Additional files

Supplementary file 1

Immunohistochemistry methodology.

For each of the immunohistochemical targets, respective source, clonalities and hosts, and methods for antigen retrieval, blocking, and primary and secondary antibody inoculation are listed. HIER = heat-induced epitope retrieval; RT = room temperature; P=polyclonal antibody; M=monoclonal antibody; Rb = rabbit; Ms = mouse

https://cdn.elifesciences.org/articles/73786/elife-73786-supp1-v2.docx
Transparent reporting form
https://cdn.elifesciences.org/articles/73786/elife-73786-transrepform1-v2.docx
Source data 1

First part of serial in vivo magnetic resonance images of animal #1, including proton density-weighted sequence, used for analysis and figure generation.

Uploaded into two separate source data for size of files. Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data1-v2.zip
Source data 2

Second part of serial in vivo magnetic resonance images of animal #1, including proton density-weighted sequence, used for analysis and figure generation.

Uploaded into two separate source data for size of files. Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data2-v2.zip
Source data 3

First part of serial in vivo magnetic resonance images of animal #2, including proton density-weighted sequence, used for analysis and figure generation.

Uploaded into two separate source data for size of files. Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data3-v2.zip
Source data 4

Second part of serial in vivo magnetic resonance images of animal #2, including proton density-weighted sequence, used for analysis and figure generation.

Uploaded into two separate source data for size of files. Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data4-v2.zip
Source data 5

Serial in vivo magnetic resonance images of animal #3, including proton density-weighted sequence, used for analysis and figure generation.

Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data5-v2.zip
Source data 6

Serial in vivo magnetic resonance images of animal #4, including proton density-weighted sequence, used for analysis and figure generation.

Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data6-v2.zip
Source data 7

Serial in vivo magnetic resonance images of animal #5, including proton density-weighted sequence, used for analysis and figure generation.

Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data7-v2.zip
Source data 8

Serial in vivo magnetic resonance images of animal #6, including proton density-weighted sequence, used for analysis and figure generation.

Files are in neuroimaging informatics technology initiative (NIFTI) format.

https://cdn.elifesciences.org/articles/73786/elife-73786-data8-v2.zip
Source data 9

Iba1 and proteolipid protein (PLP) immunohistochemistry images of animals #2–6 used for analysis and figure generation.

https://cdn.elifesciences.org/articles/73786/elife-73786-data9-v2.pdf

Download links