Search

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.

With colorectal cancer (CRC) being the second most common cause of cancer-related deaths worldwide1, there is an urgent need for better diagnostic tools and new, more targeted therapies. Here we used genetically engineered mouse models (GEMMs), 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 unsupervised metabolic profiling can stratify intestinal tissues according to underlying genetic alterations, and use mass spectrometry imaging (MSI) 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-mutant CRC, and propose one of its enzymes, i.e. adenosylhomocysteinase (AHCY), as a new therapeutic target. Collectively, we show that the profound genotype-dependent alterations in both lipid and small molecule metabolism in CRC may be exploited for tissue classification with no need for ion identification, and we applied further data analysis to expose a novel metabolic vulnerability of CRC.

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.

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.

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.

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.

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.

High dimensionality omics and hyperspectral imaging datasets present difficult challenges for feature extraction and data mining due to huge numbers of features that cannot be simultaneously examined. The sample numbers and variables of these methods are constantly growing as new technologies are developed, and computational analysis needs to evolve to keep up with growing demand. Current state of the art algorithms can handle some routine datasets but struggle when datasets grow above a certain size. We present a training deep learning via neural networks on non-linear dimensionality reduction, in particular t-distributed stochastic neighbour embedding (t-SNE), to overcome prior limitations of these methods.One Sentence SummaryAnalysis of prohibitively large datasets by combining deep learning via neural networks with non-linear dimensionality reduction.

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.

There are considerable gaps in our understanding of the relationship between human brain activity measured at different temporal and spatial scales by intracranial electroencephalography and fMRI. By comparing individual features and summary descriptions of intracranial EEG activity we determined which best predict fMRI changes in the sensorimotor cortex in two brain states: at rest and during motor performance. We also then examine the specificity of this relationship to spatial colocalisation of the two signals.We acquired electrocorticography and fMRI simultaneously (ECoG-fMRI) in the sensorimotor cortex of 3 patients with epilepsy. During motor activity, high gamma power was the only frequency band where the electrophysiological response was colocalised with fMRI measures across all subjects. The best model of fMRI changes was its principal components, a parsimonious description of the entire ECoG spectrogram. This model performed much better than a model based on the classical frequency bands both during task and rest periods or models derived on a summary of cross spectral changes (e.g. ‘root mean squared EEG frequency’). This suggests that the region specific fMRI signal is reflected in spatially and spectrally distributed EEG activity.

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.

Highlights • We created a validation method for the evaluation of automated classification of interictal spikes. • We used a modified version of Wave_clus (WC) to automatically classify the data of 5 patients. • WC classification was similar to EEG reviewers providing an unbiased evaluation of the clinical data., 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.

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.

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.

A study of the epileptic seizure dynamics using simultaneous recording of electroencephalography correlated functional magnetic resonance imaging (EEG-fMRI) data from 5 focal epilepsy patients undergoing pre-surgical evaluation was realised with the aim of identifying the focus and propagation network seizure-involved. A method based on the General Linear Model (GLM) at different neurophysiology regressor lags (LasgM), a connectivity model-based method, Dynamic Causal Modelling (DCM), and a data-driven method, Granger Causality (GC) were investigated. DCM analysis provided meaningful and significant results when a sufficient number of seizure events was recorded, but suffered from the generally poor data signal-to-noise ratio (SNR). The concordance between the LagsM results and the clinical expectation suggests that LagsM can be useful as a complementary approach. Intending to establish the validity of a GC approach for the problem addressed, a simulation study was performed and the results showed that it is not appropriated.