Investigating the relationship between thalamic iron concentration and disease severity in secondary progressive multiple sclerosis using quantitative susceptibility mapping: Cross-sectional analysis from the MS-STAT2 randomised controlled trial.

Background: Deep grey matter pathology is a key driver of disability worsening in people with multiple sclerosis. Quantitative susceptibility mapping (QSM) is an advanced magnetic resonance imaging (MRI) technique which quantifies local magnetic susceptibility from variations in phase produced by changes in the local magnetic field. In the deep grey matter, susceptibility has previously been validated against tissue iron concentration. However, it currently remains unknown whether susceptibility is abnormal in older progressive MS cohorts, and whether it correlates with disability. Objectives: To investigate differences in mean regional susceptibility in deep grey matter between people with secondary progressive multiple sclerosis (SPMS) and healthy controls; to examine in patients the relationships between deep grey matter susceptibility and clinical and imaging measures of disease severity. Methods: Baseline data from a subgroup of the MS-STAT2 trial (simvastatin vs. placebo in SPMS, NCT03387670) were included. The subgroup underwent clinical assessments and an advanced MRI protocol at 3T. A cohort of age-matched healthy controls underwent the same MRI protocol. Susceptibility maps were reconstructed using a robust QSM pipeline from multi-echo 3D gradient-echo sequence. Regions of interest (ROIs) in the thalamus, globus pallidus and putamen were segmented from 3D T1-weighted images, and lesions segmented from 3D fluid-attenuated inversion recovery images. Linear regression was used to compare susceptibility from ROIs between patients and controls, adjusting for age and sex. Where significant differences were found, we further examined the associations between ROI susceptibility and clinical and imaging measures of MS severity. Results: 149 SPMS (77% female; mean age: 53 yrs; median Expanded Disability Status Scale (EDSS): 6.0 [interquartile range 4.5-6.0]) and 33 controls (52% female, mean age: 57) were included.Thalamic susceptibility was significantly lower in SPMS compared to controls: mean (SD) 28.6 (12.8) parts per billion (ppb) in SPMS vs. 39.2 (12.7) ppb in controls; regression coefficient: -12.0 [95% confidence interval: -17.0 to -7.1], p < 0.001. In contrast, globus pallidus and putamen susceptibility were similar between both groups.In SPMS, a 10 ppb lower thalamic susceptibility was associated with a +0.13 [+0.01 to +0.24] point higher EDSS (p < 0.05), a -2.4 [-3.8 to -1.0] point lower symbol digit modality test (SDMT, p = 0.001), and a -2.4 [-3.7 to -1.1] point lower Sloan low contrast acuity, 2.5% (p < 0.01).Lower thalamic susceptibility was also strongly associated with a higher T2 lesion volume (T2LV, p < 0.001) and lower normalised whole brain, deep grey matter and thalamic volumes (all p < 0.001). Conclusions: The reduced thalamic susceptibility found in SPMS compared to controls suggests that thalamic iron concentrations are lower at this advanced stage of the disease. The observed relationships between lower thalamic susceptibility and more severe physical, cognitive and visual disability suggests that reductions in thalamic iron may correlate with important mechanisms of clinical disease progression. Such mechanisms appear to intimately link reductions in thalamic iron with higher T2LV and the development of thalamic atrophy, encouraging further research into QSM-derived thalamic susceptibility as a biomarker of disease severity in SPMS.

Body size interacts with the structure of the central nervous system: A multi-center in vivo neuroimaging study.

Clinical research emphasizes the implementation of rigorous and reproducible study designs that rely on between-group matching or controlling for sources of biological variation such as subject's sex and age. However, corrections for body size (i.e. height and weight) are mostly lacking in clinical neuroimaging designs. This study investigates the importance of body size parameters in their relationship with spinal cord (SC) and brain magnetic resonance imaging (MRI) metrics. Data were derived from a cosmopolitan population of 267 healthy human adults (age 30.1±6.6 years old, 125 females). We show that body height correlated strongly or moderately with brain gray matter (GM) volume, cortical GM volume, total cerebellar volume, brainstem volume, and cross-sectional area (CSA) of cervical SC white matter (CSA-WM; 0.44≤r≤0.62). In comparison, age correlated weakly with cortical GM volume, precentral GM volume, and cortical thickness (-0.21≥r≥-0.27). Body weight correlated weakly with magnetization transfer ratio in the SC WM, dorsal columns, and lateral corticospinal tracts (-0.20≥r≥-0.23). Body weight further correlated weakly with the mean diffusivity derived from diffusion tensor imaging (DTI) in SC WM (r=-0.20) and dorsal columns (-0.21), but only in males. CSA-WM correlated strongly or moderately with brain volumes (0.39≤r≤0.64), and weakly with precentral gyrus thickness and DTI-based fractional anisotropy in SC dorsal columns and SC lateral corticospinal tracts (-0.22≥r≥-0.25). Linear mixture of sex and age explained 26±10% of data variance in brain volumetry and SC CSA. The amount of explained variance increased at 33±11% when body height was added into the mixture model. Age itself explained only 2±2% of such variance. In conclusion, body size is a significant biological variable. Along with sex and age, body size should therefore be included as a mandatory variable in the design of clinical neuroimaging studies examining SC and brain structure.

An accessible deep learning tool for voxel-wise classification of brain malignancies from perfusion MRI.

Noninvasive differential diagnosis of brain tumors is currently based on the assessment of magnetic resonance imaging (MRI) coupled with dynamic susceptibility contrast (DSC). However, a definitive diagnosis often requires neurosurgical interventions that compromise patients' quality of life. We apply deep learning on DSC images from histology-confirmed patients with glioblastoma, metastasis, or lymphoma. The convolutional neural network trained on ∼50,000 voxels from 40 patients provides intratumor probability maps that yield clinical-grade diagnosis. Performance is tested in 400 additional cases and an external validation cohort of 128 patients. The tool reaches a three-way accuracy of 0.78, superior to the conventional MRI metrics cerebral blood volume (0.55) and percentage of signal recovery (0.59), showing high value as a support diagnostic tool. Our open-access software, Diagnosis In Susceptibility Contrast Enhancing Regions for Neuro-oncology (DISCERN), demonstrates its potential in aiding medical decisions for brain tumor diagnosis using standard-of-care MRI.

What contributes to disability in progressive MS? A brain and cervical cord-matched quantitative MRI study.

BACKGROUND: We assessed the ability of a brain-and-cord-matched quantitative magnetic resonance imaging (qMRI) protocol to differentiate patients with progressive multiple sclerosis (PMS) from controls, in terms of normal-appearing (NA) tissue abnormalities, and explain disability. METHODS: A total of 27 patients and 16 controls were assessed on the Expanded Disability Status Scale (EDSS), 25-foot timed walk (TWT), 9-hole peg (9HPT) and symbol digit modalities (SDMT) tests. All underwent 3T brain and (C2-C3) cord structural imaging and qMRI (relaxometry, quantitative magnetisation transfer, multi-shell diffusion-weighted imaging), using a fast brain-and-cord-matched protocol with brain-and-cord-unified imaging readouts. Lesion and NA-tissue volumes and qMRI metrics reflecting demyelination and axonal loss were obtained. Random forest analyses identified the most relevant volumetric/qMRI measures to clinical outcomes. Confounder-adjusted linear regression estimated the actual MRI-clinical associations. RESULTS: Several qMRI/volumetric differences between patients and controls were observed (p < 0.01). Higher NA-deep grey matter quantitative-T1 (EDSS: beta = 7.96, p = 0.006; 9HPT: beta = -0.09, p = 0.004), higher NA-white matter orientation dispersion index (TWT: beta = -3.21, p = 0.005; SDMT: beta = -847.10, p < 0.001), lower whole-cord bound pool fraction (9HPT: beta = 0.79, p = 0.001) and higher NA-cortical grey matter quantitative-T1 (SDMT = -94.31, p < 0.001) emerged as particularly relevant predictors of greater disability. CONCLUSION: Fast brain-and-cord-matched qMRI protocols are feasible and identify demyelination - combined with other mechanisms - as key for disability accumulation in PMS.

Identification of Precise 3D CT Radiomics for Habitat Computation by Machine Learning in Cancer.

Purpose To identify precise three-dimensional radiomics features in CT images that enable computation of stable and biologically meaningful habitats with machine learning for cancer heterogeneity assessment. Materials and Methods This retrospective study included 2436 liver or lung lesions from 605 CT scans (November 2010-December 2021) in 331 patients with cancer (mean age, 64.5 years ± 10.1 [SD]; 185 male patients). Three-dimensional radiomics were computed from original and perturbed (simulated retest) images with different combinations of feature computation kernel radius and bin size. The lower 95% confidence limit (LCL) of the intraclass correlation coefficient (ICC) was used to measure repeatability and reproducibility. Precise features were identified by combining repeatability and reproducibility results (LCL of ICC ≥ 0.50). Habitats were obtained with Gaussian mixture models in original and perturbed data using precise radiomics features and compared with habitats obtained using all features. The Dice similarity coefficient (DSC) was used to assess habitat stability. Biologic correlates of CT habitats were explored in a case study, with a cohort of 13 patients with CT, multiparametric MRI, and tumor biopsies. Results Three-dimensional radiomics showed poor repeatability (LCL of ICC: median [IQR], 0.442 [0.312-0.516]) and poor reproducibility against kernel radius (LCL of ICC: median [IQR], 0.440 [0.33-0.526]) but excellent reproducibility against bin size (LCL of ICC: median [IQR], 0.929 [0.853-0.988]). Twenty-six radiomics features were precise, differing in lung and liver lesions. Habitats obtained with precise features (DSC: median [IQR], 0.601 [0.494-0.712] and 0.651 [0.52-0.784] for lung and liver lesions, respectively) were more stable than those obtained with all features (DSC: median [IQR], 0.532 [0.424-0.637] and 0.587 [0.465-0.703] for lung and liver lesions, respectively; P < .001). In the case study, CT habitats correlated quantitatively and qualitatively with heterogeneity observed in multiparametric MRI habitats and histology. Conclusion Precise three-dimensional radiomics features were identified on CT images that enabled tumor heterogeneity assessment through stable tumor habitat computation. Keywords: CT, Diffusion-weighted Imaging, Dynamic Contrast-enhanced MRI, MRI, Radiomics, Unsupervised Learning, Oncology, Liver, Lung Supplemental material is available for this article. © RSNA, 2024 See also the commentary by Sagreiya in this issue.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology.

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

Feasibility of in vivo multi-parametric quantitative magnetic resonance imaging of the healthy sciatic nerve with a unified signal readout protocol.

Magnetic resonance neurography (MRN) has been used successfully over the years to investigate the peripheral nervous system (PNS) because it allows early detection and precise localisation of neural tissue damage. However, studies demonstrating the feasibility of combining MRN with multi-parametric quantitative magnetic resonance imaging (qMRI) methods, which provide more specific information related to nerve tissue composition and microstructural organisation, can be invaluable. The translation of emerging qMRI methods previously validated in the central nervous system to the PNS offers real potential to characterise in patients in vivo the underlying pathophysiological mechanisms involved in a plethora of conditions of the PNS. The aim of this study was to assess the feasibility of combining MRN with qMRI to measure diffusion, magnetisation transfer and relaxation properties of the healthy sciatic nerve in vivo using a unified signal readout protocol. The reproducibility of the multi-parametric qMRI protocol as well as normative qMRI measures in the healthy sciatic nerve are reported. The findings presented herein pave the way to the practical implementation of joint MRN-qMRI in future studies of pathological conditions affecting the PNS.

Patterns of inflammation, microstructural alterations, and sodium accumulation define multiple sclerosis subtypes after 15 years from onset.

Introduction: Conventional MRI is routinely used for the characterization of pathological changes in multiple sclerosis (MS), but due to its lack of specificity is unable to provide accurate prognoses, explain disease heterogeneity and reconcile the gap between observed clinical symptoms and radiological evidence. Quantitative MRI provides measures of physiological abnormalities, otherwise invisible to conventional MRI, that correlate with MS severity. Analyzing quantitative MRI measures through machine learning techniques has been shown to improve the understanding of the underlying disease by better delineating its alteration patterns. Methods: In this retrospective study, a cohort of healthy controls (HC) and MS patients with different subtypes, followed up 15 years from clinically isolated syndrome (CIS), was analyzed to produce a multi-modal set of quantitative MRI features encompassing relaxometry, microstructure, sodium ion concentration, and tissue volumetry. Random forest classifiers were used to train a model able to discriminate between HC, CIS, relapsing remitting (RR) and secondary progressive (SP) MS patients based on these features and, for each classification task, to identify the relative contribution of each MRI-derived tissue property to the classification task itself. Results and discussion: Average classification accuracy scores of 99 and 95% were obtained when discriminating HC and CIS vs. SP, respectively; 82 and 83% for HC and CIS vs. RR; 76% for RR vs. SP, and 79% for HC vs. CIS. Different patterns of alterations were observed for each classification task, offering key insights in the understanding of MS phenotypes pathophysiology: atrophy and relaxometry emerged particularly in the classification of HC and CIS vs. MS, relaxometry within lesions in RR vs. SP, sodium ion concentration in HC vs. CIS, and microstructural alterations were involved across all tasks.

Multimodal Analysis of the Visual Pathways in Friedreich's Ataxia Reveals Novel Biomarkers.

BACKGROUND: Optic neuropathy is a near ubiquitous feature of Friedreich's ataxia (FRDA). Previous studies have examined varying aspects of the anterior and posterior visual pathways but none so far have comprehensively evaluated the heterogeneity of degeneration across different areas of the retina, changes to the macula layers and combined these with volumetric MRI studies of the visual cortex and frataxin level. METHODS: We investigated 62 genetically confirmed FRDA patients using an integrated approach as part of an observational cohort study. We included measurement of frataxin protein levels, clinical evaluation of visual and neurological function, optical coherence tomography to determine retinal nerve fibre layer thickness and macular layer volume and volumetric brain MRI. RESULTS: We demonstrate that frataxin level correlates with peripapillary retinal nerve fibre layer thickness and that retinal sectors differ in their degree of degeneration. We also shown that retinal nerve fibre layer is thinner in FRDA patients than controls and that this thinning is influenced by the AAO and GAA1. Furthermore we show that the ganglion cell and inner plexiform layers are affected in FRDA. Our MRI data indicate that there are borderline correlations between retinal layers and areas of the cortex involved in visual processing. CONCLUSION: Our study demonstrates the uneven distribution of the axonopathy in the retinal nerve fibre layer and highlight the relative sparing of the papillomacular bundle and temporal sectors. We show that thinning of the retinal nerve fibre layer is associated with frataxin levels, supporting the use the two biomarkers in future clinical trials design. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Differentiating Multiple Sclerosis From AQP4-Neuromyelitis Optica Spectrum Disorder and MOG-Antibody Disease With Imaging.

BACKGROUND AND OBJECTIVES: Relapsing-remitting multiple sclerosis (RRMS), aquaporin-4 antibody-positive neuromyelitis optica spectrum disorder (AQP4-NMOSD), and myelin oligodendrocyte glycoprotein antibody-associated disease (MOGAD) may have overlapping clinical features. There is an unmet need for imaging markers that differentiate between them when serologic testing is unavailable or ambiguous. We assessed whether imaging characteristics typical of MS discriminate RRMS from AQP4-NMOSD and MOGAD, alone and in combination. METHODS: Adult, nonacute patients with RRMS, APQ4-NMOSD, and MOGAD and healthy controls were prospectively recruited at the National Hospital for Neurology and Neurosurgery (London, United Kingdom) and the Walton Centre (Liverpool, United Kingdom) between 2014 and 2019. They underwent conventional and advanced brain, cord, and optic nerve MRI and optical coherence tomography (OCT). RESULTS: A total of 91 consecutive patients (31 RRMS, 30 APQ4-NMOSD, and 30 MOGAD) and 34 healthy controls were recruited. The most accurate measures differentiating RRMS from AQP4-NMOSD were the proportion of lesions with the central vein sign (CVS) (84% vs 33%, accuracy/specificity/sensitivity: 91/88/93%, p < 0.001), followed by cortical lesions (median: 2 [range: 1-14] vs 1 [0-1], accuracy/specificity/sensitivity: 84/90/77%, p = 0.002) and white matter lesions (mean: 39.07 [±25.8] vs 9.5 [±14], accuracy/specificity/sensitivity: 78/84/73%, p = 0.001). The combination of higher proportion of CVS, cortical lesions, and optic nerve magnetization transfer ratio reached the highest accuracy in distinguishing RRMS from AQP4-NMOSD (accuracy/specificity/sensitivity: 95/92/97%, p < 0.001). The most accurate measures favoring RRMS over MOGAD were white matter lesions (39.07 [±25.8] vs 1 [±2.3], accuracy/specificity/sensitivity: 94/94/93%, p = 0.006), followed by cortical lesions (2 [1-14] vs 1 [0-1], accuracy/specificity/sensitivity: 84/97/71%, p = 0.004), and retinal nerve fiber layer thickness (RNFL) (mean: 87.54 [±13.83] vs 75.54 [±20.33], accuracy/specificity/sensitivity: 80/79/81%, p = 0.009). Higher cortical lesion number combined with higher RNFL thickness best differentiated RRMS from MOGAD (accuracy/specificity/sensitivity: 84/92/77%, p < 0.001). DISCUSSION: Cortical lesions, CVS, and optic nerve markers achieve a high accuracy in distinguishing RRMS from APQ4-NMOSD and MOGAD. This information may be useful in clinical practice, especially outside the acute phase and when serologic testing is ambiguous or not promptly available. CLASSIFICATION OF EVIDENCE: This study provides Class II evidence that selected conventional and advanced brain, cord, and optic nerve MRI and OCT markers distinguish adult patients with RRMS from AQP4-NMOSD and MOGAD.

SENSE EPI reconstruction with 2D phase error correction and channel-wise noise removal.

PURPOSE: To develop a robust reconstruction pipeline for EPI data that enables 2D Nyquist phase error correction using sensitivity encoding without incurring major noise artifacts in low SNR data. METHODS: SENSE with 2D phase error correction (PEC-SENSE) was combined with channel-wise noise removal using Marcenko-Pastur principal component analysis (MPPCA) to simultaneously eliminate Nyquist ghost artifacts in EPI data and mitigate the noise amplification associated with phase correction using parallel imaging. The proposed pipeline (coined SPECTRE) was validated in phantom DW-EPI data using the accuracy and precision of diffusion metrics; ground truth values were obtained from data acquired with a spin echo readout. Results from the SPECTRE pipeline were compared against PEC-SENSE reconstructions with three alternate denoising strategies: (i) no denoising; (ii) denoising of magnitude data after image formation; (iii) denoising of complex data after image formation. SPECTRE was then tested using high b $$ b $$ -value (i.e., low SNR) diffusion data (up to b = 3000 $$ b=3000 $$ s/mm 2 $$ {}^2 $$ ) in four healthy subjects. RESULTS: Noise amplification associated with phase error correction incurred a 23% bias in phantom mean diffusivity (MD) measurements. Phantom MD estimates using the SPECTRE pipeline were within 8% of the ground truth value. In healthy volunteers, the SPECTRE pipeline visibly corrected Nyquist ghost artifacts and reduced associated noise amplification in high b $$ b $$ -value data. CONCLUSION: The proposed reconstruction pipeline is effective in correcting low SNR data, and improves the accuracy and precision of derived diffusion metrics.

Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla.

PURPOSE: To compare different multi-echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi-echo gradient-recalled echo signal phase information, either before or after applying Laplacian-base methods (LBMs) for phase unwrapping or background field removal. METHODS: Multi-echo gradient-recalled echo data were simulated in a numerical head phantom, and multi-echo gradient-recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image-based estimation of gradient-recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE-dependent phase and then combined multiple TEs via either TE-weighted or SNR-weighted averaging; and 2 calculated TE-dependent susceptibility maps via either multi-step or single-step QSM and then combined multiple TEs via magnitude-weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland-Altman) and regional differences between the means (nonparametric tests). RESULTS: Nonlinearly fitting the multi-echo phase over TEs before applying LBMs provided the highest regional accuracy of χ $$ \chi $$ and the lowest phase noise propagation compared to averaging the LBM-processed TE-dependent phase. This result was especially pertinent in high-susceptibility venous regions. CONCLUSION: For multi-echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

Histo-MRI map study protocol: a prospective cohort study mapping MRI to histology for biomarker validation and prediction of prostate cancer.

INTRODUCTION: Multiparametric MRI (mpMRI) is now widely used to risk stratify men with a suspicion of prostate cancer and identify suspicious regions for biopsy. However, the technique has modest specificity and a high false-positive rate, especially in men with mpMRI scored as indeterminate (3/5) or likely (4/5) to have clinically significant cancer (csPCa) (Gleason ≥3+4). Advanced MRI techniques have emerged which seek to improve this characterisation and could predict biopsy results non-invasively. Before these techniques are translated clinically, robust histological and clinical validation is required. METHODS AND ANALYSIS: This study aims to clinically validate two advanced MRI techniques in a prospectively recruited cohort of men suspected of prostate cancer. Histological analysis of men undergoing biopsy or prostatectomy will be used for biological validation of biomarkers derived from Vascular and Extracellular Restricted Diffusion for Cytometry in Tumours and Luminal Water imaging. In particular, prostatectomy specimens will be processed using three-dimension printed patient-specific moulds to allow for accurate MRI and histology mapping. The index tests will be compared with the histological reference standard to derive false positive rate and true positive rate for men with mpMRI scores which are indeterminate (3/5) or likely (4/5) to have clinically significant prostate cancer (csPCa). Histopathological validation from both biopsy and prostatectomy samples will provide the best ground truth in validating promising MRI techniques which could predict biopsy results and help avoid unnecessary biopsies in men suspected of prostate cancer. ETHICS AND DISSEMINATION: Ethical approval was granted by the London-Queen Square Research Ethics Committee (19/LO/1803) on 23 January 2020. Results from the study will be presented at conferences and submitted to peer-reviewed journals for publication. Results will also be available on ClinicalTrials.gov. TRIAL REGISTRATION NUMBER: NCT04792138.

Comparison of multicenter MRI protocols for visualizing the spinal cord gray matter.

PURPOSE: Spinal cord gray-matter imaging is valuable for a number of applications, but remains challenging. The purpose of this work was to compare various MRI protocols at 1.5 T, 3 T, and 7 T for visualizing the gray matter. METHODS: In vivo data of the cervical spinal cord were collected from nine different imaging centers. Data processing consisted of automatically segmenting the spinal cord and its gray matter and co-registering back-to-back scans. We computed the SNR using two methods (SNR_single using a single scan and SNR_diff using the difference between back-to-back scans) and the white/gray matter contrast-to-noise ratio per unit time. Synthetic phantom data were generated to evaluate the metrics performance. Experienced radiologists qualitatively scored the images. We ran the same processing on an open-access multicenter data set of the spinal cord MRI (N = 267 participants). RESULTS: Qualitative assessments indicated comparable image quality for 3T and 7T scans. Spatial resolution was higher at higher field strength, and image quality at 1.5 T was found to be moderate to low. The proposed quantitative metrics were found to be robust to underlying changes to the SNR and contrast; however, the SNR_single method lacked accuracy when there were excessive partial-volume effects. CONCLUSION: We propose quality assessment criteria and metrics for gray-matter visualization and apply them to different protocols. The proposed criteria and metrics, the analyzed protocols, and our open-source code can serve as a benchmark for future optimization of spinal cord gray-matter imaging protocols.

Voxel-level analysis of normalized DSC-PWI time-intensity curves: a potential generalizable approach and its proof of concept in discriminating glioblastoma and metastasis.

OBJECTIVE: Standard DSC-PWI analyses are based on concrete parameters and values, but an approach that contemplates all points in the time-intensity curves and all voxels in the region-of-interest may provide improved information, and more generalizable models. Therefore, a method of DSC-PWI analysis by means of normalized time-intensity curves point-by-point and voxel-by-voxel is constructed, and its feasibility and performance are tested in presurgical discrimination of glioblastoma and metastasis. METHODS: In this retrospective study, patients with histologically confirmed glioblastoma or solitary-brain-metastases and presurgical-MR with DSC-PWI (August 2007-March 2020) were retrieved. The enhancing tumor and immediate peritumoral region were segmented on CE-T1wi and coregistered to DSC-PWI. Time-intensity curves of the segmentations were normalized to normal-appearing white matter. For each participant, average and all-voxel-matrix of normalized-curves were obtained. The 10 best discriminatory time-points between each type of tumor were selected. Then, an intensity-histogram analysis on each of these 10 time-points allowed the selection of the best discriminatory voxel-percentile for each. Separate classifier models were trained for enhancing tumor and peritumoral region using binary logistic regressions. RESULTS: A total of 428 patients (321 glioblastomas, 107 metastases) fulfilled the inclusion criteria (256 men; mean age, 60 years; range, 20-86 years). Satisfactory results were obtained to segregate glioblastoma and metastases in training and test sets with AUCs 0.71-0.83, independent accuracies 65-79%, and combined accuracies up to 81-88%. CONCLUSION: This proof-of-concept study presents a different perspective on brain MR DSC-PWI evaluation by the inclusion of all time-points of the curves and all voxels of segmentations to generate robust diagnostic models of special interest in heterogeneous diseases and populations. The method allows satisfactory presurgical segregation of glioblastoma and metastases. KEY POINTS: • An original approach to brain MR DSC-PWI analysis, based on a point-by-point and voxel-by-voxel assessment of normalized time-intensity curves, is presented. • The method intends to extract optimized information from MR DSC-PWI sequences by impeding the potential loss of information that may represent the standard evaluation of single concrete perfusion parameters (cerebral blood volume, percentage of signal recovery, or peak height) and values (mean, maximum, or minimum). • The presented approach may be of special interest in technically heterogeneous samples, and intrinsically heterogeneous diseases. Its application enables satisfactory presurgical differentiation of GB and metastases, a usual but difficult diagnostic challenge for neuroradiologist with vital implications in patient management.

Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology.

PURPOSE: Relationships between diffusion-weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange? METHODS: Monte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L / D0 , and clinical protocols (single diffusion encoding; maximum b-value: {1000, 1500, 2000} s/mm2 ; 5 unique gradient duration/separation pairs; SNR = { ∞ , 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion ( D ) and kurtosis ( K ) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell-size contrasts, and deployed on 9.4T ex vivo MRI-histology data of fixed mouse livers RESULTS: Relationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three-class cell-size classification task (for SNR = 20, bmax = 1500 s/mm2 , δ = 20 ms, and Δ = 75 ms). MRI detects cell-size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes. CONCLUSION: Salient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.

Spatial patterns of brain lesions assessed through covariance estimations of lesional voxels in multiple Sclerosis: The SPACE-MS technique.

Predicting disability in progressive multiple sclerosis (MS) is extremely challenging. Although there is some evidence that the spatial distribution of white matter (WM) lesions may play a role in disability accumulation, the lack of well-established quantitative metrics that characterise these aspects of MS pathology makes it difficult to assess their relevance for clinical progression. This study introduces a novel approach, called SPACE-MS, to quantitatively characterise spatial distributional features of brain MS lesions, so that these can be assessed as predictors of disability accumulation. In SPACE-MS, the covariance matrix of the spatial positions of each patient's lesional voxels is computed and its eigenvalues extracted. These are combined to derive rotationally-invariant metrics known to be common and robust descriptors of ellipsoid shape such as anisotropy, planarity and sphericity. Additionally, SPACE-MS metrics include a neuraxis caudality index, which we defined for the whole-brain lesion mask as well as for the most caudal brain lesion. These indicate how distant from the supplementary motor cortex (along the neuraxis) the whole-brain mask or the most caudal brain lesions are. We applied SPACE-MS to data from 515 patients involved in three studies: the MS-SMART (NCT01910259) and MS-STAT1 (NCT00647348) secondary progressive MS trials, and an observational study of primary and secondary progressive MS. Patients were assessed on motor and cognitive disability scales and underwent structural brain MRI (1.5/3.0 T), at baseline and after 2 years. The MRI protocol included 3DT1-weighted (1x1x1mm3) and 2DT2-weighted (1x1x3mm3) anatomical imaging. WM lesions were semiautomatically segmented on the T2-weighted scans, deriving whole-brain lesion masks. After co-registering the masks to the T1 images, SPACE-MS metrics were calculated and analysed through a series of multiple linear regression models, which were built to assess the ability of spatial distributional metrics to explain concurrent and future disability after adjusting for confounders. Patients whose WM lesions laid more caudally along the neuraxis or were more isotropically distributed in the brain (i.e. with whole-brain lesion masks displaying a high sphericity index) at baseline had greater motor and/or cognitive disability at baseline and over time, independently of brain lesion load and atrophy measures. In conclusion, here we introduced the SPACE-MS approach, which we showed is able to capture clinically relevant spatial distributional features of MS lesions independently of the sheer amount of lesions and brain tissue loss. Location of lesions in lower parts of the brain, where neurite density is particularly high, such as in the cerebellum and brainstem, and greater spatial spreading of lesions (i.e. more isotropic whole-brain lesion masks), possibly reflecting a higher number of WM tracts involved, are associated with clinical deterioration in progressive MS. The usefulness of the SPACE-MS approach, here demonstrated in MS, may be explored in other conditions also characterised by the presence of brain WM lesions.