|Home | About | Journals | Submit | Contact Us | Français|
Cuprizone-induced demyelination in mice is a frequently used model in preclinical multiple sclerosis research. A recent quantitative clinically-targeted MRI method, fast macromolecular proton fraction (MPF) mapping demonstrated a promise as a myelin biomarker in human and animal studies with a particular advantage of sensitivity to both white matter (WM) and gray matter (GM) demyelination. This study aimed to histologically validate the capability of MPF mapping to quantify myelin loss in brain tissues using the cuprizone demyelination model. Whole-brain MPF maps were obtained in vivo on an 11.7T animal MRI scanner from 7 cuprizone-treated and 7 control С57BL/6 mice using the fast single-point synthetic-reference method. Brain sections were histologically stained with Luxol Fast Blue (LFB) for myelin quantification. Significant (p<0.05) demyelination in cuprizone-treated animals was found according to both LFB staining and MPF in all anatomical structures (corpus callosum, anterior commissure, internal capsule, thalamus, caudoputamen, and cortex). MPF strongly correlated with quantitative histology in all animals (r=0.95, p<0.001) as well as in treatment and control groups taken separately (r=0.96, p=0.002 and r=0.93, p=0.007, respectively). Close agreement between histological myelin staining and MPF suggests that fast MPF mapping enables robust and accurate quantitative assessment of demyelination in both WM and GM.
Macromolecular proton fraction (MPF) is a biophysical parameter defined within the two-pool model of magnetization transfer (MT)1 that describes the amount of macromolecular protons characterized by solid-like nuclear magnetic resonance (NMR) spin dynamics and involved into magnetization exchange with free water protons in biological systems. During past decade, MPF has attracted remarkable attention as a quantitative biomarker of myelin due to high sensitivity to demyelination and strong correlations between MPF and myelin content in neural tissues reported in a number of studies2,3,4,5,6,7,8,9,10,11,12. Associations between MPF and myelination were established in a variety of models and objects including excised demyelinated peripheral nerve2, postmortem multiple sclerosis (MS) brain3, rat experimental autoimmune encephalomyelitis4, shiverer mouse5,6, lipopolysacharide-induced focal lesion in the rat brain7, shaking pup8, cuprizone-induced demyelination in mice9,10, as well as normal rat brain11 and spinal cord12.
MPF or its associated quantity, pool size ratio (PSR) can be measured by a variety of quantitative MT (qMT) techniques based on semi-selective perturbation of either macromolecular or water proton populations13. Common approaches allowing reconstruction of MPF maps from magnetic resonance imaging (MRI) data include Z-spectroscopic imaging based on the analysis of signal dependence on the offset frequency and power of an off-resonance saturation pulse14,15 and on-resonance techniques based on the analysis of bi-exponential longitudinal relaxation16,17. These methods allow measurements of several parameters of the two-pool model including MPF and require collection of multiple image sets that results in time-consuming data acquisition rendering them inapplicable in clinics. An alternative fast method allowing whole-brain MPF mapping based on a single MT-weighted image, reference image, and a T1 map was recently proposed18. This technique was tested in multiple sclerosis (MS)19 and mild traumatic brain injury20 studies and demonstrated a promise as a robust clinically-targeted myelin imaging approach. In the above clinical studies, the fast single-point MPF mapping method18 demonstrated the capability to detect subtle microscopic demyelination in normal-appearing brain tissues including both white matter (WM) and gray matter (GM) and revealed the primary clinical relevance of cortical GM demyelination in MS19. Time efficiency of single-point MPF mapping has been further improved based on the recent synthetic-reference image reconstruction algorithm21,22, which requires only three source images and enables generation of three-dimensional (3D) MPF maps with high spatial resolution, whole-brain coverage, and clinically acceptable scan time. While the fast single-point MPF mapping method based on either synthetic21,22 or actually acquired18 reference image showed highly consistent MPF measurements with the standard Z-spectroscopic method in both human18 and animal22 brain in vivo, this approach is still lacking direct histological validation as a quantitative myelin imaging modality.
Cuprizone-induced demyelination in mice is a widely used model of MS that reproduces certain pathological aspects of the disease, such as depletion of oligodendrocytes, myelin degradation, microglial activation, and astrogliosis accompanied with inflammation, hydrocephalus, and brain atrophy23,24,25,26. Cuprizone administration typically results in global demyelination of WM and GM with a variable extent across anatomical structures and without formation of macroscopic focal lesions23,24,25,26. This feature makes the cuprizone model a very suitable tool for evaluation and validation of quantitative imaging methods intended to measure microscopic demyelination in normal-appearing brain tissues. The overall goal of this study was to validate the fast synthetic-reference MPF mapping method using quantitative histology in the murine cuprizone demyelination model.
Example MPF maps and LFB stained sections of control and cuprizone-treated murine brains are presented in Fig. 1. MPF maps of control animals demonstrated excellent tissue contrast between WM and GM and clear delineation of anatomical structures. Cuprizone-treated mice had visually reduced tissue contrast on MPF maps, especially in the corpus callosum and a slight overall reduction of the image intensity that can be appreciated when compared to control maps presented in the same absolute scale (Fig. 1). A visible reduction in the optical density on LFB stained brain sections of mice treated with cuprizone also can be appreciated in both WM and GM structures (Fig. 1). Quantitative comparison of MPF values and LFB optical densities in a series of anatomical regions between cuprizone-treated and control animals is presented in Table 1. Both MPF and LFB staining intensity in all WM and GM anatomical structures demonstrated a statistically significant decrease in the cuprizone group. Figure 2 compares the effect sizes characterizing group differences across anatomical structures. According to both MPF and LFB data, the largest effect of demyelination was found in the corpus callosum. Despite smaller relative changes in MPF, MPF mapping and LFB staining demonstrated a comparable performance in discriminating cuprizone-treated animals from the control group in terms of the effect size due to smaller variability of MPF measurements.
Results of linear regression analyses of MPF values as a function of LFB optical density for individual and averaged across anatomical structures data are presented in Table 2, and corresponding plots are shown in Fig. 3. In the pooled sample including data for each anatomical structure in each animal at each location from bregma, MPF values and LFB optical densities were strongly correlated (r=0.88, p<0.001) (Fig. 3a). This correlation became even stronger (r=0.95, p<0.001) after averaging data across animals for each structure (Fig. 3b). Regression analyses for separate groups of cuprizone-treated and control animals demonstrated remarkably similar results in terms of correlation strengths and coefficients of regression equations (Table 2). Particularly, for the subsamples of individual animals and structures, correlation coefficients were r=0.88 (p<0.001) for the cuprizone group and r=0.87 (p<0.001) for the control group. No significant differences between groups were found for correlation coefficients (p=0.78) and slopes (p=0.80) and intercepts (p=0.79) of the regression lines. Very strong correlations were also observed for average structure data in both groups (r=0.96, p=0.002 for cuprizone-treated animals and r=0.93, p=0.007 for control animals) with no significant differences in regression coefficients (p=0.58 for slopes and 0.56 for intercepts) and correlation strengths (p=0.72).
The extent of demyelination in the corpus callosum was evaluated additionally in the rostral-caudal and medial-lateral directions. The results of the analysis of variance (ANOVA) for different regions of the corpus callosum are presented in Fig. 4. The patterns of the cuprizone treatment effect on regional MPF values and LFB optical density appeared substantially similar (Fig. 4). The significant main effect of treatment was observed for both MPF (F=39.67, p<0.001) and LFB optical density (F=13.76, p=0.006) measurements in all regions (Fig. 4). In the medial-lateral direction, no significant main effect of the region (F=0.18, p=0.68 for MPF and F=0.03, p=0.86 for LFB optical density) was observed, but significant interactions between the effects of treatment and region were found for both MPF (F=12.26, p=0.004) and LFB (F=13.76, p=0.006) data. In the rostral-caudal direction, the significant main effect of the region was identified for both MPF (F=144.37, p<0.001) and LFB (F=7.64, p<0.02), while interactions between the effects of treatment and region did not reach statistical significance (F=1.74, p=0.21 for MPF and F=2.84, p=0.13 for LFB). Interactions between factors of the rostral-caudal and medial-lateral direction were also insignificant (F=0.96, p=0.35 for MPF and F=0.03, p=0.87 for LFB).
Bland–Altman plots for repeated MPF measurements in a series of anatomic structures obtained from two consecutive MRI scans are presented in Fig. 5. Detailed report of the statistical analysis of repeatability is provided as a Supplementary Material (Supplementary Table S1). MPF values in all anatomic structures demonstrated close agreement with no significant bias between two measurements (p>0.1 for all structures) and equal variances within each measurement (p>0.1 for all structures). Within-subject coefficients of variation (CoV) for specific structures were as follows: 1.2% for the corpus callosum, 4.9% for the anterior commissure, 1.8% for the internal capsule, 1.0% for the thalamus, 1.9% for the caudoputamen, and 1.9% for the cortex.
The key finding of this study is very close agreement between MPF measured by the fast single-point method and quantitative histology in both normal and demyelinated mouse brain across a range of WM and GM anatomical structures. This agreement is confirmed by both strong correlations between MPF and histology and comparable effect sizes for detecting demyelination in both WM and GM. Fast MPF mapping and LFB staining in the corpus callosum also demonstrated very similar capability to measure regional variations in demyelination within a single anatomic structure. The results of this study corroborate earlier reports of strong correlations between MPF and LFB optical density in the normal rodent brain11 and the capability of MPF to detect cuprizone-induced demyelination not only in WM but also in cortical GM9. Of particular note, the strengths of correlations between MPF and LFB staining across WM and GM structures reported in this study appeared essentially similar to the data for the normal rat brain11, as the correlation coefficients larger than 0.9 were found in both studies. It should be emphasized that the patterns and strengths of correlations between MPF and LFB staining discovered in this study were nearly identical for the groups of control and cuprizone-treated animals. This finding suggests that demyelination provides the primary pathological substrate of an MPF decrease in the cuprizone model, whereas potential influences on MPF of other factors that might affect various quantitative MRI measurements (such as inflammation and gliosis) are minor, if not negligible. While not fully generalizable to MPF changes in focal MS lesions in humans due to a possible confounding effect of edema, this conclusion provides a uniform conceptual framework for the interpretation of MPF as a biomarker of demyelination in normal appearing brain tissues. A recent study19 have identified a significant decrease of MPF in normal-appearing WM and GM in MS and close associations between MPF in normal-appearing brain tissues and clinical status. Furthermore, MPF changes in normal appearing brain tissues demonstrated a much higher clinical relevance than MPF in lesions or lesion volume19. In accordance with a commonly accepted view of MS as a whole-brain disease27, the development and validation of sensitive and specific quantitative biomarkers of visually undetectable demyelination in normal-appearing brain tissues may bring both studies of new therapies and individual clinical care in MS to a principally new level. Taken together with earlier pre-clinical and clinical observations9,11,19,20, the results of this study suggest that MPF provides a robust quantitative measure of the myelin content in both WM and GM, while the technical simplicity, short scan time, and whole-brain coverage enabled by the single-point synthetic-reference MPF mapping method21,22 make MPF a sensible biomarker for a variety of pre-clinical and clinical studies.
Cuprizone-induced demyelination of WM in mice has been studied by a variety of MRI methods9,10,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42. A vast majority of these studies were primarily focused on the corpus callosum, the site of most prominent myelin loss23,24,25,26, which can be easily detected by routine T2- and T1-weighted imaging28,29,30,31 and a number of quantitative techniques including diffusion tensor imaging (DTI)9,29,32,33,34,35,36,37,38,39, conventional MT ratio (MTR) imaging9,28,34,36,40,41,42, qMT9,10, and multi-component relaxation analysis38. Several MRI studies also investigated demyelination in other WM structures, such as cerebral peduncles32, anterior commissure32, optic tracts32, external capsule9,10,30,32, and cingulum30,32 which are affected by cuprizone intoxication to a lesser extent32,43. Of note, demyelination of the corpus callosum in the cuprizone model exhibits substantial spatial and temporal variability29,39,42,43,44,45,46. Both MPF mapping and histology in our study confirmed the presence of the medial-to-lateral gradient of demyelination that is in agreement with the earlier publication42. However, in contrast to previous studies29,39,42,43,44,45,46, we did not find significant distinctions in the myelin loss in the rostral-caudal direction. The most probable explanation of this discrepancy is a substantially longer cuprizone treatment period (8 weeks) in our study compared to earlier ones where the rostral-caudal gradient of demyelination was detected after 4–6 weeks of cuprizone administration29,39,42,43,44,45,46. As such, regional distinction in the myelin content observable within a certain time frame might have been equalized due to either more extensive demyelination or concurrent remyelination23,39,46. In agreement with this explanation, a reduction of rostral-caudal differences in demyelination was observed during prolonged cuprizone treatment (up to12 weeks)39. Another possible reason is inherent spatial variations of the myelin content in the normal corpus callosum as indicated by both MPF mapping and LFB staining (Fig. 4), which might mask relatively subtle regional differences in the demyelination. While spatial patterns of cuprizone-induced demyelination generally have a complex dependence on the treatment dose, duration, and mouse strain29,39,42,43,44,45,46, our study indicates that MPF provides nearly identical to histology results in quantifying a variable extent of myelin loss in WM.
Cortical and subcortical GM demyelination is a characteristic feature of the cuprizone model according to histological studies47,48,49,50,51, though MRI data about GM in cuprizone-treated animals are scarce. Particularly, a decrease of MTR was reported in deep GM41, while MTR could not distinguish cortical demyelination between treatment and control groups9,34,41. A quantitative parameter based on multicomponent T2 relaxation, myelin water fraction52,53, which has been proposed as a myelin biomarker, was unable to detect both cortical and subcortical GM demyelination, while showing a significant effect in the corpus callosum and cerebellar WM in a recent cuprizone murine ex vivo study38. Quantitative indexes derived from DTI data, such as fractional anisotropy and radial diffusivity were sensitive to cuprizone-induced changes in WM9,29,32,33,34,35,36,37,38 but failed to identify distinctions in GM9,37. A more sophisticated diffusion-based technique, diffusion kurtosis imaging (DKI) demonstrated quantitative changes in the cortex of cuprizone-treated mice37. However, these changes appeared inconsistent with changes in the corpus callosum, thus suggesting a complex mechanism determining a behavior of DKI parameters that cannot be reduced to demyelination alone37. Significant cortical changes in cuprizone-treated mice were reported for MPF, cross-relaxation rate constant, and T1 based on the analysis of Z-spectroscopic qMT data9. To the best of our knowledge, none of previous cuprizone model studies quantitatively compared any MRI technique with histology in the context of absolute myelin content measurements across a range of anatomical structures including WM and both cortical and subcortical GM. This study demonstrates that MPF measurements based on the fast single-point synthetic-reference method21 provide a similar to quantitative histology performance in the assessment of both WM and GM demyelination. A number of post-mortem histological studies have indicated a considerable amount of cortical and subcortical GM demyelination in MS in addition to WM damage27,54,55,56,57. High clinical relevance of cortical demyelination in MS has been established in a recent study19, where MPF in cortical GM was found to be the strongest predictor of disability and the secondary progressive MS disease phenotype19. As such, the histologically validated capability of fast MPF mapping to accurately quantify GM demyelination makes this method especially attractive from the clinical standpoint.
Another useful feature of fast MPF mapping in the context of both pre-clinical and clinical studies is high reliability of results. Our study of scan-rescan repeatability has demonstrated excellent performance of the method in both WM and GM structures. For MPF measurements in all explored anatomical structures, except for the anterior commissure, within-subject CoV did not exceed 2%. The MPF values in the anterior commissure had a higher variability with CoV of about 5% that is probably due to the small size and associated partial volume effect for this structure. Repeatability of fast MPF mapping in the murine brain approaches earlier data for the human brain58 where within-subject CoV of 1.5 and 1.1% were reported for WM and GM, respectively58. Slightly larger CoV in the present work can be addressed to the manual procedure of image analysis as opposed to global brain tissue segmentation in the human study58. At the same time, both this and previous58 studies indicate that the effect of instrumental factors on variability of MPF measurements is minimal.
Compared to other quantitative MRI techniques proposed for myelin characterization, fast MPF mapping offers several practical advantages. An empirical value describing the MT effect, MTR has been widely used to characterize demyelination59. While sensitivity of MTR to myelin mainly originates from its dependence on MPF11,15,19, MTR also depends on T1. MPF and T1 contribute into MTR in opposite directions, thus reducing its sensitivity to myelin content variations15,19. Additionally, the dependence on T1 makes MTR sensitive to the magnetic field strength and the effect of paramagnetic ions, such as iron, which frequently confounds demyelination, particularly in subcortical GM57,60,61. MPF itself is independent of the presence of paramagnetic substances and field strength22. Alternative semi-empirical indexes describing the MT effect, such as MT saturation62, and MTR difference observed in the recent inhomogeneous MT experimental method63 were proposed to overcome limitations of MTR. While the maps of the above quantities provide similar to MPF brain tissue contrast, these techniques have not been histologically validated, and the feasibility of their use for absolute quantitation at different field strengths remains uncertain. In contrast to DTI parameters associated with myelination (fractional anisotropy and radial diffusivity)32,33,34, MPF is insensitive to directional properties and spatial organization of WM fiber tracts, which inevitably confound DTI-based measurements and result in high anatomical heterogeneity of corresponding indices that is unrelated to the myelin content64,65,66. In contrast to relaxation-based methods aimed to extract myelin water fraction from multi-component T252, T2*67, or combined T1 and T2 data analysis53, MPF is insensitive to errors caused by iron deposition21,68 and variability of relaxation times associated with magnetic field strength22. MPF also has been shown to provide a more accurate measure of myelin than multi-component T2, since the latter is subjected to systematic errors caused by inter-compartmental water exchange12. A new relaxation-based method, relaxation along a fictitious field (RAFF) has recently emerged and showed a promise for myelin quantitation based on the comparisons with histology69. However, it remains unclear at this point, whether quantitative data obtained with this approach depend on magnetic field and paramagnetic contributions and whether a RAFF mapping technique can provide sufficient time efficiency and image quality for clinical applications. Finally, MPF maps obtained by the fast single-point synthetic-reference method21,22 offer a substantially better spatial resolution and image quality in vivo than parametric maps typically available within DTI, multi-component relaxation, and multi-parameter qMT methods.
Fast 3D MPF mapping provides a robust clinically-targeted quantitative myelin imaging modality that showed the capability to detect demyelination in normal-appearing brain tissues in recent MS and mild traumatic brain injury studies. This study provides the first quantitative histological validation of the fast MPF mapping method in the animal demyelination model. Very strong correlations between histological myelin staining and MPF suggest that the method enables accurate quantitative assessment of the myelin content in both WM and GM and that demyelination is a major or even sole pathological substrate of MPF changes in demyelinating diseases.
All animal experiments were performed in accordance with the rules adopted by the European Convention for the Protection of Vertebrate Animals used for Experimental and other Scientific Purposes. The experimental protocol was approved by the Bioethical Committee of the Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences and the Bioethical Committee of the Biological Institute at the Tomsk State University. Fourteen five-week-old C57BL/6 male mice were obtained from the vivarium of the N.N.Vorozhtzov Institute of Organic Chemistry (Novosibirsk, Russian Federation). Animals were housed with a 12hours dark-light cycle at the temperature of 21±2°C, and humidity of 40±2%. Food and water were provided ad libitum. Seven mice were treated 0.3% cuprizone solution in drinking water during 8 weeks to induce demyelination, and seven control mice received regular vivarium chow and water. Imaging was performed under isoflurane anesthesia (1.5–2% in oxygen) with respiratory monitoring during the scan.
After scanning, mice were transcardially perfused with 4% paraformaldehyde (PFA) under ether anesthesia. The brains were removed and fixed overnight in PFA at 4°C. Then the brains were cryoprotected in graded concentration of sucrose in phosphate buffer (24h at 10% and 24h at 20%) at 4°C, frozen in liquid nitrogen, and stored at −80°C. Coronal brain sections with 10μm thickness were prepared using an HM525 cryostat (Thermo Fisher Scientific, Walldorf, Germany). Histological myelin staining was performed using Luxol Fast Blue (LFB) colorant (Bio-Optica, Milano, Italy). Whole-brain sections were photographed using an Axio Imager Z2 microscope (Carl Zeiss AG, Göttingen, Germany) with the EC PlanNeofluar lx/0.025 objective and AxioVision 4.8 software. Identical imaging parameters (voltage of microscope lamp and exposure time) were set for all photographed slices.
Additional seven eight-week-old C57BL/6 male mice were used to assess scan-rescan repeatability of MPF mapping. Animals were kept under standard vivarium conditions, as specified above, and scanned twice with one week interval between scans.
Mice were imaged on an 11.7 Tesla horizontal-bore animal MRI scanner (BioSpec 117/16 USR; Bruker BioSpin, Ettlingen, Germany) with manufacturer’s four-channel mouse brain surface coil. Imaging was performed under isoflurane anesthesia (1.5–2% in oxygen) with respiratory monitoring during the scan. A fast high-resolution 3D MPF mapping protocol was implemented according to the previously described single-point method with synthetic reference image normalization21,22, and scan parameters were adapted for mouse whole-brain imaging. The protocol included the following sequences applied in the coronal plane with 3D field-of-view of 20×20×24mm:
The GRE and AFI sequences were implemented with optimal radiofrequency and gradient spoiling based on the excitation pulse phase increments of 169° for GRE and 39° for AFI72. In all sequences, linear phase-encoding order with fractional (75%) k-space acquisition in the slab selection direction was used. The total protocol execution time was about 35min.
Reconstruction of MPF maps was carried out using custom-written C-language software based on the algorithm detailed elsewhere21. In brief, this algorithm computes MPF maps from three source images (MT-, T1-, and PD-weighted) with correction of B0 and B1 field non-uniformities by voxel-based iterative single-point solution of the matrix pulsed MT equation18 with a calculated synthetic reference image for data normalization21.
The same protocol and map reconstruction was implemented to examine reproducibility of MPF mapping.
MPF maps and microphotographs of LFB stained histological sections were analyzed using publically available ImageJ software (National Institutes of Health, Bethesda, MD, USA). Two brain locations (−1.58mm and +0.74mm from bregma) defined according to a mouse brain atlas73 were chosen for quantitative analysis. MPF values and LFB optical densities were measured for a series of WM and GM structures including corpus callosum, anterior commissure, internal capsule, thalamus, caudoputamen, and cortex. Regions-of-interest (ROIs) of a standard size and shape were manually placed on MPF maps and photographs of LFB stained sections within investigated brain structures using a mouse brain atlas73 as reference. The scheme of ROI measurements is exemplified in Fig. 6. MPF values for each anatomical structure were quantified by averaging ROI measurements obtained from 3D MPF maps corresponding to each location from bregma. Myelin densities on LFB stained photographs were quantified using the procedure modified from that described by Underhill, et al.11. The mean intensity of the red channel (IR) in a ROI was measured from RGB images as a quantity characterizing the complementary blue channel saturation. Specifically, visual perception of a more saturated blue color reflects stronger light absorption in the red region of the spectrum associated with stained myelin that translates into a lower intensity of the red channel in RGB images. Additionally, the mean intensity of the red channel in the background was measured on each photograph in four ROIs outside the brain. The intensities from the four above ROIs were averaged and used as the background correction factor (IB). LFB optical density (in %) was subsequently calculated for each ROI as 100×[1−(IR/IB)], which yields a higher myelin content for increasing values. Individual ROI measurements were obtained from three adjacent sections per each location from bregma and averaged for each anatomical structure.
All statistical analyses were carried out in Statistica 10.0 for Windows (StatSoft Inc, Tulsa, OK, USA). Mean values and standard deviations (SD) of MPF and LFB optical density were calculated for each anatomical structure. Normality of the data within animal groups and residuals in regression analyses was assessed using Shapiro-Wilk test. No significant deviations from the normal distribution were found, and therefore, parametric analyses were used. MPF values and quantitative histology data were compared between cuprizone-treated and control animals using one-way multivariate analysis of variance (MANOVA) followed by independent-samples t-tests for individual anatomical structures. Distinctions in variables between animal groups were described by the effect size (Cohen’s d) calculated as the ratio of the mean difference to the pooled standard deviation.
Pearson’s correlation coefficient (r) and linear regression analysis were used to determine associations between MPF values and LFB optical densities across anatomical structures and animals. Significance of differences between Pearson’s correlation coefficients for the two groups of animals was tested using Fisher’s r-to-z transformation. Slopes and intercepts of regression equations were compared between groups using analysis of covariance (ANCOVA)74 where MPF values were used as a dependent variable, LFB optical density was used as a covariate, and group membership (control or cuprizone-treated) was used as a categorical factor.
Regional distinctions in the corpus callosum demyelination were assessed using three-way analysis of variance (ANOVA) with the factors of treatment, ROI location in the rostral-caudal direction, and ROI location in the medial-lateral direction. The model included the main effects of factors and their pairwise interactions.
An agreement between MPF maps obtained from two repeated scans was assessed using Bland–Altman plots for each anatomical structure. Significance of the bias between repeated MPF measurements was examined using the one sample t-test for the mean difference. The limits of agreement were calculated as the mean difference ±1.96SD of the mean difference. To estimate variability between repeated measurements, the within-subject coefficients of variation was calculated as the percentage ratio of SD to the mean of paired measurements. The Levene’s test was used to assess equality of variances. Two-tailed tests were used in all analyses. Statistical significance was defined as a p value less than 0.05.
How to cite this article: Yu. Khodanovich, M. et al. Histological validation of fast macromolecular proton fraction mapping as a quantitative myelin imaging method in the cuprizone demyelination model. Sci. Rep. 7, 46686; doi: 10.1038/srep46686 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors acknowledge financial support from the Russian Science Foundation (project № 14-45-00040) for experimental studies, data analysis, and manuscript preparation and from the National Institutes of Health (grant R21EB016135) for software development. V.L.Y. received partial salary support from the National Multiple Sclerosis Society grant RG 4864A1/1. Animal procedures were carried out at the Center for Genetic Resources of Laboratory Animals at the Institute of Cytology and Genetics (Siberian Branch, Russian Academy of Sciences) operating under support of the Ministry of Education and Science projects RFMEFI61914X0005 and RFMEFI62114X0010.
The authors declare no competing financial interests.
Author Contributions M.Y.K., I.V.S., and V.L.Y. designed the study. M.Y.K. performed statistical analysis and wrote the manuscript. I.V.S. and T.G.T. were responsible for the development of the animal model. A.E.A. and A.V.R. performed MRI studies. V.Y.G. performed image analysis. N.M.N.-D. and L.R.M. performed histological processing and staining. V.L.Y. developed an MRI protocol and image processing software and critically revised the manuscript. All authors reviewed and approved the final version of the manuscript.