Measurement of the Mapping between Intracranial EEG and fMRI Recordings in the Human Brain.

There are considerable gaps in our understanding of the relationship between human brain activity measured at different temporal and spatial scales. Here, electrocorticography (ECoG) measures were used to predict functional MRI changes in the sensorimotor cortex in two brain states: at rest and during motor performance. The specificity of this relationship to spatial co-localisation of the two signals was also investigated. We acquired simultaneous ECoG-fMRI in the sensorimotor cortex of three patients with epilepsy. During motor activity, high gamma power was the only frequency band where the electrophysiological response was co-localised with fMRI measures across all subjects. The best model of fMRI changes across states was its principal components, a parsimonious description of the entire ECoG spectrogram. This model performed much better than any others that were based either on the classical frequency bands or on summary measures of cross-spectral changes. The region-specific fMRI signal is reflected in spatially and spectrally distributed EEG activity.

Visualisation of drug distribution in skin using correlative optical spectroscopy and mass spectrometry imaging.

A correlative methodology for label-free chemical imaging of soft tissue has been developed, combining non-linear optical spectroscopies and mass spectrometry to achieve sub-micron spatial resolution and critically improved drug detection sensitivity. The approach was applied to visualise the kinetics of drug reservoir formation within human skin following in vitro topical treatment with a commercial diclofenac gel. Non-destructive optical spectroscopic techniques, namely stimulated Raman scattering, second harmonic generation and two photon fluorescence microscopies, were used to provide chemical and structural contrast. The same tissue sections were subsequently analysed by secondary ion mass spectrometry, which offered higher sensitivity for diclofenac detection throughout the epidermis and dermis. A method was developed to combine the optical and mass spectrometric datasets using image registration techniques. The label-free, high-resolution visualisation of tissue structure coupled with sensitive chemical detection offers a powerful method for drug biodistribution studies in the skin that impact directly on topical pharmaceutical product development.

Metabolic profiling stratifies colorectal cancer and reveals adenosylhomocysteinase as a therapeutic target.

The genomic landscape of colorectal cancer (CRC) is shaped by inactivating mutations in tumour suppressors such as APC, and oncogenic mutations such as mutant KRAS. Here we used genetically engineered mouse models, and multimodal mass spectrometry-based metabolomics to study the impact of common genetic drivers of CRC on the metabolic landscape of the intestine. We show that untargeted metabolic profiling can be applied to stratify intestinal tissues according to underlying genetic alterations, and use mass spectrometry imaging to identify tumour, stromal and normal adjacent tissues. By identifying ions that drive variation between normal and transformed tissues, we found dysregulation of the methionine cycle to be a hallmark of APC-deficient CRC. Loss of Apc in the mouse intestine was found to be sufficient to drive expression of one of its enzymes, adenosylhomocysteinase (AHCY), which was also found to be transcriptionally upregulated in human CRC. Targeting of AHCY function impaired growth of APC-deficient organoids in vitro, and prevented the characteristic hyperproliferative/crypt progenitor phenotype driven by acute deletion of Apc in vivo, even in the context of mutant Kras. Finally, pharmacological inhibition of AHCY reduced intestinal tumour burden in ApcMin/+ mice indicating its potential as a metabolic drug target in CRC.

A Multimodal Desorption Electrospray Ionisation Workflow Enabling Visualisation of Lipids and Biologically Relevant Elements in a Single Tissue Section.

The colocation of elemental species with host biomolecules such as lipids and metabolites may shed new light on the dysregulation of metabolic pathways and how these affect disease pathogeneses. Alkali metals have been the subject of extensive research, are implicated in various neurodegenerative and infectious diseases and are known to disrupt lipid metabolism. Desorption electrospray ionisation (DESI) is a widely used approach for molecular imaging, but previous work has shown that DESI delocalises ions such as potassium (K) and chlorine (Cl), precluding the subsequent elemental analysis of the same section of tissue. The solvent typically used for the DESI electrospray is a combination of methanol and water. Here we show that a novel solvent system, (50:50 (%v/v) MeOH:EtOH) does not delocalise elemental species and thus enables elemental mapping to be performed on the same tissue section post-DESI. Benchmarking the MeOH:EtOH electrospray solvent against the widely used MeOH:H2O electrospray solvent revealed that the MeOH:EtOH solvent yielded increased signal-to-noise ratios for selected lipids. The developed multimodal imaging workflow was applied to a lung tissue section containing a tuberculosis granuloma, showcasing its applicability to elementally rich samples displaying defined structural information.

Exploring New Methods to Study and Moderate Proton Beam Damage for Multimodal Imaging on a Single Tissue Section.

Characterizing proton beam damage in biological materials is of interest to enable the integration of proton microprobe elemental mapping techniques with other imaging modalities. It is also of relevance to obtain a deeper understanding of mechanical damage to lipids in tissues during proton beam cancer therapy. We have developed a novel strategy to characterize proton beam damage to lipids in biological tissues based on mass spectrometry imaging. This methodology is applied to characterize changes to lipids in tissues ex vivo, irradiated under different conditions designed to mitigate beam damage. This work shows that performing proton beam irradiation at ambient pressure, as well as including the application of an organic matrix prior to irradiation, can reduce damage to lipids in tissues. We also discovered that, irrespective of proton beam irradiation, placing a sample in a vacuum prior to desorption electrospray ionization imaging can enhance lipid signals, a conclusion that may be of future benefit to the mass spectrometry imaging community.

The amino acid transporter SLC7A5 is required for efficient growth of KRAS-mutant colorectal cancer.

Oncogenic KRAS mutations and inactivation of the APC tumor suppressor co-occur in colorectal cancer (CRC). Despite efforts to target mutant KRAS directly, most therapeutic approaches focus on downstream pathways, albeit with limited efficacy. Moreover, mutant KRAS alters the basal metabolism of cancer cells, increasing glutamine utilization to support proliferation. We show that concomitant mutation of Apc and Kras in the mouse intestinal epithelium profoundly rewires metabolism, increasing glutamine consumption. Furthermore, SLC7A5, a glutamine antiporter, is critical for colorectal tumorigenesis in models of both early- and late-stage metastatic disease. Mechanistically, SLC7A5 maintains intracellular amino acid levels following KRAS activation through transcriptional and metabolic reprogramming. This supports the increased demand for bulk protein synthesis that underpins the enhanced proliferation of KRAS-mutant cells. Moreover, targeting protein synthesis, via inhibition of the mTORC1 regulator, together with Slc7a5 deletion abrogates the growth of established Kras-mutant tumors. Together, these data suggest SLC7A5 as an attractive target for therapy-resistant KRAS-mutant CRC.

Implications of Peak Selection in the Interpretation of Unsupervised Mass Spectrometry Imaging Data Analyses.

Mass spectrometry imaging can produce large amounts of complex spectral and spatial data. Such data sets are often analyzed with unsupervised machine learning approaches, which aim at reducing their complexity and facilitating their interpretation. However, choices made during data processing can impact the overall interpretation of these analyses. This work investigates the impact of the choices made at the peak selection step, which often occurs early in the data processing pipeline. The discussion is done in terms of visualization and interpretation of the results of two commonly used unsupervised approaches: t-distributed stochastic neighbor embedding and k-means clustering, which differ in nature and complexity. Criteria considered for peak selection include those based on hypotheses (exemplified herein in the analysis of metabolic alterations in genetically engineered mouse models of human colorectal cancer), particular molecular classes, and ion intensity. The results suggest that the choices made at the peak selection step have a significant impact in the visual interpretation of the results of either dimensionality reduction or clustering techniques and consequently in any downstream analysis that relies on these. Of particular significance, the results of this work show that while using the most abundant ions can result in interesting structure-related segmentation patterns that correlate well with histological features, using a smaller number of ions specifically selected based on prior knowledge about the biochemistry of the tissues under investigation can result in an easier-to-interpret, potentially more valuable, hypothesis-confirming result. Findings presented will help researchers understand and better utilize unsupervised machine learning approaches to mine high-dimensionality data.

Direct Tissue Mass Spectrometry Imaging by Atmospheric Pressure UV-Laser Desorption Plasma Postionization.

Matrix-assisted laser desorption ionization (MALDI) operated at atmospheric pressure has been shown to be a promising technique for mass spectrometry imaging of biological tissues at high spatial resolution. Recent studies have shown several orders of magnitude improvement in sensitivity afforded by coupling with a low-temperature plasma (LTP) for postionization. In this work we report the first results from "matrix-free" imaging using our atmospheric pressure (AP) transmission mode (TM) (MA)LDI source with LTP postionization. Direct MSI analysis of murine testis with no sample preparation after tissue sectioning enabled imaging of a range of lipid classes at pixel sizes of 25 µm. We compared results from the matrix-free methods with MALDI experiments in which the matrix was applied on top, underneath, or layered as a sandwich. The sandwich preparation was found to lead to ion yields approximately 2- or 3-fold higher than the other methods, indicating that the addition of a light absorbing matrix remains beneficial. Nonetheless, LDI methods confer a range of advantages, and the sensitivity improvements provided by postionization strategies are a promising step toward high-efficiency laser sampling under ambient conditions.

MALDI-2 at Atmospheric Pressure-Parameter Optimization and First Imaging Experiments.

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) is a powerful label-free technique for mapping the spatial distribution of biomolecules directly from tissue. However, like most other MSI techniques, it suffers from low ionization yields and ion suppression effects for biomolecules that might be of interest for a specific application at hand. Recently, a form of laser postionization was introduced (coined MALDI-2) that critically boosts the ion yield for many glyco- and phospholipids by several orders of magnitude and makes the detection of further biomolecular species possible. While the MALDI-2 technique is being increasingly applied by the MSI community, it is still only implemented in fine vacuum ion sources in a pressure range of about 1-10 mbar. Here, we show the first implementation of the technique to a custom-built atmospheric pressure ion source coupled to an Orbitrap Elite system. We present results from parameter optimization of MALDI-2 at atmospheric pressure, compare our findings to previously published fine vacuum data, and show first imaging results from mouse cerebellum with a 20 µm pixel size. Our findings broaden the feasibility of the technique to overall more flexible atmospheric pressure ion sources.

Construction and testing of an atmospheric-pressure transmission-mode matrix assisted laser desorption ionisation mass spectrometry imaging ion source with plasma ionisation enhancement.

Matrix assisted laser desorption ionisation mass spectrometry (MALDI-MS) at atmospheric pressure (AP) is, with a few notable exceptions, overshadowed by its vacuum based forms and AP transmission mode (TM) MALDI-MS lacks the up-take its potential benefits might suggest. The reasons for this are not fully understood and it is clear further development is required to realise the flexibility and power of this ionisation method and geometry. Here we report the build of a new AP-TM-MALDI-MSI ion source with plasma ionisation enhancement. This novel ion source is used to analyse a selection of increasingly complex systems from molecular standards to murine brain tissue sections. Significant enhancement of detected ion intensity is observed in both positive and negative ion mode in all systems, with up to 2000 fold increases observed for a range of tissue endogenous species. The substantial improvements conferred by the plasma enhancement are then employed to demonstrate the acquisition of proof of concept tissue images, with high quality spectra obtained down to 10 × 10 µm pixel size.

BOLD mapping of human epileptic spikes recorded during simultaneous intracranial EEG-fMRI: The impact of automated spike classification.

OBJECTIVES: Simultaneous intracranial EEG and functional MRI (icEEG-fMRI) can be used to map the haemodynamic (BOLD) changes associated with the generation of IEDs. Unlike scalp EEG-fMRI, in most patients who undergo icEEG-fMRI, IEDs recorded intracranially are numerous and show variability in terms of field amplitude and morphology. Therefore, visual marking can be highly subjective and time consuming. In this study, we applied an automated spike classification algorithm, Wave_clus (WC), to IEDs marked visually on icEEG data acquired during simultaneous fMRI acquisition. The motivation of this work is to determine whether using a potentially more consistent and unbiased automated approach can produce more biologically meaningful BOLD patterns compared to the BOLD patterns obtained based on the conventional, visual classification. METHODS: We analysed simultaneous icEEG-fMRI data from eight patients with severe drug resistant epilepsy, and who subsequently underwent resective surgery that resulted in a good outcome: confirmed epileptogenic zone (EZ). For each patient two fMRI analyses were performed: one based on the conventional visual IED classification and the other based on the automated classification. We used the concordance of the IED-related BOLD maps with the confirmed EZ as an indication of their biological meaning, which we compared for the automated and visual classifications for all IED originating in the EZ. RESULTS: Across the group, the visual and automated classifications resulted in 32 and 24 EZ IED classes respectively, for which 75% vs 83% of the corresponding BOLD maps were concordant. At the single-subject level, the BOLD maps for the automated approach had greater concordance in four patients, and less concordance in one patient, compared to those obtained using the conventional visual classification, and equal concordance for three remaining patients. These differences did not reach statistical significance. CONCLUSION: We found automated IED classification on icEEG data recorded during fMRI to be feasible and to result in IED-related BOLD maps that may contain similar or greater biological meaning compared to the conventional approach in the majority of the cases studied. We anticipate that this approach will help to gain significant new insights into the brain networks associated with IEDs and in relation to postsurgical outcome.

A novel scheme for the validation of an automated classification method for epileptic spikes by comparison with multiple observers.

OBJECTIVE: To validate the application of an automated neuronal spike classification algorithm, Wave_clus (WC), on interictal epileptiform discharges (IED) obtained from human intracranial EEG (icEEG) data. METHOD: Five 10-min segments of icEEG recorded in 5 patients were used. WC and three expert EEG reviewers independently classified one hundred IED events into IED classes or non-IEDs. First, we determined whether WC-human agreement variability falls within inter-reviewer agreement variability by calculating the variation of information for each classifier pair and quantifying the overlap between all WC-reviewer and all reviewer-reviewer pairs. Second, we compared WC and EEG reviewers' spike identification and individual spike class labels visually and quantitatively. RESULTS: The overlap between all WC-human pairs and all human pairs was >80% for 3/5 patients and >58% for the other 2 patients demonstrating WC falling within inter-human variation. The average sensitivity of spike marking for WC was 91% and >87% for all three EEG reviewers. Finally, there was a strong visual and quantitative similarity between WC and EEG reviewers. CONCLUSIONS: WC performance is indistinguishable to that of EEG reviewers' suggesting it could be a valid clinical tool for the assessment of IEDs. SIGNIFICANCE: WC can be used to provide quantitative analysis of epileptic spikes.

Simultaneous Intracranial EEG-fMRI Shows Inter-Modality Correlation in Time-Resolved Connectivity Within Normal Areas but Not Within Epileptic Regions.

For the first time in research in humans, we used simultaneous icEEG-fMRI to examine the link between connectivity in haemodynamic signals during the resting-state (rs) and connectivity derived from electrophysiological activity in terms of the inter-modal connectivity correlation (IMCC). We quantified IMCC in nine patients with drug-resistant epilepsy (i) within brain networks in 'healthy' non-involved cortical zones (NIZ) and (ii) within brain networks involved in generating seizures and interictal spikes (IZ1) or solely spikes (IZ2). Functional connectivity (h 2 ) estimates for 10 min of resting-state data were obtained between each pair of electrodes within each clinical zone for both icEEG and fMRI. A sliding window approach allowed us to quantify the variability over time of h 2 (vh 2) as an indicator of connectivity dynamics. We observe significant positive IMCC for h 2 and vh 2, for multiple bands in the NIZ only, with the strongest effect in the lower icEEG frequencies. Similarly, intra-modal h 2 and vh 2 were found to be differently modified as a function of different epileptic processes: compared to NIZ, [Formula: see text] was higher in IZ1, but lower in IZ2, while [Formula: see text] showed the inverse pattern. This corroborates previous observations of inter-modal connectivity discrepancies in pathological cortices, while providing the first direct invasive and simultaneous comparison in humans. We also studied time-resolved FC variability multimodally for the first time, finding that IZ1 shows both elevated internal [Formula: see text] and less rich dynamical variability, suggesting that its chronic role in epileptogenesis may be linked to greater homogeneity in self-sustaining pathological oscillatory states.

Electrophysiological correlates of the BOLD signal for EEG-informed fMRI.

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are important tools in cognitive and clinical neuroscience. Combined EEG-fMRI has been shown to help to characterise brain networks involved in epileptic activity, as well as in different sensory, motor and cognitive functions. A good understanding of the electrophysiological correlates of the blood oxygen level-dependent (BOLD) signal is necessary to interpret fMRI maps, particularly when obtained in combination with EEG. We review the current understanding of electrophysiological-haemodynamic correlates, during different types of brain activity. We start by describing the basic mechanisms underlying EEG and BOLD signals and proceed by reviewing EEG-informed fMRI studies using fMRI to map specific EEG phenomena over the entire brain (EEG-fMRI mapping), or exploring a range of EEG-derived quantities to determine which best explain colocalised BOLD fluctuations (local EEG-fMRI coupling). While reviewing studies of different forms of brain activity (epileptic and nonepileptic spontaneous activity; cognitive, sensory and motor functions), a significant attention is given to epilepsy because the investigation of its haemodynamic correlates is the most common application of EEG-informed fMRI. Our review is focused on EEG-informed fMRI, an asymmetric approach of data integration. We give special attention to the invasiveness of electrophysiological measurements and the simultaneity of multimodal acquisitions because these methodological aspects determine the nature of the conclusions that can be drawn from EEG-informed fMRI studies. We emphasise the advantages of, and need for, simultaneous intracranial EEG-fMRI studies in humans, which recently became available and hold great potential to improve our understanding of the electrophysiological correlates of BOLD fluctuations.

Characterisation and reduction of the EEG artefact caused by the helium cooling pump in the MR environment: validation in epilepsy patient data.

The EEG acquired simultaneously with fMRI is distorted by a number of artefacts related to the presence of strong magnetic fields, which must be reduced in order to allow for a useful interpretation and quantification of the EEG data. For the two most prominent artefacts, associated with magnetic field gradient switching and the heart beat, reduction methods have been developed and applied successfully. However, a number of artefacts related to the MR-environment can be found to distort the EEG data acquired even without ongoing fMRI acquisition. In this paper, we investigate the most prominent of those artefacts, caused by the Helium cooling pump, and propose a method for its reduction and respective validation in data collected from epilepsy patients. Since the Helium cooling pump artefact was found to be repetitive, an average template subtraction method was developed for its reduction with appropriate adjustments for minimizing the degradation of the physiological part of the signal. The new methodology was validated in a group of 15 EEG-fMRI datasets collected from six consecutive epilepsy patients, where it successfully reduced the amplitude of the artefact spectral peaks by 95 ± 2 % while the background spectral amplitude within those peaks was reduced by only -5 ± 4 %. Although the Helium cooling pump should ideally be switched off during simultaneous EEG-fMRI acquisitions, we have shown here that in cases where this is not possible the associated artefact can be effectively reduced in post processing.

Dynamic Causal Modelling of epileptic seizure propagation pathways: a combined EEG-fMRI study.

Simultaneous EEG-fMRI offers the possibility of non-invasively studying the spatiotemporal dynamics of epileptic activity propagation from the focus towards an extended brain network, through the identification of the haemodynamic correlates of ictal electrical discharges. In epilepsy associated with hypothalamic hamartomas (HH), seizures are known to originate in the HH but different propagation pathways have been proposed. Here, Dynamic Causal Modelling (DCM) was employed to estimate the seizure propagation pathway from fMRI data recorded in a HH patient, by testing a set of clinically plausible network connectivity models of discharge propagation. The model consistent with early propagation from the HH to the temporal-occipital lobe followed by the frontal lobe was selected as the most likely model to explain the data. Our results demonstrate the applicability of DCM to investigate patient-specific effective connectivity in epileptic networks identified with EEG-fMRI. In this way, it is possible to study the propagation pathway of seizure activity, which has potentially great impact in the decision of the surgical approach for epilepsy treatment.