Aneurin J. Kennerley

Negative Blood Oxygen Level Dependence in the Rat:A Model for Investigating the Role of Suppression in Neurovascular Coupling

Modern neuroimaging techniques rely on neurovascular coupling to show regions of increased brain activation. However, little is known of the neurovascular coupling relationships that exist for inhibitory signals. To address this issue directly we developed a preparation to investigate the signal sources of one of these proposed inhibitory neurovascular signals, the negative blood oxygen level-dependent (BOLD) response (NBR), in rat somatosensory cortex. We found a reliable NBR measured in rat somatosensory cortex in response to unilateral electrical whisker stimulation, which was located in deeper cortical layers relative to the positive BOLD response. Separate optical measurements (two-dimensional optical imaging spectroscopy and laser Doppler flowmetry) revealed that the NBR was a result of decreased blood volume and flow and increased levels of deoxyhemoglobin. Neural activity in the NBR region, measured by multichannel electrodes, varied considerably as a function of cortical depth. There was a decrease in neuronal activity in deep cortical laminae. After cessation of whisker stimulation there was a large increase in neural activity above baseline. Both the decrease in neuronal activity and increase above baseline after stimulation cessation correlated well with the simultaneous measurement of blood flow suggesting that the NBR is related to decreases in neural activity in deep cortical layers. Interestingly, the magnitude of the neural decrease was largest in regions showing stimulus-evoked positive BOLD responses. Since a similar type of neural suppression in surround regions was associated with a negative BOLD signal, the increased levels of suppression in positive BOLD regions could importantly moderate the size of the observed BOLD response.

Concurrent fMRI and optical measures for the investigation of the hemodynamic response function

Functional magnetic resonance imaging (fMRI) signal variations are based on a combination of changes in cerebral blood flow (CBF) and volume (CBV), and blood oxygenation. We investigated the relationship between these hemodynamic parameters in the rodent barrel cortex by performing fMRI concurrently with laser Doppler flowmetry (LDF) or optical imaging spectroscopy (OIS), following whisker stimulation and hypercapnic challenge. A difference between the positions of the maximum blood oxygenation level‐dependent (BOLD) and CBV changes was observed in coronal fMRI maps, with the BOLD region being more superficial. A 6.5% baseline blood volume fraction in this superficial region dropped to 4% in deeper cortical layers (corresponding to total hemoglobin baseline volumes Hbt0 = 110 μM and 67 μM, respectively), as inferred from maps of ΔR  2* . Baseline volume profiles were used to parameterize the Monte Carlo simulations (MCS) to interpret the 2D OIS. From this it was found that the optical blood volume measurements (i.e., changes in total hemoglobin) equated with CBV‐MRI measurements when the MRI data were taken from superficial cortical layers. Optical measures of activation showed a good spatial overlap with fMRI measurements taken in the same plane (covering the right hemisphere surface). Changes in CBV and CBF followed the scaling relationship CBV = CBFα, with mean α = 0.38 ± 0.06. Magn Reson Med 54:354–365, 2005.

Neurovascular coupling is brain region-dependent

Despite recent advances in alternative brain imaging technologies, functional magnetic resonance imaging (fMRI) remains the workhorse for both medical diagnosis and primary research. Indeed, the number of research articles that utilise fMRI have continued to rise unabated since its conception in 1991, despite the limitation that recorded signals originate from the cerebral vasculature rather than neural tissue. Consequently, understanding the relationship between brain activity and the resultant changes in metabolism and blood flow (neurovascular coupling) remains a vital area of research. In the past, technical constraints have restricted investigations of neurovascular coupling to cortical sites and have led to the assumption that coupling in non-cortical structures is the same as in the cortex, despite the lack of any evidence. The current study investigated neurovascular coupling in the rat using whole-brain blood oxygenation level-dependent (BOLD) fMRI and multi-channel electrophysiological recordings and measured the response to a sensory stimulus as it proceeded through brainstem, thalamic and cortical processing sites - the so-called whisker-to-barrel pathway. We found marked regional differences in the amplitude of BOLD activation in the pathway and non-linear neurovascular coupling relationships in non-cortical sites. The findings have important implications for studies that use functional brain imaging to investigate sub-cortical function and caution against the use of simple, linear mapping of imaging signals onto neural activity.

Cortical lamina-dependent blood volume changes in human brain at 7 T

Cortical layer-dependent high (sub-millimeter) resolution functional magnetic resonance imaging (fMRI) in human or animal brain can be used to address questions regarding the functioning of cortical circuits, such as the effect of different afferent and efferent connectivities on activity in specific cortical layers. The sensitivity of gradient echo (GE) blood oxygenation level-dependent (BOLD) responses to large draining veins reduces its local specificity and can render the interpretation of the underlying laminar neural activity impossible. The application of the more spatially specific cerebral blood volume (CBV)-based fMRI in humans has been hindered by the low sensitivity of the noninvasive modalities available. Here, a vascular space occupancy (VASO) variant, adapted for use at high field, is further optimized to capture layer-dependent activity changes in human motor cortex at sub-millimeter resolution. Acquired activation maps and cortical profiles show that the VASO signal peaks in gray matter at 0.8-1.6mm depth, and deeper compared to the superficial and vein-dominated GE-BOLD responses. Validation of the VASO signal change versus well-established iron-oxide contrast agent based fMRI methods in animals showed the same cortical profiles of CBV change, after normalization for lamina-dependent baseline CBV. In order to evaluate its potential of revealing small lamina-dependent signal differences due to modulations of the input-output characteristics, layer-dependent VASO responses were investigated in the ipsilateral hemisphere during unilateral finger tapping. Positive activation in ipsilateral primary motor cortex and negative activation in ipsilateral primary sensory cortex were observed. This feature is only visible in high-resolution fMRI where opposing sides of a sulcus can be investigated independently because of a lack of partial volume effects. Based on the results presented here, we conclude that VASO offers good reproducibility, high sensitivity and lower sensitivity than GE-BOLD to changes in larger vessels, making it a valuable tool for layer-dependent fMRI studies in humans.

Optimised and rapid pre-clinical screening in the SOD1(G93A) transgenic mouse model of amyotrophic lateral sclerosis (ALS).

The human SOD1G93A transgenic mouse has been used extensively since its development in 1994 as a model for amyotrophic lateral sclerosis (ALS). In that time, a great many insights into the toxicity of mutant SOD1 have been gained using this and other mutant SOD transgenic mouse models. They all demonstrate a selective toxicity towards motor neurons and in some cases features of the pathology seen in the human disease. These models have two major drawbacks. Firstly the generation of robust preclinical data in these models has been highlighted as an area for concern. Secondly, the amount of time required for a single preclinical experiment in these models (3–4 months) is a hurdle to the development of new therapies. We have developed an inbred C57BL/6 mouse line from the original mixed background (SJLxC57BL/6) SOD1G93A transgenic line and show here that the disease course is remarkably consistent and much less prone to background noise, enabling reduced numbers of mice for testing of therapeutics. Secondly we have identified very early readouts showing a large decline in motor function compared to normal mice. This loss of motor function has allowed us to develop an early, sensitive and rapid screening protocol for the initial phases of denervation of muscle fibers, observed in this model. We describe multiple, quantitative readouts of motor function that can be used to interrogate this early mechanism. Such an approach will increase throughput for reduced costs, whilst reducing the severity of the experimental procedures involved.

Theory and generalization of Monte Carlo models of the BOLD signal source.

The dependency of the blood oxygenation level dependent (BOLD) signal on underlying hemodynamics is not well understood. Building a forward biophysical model of this relationship is important for the quantitative estimation of the hemodynamic changes and neural activity underlying functional magnetic resonance imaging (fMRI) signals. We have developed a general model of the BOLD signal which can model both intra‐ and extravascular signals for an arbitrary tissue model across a wide range of imaging parameters. The model of the BOLD signal was instantiated as a look‐up‐table (LuT), and was verified against concurrent fMRI and optical imaging measurements of activation induced hemodynamics. Magn Reson Med, 2008.

Kinetic modeling of hyperpolarized 13C pyruvate metabolism in tumors using a measured arterial input function

Mathematical models are required to estimate kinetic parameters of [1‐13C] pyruvate‐lactate interconversion from magnetic resonance spectroscopy data. One‐ or two‐way exchange models utilizing a hypothetical approximation to the true arterial input function (AIF), (e.g. an ideal ‘box‐car’ function) have been used previously. We present a method for direct measurement of the AIF in the rat. The hyperpolarized [1‐13C] pyruvate signal was measured in arterial blood as it was continuously withdrawn through a small chamber. The measured signal was corrected for T1 relaxation of pyruvate, RF pulses and dispersion of blood in the chamber to allow for the estimation of the direct AIF. Using direct AIF, rather than the commonly used box‐car AIF, provided realistic estimates of the rate constant of conversion of pyruvate to lactate, kpl, the rate constant of conversion of lactate to pyruvate klp, the clearance rate constant of pyruvate from blood to tissue, Kip, and the relaxation rate of lactate T1la. Since no lactate signal was present in blood, it was possible to use a simple precursor‐product relationship, with the tumor tissue pyruvate time‐course as the input for the lactate time‐course. This provided a robust estimate of kpl, similar to that obtained using a directly measured AIF. Magn Reson Med, 70:943–953, 2013.

Directing cell therapy to anatomic target sites in vivo with magnetic resonance targeting

Cell-based therapy exploits modified human cells to treat diseases but its targeted application in specific tissues, particularly those lying deep in the body where direct injection is not possible, has been problematic. Here we use a magnetic resonance imaging (MRI) system to direct macrophages carrying an oncolytic virus, Seprehvir, into primary and metastatic tumour sites in mice. To achieve this, we magnetically label macrophages with super-paramagnetic iron oxide nanoparticles and apply pulsed magnetic field gradients in the direction of the tumour sites. Magnetic resonance targeting guides macrophages from the bloodstream into tumours, resulting in increased tumour macrophage infiltration and reduction in tumour burden and metastasis. Our study indicates that clinical MRI scanners can not only track the location of magnetically labelled cells but also have the potential to steer them into one or more target tissues.

Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling.

Understanding neurovascular coupling is a prerequisite for the interpretation of results obtained from modern neuroimaging techniques. This study investigated the hemodynamic and neural responses in rat somatosensory cortex elicited by 16 seconds electrical whisker stimuli. Hemodynamics were measured by optical imaging spectroscopy and neural activity by multichannel electrophysiology. Previous studies have suggested that the whisker-evoked hemodynamic response contains two mechanisms, a transient ‘backwards’ dilation of the middle cerebral artery, followed by an increase in blood volume localized to the site of neural activity. To distinguish between the mechanisms responsible for these aspects of the response, we presented whisker stimuli during normocapnia (‘control’), and during a high level of hypercapnia. Hypercapnia was used to ‘predilate’ arteries and thus possibly ‘inhibit’ aspects of the response related to the ‘early’ mechanism. Indeed, hemodynamic data suggested that the transient stimulus-evoked response was absent under hypercapnia. However, evoked neural responses were also altered during hypercapnia and convolution of the neural responses from both the normocapnic and hypercapnic conditions with a canonical impulse response function, suggested that neurovascular coupling was similar in both conditions. Although data did not clearly dissociate early and late vascular responses, they suggest that the neurovascular coupling relationship is neurogenic in origin.

A dynamic model of neurovascular coupling: Implications for blood vessel dilation and constriction

Neurovascular coupling in response to stimulation of the rat barrel cortex was investigated using concurrent multichannel electrophysiology and laser Doppler flowmetry. The data were used to build a linear dynamic model relating neural activity to blood flow. Local field potential time series were subject to current source density analysis, and the time series of a layer IV sink of the barrel cortex was used as the input to the model. The model output was the time series of the changes in regional cerebral blood flow (CBF). We show that this model can provide excellent fit of the CBF responses for stimulus durations of up to 16 s. The structure of the model consisted of two coupled components representing vascular dilation and constriction. The complex temporal characteristics of the CBF time series were reproduced by the relatively simple balance of these two components. We show that the impulse response obtained under the 16-s duration stimulation condition generalised to provide a good prediction to the data from the shorter duration stimulation conditions. Furthermore, by optimising three out of the total of nine model parameters, the variability in the data can be well accounted for over a wide range of stimulus conditions. By establishing linearity, classic system analysis methods can be used to generate and explore a range of equivalent model structures (e.g., feed-forward or feedback) to guide the experimental investigation of the control of vascular dilation and constriction following stimulation.