Age-Related Decline in Brain Myelination: Quantitative Macromolecular Proton Fraction Mapping, T2-Flair Hyperintensity Volume, and Anti-myelin Antibodies 7 Years Apart

1. Introduction

Identifying a pathological substrate of cognitive decline in normal aging and age-related dementias becomes increasingly important with the continuous growth of life expectancy. It is commonly recognized that normal brain aging is a multifaceted process primarily affecting the microstructure of neuropil, while keeping the global neuronal population nearly intact [1,2]. Age-related decline in the brain function suggests that disruption of normal brain connectivity may be an important underlying reason. This circumstance motivated a substantial attention to investigations of myelin changes in the aging brain.

A number of histological studies in rodents [3] and non-human primates [4] have shown the effect of age on the ultrastructure of myelinated nerve fibers, such as “balloonsˮ in myelin sheaths [4,5], increase in the number of lamellae [6], accumulation of dense cytoplasm with vesicular inclusions [5]. Paler staining of white matter was observed in the brains of old primates compared to young animals [7]. Histological studies on post-mortem human brain have also shown a decrease in the intensity of staining of nerve fibers in the cerebral cortex [8]. Stereological analyses [9] indicated a 15-17% age-related loss of white matter in the cerebral hemispheres and 27% reduction in the total length of myelinated nerve fibers.

Histological studies of age-related myelin changes in humans are rare, highly dependent on the quality of the biomaterial fixation, and can only be carried out post-mortem. An alternative way that allows for longitudinal study of brain myelination is the use of non-invasive imaging. A number of magnetic resonance imaging (MRI) methods are sensitive to myelin, though their specificity varies greatly [10]. Nevertheless, modern neuroimaging studies collectively suggest a progressive decline in myelination in the aging brain. Diffusion tensor imaging (DTI) has been the most common approach to study age-related changes in white matter (WM) [11,12,13,14,15,16,17]. DTI studies demonstrated age dependence of myelin-associated indices (fractional anisotropy (FA) and radial diffusivity (RD)) consistent with demyelination in various WM tracts including frontal WM, internal capsule (IC), and corpus callosum (CC) [16,17]. A major limitation of DTI for myelin assessment is the sensitivity of FA and RD to the direction of diffusion in regions containing multidirectional fibers and changes in tissue microstructure unrelated to myelination [15]. Another, more specific to myelin approach, myelin water fraction (MWF) mapping has showed age-related decrease in myelination for most WM regions [13,14,16,17,18,19]. A semi-quantitative technique for studies of cortical myelination based on the signal ratio obtained from T1- and T2-weighted images revealed inverted U-shaped aging profiles in the cortex suggesting protracted myelination in younger ages followed by a decline in elderly subjects [20]. Similar profiles were identified later for the age dependence of R1=1/T1 [21] and MWF [13] in WM. Quantitative interpretation of the above findings strictly in terms of myelin is difficult because changes in relaxation times T1 and T2 affect both measured signal intensities and MWF [22]. Such changes are non-specific to myelin and may be caused by age-related iron deposition, water content alterations, and gliosis [22,23].

The magnetization transfer (MT) effect provides an alternative mechanism for probing myelination changes via interaction between water and macromolecules in tissues [24,25]. In early brain aging studies (detailed review can be found in [24], the MT effect was assessed by the empirical semi-quantitative index MT ratio (MTR). Several studies have shown negative associations between MTR and age, however, MTR provided lower sensitivity to age-related changes as compared with MWF and DTI [14]. The main limitation of MTR, which likely explains the lower specificity and sensitivity for myelin, is its dependence on longitudinal relaxation [26,27]. This dependence partially offsets the effect of macromolecular composition changes and additionally introduces unwanted sensitivity to paramagnetic effects. A newer MT quantitation approach, MT saturation index (MTsat) is independent of relaxation and offers improved specificity and sensitivity to age-related changes in myelination [28,29,30,31]. An alternative and more comprehensive way to characterize the MT effect is based on quantitative mapping of the fundamental parameters of the two-pool MT model [23] using specialized quantitative MT (qMT) techniques [25]. While applications of qMT in aging research are limited to date, particular parameters of the two-pool model (specifically, the macromolecular proton fraction (MPF) and cross-relaxation rate constant) showed negative correlations with age at the whole-brain level and in certain WM structures [32].

Among the parameters of the two-pool model, MPF has attracted major interest as a biomarker of myelin [25]. Multiple studies in animal models evidenced strong correlations between MPF and myelin histology and high sensitivity of this parameter to myelin damage as detailed in the recent review [25]. A recently emerged fast single-point MPF mapping method [33,34] allows reconstruction of MPF maps in isolation from other two-pool model parameters and overcomes the key limitations of the most qMT techniques related to extremely time-consuming data acquisition and high sensitivity to noise. The fast MPF mapping method [29,30] offers several practical advantages including independence of magnetic field strength [35], insensitivity to major confounders of alternative myelin imaging methods, such as myelinated fiber orientation [15] and iron deposition [36], high reproducibility [37,38,39], and easy implementation based on routine MRI equipment without modification of original manufacturers’ pulse sequences [47,48,49]. This method has been extensively validated by histology in the normal animal brain [26,37] and the animal models of cuprizone-induced demyelination and remyelination [37,40], ischemic stroke [41,42], and neonatal development [43].

Recent applications of fast MPF mapping in the studies of multiple sclerosis [27,36], mild traumatic brain injury [44], schizophrenia [39], and normal brain development in fetuses [45,46], children [47,48], and adolescents [49] have shown the feasibility of using this method as a robust and versatile tool for the quantitative assessment of demyelination and myelin development in humans.

The primary goal of this study was quantification of age-related changes in global and regional myelination of the normal adult human brain using the fast MPF mapping method. Unlike the previous neuroimaging studies where the effect of age on myelination was examined cross-sectionally at a single time point per subject, the present study is based on the longitudinal assessment of myelination in the same subjects with the seven-year time interval. Additionally, we investigated age-related changes in the concentrations of autoantibodies to major myelin proteins in blood plasma, which served as peripheral markers of myelin damage and white matter hyperintensities on T2-FLAIR images, which are commonly viewed as signs of focal demyelination.

2. Materials and Methods

2.1. Subjects

Eleven healthy volunteers participated in the study. The inclusion criteria were the following: age from 33 to 60 years, the absence of the history of traumatic brain injury and the absence of any diagnosed neurologic or psychiatric condition. The exclusion criteria: contraindications to MRI, inability to tolerate the MRI procedure, and self-withdrawal from the study. Written informed consent was obtained from all participants.

For the first MR imaging, participants were recruited from staff at the Mental Health Research Institute between April and November 2015. The second MRI imaging was performed on the same participants between September and November 2022. The demographic characteristics of participants are shown in Table 1.

2.2. MRI data acquisition

All participants underwent the first (2015) and repeated (2022) MR imaging using 1.5T clinical scanner Magnetom Essenza (Siemens, Erlangen, Germany). The MRI equipment did not undergo any software and hardware upgrades or major repairs during the study. The fast MPF mapping protocol [39] included three 3D spoiled gradient-echo pulse sequences with following acquisition parameters:

• MT-weighted: TR = 20 ms, echo time (TE) = 4.76 ms, flip angle (FA) = 8°, scan time 5 min 40 s;
• T1-weighted: TR =16 ms, TE = 4.76 ms, FA =18°, scan time 4 min 32 s;
• Proton-density-weighted: TR= 16 ms, TE = 4.76 ms, FA= 3°, scan time 4 min 32 s.

In addition, the following sequences were included in the protocol:

• 3D FLAIR-SPACE-FSE: TR = 5000 ms, TE = 390 ms, TI = 1800 ms and
• 3D T2-SPACE-FSE: TR=3000ms, TE=335ms.

All scans were acquired in the sagittal plane with a voxel size of 1.25 × 1.25 × 1.25 mm3 (matrix 192 × 192 × 160, field of view 240 × 240 × 200 mm3), single signal averaging.

The total scanning time was about 35 minutes.

2.3. Image processing

MPF maps were reconstructed using the previously developed software in the C++ language (available at https://www.macromolecularmri.org/), which implements a single-point algorithm with a synthetic reference image [33,50]. MPF maps obtained from the first and repeat MRI were examined for 7-year global and regional changes of WM and GM. In addition, age-related changes in the volume of T2-FLAIR hyperintensities and MPF values in the hyperintense areas were assessed.

Global WM and GM changes on MPF maps were assessed using the freely available MRIcro [51], ImageJ (Fiji) [52], ITK-snap [53], and FSL [54] software tools. The first stage involved skull stripping using a mask, which was obtained by applying the BET algorithm to the PD-weighted images in the MRIcro application. The mask was then converted to a binary image in ImageJ using the Threshold function and applied to the MPF maps to remove extracerebral tissue. Automatic global segmentation of MPF maps was performed using the FSL package to obtain masks of WM, GM, mixed WM-GM and mixed GM with cerebrospinal fluid (CSF) as detailed earlier [27,39]. Masks were used to measure the mean MPF values of these compartments. The measurements were carried out using ITK-SNAP application.

Regional WM and GM segmentation was carried out using Advanced Normalization Tools (ANTs) [55,56] and Eve anatomical atlas [57]. T1 template image of Eve atlas was registered to individual MPF maps using antsRegistrationSyNQuick algorithm. Then the obtained deformation field was applied to Type-III Eve atlas segmentation [57] to register the template atlas labels to individual MPF maps (Figure 1).

The measurements on MPF maps were performed for 118 GM and WM structures of right and left hemispheres using ITK-snap software. The list of structures included:

1. Juxtacortical (superficial) WM: superior parietal, superior, middle, and inferior frontal; precentral; postcentral; angular; pre-cuneus; cuneus; lingual; fusiform; superior, inferior, and middle occipital; superior, inferior, and middle temporal; lateral and middle fronto-orbital, supramarginal, rectus, cingulum (parts of cingulate gyrus and hippocampus);

2. WM pathways and fasciculi: corticospinal tract (CST); medial lemniscus; anterior limb, posterior limb, and retrolenticular part of internal capsule (IC); inferior, superior, and middle cerebellar peduncles (CP); cerebral peduncles; posterior thalamic radiation; anterior, superior, and posterior corona radiata (CR); fornix (FX) (stria terminalis, column and body); superior longitudinal (SL) fasciculus; superior (SFO) and inferior fronto-occipital (IFO) fasciculi; uncinate fasciculus; sagittal stratum; external capsule; pontine crossing tract; genu, body, and splenium of corpus callosum (CC); tapatum;

3. Subcortical and allocortical GM structures: amygdala; hippocampus; entorhinal area; caudate nucleus; putamen; globus pallidus; thalamus;

4. Brainstem structures: midbrain; pons; medulla.

The measurements for left and right hemispheres were averaged for all brain structures except for juxtacortical WM. A series of juxtacortical WM labels corresponding to specific gyri were analyzed separately from other WM structures, taking into account whether they belonged to the left or right hemisphere.

WM hyperintensities were outlined manually by two operators blinded to the subject information and scan time point on T2-FLAIR images with the guidance of T2-weighted images. Then, T2-FLAIR images were registered to MPF maps using ITK-snap software to measure mean MPF values in the outlined areas.

2.4. ELISA

Whole blood samples from each subject were taken from the cubital vein after 12-h over-night fasting and centrifuged for 20 min at 2000× g at 4 ◦C. The serum was isolated and stored at −80°C. Quantitative analyses of IgG antibodies against myelin basic protein (MBP) and proteolipid protein (PLP) in the serum were executed using respective ELISA kits by Cloud-Clone Corp. (CCC, USA) according to the manufacturer’s instructions. The absorbance of each well at 450 nm was measured on a Varioskan LUX spectrophotometer (Thermo Scientific, Waltham, MA, USA) located at core facility Medical Genomics at Tomsk National Research Center. Measurement results are presented in terms of optical density units (ODU).

2.5. Statistical analysis

Statistical analysis was performed using Statistica 10.0 software. Differences between time points, inter-hemispheric, and gender differences were analyzed using the repeated measures analysis of variance (ANOVA) followed by post-hoc Fisher LSD tests. In the analyses of multiple brain structures, p-values were adjusted using the Benjamini–Hochberg procedure for false discovery rate (FDR) correction to prevent false positive results in multiple comparisons. FDR level was set to 0.05. Associations between variables were assessed using the Pearson correlation coefficient. Statistical significance for all analyses was taken as less than 0.05.

3. Results

3.1. Age-related global changes in the brain myelination

Figure 2 demonstrates global changes in MPF and volumes of WM, GM, and mixed WM-GM compartments. Age-related decrease of myelination was significant in WM, GM, and mixed WM-GM (Figure 2b). However, the percentage MPF decrease in WM (1.30±0.55%) was twice greater than that in GM (0.75±0.72%) and mixed WM-GM (0.75±0.72%) (Figure 2с). A decrease in volume was significant only for GM with the relative change of 2.6% (Figure 2d). No gender differences in MPF were found for both WM, mixed WM-GM, and GM.

3.2. Age-related changes in separate WM and GM structures

Figure 3 shows age-related MPF differences in the explored WM and GM structures in detail. Additionally, Figure 3(a) highlights interhemispheric differences in MPF for juxtacortical WM for both time points. Interhemispheric differences, with a few exceptions, coincide for 2015 and 2022. In the frontal lobe, the MPF was higher in the right hemisphere compared to the left one (differences were significant for the rectus, middle frontal, lateral and middle fronto-orbital cortical WM for both time points, for precentral WM – only for 2015). Higher MPF values in the right hemisphere were also found for the fusiform WM, for both time points. For the parietal occipital lobes, and singulum, on the contrary, the MPF values were higher in the left hemisphere (differences are significant for precuneus, cuneus, supramarginal WM for both time points, for lingual, occipital and hippocampal part of cingulum WM – only for 2022). Temporal juxtacortical WM didn’t show significant interhemispheric differences.

In almost all structures of both WM and GM, a decrease in MPF was observed, but not in all structures it was statistically significant. All frontal juxtacortical WM areas showed significant MPF decrease both in left and right hemispheres (except for precentral WM, significant differences were only in the right hemisphere), whereas lateral occipital WM did not show significant decrease over 7 years (Figure 3a). Lateral temporal, precentral, and postcentral WM showed significant differences mostly in the right hemisphere.

Medial and inferior part of juxtacortical WM revealed some heterogeneity of age-related changes. The cingulate part of cingulum, fusiform, and rectus WM showed significant decrease in both hemispheres, whereas in the cuneus, pre-cuneus, and hippocampal part of cingulum WM differences were not observed. Lingual WM significantly decreased only in the left hemisphere. A significant decrease in MPF was also found for the majority of the investigated WM pathways (Figure 3c) including CC (except for the splenum), posterior thalamic radiations, anterior and superior CR, SL, uncinate, SFO and IFO fasciculi, anterior and posterior limbs of IC, external capsule, column and body of FX, sagittal stratum, and tapatum. No differences were found in the cerebellar pathways, cerebral peduncles, pontine crossing tract, stria terminalis, posterior part of CR, retrolenticular part of IC, and splenum of CC. In addition, MPF significantly decreased in all investigated GM and mixed WM-GM structures (Figure 3b) including all basal ganglia, thalamus, hippocampus, amygdala, entorhynal area, and brainstem (except for the pons). No significant gender effect on the amount of the MPF decrease was found.

The percentage changes between the mean MPF values obtained in 2015 and 2022 for the same brain structures are shown in Figure 4. The most apparent decrease in MPF in the juxtacortical WM (Figure 4a), by 1.8-5.4%, was observed in the frontal lobe. In other lobes of the brain, the decrease in MPF was much smaller. Only for fusiform and singular juxtacortical WM the decrease reached 2% in both hemispheres. The same magnitude of a decrease in MPF was found for the precentral and central, as well as inferotemporal juxtacortical WM of the right hemisphere. Among the WM tracts (Figure 4b), the genu of CC achieved the largest decrease of 4%. Anterior WM regions (anterior CR, anterior limb of IC, genu of CC) generally showed a higher percentage reduction in MPF (exceeding 2%) as compared to posterior WM regions. In the segmented tracts, MPF reductions were mostly significant with the largest effect in the uncinate fasciculus (2.8%). Among all brain regions studied, the greatest decrease of 5.9% was observed for the caudate nucleus (Figure 4c). In other deep GM structures, the differences were smaller and amounted to 1.6% for the globus pallidus, 1.9% for the putamen, and 2.2% for the thalamus. Allocortical structures showed a decrease of around 2% (2.4% for the amygdala, 2.0% for the hippocampus, and 2.2% for the entorhinal area). MPF in the brainstem structures decreased to the smallest extent (1.3% for midbrain, 1.2% for medulla, and 0.6% for pons).

3.3. Age-related changes in the volume of T2-FLAIR hyperintensities

Typical age-related changes of the brain reflected in T2-FLAIR images and MPF maps are shown in Figure 5. Specifically, zones of hyperintensities were detected on images obtained in 2022; those hyperintensity zones were mainly located periventricularly (Figure 5a). Over 7 years, the volume of hyperintense zones increased by an average of 35.2% (Figure 5b), both due to an increase in the volume of lesions adjacent to the lateral ventricles and due to the appearance of new lesions in bulk WM. (Figure 5a). Average MPF within T2-FLAIR hyperintensities slightly decreased but the differences between 2015 and 2022 were not statistically significant (Figure 5b).

The percentage changes in the volume of T2-FLAIR hyperintensities between men and women was close to significant (p = 0.053). In men, the percentage changes in the volume of hyperintensities from 2015 to 2022 was 20.3±20.0% in average, while in women the averaged changes reached 50.0±24.4%.

3.4. Age-related changes in myelin-related autoantibodies

In 2022, the levels of antibodies against MBP tended to increase compared to 2015 (0.146±0.109 ODU in 2015, 0.158±0.154 ODU in 2022), but this difference (8.8±1.6%) was not statistically significant. The levels of antibodies against PLP in 2022 and 2015 were near identical (0.042±0.007 ODU in 2015, 0.043±0.006 ODU in 2022, percentage changes 1.3±0.8%). No gender differences were found for both antibodies.

Significant negative correlations were found between the percentage differences in anti-MBP antibody concentration and the percentage differences in global WM (r = 0.69, p < 0.05) and mixed WM-GM (r = 0.63, p < 0.05). The correlation of percentage changes in anti-MBP antibodies with the percentage changes in the volume of WM hyperintensities was insignificant. No significant correlations were found for the age-related percentage differences in antibodies to PLP.

4. Discussion

This study for the first time has demonstrated an age-related decrease in WM and GM myelination measured by quantitative MPF mapping. Specifically, a significant decrease in MPF between repeated MRI scans 7 years apart was found for both global WM and mixed WM-GM compartments, as well as for the majority of WM and GM structures studied. The study was conducted on the same subjects using the same MRI scanner, which previously has shown highly reproducible results [39].

To the best of our knowledge, no other longitudinal studies were published so far describing age-related changes in myelin content of the same subjects examined in two consecutive scans several years apart. Typically, the publications of the age-related alterations in brain myelination report the results of the cross-sectional studies on a fairly large sample of healthy subjects of different ages [11,12,13,14,18,19,30,58,59,60]. Therefore, our current study is unique in this aspect.

Several earlier quantitative neuroimaging studies reported findings consistent with the decline of myelin content during normal aging. Particularly, DTI studies showed age-related changes in FA and RD in various WM regions and fiber tracts [11,12,61]. O’Sullivan et al. [12] reported age-related decrease in FA in the frontal WM. In the work of Salat et al. [11] age-related FA decrease was found only in the frontal WM, the posterior limb of the IC, and the genu of the CC, out of 9 studies brain areas; the FA in the temporal and posterior WM was relatively preserved. Most recent publications [13,14], reported stronger correlation of MWF with age but weaker correspondence to age of FA and RD parameters. Arshad et al. [13] investigated the relationship of DTI and MWF parameters with age within two models - linear and quadratic (inverted-U). It was found that FA and RD correlate with age only within linear model, while MWF correlates with age within inverted-U (quadratic) model, which better describes age-myelin association in the wide range of ages [13,61]. Similar results were reported by Faizy et al [14]. Specifically, it was found that MWF significantly decreased with age in most WM regions (except corticospinal tract). FA and MTR were associated with lower MWFs in the commissural fiber tracts only; RD had no regional effects on MWF. Authors concluded that DTI and MTI methods have limited specificity to myelin. Other publications that used MTI and MTR methods for the evaluation of age-dependent myelin changes [14,30,59] also showed weak dependence of MTR from age.

Our studies of myelination using MPF mapping to assess myelination are consistent with the above findings regarding myelin decline with age.

Our results highlight the capability of fast MPF mapping to capture very subtle changes in brain myelination with the use of a relatively small sample size. This feature is a direct consequence of the high precision of the fast single point MPF mapping method. Earlier studies demonstrated excellent scan-rescan repeatability of MPF measurements in brain tissues confirmed by small within-subject coefficients of variation (CV) in ranges of 1-2% in animals [37] and 0.8-1.6% in humans [38,39]. Of particular importance, the recent study performed on the same MRI unit with the same imaging protocol and global segmentation pipeline [39] reported the coefficients of variations of 0.8% in WM and 1.0% in GM in a short-term scan-rescan setting. These results corroborate our findings of a significant decrease in MPF by 0.75% in GM and 1.3% in WM. In summary, the current study affirms the feasibility of tracing a small decline in the myelin content associated with normal brain aging, which has a translation potential of detecting larger effects of myelin disruption in various neurological conditions at the individual level.

Unlike MTR, MPF does not depend on T1 in tissues [33,34], which is especially important for studying age-related changes in myelin due to iron accumulation in the older brain [31]. Our recent study on patients with multiple sclerosis (MS) proved that MPF is an iron-insensitive measure of demyelination [36]. Our current study applied MPF to longitudinal repetitive measurements of myelin in the brain of healthy adults. A significant global decrease by 1.3% (p<0.001) over 7 years was found in WM myelination in the same subjects. A two-fold smaller, but statistically significant, decrease in MPF was also shown for GM (0.75% change, p<0.05), and mixed WM-GM (0.75% change, p<0.01) compartments. Similar to other studies [11,12,19,30], we found anterior-to-posterior gradient in brain myelination changes. The greatest regional age-related changes in WM, reaching 5% over 7 years, were detected in the frontal cortical WM. A decrease in myelin content was less pronounced, but statistically significant, in the parietal and temporal cortical WM, and in the majority (15 of 26) subcortical WM regions. We also had shown an age-related myelination decrease in mixed WM-GM brainstem structures (except the pons), and all explored allocortical and deep GM structures. Myelination decrease was most significant in caudate putamen, 2-6% over 7 years. To our knowledge, age-related decrease in GM myelination has not been previously reported.

It is important to outline that the decrease in MPF in global WM and mixed WM-GM brain compartments significantly correlated with increased levels of MBP autoantibodies and major myelin proteins in the same subjects. Studies of autoantibodies of myelin-specific proteins in serum have mainly been carried out in the context of the hypothesis of their participation in the pathogenesis of multiple sclerosis [62,63,64,65]. The presence of antibodies to MBP was detected both in the serum of patients with MS and in the serum of healthy controls, although in lower concentrations [62,65]. For example, the study by Vojdani et al. [62], conducted with a control group of similar age (32-48 years), detected similar levels of antibodies to MBP as was shown in our current study, in contrast to its threefold increase in MS patients. Greer et al. [64] provided more compelling evidence for the involvement of anti-PLP antibodies in the pathogenesis of MS, showing the association of these antibodies with clinical scores and disease severity. In the same study, the level of anti-PLP antibodies in healthy controls was extremely low, which is similar to our results. To the best of our knowledge, no studies investigated anti-myelin plasma antibodies in normal aging.

We did not find correlations between age-related changes in myelin antibodies and changes in the volume of T2-FLAIR hyperintensities, although volume of those hyperintensity zones in WM increased significantly with age. MPF values within these lesions were substantially lower than those in the surrounding WM, thus suggesting the loss of myelin. Based on the literature, T2-FLAIR hyperintensities often correspond to the areas of demyelination but may also represent tissue damage associated with disruption of the structure of the extracellular matrix and the outflow of tissue fluids [66]. Such foci of hyperintensity frequently occur in elderly healthy people [67] and are arguably considered as risk factors for stroke and dementia [68]. Although the T2-FLAIR sequences are typically used for identifying of demyelination foci in the clinics, studies showed that hyperintensity of the T2-FLAIR signal cannot be used for unambiguous detection of myelin damage and, especially, its quantitative assessment [66]. Our results indicate that visible abnormalities of the T2-weighted signal cannot explain the global trend of age-related myelin loss, which is widespread in radiologically normal brain tissue.

The discovered interhemispheric differences in cortical WM are particularly interesting. These differences were found independently for most brain regions in different years of the study. In other words, MPF maps obtained in 2015 and 2022 were not registered to each other but were registered to the brain atlas separately.

For example, MPF in the frontal WM (rectus, precentral, middle frontal, lateral and middle fronto-orbital WM) was significantly higher in right hemisphere; however, MPF in the parietal and occipital lobe (precuneus, cuneus, supramarginal, lingual, and occipital WM) was higher for the left hemisphere. The exception is fusiform WM, where myelination is higher in the right hemisphere. Age-related changes were associated with interhemispheric differences in more regions of juxtacortical WM: in 2022 unlike 2015 interhemispheric differences were found in the lingual, occipital and hippocampal part of cingulum WM whereas only the precentral WM showed interhemispheric differences in 2015 unlike 2022.

Previously published studies of age-related changes in brain myelination are contradictory. Faizy et al. [14] did not find interhemispheric differences in the frontal, parietal, and occipital WM and CST despite the use of several measurement methods (MWF, FA, MD, RD, and MTR). Other studies that used MTR method to assess myelination, have found larger values for this parameter [59] and large age-related MTR differences in the left hemisphere [30]. MPF mapping and detailed segmentation used in the present study were able to find consistent interhemispheric differences in juxtacortical WM myelination and their dependence on age.

5. Conclusions

This is the first longitudinal study of age-related changes in brain myelination assessed with quantitative MPF mapping on the same subjects seven years apart. A significant age-related decrease in MPF was found in global WM and GM, as well as in 48 out of 82 WM and GM regions. The most pronounced MPF decrease was observed for the frontal juxtacortical WM (2-5%), genu of CC (4.0%), and caudate nucleus (5.9%). The smallest age-related changes in MPF were found in the brainstem (0.6-1.3%). The age-related decrease in MPF significantly correlated with an increase in the level of anti-MBP antibodies in the blood plasma. The volume of WM hyperintensities increased with age but did not correlate with either changes in MPF or the concentration of myelin-specific antibodies.

Our results demonstrate high sensitivity of MPF to age-related changes in brain myelination which confirms the feasibility of longitudinal studies based on the fast MPF mapping method in a clinical setting. Information about longitudinal MPF changes during normal aging provides a methodological background for future neuroimaging studies of myelin in age-related brain diseases.

6. Study limitations

The study was conducted on a small sample and a small age range of the studied subjects. Age-related decline in myelination has not been studied for subjects older than 67 years. Statistical analysis was performed within a linear model, whereas age-related changes in myelin were shown to be more consistent with a quadratic (inverted-U) model. However, a quadratic model would be less suitable for our study due to a narrow range of ages and small sample size.