Psilocybin desynchronizes the human brain.

A single dose of psilocybin, a psychedelic that acutely causes distortions of space-time perception and ego dissolution, produces rapid and persistent therapeutic effects in human clinical trials1-4. In animal models, psilocybin induces neuroplasticity in cortex and hippocampus5-8. It remains unclear how human brain network changes relate to subjective and lasting effects of psychedelics. Here we tracked individual-specific brain changes with longitudinal precision functional mapping (roughly 18 magnetic resonance imaging visits per participant). Healthy adults were tracked before, during and for 3 weeks after high-dose psilocybin (25 mg) and methylphenidate (40 mg), and brought back for an additional psilocybin dose 6-12 months later. Psilocybin massively disrupted functional connectivity (FC) in cortex and subcortex, acutely causing more than threefold greater change than methylphenidate. These FC changes were driven by brain desynchronization across spatial scales (areal, global), which dissolved network distinctions by reducing correlations within and anticorrelations between networks. Psilocybin-driven FC changes were strongest in the default mode network, which is connected to the anterior hippocampus and is thought to create our sense of space, time and self. Individual differences in FC changes were strongly linked to the subjective psychedelic experience. Performing a perceptual task reduced psilocybin-driven FC changes. Psilocybin caused persistent decrease in FC between the anterior hippocampus and default mode network, lasting for weeks. Persistent reduction of hippocampal-default mode network connectivity may represent a neuroanatomical and mechanistic correlate of the proplasticity and therapeutic effects of psychedelics.

Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

PARE: A framework for removal of confounding effects from any distance-based dimension reduction method.

Dimension reduction tools preserving similarity and graph structure such as t-SNE and UMAP can capture complex biological patterns in high-dimensional data. However, these tools typically are not designed to separate effects of interest from unwanted effects due to confounders. We introduce the partial embedding (PARE) framework, which enables removal of confounders from any distance-based dimension reduction method. We then develop partial t-SNE and partial UMAP and apply these methods to genomic and neuroimaging data. For lower-dimensional visualization, our results show that the PARE framework can remove batch effects in single-cell sequencing data as well as separate clinical and technical variability in neuroimaging measures. We demonstrate that the PARE framework extends dimension reduction methods to highlight biological patterns of interest while effectively removing confounding effects.

Imaging chronic active lesions in multiple sclerosis: a consensus statement.

Chronic active lesions (CAL) are an important manifestation of chronic inflammation in multiple sclerosis (MS) and have implications for non-relapsing biological progression. In recent years, the discovery of innovative magnetic resonance imaging (MRI) and PET derived biomarkers has made it possible to detect CAL, and to some extent quantify them, in the brain of persons with MS, in vivo. Paramagnetic rim lesions on susceptibility-sensitive MRI sequences, MRI-defined slowly expanding lesions on T1-weighted (T1-w) and T2-w scans, and 18-kDa translocator protein-positive lesions on PET are promising candidate biomarkers of CAL. While partially overlapping, these biomarkers do not have equivalent sensitivity and specificity to histopathological CAL. Standardization in the use of available imaging measures for CAL identification, quantification, and monitoring is lacking. To fast-forward clinical translation of CAL, the North American Imaging in Multiple Sclerosis Cooperative developed a Consensus Statement, which provides guidance for the radiological definition and measurement of CAL. The proposed manuscript presents this Consensus Statement, summarizes the multistep process leading to it, and identifies the remaining major gaps in knowledge.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging.

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

Penalized decomposition using residuals (PeDecURe) for feature extraction in the presence of nuisance variables.

Neuroimaging data are an increasingly important part of etiological studies of neurological and psychiatric disorders. However, mitigating the influence of nuisance variables, including confounders, remains a challenge in image analysis. In studies of Alzheimer's disease, for example, an imbalance in disease rates by age and sex may make it difficult to distinguish between structural patterns in the brain (as measured by neuroimaging scans) attributable to disease progression and those characteristic of typical human aging or sex differences. Concerningly, when not properly accounted for, nuisance variables pose threats to the generalizability and interpretability of findings from these studies. Motivated by this critical issue, in this work, we examine the impact of nuisance variables on feature extraction methods and propose Penalized Decomposition Using Residuals (PeDecURe), a new method for obtaining nuisance variable-adjusted features. PeDecURe estimates primary directions of variation which maximize covariance between partially residualized imaging features and a variable of interest (e.g., Alzheimer's diagnosis) while simultaneously mitigating the influence of nuisance variation through a penalty on the covariance between partially residualized imaging features and those variables. Using features derived using PeDecURe's first direction of variation, we train a highly accurate and generalizable predictive model, as evidenced by its robustness in testing samples with different underlying nuisance variable distributions. We compare PeDecURe to commonly used decomposition methods (principal component analysis (PCA) and partial least squares) as well as a confounder-adjusted variation of PCA. We find that features derived from PeDecURe offer greater accuracy and generalizability and lower correlations with nuisance variables compared with the other methods. While PeDecURe is primarily motivated by challenges that arise in the analysis of neuroimaging data, it is broadly applicable to data sets with highly correlated features, where novel methods to handle nuisance variables are warranted.

Similarity-based multimodal regression.

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

Effect of patient factors, center, and era on Fontan timing: An observational study using the Pediatric Health Information Systems Database.

BACKGROUND: There are no consensus guidelines defining optimal timing for the Fontan operation, the last planned surgery in staged palliation for single-ventricle heart disease. OBJECTIVES: Identify patient-level characteristics, center-level variation, and secular trends driving Fontan timing. METHODS: A retrospective observational study of subjects who underwent Fontan from 2007 to 2021 at centers in the Pediatric Health Information Systems database was performed using linear mixed-effects modeling in which age at Fontan was regressed on patient characteristics and date of operation with center as random effect. RESULTS: We included 10,305 subjects (40.4% female, 44% non-white) at 47 centers. Median age at Fontan was 3.4 years (IQR 2.6-4.4). Hypoplastic left heart syndrome (-4.4 months, 95%CI -5.5 to -3.3) and concomitant conditions (-2.6 months, 95%CI -4.1 to -1.1) were associated with younger age at Fontan. Subjects with technology-dependence (+4.6 months, 95%CI 3.1-6.1) were older at Fontan. Black (+4.1 months, 95%CI 2.5-5.7) and Asian (+8.3 months, 95%CI 5.4-11.2) race were associated with older age at Fontan. There was significant variation in Fontan timing between centers. Center accounted for 10% of variation (ICC 0.10, 95%CI 0.07-0.14). Center surgical volume was not associated with Fontan timing (P = .21). Operation year was associated with age at Fontan, with a 3.1 month increase in age for every 5 years (+0.61 months, 95%CI 0.48-0.75). CONCLUSIONS: After adjusting for patient-level characteristics there remains significant inter-center variation in Fontan timing. Age at Fontan has increased. Future studies addressing optimal Fontan timing are warranted.

SpaGCN: Integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network.

Recent advances in spatially resolved transcriptomics (SRT) technologies have enabled comprehensive characterization of gene expression patterns in the context of tissue microenvironment. To elucidate spatial gene expression variation, we present SpaGCN, a graph convolutional network approach that integrates gene expression, spatial location and histology in SRT data analysis. Through graph convolution, SpaGCN aggregates gene expression of each spot from its neighboring spots, which enables the identification of spatial domains with coherent expression and histology. The subsequent domain guided differential expression (DE) analysis then detects genes with enriched expression patterns in the identified domains. Analyzing seven SRT datasets using SpaGCN, we show it can detect genes with much more enriched spatial expression patterns than competing methods. Furthermore, genes detected by SpaGCN are transferrable and can be utilized to study spatial variation of gene expression in other datasets. SpaGCN is computationally fast, platform independent, making it a desirable tool for diverse SRT studies.

Early Magnetic Resonance Imaging Features of New Paramagnetic Rim Lesions in Multiple Sclerosis.

OBJECTIVE: To determine early magnetic resonance imaging (MRI) features of new multiple sclerosis (MS) lesions that will develop into paramagnetic rim lesions (PRLs), which have been associated with progressive tissue injury in MS. METHODS: New contrast-enhancing lesions observed on routine clinical MRI were imaged at 7 T within 4 weeks of observation, and 3 and 6 months later. The 6-month MRI was used to classify PRL status (PRL or non-PRL). The relationship between early lesion characteristics and subsequent PRL status was assessed using generalized linear mixed effects models. Random forest classification was performed to classify early predictors of subsequent PRL status. RESULTS: From 93 contrast-enhancing lesions in 23 MS patients, 37 lesions developed into a PRL. In lesions that developed into PRLs compared with those that did not, the average lesion T1 on the initial 7 T MRI was 1994 ms compared with 1,670 ms (p-value <0.001), and the average volume was 168.7 mL compared with 44 mL (p-value <0.001) in lesions that did not. These volume differences were also found on 3 T scans (p-value <0.001), and for intensity-normalized T1 -w (p-value = 0.011) and fluid-attenuated inversion recovery (p-value = 0.005). The area under the receiver operating characteristic curve for the random forest classification with leave-one-out cross-validation was found to be 0.86 using initial 7 T features. INTERPRETATION: New MS lesions that evolve into PRLs can be identified early in lesion evolution. These findings suggest that biological mechanisms underlying PRL development begin early, which has important implications for clinical trials targeting PRLs development and subsequent therapeutics. ANN NEUROL 2023;94:736-744.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium.

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Choroid plexus volume differentiates MS from its mimics.

This study aimed to determine whether choroid plexus volume (CPV) could differentiate multiple sclerosis (MS) from its mimics. A secondary analysis of two previously enrolled studies, 50 participants with MS and 64 with alternative diagnoses were included. CPV was automatically segmented from 3T magnetic resonance imaging (MRI), followed by manual review to remove misclassified tissue. Mean normalized choroid plexus volume (nCPV) to intracranial volume demonstrated relatively high specificity for MS participants in each cohort (0.80 and 0.76) with an area under the receiver-operator characteristic curve of 0.71 (95% confidence interval (CI) = 0.55-0.87) and 0.65 (95% CI = 0.52-0.77). In this preliminary study, nCPV differentiated MS from its mimics.

A multicenter pilot study evaluating simplified central vein assessment for the diagnosis of multiple sclerosis.

BACKGROUND: The central vein sign (CVS) is a proposed magnetic resonance imaging (MRI) biomarker for multiple sclerosis (MS); the optimal method for abbreviated CVS scoring is not yet established. OBJECTIVE: The aim of this study was to evaluate the performance of a simplified approach to CVS assessment in a multicenter study of patients being evaluated for suspected MS. METHODS: Adults referred for possible MS to 10 sites were recruited. A post-Gd 3D T2*-weighted MRI sequence (FLAIR*) was obtained in each subject. Trained raters at each site identified up to six CVS-positive lesions per FLAIR* scan. Diagnostic performance of CVS was evaluated for a diagnosis of MS which had been confirmed using the 2017 McDonald criteria at thresholds including three positive lesions (Select-3*) and six positive lesions (Select-6*). Inter-rater reliability assessments were performed. RESULTS: Overall, 78 participants were analyzed; 37 (47%) were diagnosed with MS, and 41 (53%) were not. The mean age of participants was 45 (range: 19-64) years, and most were female (n = 55, 71%). The area under the receiver operating characteristic curve (AUROC) for the simplified counting method was 0.83 (95% CI: 0.73-0.93). Select-3* and Select-6* had sensitivity of 81% and 65% and specificity of 68% and 98%, respectively. Inter-rater agreement was 78% for Select-3* and 83% for Select-6*. CONCLUSION: A simplified method for CVS assessment in patients referred for suspected MS demonstrated good diagnostic performance and inter-rater agreement.

Alcohol for seizure induction in the epilepsy monitoring unit.

RATIONALE: Seizure induction techniques are used in the epilepsy monitoring unit (EMU) to increase diagnostic yield and reduce length of stay. There are insufficient data on the efficacy of alcohol as an induction technique. METHODS: We performed a retrospective cohort study using six years of EMU data at our institution. We compared cases who received alcohol for seizure induction to matched controls who did not. The groups were matched on the following variables: age, reason for admission, length of stay, number of antiseizure medications (ASM) at admission, whether ASMs were tapered during admission, and presence of interictal epileptiform discharges. We used both propensity score and exact matching strategies. We compared the likelihood of epileptic seizures and nonepileptic events in cases versus controls using Kaplan-Meier time-to-event analysis, as well as odds ratios for these outcomes occurring at any time during the admission. RESULTS: We analyzed 256 cases who received alcohol (median dose 2.5 standard drinks) and 256 propensity score-matched controls. Cases who received alcohol were no more likely than controls to have an epileptic seizure (X2(1) = 0.01, p = 0.93) or nonepileptic event (X2(1) = 2.1, p = 0.14) in the first 48 h after alcohol administration. For the admission overall, cases were no more likely to have an epileptic seizure (OR 0.89, 95 % CI 0.61-1.28, p = 0.58), nonepileptic event (OR 0.97, CI 0.62-1.53, p = 1.00), nor require rescue benzodiazepine (OR 0.63, CI 0.35-1.12, p = 0.15). Stratified analyses revealed no increased risk of epileptic seizure in any subgroups. Sensitivity analysis using exact matching showed that results were robust to matching strategy. CONCLUSIONS: Alcohol was not an effective induction technique in the EMU. This finding has implications for counseling patients with epilepsy about the risks of drinking alcohol in moderation in their daily lives.

Neural Responses to Auditory Food Stimuli Following Cognitive Behavioral Therapy for Binge-Eating Disorder.

OBJECTIVE: Adults with binge-eating disorder (BED), compared with those without BED, demonstrate higher blood-oxygen-level-dependent (BOLD) response to food cues in reward-related regions of the brain. It is not known whether cognitive behavioral therapy (CBT) can reverse this reward system hyperactivation. This randomized controlled trial (RCT) assessed changes in BOLD response to binge-eating cues following CBT versus wait-list control (WLC). METHOD: Females with BED (N = 40) were randomized to CBT or WLC. Participants completed assessments at baseline and 16 weeks including measures of eating and appetite and functional magnetic resonance imaging (fMRI) to measure BOLD response while listening to personalized scripts of binge-eating and neutral-relaxing cues. Data were analyzed using general linear models with mixed effects. RESULTS: Overall retention rate was 87.5%. CBT achieved significantly greater reductions in binge-eating episodes than WLC (mean ± standard error decline of 14.6 ± 2.7 vs. 5.7 ± 2.8 episodes in the past 28 days, respectively; p = 0.03). CBT and WLC did not differ significantly in changes in neural responses to binge-eating stimuli during the fMRI sessions. Compared with WLC, CBT had significantly greater improvements in reward-based eating drive, disinhibition, and hunger as assessed by questionnaires (ps < 0.05). DISCUSSION: CBT was effective in reducing binge eating, but, contrary to our hypothesis, CBT did not improve BOLD response to auditory binge-eating stimuli in reward regions of the brain. Further studies are needed to assess mechanisms underlying improvements with CBT for BED. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT03604172.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization.

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

ModelArray: An R package for statistical analysis of fixel-wise data.

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Matrix decomposition for modeling lesion development processes in multiple sclerosis.

Our main goal is to study and quantify the evolution of multiple sclerosis lesions observed longitudinally over many years in multi-sequence structural magnetic resonance imaging (sMRI). To achieve that, we propose a class of functional models for capturing the temporal dynamics and spatial distribution of the voxel-specific intensity trajectories in all sMRI sequences. To accommodate the hierarchical data structure (observations nested within voxels, which are nested within lesions, which, in turn, are nested within study participants), we use structured functional principal component analysis. We propose and evaluate the finite sample properties of hypothesis tests of therapeutic intervention effects on lesion evolution while accounting for the multilevel structure of the data. Using this novel testing strategy, we found statistically significant differences in lesion evolution between treatment groups.