Differential frequency in imaging-based outcome measurement: Bias in real-world oncology comparative-effectiveness studies.

BACKGROUND: Comparative-effectiveness studies using real-world data (RWD) can be susceptible to surveillance bias. In solid tumor oncology studies, analyses of endpoints such as progression-free survival (PFS) are based on progression events detected by imaging assessments. This study aimed to evaluate the potential bias introduced by differential imaging assessment frequency when using electronic health record (EHR)-derived data to investigate the comparative effectiveness of cancer therapies. METHODS: Using a nationwide de-identified EHR-derived database, we first analyzed imaging assessment frequency patterns in patients diagnosed with advanced non-small cell lung cancer (aNSCLC). We used those RWD inputs to develop a discrete event simulation model of two treatments where disease progression was the outcome and PFS was the endpoint. Using this model, we induced bias with differential imaging assessment timing and quantified its effect on observed versus true treatment effectiveness. We assessed percent bias in the estimated hazard ratio (HR). RESULTS: The frequency of assessments differed by cancer treatment types. In simulated comparative-effectiveness studies, PFS HRs estimated using real-world imaging assessment frequencies differed from the true HR by less than 10% in all scenarios (range: 0.4% to -9.6%). The greatest risk of biased effect estimates was found comparing treatments with widely different imaging frequencies, most exaggerated in disease settings where time to progression is very short. CONCLUSIONS: This study provided evidence that the frequency of imaging assessments to detect disease progression can differ by treatment type in real-world patients with cancer and may induce some bias in comparative-effectiveness studies in some situations.

Functional connectome reorganization relates to post-stroke motor recovery and structural and functional disconnection.

Motor recovery following ischemic stroke is contingent on the ability of surviving brain networks to compensate for damaged tissue. In rodent models, sensory and motor cortical representations have been shown to remap onto intact tissue around the lesion site, but remapping to more distal sites (e.g. in the contralesional hemisphere) has also been observed. Resting state functional connectivity (FC) analysis has been employed to study compensatory network adaptations in humans, but mechanisms and time course of motor recovery are not well understood. Here, we examine longitudinal FC in 23 first-episode ischemic pontine stroke patients and utilize a graph matching approach to identify patterns of functional connectivity reorganization during recovery. We quantified functional reorganization between several intervals ranging from 1 week to 6 months following stroke, and demonstrated that the areas that undergo functional reorganization most frequently are in cerebellar/subcortical networks. Brain regions with more structural and functional connectome disruption due to the stroke also had more remapping over time. Finally, we show that functional reorganization is correlated with the extent of motor recovery in the early to late subacute phases, and furthermore, individuals with greater baseline motor impairment demonstrate more extensive early subacute functional reorganization (from one to two weeks post-stroke) and this reorganization correlates with better motor recovery at 6 months. Taken together, these results suggest that our graph matching approach can quantify recovery-relevant, whole-brain functional connectivity network reorganization after stroke.

Estimation of Multiple Sclerosis lesion age on magnetic resonance imaging.

We introduce the first-ever statistical framework for estimating the age of Multiple Sclerosis (MS) lesions from magnetic resonance imaging (MRI). Estimating lesion age is an important step when studying the longitudinal behavior of MS lesions and can be used in applications such as studying the temporal dynamics of chronic active MS lesions. Our lesion age estimation models use first order radiomic features over a lesion derived from conventional T1 (T1w) and T2 weighted (T2w) and fluid attenuated inversion recovery (FLAIR), T1w with gadolinium contrast (T1w+c), and Quantitative Susceptibility Mapping (QSM) MRI sequences as well as demographic information. For this analysis, we have a total of 32 patients with 53 new lesions observed at 244 time points. A one or two step random forest model for lesion age is fit on a training set using a lesion volume cutoff of 15 mm3 or 50 mm3. We explore the performance of nine different modeling scenarios that included various combinations of the MRI sequences and demographic information and a one or two step random forest models, as well as simpler models that only uses the mean radiomic feature from each MRI sequence. The best performing model on a validation set is a model that uses a two-step random forest model on the radiomic features from all of the MRI sequences with demographic information using a lesion volume cutoff of 50 mm3. This model has a mean absolute error of 7.23 months (95% CI: [6.98, 13.43]) and a median absolute error of 5.98 months (95% CI: [5.26, 13.25]) in the validation set. For this model, the predicted age and actual age have a statistically significant association (p-value <0.001) in the validation set.

Removing inter-subject technical variability in magnetic resonance imaging studies.

Magnetic resonance imaging (MRI) intensities are acquired in arbitrary units, making scans non-comparable across sites and between subjects. Intensity normalization is a first step for the improvement of comparability of the images across subjects. However, we show that unwanted inter-scan variability associated with imaging site, scanner effect, and other technical artifacts is still present after standard intensity normalization in large multi-site neuroimaging studies. We propose RAVEL (Removal of Artificial Voxel Effect by Linear regression), a tool to remove residual technical variability after intensity normalization. As proposed by SVA and RUV [Leek and Storey, 2007, 2008, Gagnon-Bartsch and Speed, 2012], two batch effect correction tools largely used in genomics, we decompose the voxel intensities of images registered to a template into a biological component and an unwanted variation component. The unwanted variation component is estimated from a control region obtained from the cerebrospinal fluid (CSF), where intensities are known to be unassociated with disease status and other clinical covariates. We perform a singular value decomposition (SVD) of the control voxels to estimate factors of unwanted variation. We then estimate the unwanted factors using linear regression for every voxel of the brain and take the residuals as the RAVEL-corrected intensities. We assess the performance of RAVEL using T1-weighted (T1-w) images from more than 900 subjects with Alzheimer's disease (AD) and mild cognitive impairment (MCI), as well as healthy controls from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database. We compare RAVEL to two intensity-normalization-only methods: histogram matching and White Stripe. We show that RAVEL performs best at improving the replicability of the brain regions that are empirically found to be most associated with AD, and that these regions are significantly more present in structures impacted by AD (hippocampus, amygdala, parahippocampal gyrus, enthorinal area, and fornix stria terminals). In addition, we show that the RAVEL-corrected intensities have the best performance in distinguishing between MCI subjects and healthy subjects using the mean hippocampal intensity (AUC=67%), a marked improvement compared to results from intensity normalization alone (AUC=63% and 59% for histogram matching and White Stripe, respectively). RAVEL is promising for many other imaging modalities.

Statistical estimation of T1 relaxation times using conventional magnetic resonance imaging.

Quantitative T1 maps estimate T1 relaxation times and can be used to assess diffuse tissue abnormalities within normal-appearing tissue. T1 maps are popular for studying the progression and treatment of multiple sclerosis (MS). However, their inclusion in standard imaging protocols remains limited due to the additional scanning time and expert calibration required and susceptibility to bias and noise. Here, we propose a new method of estimating T1 maps using four conventional MR images, which are intensity-normalized using cerebellar gray matter as a reference tissue and related to T1 using a smooth regression model. Using cross-validation, we generate statistical T1 maps for 61 subjects with MS. The statistical maps are less noisy than the acquired maps and show similar reproducibility. Tests of group differences in normal-appearing white matter across MS subtypes give similar results using both methods.

Clinical 3-tesla FLAIR* MRI improves diagnostic accuracy in multiple sclerosis.

OBJECTIVE: To evaluate clinical fluid-attenuated inversion recovery (FLAIR)* 3T magnetic resonance imaging (MRI), which is sensitive to perivenular inflammatory demyelinating lesions, in diagnosing multiple sclerosis (MS). BACKGROUND: Central veins may be a distinguishing feature of MS lesions. FLAIR*, a combined contrast derived from clinical MRI scans, has not been studied as a clinical tool for diagnosing MS. METHODS: Two experienced MS neurologists evaluated 87 scan pairs (T2-FLAIR/FLAIR*), separately and side-by-side, from 68 MS cases, 8 healthy volunteers, and 11 individuals with other neurological diseases. Raters judged cases based on experience, published criteria, and a visual assessment of the "40% rule," whereby MS is favored if >40% of lesions demonstrate a central vein. Diagnostic accuracy was determined with area under the receiver operating characteristic curve (AUC), and inter-rater reliability was assessed with Cohen's kappa (κ). RESULTS: Diagnostic accuracy was high: rater 1, AUC 0.94 (95% confidence interval: 0.89, 0.97) for T2-FLAIR, 0.95 (0.92, 0.98) for FLAIR*; rater 2, 0.94 (0.90, 0.98) and 0.90 (0.85, 0.95). AUC improved when images were considered together: rater 1, 0.99 (0.98, 1.00); rater 2, 0.98 (0.96, 0.99). Inter-rater agreement was substantial for T2-FLAIR (κ = 0.68) and FLAIR* (κ = 0.74), despite low agreement on the 40% rule (κ = 0.47) ([Formula: see text] in all cases). CONCLUSIONS: Joint clinical evaluation of T2-FLAIR and FLAIR* images modestly improves diagnostic accuracy for MS and does not require counting lesions with central veins.

Quantitative Intracerebral Hemorrhage Localization.

BACKGROUND AND PURPOSE: The location of intracerebral hemorrhage (ICH) is currently described in a qualitative way; we provide a quantitative framework for estimating ICH engagement and its relevance to stroke outcomes. METHODS: We analyzed 111 patients with ICH from the Minimally Invasive Surgery Plus Recombinant-Tissue Plasminogen Activator for Intracerebral Evacuation (MISTIE) II clinical trial. We estimated ICH engagement at a population level using image registration of computed tomographic scans to a template and a previously labeled atlas. Predictive regions of National Institutes of Health Stroke Scale and Glasgow Coma Scale stroke severity scores, collected at enrollment, were estimated. RESULTS: The percent coverage of the ICH by these regions strongly outperformed the reader-labeled locations. The adjusted R(2) almost doubled from 0.129 (reader-labeled model) to 0.254 (quantitative location model) for National Institutes of Health Stroke Scale and more than tripled from 0.069 (reader-labeled model) to 0.214 (quantitative location model). A permutation test confirmed that the new predictive regions are more predictive than chance: P<0.001 for National Institutes of Health Stroke Scale and P<0.01 for Glasgow Coma Scale. CONCLUSIONS: Objective measures of ICH location and engagement using advanced computed tomographic imaging processing provide finer, objective, and more quantitative anatomic information than that provided by human readers. CLINICAL TRIAL REGISTRATION: URL: http://www.clinicaltrials.gov. Unique identifier: NCT00224770.

Scan-stratified case-control sampling for modeling blood-brain barrier integrity in multiple sclerosis.

Multiple sclerosis (MS) is an immune-mediated neurological disease that causes morbidity and disability. In patients with MS, the accumulation of lesions in the white matter of the brain is associated with disease progression and worse clinical outcomes. Breakdown of the blood-brain barrier in newer lesions is indicative of more active disease-related processes and is a primary outcome considered in clinical trials of treatments for MS. Such abnormalities in active MS lesions are evaluated in vivo using contrast-enhanced structural MRI, during which patients receive an intravenous infusion of a costly magnetic contrast agent. In some instances, the contrast agents can have toxic effects. Recently, local image regression techniques have been shown to have modest performance for assessing the integrity of the blood-brain barrier based on imaging without contrast agents. These models have centered on the problem of cross-sectional classification in which patients are imaged at a single study visit and pre-contrast images are used to predict post-contrast imaging. In this paper, we extend these methods to incorporate historical imaging information, and we find the proposed model to exhibit improved performance. We further develop scan-stratified case-control sampling techniques that reduce the computational burden of local image regression models, while respecting the low proportion of the brain that exhibits abnormal vascular permeability.

Increased susceptibility to vaginal simian/human immunodeficiency virus transmission in pig-tailed macaques coinfected with Chlamydia trachomatis and Trichomonas vaginalis.

BACKGROUND: Sexually transmitted infections (STIs) are associated with an increased risk of human immunodeficiency virus (HIV) infection, but their biological effect on HIV susceptibility is not fully understood. METHODS: Female pig-tailed macaques inoculated with Chlamydia trachomatis and Trichomonas vaginalis (n = 9) or medium (controls; n = 7) were repeatedly challenged intravaginally with SHIVSF162p3. Virus levels were evaluated by real-time polymerase chain reaction, plasma and genital cytokine levels by Luminex assays, and STI clinical signs by colposcopy. RESULTS: Simian/HIV (SHIV) susceptibility was enhanced in STI-positive macaques (P = .04, by the log-rank test; relative risk, 2.5 [95% confidence interval, 1.1-5.6]). All STI-positive macaques were SHIV infected, whereas 3 controls (43%) remained uninfected. Moreover, relative to STI-negative animals, SHIV infections occurred earlier in the menstrual cycle in STI-positive macaques (P = .01, by the Wilcoxon test). Levels of inflammatory cytokines (interferon γ, interleukin 6, and granulocyte colony-stimulating factor [G-CSF]) were higher in STI-positive macaques during STI inoculation and SHIV exposure periods (P ≤ .05, by the Wilcoxon test). CONCLUSIONS: C. trachomatis and T. vaginalis infection increase the susceptibility to SHIV, likely because of prolonged genital tract inflammation. These novel data demonstrate a biological link between these nonulcerative STIs and the risk of SHIV infection, supporting epidemiological associations of HIV and STIs. This study establishes a macaque model for studies of high-risk HIV transmission and prevention.

A multi-institutional machine learning algorithm for prognosticating facial nerve injury following microsurgical resection of vestibular schwannoma.

Vestibular schwannomas (VS) are the most common tumor of the skull base with available treatment options that carry a risk of iatrogenic injury to the facial nerve, which can significantly impact patients' quality of life. As facial nerve outcomes remain challenging to prognosticate, we endeavored to utilize machine learning to decipher predictive factors relevant to facial nerve outcomes following microsurgical resection of VS. A database of patient-, tumor- and surgery-specific features was constructed via retrospective chart review of 242 consecutive patients who underwent microsurgical resection of VS over a 7-year study period. This database was then used to train non-linear supervised machine learning classifiers to predict facial nerve preservation, defined as House-Brackmann (HB) I vs. facial nerve injury, defined as HB II-VI, as determined at 6-month outpatient follow-up. A random forest algorithm demonstrated 90.5% accuracy, 90% sensitivity and 90% specificity in facial nerve injury prognostication. A random variable (rv) was generated by randomly sampling a Gaussian distribution and used as a benchmark to compare the predictiveness of other features. This analysis revealed age, body mass index (BMI), case length and the tumor dimension representing tumor growth towards the brainstem as prognosticators of facial nerve injury. When validated via prospective assessment of facial nerve injury risk, this model demonstrated 84% accuracy. Here, we describe the development of a machine learning algorithm to predict the likelihood of facial nerve injury following microsurgical resection of VS. In addition to serving as a clinically applicable tool, this highlights the potential of machine learning to reveal non-linear relationships between variables which may have clinical value in prognostication of outcomes for high-risk surgical procedures.

Quantifying interictal intracranial EEG to predict focal epilepsy.

Introduction: Intracranial EEG (IEEG) is used for 2 main purposes, to determine: (1) if epileptic networks are amenable to focal treatment and (2) where to intervene. Currently these questions are answered qualitatively and sometimes differently across centers. There is a need for objective, standardized methods to guide surgical decision making and to enable large scale data analysis across centers and prospective clinical trials. Methods: We analyzed interictal data from 101 patients with drug resistant epilepsy who underwent presurgical evaluation with IEEG. We chose interictal data because of its potential to reduce the morbidity and cost associated with ictal recording. 65 patients had unifocal seizure onset on IEEG, and 36 were non-focal or multi-focal. We quantified the spatial dispersion of implanted electrodes and interictal IEEG abnormalities for each patient. We compared these measures against the "5 Sense Score (5SS)," a pre-implant estimate of the likelihood of focal seizure onset, and assessed their ability to predict the clinicians' choice of therapeutic intervention and the patient outcome. Results: The spatial dispersion of IEEG electrodes predicted network focality with precision similar to the 5SS (AUC = 0.67), indicating that electrode placement accurately reflected pre-implant information. A cross-validated model combining the 5SS and the spatial dispersion of interictal IEEG abnormalities significantly improved this prediction (AUC = 0.79; p<0.05). The combined model predicted ultimate treatment strategy (surgery vs. device) with an AUC of 0.81 and post-surgical outcome at 2 years with an AUC of 0.70. The 5SS, interictal IEEG, and electrode placement were not correlated and provided complementary information. Conclusions: Quantitative, interictal IEEG significantly improved upon pre-implant estimates of network focality and predicted treatment with precision approaching that of clinical experts. We present this study as an important step in building standardized, quantitative tools to guide epilepsy surgery.

Association of brain tissue cerebrospinal fluid fraction with age in healthy cognitively normal adults.

Background and purpose: Our objective was to apply multi-compartment T2 relaxometry in cognitively normal individuals aged 20-80 years to study the effect of aging on the parenchymal CSF fraction (CSFF), a potential measure of the subvoxel CSF space. Materials and methods: A total of 60 volunteers (age range, 22-80 years) were enrolled. Voxel-wise maps of short-T2 myelin water fraction (MWF), intermediate-T2 intra/extra-cellular water fraction (IEWF), and long-T2 CSFF were obtained using fast acquisition with spiral trajectory and adiabatic T2prep (FAST-T2) sequence and three-pool non-linear least squares fitting. Multiple linear regression analyses were performed to study the association between age and regional MWF, IEWF, and CSFF measurements, adjusting for sex and region of interest (ROI) volume. ROIs include the cerebral white matter (WM), cerebral cortex, and subcortical deep gray matter (GM). In each model, a quadratic term for age was tested using an ANOVA test. A Spearman's correlation between the normalized lateral ventricle volume, a measure of organ-level CSF space, and the regional CSFF, a measure of tissue-level CSF space, was computed. Results: Regression analyses showed that there was a statistically significant quadratic relationship with age for CSFF in the cortex (p = 0.018), MWF in the cerebral WM (p = 0.033), deep GM (p = 0.017) and cortex (p = 0.029); and IEWF in the deep GM (p = 0.033). There was a statistically highly significant positive linear relationship between age and regional CSFF in the cerebral WM (p < 0.001) and deep GM (p < 0.001). In addition, there was a statistically significant negative linear association between IEWF and age in the cerebral WM (p = 0.017) and cortex (p < 0.001). In the univariate correlation analysis, the normalized lateral ventricle volume correlated with the regional CSFF measurement in the cerebral WM (ρ = 0.64, p < 0.001), cortex (ρ = 0.62, p < 0.001), and deep GM (ρ = 0.66, p < 0.001). Conclusion: Our cross-sectional data demonstrate that brain tissue water in different compartments shows complex age-dependent patterns. Parenchymal CSFF, a measure of subvoxel CSF-like water in the brain tissue, is quadratically associated with age in the cerebral cortex and linearly associated with age in the cerebral deep GM and WM.

Mapping the relationship of white matter lesions to depression in multiple sclerosis.

Importance: Multiple sclerosis (MS) is an immune-mediated neurological disorder that affects nearly one million people in the United States. Up to 50% of patients with MS experience depression. Objective: To investigate how white matter network disruption is related to depression in MS. Design: Retrospective case-control study of participants who received research-quality 3-tesla neuroimaging as part of MS clinical care from 2010-2018. Analyses were performed from May 1 to September 30, 2022. Setting: Single-center academic medical specialty MS clinic. Participants: Participants with MS were identified via the electronic health record (EHR). All participants were diagnosed by an MS specialist and completed research-quality MRI at 3T. After excluding participants with poor image quality, 783 were included. Inclusion in the depression group (MS+Depression) required either: 1) ICD-10 depression diagnosis (F32-F34.*); 2) prescription of antidepressant medication; or 3) screening positive via Patient Health Questionnaire-2 (PHQ-2) or -9 (PHQ-9). Age- and sex-matched nondepressed comparators (MS-Depression) included persons with no depression diagnosis, no psychiatric medications, and were asymptomatic on PHQ-2/9. Exposure: Depression diagnosis. Main Outcomes and Measures: We first evaluated if lesions were preferentially located within the depression network compared to other brain regions. Next, we examined if MS+Depression patients had greater lesion burden, and if this was driven by lesions specifically in the depression network. Outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. Secondary measures included between-diagnosis lesion burden, stratified by brain network. Linear mixed-effects models were employed. Results: Three hundred-eighty participants met inclusion criteria, (232 MS+Depression: age[SD]=49[12], %females=86; 148 MS-Depression: age[SD]=47[13], %females=79). MS lesions preferentially affected fascicles within versus outside the depression network (ß=0.09, 95% CI=0.08-0.10, P<0.001). MS+Depression had more white matter lesion burden (ß=0.06, 95% CI=0.01-0.10, P=0.015); this was driven by lesions within the depression network (ß=0.02, 95% CI 0.003-0.040, P=0.020). Conclusions and Relevance: We provide new evidence supporting a relationship between white matter lesions and depression in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression had more disease than MS-Depression, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

Mapping the Relationship of White Matter Lesions to Depression in Multiple Sclerosis.

BACKGROUND: Multiple sclerosis (MS) is an immune-mediated neurological disorder, and up to 50% of patients experience depression. We investigated how white matter network disruption is related to depression in MS. METHODS: Using electronic health records, 380 participants with MS were identified. Depressed individuals (MS+Depression group; n = 232) included persons who had an ICD-10 depression diagnosis, had a prescription for antidepressant medication, or screened positive via Patient Health Questionnaire (PHQ)-2 or PHQ-9. Age- and sex-matched nondepressed individuals with MS (MS-Depression group; n = 148) included persons who had no prior depression diagnosis, had no psychiatric medication prescriptions, and were asymptomatic on PHQ-2 or PHQ-9. Research-quality 3T structural magnetic resonance imaging was obtained as part of routine care. We first evaluated whether lesions were preferentially located within the depression network compared with other brain regions. Next, we examined if MS+Depression patients had greater lesion burden and if this was driven by lesions in the depression network. Primary outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. RESULTS: MS lesions preferentially affected fascicles within versus outside the depression network (ß = 0.09, 95% CI = 0.08 to 0.10, p < .001). MS+Depression patients had more lesion burden (ß = 0.06, 95% CI = 0.01 to 0.10, p = .015); this was driven by lesions within the depression network (ß = 0.02, 95% CI = 0.003 to 0.040, p = .020). CONCLUSIONS: We demonstrated that lesion location and burden may contribute to depression comorbidity in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression patients had more disease than MS-Depression patients, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

Diagnostic accuracy of shuttle CT angiography (CTA) and helical CTA in the diagnosis of vasospasm.

PURPOSE: To evaluate the diagnostic accuracy of computed tomography angiography (CTA) acquired with shuttle technique (CTAs) and helical CTA (CTAh) for vasospasm, using digital subtraction angiography (DSA) obtained within 24 h as reference standard. METHODS: Thirty-six patients with suspected vasospasm in the setting of aneurysmal subarachnoid hemorrhage (ASAH, 30/36) or acute inflammatory/infectious conditions (6/36) who underwent CTAs (17/36) or CTAh (19/36) followed by DSA within 24 h were identified retrospectively. Presence of vasospasm in the proximal cerebral arterial segments was assessed qualitatively and semi-quantitatively on CTA and subsequent DSA. Sensitivity, specificity, and receiver operating characteristic (ROC) curves were calculated. Inter-rater variability was assessed using Cohen's kappa. RESULTS: On CTAs, 35% of patients had low and 65% had high vasospasm burden. On CTAh, 37% had low and 63% had high vasospasm burden. ROC analysis demonstrated an AUC of 0.87 for CTAs (95%CI 0.67-1.00, p = 0.015) and 0.88 for CTAh (0.72-1.00, p = 0.028). Cohen's kappa was 0.843 (95%CI 0.548-1.000). Thresholding with Youden's J index, CTAs had a sensitivity of 85.71% (95%CI 48.69 to 99.27) and specificity of 66.67% (35.42 to 87.94). CTAh had sensitivity of 100% (56.55 to 100.00) and specificity of 78.57% (52.41 to 92.43). CONCLUSION: CTAs and CTAh yielded similar sensitivity, specificity, and AUC values on ROC analysis for the detection of vasospasm using DSA as reference standard. Our findings suggest that CTAs is a promising alternative to CTAh especially in patients requiring serial imaging, given its potential advantages of decreased radiation exposure, contrast dose, and cost-effectiveness.

Quantitative susceptibility mapping versus phase imaging to identify multiple sclerosis iron rim lesions with demyelination.

BACKGROUND AND PURPOSE: To compare quantitative susceptibility mapping (QSM) and high-pass-filtered (HPF) phase imaging for (1) identifying chronic active rim lesions with more myelin damage and (2) distinguishing patients with increased clinical disability in multiple sclerosis. METHODS: Eighty patients were scanned with QSM for paramagnetic rim detection and Fast Acquisition with Spiral Trajectory and T2prep for myelin water fraction (MWF). Chronic lesions were classified based on the presence/absence of rim on HPF and QSM images. A lesion-level linear mixed-effects model with MWF as the outcome was used to compare myelin damage among the lesion groups. A multiple patient-level linear regression model was fit to establish the association between Expanded Disease Status Scale (EDSS) and the log of the number of rim lesions. RESULTS: Of 2062 lesions, 188 (9.1%) were HPF rim+/QSM rim+, 203 (9.8%) were HPF rim+/QSM rim-, and the remainder had no rim. In the linear mixed-effects model, HPF rim+/QSM rim+ lesions had significantly lower MWF than both HPF rim+/QSM rim- (p < .001) and HPF rim-/QSM rim- (p < .001) lesions, while the MWF difference between HPF rim+/QSM rim- and HPF rim-/QSM rim- lesions was not statistically significant (p = .130). Holding all other factors constant, the log number of QSM rim+ lesion was associated with EDSS increase (p = .044). The association between the log number of HPF rim+ lesions and EDSS was not statistically significant (p = .206). CONCLUSIONS: QSM identifies paramagnetic rim lesions that on average have more myelin damage and stronger association with clinical disability than those detected by phase imaging.

Quantitative Water Permeability Mapping of Blood-Brain-Barrier Dysfunction in Aging.

Blood-brain-barrier (BBB) dysfunction is a hallmark of aging and aging-related disorders, including cerebral small vessel disease and Alzheimer's disease. An emerging biomarker of BBB dysfunction is BBB water exchange rate (kW) as measured by diffusion-weighted arterial spin labeling (DW-ASL) MRI. We developed an improved DW-ASL sequence for Quantitative Permeability Mapping and evaluated whole brain and region-specific kW in a cohort of 30 adults without dementia across the age spectrum. In this cross-sectional study, we found higher kW values in the cerebral cortex (mean = 81.51 min-1, SD = 15.54) compared to cerebral white matter (mean = 75.19 min-1, SD = 13.85) (p < 0.0001). We found a similar relationship for cerebral blood flow (CBF), concordant with previously published studies. Multiple linear regression analysis with kW as an outcome showed that age was statistically significant in the cerebral cortex (p = 0.013), cerebral white matter (p = 0.033), hippocampi (p = 0.043), orbitofrontal cortices (p = 0.042), and precunei cortices (p = 0.009), after adjusting for sex and number of vascular risk factors. With CBF as an outcome, age was statistically significant only in the cerebral cortex (p = 0.026) and precunei cortices (p = 0.020). We further found moderate negative correlations between white matter hyperintensity (WMH) kW and WMH volume (r = -0.51, p = 0.02), and normal-appearing white matter (NAWM) and WMH volume (r = -0.44, p = 0.05). This work illuminates the relationship between BBB water exchange and aging and may serve as the basis for BBB-targeted therapies for aging-related brain disorders.

Global cerebrospinal fluid as a zero-reference regularization for brain quantitative susceptibility mapping.

BACKGROUND AND PURPOSE: The objective ofthis study was to demonstrate a global cerebrospinal fluid (CSF) method for a consistent and automated zero referencing of brain quantitative susceptibility mapping (QSM). METHODS: Whole brain CSF mask was automatically segmented by thresholding the gradient echo transverse relaxation ( R2∗) map, and regularization was employed to enforce uniform susceptibility distribution within the CSF volume in the field-to-susceptibility inversion. This global CSF regularization method was compared with a prior ventricular CSF regularization. Both reconstruction methods were compared in a repeatability study of 12 healthy subjects using t-test on susceptibility measurements, and in patient studies of 17 multiple sclerosis (MS) and 10 Parkinson's disease (PD) patients using Wilcoxon rank-sum test on radiological scores. RESULTS: In scan-rescan experiments, global CSF regularization provided more consistent CSF volume as well as higher repeatability of QSM measurements than ventricular CSF regularization with a smaller bias: -2.7 parts per billion (ppb) versus -0.13 ppb (t-test p<0.05) and a narrower 95% limits of agreement: [-7.25, 6.99] ppb versus [-16.60, 11.19 ppb] (f-test p<0.05). In PD and MS patients, global CSF regularization reduced smoothly varying shadow artifacts and significantly improved the QSM quality score (p<0.001). CONCLUSIONS: The proposed whole brain CSF method for QSM zero referencing improves repeatability and image quality of brain QSM compared to the ventricular CSF method.

QSMRim-Net: Imbalance-aware learning for identification of chronic active multiple sclerosis lesions on quantitative susceptibility maps.

BACKGROUND AND PURPOSE: Chronic active multiple sclerosis (MS) lesions are characterized by a paramagnetic rim at the edge of the lesion and are associated with increased disability in patients. Quantitative susceptibility mapping (QSM) is an MRI technique that is sensitive to chronic active lesions, termed rim + lesions on the QSM. We present QSMRim-Net, a data imbalance-aware deep neural network that fuses lesion-level radiomic and convolutional image features for automated identification of rim + lesions on QSM. METHODS: QSM and T2-weighted-Fluid-Attenuated Inversion Recovery (T2-FLAIR) MRI of the brain were collected at 3 T for 172 MS patients. Rim + lesions were manually annotated by two human experts, followed by consensus from a third expert, for a total of 177 rim + and 3986 rim negative (rim-) lesions. Our automated rim + detection algorithm, QSMRim-Net, consists of a two-branch feature extraction network and a synthetic minority oversampling network to classify rim + lesions. The first network branch is for image feature extraction from the QSM and T2-FLAIR, and the second network branch is a fully connected network for QSM lesion-level radiomic feature extraction. The oversampling network is designed to increase classification performance with imbalanced data. RESULTS: On a lesion-level, in a five-fold cross validation framework, the proposed QSMRim-Net detected rim + lesions with a partial area under the receiver operating characteristic curve (pROC AUC) of 0.760, where clinically relevant false positive rates of less than 0.1 were considered. The method attained an area under the precision recall curve (PR AUC) of 0.704. QSMRim-Net out-performed other state-of-the-art methods applied to the QSM on both pROC AUC and PR AUC. On a subject-level, comparing the predicted rim + lesion count and the human expert annotated count, QSMRim-Net achieved the lowest mean square error of 0.98 and the highest correlation of 0.89 (95% CI: 0.86, 0.92). CONCLUSION: This study develops a novel automated deep neural network for rim + MS lesion identification using T2-FLAIR and QSM images.

Brain oxygen extraction fraction mapping in patients with multiple sclerosis.

We aimed to demonstrate the feasibility of whole brain oxygen extraction fraction (OEF) mapping for measuring lesion specific and regional OEF abnormalities in multiple sclerosis (MS) patients. In 22 MS patients and 11 healthy controls (HC), OEF and neural tissue susceptibility (χn) maps were computed from MRI multi-echo gradient echo data. In MS patients, 80 chronic active lesions with hyperintense rim on quantitative susceptibility mapping were identified, and the mean OEF and χn within the rim and core were compared using linear mixed-effect model analysis. The rim showed higher OEF and χn than the core: relative to their adjacent normal appearing white matter, OEF contrast = -6.6 ± 7.0% vs. -9.8 ± 7.8% (p < 0.001) and χn contrast = 33.9 ± 20.3 ppb vs. 25.7 ± 20.5 ppb (p = 0.017). Between MS and HC, OEF and χn were compared using a linear regression model in subject-based regions of interest. In the whole brain, compared to HC, MS had lower OEF, 30.4 ± 3.3% vs. 21.4 ± 4.4% (p < 0.001), and higher χn, -23.7 ± 7.0 ppb vs. -11.3 ± 7.7 ppb (p = 0.018). Our feasibility study suggests that OEF may serve as a useful quantitative marker of tissue oxygen utilization in MS.