Contrast-agent-based perfusion MRI code repository and testing framework: ISMRM Open Science Initiative for Perfusion Imaging (OSIPI).

PURPOSE: Software has a substantial impact on quantitative perfusion MRI values. The lack of generally accepted implementations, code sharing and transparent testing reduces reproducibility, hindering the use of perfusion MRI in clinical trials. To address these issues, the ISMRM Open Science Initiative for Perfusion Imaging (OSIPI) aimed to establish a community-led, centralized repository for sharing open-source code for processing contrast-based perfusion imaging, incorporating an open-source testing framework. METHODS: A repository was established on the OSIPI GitHub website. Python was chosen as the target software language. Calls for code contributions were made to OSIPI members, the ISMRM Perfusion Study Group, and publicly via OSIPI websites. An automated unit-testing framework was implemented to evaluate the output of code contributions, including visual representation of the results. RESULTS: The repository hosts 86 implementations of perfusion processing steps contributed by 12 individuals or teams. These cover all core aspects of DCE- and DSC-MRI processing, including multiple implementations of the same functionality. Tests were developed for 52 implementations, covering five analysis steps. For T1 mapping, signal-to-concentration conversion and population AIF functions, different implementations resulted in near-identical output values. For the five pharmacokinetic models tested (Tofts, extended Tofts-Kety, Patlak, two-compartment exchange, and two-compartment uptake), differences in output parameters were observed between contributions. CONCLUSIONS: The OSIPI DCE-DSC code repository represents a novel community-led model for code sharing and testing. The repository facilitates the re-use of existing code and the benchmarking of new code, promoting enhanced reproducibility in quantitative perfusion imaging.

A community-endorsed open-source lexicon for contrast agent-based perfusion MRI: A consensus guidelines report from the ISMRM Open Science Initiative for Perfusion Imaging (OSIPI).

This manuscript describes the ISMRM OSIPI (Open Science Initiative for Perfusion Imaging) lexicon for dynamic contrast-enhanced and dynamic susceptibility-contrast MRI. The lexicon was developed by Taskforce 4.2 of OSIPI to provide standardized definitions of commonly used quantities, models, and analysis processes with the aim of reducing reporting variability. The taskforce was established in February 2020 and consists of medical physicists, engineers, clinicians, data and computer scientists, and DICOM (Digital Imaging and Communications in Medicine) standard experts. Members of the taskforce collaborated via a slack channel and quarterly virtual meetings. Members participated by defining lexicon items and reporting formats that were reviewed by at least two other members of the taskforce. Version 1.0.0 of the lexicon was subject to open review from the wider perfusion imaging community between January and March 2022, and endorsed by the Perfusion Study Group of the ISMRM in the summer of 2022. The initial scope of the lexicon was set by the taskforce and defined such that it contained a basic set of quantities, processes, and models to enable users to report an end-to-end analysis pipeline including kinetic model fitting. We also provide guidance on how to easily incorporate lexicon items and definitions into free-text descriptions (e.g., in manuscripts and other documentation) and introduce an XML-based pipeline encoding format to encode analyses using lexicon definitions in standardized and extensible machine-readable code. The lexicon is designed to be open-source and extendable, enabling ongoing expansion of its content. We hope that widespread adoption of lexicon terminology and reporting formats described herein will increase reproducibility within the field.

The future of MRI in radiation therapy: Challenges and opportunities for the MR community.

Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field.

Eye-specific quantitative dynamic contrast-enhanced MRI analysis for patients with intraocular masses.

OBJECTIVE: Dynamic contrast enhanced (DCE)-MRI is currently not generally used for intraocular masses as lesions are small, have an inhomogeneous T1 and the eye is prone to motion. The aim of this paper is to address these eye-specific challenges, enabling accurate ocular DCE-MRI. MATERIALS & METHODS: DCE-MRI of 19 uveal melanoma (UM) patients was acquired using a fat-suppressed 3D spoiled gradient echo sequence with TWIST (time-resolved angiography with stochastic trajectories sequence). The analysis consisted of a two-step registration method to correct for both head and eye motion. A T1 map was calculated to convert signal intensities to concentrations. Subsequently, the Tofts model was fitted voxel wise to obtain Ktrans and ve. RESULTS: Registration significantly improved the concentration curve quality (p < 0.001). The T1 of melanotic lesions was significantly lower than amelanotic lesions (888 ms vs 1350 ms, p = 0.03). The average achieved B1+ in the lesions was 91%. The average Ktrans was 0.46 min-1 (range 0.13-1.0) and the average ve was 0.22 (range 0.10-0.51). CONCLUSION: Using this eye-specific analysis, DCE of intraocular masses is possible which might aid in the diagnosis, prognosis and follow-up of UM.

Histopathological Features of MRI-Invisible Regions of Prostate Cancer Lesions.

BACKGROUND: Previous studies have reported tumor volume underestimation with multiparametric (mp)MRI in prostate cancer diagnosis. PURPOSE: To investigate why some parts of lesions are not visible on mpMRI by comparing their histopathology features to those of visible regions. STUDY TYPE: Retrospective. POPULATION: Thirty-four patients with biopsy-proven prostate cancer scheduled for prostatectomy (median 68.7 years). FIELD STRENGTH/SEQUENCE: T2 -weighted, diffusion-weighted imaging, T2 mapping, and dynamic contrast-enhanced MRI on two 3T systems and one 1.5T system. ASSESSMENT: Two readers delineated suspicious lesions on mpMRI. A pathologist delineated the lesions on histopathology. A patient-customized mold enabled the registration of histopathology and MRI. On histopathology we identified mpMRI visible and invisible lesions. Subsequently, within the visible lesions we identified regions that were visible and regions that were invisible on mpMRI. For each lesion and region the following characteristics were determined: size, location, International Society of Urological Pathology (ISUP) grade, and Gleason subpatterns (density [dense/intermediate], tumor morphology [homogeneous/heterogeneous], cribriform growth [yes/no]). STATISTICAL TESTS: With generalized linear mixed-effect modeling we investigated which features explain why a lesion or a region was invisible on MRI. We compared imaging values (T2 , ADC, and Ktrans ) for these features with one-way analysis of variance (ANOVA). RESULTS: Small, anterior, and ISUP grade 1-2 lesions (n = 34) were missed more frequent than large, posterior, ISUP grade ≥ 3 lesions (n = 35). Invisible regions on mpMRI had lower tumor density, heterogeneous tumor morphology, and were located in the transition zone. Both T2 and ADC values were higher in "intermediate" compared with "dense" regions (P = 0.002 and < 0.001) and in regions with heterogeneous compared with homogeneous morphology (P < 0.001 and 0.03). Ktrans was not significantly different (P = 0.24 and 0.99). DATA CONCLUSION: Regions of prostate cancer lesions that are invisible on mpMRI have different histopathology features than visible regions. This may have implications for monitoring during active surveillance and focal treatment strategies. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2020;51:1235-1246.

Improved repeatability of dynamic contrast-enhanced MRI using the complex MRI signal to derive arterial input functions: a test-retest study in prostate cancer patients.

PURPOSE: The arterial input function (AIF) is a major source of uncertainty in tracer kinetic (TK) analysis of dynamic contrast-enhanced (DCE)-MRI data. The aim of this study was to investigate the repeatability of AIFs extracted from the complex signal and of the resulting TK parameters in prostate cancer patients. METHODS: Twenty-two patients with biopsy-proven prostate cancer underwent a 3T MRI exam twice. DCE-MRI data were acquired with a 3D spoiled gradient echo sequence. AIFs were extracted from the magnitude of the signal (AIFMAGN ), phase (AIFPHASE ), and complex signal (AIFCOMPLEX ). The Tofts model was applied to extract Ktrans , kep and ve . Repeatability of AIF curve characteristics and TK parameters was assessed with the within-subject coefficient of variation (wCV). RESULTS: The wCV for peak height and full width at half maximum for AIFCOMPLEX (7% and 8%) indicated an improved repeatability compared to AIFMAGN (12% and 12%) and AIFPHASE (12% and 7%). This translated in lower wCV values for Ktrans (11%) with AIFCOMPLEX in comparison to AIFMAGN (24%) and AIFPHASE (15%). For kep , the wCV was 16% with AIFMAGN , 13% with AIFPHASE , and 13% with AIFCOMPLEX . CONCLUSION: Repeatability of AIFPHASE and AIFCOMPLEX is higher than for AIFMAGN , resulting in a better repeatability of TK parameters. Thus, use of either AIFPHASE or AIFCOMPLEX improves the robustness of quantitative analysis of DCE-MRI in prostate cancer.

Quantitative 3T multiparametric MRI of benign and malignant prostatic tissue in patients with and without local recurrent prostate cancer after external-beam radiation therapy.

BACKGROUND: Post-radiotherapy locally recurrent prostate cancer (PCa) patients are candidates for focal salvage treatment. Multiparametric MRI (mp-MRI) is attractive for tumor localization. However, radiotherapy-induced tissue changes complicate image interpretation. To develop focal salvage strategies, accurate tumor localization and distinction from benign tissue is necessary. PURPOSE: To quantitatively characterize radio-recurrent tumor and benign radiation-induced changes using mp-MRI, and investigate which sequences optimize the distinction between tumor and benign surroundings. STUDY TYPE: Prospective case-control. SUBJECTS: Thirty-three patients with biochemical failure after external-beam radiotherapy (cases), 35 patients without post-radiotherapy recurrent disease (controls), and 13 patients with primary PCa (untreated). FIELD STRENGTH/SEQUENCES: 3T; quantitative mp-MRI: T2 -mapping, ADC, and Ktrans and kep maps. ASSESSMENT: Quantitative image-analysis of prostatic regions, within and between cases, controls, and untreated patients. STATISTICAL TESTS: Within-groups: nonparametric Friedman analysis of variance with post-hoc Wilcoxon signed-rank tests; between-groups: Mann-Whitney tests. All with Bonferroni corrections. Generalized linear mixed modeling to ascertain the contribution of each map and location to tumor likelihood. RESULTS: Benign imaging values were comparable between cases and controls (P = 0.15 for ADC in the central gland up to 0.91 for kep in the peripheral zone), both with similarly high peri-urethral Ktrans and kep values (min-1 ) (median [range]: Ktrans = 0.22 [0.14-0.43] and 0.22 [0.14-0.36], P = 0.60, kep = 0.43 [0.24-0.57] and 0.48 [0.32-0.67], P = 0.05). After radiotherapy, benign central gland values were significantly decreased for all maps (P ≤ 0.001) as well as T2 , Ktrans , and kep of benign peripheral zone (all with P ≤ 0.002). All imaging maps distinguished recurrent tumor from benign peripheral zone, but only ADC, Ktrans , and kep were able to distinguish it from benign central gland. Recurrent tumor and peri-urethral Ktrans values were not significantly different (P = 0.81), but kep values were (P < 0.001). Combining all quantitative maps and voxel location resulted in an optimal distinction between tumor and benign voxels. DATA CONCLUSION: Mp-MRI can distinguish recurrent tumor from benign tissue. LEVEL OF EVIDENCE: 2 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:269-278.

Quantitative 3-T multi-parametric MRI and step-section pathology of recurrent prostate cancer patients after radiation therapy.

OBJECTIVES: Diagnosis of radio-recurrent prostate cancer using multi-parametric MRI (mp-MRI) can be challenging due to the presence of radiation effects. We aim to characterize imaging of prostate tissue after radiation therapy (RT), using histopathology as ground truth, and to investigate the visibility of tumor lesions on mp-MRI. METHODS: Tumor delineated histopathology slides from salvage radical prostatectomy patients, primarily treated with RT, were registered to MRI. Median T2-weighted, ADC, Ktrans, and kep values in tumor and other regions were calculated. Two radiologists independently performed mp-MRI-based tumor delineations which were compared with the true pathological extent. General linear mixed-effect modeling was used to establish the contribution of each imaging modality and combinations thereof in distinguishing tumor and benign voxels. RESULTS: Nineteen of the 21 included patients had tumor in the available histopathology slides. Recurrence was predominantly multifocal with large tumor foci seen after external beam radiotherapy, whereas these were small and sparse after low-dose-rate brachytherapy. MRI-based delineations missed small foci and slightly underestimated tumor extent. The combination of T2-weighted, ADC, Ktrans, and kep had the best performance in distinguishing tumor and benign voxels. CONCLUSIONS: Using high-resolution histopathology delineations, the real tumor extent and size were found to be underestimated on MRI. mp-MRI obtained the best performance in identifying tumor voxels. Appropriate margins around the visible tumor-suspected region should be included when designing focal salvage strategies. Recurrent tumor delineation guidelines are warranted. KEY POINTS: • Compared to the use of individual sequences, multi-parametric MRI obtained the best performance in distinguishing recurrent tumor from benign voxels. • Delineations based on mp-MRI miss smaller foci and slightly underestimate tumor volume of local recurrent prostate cancer. • Focal salvage strategies should include appropriate margins around the visible tumor.

Performance of a fast and high-resolution multi-echo spin-echo sequence for prostate T2 mapping across multiple systems.

PURPOSE: To evaluate the performance of a multi-echo spin-echo sequence with k-t undersampling scheme (k-t T2 ) in prostate cancer. METHODS: Phantom experiments were performed at five systems to estimate the bias, short-term repeatability, and reproducibility across all systems expressed with the within-subject coefficient of variation (wCV). Monthly measurements were performed on two systems for long-term repeatability estimation. To evaluate clinical repeatability, two T2 maps (voxel size 0.8 × 0.8 × 3 mm3 ; 5 min) were acquired at separate visits on one system for 13 prostate cancer patients. Repeatability was assessed per patient in relation to spatial resolution. T2 values were compared for tumor, peripheral zone, and transition zone. RESULTS: Phantom measurements showed a small bias (median = -0.9 ms) and good short-term repeatability (median wCV = 0.5%). Long-term repeatability was 0.9 and 1.1% and reproducibility between systems was 1.7%. The median bias observed in patients was -1.1 ms. At voxel level, the median wCV was 15%, dropping to 4% for structures of 0.5 cm3 . The median tumor T2 values (79 ms) were significantly lower (P < 0.001) than in the peripheral zone (149 ms), but overlapped with the transition zone (91 ms). CONCLUSIONS: Reproducible T2 mapping of the prostate is feasible with good spatial resolution in a clinically reasonable scan time, allowing reliable measurement of T2 in structures as small as 0.5 cm3 . Magn Reson Med 79:1586-1594, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Impact of contrast agent injection duration on dynamic contrast-enhanced MRI quantification in prostate cancer.

The volume transfer constant Ktrans , which describes the leakage of contrast agent (CA) from vasculature into tissue, is the most commonly reported quantitative parameter for dynamic contrast-enhanced (DCE-) MRI. However, the variation in reported Ktrans values between studies from different institutes is large. One of the primary sources of uncertainty is quantification of the arterial input function (AIF). The aim of this study is to determine the influence of the CA injection duration on the AIF and tracer kinetic analysis (TKA) parameters (i.e. Ktrans , kep and ve ). Thirty-one patients with prostate cancer received two DCE-MRI examinations with an injection duration of 5 s in the first examination and a prolonged injection duration in the second examination, varying between 7.5 s and 30 s. The DCE examination was carried out on a 3.0 T MRI scanner using a transversal T1 -weighted 3D spoiled gradient echo sequence (300 s duration, dynamic scan time of 2.5 s). Data of 29 of the 31 were further analysed. AIFs were determined from the phase signal in the left and right femoral arteries. Ktrans , kep and ve were estimated with the standard Tofts model for regions of healthy peripheral zone and tumour tissue. We observed a significantly smaller peak height and increased width in the AIF for injection durations of 15 s and longer. However, we did not find significant differences in Ktrans , kep or ve for the studied injection durations. The study demonstrates that the TKA parameters Ktrans , kep and ve , measured in the prostate, do not show a significant change as a function of injection duration.

Are Epilepsy-Related fMRI Components Dependent on the Presence of Interictal Epileptic Discharges in Scalp EEG?

Spatial independent component analysis (ICA) is increasingly being used to extract resting-state networks from fMRI data. Previous studies showed that ICA also reveals independent components (ICs) related to the seizure onset zone. However, it is currently unknown how these epileptic ICs depend on the presence of interictal epileptic discharges (IEDs) in the EEG. The goal of this study was to explore the relation between ICs obtained from fMRI epochs during the occurrence of IEDs in the EEG and those without IEDs. fMRI data sets with co-registered EEG were retrospectively selected of patients from whom the location of the epileptogenic zone was confirmed by outcome of surgery (n = 8). The fMRI data were split into two epochs: one with IEDs visible in scalp EEG and one without. Spatial ICA was applied to the fMRI data of each part separately. The maps of all resulting components were compared to the resection area and the EEG-fMRI correlation pattern by computing a spatial correlation coefficient to detect the epilepsy-related component. For all patients, except one, there was a remarkable resemblance between the epilepsy-related components selected during epochs with IEDs and those without IEDs. These findings suggest that epilepsy-related ICs are not dependent on the presence of IEDs in scalp EEG. Since these epileptic ICs showed partial overlap with resting-state networks of healthy volunteers (n = 10), our study supports the need for new ways to classify epileptic ICs.

Quantitative MRI on MR-Linacs: Towards Biological Image-Guided Adaptive Radiotherapy.

Recognizing the potential of quantitative imaging biomarkers (QIBs) in radiotherapy, many studies have investigated the prognostic value of quantitative MRI (qMRI). With the introduction of MRI-guided radiotherapy systems, the practical challenges of repeated imaging have been substantially reduced. Since patients are treated inside an MRI scanner, acquisition of qMRI can be done during each fraction with limited or no prolongation of the fraction duration. In this review paper, we identify the steps that need been taken to move from MR as an imaging technique to a useful biomarker for MRI-guided radiotherapy (MRgRT).

Novel artefact removal algorithms for co-registered EEG/fMRI based on selective averaging and subtraction.

Co-registered EEG and functional MRI (EEG/fMRI) is a potential clinical tool for planning invasive EEG in patients with epilepsy. In addition, the analysis of EEG/fMRI data provides a fundamental insight into the precise physiological meaning of both fMRI and EEG data. Routine application of EEG/fMRI for localization of epileptic sources is hampered by large artefacts in the EEG, caused by switching of scanner gradients and heartbeat effects. Residuals of the ballistocardiogram (BCG) artefacts are similarly shaped as epileptic spikes, and may therefore cause false identification of spikes. In this study, new ideas and methods are presented to remove gradient artefacts and to reduce BCG artefacts of different shapes that mutually overlap in time. Gradient artefacts can be removed efficiently by subtracting an average artefact template when the EEG sampling frequency and EEG low-pass filtering are sufficient in relation to MR gradient switching (Gonçalves et al., 2007). When this is not the case, the gradient artefacts repeat themselves at time intervals that depend on the remainder between the fMRI repetition time and the closest multiple of the EEG acquisition time. These repetitions are deterministic, but difficult to predict due to the limited precision by which these timings are known. Therefore, we propose to estimate gradient artefact repetitions using a clustering algorithm, combined with selective averaging. Clustering of the gradient artefacts yields cleaner EEG for data recorded during scanning of a 3T scanner when using a sampling frequency of 2048 Hz. It even gives clean EEG when the EEG is sampled with only 256 Hz. Current BCG artefacts-reduction algorithms based on average template subtraction have the intrinsic limitation that they fail to deal properly with artefacts that overlap in time. To eliminate this constraint, the precise timings of artefact overlaps were modelled and represented in a sparse matrix. Next, the artefacts were disentangled with a least squares procedure. The relevance of this approach is illustrated by determining the BCG artefacts in a data set consisting of 29 healthy subjects recorded in a 1.5 T scanner and 15 patients with epilepsy recorded in a 3 T scanner. Analysis of the relationship between artefact amplitude, duration and heartbeat interval shows that in 22% (1.5T data) to 30% (3T data) of the cases BCG artefacts show an overlap. The BCG artefacts of the EEG/fMRI data recorded on the 1.5T scanner show a small negative correlation between HBI and BCG amplitude. In conclusion, the proposed methodology provides a substantial improvement of the quality of the EEG signal without excessive computer power or additional hardware than standard EEG-compatible equipment.

EEG-fMRI correlation patterns in the presurgical evaluation of focal epilepsy: a comparison with electrocorticographic data and surgical outcome measures.

EEG-correlated functional MRI (EEG-fMRI) visualizes brain regions associated with interictal epileptiform discharges (IEDs). This technique images the epileptiform network, including multifocal, superficial and deeply situated cortical areas. To understand the role of EEG-fMRI in presurgical evaluation, its results should be validated relative to a gold standard. For that purpose, EEG-fMRI data were acquired for a heterogeneous group of surgical candidates (n=16) who were later implanted with subdural grids and strips (ECoG). The EEG-fMRI correlation patterns were systematically compared with brain areas involved in IEDs ECoG, using a semi-automatic analysis method, as well as to the seizure onset zone, resected area, and degree of seizure freedom. In each patient at least one of the EEG-fMRI areas was concordant with an interictally active ECoG area, always including the early onset area of IEDs in the ECoG data. This confirms that EEG-fMRI reflects a pattern of onset and propagation of epileptic activity. At group level, 76% of the BOLD regions that were covered with subdural grids, were concordant with interictally active ECoG electrodes. Due to limited spatial sampling, 51% of the BOLD regions were not covered with electrodes and could, therefore, not be validated. From an ECoG perspective it appeared that 29% of the interictally active ECoG regions were missed by EEG-fMRI and that 68% of the brain regions were correctly identified as inactive with EEG-fMRI. Furthermore, EEG-fMRI areas included the complete seizure onset zone in 83% and resected area in 93% of the data sets. No clear distinction was found between patients with a good or poor surgical outcome: in both patient groups, EEG-fMRI correlation patterns were found that were either focal or widespread. In conclusion, by comparison of EEG-fMRI with interictal invasive EEG over a relatively large patient population we were able to show that the EEG-fMRI correlation patterns are spatially accurate at the level of neurosurgical units (i.e. anatomical brain regions) and reflect the underlying network of IEDs. Therefore, we expect that EEG-fMRI can play an important role for the determination of the implantation strategy.

Recommendations for improved reproducibility of ADC derivation on behalf of the Elekta MRI-linac consortium image analysis working group.

BACKGROUND AND PURPOSE: The apparent diffusion coefficient (ADC), a potential imaging biomarker for radiotherapy response, needs to be reproducible before translation into clinical use. The aim of this study was to evaluate the multi-centre delineation- and calculation-related ADC variation and give recommendations to minimize it. MATERIALS AND METHODS: Nine centres received identical diffusion-weighted and anatomical magnetic resonance images of different cancerous tumours (adrenal gland, pelvic oligo metastasis, pancreas, and prostate). All centres delineated the gross tumour volume (GTV), clinical target volume (CTV), and viable tumour volume (VTV), and calculated ADCs using both their local calculation methods and each of the following calculation conditions: b-values 0-500 vs. 150-500 s/mm2, region-of-interest (ROI)-based vs. voxel-based calculation, and mean vs. median. ADC variation was assessed using the mean coefficient of variation across delineations (CVD) and calculation methods (CVC). Absolute ADC differences between calculation conditions were evaluated using Friedman's test. Recommendations for ADC calculation were formulated based on observations and discussions within the Elekta MRI-linac consortium image analysis working group. RESULTS: The median (range) CVD and CVC were 0.06 (0.02-0.32) and 0.17 (0.08-0.26), respectively. The ADC estimates differed 18% between b-value sets and 4% between ROI/voxel-based calculation (p-values < 0.01). No significant difference was observed between mean and median (p = 0.64). Aligning calculation conditions between centres reduced CVC to 0.04 (0.01-0.16). CVD was comparable between ROI types. CONCLUSION: Overall, calculation methods had a larger impact on ADC reproducibility compared to delineation. Based on the results, significant sources of variation were identified, which should be considered when initiating new studies, in particular multi-centre investigations.

A semi-automatic method to determine electrode positions and labels from gel artifacts in EEG/fMRI-studies.

The analysis of simultaneous EEG and fMRI data is generally based on the extraction of regressors of interest from the EEG, which are correlated to the fMRI data in a general linear model setting. In more advanced approaches, the spatial information of EEG is also exploited by assuming underlying dipole models. In this study, we present a semi automatic and efficient method to determine electrode positions from electrode gel artifacts, facilitating the integration of EEG and fMRI in future EEG/fMRI data models. In order to visualize all electrode artifacts simultaneously in a single view, a surface rendering of the structural MRI is made using a skin triangular mesh model as reference surface, which is expanded to a "pancake view". Then the electrodes are determined with a simple mouse click for each electrode. Using the geometry of the skin surface and its transformation to the pancake view, the 3D coordinates of the electrodes are reconstructed in the MRI coordinate frame. The electrode labels are attached to the electrode positions by fitting a template grid of the electrode cap in which the labels are known. The correspondence problem between template and sample electrodes is solved by minimizing a cost function over rotations, shifts and scalings of the template grid. The crucial step here is to use the solution of the so-called "Hungarian algorithm" as a cost function, which makes it possible to identify the electrode artifacts in arbitrary order. The template electrode grid has to be constructed only once for each cap configuration. In our implementation of this method, the whole procedure can be performed within 15 min including import of MRI, surface reconstruction and transformation, electrode identification and fitting to template. The method is robust in the sense that an electrode template created for one subject can be used without identification errors for another subject for whom the same EEG cap was used. Furthermore, the method appears to be robust against spurious or missing artifacts. We therefore consider the proposed method as a useful and reliable tool within the larger toolbox required for the analysis of co-registered EEG/fMRI data.

Longitudinal Correlations Between Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced (DCE) MRI During Radiotherapy in Prostate Cancer Patients.

Purpose: Intravoxel incoherent motion (IVIM) is a promising technique that can acquire perfusion information without the use of contrast agent, contrary to the more established dynamic contrast-enhanced (DCE) technique. This is of interest for treatment response monitoring, where patients can be imaged on each treatment fraction. In this study, longitudinal correlations between IVIM- and DCE parameters were assessed in prostate cancer patients receiving radiation treatment. Materials and Methods: 20 prostate cancer patients were treated on a 1.5 T MR-linac with 20 x 3 or 3.1 Gy. Weekly IVIM and DCE scans were acquired. Tumors, the peripheral zone (PZ), and the transition zone (TZ) were delineated on a T2-weighted scan acquired on the first fraction. IVIM and DCE scans were registered to this scan and the delineations were propagated. Median values from these delineations were used for further analysis. The IVIM parameters D, f, D* and the product fD* were calculated. The Tofts model was used to calculate the DCE parameters Ktrans, kep and ve. Pearson correlations were calculated for the IVIM and DCE parameters on values from the first fraction for each region of interest (ROI). For longitudinal analysis, the repeated measures correlation coefficient was used to determine correlations between IVIM and DCE parameters in each ROI. Results: When averaging over patients, an increase during treatment in all IVIM and DCE parameters was observed in all ROIs, except for D in the PZ and TZ. No significant Pearson correlations were found between any pair of IVIM and DCE parameters measured on the first fraction. Significant but low longitudinal correlations were found for some combinations of IVIM and DCE parameters in the PZ and TZ, while no significant longitudinal correlations were found in the tumor. Notably in the TZ, for both f and fD*, significant longitudinal correlations with all DCE parameters were found. Conclusions: The increase in IVIM- and DCE parameters when averaging over patients indicates a measurable response to radiation treatment with both techniques. Although low, significant longitudinal correlations were found which suggests that IVIM could potentially be used as an alternative to DCE for treatment response monitoring.

Visualising spatial heterogeneity in glioblastoma using imaging habitats.

Glioblastoma is a high-grade aggressive neoplasm characterised by significant intra-tumoral spatial heterogeneity. Personalising therapy for this tumour requires non-invasive tools to visualise its heterogeneity to monitor treatment response on a regional level. To date, efforts to characterise glioblastoma's imaging features and heterogeneity have focussed on individual imaging biomarkers, or high-throughput radiomic approaches that consider a vast number of imaging variables across the tumour as a whole. Habitat imaging is a novel approach to cancer imaging that identifies tumour regions or 'habitats' based on shared imaging characteristics, usually defined using multiple imaging biomarkers. Habitat imaging reflects the evolution of imaging biomarkers and offers spatially preserved assessment of tumour physiological processes such perfusion and cellularity. This allows for regional assessment of treatment response to facilitate personalised therapy. In this review, we explore different methodologies to derive imaging habitats in glioblastoma, strategies to overcome its technical challenges, contrast experiences to other cancers, and describe potential clinical applications.