Search

Purpose Characterizing vessel territories can provide crucial information for evaluation of cerebrovascular disorders. In this study, we present a novel postprocessing pipeline for vascular territorial imaging of cerebral arteries based on a noncontrast enhanced time-resolved 4D magnetic resonance angiography (MRA). Methods Eight healthy participants, 1 Moyamoya patient, and 1 arteriovenous malformations patient were recruited. Territorial segmentation and relative blood flow rate calculations of cerebral arteries including left and right middle cerebral arteries and left and right posterior cerebral arteries were carried out based on the 4D MRA-derived arterial arrival time maps of intracranial vessels. Results Among healthy young subjects, the average relative blood flow rate values corresponding to left and right middle cerebral arteries and left and right posterior cerebral arteries were 35.9 ± 5.9%, 32.9 ± 7.5%, 15.4 ± 3.8%, and 15.9 ± 2.5%, respectively. Excellent agreement was observed between relative blood flow rate values obtained from the proposed 4D MRA-based method and reference 2D phase contrast MRI. Abnormal cerebral circulations were visualized and quantified on both patients using the developed technique. Conclusion The vascular territorial imaging technique developed in this study allowed for the generation of territorial maps with user-defined level of details within a clinically feasible scan time, and as such may provide useful information to assess cerebral circulation balance in different pathologies.

Purpose Optic pathway gliomas (OPG) are low-grade pilocytic astrocytomas accounting for 3-5% of pediatric intracranial tumors. Accurate and quantitative follow-up of OPG using magnetic resonance imaging (MRI) is crucial for therapeutic decision making, yet is challenging due to the complex shape and heterogeneous tissue pattern which characterizes these tumors. The aim of this study was to implement automatic methods for segmentation and classification of OPG and its components, based on MRI. Methods A total of 202 MRI scans from 29 patients with chiasmatic OPG scanned longitudinally were retrospectively collected and included in this study. Data included T2 and post-contrast T1 weighted images. The entire tumor volume and its components were manually annotated by a senior neuro-radiologist, and inter- and intra-rater variability of the entire tumor volume was assessed in a subset of scans. Automatic tumor segmentation was performed using deep-learning method with U-Net+ResNet architecture. A fivefold cross-validation scheme was used to evaluate the automatic results relative to manual segmentation. Voxel-based classification of the tumor into enhanced, non-enhanced, and cystic components was performed using fuzzy c-means clustering. Results The results of the automatic tumor segmentation were: mean dice score = 0.736 ± 0.025, precision = 0.918 ± 0.014, and recall = 0.635 ± 0.039 for the validation data, and dice score = 0.761 ± 0.011, precision = 0.794 ± 0.028, and recall = 0.742 ± 0.012 for the test data. The accuracy of the voxel-based classification of tumor components was 0.94, with precision = 0.89, 0.97, and 0.85, and recall = 1.00, 0.79, and 0.94 for the non-enhanced, enhanced, and cystic components, respectively. Conclusion This study presents methods for automatic segmentation of chiasmatic OPG tumors and classification into the different components of the tumor, based on conventional MRI. Automatic quantitative longitudinal assessment of these tumors may improve radiological monitoring, facilitate early detection of disease progression and optimize therapy management.

BACKGROUND Differentiation between glioblastoma and brain metastasis is highly important due to differing medical treatment strategies. While MRI is the modality of choice for the assessment of patients with brain tumors, differentiation between glioblastoma and solitary brain metastasis may be challenging due to their similar appearance on MRI. PURPOSE To differentiate between glioblastoma and brain metastasis subtypes using radiomics analysis based on conventional post-contrast T1 -weighted (T1 W) MRI. STUDY TYPE Retrospective. SUBJECTS Data were acquired from 439 patients: 212 patients with glioblastoma and 227 patients with brain metastasis (breast, lung, and others). FIELD STRENGTH/SEQUENCE Post-contrast 3D T1 W gradient echo images, acquired with 1.5 and 3.0 T MR systems. ASSESSMENT Analysis included image preprocessing, segmentation of tumor area, and features extraction including: patients' clinical information, tumor location, first- and second-order statistical, morphological, wavelet features, and bag-of-features. Following dimension reduction, classification was performed using various machine-learning algorithms including support-vector machine (SVM), k-nearest neighbor, decision trees, and ensemble classifiers. STATISTICAL TESTS For classification, the data were divided into training (80%) and testing datasets (20%). Following optimization of the classifiers, mean sensitivity, specificity, accuracy, and area under the receiver operating characteristic curve (AUC) were calculated. RESULTS For the testing dataset, the best results for differentiation of glioblastoma from brain metastasis were obtained using the SVM classifier with mean accuracy = 0.85, sensitivity = 0.86, specificity = 0.85, and AUC = 0.96. The best classification results between glioblastoma and brain metastasis subtypes were obtained using SVM classifier with mean accuracy = 0.85, 0.89, 0.75, 0.90; sensitivity = 1.00, 0.60, 0.57, 0.11; specificity = 0.76, 0.92, 0.87, 0.99; and AUC = 0.98, 0.81, 0.83, 0.57 for the glioblastoma, breast, lung, and other brain metastases, respectively. DATA CONCLUSION Differentiation between glioblastoma and brain metastasis showed a high success rate based on postcontrast T1 W MRI. Classification between glioblastoma and brain metastasis subtypes may require additional MR sequences with other tissue contrasts. LEVEL OF EVIDENCE 1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:519-528.

BACKGROUND High-grade gliomas (HGGs) induce both vasogenic edema and extensive infiltration of tumor cells, both of which present with similar appearance on conventional MRI. Using current radiological criteria, differentiation between these tumoral and nontumoral areas within the nonenhancing lesion area remains challenging. PURPOSE To use radiomics patch-based analysis, based on conventional MRI, for the classification of the nonenhancing lesion area in patients with HGG into tumoral and nontumoral components. STUDY TYPE Prospective. SUBJECTS In all, 179 MRI scans were obtained from 102 patients: 67 patients with HGG and 35 patients with brain metastases. A subgroup of 15 patients with HGG were scanned before and following administration of bevacizumab. FIELD STRENGTH/SEQUENCE Pre and postcontrast agent T1 -weighted-imaging (WI), T2 WI, FLAIR, diffusion-tensor-imaging (DTI), and dynamic-contrast-enhanced (DCE)-MRI at 3T. ASSESSMENT A total of 225 histograms and gray-level-co-occurrence matrix-based features were extracted from the nonenhancing lesion area. Tumoral volumes of interest (VOIs) were defined at the peritumoral area in patients with HGG; nontumoral VOIs were defined in patients with brain metastasis. Twenty machine-learning algorithms including support-vector-machine (SVM), k-nearest neighbor, decision-trees, and ensemble classifiers were tested. The best classifier was trained on the entire labeled data, and was used to classify the entire data. STATISTICAL TESTS Dimensional reduction was performed on the 225 features using principal component analysis. Classification results were evaluated based on the sensitivity, specificity, and accuracy of each of the 20 classifiers, first based on a training and testing dataset (80% of the labeled data) in a 5-fold manner, and next by applying the best classifier to the validation data (the remaining 20% of the labeled data). Results were additionally evaluated by assessing differences in dynamic-contrast-enhanced plasma-volume (vp ) and volume-transfer-constant (ktrans ) values between the two components using Mann-Whitney U-test/t-test. RESULTS The best classification into tumoral and nontumoral lesion components was obtained using a linear SVM classifier, with average accuracy of 87%, sensitivity 86%, and specificity of 89% (for the training and testing data). Significantly higher vp and ktrans values (P < 0.0001) were detected in the tumoral compared to the nontumoral component. Preliminary classification results in a subgroup of patients treated with bevacizumab demonstrated a reduction mainly in the nontumoral component following administration of bevacizumab, enabling early assessment of disease progression in some patients. DATA CONCLUSION A radiomics patch-based analysis enables classification of the nonenhancing lesion area in patients with HGG. Preliminary results were promising and the proposed method has the potential to assist in clinical decision-making and to improve therapy response assessment in patients with HGG. LEVEL OF EVIDENCE 1 Technical Efficacy Stage 4 J. Magn. Reson. Imaging 2018.

Purpose In many clinical MRI scenarios, existing imaging information can be used to significantly shorten acquisition time or to improve Signal to Noise Ratio (SNR). In this paper the authors present a framework, referred to as FASTMER, for fast MRI by exploiting a reference image. Methods The proposed approach utilizes the possible similarity of the reference image to the acquired image, which exists in many clinical MRI imaging scenarios. Examples include similarity between adjacent slices in high resolution MRI, similarity between various contrasts in the same scan and similarity between different scans of the same patient. To account for the fact that the reference image may exhibit low similarity with the acquired image the authors develop an iterative weighted reconstruction approach, which tunes the weights according to the degree of similarity. Results Experimental results demonstrate the performance of the method in three different clinical MRI scenarios: The first example demonstrates SNR improvement in high resolution brain MRI, the second scenario exploits similarity between T2-weighted and fluid-attenuated inversion recovery (FLAIR) for fast FLAIR scanning and the last application utilizes similarity between baseline and follow-up scans for fast follow-up acquisition. The results show that FASTMER outperforms image reconstruction of existing state-of-the-art methods. Conclusions The authors present a framework for fast MRI by exploiting a reference image. Recovery is based on an iterative algorithm that supports cases in which similarity to the reference scan is not guaranteed. This extends the applicability of the FASTMER to different MRI scanning scenarios. Thanks to the existence of reference images in various clinical imaging tasks, the proposed framework can play a major role in improving reconstruction in many MR applications.

PURPOSE To optimize the analysis of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) under the two-compartment-exchange-model (2CXM) and to incorporate voxelwise bolus-arrival-time (BAT). MATERIALS AND METHODS The accuracy of the pharmacokinetic (PK) parameters, extracted from 3T DCE-MRI using 2CXM, was tested under several conditions: eight algorithms for data estimation; correction for BAT; using model selection; different temporal resolution and scan duration. Comparisons were performed on simulated data. The best algorithm was applied to seven patients with brain tumors or following stroke. The extracted perfusion parameters were compared to those of dynamic susceptibility contrast MRI (DSC-MRI). RESULTS ACoPeD (AIF-corrected-perfusion-DCE-MRI), an analysis using a 2nd derivative regularized-spline and incorporating BAT, achieved the most accurate estimation in simulated data, mean-relative-error: Fp , F, vp , ve : 24.8%, 41.7%, 26.4%, 27.2% vs. 76.5%, 190.8%, 78.8%, 82.39% of the direct four parameters estimation (one-sided two-sample t-test, P < 0.001). Correction for BAT increased the estimation accuracy of the PK parameters by more than 30% and provided a supertemporal resolution estimation of the BAT (higher than the acquired resolution, mean-absolute-error 0.2 sec). High temporal resolution (∼2 sec) is required to avoid biased estimation of PK parameters, and long scan duration (∼20 min) is important for reliable permeability but not for perfusion estimations, mean-error-reduction: E: ∼12%, ve : ∼6%. Using ACoPeD, PK values from normal-appearing white matter, gray matter, and lesion were extracted from patients. Preliminary results showed significant voxelwise correlations to DSC-MRI, between flow values in a patient following stroke (r = 0.49, P < 0.001), and blood volume in a patient with a brain tumor (r = 0.62, P < 0.001). CONCLUSION This study proposes an optimized analysis method, ACoPeD, for tissue perfusion and permeability estimation using DCE-MRI, to be used in clinical settings. LEVEL OF EVIDENCE 1 J. Magn. Reson. Imaging 2017;45:237-249.

We present a novel method for fiber-based comparison of diffusion tensor imaging (DTI) scans of groups of subjects. The method entails initial preprocessing and fiber reconstruction by tractography of each brain in its native coordinate system. Several diffusion parameters are sampled along each fiber and used in subsequent comparisons. A spatial correspondence between subjects is established based on geometric similarity between fibers in a template set (several choices for template are explored), and fibers in all other subjects. Diffusion parameters between groups are compared statistically for each template fiber. Results are presented at single fiber resolution. As an initial exploratory step in neurological population studies this method points to the locations affected by the pathology of interest, without requiring a hypothesis. It does not make any grouping assumptions on the fibers and no manual intervention is needed. The framework was applied here to 18 healthy subjects and 23 amyotrophic lateral sclerosis (ALS) patients. The results are compatible with previous findings and with the tract based spatial statistics (TBSS) method. Hum Brain Mapp 37:477-490, 2016. © 2015 Wiley Periodicals, Inc.

The hippocampus is known to play a vital role in learning and memory and was demonstrated as an early imaging marker for Alzheimer's disease (AD). However, its role as a predictor for mild cognitive impairment and dementia following stroke is unclear. The main purpose of this study was to examine the associations between hippocampal volume, mean diffusivity (MD) and connectivity and cognitive state following stroke. Eighty three consecutive first ever mild to moderate stroke or transient ischemic attack (TIA) survivors from our ongoing prospective TABASCO (Tel Aviv Brain Acute Stroke Cohort) study underwent magnetic resonance imaging scans within 7 days of stroke onset. Hippocampal volume was measured from T1 weighted images, hippocampal mean diffusivity was calculated from diffusion tensor imaging and connectivity was calculated from resting state fMRI. Global cognitive assessments were evaluated during hospitalization and 6 and 12 months later using a computerized neuropsychological battery. Multiple linear regression analysis was used to test which of the hippocampi measurements best predict cognitive state. All three imaging parameters were significantly correlated to each other (|r's| >0.3, P's

Purpose: Repeated brain MRI scans are performed in many clinical scenarios, such as follow up of patients with tumors and therapy response assessment. In this paper, the authors show an approach to utilize former scans of the patient for the acceleration of repeated MRI scans. Methods: The proposed approach utilizes the possible similarity of the repeated scans in longitudinal MRI studies. Since similarity is not guaranteed, sampling and reconstruction are adjusted during acquisition to match the actual similarity between the scans. The baseline MR scan is utilized both in the sampling stage, via adaptive sampling, and in the reconstruction stage, with weighted reconstruction. In adaptive sampling, k-space sampling locations are optimized during acquisition. Weighted reconstruction uses the locations of the nonzero coefficients in the sparse domains as a prior in the recovery process. The approach was tested on 2D and 3D MRI scans of patients with brain tumors. Results: The longitudinal adaptive compressed sensing MRI (LACS-MRI) scheme provides reconstruction quality which outperforms other CS-based approaches for rapid MRI. Examples are shown on patients with brain tumors and demonstrate improved spatial resolution. Compared with data sampled at the Nyquist rate, LACS-MRI exhibits signal-to-error ratio (SER) of 24.8 dB with undersampling factor of 16.6 in 3D MRI. Conclusions: The authors presented an adaptive method for image reconstruction utilizing similarity of scans in longitudinal MRI studies, where possible. The proposed approach can significantly reduce scanning time in many applications that consist of disease follow-up and monitoring of longitudinal changes in brain MRI.

Background Optic pathway gliomas (OPG) represent 5% of pediatric brain tumors and compose a major therapeutic dilemma to the treating physicians. While chemotherapy is widely used for these tumors, our ability to predict radiological response is still lacking. In this study, we use volumetric imaging to examine in detail the long-term effect of chemotherapy on the tumor as well as its various sub-components. Procedure The tumors of 15 patients with OPG, treated with chemotherapy, were longitudinally measured using our novel, previously described volumetric method. Patients were treated with up to five lines of chemotherapy. Sufficient follow-up imaging data, and patient's numbers, allowed for analysis of two treatment lines. Volumetric measurements of the tumors were segmented into solid-non-enhancing, solid-enhancing, and cystic components. Outcome analysis was done per specific treatment line and for the overall follow-up period. Results An average reduction of 9.7% (±23%) in the gross-total-solid volume (GTSV) was noted following treatment with vincristine and carboplatin. The cystic component grew under therapy by an average of 12.6% (±39%). When measured over the course of the whole study period, the cystic component grew by an average of 35% (±100%) and the GTSV increased by 12% (±35%). Conclusion Initial treatment with vincristine and carboplatin seems to have a minimal initial effect, mostly on the solid components. The cystic component in itself seems to be unaffected by chemotherapy, and contributes to the subsequent growth of the total volume. During the overall treatment period, both solid and cystic components grew regardless of combined treatment methods. Pediatr Blood Cancer 2015;62:1353–1359. © 2015 Wiley Periodicals, Inc.

Background: Patients with Parkinson’s disease have reduced gray matter volume and fractional anisotropy in both cortical and sub-cortical structures, yet changes in the pre-motor phase of the disease are unknown. Methods: A comprehensive imaging study using voxelbased morphometry and diffusion tensor imaging tract-based spatial statistics analysis was performed on 64 Ashkenazi Jewish asymptomatic first degree relatives of patients with Parkinson’s disease (30 mutation carriers), who carry the G2019S mutation in the leucine-rich repeat kinase 2 (LRRK2) gene. Results: No between-group differences in gray matter volume could be noted in either whole-brain or volume-of-interest analysis. Diffusion tensor imaging analysis did not identify group differences in white matter areas, and volume-of-interest analysis identified no differences in diffusivity parameters in Parkinson’s disease-related structures. Conclusions: G2019S carriers do not manifest changes in gray matter volume or diffusivity parameters in Parkinson’s disease-related structures prior to the appearance of motor symptoms. V C 2014 Interna

relaxationtime from a set of variable flip-angle (FA) spoiled gradientrecalled echo images.Materials and Methods: The proposed method: (a) usesa uniform weighting of all FAs, (b) takes into accountglobal inaccuracies in the generation of the prescribedFAs by estimating the actual FAs, and (c) incorporatesdata-driven local B

Purpose Stroke imaging studies during the acute phase are likely to precede several vascular brain mechanisms, which have an important role in patient outcome. The aim of this study was to identify within the lesion area during the subacute phase (≥1 day) reactive tissue, which may have the potential for recovery. Methods Twenty seven stroke patients from two cohorts were included. MRI performed during the subacute phase included conventional, perfusion and diffusion imaging. In cohort I, unsupervised multiparametric classification of the lesion area was performed. In cohort II threshold based classification was performed during the subacute phase, and radiological outcome was assessed at follow-up scan. Results Three tissue classes were identified in cohort I, referred to as irreversibly damaged, intermediary, and reactive tissue. Based on threshold values defined in cohort I, the reactive tissue was identified in 11/13 patients in cohort II, and showed tissue preservation/partial recovery in 9/11 patients at follow-up scan. The irreversibly damaged tissue was identified in 7/13 patients in cohort II, and predicted tissue necrosis in all cases. Conclusion Identification of reactive tissue following stroke during the subacute phase can improve radiological assessment, contribute to the understanding of brain recovery processes and has implications for new therapeutic approaches. Magn Reson Med 72:1381–1388, 2014. © 2013 Wiley Periodicals, Inc.

Purpose: To evaluate the sensitivity of high b value diffusion weight magnetic resonance imaging (DWI) in detecting normal white matter maturation, compare it to conventional diffusion tensor imaging (DTI), and to obtain normative quantitative data using this method. Materials and Methods: High b value DWI (bmax 6000 sec/mm 2 ) using q-space analysis and conventional DTI (b 1000 sec/mm 2 ) were performed on 36 healthy subjects aged 4 months to 23 years. Fractional-anisotropy (FA), apparent-displacement, and apparent-probability values were measured in all slices and in six regions of interest (ROIs) of large fiber tracks. Values were correlated with each other and with age using regression analysis. Results: FA, displacement, and probability indices from all slices were highly correlated with each other (r 0.87, P 0.0001) and with age (r 0.82, P 0.0001). All age-related changes in the six pre-determined ROIs were best fitted by mono-exponential functions. Changes in the splenium extended to a later age when compared with the genu of the corpus-callosum, while the centrum semi-ovale demonstrated the latest changes with age. Conclusions: High b-value DWI and DTI showed changes in white matter from infancy through adulthood. However, high b-value detects a signal that is likely to originate mainly from the intra-axonal water population, and thus may represent different aspects of development and different sensitivity to pathology.

Multiple sclerosis (MS) is an inflammatory disease of the central nervous system (CNS) which affects nearly one million people worldwide, leading to a progressive decline of motor and sensory functions, and permanent disability. High b-value diffusion-weighted MR images (b of up to 14000 s/mm 2 ) were acquired from the brains of controls and MS patients. These diffusion MR images, in which signal decay is not monoexponential, were analyzed using the q-space approach that emphasizes the diffusion characteristics of the slow-diffusing component. From this analysis, displacement and probability maps were constructed. The computed q-space analyzed MR images that were compared with conventional T1, T2 (fluid attenuated inversion recovery (FLAIR)), and diffusion tensor imaging (DTI) images were found to be sensitive to the pathophysiological state of white matter. The indices used to construct this qspace analyzed MR maps, provided a pronounced differentiation between normal tissue and tissues classified as MS plaques by the FLAIR images. More importantly, a pronounced differentiation was also observed between tissues classified by the FLAIR MR images as normal-appearing white matter (NAWM) in the MS brains, which are known to be abnormal, and the respective control tissues. The potential diagnostic capacity of high b-value diffusion q-space analyzed MR images is discussed, and experimental data that explains the consequences of using the q-space approach once the short pulse gradient approximation is violated are presented. Magn Reson Med 47:115‐126, 2002. © 2002 Wiley-Liss, Inc.

E. Kliper, E. Ben Assayag, S. Shenhar-Tsarfaty, L. Shopin, H. Hallevi, E. Uriel, A. Korczyn A, N. M. Bornstein, T. Hendler, O. Aizenstein, and D. Ben Bashat D The Wohl institute for Advanced Imaging, Tel Aviv Sourasky Medical Center, Tel Aviv, Israel, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel, Departments of Neurology, Tel Aviv Sourasky Medical Center, Israel, Israel, Departments of Radiology, Tel Aviv Sourasky Medical Center, Israel, Israel, The Wohl Institute for Advanced Imaging, Tel Aviv Sourasky Medical Center, Israel, Israel

Purpose: Tracking the progression of low grade tumors (LGTs) is a challenging task, due to their slow growth rate and associated complex internal tumor components, such as heterogeneous enhancement, hemorrhage, and cysts. In this paper, the authors show a semiautomatic method to reliably track the volume of LGTs and the evolution of their internal components in longitudinal MRI scans. Methods: The authors' method utilizes a spatiotemporal evolution modeling of the tumor and its internal components. Tumor components gray level parameters are estimated from the follow-up scan itself, obviating temporal normalization of gray levels. The tumor delineation procedure effectively incorporates internal classification of the baseline scan in the time-series as prior data to segment and classify a series of follow-up scans. The authors applied their method to 40 MRI scans of ten patients, acquired at two different institutions. Two types of LGTs were included: Optic pathway gliomas and thalamic astrocytomas. For each scan, a “gold standard” was obtained manually by experienced radiologists. The method is evaluated versus the gold standard with three measures: gross total volume error, total surface distance, and reliability of tracking tumor components evolution. Results: Compared to the gold standard the authors' method exhibits a mean Dice similarity volumetric measure of 86.58% and a mean surface distance error of 0.25 mm. In terms of its reliability in tracking the evolution of the internal components, the method exhibits strong positive correlation with the gold standard. Conclusions: The authors' method provides accurate and repeatable delineation of the tumor and its internal components, which is essential for therapy assessment of LGTs. Reliable tracking of internal tumor components over time is novel and potentially will be useful to streamline and improve follow-up of brain tumors, with indolent growth and behavior.