Multi-centre reproducibility of diffusion MRI parameters for clinical sequences in the brain

Inter-visit and inter-reader reproducibility of multi-parametric diffusion-weighted MR imaging in longitudinally imaged patients with metabolic dysfunction-associated fatty liver disease and healthy volunteers

BACKGROUND

Despite the widespread use of diffusion-weighted imaging (DWI) in metabolic dysfunction-associated fatty liver disease (MAFLD), MRI acquisition and quantification techniques vary in the literature suggesting the need for established and reproducible protocols. The goal of this study was to assess inter-visit and inter-reader reproducibility of DWI- and IVIM-derived parameters in patients with MAFLD and healthy volunteers using extensive sampling of the "fast" compartment, non-rigid registration, and exclusion voxels with poor fit quality.

METHODS

From June 2019 to April 2023, 31 subjects (20 patients with biopsy-proven MAFLD and 11 healthy volunteers) were included in this IRB-approved study. Subjects underwent MRI examinations twice within 40 days. 3.0 T DWI was acquired using a respiratory-triggered spin-echo diffusion-weighted echo-planar imaging sequence (b-values of 0, 10, 20, 30, 40, 50, 100, 200, 400, 800 s/mm2). DWI series were co-registered prior to voxel-wise non-linear regression of the IVIM model and voxels with poor fit quality were excluded (normalized root mean squared error ≥ 0.05). IVIM parameters (perfusion fraction, f; diffusion coefficient, D; and pseudo-diffusion coefficient, D*), and apparent diffusion coefficients (ADC) were computed from manual segmentation of the right liver lobe performed by two analysts on two MRI examinations.

RESULTS

All results are reported for f, D, D*, and ADC respectively. For inter-reader agreement on the first visit, ICC were of 0.985, 0.994, 0.986, and 0.993 respectively. For intra-reader agreement of analyst 1 assessed on both imaging examinations, ICC between visits were of 0.805, 0.759, 0.511, and 0.850 respectively. For inter-reader agreement on the first visit, mean bias and 95 % limits of agreement were (0.00 ± 0.03), (-0.01 ± 0.03) × 10-3 mm2/s, (0.70 ± 10.40) × 10-3 mm2/s, and (-0.02 ± 0.04) × 10-3 mm2/s respectively. For intra-reader agreement of analyst 1, mean bias and 95 % limits of agreement were (0.01 ± 0.09) × 10-3 mm2/s, (-0.01 ± 0.21) × 10-3 mm2/s, (-13.37 ± 56.19) × 10-3 mm2/s, and (-0.01 ± 0.16) × 10-3 mm2/s respectively. Except for parameter D* that was associated with between-subjects parameter variability (P = 0.009), there was no significant variability between subjects, examinations, or readers.

CONCLUSION

With our approach, IVIM parameters f, D, D*, and ADC provided excellent inter-reader agreement and good to very good inter-visit or intra-reader agreement, thus showing the reproducibility of IVIM-DWI of the liver in MAFLD patients and volunteers.

Machine Learning Analysis in Diffusion Kurtosis Imaging for Discriminating Pediatric Posterior Fossa Tumors: A Repeatability and Accuracy Pilot Study

Background and purpose: Differentiating pediatric posterior fossa (PF) tumors such as medulloblastoma (MB), ependymoma (EP), and pilocytic astrocytoma (PA) remains relevant, because of important treatment and prognostic implications. Diffusion kurtosis imaging (DKI) has not yet been investigated for discrimination of pediatric PF tumors. Estimating diffusion values from whole-tumor-based (VOI) segmentations may improve diffusion measurement repeatability compared to conventional region-of-interest (ROI) approaches. Our purpose was to compare repeatability between ROI and VOI DKI-derived diffusion measurements and assess DKI accuracy in discriminating among pediatric PF tumors. Materials and methods: We retrospectively analyzed 34 children (M, F, mean age 7.48 years) with PF tumors who underwent preoperative examination on a 3 Tesla magnet, including DKI. For each patient, two neuroradiologists independently segmented the whole solid tumor, the ROI of the area of maximum tumor diameter, and a small 5 mm ROI. The automated analysis pipeline included inter-observer variability, statistical, and machine learning (ML) analyses. We evaluated inter-observer variability with coefficient of variation (COV) and Bland-Altman plots. We estimated DKI metrics accuracy in discriminating among tumor histology with MANOVA analysis. In order to account for class imbalances, we applied SMOTE to balance the dataset. Finally, we performed a Random Forest (RF) machine learning classification analysis based on all DKI metrics from the SMOTE dataset by partitioning 70/30 the training and testing cohort. Results: Tumor histology included medulloblastoma (15), pilocytic astrocytoma (14), and ependymoma (5). VOI-based measurements presented lower variability than ROI-based measurements across all DKI metrics and were used for the analysis. DKI-derived metrics could accurately discriminate between tumor subtypes (Pillai's trace: p < 0.001). SMOTE generated 11 synthetic observations (10 EP and 1 PA), resulting in a balanced dataset with 45 instances (34 original and 11 synthetic). ML analysis yielded an accuracy of 0.928, which correctly predicted all but one lesion in the testing set. Conclusions: VOI-based measurements presented improved repeatability compared to ROI-based measurements across all diffusion metrics. An ML classification algorithm resulted accurate in discriminating PF tumors on a SMOTE-generated dataset. ML techniques based on DKI-derived metrics are useful for the discrimination of pediatric PF tumors.

Tractometry of Human Visual White Matter Pathways in Health and Disease

Diffusion-weighted MRI (dMRI) provides a unique non-invasive view of human brain tissue properties. The present review article focuses on tractometry analysis methods that use dMRI to assess the properties of brain tissue within the long-range connections comprising brain networks. We focus specifically on the major white matter tracts that convey visual information. These connections are particularly important because vision provides rich information from the environment that supports a large range of daily life activities. Many of the diseases of the visual system are associated with advanced aging, and tractometry of the visual system is particularly important in the modern aging society. We provide an overview of the tractometry analysis pipeline, which includes a primer on dMRI data acquisition, voxelwise model fitting, tractography, recognition of white matter tracts, and calculation of tract tissue property profiles. We then review dMRI-based methods for analyzing visual white matter tracts: the optic nerve, optic tract, optic radiation, forceps major, and vertical occipital fasciculus. For each tract, we review background anatomical knowledge together with recent findings in tractometry studies on these tracts and their properties in relation to visual function and disease. Overall, we find that measurements of the brain's visual white matter are sensitive to a range of disorders and correlate with perceptual abilities. We highlight new and promising analysis methods, as well as some of the current barriers to progress toward integration of these methods into clinical practice. These barriers, such as variability in measurements between protocols and instruments, are targets for future development.

Dynamic Evolution of Infarct Volumes at MRI in Ischemic Stroke Due to Large Vessel Occlusion

BACKGROUND AND OBJECTIVES

The typical infarct volume trajectories in stroke patients, categorized as slow or fast progressors, remain largely unknown. This study aimed to reveal the characteristic spatiotemporal evolutions of infarct volumes caused by large vessel occlusion (LVO) and show that such growth charts help anticipate clinical outcomes.

METHODS

We conducted a secondary analysis from prospectively collected databases (FRAME, 2017-2019; ETIS, 2015-2022). We selected acute MRI data from anterior LVO stroke patients with witnessed onset, which were divided into training and independent validation datasets. In the training dataset, using Gaussian mixture analysis, we classified the patients into 3 growth groups based on their rate of infarct growth (diffusion volume/time-to-imaging). Subsequently, we extrapolated pseudo-longitudinal models of infarct growth for each group and generated sequential frequency maps to highlight the spatial distribution of infarct growth. We used these charts to attribute a growth group to the independent patients from the validation dataset. We compared their 3-month modified Rankin scale (mRS) with the predicted values based on a multivariable regression model from the training dataset that used growth group as an independent variable.

RESULTS

We included 804 patients (median age 73.0 years [interquartile range 61.2-82.0 years]; 409 men). The training dataset revealed nonsupervised clustering into 11% (74/703) slow, 62% (437/703) intermediate, and 27% (192/703) fast progressors. Infarct volume evolutions were best fitted with a linear (r = 0.809; p < 0.001), cubic (r = 0.471; p < 0.001), and power (r = 0.63; p < 0.001) function for the slow, intermediate, and fast progressors, respectively. Notably, the deep nuclei and insular cortex were rapidly affected in the intermediate and fast groups with further cortical involvement in the fast group. The variable growth group significantly predicted the 3-month mRS (multivariate odds ratio 0.51; 95% CI 0.37-0.72, p < 0.0001) in the training dataset, yielding a mean area under the receiver operating characteristic curve of 0.78 (95% CI 0.66-0.88) in the independent validation dataset.

DISCUSSION

We revealed spatiotemporal archetype dynamic evolutions following LVO stroke according to 3 growth phenotypes called slow, intermediate, and fast progressors, providing insight into anticipating clinical outcome. We expect this could help in designing neuroprotective trials aiming at modulating infarct growth before EVT.

Implementing ABCD studyⓇ MRI sequences for multi-site cohort studies: Practical guide to necessary steps, preprocessing methods, and challenges

Large multi-site studies that combine magnetic resonance imaging (MRI) data across research sites present exceptional opportunities to advance neuroscience research. However, scanner or site variability and non-standardised image acquisition protocols, data processing and analysis pipelines can adversely affect the reliability and repeatability of MRI derived brain measures. We implemented a standardised MRI protocol based on that used in the Adolescent Brain Cognition Development (ABCD)Ⓡ study in two sites, and across four MRI scanners. Twice repeated measurements of a single healthy volunteer were obtained in two sites and in four 3T MRI scanners (vendors: Siemens, Philips, and GE). Imaging data included anatomical scans (T1 weighted, T2 weighted), diffusion weighted imaging (DWI) and resting state functional MRI (rs-fMRI). Standardised containerized pipelines were utilised to pre-process the data and different image quality metrics and test-retest variability of different brain metrics were evaluated. The implementation of the MRI protocols was possible with minor adjustments in acquisition (e.g. repetition time (TR), higher b-values) and exporting (DICOM formats) of images due to different technical performance of the scanners. This study provides practical insights into the implementation of standardised sequences and data processing for multisite studies, showcase the benefits of containerised preprocessing tools, and highlights the need for careful optimisation of multisite image acquisition.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

Quantitative analysis of apparent diffusion coefficients to predict neurological prognosis in cardiac arrest survivors: an observational derivation and internal-external validation study

BACKGROUND

This study aimed to validate apparent diffusion coefficient (ADC) values and thresholds to predict poor neurological outcomes in out-of-hospital cardiac arrest (OHCA) survivors by quantitatively analysing the ADC values via brain magnetic resonance imaging (MRI).

METHODS

This observational study used prospectively collected data from two tertiary academic hospitals. The derivation cohort comprised 70% of the patients randomly selected from one hospital, whereas the internal validation cohort comprised the remaining 30%. The external validation cohort used the data from another hospital, and the MRI data were restricted to scans conducted at 3 T within 72-96 h after an OHCA experience. We analysed the percentage of brain volume below a specific ADC value at 50-step intervals ranging from 200 to 1200 × 10-6 mm2/s, identifying thresholds that differentiate between good and poor outcomes. Poor neurological outcomes were defined as cerebral performance categories 3-5, 6 months after experiencing an OHCA.

RESULTS

A total of 448 brain MRI scans were evaluated, including a derivation cohort (n = 224) and internal/external validation cohorts (n = 96/128, respectively). The proportion of brain volume with ADC values below 450, 500, 550, 600, and 650 × 10-6 mm2/s demonstrated good to excellent performance in predicting poor neurological outcomes in the derivation group (area under the curve [AUC] 0.89-0.91), and there were no statistically significant differences in performances among the derivation, internal validation, and external validation groups (all P > 0.5). Among these, the proportion of brain volume with an ADC below 600 × 10-6 mm2/s predicted a poor outcome with a 0% false-positive rate (FPR) and 76% (95% confidence interval [CI] 68-83) sensitivity at a threshold of > 13.2% in the derivation cohort. In both the internal and external validation cohorts, when using the same threshold, a specificity of 100% corresponded to sensitivities of 71% (95% CI 58-81) and 78% (95% CI 66-87), respectively.

CONCLUSIONS

In this validation study, by consistently restricting the MRI types and timing during quantitative analysis of ADC values in brain MRI, we observed high reproducibility and sensitivity at a 0% FPR. Prospective multicentre studies are necessary to validate these findings.

SARS-CoV2 evokes structural brain changes resulting in declined executive function

BACKGROUND

Several research has underlined the multi-system character of COVID-19. Though effects on the Central Nervous System are mainly discussed as disease-specific affections due to the virus' neurotropism, no comprehensive disease model of COVID-19 exists on a neurofunctional base by now. We aimed to investigate neuroplastic grey- and white matter changes related to COVID-19 and to link these changes to neurocognitive testings leading towards a multi-dimensional disease model.

METHODS

Groups of acutely ill COVID-19 patients (n = 16), recovered COVID-19 patients (n = 21) and healthy controls (n = 13) were prospectively included into this study. MR-imaging included T1-weighted sequences for analysis of grey matter using voxel-based morphometry and diffusion-weighted sequences to investigate white matter tracts using probabilistic tractography. Comprehensive neurocognitive testing for verbal and non-verbal domains was performed.

RESULTS

Alterations strongly focused on grey matter of the frontal-basal ganglia-thalamus network and temporal areas, as well as fiber tracts connecting these areas. In acute COVID-19 patients, a decline of grey matter volume was found with an accompanying diminution of white matter tracts. A decline in executive function and especially verbal fluency was found in acute patients, partially persisting in recovered.

CONCLUSION

Changes in gray matter volume and white matter tracts included mainly areas involved in networks of executive control and language. Deeper understanding of these alterations is necessary especially with respect to long-term impairments, often referred to as 'Post-COVID'.

Repeatability of diffusion-based stiffness prediction - A healthy volunteer study

INTRODUCTION

The study investigated the repeatability of brain diffusion-based stiffness prediction (DWIstiff) in healthy volunteers.

METHODS

Thirty-one healthy volunteers were examined with DWIstiff using two different sets of b-values: b200-1500 s/mm2 (DWIstiff, 1500) and b200-1000 s/mm2 (DWIstiff, 1000). Each b-value set was scanned twice per imaging session without repositioning the participants. DWIstiff images were reconstructed from each set. Two observers delineated regions of interest (ROIs) on each DWIstiff image. The repeatability coefficient (RC), coefficient of variation (CV), inter- and intraobserver agreement were calculated.

RESULTS

After excluding three participants due to image artifacts, the study included twenty-eight volunteers (mean age (range)) 37 years (24-62), 10 males, 18 females). For DWIstiff, 1500, the lowest and the highest RCs were in the parietal lobe (0.52) and respectively the brain stem (1.17). The lowest RC for DWIstiff, 1000 was in the frontal lobe (0.42) and the highest in the brain stem (1.58). The CV for whole brain measurements was 3.83 % for DWIstiff, 1500 and 4.93 % for DWIstiff, 1000. The Bland‒Altman (BA) limits of agreement (LoA) for the intraobserver agreement of DWIstiff, 1500 were -0.90 to 1.06 and respectively -0.78 to 0.88 for DWIstiff, 1000. Regarding interobserver agreement, the LoA were -0.85 to 0.94 for DWIstiff, 1500 and -0.61 to 0.66 for DWIstiff, 1000.

CONCLUSION

DWIstiff is a precise technique with some observer dependence. Repeatability is higher for DWIstiff, 1000 s/mm2 than for DWIstiff 1500 s/mm2.

IMPLICATIONS FOR PRACTICE

Our findings suggest that DWIstiff can reliably detect stiffness changes larger than 4.93 % in healthy volunteers. Further studies should investigate whether the repeatability of DWIstiff may be affected by the presence of pathology such as a brain tumor.

Incremental value of volumetric multiparametric MRI over Fudan score for prognosis of unresectable intrahepatic cholangiocarcinoma treated with systemic chemotherapy

BACKGROUND

Individualized patient care requires prognostic models customized to a tumor and an individual's disease profile for reliable survival prediction. MRI has prognostic value for intrahepatic cholangiocarcinoma (ICCA). Existing prognostic models for ICCA exclude imaging-based information about an individual's tumor that may reflect important aspects of tumor's biology. Fudan score, a prognostic model applicable to unresectable ICCA, is limited by subjective morphologic imaging parameters.

OBJECTIVES

To assess the prognostic value of baseline volumetric multiparametric MRI in unresectable intrahepatic cholangiocarcinoma (ICCA) treated with systemic chemotherapy and the incremental value of MRI over the Fudan score.

METHODS

This retrospective study included 114 ICCA patients treated with systemic chemotherapy between 2007 and 2021 after a baseline MRI. The single largest tumor was volumetrically assessed for anatomic (total tumor volume and diameter) and functional parameters (viable tumor volume, percentage-viable tumor volume, viable tumor burden, and ADC). A derivation cohort of 30 patients was utilized to identify MRI parameters associated with overall survival (OS) using Cox regression analysis. The incremental value of MRI over Fudan score was assessed on an independent sub-cohort of 84 patients using Kaplan-Meier analysis and C-index.

RESULTS

114 patients (64 years +/- 11; 61 women) were evaluated. Pre-treatment high (>1350x10-6 mm2/sec) ADC was the only independent predictor of OS (HR, 8.07; P < 0.001). Replacing subjective tumor boundary with objective ADC value, and using modified biochemical thresholds increased the prognostic stratification for the risk groups in the modified ADC-Fudan model compared to the original Fudan model (median survival 12 and 4.5 months; P = 0.055; vs. 11 and 3 months; P < 0.001). The modified ADC-Fudan model demonstrated an 11 % improvement over the original Fudan model (c-index: 0.80 vs. 0.69; P = 0.044) for survival prediction.

CONCLUSIONS

High pre-treatment volumetric ADC was associated with unfavorable prognosis in patients with unresectable intrahepatic cholangiocarcinoma treated with systemic chemotherapy. Supplementing the original Fudan model with ADC and modified serum marker thresholds improved the survival prediction performance by 11% in the resulting modified ADC-Fudan model.

CLINICAL IMPACT

Volumetric MRI could improve the survival prediction among ICCA patients prior to receiving potentially toxic and expensive palliative chemotherapies. This could potentially guide individualized therapy for this patient cohort.

Repeatability and reproducibility of apparent exchange rate measurements in yeast cell phantoms using filter-exchange imaging

OBJECTIVES

Development of a protocol for validation and quality assurance of filter-exchange imaging (FEXI) pulse sequences with well-defined and reproducible phantoms.

MATERIALS AND METHODS

A FEXI pulse sequence was implemented on a 7 T preclinical MRI scanner. Six experiments in three different test categories were established for sequence validation, demonstration of the reproducibility of phantoms and the measurement of induced changes in the apparent exchange rate (AXR). First, an ice-water phantom was used to investigate the consistency of apparent diffusion coefficient (ADC) measurements with different diffusion filters. Second, yeast cell phantoms were utilized to validate the determination of the AXR in terms of repeatability (same phantom and session), reproducibility (separate but comparable phantoms in different sessions) and directionality of diffusion encodings. Third, the yeast cell phantoms were, furthermore, used to assess potential AXR bias because of altered cell density and temperature. In addition, a treatment experiment with aquaporin inhibitors was performed to evaluate the influence of these compounds on the cell membrane permeability in yeast cells.

RESULTS

FEXI-based ADC measurements of an ice-water phantom were performed for three different filter strengths, showed good agreement with the literature value of 1.099 × 10-3 mm2/s and had a maximum coefficient of variation (CV) of 0.55% within the individual filter strengths. AXR estimation in a single yeast cell phantom and imaging session with five repetitions resulted in an overall mean value of (1.49 ± 0.05) s-1 and a CV of 3.4% between the chosen regions of interest. For three separately prepared phantoms, AXR measurements resulted in a mean value of (1.50 ± 0.04) s-1 and a CV of 2.7% across the three phantoms, demonstrating high reproducibility. Across three orthogonal diffusion directions, a mean value of (1.57 ± 0.03) s-1 with a CV of 1.9% was detected, consistent with isotropy of AXR in yeast cells. Temperature and AXR were linearly correlated (R2 = 0.99) and an activation energy EA of 37.7 kJ/mol was determined by Arrhenius plot. Furthermore, a negative correlation was found between cell density (as determined by the reference ADC/fe) and AXR (R2 = 0.95). The treatment experiment resulted in significantly decreased AXR values at different temperatures in the treated sample compared to the untreated control indicating an inhibiting effect.

CONCLUSIONS

Using ice-water and yeast cell-based phantoms, a protocol for the validation of FEXI pulse sequences was established for the assessment of stability, repeatability, reproducibility and directionality. In addition, a strong dependence of AXR on cell density and temperature was shown. As AXR is an emerging novel imaging biomarker, the suggested protocol will be useful for quality assurance of AXR measurements within a study and potentially across multiple sites.

Traveling Subject-Informed Harmonization Increases Reliability of Brain Diffusion Tensor and Neurite Mapping

Diffusion-weighted magnetic resonance imaging (dMRI) of brain has helped elucidate the microstructural changes of psychiatric and neurodegenerative disorders. Inconsistency between MRI models has hampered clinical application of dMRI-based metrics. Using harmonized dMRI data of 300 scans from 69 traveling subjects (TS) scanning the same individuals at multiple conditions with 13 MRI models and 2 protocols, the widely-used metrics such as diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) were evaluated before and after harmonization with a combined association test (ComBat) or TS-based general linear model (TS-GLM). Results showed that both ComBat and TS-GLM significantly reduced the effects of the MRI site, model, and protocol for diffusion metrics while maintaining the intersubject biological effects. The harmonization power of TS-GLM based on TS data model is more powerful than that of ComBat. In conclusion, our research demonstrated that although ComBat and TS-GLM harmonization approaches were effective at reducing the scanner effects of the site, model, and protocol for DTI and NODDI metrics in WM, they exhibited high retainability of biological effects. Therefore, we suggest that, after harmonizing DTI and NODDI metrics, a multisite study with large cohorts can accurately detect small pathological changes by retaining pathological effects.

Diffusion-weighted imaging on an MRI-linear accelerator to identify adversely prognostic tumour regions in glioblastoma during chemoradiation

BACKGROUND AND PURPOSE

Survival in glioblastoma might be extended by escalating the radiotherapy dose to treatment-resistant tumour and adapting to tumour changes. Diffusion-weighted imaging (DWI) on MRI-linear accelerators (MR-Linacs) could be used to identify a dose escalation target, but its prognostic value must be demonstrated. The purpose of this study was to determine whether MR-Linac DWI can assess treatment response in glioblastoma and whether changes in DWI show greater prognostic value than changes in the contrast-enhancing gross tumour volume (GTV).

MATERIALS AND METHODS

Seventy-five patients with glioblastoma were treated with chemoradiotherapy, of which 32 were treated on a 1.5 T MRI-linear accelerator (MR-Linac). Patients were imaged with simulation MRI scanners (MR-sim) at treatment planning and weeks 2, 4, and 10 after treatment start. Twenty-eight patients had additional MR-Linac DWI sequences. Cox modelling was used to evaluate the correlation of overall and progression-free survival (OS and PFS) with clinical variables and volumetric changes in the GTV and low-ADC regions (ADC < 1.25 µm2/ms within GTV).

RESULTS

In total, 479 MR-Linac DWI and 289 MR-sim DWI datasets were analyzed. MR-Linac low-ADC changes between weeks 2 and 5 inclusive were prognostic for OS (hazard ratio lower limits ≥ 1.2, p-values ≤ 0.02). MR-sim low-ADC changes showed greater correlation with OS and PFS than GTV changes (e.g., OS hazard ratio at week 2 was 3.4 (p <0.001) for low-ADC versus 2.0 (p = 0.022) for GTV).

CONCLUSION

MR-Linac DWI can measure low-ADC tumour volumes that correlate with OS and PFS better than contrast-enhancing GTV. Low-ADC regions could serve as dose escalation targets.

Toward vendor-independent measurement of cerebral venous oxygenation: Comparison of TRUST MRI across three major MRI manufacturers and association with end-tidal CO2

Cerebral venous oxygenation (Yv ) is a valuable biomarker for a variety of brain diseases. T2 relaxation under spin tagging (TRUST) MRI is a widely used method for Yv quantification. In this work, there were two main objectives. The first was to evaluate the reproducibility of TRUST Yv measurements across MRI scanners from different vendors. The second was to examine the correlation between Yv and end-tidal CO2 (EtCO2 ) in a multisite, multivendor setting and determine the usefulness of this correlation to account for variations in Yv caused by normal variations and physiological fluctuations. Standardized TRUST pulse sequences were implemented on three scanners from major MRI vendors (GE, Siemens, Philips). These scanners were located at two research institutions. Ten healthy subjects were scanned. On each scanner, the subject underwent two scan sessions, each of which included three TRUST scans, to evaluate the intrasession and intersession reproducibility of Yv . Each scanner was also equipped with a capnograph device to record the EtCO2 of the subject during the MRI scan. We found no significant bias in Yv measurements across the three scanners (P = 0.18). The measured Yv values on the three scanners were also strongly correlated with each other (intraclass correlation coefficients > 0.85, P < 0.001). The intrasession and intersession coefficients of variation of Yv were less than 4% and showed no significant difference among the scanners. In addition, our results revealed that (1) within the same subject, Yv increased with EtCO2 at a rate of 1.24 ± 0.17%/mmHg (P < 0.0001), and (2) across different subjects, individuals with a higher EtCO2 had a higher Yv , at a rate of 0.94 ± 0.36%/mmHg (P = 0.01). These results suggest that (1) the standardized TRUST sequences had similar accuracies and reproducibilities for the quantification of Yv across the scanners, and (2) recording of EtCO2 may be a useful complement to Yv measurement to account for CO2 -related physiological fluctuations in Yv in multisite, multivendor studies.

Immune checkpoint inhibitor treatment of brain metastasis associated with a less invasive growth pattern, higher T-cell infiltration and raised tumor ADC on diffusion weighted MRI

BACKGROUND

Brain metastases are the most common intracranial tumors with an increasing incidence. They are an important cause of morbidity and mortality in patients with solid organ cancer and a focus of recent clinical research and experimental interest. Immune checkpoint inhibitors are being increasingly used to treat solid organ cancers.

METHODS

To determine whether immune checkpoint inhibitors were biologically effective in the brain, we compared melanoma brain metastasis samples where treatment with ipilimumab had occurred preoperatively to those who had not received any immune modulating therapy and looked for histopathological (invasion, vascularity, metastasis inducing proteins, matrix metalloproteinases, immune cell infiltration, tissue architecture) and advanced MRI differences (diffusion weighted imaging).

RESULTS

Co-localized tissue samples from the same regions as MRI regions of interest showed significantly lower vascularity (density of CD34 + vessels) in the core and higher T-cell infiltration (CD3 + cells) in the leading edge for ipilimumab-treated brain metastasis samples than for untreated cases and this correlated with a higher tumor ADC signal at post-treatment/preoperative MRI brain.

CONCLUSIONS

Treatment of a melanoma brain metastasis with ipilimumab appears to cause measurable biological changes in the tumor that can be correlated with post-treatment diffusion weighted MRI imaging, suggesting both a mechanism of action and a possible surrogate marker of efficacy.

Diffusion-Weighted Imaging of the Orbit: A Case Series and Systematic Review

PURPOSE

To describe the findings of diffusion-weighted imaging (DWI) for a series of orbital lesions and provide a systematic review of relevant literature.

METHODS

A retrospective review of 20 patients with orbital lesions who underwent MRI with DWI at two academic institutions between 2015 and 2020 was performed. Lesion diagnosis was histopathologically confirmed except a presumed cavernous hemangioma. Echoplanar diffusion-weighted images had been acquired using 2 or 3 b values (b=0 and 1000 or b=0, 500, and 1000) at 1.5T or 3T. Lesions with significant artifacts were excluded. DWI sequences were analyzed by neuro-radiologists blinded to the diagnosis. Mean ADC values of lesions were calculated from a single region of interest. An independent two-tailed t test was used to compare categories of lesions with p < 0.05 considered significant. A systematic review of the literature was performed.

RESULTS

Our study included 21 lesions. ADC values were significantly lower for malignant lesions (0.628 ± 0.125 × 10 -3 mm 2 /s) than inflammatory lesions (1.167 ± 0.381 × 10 -3 mm 2 /s) ( p < 0.001). ADC values were significantly lower for orbital lymphoma (mean 0.621 ± 0.147 × 10 -3 mm 2 /s) than idiopathic orbital inflammation (mean 1.188 ± 0.269 × 10 -3 mm 2 /s) with no overlap ( p < 0.001).

CONCLUSIONS

Orbital malignancies demonstrated lower ADC values, while inflammatory processes demonstrated higher ADC values, except IgG4-related disease. DWI and ADC values differentiated idiopathic orbital inflammation from orbital lymphoma. This study highlights the role of DWI in evaluating orbital pathology.

Diffusion MRI-based connectomics features improve the noninvasive prediction of H3K27M mutation in brainstem gliomas

PURPOSE

To establish an individualized predictive model to identify patients with brainstem gliomas (BSGs) at high risk of H3K27M mutation, with the inclusion of brain structural connectivity analysis based on diffusion MRI (dMRI).

MATERIALS AND METHODS

A primary cohort of 133 patients with BSGs (80 H3K27M-mutant) were retrospectively included. All patients underwent preoperative conventional MRI and dMRI. Tumor radiomics features were extracted from conventional MRI, while two kinds of global connectomics features were extracted from dMRI. A machine learning-based individualized H3K27M mutation prediction model combining radiomics and connectomics features was generated with a nested cross validation strategy. Relief algorithm and SVM method were used in each outer LOOCV loop to select the most robust and discriminative features. Additionally, two predictive signatures were established using the LASSO method, and simplified logistic models were built using multivariable logistic regression analysis. An independent cohort of 27 patients was used to validate the best model.

RESULTS

35 tumor-related radiomics features, 51 topological properties of brain structural connectivity networks, and 11 microstructural measures along white matter tracts were selected to construct a machine learning-based H3K27M mutation prediction model, which achieved an AUC of 0.9136 in the independent validation set. Radiomics- and connectomics-based signatures were generated and simplified combined logistic model was built, upon which derived nomograph achieved an AUC of 0.8827 in the validation cohort.

CONCLUSION

dMRI is valuable in predicting H3K27M mutation in BSGs, and connectomics analysis is a promising approach. Combining multiple MRI sequences and clinical features, the established models have good performance.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

Direct jet coaxial electrospinning of axon-mimicking fibers for diffusion tensor imaging

Hollow polymer microfibers with variable microstructural and hydrophilic properties were proposed as building elements to create axon-mimicking phantoms for validation of diffusion tensor imaging (DTI). The axon-mimicking microfibers were fabricated in a mm-thick 3D anisotropic fiber strip, by direct jet coaxial electrospinning of PCL/polysiloxane-based surfactant (PSi) mixture as shell and polyethylene oxide (PEO) as core. Hydrophilic PCL-PSi fiber strips were first obtained by carefully selecting appropriate solvents for the core and appropriate fiber collector rotating and transverse speeds. The porous cross-section and anisotropic orientation of axon-mimicking fibers were then quantitatively evaluated using two ImageJ plugins-nearest distance (ND) and directionality based on their scanning electron microscopy (SEM) images. Third, axon-mimicking phantom was constructed from PCL-PSi fiber strips with variable porous-section and fiber orientation and tested on a 3T clinical MR scanner. The relationship between DTI measurements (mean diffusivity [MD] and fractional anisotropy [FA]) of phantom samples and their pore size and fiber orientation was investigated. Two key microstructural parameters of axon-mimicking phantoms including normalized pore distance and dispersion of fiber orientation could well interpret the variations in DTI measurements. Two PCL-PSi phantom samples made from different regions of the same fiber strips were found to have similar MD and FA values, indicating that the direct jet coaxial electrospun fiber strips had consistent microstructure. More importantly, the MD and FA values of the developed axon-mimicking phantoms were mostly in the biologically relevant range.

A systematic review of the relationship between magnetic resonance imaging based resting-state and structural networks in the rodent brain

Recent developments in rodent brain imaging have enabled translational characterization of functional and structural connectivity at the whole brain level in vivo. Nevertheless, fundamental questions about the link between structural and functional networks remain unsolved. In this review, we systematically searched for experimental studies in rodents investigating both structural and functional network measures, including studies correlating functional connectivity using resting-state functional MRI with diffusion tensor imaging or viral tracing data. We aimed to answer whether functional networks reflect the architecture of the structural connectome, how this reciprocal relationship changes throughout a disease, how structural and functional changes relate to each other, and whether changes follow the same timeline. We present the knowledge derived exclusively from studies that included in vivo imaging of functional and structural networks. The limited number of available reports makes it difficult to draw general conclusions besides finding a spatial and temporal decoupling between structural and functional networks during brain disease. Data suggest that when overcoming the currently limited evidence through future studies with combined imaging in various disease models, it will be possible to explore the interaction between both network systems as a disease or recovery biomarker.

Early MRI Predictors of Relapse in Primary Central Nervous System Lymphoma Treated with MATRix Immunochemotherapy

Primary Central Nervous System Lymphoma (PCNSL) is a highly malignant brain tumour. We investigated dynamic changes in tumour volume and apparent diffusion coefficient (ADC) measurements for predicting outcome following treatment with MATRix chemotherapy in PCNSL. Patients treated with MATRix (n = 38) underwent T1 contrast-enhanced (T1CE) and diffusion-weighted imaging (DWI) before treatment, after two cycles and after four cycles of chemotherapy. Response was assessed using the International PCNSL Collaborative Group (IPCG) imaging criteria. ADC histogram parameters and T1CE tumour volumes were compared among response groups, using one-way ANOVA testing. Logistic regression was performed to examine those imaging parameters predictive of response. Response after two cycles of chemotherapy differed from response after four cycles; of the six patients with progressive disease (PD) after four cycles of treatment, two (33%) had demonstrated a partial response (PR) or complete response (CR) after two cycles. ADC_mean at baseline, T1CE at baseline and T1CE percentage volume change differed between response groups (0.005 < p < 0.038) and were predictive of MATRix treatment response (area under the curve: 0.672-0.854). Baseline ADC and T1CE metrics are potential biomarkers for risk stratification of PCNSL patients early during remission induction therapy with MATRix. Standard interim response assessment (after two cycles) according to IPCG imaging criteria does not reliably predict early disease progression in the context of a conventional treatment approach.

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

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

White matter microstructural perturbations after total sleep deprivation in depression

Background

Total sleep deprivation (TSD) transiently reverses depressive symptoms in a majority of patients with depression. How TSD modulates diffusion tensor imaging (DTI) measures of white matter (WM) microstructure, which may be linked with TSD's rapid antidepressant effects, remains uncharacterized.

Methods

Patients with depression (N = 48, mean age = 33, 26 women) completed diffusion-weighted imaging and Hamilton Depression Rating (HDRS) and rumination scales before and after >24 h of TSD. Healthy controls (HC) (N = 53, 23 women) completed the same assessments at baseline, and after receiving TSD in a subset of HCs (N = 15). Tract based spatial statistics (TBSS) investigated voxelwise changes in fractional anisotropy (FA) across major WM pathways pre-to-post TSD in patients and HCs and between patients and HCs at baseline. Post hoc analyses tested for TSD effects for other diffusion metrics, and the relationships between change in diffusion measures with change in mood and rumination symptoms.

Results

Significant improvements in mood and rumination occurred in patients with depression (both p < 0.001), but not in HCs following TSD. Patients showed significant (p < 0.05, corrected) decreases in FA values in multiple WM tracts, including the body of the corpus callosum and anterior corona radiata post-TSD. Significant voxel-level changes in FA were not observed in HCs who received TSD (p > 0.05). However, differential effects of TSD between HCs and patients were found in the superior corona radiata, frontal WM and the posterior thalamic radiation (p < 0.05, corrected). A significant (p < 0.05) association between change in FA and axial diffusivity within the right superior corona radiata and improvement in rumination was found post-TSD in patients.

Conclusion

Total sleep deprivation leads to rapid microstructural changes in WM pathways in patients with depression that are distinct from WM changes associated with TSD observed in HCs. WM tracts including the superior corona radiata and posterior thalamic radiation could be potential biomarkers of the rapid therapeutic effects of TSD. Changes in superior corona radiata FA, in particular, may relate to improvements in maladaptive rumination.

Diagnostic accuracy of apparent diffusion coefficient to differentiate intrapancreatic accessory spleen from pancreatic neuroendocrine tumors

BACKGROUND

Intrapancreatic accessory spleen (IPAS) shares similar imaging findings with hypervascular pancreatic neuroendocrine tumors (PNETs), which may lead to unnecessary surgery.

AIM

To investigate and compare the diagnostic performance of absolute apparent diffusion coefficient (ADC) and normalized ADC (lesion-to-spleen ADC ratios) in the differential diagnosis of IPAS from PNETs.

METHODS

A retrospective study consisting of 29 patients (16 PNET patients vs 13 IPAS patients) who underwent preoperative contrast-enhanced magnetic resonance imaging together with diffusion-weighted imaging/ADC maps between January 2017 and July 2020 was performed. Two independent reviewers measured ADC on all lesions and spleens, and normalized ADC was calculated for further analysis. The receiver operating characteristics analysis was carried out for evaluating the diagnostic performance of both absolute ADC and normalized ADC values in the differential diagnosis between IPAS and PNETs by clarifying sensitivity, specificity, and accuracy. Inter-reader reliability for the two methods was evaluated.

RESULTS

IPAS had a significantly lower absolute ADC (0.931 ± 0.773 × 10^-3 mm²/s vs 1.254 ± 0.219 × 10^-3 mm²/s) and normalized ADC value (1.154 ± 0.167 vs 1.591 ± 0.364) compared to PNET. A cutoff value of 1.046 × 10^-3 mm²/s for absolute ADC was associated with 81.25% sensitivity, 100% specificity, and 89.66% accuracy with an area under the curve of 0.94 (95% confidence interval: 0.8536-1.000) for the differential diagnosis of IPAS from PNET. Similarly, a cutoff value of 1.342 for normalized ADC was associated with 81.25% sensitivity, 92.31% specificity, and 86.21% accuracy with an area under the curve of 0.91 (95% confidence interval: 0.8080-1.000) for the differential diagnosis of IPAS from PNET. Both methods showed excellent inter-reader reliability with intraclass correlation coefficients for absolute ADC and ADC ratio being 0.968 and 0.976, respectively.

CONCLUSION

Both absolute ADC and normalized ADC values can facilitate the differentiation between IPAS and PNET.

Individualised profiling of white matter organisation in moderate-to-severe traumatic brain injury patients

BACKGROUND AND PURPOSE

Approximately 65% of moderate-to-severe traumatic brain injury (m-sTBI) patients present with poor long-term behavioural outcomes, which can significantly impair activities of daily living. Numerous diffusion-weighted MRI studies have linked these poor outcomes to decreased white matter integrity of several commissural tracts, association fibres and projection fibres in the brain. However, most studies have focused on group-based analyses, which are unable to deal with the substantial between-patient heterogeneity in m-sTBI. As a result, there is increasing interest and need in conducting individualised neuroimaging analyses.

MATERIALS AND METHODS

Here, we generated a detailed subject-specific characterisation of microstructural organisation of white matter tracts in 5 chronic patients with m-sTBI (29 - 49y, 2 females), presented as a proof-of-concept. We developed an imaging analysis framework using fixel-based analysis and TractLearn to determine whether the values of fibre density of white matter tracts at the individual patient level deviate from the healthy control group (n = 12, 8F, M_age = 35.7y, age range 25 - 64y).

RESULTS

Our individualised analysis revealed unique white matter profiles, confirming the heterogenous nature of m-sTBI and the need of individualised profiles to properly characterise the extent of injury. Future studies incorporating clinical data, as well as utilising larger reference samples and examining the test-retest reliability of the fixel-wise metrics are warranted.

CONCLUSIONS

Individualised profiles may assist clinicians in tracking recovery and planning personalised training programs for chronic m-sTBI patients, which is necessary to achieve optimal behavioural outcomes and improved quality of life.

Meta-analysis of the normal diffusion tensor imaging values of the peripheral nerves in the upper limb

Peripheral neuropathy affects 1 in 10 adults over the age of 40 years. Given the absence of a reliable diagnostic test for peripheral neuropathy, there has been a surge of research into diffusion tensor imaging (DTI) because it characterises nerve microstructure and provides reproducible proxy measures of myelination, axon diameter, fibre density and organisation. Before researchers and clinicians can reliably use diffusion tensor imaging to assess the 'health' of the major nerves of the upper limb, we must understand the "normal" range of values and how they vary with experimental conditions. We searched PubMed, Embase, medRxiv and bioRxiv for studies which reported the findings of DTI of the upper limb in healthy adults. Four review authors independently triple extracted data. Using the meta suite of Stata 17, we estimated the normal fractional anisotropy (FA) and diffusivity (mean, MD; radial, RD; axial AD) values of the median, radial and ulnar nerve in the arm, elbow and forearm. Using meta-regression, we explored how DTI metrics varied with age and experimental conditions. We included 20 studies reporting data from 391 limbs, belonging to 346 adults (189 males and 154 females, ~ 1.2 M:1F) of mean age 34 years (median 31, range 20-80). In the arm, there was no difference in the FA (pooled mean 0.59 mm²/s [95% CI 0.57, 0.62]; I² 98%) or MD (pooled mean 1.13 × 10^-3 mm²/s [95% CI 1.08, 1.18]; I² 99%) of the median, radial and ulnar nerves. Around the elbow, the ulnar nerve had a 12% lower FA than the median and radial nerves (95% CI - 0.25, 0.00) and significantly higher MD, RD and AD. In the forearm, the FA (pooled mean 0.55 [95% CI 0.59, 0.64]; I² 96%) and MD (pooled mean 1.03 × 10^-3 mm²/s [95% CI 0.94, 1.12]; I² 99%) of the three nerves were similar. Multivariable meta regression showed that the b-value, TE, TR, spatial resolution and age of the subject were clinically important moderators of DTI parameters in peripheral nerves. We show that subject age, as well as the b-value, TE, TR and spatial resolution are important moderators of DTI metrics from healthy nerves in the adult upper limb. The normal ranges shown here may inform future clinical and research studies.

Vacuoles related to tissue neuron-astrocyte ratio and infiltration of macrophages/monocytes contribute to hyperintense brain signals on diffusion-weighted magnetic resonance imaging in sporadic Creutzfeldt-Jakob disease

BACKGROUND

Radiological features in patients with sporadic Creutzfeldt-Jakob disease (sCJD) are hyperintensity of the cerebral cortex and the basal ganglia displayed by diffusion-weighted magnetic resonance imaging (DW-MRI). We performed a quantitative study on neuropathological and radiological findings.

METHODS

Patient 1 received a definite diagnosis of MM1-type sCJD, while patient 2 received a definite diagnosis of MM1 + 2-type sCJD. Two DW-MRI scans were performed on each patient. DW-MRI was either taken the day before or on the day of the patient's death, and several hyperintense or isointense areas were marked as a region of interest (ROI). Mean signal intensity of the ROI was measured. Pathological quantitative assessments of the vacuoles, astrocytosis, infiltration of monocytes/macrophages, and proliferation of microglia was performed. Vacuole load (% area vacuole), glial fibrillary acidic protein (GFAP), CD68, and Iba-1 load were calculated. We defined spongiform change index (SCI) to indicate vacuoles related to a tissue neuron-astrocyte ratio. We assessed the correlation of intensity of the last DW-MRI and the pathological findings as well as association of changes of the signal intensity on the sequential images and the pathological findings.

RESULT

We observed a strong positive correlation between SCI and DW-MRI intensity. In the analysis using serial DW-MRI and pathological findings, we found that CD68 load was significantly larger in areas where signal intensity decreased, as compared to those areas where hyperintensity remained unchanged.

CONCLUSION

DW-MRI intensity in sCJD is associated with the ratio of neuron to astrocyte in the vacuoles and the infiltration of macrophages and/or monocytes.

Peak Width of Skeletonized Mean Diffusivity: A Neuroimaging Marker for White Matter Injury

A leading cause of white matter (WM) injury in older individuals is cerebral small vessel disease (SVD). Cerebral SVD is the most prevalent vascular contributor to cognitive impairment and dementia. Therapeutic progress for cerebral SVD and other WM disorders depends on the development and validation of neuroimaging markers suitable as outcome measures in future interventional trials. Diffusion-tensor imaging (DTI) is one of the best-suited MRI techniques for assessing the extent of WM damage in the brain. But the optimal method to analyze individual DTI data remains hindered by labor-intensive and time-consuming processes. Peak width of skeletonized mean diffusivity (PSMD), a recently developed fast, fully automated DTI marker, was designed to quantify the WM damage secondary to cerebral SVD and reflect related cognitive impairment. Despite its promising results, knowledge about PSMD is still limited in the radiologic community. This focused review provides an overview of the technical details of PSMD while synthesizing the available data on its clinical and neuroimaging associations. From a critical expert viewpoint, the authors discuss the limitations of PSMD and its current validation status as a neuroimaging marker for vascular cognitive impairment. Finally, they point out the gaps to be addressed to further advance the field.

Diffusion Tensor and Kurtosis Imaging Findings the First Year following Mild Traumatic Brain Injury

Despite enormous research interest in diffusion tensor imaging and diffusion kurtosis imaging (DTI; DKI) following mild traumatic brain injury (MTBI), it remains unknown how diffusion in white matter evolves post-injury and relates to acute MTBI characteristics. This prospective cohort study aimed to characterize diffusion changes in white matter the first year after MTBI. Patients with MTBI (n = 193) and matched controls (n = 83) underwent 3T magnetic resonance imaging (MRI) within 72 h and 3- and 12-months post-injury. Diffusion data were analyzed in three steps: 1) voxel-wise comparisons between the MTBI and control group were performed with tract-based spatial statistics at each time-point; 2) clusters of significant voxels identified in step 1 above were evaluated longitudinally with mixed-effect models; 3) the MTBI group was divided into: (A) complicated (with macrostructural findings on MRI) and uncomplicated MTBI; (B) long (1-24 h) and short (< 1 h) post-traumatic amnesia (PTA); and (C) other and no other concurrent injuries to investigate if findings in step 1 were driven mainly by aberrant diffusion in patients with a more severe injury. At 72 h, voxel-wise comparisons revealed significantly lower fractional anisotropy (FA) in one tract and significantly lower mean kurtosis (Kmean) in 11 tracts in the MTBI compared with control group. At 3 months, the MTBI group had significantly higher mean diffusivity in eight tracts compared with controls. At 12 months, FA was significantly lower in four tracts and Kmean in 10 tracts in patients with MTBI compared with controls. There was considerable overlap in affected tracts across time, including the corpus callosum, corona radiata, internal and external capsule, and cerebellar peduncles. Longitudinal analyses revealed that the diffusion metrics remained relatively stable throughout the first year after MTBI. The significant group*time interactions identified were driven by changes in the control rather than the MTBI group. Further, differences identified in step 1 did not result from greater diffusion abnormalities in patients with complicated MTBI, long PTA, or other concurrent injuries, as standardized mean differences in diffusion metrics between the groups were small (0.07 ± 0.11) and non-significant. However, follow-up voxel-wise analyses revealed that other concurrent injuries had effects on diffusion metrics, but predominantly in other metrics and at other time-points than the effects observed in the MTBI versus control group analysis. In conclusion, patients with MTBI differed from controls in white matter integrity already 72 h after injury. Diffusion metrics remained relatively stable throughout the first year after MTBI and were not driven by deviating diffusion in patients with a more severe MTBI.

Inter-vendor and inter-observer reliability of diffusion tensor imaging in the musculoskeletal system: a multiscanner MR study

BACKGROUND

To evaluate the inter-observer and inter-vendor reliability of diffusion tensor imaging parameters in the musculoskeletal system.

METHODS

This prospective study included six healthy volunteers three men (mean age: 42; range: 31-52 years) and three women (mean age: 36; range: 30-44 years). Each subject was scanned using different 3 Tesla magnetic resonance scanners from three different vendors at three different sites bilaterally. First, the intra-class correlation coefficient was used to determine between-observers agreement for overall measurements and clinical sites. Next, between-group comparisons were made through the nonparametric Friedman's test. Finally, the Bland-Altman method was used to determine agreement among the three scanner measurements, comparing them two by two.

RESULTS

A total of 792 measurement were calculated. ICC reported high levels of agreement between the two observers. ICC related to MD, FA, and RD measurements ranged from 0.88 (95% CI 0.85-0.90) to 0.95 (95% CI 0.94-0.96), from 0.85 (95% CI 0.81-0.88) to 0.95 (95% CI 0.93-0.96), and from 0.89 (0.85-0.90) to 0.92 (0.90-0.94). No statistically significant inter-vendor differences were observed. The Bland-Altmann method confirmed a high correlation between parameter values.

CONCLUSION

An excellent inter-observer and inter-vendor reliability was found in our study.

Association Between Quantitative Diffusion-Weighted Magnetic Resonance Neuroimaging and Outcome After Pediatric Cardiac Arrest

BACKGROUND AND OBJECTIVES

Diffusion MRI can quantify the extent of hypoxic-ischemic brain injury after cardiac arrest. Our objective was to determine the association between the adult-derived threshold of apparent diffusion coefficient (ADC) <650 × 10-6 mm2/s in >10% of brain tissue and an unfavorable outcome after pediatric cardiac arrest. Since ADC decreases exponentially as a function of increasing age, we determined the association between (1) having >10% of brain tissue below a novel age-dependent ADC threshold, and (2) age-normalized whole-brain mean ADC and unfavorable outcome.

METHODS

This was a retrospective study of patients aged ≤18 years who had cardiac arrest and a clinically obtained brain MRI within 7 days. The primary outcome was unfavorable neurologic status at hospital discharge based on the Pediatric Cerebral Performance Category score. ADC images were extracted from 3-direction diffusion imaging. We determined whether each patient had >10% of voxels with an ADC below prespecified thresholds. We computed the whole-brain mean ADC for each patient.

RESULTS

One hundred thirty-four patients were analyzed. Patients with ADC <650 × 10-6 mm2/s in >10% of voxels had 15 times higher odds (95% CI 5-65) of an unfavorable outcome compared with patients with ADC <650 × 10-6 mm2/s (area under the receiver operating characteristic curve [AUROC] 0.72 [95% CI 0.63-0.80]). These ADC criteria had a sensitivity and specificity of 0.49 and 0.94, respectively, and positive and negative predictive values of 0.93 and 0.52, respectively, for an unfavorable outcome. The age-dependent ADC threshold that yielded optimal sensitivity and specificity for unfavorable outcomes was <300 × 10-6 mm2/s below each patient's predicted whole-brain mean ADC. The sensitivity, specificity, and positive and negative predictive values for this ADC threshold were 0.53, 0.96, 0.96, and 0.54, respectively (odds ratio [OR] 26.4 [95% CI 7.5-168.3]; AUROC 0.74 [95% CI 0.66-0.83]). Lower age-normalized whole-brain mean ADC was also associated with an unfavorable outcome (OR 0.42 [0.24-0.64], AUROC 0.76 [95% CI 0.66-0.82]).

DISCUSSION

Quantitative diffusion thresholds on MRI within 7 days after cardiac arrest were associated with an unfavorable outcome in children. The age-independent ADC threshold was highly specific for predicting an unfavorable outcome. However, the specificity and sensitivity increased when using age-dependent ADC thresholds. Age-dependent ADC thresholds may improve prognostic accuracy and require further investigation in larger cohorts.

CLASSIFICATION OF EVIDENCE

This study provides Class III evidence that quantitative diffusion-weighted imaging within 7 days postarrest can predict an unfavorable clinical outcome in children.

Validation of cross-sectional and longitudinal ComBat harmonization methods for magnetic resonance imaging data on a travelling subject cohort

Background

The growth in multi-center neuroimaging studies generated a need for methods that mitigate the differences in hardware and acquisition protocols across sites i.e., scanner effects. ComBat harmonization methods have shown promise but have not yet been tested on all the data types commonly studied with magnetic resonance imaging (MRI). This study aimed to validate neuroCombat, longCombat and gamCombat on both structural and diffusion metrics in both cross-sectional and longitudinal data.

Methods

We used a travelling subject design whereby 73 healthy volunteers contributed 161 scans across two sites and four machines using one T1 and five diffusion MRI protocols. Scanner was defined as a composite of site, machine and protocol. A common pipeline extracted two structural metrics (volumes and cortical thickness) and two diffusion tensor imaging metrics (mean diffusivity and fractional anisotropy) for seven regions of interest including gray and (except for cortical thickness) white matter regions.

Results

Structural data exhibited no significant scanner effect and therefore did not benefit from harmonization in our particular cohort. Indeed, attempting harmonization obscured the true biological effect for some regions of interest. Diffusion data contained marked scanner effects and was successfully harmonized by all methods, resulting in smaller scanner effects and better detection of true biological effects. LongCombat less effectively reduced the scanner effect for cross-sectional white matter data but had a slightly lower probability of incorrectly finding group differences in simulations, compared to neuroCombat and gamCombat. False positive rates for all methods and all metrics did not significantly exceed 5%.

Conclusions

Statistical harmonization of structural data is not always necessary and harmonization in the absence of a scanner effect may be harmful. Harmonization of diffusion MRI data is highly recommended with neuroCombat, longCombat and gamCombat performing well in cross-sectional and longitudinal settings.

Evaluating Magnetic Resonance Diffusion Properties Together with Brain Volumetry May Predict Progression to Multiple Sclerosis

RATIONALE AND OBJECTIVES

Although the gold standard in predicting future progression from clinically isolated syndrome (CIS) to clinically definite multiple sclerosis (CDMS) consists in the McDonald criteria, efforts are being made to employ various advanced MRI techniques for predicting clinical progression. This study's main aim was to evaluate the predictive power of diffusion tensor imaging (DTI) of the brain and brain volumetry to distinguish between patients having CIS with future progression to CDMS from those without progression during the following 2 years and to compare those parameters with conventional MRI evaluation.

MATERIALS AND METHODS

All participants underwent an MRI scan of the brain. DTI and volumetric data were processed and various parameters were compared between the study groups.

RESULTS

We found significant differences between the subgroups of patients differing by future progression to CDMS in most of those DTI and volumetric parameters measured. Fractional anisotropy of water diffusion proved to be the strongest predictor of clinical conversion among all parameters evaluated, demonstrating also higher specificity compared to evaluation of conventional MRI images according to McDonald criteria.

CONCLUSION

Conclusion: Our results provide evidence that the evaluation of DTI parameters together with brain volumetry in patients with early-stage CIS may be useful in predicting conversion to CDMS within the following 2 years of the disease course.

Registration and quantification network (RQnet) for IVIM-DKI analysis in MRI

PURPOSE

A deep learning method is proposed for aligning diffusion weighted images (DWIs) and estimating intravoxel incoherent motion-diffusion kurtosis imaging parameters simultaneously.

METHODS

We propose an unsupervised deep learning method that performs 2 tasks: registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. A common registration method in diffusion MRI is based on minimizing dissimilarity between various DWIs, which may result in registration errors due to different contrasts in different DWIs. We designed a novel unsupervised deep learning method for both accurate registration and quantification of various diffusion parameters. In order to generate motion-simulated training data and test data, 17 volunteers were scanned without moving their heads, and 4 volunteers moved their heads during the scan in a 3 Tesla MRI. In order to investigate the applicability of the proposed method to other organs, kidney images were also obtained. We compared the registration accuracy of the proposed method, statistical parametric mapping, and a deep learning method with a normalized cross-correlation loss. In the quantification part of the proposed method, a deep learning method that considered the diffusion gradient direction was used.

RESULTS

Simulations and experimental results showed that the proposed method accurately performed registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. The registration accuracy of the proposed method was high for all b values. Furthermore, quantification performance was analyzed through simulations and in vivo experiments, where the proposed method showed the best performance among the compared methods.

CONCLUSION

The proposed method aligns the DWIs and accurately quantifies the intravoxel incoherent motion-diffusion kurtosis imaging parameters.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

Does Loss of Integrity of the Cingulum Bundle Link Amyloid-β Accumulation and Neurodegeneration in Alzheimer’s Disease?

Background: Alzheimer’s disease is characterized by the accumulation of amyloid-β (Aβ) into plaques, aggregation of tau into neurofibrillary tangles, and neurodegenerative processes including atrophy. However, there is a poorly understood spatial discordance between initial Aβ deposition and local neurodegeneration. Objective: Here, we test the hypothesis that the cingulum bundle links Aβ deposition in the cingulate cortex to medial temporal lobe (MTL) atrophy. Methods: 21 participants with mild cognitive impairment (MCI) from the UMC Utrecht memory clinic (UMCU, discovery sample) and 37 participants with MCI from Alzheimer’s Disease Neuroimaging Initiative (ADNI, replication sample) with available Aβ-PET scan, T1-weighted and diffusion-weighted MRI were included. Aβ load of the cingulate cortex was measured by the standardized uptake value ratio (SUVR), white matter integrity of the cingulum bundle was assessed by mean diffusivity and atrophy of the MTL by normalized MTL volume. Relationships were tested with linear mixed models, to accommodate multiple measures for each participant. Results: We found at most a weak association between cingulate Aβ and MTL volume (added R2 <0.06), primarily for the posterior hippocampus. In neither sample, white matter integrity of the cingulum bundle was associated with cingulate Aβ or MTL volume (added R2 <0.01). Various sensitivity analyses (Aβ-positive individuals only, posterior cingulate SUVR, MTL sub region volume) provided similar results. Conclusion: These findings, consistent in two independent cohorts, do not support our hypothesis that loss of white matter integrity of the cingulum is a connecting factor between cingulate gyrus Aβ deposition and MTL atrophy.

Early Onset Diffusion Abnormalities in Refractory Headache Disorders

Objective

This study sought to determine if individuals with medically refractory migraine headache have volume or diffusion abnormalities on neuroimaging compared to neurotypical individuals.

Background

Neuroimaging biomarkers in headache medicine continue to be limited. Early prediction of medically refractory headache and migraine disorders could result in earlier administration of high efficacy therapeutics.

Methods

A single-center, retrospective, case control study was performed. All patients were evaluated clinically between 2014 and 2018. Individuals with medically refractory migraine headache (defined by ICDH-3 criteria) without any other chronic medical diseases were enrolled. Patients had to have failed more than two therapeutics and aura was not exclusionary. The initial MRI study for each patient was reviewed. Multiple brain regions were analyzed for volume and apparent diffusion coefficient values. These were compared to 81 neurotypical control patients.

Results

A total of 79 patients with medically refractory migraine headache were included and compared to 74 neurotypical controls without headache disorders. Time between clinical diagnosis and neuroimaging was a median of 24 months (IQR: 12.0–37.0). Comparison of individuals with medically refractory migraine headache to controls revealed statistically significant differences in median apparent diffusion coefficient (ADC) in multiple brain subregions (p < 0.001). Post-hoc pair-wise analysis comparing individuals with medically refractory migraine headache to control patients revealed significantly decreased median ADC values for the thalamus, caudate, putamen, pallidum, amygdala, brainstem, and cerebral white matter. No volumetric differences were observed between groups.

Conclusions

In individuals with medically refractory MH, ADC changes are measurable in multiple brain structures at an early age, prior to the failure of multiple pharmacologic interventions and the diagnosis of medically refractory MH. This data supports the hypothesis that structural connectivity issues may predispose some patients toward more medically refractory pain disorders such as MH.

WA, Boyle J, Coll‐Font J, Dall'Armellina E, Ennis DB, Ferreira PF, Kalra P, Kolipaka A, Kozerke S, Lohr D, Mongeon F, Moulin K, Nguyen C, Nielles‐Vallespin S, Raterman B, Schreiber LM, Scott AD, Sosnovik DE, Stoeck CT, Tous C, Tunnicliffe EM, Weng AM, Croisille P, Viallon M, Schneider JE. Validation of cardiac diffusion tensor imaging sequences: A multicentre test-retest phantom study

Cardiac diffusion tensor imaging (DTI) is an emerging technique for the in vivo characterisation of myocardial microstructure, and there is a growing need for its validation and standardisation. We sought to establish the accuracy, precision, repeatability and reproducibility of state-of-the-art pulse sequences for cardiac DTI among 10 centres internationally. Phantoms comprising 0%-20% polyvinylpyrrolidone (PVP) were scanned with DTI using a product pulsed gradient spin echo (PGSE; N = 10 sites) sequence, and a custom motion-compensated spin echo (SE; N = 5) or stimulated echo acquisition mode (STEAM; N = 5) sequence suitable for cardiac DTI in vivo. A second identical scan was performed 1-9 days later, and the data were analysed centrally. The average mean diffusivities (MDs) in 0% PVP were (1.124, 1.130, 1.113) x 10^-3  mm² /s for PGSE, SE and STEAM, respectively, and accurate to within 1.5% of reference data from the literature. The coefficients of variation in MDs across sites were 2.6%, 3.1% and 2.1% for PGSE, SE and STEAM, respectively, and were similar to previous studies using only PGSE. Reproducibility in MD was excellent, with mean differences in PGSE, SE and STEAM of (0.3 ± 2.3, 0.24 ± 0.95, 0.52 ± 0.58) x 10^-5  mm² /s (mean ± 1.96 SD). We show that custom sequences for cardiac DTI provide accurate, precise, repeatable and reproducible measurements. Further work in anisotropic and/or deforming phantoms is warranted.

Challenges in ensuring the generalizability of image quantitation methods for MRI

Image quantitation methods including quantitative MRI, multiparametric MRI, and radiomics offer great promise for clinical use. However, many of these methods have limited clinical adoption, in part due to issues of generalizability, that is, the ability to translate methods and models across institutions. Researchers can assess generalizability through measurement of repeatability and reproducibility, thus quantifying different aspects of measurement variance. In this article, we review the challenges to ensuring repeatability and reproducibility of image quantitation methods as well as present strategies to minimize their variance to enable wider clinical implementation. We present possible solutions for achieving clinically acceptable performance of image quantitation methods and briefly discuss the impact of minimizing variance and achieving generalizability towards clinical implementation and adoption.

Regional Distribution of Brain Injury After Cardiac Arrest: Clinical and Electrographic Correlates

BACKGROUND AND OBJECTIVES

Disorders of consciousness, EEG background suppression, and epileptic seizures are associated with poor outcome after cardiac arrest. Our objective was to identify the distribution of diffusion MRI-measured anoxic brain injury after cardiac arrest and to define the regional correlates of disorders of consciousness, EEG background suppression, and seizures.

METHODS

We analyzed patients from a single-center database of unresponsive patients who underwent diffusion MRI after cardiac arrest (n = 204). We classified each patient according to recovery of consciousness (command following) before discharge, the most continuous EEG background (burst suppression vs continuous), and the presence or absence of seizures. Anoxic brain injury was measured with the apparent diffusion coefficient (ADC) signal. We identified ADC abnormalities relative to controls without cardiac arrest (n = 48) and used voxel lesion symptom mapping to identify regional associations with disorders of consciousness, EEG background suppression, and seizures. We then used a bootstrapped lasso regression procedure to identify robust, multivariate regional associations with each outcome variable. Last, using area under receiver operating characteristic curves, we then compared the classification ability of the strongest regional associations to that of brain-wide summary measures.

RESULTS

Compared to controls, patients with cardiac arrest demonstrated ADC signal reduction that was most significant in the occipital lobes. Disorders of consciousness were associated with reduced ADC most prominently in the occipital lobes but also in deep structures. Regional injury more accurately classified patients with disorders of consciousness than whole-brain injury. Background suppression mapped to a similar set of brain regions, but regional injury could no better classify patients than whole-brain measures. Seizures were less common in patients with more severe anoxic injury, particularly in those with injury to the lateral temporal white matter.

DISCUSSION

Anoxic brain injury was most prevalent in posterior cerebral regions, and this regional pattern of injury was a better predictor of disorders of consciousness than whole-brain injury measures. EEG background suppression lacked a specific regional association, but patients with injury to the temporal lobe were less likely to have seizures. Regional patterns of anoxic brain injury are relevant to the clinical and electrographic sequelae of cardiac arrest and may hold importance for prognosis.

CLASSIFICATION OF EVIDENCE

This study provides Class IV evidence that disorders of consciousness after cardiac arrest are associated with widely lower ADC values on diffusion MRI and are most strongly associated with reductions in occipital ADC.

Predictive MRI Biomarkers in MS—A Critical Review

Background and Objectives: In this critical review, we explore the potential use of MRI measurements as prognostic biomarkers in multiple sclerosis (MS) patients, for both conventional measurements and more novel techniques such as magnetization transfer, diffusion tensor, and proton spectroscopy MRI. Materials and Methods: All authors individually and comprehensively reviewed each of the aspects listed below in PubMed, Medline, and Google Scholar. Results: There are numerous MRI metrics that have been proven by clinical studies to hold important prognostic value for MS patients, most of which can be readily obtained from standard 1.5T MRI scans. Conclusions: While some of these parameters have passed the test of time and seem to be associated with a reliable predictive power, some are still better interpreted with caution. We hope this will serve as a reminder of how vast a resource we have on our hands in this versatile tool—it is up to us to make use of it.

Initial Experience on Hyperpolarized [1-13C]Pyruvate MRI Multicenter Reproducibility—Are Multicenter Trials Feasible?

Background: Magnetic resonance imaging (MRI) with hyperpolarized [1-¹³C]pyruvate allows real-time and pathway specific clinical detection of otherwise unimageable in vivo metabolism. However, the comparability between sites and protocols is unknown. Here, we provide initial experiences on the agreement of hyperpolarized MRI between sites and protocols by repeated imaging of same healthy volunteers in Europe and the US. Methods: Three healthy volunteers traveled for repeated multicenter brain MRI exams with hyperpolarized [1-¹³C]pyruvate within one year. First, multisite agreement was assessed with the same echo-planar imaging protocol at both sites. Then, this was compared to a variable resolution echo-planar imaging protocol. In total, 12 examinations were performed. Common metrics of ¹³C-pyruvate to ¹³C-lactate conversion were calculated, including the k_PL, a model-based kinetic rate constant, and its model-free equivalents. Repeatability was evaluated with intraclass correlation coefficients (ICC) for absolute agreement computed using two-way random effects models. Results: The mean k_PL across all examinations in the multisite comparison was 0.024 ± 0.0016 s⁻¹. The ICC of the k_PL was 0.83 (p = 0.14) between sites and 0.7 (p = 0.09) between examinations of the same volunteer at any of the two sites. For the model-free metrics, the lactate Z-score had similar site-to-site ICC, while it was considerably lower for the lactate-to-pyruvate ratio. Conclusions: Estimation of metabolic conversion from hyperpolarized [1-¹³C]pyruvate to lactate using model-based metrics such as k_PL suggests close agreement between sites and examinations in volunteers. Our initial results support harmonization of protocols, support multicenter studies, and inform their design.

Reproducibility of diffusion tensor image analysis along the perivascular space (DTI-ALPS) for evaluating interstitial fluid diffusivity and glymphatic function: CHanges in Alps index on Multiple conditiON acquIsition eXperiment (CHAMONIX) study

PURPOSE

The diffusion tensor image analysis along the perivascular space (DTI-ALPS) method was developed to evaluate the brain's glymphatic function or interstitial fluid dynamics. This study aimed to evaluate the reproducibility of the DTI-ALPS method and the effect of modifications in the imaging method and data evaluation.

MATERIALS AND METHODS

Seven healthy volunteers were enrolled in this study. Image acquisition was performed for this test-retest study using a fixed imaging sequence and modified imaging methods which included the placement of region of interest (ROI), imaging plane, head position, averaging, number of motion-proving gradients, echo time (TE), and a different scanner. The ALPS-index values were evaluated for the change of conditions listed above.

RESULTS

This test-retest study by a fixed imaging sequence showed very high reproducibility (intraclass coefficient = 0.828) for the ALPS-index value. The bilateral ROI placement showed higher reproducibility. The number of averaging and the difference of the scanner did not influence the ALPS-index values. However, modification of the imaging plane and head position impaired reproducibility, and the number of motion-proving gradients affected the ALPS-index value. The ALPS-index values from 12-axis DTI and 3-axis diffusion-weighted image (DWI) showed good correlation (r = 0.86). Also, a shorter TE resulted in a larger value of the ALPS-index.

CONCLUSION

ALPS index was robust under the fixed imaging method even when different scanners were used. ALPS index was influenced by the imaging plane, the number of motion-proving gradient axes, and TE in the imaging sequence. These factors should be uniformed in the planning ALPS method studies. The possibility to develop a 3-axis DWI-ALPS method using three axes of the motion-proving gradient was also suggested.

Differences in the choroid plexus volume and microstructure are associated with body adiposity

The choroid plexus (CP) is a cerebral structure located in the ventricles that functions in producing most of the brain's cerebrospinal fluid (CSF) and transporting proteins and immune cells. Alterations in CP structure and function has been implicated in several pathologies including aging, multiple sclerosis, Alzheimer's disease, and stroke. However, identification of changes in the CP remains poorly characterized in obesity, one of the main risk factors of neurodegeneration, including in the absence of frank central nervous system alterations. Our goal here was to characterize the association between obesity, measured by the body mass index (BMI) or waist circumference (WC) metrics, and CP microstructure and volume, assessed using advanced magnetic resonance imaging (MRI) methodology. This cross-sectional study was performed in the clinical unit of the National Institute on Aging and included a participant population of 123 cognitively unimpaired individuals spanning the age range of 22 - 94 years. Automated segmentation methods from FreeSurfer were used to identify the CP structure. Our analysis included volumetric measurements, quantitative relaxometry measures (T ₁ and T ₂), and the diffusion tensor imaging (DTI) measure of mean diffusivity (MD). Strong positive associations were observed between WC and all MRI metrics, as well as CP volume. When comparing groups based on the established cutoff point by the National Institutes of Health for WC, a modest difference in MD and a significant difference in T ₁ values were observed between obese and lean individuals. We also found differences in T1 and MD between obese and overweight individuals as defined using the BMI cutoff. We conjecture that these observations in CP volume and microstructure are due to obesity-induced inflammation, diet, or, very likely, dysregulations in leptin binding and transport. These findings demonstrate that obesity is strongly associated with a decline in CP microstructural integrity. We expect that this work will lay the foundation for further investigations on obesity-induced alterations in CP structure and function.

MASiVar: Multisite, multiscanner, and multisubject acquisitions for studying variability in diffusion weighted MRI

PURPOSE

Diffusion-weighted imaging allows investigators to identify structural, microstructural, and connectivity-based differences between subjects, but variability due to session and scanner biases is a challenge.

METHODS

To investigate DWI variability, we present MASiVar, a multisite data set consisting of 319 diffusion scans acquired at 3 T from b = 1000 to 3000 s/mm2 across 14 healthy adults, 83 healthy children (5 to 8 years), three sites, and four scanners as a publicly available, preprocessed, and de-identified data set. With the adult data, we demonstrate the capacity of MASiVar to simultaneously quantify the intrasession, intersession, interscanner, and intersubject variability of four common DWI processing approaches: (1) a tensor signal representation, (2) a multi-compartment neurite orientation dispersion and density model, (3) white-matter bundle segmentation, and (4) structural connectomics. Respectively, we evaluate region-wise fractional anisotropy, mean diffusivity, and principal eigenvector; region-wise CSF volume fraction, intracellular volume fraction, and orientation dispersion index; bundle-wise shape, volume, fractional anisotropy, and length; and whole connectome correlation and maximized modularity, global efficiency, and characteristic path length.

RESULTS

We plot the variability in these measures at each level and find that it consistently increases with intrasession to intersession to interscanner to intersubject effects across all processing approaches and that sometimes interscanner variability can approach intersubject variability.

CONCLUSIONS

This study demonstrates the potential of MASiVar to more globally investigate DWI variability across multiple levels and processing approaches simultaneously and suggests harmonization between scanners for multisite analyses should be considered before inference of group differences on subjects.

Functional Connectivity-Based Prediction of Autism on Site Harmonized ABIDE Dataset

OBJECTIVE

The larger sample sizes available from multi-site publicly available neuroimaging data repositories makes machine-learning based diagnostic classification of mental disorders more feasible by alleviating the curse of dimensionality. However, since multi-site data are aggregated post-hoc, i.e. they were acquired from different scanners with different acquisition parameters, non-neural inter-site variability may mask inter-group differences that are at least in part neural in origin. Hence, the advantages gained by the larger sample size in the context of machine-learning based diagnostic classification may not be realized.

METHODS

We address this issue using harmonization of multi-site neuroimaging data using the ComBat technique, which is based on an empirical Bayes formulation to remove inter-site differences in data distributions, to improve diagnostic classification accuracy. Specifically, we demonstrate this using ABIDE (Autism Brain Imaging Data Exchange) multi-site data for classifying individuals with Autism from healthy controls using resting state fMRI-based functional connectivity data.

RESULTS

Our results show that higher classification accuracies across multiple classification models can be obtained (especially for models based on artificial neural networks) from multi-site data post harmonization with the ComBat technique as compared to without harmonization, outperforming earlier results from existing studies using ABIDE. Furthermore, our network ablation analysis facilitated important insights into autism spectrum disorder pathology and the connectivity in networks shown to be important for classification covaried with verbal communication impairments in Autism.

CONCLUSION

Multi-site data harmonization using ComBat improves neuroimaging-based diagnostic classification of mental disorders.

SIGNIFICANCE

ComBat has the potential to make AI-based clinical decision-support systems more feasible in psychiatry.