VesselBoost: A Python Toolbox for Small Blood Vessel Segmentation in Human Magnetic Resonance Angiography Data.

Magnetic resonance angiography (MRA) performed at ultra-high magnetic field provides a unique opportunity to study the arteries of the living human brain at the mesoscopic level. From this, we can gain new insights into the brain's blood supply and vascular disease affecting small vessels. However, for quantitative characterization and precise representation of human angioarchitecture to, for example, inform blood-flow simulations, detailed segmentations of the smallest vessels are required. Given the success of deep learning-based methods in many segmentation tasks, we here explore their application to high-resolution MRA data, and address the difficulty of obtaining large data sets of correctly and comprehensively labelled data. We introduce VesselBoost, a vessel segmentation package, which utilizes deep learning and imperfect training labels for accurate vasculature segmentation. Combined with an innovative data augmentation technique, which leverages the resemblance of vascular structures, VesselBoost enables detailed vascular segmentations.

Three-dimensional EPI with shot-selective CAIPIRIHANA for rapid high-resolution quantitative susceptibility mapping at 3 T.

PURPOSE: QSM provides insight into healthy brain aging and neuropathologies such as multiple sclerosis (MS), traumatic brain injuries, brain tumors, and neurodegenerative diseases. Phase data for QSM are usually acquired from 3D gradient-echo (3D GRE) scans with long acquisition times that are detrimental to patient comfort and susceptible to patient motion. This is particularly true for scans requiring whole-brain coverage and submillimeter resolutions. In this work, we use a multishot 3D echo plannar imaging (3D EPI) sequence with shot-selective 2D CAIPIRIHANA to acquire high-resolution, whole-brain data for QSM with minimal distortion and blurring. METHODS: To test clinical viability, the 3D EPI sequence was used to image a cohort of MS patients at 1-mm isotropic resolution at 3 T. Additionally, 3D EPI data of healthy subjects were acquired at 1-mm, 0.78-mm, and 0.65-mm isotropic resolution with varying echo train lengths (ETLs) and compared with a reference 3D GRE acquisition. RESULTS: The appearance of the susceptibility maps and the susceptibility values for segmented regions of interest were comparable between 3D EPI and 3D GRE acquisitions for both healthy and MS participants. Additionally, all lesions visible in the MS patients on the 3D GRE susceptibility maps were also visible on the 3D EPI susceptibility maps. The interplay among acquisition time, resolution, echo train length, and the effect of distortion on the calculated susceptibility maps was investigated. CONCLUSION: We demonstrate that the 3D EPI sequence is capable of rapidly acquiring submillimeter resolutions and providing high-quality, clinically relevant susceptibility maps.

Adiposity associates with lower plasma resolvin E1 (Rve1): a population study.

BACKGROUND: Inadequate inflammation resolution may contribute to persistent low-grade inflammation that accompanies many chronic conditions. Resolution of inflammation is an active process driven by Specialized Pro-resolving Mediators (SPM) that derive from long chain n-3 and n-6 fatty acids. This study examined plasma SPM in relation to sex differences, lifestyle and a broad range cardiovascular disease (CVD) risk factors in 978, 27-year olds from the Australian Raine Study. METHODS: Plasma SPM pathway intermediates (18-HEPE, 17-HDHA and 14-HDHA), and SPM (E- and D-series resolvins, PD1, MaR1) and LTB4 were measured by liquid chromatography-tandem mass spectrometry (LCMSMS). Pearson correlations and multiple regression analyses assessed relationships between SPM and CVD risk factors. Unpaired t-tests or ANOVA assessed the effect of sex, smoking, unhealthy alcohol consumption and obesity on SPM. RESULTS: Women had higher 17-HDHA (p = 0.01) and lower RvE1 (p < 0.0001) and RvD1 (p = 0.05) levels compared with men. In univariate analysis, obesity associated with lower RvE1 (p = 0.002), whereas smoking (p < 0.001) and higher alcohol consumption (p < 0.001) associated with increased RvE1. In multiple regression analysis, plasma RvE1 was negatively associated with a range of measures of adiposity including BMI, waist circumference, waist-to-height ratio, abdominal subcutaneous fat volume, and skinfold thicknesses in both men and women. CONCLUSION: This population study suggests that a deficiency in plasma RvE1 may occur in response to increasing adiposity. This observation could be relevant to ongoing inflammation that associates with CVD and other chronic diseases.

Improving quantitative susceptibility mapping for the identification of traumatic brain injury neurodegeneration at the individual level.

BACKGROUND: Emerging evidence suggests that traumatic brain injury (TBI) is a major risk factor for developing neurodegenerative disease later in life. Quantitative susceptibility mapping (QSM) has been used by an increasing number of studies in investigations of pathophysiological changes in TBI. However, generating artefact-free quantitative susceptibility maps in brains with large focal lesions, as in the case of moderate-to-severe TBI (ms-TBI), is particularly challenging. To address this issue, we utilized a novel two-pass masking technique and reconstruction procedure (two-pass QSM) to generate quantitative susceptibility maps (QSMxT; Stewart et al., 2022, Magn Reson Med.) in combination with the recently developed virtual brain grafting (VBG) procedure for brain repair (Radwan et al., 2021, NeuroImage) to improve automated delineation of brain areas. We used QSMxT and VBG to generate personalised QSM profiles of individual patients with reference to a sample of healthy controls. METHODS: Chronic ms-TBI patients (N = 8) and healthy controls (N = 12) underwent (multi-echo) GRE, and anatomical MRI (MPRAGE) on a 3T Siemens PRISMA scanner. We reconstructed the magnetic susceptibility maps using two-pass QSM from QSMxT. We then extracted values of magnetic susceptibility in grey matter (GM) regions (following brain repair via VBG) across the whole brain and determined if they deviate from a reference healthy control group [Z-score < -3.43 or > 3.43, relative to the control mean], with the aim of obtaining personalised QSM profiles. RESULTS: Using two-pass QSM, we achieved susceptibility maps with a substantial increase in quality and reduction in artefacts irrespective of the presence of large focal lesions, compared to single-pass QSM. In addition, VBG minimised the loss of GM regions and exclusion of patients due to failures in the region delineation step. Our findings revealed deviations in magnetic susceptibility measures from the HC group that differed across individual TBI patients. These changes included both increases and decreases in magnetic susceptibility values in multiple GM regions across the brain. CONCLUSIONS: We illustrate how to obtain magnetic susceptibility values at the individual level and to build personalised QSM profiles in ms-TBI patients. Our approach opens the door for QSM investigations in more severely injured patients. Such profiles are also critical to overcome the inherent heterogeneity of clinical populations, such as ms-TBI, and to characterize the underlying mechanisms of neurodegeneration at the individual level more precisely. Moreover, this new personalised QSM profiling could in the future assist clinicians in assessing recovery and formulating a neuroscience-guided integrative rehabilitation program tailored to individual TBI patients.

Super-resolution QSM in little or no additional time for imaging (NATIve) using 2D EPI imaging in 3 orthogonal planes.

Quantitative Susceptibility Mapping has the potential to provide additional insights into neurological diseases but is typically based on a quite long (5-10 min) 3D gradient-echo scan which is highly sensitive to motion. We propose an ultra-fast acquisition based on three orthogonal (sagittal, coronal and axial) 2D simultaneous multi-slice EPI scans with 1 mm in-plane resolution and 3 mm thick slices. Images in each orientation are corrected for susceptibility-related distortions and co-registered with an iterative non-linear Minimum Deformation Averaging (Volgenmodel) approach to generate a high SNR, super-resolution data set with an isotropic resolution of close to 1 mm. The net acquisition time is 3 times the volume acquisition time of EPI or about 12 s, but the three volumes could also replace "dummy scans" in fMRI, making it feasible to acquire QSM in little or No Additional Time for Imaging (NATIve). NATIve QSM values agreed well with reference 3D GRE QSM in the basal ganglia in healthy subjects. In patients with multiple sclerosis, there was also a good agreement between the susceptibility values within lesions and control ROIs and all lesions which could be seen on 3D GRE QSMs could also be visualized on NATIve QSMs. The approach is faster than conventional 3D GRE by a factor of 25-50 and faster than 3D EPI by a factor of 3-5. As a 2D technique, NATIve QSM was shown to be much more robust to motion than the 3D GRE and 3D EPI, opening up the possibility of studying neurological diseases involving iron accumulation and demyelination in patients who find it difficult to lie still for long enough to acquire QSM data with conventional methods.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL).

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

Variability of visual field maps in human early extrastriate cortex challenges the canonical model of organization of V2 and V3.

Visual field maps in human early extrastriate areas (V2 and V3) are traditionally thought to form mirror-image representations which surround the primary visual cortex (V1). According to this scheme, V2 and V3 form nearly symmetrical halves with respect to the calcarine sulcus, with the dorsal halves representing lower contralateral quadrants, and the ventral halves representing upper contralateral quadrants. This arrangement is considered to be consistent across individuals, and thus predictable with reasonable accuracy using templates. However, data that deviate from this expected pattern have been observed, but mainly treated as artifactual. Here, we systematically investigate individual variability in the visual field maps of human early visual cortex using the 7T Human Connectome Project (HCP) retinotopy dataset. Our results demonstrate substantial and principled inter-individual variability. Visual field representation in the dorsal portions of V2 and V3 was more variable than in their ventral counterparts, including substantial departures from the expected mirror-symmetrical patterns. In addition, left hemisphere retinotopic maps were more variable than those in the right hemisphere. Surprisingly, only one-third of individuals had maps that conformed to the expected pattern in the left hemisphere. Visual field sign analysis further revealed that in many individuals the area conventionally identified as dorsal V3 shows a discontinuity in the mirror-image representation of the retina, associated with a Y-shaped lower vertical representation. Our findings challenge the current view that inter-individual variability in early extrastriate cortex is negligible, and that the dorsal portions of V2 and V3 are roughly mirror images of their ventral counterparts.

Neurodesk: An accessible, flexible, and portable data analysis environment for reproducible neuroimaging.

Neuroimaging data analysis often requires purpose-built software, which can be challenging to install and may produce different results across computing environments. Beyond being a roadblock to neuroscientists, these issues of accessibility and portability can hamper the reproducibility of neuroimaging data analysis pipelines. Here, we introduce the Neurodesk platform, which harnesses software containers to support a comprehensive and growing suite of neuroimaging software (https://www.neurodesk.org/). Neurodesk includes a browser-accessible virtual desktop environment and a command line interface, mediating access to containerized neuroimaging software libraries on various computing platforms, including personal and high-performance computers, cloud computing and Jupyter Notebooks. This community-oriented, open-source platform enables a paradigm shift for neuroimaging data analysis, allowing for accessible, flexible, fully reproducible, and portable data analysis pipelines.

NeXtQSM-A complete deep learning pipeline for data-consistent Quantitative Susceptibility Mapping trained with hybrid data.

Deep learning based Quantitative Susceptibility Mapping (QSM) has shown great potential in recent years, obtaining similar results to established non-learning approaches. Many current deep learning approaches are not data consistent, require in vivo training data or solve the QSM problem in consecutive steps resulting in the propagation of errors. Here we aim to overcome these limitations and developed a framework to solve the QSM processing steps jointly. We developed a new hybrid training data generation method that enables the end-to-end training for solving background field correction and dipole inversion in a data-consistent fashion using a variational network that combines the QSM model term and a learned regularizer. We demonstrate that NeXtQSM overcomes the limitations of previous deep learning methods. NeXtQSM offers a new deep learning based pipeline for computing quantitative susceptibility maps that integrates each processing step into the training and provides results that are robust and fast.

Spikes with and without concurrent high-frequency oscillations: Topographic relationship and neural correlates using EEG-fMRI.

OBJECTIVE: Epilepsy surgery is the best therapeutic option for patients with drug-resistant focal epilepsy. During presurgical investigation, interictal spikes can provide important information on eligibility, lateralisation and localisation of the surgical target. However, their relationship to epileptogenic tissue is variable. Interictal spikes with concurrent high-frequency oscillations (HFOs) have been postulated to reflect epileptogenic tissue more reliably. Here, we studied the voltage distribution of scalp-recorded spikes with and without concurrent HFO and identified their respective haemodynamic correlates using simultaneous electroencephalography and functional Magnetic Resonance Imaging (EEG-fMRI). METHODS: The scalp topography of spikes with and without concurrent HFOs were assessed in 31 consecutive patients with focal drug-resistant epilepsy who showed interictal spikes during presurgical evaluation. Simultaneous EEG-fMRI was then used in 17 patients with spikes and concurrent HFOs. Haemodynamic changes were obtained from the spatial correlation between the patient-specific voltage map of each spike population and the intra-scanner EEG. The haemodynamic response of spikes with and without HFOs were compared in terms of their spatial similarity, strength, the distance between activation peaks and concordance with interictal localisation. RESULTS: Twenty-five patients showed spikes with and without concurrent HFOs. Among patients with both types of spikes, most spikes were not associated with HFOs (p < 0.0001, Mann-Whitney test). Twenty of the 25 patients showed an average of 8 ± 6 (standard deviation) electrodes with significant voltage differences (p = 0.025, permutation test corrected for multiple comparisons) on scalp electrodes within and distant to the spike field. Comparing the haemodynamic response between both spike populations, we found no significant differences in the peak strength (p = 0.71, Mann-Whitney test), spatial distribution (p = 0.113, One-sample Wilcoxon test) and distance between activation peaks (p = 0.5, One-sample Wilcoxon test), with all peaks being co-localised in the same lobe. SIGNIFICANCE: Our data showed that spikes with and without HFOs have different scalp voltage distributions. However, when assessing the haemodynamic changes of each spike type, we found that both elicit similar haemodynamic changes and share high spatial similarity suggesting that the epileptic networks of spikes with and without HFOs have the same underlying neural substrate.

Deep learning-based quantitative susceptibility mapping (QSM) in the presence of fat using synthetically generated multi-echo phase training data.

PURPOSE: To enable a fast and automatic deep learning-based QSM reconstruction of tissues with diverse chemical shifts, relevant to most regions outside the brain. METHODS: A UNET was trained to reconstruct susceptibility maps using synthetically generated, unwrapped, multi-echo phase data as input. The RMS error with respect to synthetic validation data was computed. The method was tested on two in vivo knee and two pelvis data sets. Comparisons were made to a conventional fat-water separation pipeline by applying a commonly used graph-cut algorithm, both without and with an extended mask for background field removal (FWS-CONV-QSM and FWS-MASK-CONV-QSM, respectively). Several regions of interest were segmented and compared. Furthermore, the approach was tested on a prostate cancer patient receiving low-dose-rate brachytherapy, to detect and localize the seeds by MRI. RESULTS: The RMS error was 0.292 ppm with FWS-CONV-QSM and 0.123 ppm for the UNET approach. Susceptibility maps were reconstructed much faster (< 10 s) and completely automatically (no background masking needed) by the UNET compared with the other applied techniques (5 min 51 s and 22 min 44 s for CONV-QSM and FWS-MASK-CONV-QSM, respectively. Background artifacts, fat-water swaps, and hypointense artifacts between I-125 seeds of a patient receiving low-dose brachytherapy in the prostate were largely reduced in the UNET approach. CONCLUSIONS: Deep learning-based QSM reconstruction, trained solely with synthetic data, is well-suited to rapidly reconstructing high-quality susceptibility maps in the presence of fat without needing masking for background field removal.

Overview of quantitative susceptibility mapping using deep learning: Current status, challenges and opportunities.

Quantitative susceptibility mapping (QSM) has gained broad interest in the field by extracting bulk tissue magnetic susceptibility, predominantly determined by myelin, iron and calcium from magnetic resonance imaging (MRI) phase measurements in vivo. Thereby, QSM can reveal pathological changes of these key components in a variety of diseases. QSM requires multiple processing steps such as phase unwrapping, background field removal and field-to-source inversion. Current state-of-the-art techniques utilize iterative optimization procedures to solve the inversion and background field correction, which are computationally expensive and require a careful choice of regularization parameters. With the recent success of deep learning using convolutional neural networks for solving ill-posed reconstruction problems, the QSM community also adapted these techniques and demonstrated that the QSM processing steps can be solved by efficient feed forward multiplications not requiring either iterative optimization or the choice of regularization parameters. Here, we review the current status of deep learning-based approaches for processing QSM, highlighting limitations and potential pitfalls, and discuss the future directions the field may take to exploit the latest advances in deep learning for QSM.

QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping.

PURPOSE: Quantitative susceptibility mapping (QSM) estimates the spatial distribution of tissue magnetic susceptibilities from the phase of a gradient-echo signal. QSM algorithms require a signal mask to delineate regions with reliable phase for subsequent susceptibility estimation. Existing masking techniques used in QSM have limitations that introduce artifacts, exclude anatomical detail, and rely on parameter tuning and anatomical priors that narrow their application. Here, a robust masking and reconstruction procedure is presented to overcome these limitations and enable automated QSM processing. Moreover, this method is integrated within an open-source software framework: QSMxT. METHODS: A robust masking technique that automatically separates reliable from less reliable phase regions was developed and combined with a two-pass reconstruction procedure that operates on the separated sources before combination, extracting more information and suppressing streaking artifacts. RESULTS: Compared with standard masking and reconstruction procedures, the two-pass inversion reduces streaking artifacts caused by unreliable phase and high dynamic ranges of susceptibility sources. It is also robust across a range of acquisitions at 3 T in volunteers and phantoms, at 7 T in tumor patients, and in an in silico head phantom, with significant artifact and error reductions, greater anatomical detail, and minimal parameter tuning. CONCLUSION: The two-pass masking and reconstruction procedure separates reliable from less reliable phase regions, enabling a more accurate QSM reconstruction that mitigates artifacts, operates without anatomical priors, and requires minimal parameter tuning. The technique and its integration within QSMxT makes QSM processing more accessible and robust to streaking artifacts.

Predicting the retinotopic organization of human visual cortex from anatomy using geometric deep learning.

Whether it be in a single neuron or a more complex biological system like the human brain, form and function are often directly related. The functional organization of human visual cortex, for instance, is tightly coupled with the underlying anatomy with cortical shape having been shown to be a useful predictor of the retinotopic organization in early visual cortex. Although the current state-of-the-art in predicting retinotopic maps is able to account for gross individual differences, such models are unable to account for any idiosyncratic differences in the structure-function relationship from anatomical information alone due to their initial assumption of a template. Here we developed a geometric deep learning model capable of exploiting the actual structure of the cortex to learn the complex relationship between brain function and anatomy in human visual cortex such that more realistic and idiosyncratic maps could be predicted. We show that our neural network was not only able to predict the functional organization throughout the visual cortical hierarchy, but that it was also able to predict nuanced variations across individuals. Although we demonstrate its utility for modeling the relationship between structure and function in human visual cortex, our approach is flexible and well-suited for a range of other applications involving data structured in non-Euclidean spaces.

Centering inclusivity in the design of online conferences-An OHBM-Open Science perspective.

As the global health crisis unfolded, many academic conferences moved online in 2020. This move has been hailed as a positive step towards inclusivity in its attenuation of economic, physical, and legal barriers and effectively enabled many individuals from groups that have traditionally been underrepresented to join and participate. A number of studies have outlined how moving online made it possible to gather a more global community and has increased opportunities for individuals with various constraints, e.g., caregiving responsibilities. Yet, the mere existence of online conferences is no guarantee that everyone can attend and participate meaningfully. In fact, many elements of an online conference are still significant barriers to truly diverse participation: the tools used can be inaccessible for some individuals; the scheduling choices can favour some geographical locations; the set-up of the conference can provide more visibility to well-established researchers and reduce opportunities for early-career researchers. While acknowledging the benefits of an online setting, especially for individuals who have traditionally been underrepresented or excluded, we recognize that fostering social justice requires inclusivity to actively be centered in every aspect of online conference design. Here, we draw from the literature and from our own experiences to identify practices that purposefully encourage a diverse community to attend, participate in, and lead online conferences. Reflecting on how to design more inclusive online events is especially important as multiple scientific organizations have announced that they will continue offering an online version of their event when in-person conferences can resume.

Deep learning in magnetic resonance image reconstruction.

Magnetic resonance (MR) imaging visualises soft tissue contrast in exquisite detail without harmful ionising radiation. In this work, we provide a state-of-the-art review on the use of deep learning in MR image reconstruction from different image acquisition types involving compressed sensing techniques, parallel image acquisition and multi-contrast imaging. Publications with deep learning-based image reconstruction for MR imaging were identified from the literature (PubMed and Google Scholar), and a comprehensive description of each of the works was provided. A detailed comparison that highlights the differences, the data used and the performance of each of these works were also made. A discussion of the potential use cases for each of these methods is provided. The sparse image reconstruction methods were found to be most popular in using deep learning for improved performance, accelerating acquisitions by around 4-8 times. Multi-contrast image reconstruction methods rely on at least one pre-acquired image, but can achieve 16-fold, and even up to 32- to 50-fold acceleration depending on the set-up. Parallel imaging provides frameworks to be integrated in many of these methods for additional speed-up potential. The successful use of compressed sensing techniques and multi-contrast imaging with deep learning and parallel acquisition methods could yield significant MR acquisition speed-ups within clinical routines in the near future.

Improving FLAIR SAR efficiency at 7T by adaptive tailoring of adiabatic pulse power through deep learning B1+ estimation.

PURPOSE: The purpose of this study is to demonstrate a method for specific absorption rate (SAR) reduction for 2D T2 -FLAIR MRI sequences at 7 T by predicting the required adiabatic radiofrequency (RF) pulse power and scaling the RF amplitude in a slice-wise fashion. METHODS: We used a time-resampled frequency-offset corrected inversion (TR-FOCI) adiabatic pulse for spin inversion in a T2 -FLAIR sequence to improve B1+ homogeneity and calculated the pulse power required for adiabaticity slice-by-slice to minimize the SAR. Drawing on the implicit B1+ inhomogeneity in a standard localizer scan, we acquired 3D AutoAlign localizers and SA2RAGE B1+ maps in 28 volunteers. Then, we trained a convolutional neural network (CNN) to estimate the B1+ profile from the localizers and calculated pulse scale factors for each slice. We assessed the predicted B1+ profiles and the effect of scaled pulse amplitudes on the FLAIR inversion efficiency in oblique transverse, sagittal, and coronal orientations. RESULTS: The predicted B1+ amplitude maps matched the measured ones with a mean difference of 9.5% across all slices and participants. The slice-by-slice scaling of the TR-FOCI inversion pulse was most effective in oblique transverse orientation and resulted in a 1 min and 30 s reduction in SAR induced delay time while delivering identical image quality. CONCLUSION: We propose a SAR reduction technique based on the estimation of B1+ profiles from standard localizer scans using a CNN and show that scaling the inversion pulse power slice-by-slice for FLAIR sequences at 7T reduces SAR and scan time without compromising image quality.

Functional connectivity of the irritative zone identified by electrical source imaging, and EEG-correlated fMRI analyses.

OBJECTIVE: The irritative zone - the area generating epileptic spikes - can be studied non-invasively during the interictal period using Electrical Source Imaging (ESI) and simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI). Although the techniques yield results which may overlap spatially, differences in spatial localization of the irritative zone within the same patient are consistently observed. To investigate this discrepancy, we used Blood Oxygenation Level Dependent (BOLD) functional connectivity measures to examine the underlying relationship between ESI and EEG-fMRI findings. METHODS: Fifteen patients (age 20-54), who underwent presurgical epilepsy investigation, were scanned using a single-session resting-state EEG-fMRI protocol. Structural MRI was used to obtain the electrode localisation of a high-density 64-channel EEG cap. Electrical generators of interictal epileptiform discharges were obtained using a distributed local autoregressive average (LAURA) algorithm as implemented in Cartool EEG software. BOLD activations were obtained using both spike-related and voltage-map EEG-fMRI analysis. The global maxima of each method were used to investigate the temporal relationship of BOLD time courses and to assess the spatial similarity using the Dice similarity index between functional connectivity maps. RESULTS: ESI, voltage-map and spike-related EEG-fMRI methods identified peaks in 15 (100%), 13 (67%) and 8 (53%) of the 15 patients, respectively. For all methods, maxima were localised within the same lobe, but differed in sub-lobar localisation, with a median distance of 22.8 mm between the highest peak for each method. The functional connectivity analysis showed that the temporal correlation between maxima only explained 38% of the variance between the time course of the BOLD response at the maxima. The mean Dice similarity index between seed-voxel functional connectivity maps showed poor spatial agreement. SIGNIFICANCE: Non-invasive methods for the localisation of the irritative zone have distinct spatial and temporal sensitivity to different aspects of the local cortical network involved in the generation of interictal epileptiform discharges.

MRI phase offset correction method impacts quantitative susceptibility mapping.

Individual channel ultra-high field (7T) phase images have to be phase offset corrected prior to the mapping of magnetic susceptibility of tissue. Whilst numerous methods have been proposed for gradient recalled echo MRI phase offset correction, it remains unclear how they affect quantitative magnetic susceptibility values derived from phase images. Methods already proposed either employ a single or multiple echo time MRI data. In terms of the latter, offsets can be derived using an ultra-short echo time acquisition, or by estimating the offset based on two echo points with the assumption of linear phase evolution with echo time. Our evaluation involved 32 channel multi-echo time 7T GRE (Gradient Recalled Echo) and ultra-short echo time PETRA (Pointwise Encoding Time Reduction with Radial Acquisition) MRI data collected for a susceptibility phantom and three human brains. The combined phase images generated using four established offset correction methods (two single and two multiple echo time) were analysed, followed by an assessment of quantitative susceptibility values obtained for a phantom and human brains. The effectiveness of each method in removing the offsets was shown to reduce with increased echo time, decreased signal intensity and reduced overlap in coil sensitivity profiles. Quantitative susceptibility values and how they change with echo time were found to be method specific. Phase offset correction methods based on single echo time data have a tendency to produce more accurate and less noisy quantitative susceptibility maps in comparison with methods employing multiple echo time data.