ICA-based denoising for ASL perfusion imaging

Arterial Spin Labelling (ASL) imaging derives a perfusion image by tracing the accumulation of magnetically labeled blood water in the brain. As the image generated has an intrinsically low signal to noise ratio (SNR), multiple measurements are routinely acquired and averaged, at a penalty of increased scan duration and opportunity for motion artefact. However, this strategy alone might be ineffective in clinical settings where the time available for acquisition is limited and patient motion are increased. This study investigates the use of an Independent Component Analysis (ICA) approach for denoising ASL data, and its potential for automation.72 ASL datasets (pseudo-continuous ASL; 5 different post-labeling delays: 400, 800, 1200, 1600, 2000 m s; total volumes = 60) were collected from thirty consecutive acute stroke patients. The effects of ICA-based denoising (manual and automated) where compared to two different denoising approaches, aCompCor, a Principal Component-based method, and Enhancement of Automated Blood Flow Estimates (ENABLE), an algorithm based on the removal of corrupted volumes. Multiple metrics were used to assess the changes in the quality of the data following denoising, including changes in cerebral blood flow (CBF) and arterial transit time (ATT), SNR, and repeatability. Additionally, the relationship between SNR and number of repetitions acquired was estimated before and after denoising the data.The use of an ICA-based denoising approach resulted in significantly higher mean CBF and ATT values (p [less than] 0.001), lower CBF and ATT variance (p [less than] 0.001), increased SNR (p [less than] 0.001), and improved repeatability (p [less than] 0.05) when compared to the raw data. The performance of manual and automated ICA-based denoising was comparable. These results went beyond the effects of aCompCor or ENABLE. Following ICA-based denoising, the SNR was higher using only 50% of the ASL-dataset collected than when using the whole raw data.The results show that ICA can be used to separate signal from noise in ASL data, improving the quality of the data collected. In fact, this study suggests that the acquisition time could be reduced by 50% without penalty to data quality, something that merits further study. Independent component classification and regression can be carried out either manually, following simple criteria, or automatically

Multi-parametric Imaging Using Hybrid PET/MR to Investigate the Epileptogenic Brain

Neuroimaging analysis has led to fundamental discoveries about the healthy and pathological human brain. Different imaging modalities allow garnering complementary information about brain metabolism, structure and function. To ensure that the integration of imaging data from these modalities is robust and reliable, it is fundamental to attain deep knowledge of each modality individually. Epilepsy, a neurological condition characterised by recurrent spontaneous seizures, represents a field in which applications of neuroimaging and multi-parametric imaging are particularly promising to guide diagnosis and treatment. In this PhD thesis, I focused on different imaging modalities and investigated advanced denoising and analysis strategies to improve their application to epilepsy. The first project focused on fluorodeoxyglucose (FDG) positron emission tomography (PET), a well-established imaging modality assessing brain metabolism, and aimed to develop a novel, semi-quantitative pipeline to analyse data in children with epilepsy, thus aiding presurgical planning. As pipelines for FDG-PET analysis in children are currently lacking, I developed age-appropriate templates to provide statistical parametric maps identifying epileptogenic areas on patient scans. The second and third projects focused on two magnetic resonance imaging (MRI) modalities: resting-state functional MRI (rs-fMRI) and arterial spin labelling (ASL), respectively. The aim was to i) probe the efficacy of different fMRI denoising pipelines, and ii) formally compare different ASL data acquisition strategies. In the former case, I compared different pre-processing methods and assessed their impact on fMRI signal quality and related functional connectivity analyses. In the latter case, I compared two ASL sequences to investigate their ability to quantify cerebral blood flow and interregional brain connectivity. The final project addressed the combination of rs-fMRI and ASL, and leveraged graph-theoretical analysis tools to i) compare metrics estimated via these two imaging modalities in healthy subjects and ii) assess topological changes captured by these modalities in a sample of temporal lobe epilepsy patients

ExploreASL: An image processing pipeline for multi-center ASL perfusion MRI studies

Arterial spin labeling (ASL) has undergone significant development since its inception, with a focus on improving standardization and reproducibility of its acquisition and quantification. In a community-wide effort towards robust and reproducible clinical ASL image processing, we developed the software package ExploreASL, allowing standardized analyses across centers and scanners. The procedures used in ExploreASL capitalize on published image processing advancements and address the challenges of multi-center datasets with scanner-specific processing and artifact reduction to limit patient exclusion. ExploreASL is self-contained, written in MATLAB and based on Statistical Parameter Mapping (SPM) and runs on multiple operating systems. To facilitate collaboration and data-exchange, the toolbox follows several standards and recommendations for data structure, provenance, and best analysis practice. ExploreASL was iteratively refined and tested in the analysis of >10,000 ASL scans using different pulse-sequences in a variety of clinical populations, resulting in four processing modules: Import, Structural, ASL, and Population that perform tasks, respectively, for data curation, structural and ASL image processing and quality control, and finally preparing the results for statistical analyses on both single-subject and group level. We illustrate ExploreASL processing results from three cohorts: perinatally HIV-infected children, healthy adults, and elderly at risk for neurodegenerative disease. We show the reproducibility for each cohort when processed at different centers with different operating systems and MATLAB versions, and its effects on the quantification of gray matter cerebral blood flow. ExploreASL facilitates the standardization of image processing and quality control, allowing the pooling of cohorts which may increase statistical power and discover between-group perfusion differences. Ultimately, this workflow may advance ASL for wider adoption in clinical studies, trials, and practice

Improved quantification of perfusion in patients with cerebrovascular disease.

In recent years measurements of cerebral perfusion using bolus-tracking MRI have become common clinical practice in the diagnosis and management of patients with stroke and cerebrovascular disease. An active area of research is the development of methods to identify brain tissue that is at risk of irreversible damage, but amenable to salvage using reperfusion treatments, such as thrombolysis. However, the specificity and sensitivity of these methods are limited by the inaccuracies in the perfusion data. Accurate measurements of perfusion are difficult to obtain, especially in patients with cerebrovascular diseases. In particular, if the bolus of MR contrast is delayed and/or dispersed due to cerebral arterial abnormalities, perfusion is likely to be underestimated using the standard analysis techniques. The potential for such underestimation is often overlooked when using the perfusion maps to assess stroke patients. Since thrombolysis can increase the risk of haemorrhage, a misidentification of 'at-risk' tissue has potentially dangerous clinical implications. This thesis presents several methodologies which aim to improve the accuracy and interpretation of the analysed bolus-tracking data. Two novel data analysis techniques are proposed, which enable the identification of brain regions where delay and dispersion of the bolus are likely to bias the perfusion measurements. In this way true hypoperfusion can be distinguished from erroneously low perfusion estimates. The size of the perfusion measurement errors are investigated in vivo, and a parameterised characterisation of the bolus delay and dispersion is obtained. Such information is valuable for the interpretation of in vivo data, and for further investigation into the effects of abnormal vasculature on perfusion estimates. Finally, methodology is presented to minimise the perfusion measurement errors prevalent in patients with cerebrovascular diseases. The in vivo application of this method highlights the dangers of interpreting perfusion values independently of the bolus delay and dispersion

Arterial Spin Labeling Reveals Disrupted Brain Networks and Functional Connectivity in Drug-Resistant Temporal Epilepsy

Resting-state networks (RSNs) and functional connectivity (FC) have been increasingly exploited for mapping brain activity and identifying abnormalities in pathologies, including epilepsy. The majority of studies currently available are based on bloodoxygenation- level-dependent (BOLD) contrast in combination with either independent component analysis (ICA) or pairwise region of interest (ROI) correlations. Despite its success, this approach has several shortcomings as BOLD is only an indirect and non-quantitative measure of brain activity. Conversely, promising results have recently been achieved by arterial spin labeling (ASL) MRI, primarily developed to quantify brain perfusion. However, the wide application of ASL-based FC has been hampered by its complexity and relatively low robustness to noise, leaving several aspects of this approach still largely unexplored. In this study, we firstly aimed at evaluating the effect of noise reduction on spatio-temporal ASL analyses and quantifying the impact of two ad-hoc processing pipelines (basic and advanced) on connectivity measures. Once the optimal strategy had been defined, we investigated the applicability of ASL for connectivity mapping in patients with drug-resistant temporal epilepsy vs. controls (10 per group), aiming at revealing between-group voxel-wise differences in each RSN and ROI-wise FC changes. We first found ASL was able to identify the main network (DMN) along with all the others generally detected with BOLD but never previously reported from ASL. For all RSNs, ICA-based denoising (advanced pipeline) allowed to increase their similarity with the corresponding BOLD template. ASL-based RSNs were visibly consistent with literature findings; however, group differences could be identified in the structure of some networks. Indeed, statistics revealed areas of significant FC decrease in patients within different RSNs, such as DMN and cerebellum (CER), while significant increases were found in some cases, such as the visual networks. Finally, the ROI-based analyses identified several inter-hemispheric dysfunctional links (controls > patients) mainly between areas belonging to the DMN, right-left thalamus and right-left temporal lobe. Conversely, fewer connections, predominantly intra-hemispheric, showed the opposite pattern (controls < patients). All these elements provide novel insights into the pathological modulations characterizing a “network disease” as epilepsy, shading light on the importance of perfusion-based approaches for identifying the disrupted areas and communications between brain regions

Opioid suppression of conditioned anticipatory brain responses to breathlessness

Opioid painkillers are a promising treatment for chronic breathlessness, but are associated with potentially fatal side effects. In the treatment of breathlessness, their mechanisms of action are unclear. A better understanding might help to identify safer alternatives. Learned associations between previously neutral stimuli (e.g. stairs) and repeated breathlessness induce an anticipatory threat response that may worsen breathlessness, contributing to the downward spiral of decline seen in clinical populations. As opioids are known to influence associative learning, we hypothesized that they may interfere with the brain processes underlying a conditioned anticipatory response to breathlessness in relevant brain areas, including the amygdala and the hippocampus. Healthy volunteers viewed visual cues (neutral stimuli) immediately before induction of experimental breathlessness with inspiratory resistive loading. Thus, an association was formed between the cue and breathlessness. Subsequently, this paradigm was repeated in two identical neuroimaging sessions with intravenous infusions of either low-dose remifentanil (0.7ng/ml target controlled infusion) or saline (randomised). During saline infusion, breathlessness anticipation activated the right anterior insula and the adjacent operculum. Breathlessness was associated with activity in a network including the insula, operculum, dorsolateral prefrontal cortex, anterior cingulate cortex and the primary sensory and motor cortices. Remifentanil reduced breathlessness unpleasantness but not breathlessness intensity. Remifentanil depressed anticipatory activity in the amygdala and the hippocampus that correlated with reductions in breathlessness unpleasantness. During breathlessness, remifentanil decreased activity in the anterior insula, anterior cingulate cortex and sensory motor cortices. Remifentanil-induced reduction in breathlessness unpleasantness was associated with increased activity in the rostral anterior cingulate cortex and nucleus accumbens, components of the endogenous opioid system known to decrease the perception of aversive stimuli. These findings suggest that in addition to effects on brainstem respiratory control, opioids palliate breathlessness through an interplay of altered associative learning mechanisms. These mechanisms provide potential targets for novel ways to develop and assess treatments for chronic breathlessness

Application of resting-state fMRI methods to acute ischemic stroke

Diffusion weighted imaging (DWI) and dynamic susceptibility contrast-enhanced (DSC) perfusion-weighted imaging (PWI) are commonly employed in clinical practice and in research to give pathophysiological information for patients with acute ischemic stroke. DWI is thought to roughly reflect the severely damaged infarct core, while DSC-PWI reflects the area of hypoperfusion. The volumetric difference between DWI and DSC-PWI is termed the PWI/DWI-mismatch, and has been suggested as an MRI surrogate of the ischemic penumbra. However, due to the application of a contrast agent, which has potentially severe side-effects (e.g., nephrogenic systemic fibrosis), the DSC-PWI precludes repetitive examinations for monitoring purposes. New approaches are being sought to overcome this shortcoming. BOLD (blood oxygen-level dependent) signal can reflect the metabolism of blood oxygen in the brain and hemodynamics can be assessed with resting-state fMRI. The aim of this thesis was to use resting-state fMRI as a new approach to give similar information as DSC-PWI. This thesis comprises two studies: In the first study (see Chapter 2), two resting-state fMRI methods, local methods which compare low frequency amplitudes between two hemispheres and a k-means clustering approach, were applied to investigate the functional damage of patients with acute ischemic stroke both in the time domain and frequency domain. We found that the lesion areas had lower amplitudes than contralateral homotopic healthy tissues. We also differentiated the lesion areas from healthy tissues using a k-means clustering approach. In the second study (see Chapter 3), time-shift analysis (TSA), which assesses time delays of the spontaneous low frequency fluctuations of the resting-state BOLD signal, was applied to give similar pathophysiological information as DSC-PWI in the acute phase of stroke. We found that areas which showed a pronounced time delay to the respective mean time course were very similar to the hypoperfusion area. In summary, we suggest that the resting-state fMRI methods, especially the time-shift analysis (TSA), may provide comparable information to DSC-PWI and thus serve as a useful diagnostic tool for stroke MRI without the need for the application of a contrast agent

Intelligent Imaging of Perfusion Using Arterial Spin Labelling

Arterial spin labelling (ASL) is a powerful magnetic resonance imaging technique, which can be used to noninvasively measure perfusion in the brain and other organs of the body. Promising research results show how ASL might be used in stroke, tumours, dementia and paediatric medicine, in addition to many other areas. However, significant obstacles remain to prevent widespread use: ASL images have an inherently low signal to noise ratio, and are susceptible to corrupting artifacts from motion and other sources. The objective of the work in this thesis is to move towards an "intelligent imaging" paradigm: one in which the image acquisition, reconstruction and processing are mutually coupled, and tailored to the individual patient. This thesis explores how ASL images may be improved at several stages of the imaging pipeline. We review the relevant ASL literature, exploring details of ASL acquisitions, parameter inference and artifact post-processing. We subsequently present original work: we use the framework of Bayesian experimental design to generate optimised ASL acquisitions, we present original methods to improve parameter inference through anatomically-driven modelling of spatial correlation, and we describe a novel deep learning approach for simultaneous denoising and artifact filtering. Using a mixture of theoretical derivation, simulation results and imaging experiments, the work in this thesis presents several new approaches for ASL, and hopefully will shape future research and future ASL usage

Calcium channel blockade with nimodipine reverses MRI evidence of cerebral oedema following acute hypoxia

Acute cerebral hypoxia causes rapid calcium shifts leading to neuronal damage and death. Calcium channel antagonists improve outcomes in some clinical conditions, but mechanisms remain unclear. In 18 healthy participants we: (i) quantified with multiparametric MRI the effect of hypoxia on the thalamus, a region particularly sensitive to hypoxia, and on the whole brain in general; (ii) investigated how calcium channel antagonism with the drug nimodipine affects the brain response to hypoxia. Hypoxia resulted in a significant decrease in apparent diffusion coefficient (ADC), a measure particularly sensitive to cell swelling, in a widespread network of regions across the brain, and the thalamus in particular. In hypoxia, nimodipine significantly increased ADC in the same brain regions, normalizing ADC towards normoxia baseline. There was positive correlation between blood nimodipine levels and ADC change. In the thalamus, there was a significant decrease in the amplitude of low frequency fluctuations (ALFF) in resting state functional MRI and an apparent increase of grey matter volume in hypoxia, with the ALFF partially normalized towards normoxia baseline with nimodipine. This study provides further evidence that the brain response to acute hypoxia is mediated by calcium, and importantly that manipulation of intracellular calcium flux following hypoxia may reduce cerebral cytotoxic oedem