Cortical lesions uniquely predict motor disability accrual and form rarely in the absence of new white matter lesions in multiple sclerosis =========================================================================================================================================== * Erin S Beck * W Andrew Mullins * Jonadab dos Santos Silva * Stefano Filippini * Prasanna Parvathaneni * Josefina Maranzano * Mark Morrison * Daniel J Suto * Corinne Donnay * Henry Dieckhaus * Nicholas J Luciano * Kanika Sharma * María Ines Gaitán * Jiaen Liu * Jacco A de Zwart * Peter van Gelderen * Irene Cortese * Sridar Narayanan * Jeff H Duyn * Govind Nair * Pascal Sati * Daniel S Reich ## Abstract **Background and objectives** Cortical lesions (CL) are common in multiple sclerosis (MS) and associate with disability and progressive disease. We asked whether CL continue to form in people with stable white matter lesions (WML) and whether the association of CL with worsening disability relates to pre-existing or new CL. **Methods** A cohort of adults with MS were evaluated annually with 7 tesla (T) brain magnetic resonance imaging (MRI) and 3T brain and spine MRI for 2 years, and clinical assessments for 3 years. CL were identified on 7T images at each timepoint. WML and brain tissue segmentation were performed using 3T images at baseline and year 2. **Results** 59 adults with MS had ≥1 7T follow-up visit (mean follow-up time 2±0.5 years). 9 had “active” relapsing-remitting MS (RRMS), defined as new WML in the year prior to enrollment. Of the remaining 50, 33 had “stable” RRMS, 14 secondary progressive MS (SPMS), and 3 primary progressive MS. 16 total new CL formed in the active RRMS group (median 1, range 0-10), 7 in the stable RRMS group (median 0, range 0-5), and 4 in the progressive MS group (median 0, range 0-1) (p=0.006, stable RR vs PMS p=0.88). New CL were not associated with greater change in any individual disability measure or in a composite measure of disability worsening (worsening Expanded Disability Status Scale or 9-hole peg test or 25-foot timed walk). Baseline CL volume was higher in people with worsening disability (median 1010μl, range 13-9888 vs median 267μl, range 0-3539, p=0.001, adjusted for age and sex) and in individuals with RRMS who subsequently transitioned to SPMS (median 2183μl, range 270-9888 vs median 321μl, range 0-6392 in those who remained RRMS, p=0.01, adjusted for age and sex). Baseline WML volume was not associated with worsening disability or transition from RRMS to SPMS. **Discussion** CL formation is rare in people with stable WML, even in those with worsening disability. CL but not WML burden predicts future worsening of disability, suggesting that the relationship between CL and disability progression is related to long-term effects of lesions that form in the earlier stages of disease, rather than to ongoing lesion formation. ## Introduction In progressive multiple sclerosis (MS), there is gradual accrual of disability over time, often in the absence of new white matter or spinal cord lesions. This accrual of disability may be related to long-term effects of focal demyelination that occurs early in disease, persistent inflammation either diffusely or within pre-existing lesions, and/or ongoing focal demyelination that is not detected on standard MRI, including in the cerebral cortex.1, 2 Cortical MS lesions are associated with disability and progression, possibly to a greater extent than white matter lesions.3–6 Cortical lesions, and in particular subpial cortical lesions, which involve the superficial cortical layers, have been hypothesized to form via a partially distinct mechanism from other MS lesions that is related to inflammation in the overlying meninges.7–9 Despite their acknowledged prevalence on histopathologic studies, where a median of 14% of the cortex is demyelinated,10–12 cortical lesions have until recently been difficult to visualize on in vivo MRI. This has limited our understanding of when they develop and how they contribute to disability. Specifically, it is unclear whether the association between cortical lesions and progression is due to a persistent and long-term effect of cortical lesions that form early in the disease course or to ongoing cortical lesion formation in later stages of the disease. Intriguingly, recent studies have indicated that cortical lesion formation may be higher in progressive MS than in relapsing remitting MS (RRMS),13, 14 in contrast to white matter lesion formation, the vast majority of which occurs in the relapsing phase of the disease. In this study, we followed three groups longitudinally — people in whom new white matter lesions had formed in the year prior to their enrollment (active RRMS); people with RRMS with stable white matter lesions in the year prior to enrollment (stable RRMS); and people with progressive MS with stable white matter lesions in the year prior to enrollment — with the goal of determining whether cortical lesions continue to form in people with stable white matter lesions and whether the association of cortical lesions with worsening disability is related to pre-existing or new cortical lesions. We found, surprisingly, that cortical lesion formation is rare in people with stable white matter lesions, even in those with worsening disability. Importantly, cortical but not white matter lesion burden predicts future worsening of disability, suggesting that the relationship between cortical lesions and disability progression is related to long-term effects of lesions that formed in the earlier stages of disease, rather than to ongoing lesion formation. ## Methods ### Clinical cohort Participants enrolled in an institutional review board-approved MS natural history study at the National Institutes of Health provided written, informed consent. We built a prospective subcohort of MS patients who were ≥18 years old and without 7T MRI contraindication. Between May 2017 and November 2018, we prospectively enrolled individuals with stable white matter and spinal cord lesions by 3T MRI in the year prior to enrollment (stable MS), including 36 with RRMS and 19 with progressive MS (15 secondary progressive (SPMS) and four primary progressive (PPMS)). For comparison, we also enrolled nine individuals with RRMS who had at least one new or contrast-enhancing white matter lesion in the year prior to enrollment (active MS). As we did not have any data on rate of cortical lesion accumulation using our more sensitive 7T techniques at the time the study was initiated, sample size was driven by recruitment feasibility in our center, and we prioritized recruiting participants with radiographically stable MS. Participants were evaluated annually for two years, including a clinical history and physical exam, expanded disability status scale (EDSS), paced auditory symbol addition test (PASAT), symbol digit modality test (SDMT, paper based), 25-foot timed walk (25TW), and 9-hole peg test (9HPT), a 7T brain MRI, and a 3T brain and spinal cord MRI. If obtained more than one day apart, the 25TW and 9HPT were performed at both 3T and 7T visits, and results were averaged. After the first two years, participants were followed annually with the same clinical testing for up to two additional years. Clinical relapses were determined based on chart review by an MS specialist. ### MRI acquisition All 3T scans were acquired on a single Magnetom Skyra scanner (Siemens, Erlangen, Germany), equipped with a 32-channel head coil. The first 17 baseline scans were acquired prior to a software upgrade (D13 to E11). 3T brain scans included axial 2D proton density (PD)/T2w, sagittal 3D magnetization prepared 2 rapid gradient echo (MP2RAGE), and sagittal 3D T2w fluid attenuated inversion recovery (FLAIR) before and after administration of gadobutrol (unless contraindicated or refused by the participant). Several upgrades to the MP2RAGE version on the 3T scanner occurred in the first year of the study, after which there were no further changes in MP2RAGE version. Spinal cord scans included sagittal 2D short-tau inversion recovery (STIR) and sagittal 3D T1w gradient-recalled echo (GRE), sagittal 3D MP2RAGE of the cervical spine, and axial 3D MP2RAGE of the thoracic spine. All 7T brain scans were performed on a 7T whole-body research system (Siemens) equipped with a single-channel transmit, 32-channel phased array receive head coil. A software upgrade (B17 to E12) was performed prior to the final five year-2 scans. 7T scans included axial 3D MP2RAGE (0.5mm isometric, acquired four times per scan session), sagittal 3D MP2RAGE (0.7mm isometric), and sagittal 3D segmented T2*w echo-planar imaging (EPI) (0.5mm isometric, acquired in two partially overlapping volumes for full brain coverage). At the baseline and year-1 timepoints, axial 2D T2*w multi-echo GRE was also acquired. At the year-2 timepoint, for 30 scans (all scans acquired after April 2020), an axial 3D T2*w multi-echo GRE with navigator-based B and motion correction was acquired in place of the 2D T2*w multi-echo GRE, as we found that this sequence improved image quality.15, 16 Both types of T2*w GRE scans were 0.5mm isometric and acquired in three partially overlapping volumes for near full supratentorial brain coverage. MP2RAGE data were processed into uniform denoised images (hereafter T1w MP2RAGE) and T1 maps using manufacturer-provided software. 7T images did not fully cover the inferior temporal and occipital lobes. Median interval between 3T and 7T MRI was 27 days (IQR 61, range 1–294). Sequence parameters are listed in Supplementary Table 1. ### Image processing The 7T 0.5mm T1w and T1 map MP2RAGE repetitions were coregistered and median T1w and T1 map images were generated, as described previously.17 The 7T T2*w GRE magnitude images were averaged across echo times, and within each timepoint, the averaged GRE images and the T2*w EPI images were linearly registered (FMRIB’s Linear Image Registration Tool, FLIRT18, 19) to the 7T T1w MP2RAGE median images. For the year-1 and year-2 timepoints, 7T T1 MP2RAGE median images were then registered to the baseline T1w MP2RAGE (Analysis of Functional NeuroImages, AFNI, 3DAllineate20, 21), and the transformation matrix was applied to the T2*w EPI and GRE magnitude images. Subtraction images were generated between registered T1 map images at each timepoint. Baseline 3T images were processed and registered as described previously.22 For follow-up timepoints, FLAIR images were registered to T1w MP2RAGE (FLIRT), followed by registration of all FLAIR and T1w MP2RAGE to the baseline T1w MP2RAGE (AFNI, 3DAllineate). ### Longitudinal cortical and white matter lesion segmentation Cortical lesions were manually segmented on baseline images using median 7T T1w and T1 map MP2RAGE, T2*w GRE, and T2*w EPI images by two independent raters as previously described.23 Lesions were categorized as leukocortical (type 1, involving cortex and white matter), intracortical (type 2, confined to the cortex and not touching the pial surface of the brain), and subpial (type 3, involving the cortex exclusively and touching the pial surface, and type 4, involving cortex and the white matter and touching the pial surface).10 Cortical lesions were hypointense on T1w MP2RAGE images and/or hyperintense on T2*w images and were seen on at least two consecutive axial slices. Cortical lesion volumes and median T1 within lesions, including the white matter portion of leukocortical lesions, were calculated. For each subsequent timepoint, MP2RAGE and T2*w images were evaluated by two independent raters for the presence of new cortical lesions compared to the baseline scan. Cortical lesions that were visible at year 1 or 2 and whose presence at earlier timepoints was potentially obscured by motion were not segmented. In addition, the mask for each cortical lesion present at baseline was evaluated and adjusted at each timepoint both for imperfect registration as well as for changes in lesion size. White matter lesions were segmented on baseline 3T images as previously described.23 Baseline white matter lesion masks were adjusted using registered year-2 T1w MP2RAGE and FLAIR images, and new lesions were identified and segmented. Longitudinal 7T MP2RAGE images were also assessed for new white matter lesions. Occasionally, small, new white matter lesions were identified only on 7T images but were seen in retrospect on 3T images. Once identified on 7T images, these new lesions were segmented on the corresponding 3T images. ### Spinal cord lesion assessment Spinal cord lesions were identified by two independent raters followed by consensus review using sagittal STIR and MP2RAGE images of the cervical spinal cord and sagittal STIR and axial MP2RAGE images of the thoracic spinal cord at baseline and year 2. ### Brain volume and brain atrophy measurement Brain volume measurements were obtained from 0.7mm T1w MP2RAGE images at all available 7T timepoints using Pseudo-Label Assisted “no-new-Net” version of U-Net algorithm (PLAn).24 PLAn is a transfer-learning method that uses pseudo-labels obtained from C-DEF segmentation of 3T images25 which were then fine tuned on MP2RAGE images from 7T in a separate cohort of patients. Retraining of PLAn was not required herein as the segmentation model trained on images from the scanner already existed. PLAn segmentations were assessed for quality and adjusted manually as needed. 0.7mm MP2RAGE was not available for all participants at all timepoints. For 47 participants, 0.7mm MP2RAGE was acquired at baseline and at least one follow-up timepoint. 10 scans from 5 individuals were discarded due to poor segmentation, leaving a total of 42 individuals with multiple timepoints available to analyze. Lesions identified by PLAn were filled and cortex, white matter, and deep gray matter volumes were extracted. Because of inconsistent infratentorial coverage in our 7T images, the brainstem and cerebellum were not included in the brain volume measurements. Baseline brain volumes were normalized using the MNI152 intracranial volume ([https://nist.mni.mcgill.ca/icbm-152lin/](https://nist.mni.mcgill.ca/icbm-152lin/)). To determine the rate of atrophy, annual atrophy rates were calculated from raw brain volumes at each timepoint (baseline to year 1 and year 1 to year 2), which were then averaged across timepoints for each individual. ### Statistical analyses Comparisons between stable RRMS and PMS and between people with stable vs new lesions were performed using t-tests, Fisher’s exact test, and Mann-Whitney tests as appropriate. Comparisons between active RRMS, stable RRMS, and PMS were performed using Kruskal-Wallis tests and one-way ANOVA as appropriate. To investigate the relationship between baseline MRI measures and subsequent disability change and brain atrophy, we first performed Box-Cox transformation of changes in disability measures and lesion counts and volumes, and then calculated partial correlation coefficients, adjusting for age and sex, with 5000 nonparametric bootstrapping iterations to increase analysis robustness.26, 27 Comparisons between people with progression of disability, defined as an increase in EDSS of ≥1 for baseline EDSS <6 or ≥0.5 for baseline EDSS ≥6 or 20% increase in 25TW or 20% increase in 9HPT, vs stable disability and between people who transitioned from RRMS to SPMS vs those who remained RRMS were performed using a multivariable generalized linear model adjusted for age and sex. Box-Cox transformation of lesion counts and volumes was performed prior to running these models. Adjustment for multiple comparisons was performed using the Benjamini and Hochberg procedure.28 Statistical analyses were performed on IBM SPSS Statistics, version 26 (IBM Corp., Armonk, NY, USA). ## Results ### Cortical lesion formation is rare in people with stable white matter lesions and is not associated with worsening disability 64 individuals with MS, including 9 active RRMS, 36 stable RRMS, and 19 progressive MS (15 SPMS, 4 PPMS), underwent baseline clinical evaluation, 3T brain and spinal cord MRI, and 7T brain MRI. The results of the baseline analyses have been previously described.23 59/64 participants (92%) returned for at least one 7T follow-up visit (9 active RRMS, 33 stable RRMS, and 17 progressive MS [14 SPMS, 3 PPMS]) (Table 1). Mean ± standard deviation (SD) total 7T follow-up was 2 ± 0.5 years and did not differ between groups (p=0.77). View this table: [Table 1.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/T1) Table 1. Baseline cohort characteristics 12 individuals formed at least one new cortical lesion, including 5/9 (56%) active RRMS, 3/33 (9%) stable RRMS, and 4/17 (24%) progressive MS (p=0.008, stable RR vs PMS p=0.21). A total of 16 new cortical lesions formed in the active RRMS group (median 1 per participant, IQR 2, range 0-10), 7 in the stable RRMS group (median 0, IQR 0, range 0-5), and 4 in the progressive MS group (median 0, IQR 1, range 0-1) (p=0.006, stable RR vs PMS p=0.88) (Figure 1A, Table 2). There was no difference in the subtypes of new cortical lesions (leukocortical, intracortical, and subpial) between the three groups (p=0.87). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/25/2023.09.22.23295974/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/F1) Figure 1. New cortical lesion formation is rare in people with stable white matter lesions. (A) New cortical and white matter lesion formation was frequent in people with active RRMS (defined as having at least one new or enhancing white matter lesion in the year prior to enrollment) but was infrequent in people with radiographically stable RRMS or progressive MS. New cortical lesion formation was not related to baseline cortical lesion number (gray bars, each bar represents a single participant). (B-D) Appearance of new cortical lesions can differ from the appearance of chronic cortical lesions (cortical lesions are denoted with arrows). (B) A subpial lesion that was new on the year-2 MRI appears hypointense on T2*w GRE. (C) A leukocortical lesion that was present at baseline is initially hyperintense on T2*w GRE with a band of hypointensity at the cortex-white matter junction. This lesion expands over time and becomes more hypointense on T1w MP2RAGE, while at the same time the hypointense band on T2*w GRE disappears. (D) A subpial lesion that was new on the year-1 MRI initially appears isointense on T2*w GRE and then becomes more hyperintense at year-2. (E) A white matter lesion in the same individual as (C) was new on the year-2 MRI and is hyperintense on T2*w GRE. RRMS: relapsing remitting multiple sclerosis. SPMS: secondary progressive multiple sclerosis. PPMS: primary progressive multiple sclerosis. T2*w GRE: T2* weighted gradient recalled echo. MP2RAGE: magnetization prepared 2 rapid acquisition gradient echo. View this table: [Table 2.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/T2) Table 2. Longitudinal MRI changes by MS phenotype. Individuals with new cortical lesions did not differ from those with no new cortical lesions with respect to age, sex, disease duration, or baseline cortical lesion burden (Table 3). Disease modifying therapy (categorized as none, injectable, oral, infusion, or changed during the study) did not differ between people with new vs no new cortical lesions. 3/9 people on infusion therapies formed one new cortical lesion each. Two people were on ocrelizumab: one received their first dose approximately 1 year before the baseline MRI and the other received their first dose 8 months before the baseline MRI; neither had any interruptions in treatment. One person was on natalizumab for over 1 year before the baseline MRI and did not have any known interruptions in treatment. There was no difference in change in any of the disability measures or in overall disability progression (defined as an increase in EDSS of ≥1 for baseline EDSS <6 or ≥0.5 for baseline EDSS ≥6 or 20% increase in 25TW or 20% increase in 9HPT) in people with new vs no new cortical lesions (mean clinical follow-up time 2.8 ± 0.7 years). View this table: [Table 3.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/T3) Table 3. Characteristics of individuals with new vs stable lesions. 19 individuals formed a total of 96 new white matter lesion, including 9/9 (100%) active RRMS, 8/33 (24%) stable RRMS, and 2/17 (12%) progressive MS (p<0.0001, stable RR vs PMS p=0.46) (Figure 1A, Table 2). There were no differences in age, sex, disease duration, baseline white matter lesion volume, disease modifying therapy, or change in disability between people with new vs no new white matter lesions. 1/9 people on infusion therapies (ocrelizumab) formed a new white matter lesion; this individual’s first dose of ocrelizumab was two months before their baseline MRI. There was no difference in change in any of the disability measures or in overall disability progression in people with new vs no new white matter lesions (Table 3). 6/12 people (50%) with new cortical lesions also formed new white matter lesions, while 13/47 (27%) of people without new cortical lesions formed new white matter lesions (p=0.17). There was no difference in the subtypes of cortical lesions that formed in people with vs without new white matter lesions (p=0.69). Three white matter lesions in three individuals expanded into the cortex during the 7T follow up period. No intracortical or subpial lesions expanded into the white matter. New spinal cord lesions were rare (10 total new lesions in 4 individuals) and did not differ between groups, nor did change in brain volume (Table 2). Four individuals had a total of six relapses during the follow-up period. Two individuals were in the active RRMS group, and both had new lesions in the white matter, spinal cord, and cortex. Two individuals with clinical relapses were in the stable RRMS group; neither of these individuals had new lesions in the brain or spinal cord. None of the four had overall disability progression (composite measure of EDSS, 9HPT, and 25TW), and none transitioned from RRMS to SPMS. ### New cortical lesions have a distinct appearance on T2*w images While most cortical lesions are hypointense on T1w images and hyperintense on T2*w images, we observed that 35% of the new cortical lesions, while hypointense on T1w images, were hypointense or isointense on T2*w images. In contrast, only 10% of new white matter lesions were iso- or hypointense on T2*w images (p=0.01) (Figure 1B-E). ### Baseline cortical lesion burden, but not white matter or spinal cord lesion burden, is associated with worsening disability and transition from relapsing to progressive MS To determine the association between baseline cortical lesion burden and subsequent change in disability, we considered all available timepoints at which there was clinical follow-up. Mean ± SD clinical follow up time was 2.8 ± 0.7 years. There was no difference in clinical follow-up time between active RRMS, stable RRMS, and progressive MS (p=0.07). Baseline total cortical lesion volume was correlated with percent change in 9HPT (r=0.69, p<0.001, adjusted for age and sex, significant after adjusting for multiple comparisons) and 25FTW (r=0.49, p<0.001, significant after adjusting for multiple comparisons). Baseline leukocortical, intracortical, and subpial lesion volume were each associated with percent change in 9HPT and percent change in 25FTW. Baseline cortical lesion volume was not correlated with change in EDSS, PASAT, or SDMT. Baseline white matter lesion volume was correlated with percent change in 25FTW (r=0.45, p<0.001) and PASAT (r=-0.38, p<0.01). Baseline spinal cord lesion number and normalized brain volume were not correlated with change in any of the individual disability measures (Table 4). View this table: [Table 4.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/T4) Table 4. Association between baseline MRI measures and longitudinal change in disability. We also determined whether baseline lesion burden was associated with overall progression of disability, defined as an increase in EDSS of ≥1 for baseline EDSS <6 or ≥0.5 for baseline EDSS ≥6 or 20% increase in 25TW or 20% increase in 9HPT. Baseline cortical lesion volume was higher in people with subsequent progression of disability (median 1010 ml, IQR 3162, range 13-9888) than in those without (median 267 ml, IQR 854, range 0-3539, B 0.809, p=0.001, adjusted for age and sex) (Figure 2A, table 5). Baseline volume was also higher in people with progression disability for each of the cortical lesion subtypes. In contrast, baseline white matter lesion volume did not differ between people with subsequent progression of disability (median 10.4 ml, IQR 14.2, range 1.5-63.1) vs those without (median 6.5 ml, IQR 12.0, range 0.9-50.8, B 0.206, p=0.42). There was also no difference in baseline spinal cord lesion number or baseline normalized brain volume between people with progression of disability vs those without (Figure 2A, table 5). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/25/2023.09.22.23295974/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/F2) Figure 2. People with worsening disability and transition from relapsing to progressive MS have higher cortical lesion burden at baseline. (A) Individuals with progression of disability during the study (“Prog,” defined as an increase in EDSS of ≥1 for baseline EDSS <6 or ≥0.5 for baseline EDSS ≥6 or 20% increase in 25TW or 20% increase in 9HPT) had higher baseline cortical lesion volume and subpial lesion volume, but no difference in baseline white matter lesion volume, spinal cord lesion number, or normalized brain volume, compared to people without progression of disability (“No prog”). Similar results were observed when comparing individuals who transitioned from relapsing remitting to secondary progressive MS with those who remained relapsing remitting (B). Comparisons were performed using a multivariable generalized linear model adjusted for age and sex. RR: relapsing remitting. SP: secondary progressive. * p<0.05, ** p<0.01 View this table: [Table 5.](http://medrxiv.org/content/early/2023/09/25/2023.09.22.23295974/T5) Table 5. Baseline MRI measures in people with progression of disability vs stable disability. Increase in cortical lesion volume was higher in people with progression of disability (B 0.831, p=0.007, adjusted for age and sex), however this was not the case for percent change in cortical lesion volume, nor was change in cortical lesion volume associated with progression of disability when baseline cortical lesion volume was included in the model (B -0.245, p=0.56). Change in white matter lesion volume and rate of brain atrophy did not differ between people with progression of disability vs those without (Supplementary table 2). 6/42 individuals with RRMS at baseline had transitioned to SPMS by the time of the final clinical follow-up, as determined clinically by the evaluating neurologist. At baseline, individuals who subsequently transitioned to SPMS had higher total cortical lesion volume (median 2183 μl, IQR 3673, range 270-9888 vs median 321 μl, IQR 832, range 0-6392, B 0.873, p=0.01, adjusted for age and sex), and subpial lesion volume (median 1489 μl, IQR 3850, range 143-9660 vs median 119 μl, IQR 372, range 0-5233, B=1.047, p=0.003) than those who remained RRMS. There was no difference in leukocortical lesion volume, white matter lesion volume, or spinal cord lesion number between those who transitioned from RRMS to SPMS vs those who remained relapsing remitting, nor did baseline brain volume differ between the groups (Figure 2B, Supplementary Table 3). ### Association between baseline lesion burden and rate of subsequent brain atrophy Baseline white matter lesion volume was associated with subsequent brain atrophy rate (r=0.320, p=0.04, adjusted for age and sex), as was leukocortical lesion number (r=0.375, p=0.02). There was no association between baseline cortical lesion number or volume, subpial or intracortical lesion number or volume, or spinal cord lesion number and subsequent change in brain volumes (Supplementary Table 4). ## Discussion We characterized cortical lesions longitudinally in a cohort of individuals with MS, most of whom had not formed new white matter lesions in the year prior to enrollment in the study. We found that cortical lesion formation was rare in this population, was not higher in progressive vs relapsing MS, and was not associated with worsening of disability over time. Importantly, however, baseline cortical lesion burden, but not white matter or spinal cord lesion burden, was higher in people who had subsequent worsening of motor disability and in people who transitioned from relapsing to secondary progressive MS. Cortical lesion formation in this study was lower than what has previously been reported,29, 30 especially in the participants who had stable white matter lesions prior to enrollment. Also in contrast to prior studies,13, 14 we did not find a higher rate of cortical lesion formation in progressive MS. These previous studies were published between 2009 and 2019, and participants were on lower efficacy treatments than the participants here, of whom 37/59 (63%) were on oral or infusion disease modifying therapies. Limited observation has previously demonstrated lower rates of cortical lesion formation in people on disease-modifying therapy.31, 32 However, these studies were done using 1.5T double inversion recovery MRI, which has poor sensitivity for cortical, and in particular subpial, lesions.33 Dedicated prospective studies will be needed to test the effects of modern disease-modifying therapies on cortical lesion formation. It is possible that despite evident differences in mechanisms of cortical and white matter lesion formation within the central nervous system, peripheral immune mediators targeted by disease-modifying therapies may be shared. The lower rate of cortical lesion formation in this cohort could also be due to the purposeful recruitment of individuals with stable white matter lesions, resulting in a cohort that was older (mean 49 years) than in prior studies and for the most part had longstanding disease (mean 14 years). We did not find differences in baseline demographic, clinical, or MRI measures in individuals with new vs no new cortical lesions, however this may be due to the very small number of individuals with new cortical lesions. Our finding that cortical lesion burden predicts disability progression and transition from relapsing to secondary progressive MS is in agreement with prior studies, including at least one prior 7T study.14, 30, 34, 35 The lack of new cortical lesions in this older cohort with longstanding disease, combined with the association between baseline cortical lesion burden and 2-year disability worsening, suggests that, at least in the current era, cortical lesions form early in disease12, 36 and then may exert long-term effects on disability. Cortical lesions could potentially exert these long term effects on disability by causing local and/or remote neurodegeneration and cell death, reflected as decreased brain volume on MRI. In both cross-sectional and longitudinal studies, the relationship between cortical lesion burden and whole brain and cortical atrophy has been inconsistent.6, 23, 34, 37–39 Here, we did not find a consistent relationship between baseline cortical lesion burden and subsequent brain atrophy. It is possible that there are more local effects on brain structure or microstructural changes that are not visible with the MRI methods used here. One potential limitation of this study is that there were changes to scanner software versions on both the 3T and 7T scanners during the study as well as a switch partway through year 2 to a motion and B-corrected version of the T2*w GRE sequence. Because of significant differences in image contrast on 3T MP2RAGE images acquired with different WIP sequences, which we found affected tissue segmentation and brain volume calculations, we used 7T MP2RAGE images for brain volume calculations. This meant that we were unable to include the infratentorial brain in our atrophy measurements due to inconsistent coverage. In addition, the lack of a longitudinal brain segmentation methods that works well with 7T images meant that we instead used a cross-sectional method at each timepoint, potentially decreasing our power to detect changes in brain volume over time. For longitudinal cortical lesion assessment, we were careful not to count as “new” any cortical lesion that was clearly seen on the motion corrected T2*w GRE images but could have been obscured by artifact on images from the earlier timepoints. This occurred rarely, and we did not observe a difference in new cortical lesion count between timepoints with vs without motion correction, but it is possible that there was a slight undercounting of new cortical lesions. This study is also limited by cohort size, which, although larger than most 7T studies in MS to date, limits our ability to identify factors associated with cortical lesion formation. This is an important question for future studies, as cortical and white matter lesion burden are not well correlated and the factors determining which patients develop cortical vs white matter lesions are unclear. If confirmed, our results have important implications for clinical care and research. They suggest that it is necessary to stop cortical lesion formation early in disease and provide some support for the possibility that existing disease-modifying therapies may offer at least partial control. The association we find between cortical lesion burden and subsequent worsening of disability may be useful for prognostication and for selecting participants for clinical trials in progressive MS who are likely to worsen over a relatively short follow up and might benefit most from treatment. ## Supporting information Supplemental tables [[supplements/295974_file05.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Funding This research was supported by the Intramural Research Program of the National Institute of Neurological Disorders and Stroke, National Institutes of Health. Erin Beck was supported by a Clinician Scientist Development Award and a Career Transition Fellowship Award from the National Multiple Sclerosis Society. ## Acknowledgments We acknowledge the staff of the National Institute of Neurological Disorders and Stroke Neuroimmunology Clinic for care of and collection of clinical data from the study participants and the staff of the National Institute of Mental Health Functional MRI Facility, Matthew Schindler, Martina Absinta, Hadar Kolb, and Maxime Donadieu for help with MRI acquisitions. We thank Siemens for providing access to research pulse sequences. We thank Tianxia Wu for her advice on statistical analyses. ## Footnotes * In the original version of the manuscript, one of the author's names was misspelled (Henry Dieckhaus). * Received September 22, 2023. * Revision received September 24, 2023. * Accepted September 25, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license ## References 1. 1.Kuhlmann T, Moccia M, Coetzee T, et al. Multiple sclerosis progression: time for a new mechanism-driven framework. Lancet Neurol. Jan 2023;22(1):78–88. doi:10.1016/S1474-4422(22)00289-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(22)00289-7&link_type=DOI) 2. 2.Reich DS, Lucchinetti CF, Calabresi PA. Multiple Sclerosis. N Engl J Med. Jan 11 2018;378(2):169–180. doi:10.1056/NEJMra1401483 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra1401483&link_type=DOI) 3. 3.Calabrese M, Agosta F, Rinaldi F, et al. Cortical lesions and atrophy associated with cognitive impairment in relapsing-remitting multiple sclerosis. Arch Neurol. Sep 2009;66(9):1144–50. doi:10.1001/archneurol.2009.174 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneurol.2009.174&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19752305&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269765000014&link_type=ISI) 4. 4.Roosendaal SD, Moraal B, Pouwels PJ, et al. Accumulation of cortical lesions in MS: relation with cognitive impairment. Mult Scler. Jun 2009;15(6):708–14. doi:10.1177/1352458509102907 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458509102907&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19435749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267186000008&link_type=ISI) 5. 5.Nielsen AS, Kinkel RP, Madigan N, Tinelli E, Benner T, Mainero C. Contribution of cortical lesion subtypes at 7T MRI to physical and cognitive performance in MS. Neurology. Aug 2013;81(7):641–9. doi:10.1212/WNL.0b013e3182a08ce8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3182a08ce8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23864311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 6. 6.Harrison DM, Roy S, Oh J, et al. Association of Cortical Lesion Burden on 7-T Magnetic Resonance Imaging With Cognition and Disability in Multiple Sclerosis. JAMA Neurol. Sep 2015;72(9):1004–12. doi:10.1001/jamaneurol.2015.1241 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamaneurol.2015.1241&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26192316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 7. 7.Howell OW, Reeves CA, Nicholas R, et al. Meningeal inflammation is widespread and linked to cortical pathology in multiple sclerosis. Brain. Sep 2011;134(Pt 9):2755–71. doi:10.1093/brain/awr182 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awr182&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21840891&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294959800028&link_type=ISI) 8. 8.Magliozzi R, Howell OW, Nicholas R, et al. Inflammatory intrathecal profiles and cortical damage in multiple sclerosis. Ann Neurol. Apr 2018;83(4):739–755. doi:10.1002/ana.25197 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25197&link_type=DOI) 9. 9.Magliozzi R, Howell OW, Reeves C, et al. A Gradient of neuronal loss and meningeal inflammation in multiple sclerosis. Ann Neurol. Oct 2010;68(4):477–93. doi:10.1002/ana.22230 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.22230&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20976767&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283398600009&link_type=ISI) 10. 10.Bo L, Vedeler CA, Nyland HI, Trapp BD, Mork SJ. Subpial demyelination in the cerebral cortex of multiple sclerosis patients. J Neuropathol Exp Neurol. Jul 2003;62(7):723–32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnen/62.7.723&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12901699&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184177300003&link_type=ISI) 11. 11.Brownell B, Hughes JT. The distribution of plaques in the cerebrum in multiple sclerosis. J Neurol Neurosurg Psychiatry. Nov 1962;25:315–20. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJqbm5wIjtzOjU6InJlc2lkIjtzOjg6IjI1LzQvMzE1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMjUvMjAyMy4wOS4yMi4yMzI5NTk3NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. 12.Lucchinetti CF, Popescu BF, Bunyan RF, et al. Inflammatory cortical demyelination in early multiple sclerosis. N Engl J Med. Dec 08 2011;365(23):2188–97. doi:10.1056/NEJMoa1100648 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1100648&link_type=DOI) 13. 13.Sethi V, Yousry T, Muhlert N, et al. A longitudinal study of cortical grey matter lesion subtypes in relapse-onset multiple sclerosis. J Neurol Neurosurg Psychiatry. Jul 2016;87(7):750–3. doi:10.1136/jnnp-2015-311102 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI4Ny83Lzc1MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzI1LzIwMjMuMDkuMjIuMjMyOTU5NzQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 14. 14.Treaba CA, Granberg TE, Sormani MP, et al. Longitudinal Characterization of Cortical Lesion Development and Evolution in Multiple Sclerosis with 7.0-T MRI. Radiology. 06 2019;291(3):740–749. doi:10.1148/radiol.2019181719 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1148/radiol.2019181719&link_type=DOI) 15. 15.Liu J, Beck ES, Filippini S, et al. Navigator-Guided Motion and B0 Correction of T2*-Weighted Magnetic Resonance Imaging Improves Multiple Sclerosis Cortical Lesion Detection. Invest Radiol. Jul 1 2021;56(7):409–416. doi:10.1097/RLI.0000000000000754 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/RLI.0000000000000754&link_type=DOI) 16. 16.Liu J, van Gelderen P, de Zwart JA, Duyn JH. Reducing motion sensitivity in 3D high-resolution T2*-weighted MRI by navigator-based motion and nonlinear magnetic field correction. Neuroimage. Feb 1 2020;206:116332. doi:10.1016/j.neuroimage.2019.116332 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2019.116332&link_type=DOI) 17. 17.Beck ES, Sati P, Sethi V, et al. Improved Visualization of Cortical Lesions in Multiple Sclerosis Using 7T MP2RAGE. AJNR Am J Neuroradiol. Feb 8 2018;doi:10.3174/ajnr.A5534 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo4OiIzOS8zLzQ1OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzI1LzIwMjMuMDkuMjIuMjMyOTU5NzQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 18. 18.Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal. Jun 2001;5(2):143–56. doi:10.1016/s1361-8415(01)00036-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1361-8415(01)00036-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11516708&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169672400005&link_type=ISI) 19. 19.Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage. Oct 2002;17(2):825–41. doi:10.1016/s1053-8119(02)91132-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1053-8119(02)91132-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12377157&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178642000027&link_type=ISI) 20. 20.Cox RW. AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res. Jun 1996;29(3):162–73. doi:10.1006/cbmr.1996.0014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/cbmr.1996.0014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8812068&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UV56700002&link_type=ISI) 21. 21.Cox RW, Hyde JS. Software tools for analysis and visualization of fMRI data. NMR Biomed. Jun-Aug 1997;10(4-5):171–8. doi:10.1002/(sici)1099-1492(199706/08)10:4/5<171::aid-nbm453>3.0.co;2-l [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/(sici)1099-1492(199706/08)10:4/5<171::aid-nbm453>3.0.co;2-l&link_type=DOI) 22. 22.Maranzano J, Dadar M, Rudko DA, et al. Comparison of Multiple Sclerosis Cortical Lesion Types Detected by Multicontrast 3T and 7T MRI. AJNR Am J Neuroradiol. Jul 2019;40(7):1162–1169. doi:10.3174/ajnr.A6099 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo5OiI0MC83LzExNjIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8yNS8yMDIzLjA5LjIyLjIzMjk1OTc0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 23. 23.Beck ES, Maranzano J, Luciano NJ, et al. Cortical lesion hotspots and association of subpial lesions with disability in multiple sclerosis. Mult Scler. Aug 2022;28(9):1351–1363. doi:10.1177/13524585211069167 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/13524585211069167&link_type=DOI) 24. 24.Diekhaus H, Donnay C, Gaitan MI, et al. Pseudo-Label Assisted Nnu-Net (PLAn) Enables Automatic Segmentation of 7T MRI From a Single Acquisition. medRxiv. 2022;2022.12.22.22283866 doi:10.1101/2022.12.22.22283866 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMi4xMi4yMi4yMjI4Mzg2NnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMjUvMjAyMy4wOS4yMi4yMzI5NTk3NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 25. 25.Selvaganesan K, Whitehead E, DeAlwis PM, et al. Robust, atlas-free, automatic segmentation of brain MRI in health and disease. Heliyon. Feb 2019;5(2):e01226. doi:10.1016/j.heliyon.2019.e01226 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.heliyon.2019.e01226&link_type=DOI) 26. 26.Chernick MR. Bootstrap methods: A guide for practicioners and researchers. 2 ed. Wiley-Interscience; 2008. 27. 27.Efron B. Bootstrap Methods: Another Look at the Jackknife. Annals of Statistics. 1979;7(1):1–26. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/aos/1176344552&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GL00300001&link_type=ISI) 28. 28.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society, Series B (Methodological). 1995;57(1):289–300. [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 29. 29.Calabrese M, Rocca MA, Atzori M, et al. Cortical lesions in primary progressive multiple sclerosis: a 2-year longitudinal MR study. Neurology. Apr 14 2009;72(15):1330–6. doi:10.1212/WNL.0b013e3181a0fee5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181a0fee5&link_type=DOI) 30. 30.Calabrese M, Rocca MA, Atzori M, et al. A 3-year magnetic resonance imaging study of cortical lesions in relapse-onset multiple sclerosis. Ann Neurol. Mar 2010;67(3):376–83. doi:10.1002/ana.21906 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.21906&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20373349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276216800013&link_type=ISI) 31. 31.Rinaldi F, Calabrese M, Seppi D, Puthenparampil M, Perini P, Gallo P. Natalizumab strongly suppresses cortical pathology in relapsing-remitting multiple sclerosis. Mult Scler. Dec 2012;18(12):1760–7. doi:10.1177/1352458512447704 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458512447704&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22570359&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 32. 32.Rinaldi F, Perini P, Atzori M, Favaretto A, Seppi D, Gallo P. Disease-modifying drugs reduce cortical lesion accumulation and atrophy progression in relapsing-remitting multiple sclerosis: results from a 48-month extension study. Mult Scler Int. 2015;2015:369348. doi:10.1155/2015/369348 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2015/369348&link_type=DOI) 33. 33.Kilsdonk ID, Jonkman LE, Klaver R, et al. Increased cortical grey matter lesion detection in multiple sclerosis with 7 T MRI: a post-mortem verification study. Brain. May 2016;139(Pt 5):1472–81. doi:10.1093/brain/aww037 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/aww037&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26956422&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 34. 34.Calabrese M, Poretto V, Favaretto A, et al. Cortical lesion load associates with progression of disability in multiple sclerosis. Brain. Oct 2012;135(Pt 10):2952–61. doi:10.1093/brain/aws246 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/aws246&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23065788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310156700006&link_type=ISI) 35. 35.Calabrese M, Romualdi C, Poretto V, et al. The changing clinical course of multiple sclerosis: a matter of gray matter. Ann Neurol. Jul 2013;74(1):76–83. doi:10.1002/ana.23882 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.23882&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23494723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 36. 36.Maranzano J, Till C, Assemlal HE, et al. Detection and clinical correlation of leukocortical lesions in pediatric-onset multiple sclerosis on multi-contrast MRI. Mult Scler. Jun 2019;25(7):980–986. doi:10.1177/1352458518779952 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458518779952&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29852831&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 37. 37.Treaba CA, Herranz E, Barletta VT, et al. The relevance of multiple sclerosis cortical lesions on cortical thinning and their clinical impact as assessed by 7.0-T MRI. J Neurol. Jul 2021;268(7):2473–2481. doi:10.1007/s00415-021-10400-4 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00415-021-10400-4&link_type=DOI) 38. 38.Nielsen AS, Kinkel RP, Tinelli E, Benner T, Cohen-Adad J, Mainero C. Focal cortical lesion detection in multiple sclerosis: 3 Tesla DIR versus 7 Tesla FLASH-T2. J Magn Reson Imaging. Mar 2012;35(3):537–42. doi:10.1002/jmri.22847 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmri.22847&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22045554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F25%2F2023.09.22.23295974.atom) 39. 39.Pareto D, Sastre-Garriga J, Auger C, et al. Juxtacortical Lesions and Cortical Thinning in Multiple Sclerosis. AJNR Am J Neuroradiol. Dec 2015;36(12):2270–6. doi:10.3174/ajnr.A4485 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czoxMDoiMzYvMTIvMjI3MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzI1LzIwMjMuMDkuMjIuMjMyOTU5NzQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)