Study-level wavelet cluster analysis and data-driven signal models in pharmacological MRI

Modeling the acute pharmacological response to selective serotonin reuptake inhibitors in human brain using simultaneous PET/MR imaging

Pharmacological imaging of the effects of selective serotonin reuptake inhibitors (SSRI) may aid the clarification of their mechanism of action and influence treatment of highly prevalent neuropsychiatric conditions if the detected effects could be related to patient outcomes. In a randomized double-blind design, 38 healthy participants received a constant infusion of 8 mg citalopram or saline during either their first or second of two PET/MR scans. Resting-state functional MRI (fMRI) was acquired simultaneously with PET data on the binding of serotonin transporters (5-HTT) using [¹¹C]DASB. Three different approaches for modeling of pharmacological fMRI response were tested separately. These relied on the use of regressors corresponding to (1) the drug infusion paradigm, (2) time courses of citalopram plasma concentrations and (3) changes in 5-HTT binding measured in each individual, respectively. Furthermore, the replication of results of a widely used model-free analysis method was attempted which assesses the deviation of signal in discrete time bins of fMRI data acquired after start of drug infusion. Following drug challenge, average 5-HTT occupancy was 69±7% and peak citalopram plasma levels were 111.8 ± 21.1 ng/ml. None of the applied methods could detect significant differences in the pharmacological response between SSRI and placebo scans. The failed replication of SSRI effects reported in the literature despite a threefold larger sample size highlights the importance of appropriate correction for family-wise error in order to avoid spurious results in pharmacological imaging. This calls for the development of analysis methods which take regional specialization and the dynamics of brain activity into account.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

Data collection and analysis strategies for phMRI

Although functional MRI traditionally has been applied mainly to study changes in task-induced brain function, evolving acquisition methodologies and improved knowledge of signal mechanisms have increased the utility of this method for studying responses to pharmacological stimuli, a technique often dubbed "phMRI". The proliferation of higher magnetic field strengths and the use of exogenous contrast agent have boosted detection power, a critical factor for successful phMRI due to the restricted ability to average multiple stimuli within subjects. Receptor-based models of neurovascular coupling, including explicit pharmacological models incorporating receptor densities and affinities and data-driven models that incorporate weak biophysical constraints, have demonstrated compelling descriptions of phMRI signal induced by dopaminergic stimuli. This report describes phMRI acquisition and analysis methodologies, with an emphasis on data-driven analyses. As an example application, statistically efficient data-driven regressors were used to describe the biphasic response to the mu-opioid agonist remifentanil, and antagonism using dopaminergic and GABAergic ligands revealed modulation of the mesolimbic pathway. Results illustrate the power of phMRI as well as our incomplete understanding of mechanisms underlying the signal. Future directions are discussed for phMRI acquisitions in human studies, for evolving analysis methodologies, and for interpretative studies using the new generation of simultaneous PET/MRI scanners. This article is part of the Special Issue Section entitled 'Neuroimaging in Neuropharmacology'.

Test-retest reliability of the BOLD pharmacological MRI response to ketamine in healthy volunteers

The pharmacological MRI (phMRI) technique is being increasingly used in both pre-clinical and clinical models to investigate pharmacological effects on task-free brain function. Ketamine, an N-methyl-d-aspartate receptor (NMDAR) antagonist, induces a strong phMRI response and represents a promising pharmacological model to investigate the role of glutamatergic abnormalities in psychiatric symptomatology. The aim of this study was to assess whether the brain response to ketamine is reliable in order to validate ketamine phMRI as a mechanistic marker of glutamatergic dysfunction and to determine its utility in repeated measures designs to detect the modulatory effect of other drugs. Thus we assessed the test-retest reliability of the brain response to ketamine in healthy volunteers and identified an optimal modelling approach with reliability as our selection criterion. PhMRI data were collected from 10 healthy male participants, at rest, on two separate occasions. Subanaesthetic doses of I.V. ketamine infusion (target plasma levels 50 ng/mL and 75 ng/mL) were administered in both sessions. Test-retest reliability of the ketamine phMRI response was assessed voxel-wise and on pre-defined ROIs for a range of temporal design matrices including different combinations of nuisance regressors designed to model shape variance, linear drift and head motion. Effect sizes are also reported. All models showed a significant and widespread response to low-dose ketamine in predicted cerebral networks and as expected, increasing the number of model parameters improved model fit. Reliability of the predefined ROIs differed between the different models assessed. Using reliability as the selection criterion, a model capturing subject motion and linear drift performed the best across two sessions. The anatomical distribution of effects for all models was consistent with results of previous imaging studies in humans with BOLD signal increases in regions including midline cingulate and supracingulate cortex, thalamus, insula, anterior temporal lobe and ventrolateral prefrontal structures, and BOLD signal decreases in the subgenual cingulate cortex. This study represents the first investigation of the test-retest reliability of the BOLD phMRI response to acute ketamine challenge. All models tested were effective at describing the ketamine response although the design matrix associated with the highest reliability may represent a robust and well-characterised ketamine phMRI assay more suitable for repeated-measures designs. This ketamine assay is applicable as a model of neurotransmitter dysfunction suitable as a pharmacodynamic imaging tool to test and validate modulatory interventions, as a model of NMDA hypofunction in psychiatric disorders, and may be adapted to understand potential antidepressant and analgesic effects of NMDAR antagonists.

Brain reinforcement system function is ghrelin dependent: studies in the rat using pharmacological fMRI and intracranial self-stimulation

Ghrelin (GHR) is an orexigenic gut peptide that interacts with brain ghrelin receptors (GHR-Rs) to promote food intake. Recent research suggests that GHR acts as a modulator of motivated behavior, suggesting a direct influence of GHR on brain reinforcement circuits. In the present studies, we investigated the role of GHR and GHR-Rs in brain reinforcement function. Pharmacological magnetic resonance imaging was used to spatially resolve the functional activation produced by systemic administration of an orexigenic GHR dose. The imaging data revealed a focal activation of a network of subcortical structures that comprise brain reinforcement circuits-ventral tegmental area, lateral hypothalamus and nucleus accumbens. We next analyzed whether brain reinforcement circuits require functional GHR-Rs. To this purpose, wild-type (WT) or mutant rats sustaining N-ethyl-N-nitrosourea-induced knockout of GHR-Rs (GHR-R null rats) were implanted with stimulating electrodes aimed at the lateral hypothalamus, shaped to respond for intracranial self-stimulation (ICSS) and then tested using a rate-frequency procedure to examine ICSS response patterns. WT rats were readily shaped using stimulation intensities of 75 µA, whereas GHR-R null rats required 300 µA for ICSS shaping. No differences in rate-frequency curves were noted for WT rats at 75 µA and GHR-R null rats at 300 µA. When current intensity was lowered to 100 µA, GHR-R null rats did not respond for ICSS. Taken collectively, these data suggest that systemic GHR can activate mesolimbic dopaminergic areas, and highlight a facilitative role of GHR-Rs on the activity of brain reinforcement systems.

Voxel scale complex networks of functional connectivity in the rat brain: neurochemical state dependence of global and local topological properties

Network analysis of functional imaging data reveals emergent features of the brain as a function of its topological properties. However, the brain is not a homogeneous network, and the dependence of functional connectivity parameters on neuroanatomical substrate and parcellation scale is a key issue. Moreover, the extent to which these topological properties depend on underlying neurochemical changes remains unclear. In the present study, we investigated both global statistical properties and the local, voxel-scale distribution of connectivity parameters of the rat brain. Different neurotransmitter systems were stimulated by pharmacological challenge (d-amphetamine, fluoxetine, and nicotine) to discriminate between stimulus-specific functional connectivity and more general features of the rat brain architecture. Although global connectivity parameters were similar, mapping of local connectivity parameters at high spatial resolution revealed strong neuroanatomical dependence of functional connectivity in the rat brain, with clear differentiation between the neocortex and older brain regions. Localized foci of high functional connectivity independent of drug challenge were found in the sensorimotor cortices, consistent with the high neuronal connectivity in these regions. Conversely, the topological properties and node roles in subcortical regions varied with neurochemical state and were dependent on the specific dynamics of the different functional processes elicited.

Modulation of fronto-cortical activity by modafinil: a functional imaging and fos study in the rat

Modafinil (MOD) is a wake-promoting drug with pro-cognitive properties. Despite its increasing use, the neuronal substrates of MOD action remain elusive. In particular, animal studies have highlighted a putative role of diencephalic areas as primary neuronal substrate of MOD action, with inconsistent evidence of recruitment of fronto-cortical areas despite the established pro-cognitive effects of the drug. Moreover, most animal studies have employed doses of MOD of limited clinical relevance. We used pharmacological magnetic resonance imaging (phMRI) in the anesthetized rat to map the circuitry activated by a MOD dose producing clinically relevant plasma exposure, as here ascertained by pharmacokinetic measurements. We observed prominent and sustained activation of the prefrontal and cingulate cortex, together with weaker but significant activation of the somatosensory cortex, medial thalamic domains, hippocampus, ventral striatum and dorsal raphe. Correlation analysis of phMRI data highlighted enhanced connectivity within a neural network including dopamine projections from the ventral tegmental area to the nucleus accumbens. The pro-arousing effect of MOD was assessed using electroencephalographic recording under anesthetic conditions comparable to those used for phMRI, together with the corresponding Fos immunoreactivity distribution. MOD produced electroencephalogram desynchronization, resulting in reduced delta and increased theta frequency bands, and a pattern of Fos induction largely consistent with the phMRI study. Altogether, these findings show that clinically relevant MOD doses can robustly activate fronto-cortical areas involved in higher cognitive functions and a network of pro-arousing areas, which provide a plausible substrate for the wake-promoting and pro-cognitive effects of the drug.

Neuroimaging evidence of altered fronto-cortical and striatal function after prolonged cocaine self-administration in the rat

Cocaine addiction is often modeled in experimental paradigms where rodents learn to self-administer (SA) the drug. However, the extent to which these models replicate the functional alterations observed in clinical neuroimaging studies of cocaine addiction remains unknown. We used magnetic resonance imaging (MRI) to assess basal and evoked brain function in rats subjected to a prolonged, extended-access cocaine SA scheme. Specifically, we measured basal cerebral blood volume (bCBV), an established correlate of basal metabolism, and assessed the reactivity of the dopaminergic system by mapping the pharmacological MRI (phMRI) response evoked by the dopamine-releaser amphetamine. Cocaine-exposed subjects exhibited reduced bCBV in fronto-cortical areas, nucleus accumbens, ventral hippocampus, and thalamus. The cocaine group also showed an attenuated functional response to amphetamine in ventrostriatal areas, an effect that was significantly correlated with total cocaine intake. An inverse relationship between bCBV in the reticular thalamus and the frontal response elicited by amphetamine was found in control subjects but not in the cocaine group, suggesting that the inhibitory interplay within this attentional circuit may be compromised by the drug. Importantly, histopathological analysis did not reveal significant alterations of the microvascular bed in the brain of cocaine-exposed subjects, suggesting that the imaging findings cannot be merely ascribed to cocaine-induced vascular damage. These results document that chronic, extended-access cocaine SA in the rat produces focal fronto-cortical and striatal alterations that serve as plausible neurobiological substrate for the behavioral expression of compulsive drug intake in laboratory animals.

Functional magnetic resonance imaging reveals different neural substrates for the effects of orexin-1 and orexin-2 receptor antagonists

Orexins are neuro-modulatory peptides involved in the control of diverse physiological functions through interaction with two receptors, orexin-1 (OX1R) and orexin-2 (OX2R). Recent evidence in pre-clinical models points toward a putative dichotomic role of the two receptors, with OX2R predominantly involved in the regulation of the sleep/wake cycle and arousal, and the OX1R being more specifically involved in reward processing and motivated behaviour. However, the specific neural substrates underlying these distinct processes in the rat brain remain to be elucidated. Here we used functional magnetic resonance imaging (fMRI) in the rat to map the modulatory effect of selective OXR blockade on the functional response produced by D-amphetamine, a psychostimulant and arousing drug that stimulates orexigenic activity. OXR blockade was produced by GSK1059865 and JNJ1037049, two novel OX1R and OX2R antagonists with unprecedented selectivity at the counter receptor type. Both drugs inhibited the functional response to D-amphetamine albeit with distinct neuroanatomical patterns: GSK1059865 focally modulated functional responses in striatal terminals, whereas JNJ1037049 induced a widespread pattern of attenuation characterised by a prominent cortical involvement. At the same doses tested in the fMRI study, JNJ1037049 exhibited robust hypnotic properties, while GSK1059865 failed to display significant sleep-promoting effects, but significantly reduced drug-seeking behaviour in cocaine-induced conditioned place preference. Collectively, these findings highlight an essential contribution of the OX2R in modulating cortical activity and arousal, an effect that is consistent with the robust hypnotic effect exhibited by JNJ1037049. The subcortical and striatal pattern observed with GSK1059865 represent a possible neurofunctional correlate for the modulatory role of OX1R in controlling reward-processing and goal-oriented behaviours in the rat.

Functional and pharmacological MRI in understanding brain function at a systems level

Functional magnetic resonance imaging (fMRI) methods have been extensively applied to study the human brain and its functional organization in healthy and disease states. A strong rationale exists for the extension of this approach to animal models as a translational tool to bridge clinical and preclinical research. Specifically, the development of pharmacological MRI (phMRI), i.e., the use of fMRI to map spatiotemporal patterns of brain activity induced by pharmacological agents, has provided a robust and flexible tool to resolve brain circuits and mechanism-specific functional changes produced by selective intervention in different neurotransmitter systems in vivo. This chapter describes the methodological aspects of fMRI and phMRI in preclinical species, and some of the key findings, with a special emphasis on the translational potential of these methods in neuropharmacological research.

Wavelet-based fMRI analysis: 3-D denoising, signal separation, and validation metrics

We present a novel integrated wavelet-domain based framework (w-ICA) for 3-D denoising functional magnetic resonance imaging (fMRI) data followed by source separation analysis using independent component analysis (ICA) in the wavelet domain. We propose the idea of a 3-D wavelet-based multi-directional denoising scheme where each volume in a 4-D fMRI data set is sub-sampled using the axial, sagittal and coronal geometries to obtain three different slice-by-slice representations of the same data. The filtered intensity value of an arbitrary voxel is computed as an expected value of the denoised wavelet coefficients corresponding to the three viewing geometries for each sub-band. This results in a robust set of denoised wavelet coefficients for each voxel. Given the de-correlated nature of these denoised wavelet coefficients, it is possible to obtain more accurate source estimates using ICA in the wavelet domain. The contributions of this work can be realized as two modules: First, in the analysis module we combine a new 3-D wavelet denoising approach with signal separation properties of ICA in the wavelet domain. This step helps obtain an activation component that corresponds closely to the true underlying signal, which is maximally independent with respect to other components. Second, we propose and describe two novel shape metrics for post-ICA comparisons between activation regions obtained through different frameworks. We verified our method using simulated as well as real fMRI data and compared our results against the conventional scheme (Gaussian smoothing+spatial ICA: s-ICA). The results show significant improvements based on two important features: (1) preservation of shape of the activation region (shape metrics) and (2) receiver operating characteristic curves. It was observed that the proposed framework was able to preserve the actual activation shape in a consistent manner even for very high noise levels in addition to significant reduction in false positive voxels.

A statistical framework for optimal design matrix generation with application to fMRI

The general linear model (GLM) is a well established tool for analyzing functional magnetic resonance imaging (fMRI) data. Most fMRI analyses via GLM proceed in a massively univariate fashion where the same design matrix is used for analyzing data from each voxel. A major limitation of this approach is the locally varying nature of signals of interest as well as associated confounds. This local variability results in a potentially large bias and uncontrolled increase in variance for the contrast of interest. The main contributions of this paper are two fold: 1) we develop a statistical framework that enables estimation of an optimal design matrix while explicitly controlling the bias variance decomposition over a set of potential design matrices and 2) we develop and validate a numerical algorithm for computing optimal design matrices for general fMRI data sets. The implications of this framework include the ability to match optimally the magnitude of underlying signals to their true magnitudes while also matching the "null" signals to zero size thereby optimizing both the sensitivity and specificity of signal detection. By enabling the capture of multiple profiles of interest using a single contrast (as opposed to an F-test) in a way that optimizes for both bias and variance enables the passing of first level parameter estimates and their variances to the higher level for group analysis which is not possible using F-tests. We demonstrate the application of this approach to in vivo pharmacological fMRI data capturing the acute response to a drug infusion, to task-evoked, block design fMRI and to the estimation of a haemodynamic response function (HRF) in event-related fMRI. Although developed with motivation from fMRI, our framework is quite general and has potentially wide applicability to a variety of disciplines.

Functional connectivity in the rat brain: a complex network approach

Functional connectivity analyses of fMRI data can provide a wealth of information on the brain functional organization and have been widely applied to the study of the human brain. More recently, these methods have been extended to preclinical species, thus providing a powerful translational tool. Here, we review methods and findings of functional connectivity studies in the rat. More specifically, we focus on correlation analysis of pharmacological MRI (phMRI) responses, an approach that has enabled mapping the patterns of connectivity underlying major neurotransmitter systems in vivo. We also review the use of novel statistical approaches based on a network representation of the functional connectivity and their application to the study of the rat brain functional architecture.

Antagonism at serotonin 5-HT(2A) receptors modulates functional activity of frontohippocampal circuit

RATIONALE

Several second-generation antipsychotics are characterised by a significant antagonistic effect at serotonin 5-HT(2A) receptors (5-HT(2A)R), a feature that has been associated with lower incidence of extra-pyramidal symptoms and a putative amelioration of positive and negative symptoms experienced by schizophrenic patients. However, the neurofunctional substrate of 5-HT(2A) antagonism and its exact contribution to the complex pharmacological profile of these drugs remain to be elucidated.

OBJECTIVES

Here, we used pharmacological magnetic resonance imaging to map the modulatory effects of the selective 5-HT(2A)R antagonist Ml00907 on the spatiotemporal patterns of brain activity elicited by acute phencyclidine (PCP) challenge in the rat. PCP is a non-competitive NMDA receptor antagonist that induces dysregulation of corticolimbic glutamatergic neurotransmission and produces cognitive impairment and psychotic-like symptoms reminiscent of those observed in schizophrenia.

RESULTS

Pre-administration of M100907 produced focal and region-dependent attenuation of PCP-induced response in frontoseptohippocampal areas. As early studies highlighted a permissive role of 5-HT(2A)R on frontal dopamine release, the role of post-synaptic dopamine D(1) receptors on PCP-induced response was examined by using the potent antagonist SCH23390. Interestingly, SCH23390 did not affect PCP's response in any of the regions examined. This finding rules out a significant contribution of dopamine in the functional changes mapped and, indirectly, the inhibitory effect of M100907, in favour of a glutamatergic origin.

CONCLUSIONS

Our data expand recent evidence suggesting a key role of 5-HT(2A)R in modulating glutamate-mediated cognitive performance in the prefrontal cortex and highlight the whole frontoseptohippocampal circuit as a key functional substrate of 5-HT(2A)R antagonism in normal and disease states.

Pharmacological stimulation of NMDA receptors via co-agonist site suppresses fMRI response to phencyclidine in the rat

RATIONALE

Increasing experimental evidence suggests that impaired N-methyl-D: -aspartic acid (NMDA) receptor (NMDAr) function could be a key pathophysiological determinant of schizophrenia. Agonists at the allosteric glycine (Gly) binding site of the NMDA complex can promote NMDAr activity, a strategy that could provide therapeutic efficacy for the disorder. NMDAr antagonists like phencyclidine (PCP) can induce psychotic and dissociative symptoms similar to those observed in schizophrenia and are therefore widely used experimentally to impair NMDA neurotransmission in vivo.

OBJECTIVES

In the present study, we used pharmacological magnetic resonance imaging (phMRI) to investigate the modulatory effects of endogenous and exogenous agonists at the NMDAr Gly site on the spatiotemporal patterns of brain activation induced by acute PCP challenge in the rat. The drugs investigated were D: -serine, an endogenous agonist of the NMDAr Gly site, and SSR504734, a potent Gly transporter type 1 (GlyT-1) inhibitor that can potentiate NMDAr function by increasing synaptic levels of Gly.

RESULTS

Acute administration of PCP induced robust and sustained activation of discrete cortico-limbo-thalamic circuits. Pretreatment with D: -serine (1 g/kg) or SSR504734 (10 mg/kg) completely inhibited PCP-induced functional activation. This effect was accompanied by weak but sustained deactivation particularly in cortical areas.

CONCLUSIONS

These findings suggest that agents that stimulate NMDAr via Gly co-agonist site can potentiate NMDAr activity in the living brain and corroborate the potential for this class of drugs to provide selective enhancement of NMDAr neurotransmission in schizophrenia.

Differential effects of antipsychotic and glutamatergic agents on the phMRI response to phencyclidine

Acute administration of NMDA receptor (NMDAR) antagonists such as phencyclidine (PCP) or ketamine induces symptoms that closely resemble those of schizophrenia in humans, a finding that has led to the hypothesis that a decreased NMDAR function may be a predisposing or even causative factor in schizophrenia. However, the precise neuropharmacological mechanisms underlying these effects remain to be fully elucidated. Here, we applied pharmacological MRI (phMRI) to examine the brain circuitry underlying the psychotomimetic action of PCP in the anesthetized rat, and investigated how these functional changes are modulated by drugs that possess distinct pharmacological mechanisms. Acute administration of PCP (0.5 mg/kg i.v.) produced robust and sustained positive relative cerebral blood volume (rCBV) changes in discrete cortico-limbo-thalamic regions. Pretreatment with the selective D2 dopamine antagonist raclopride (0.3 mg/kg i.p.) did not significantly affect the rCBV response to PCP, while the atypical antipsychotic clozapine (5 mg/kg i.p.) produced region-dependent effects, with complete suppression of the rCBV response in the thalamus, and weaker attenuation of the response in cortical and hippocampal structures. The response to PCP was strongly suppressed in all regions by pretreatment with two drugs that can inhibit aberrant glutamatergic activity: the anticonvulsant lamotrigine (10 mg/kg i.p.) and the mGluR2/3 agonist LY354740 (10 mg/kg i.p.). Taken together, our findings corroborate the pivotal role of dysfunctional glutamatergic neurotransmission in the functional response elicited by PCP, while the lack of effect of raclopride argues against a primary role of dopamine D2 receptor activation in this process. Finally, the thalamic effect of clozapine could be key to elucidating the functional basis of its pharmacological action.

Drug-anaesthetic interaction in phMRI: the case of the psychotomimetic agent phencyclidine

Pharmacological magnetic resonance imaging (phMRI) provides a powerful means to map the effects of drugs on brain activity, with important applications in pharmacological research. However, phMRI studies in preclinical species are often conducted under general anaesthesia as a means to avoid head motion and to minimise the stress induced by the procedure. Under these conditions, the phMRI response to the drug of interest may be affected by interactions with the anaesthetic agent, with consequences for the interpretation of the data. Here, we have investigated the phMRI response to phencyclidine (PCP), an NMDA receptor blocker, in the halothane-anaesthetised rat for varying levels of anaesthesia and different PCP challenge doses. PCP induces psychotic-like symptoms in humans and laboratory animals and is widely applied as a pharmacological model of schizophrenia. However, PCP possesses anaesthetic properties per se, and its interactions with halothane might result in significant effects on the phMRI activation patterns. We observed two qualitatively different patterns of phMRI response. At 0.5 mg/kg iv PCP and 0.8% halothane maintenance anaesthesia, the lowest doses explored, an activation of discrete cortico-limbo-thalamic structures was observed, consistent with neuroimaging studies in humans and 2-deoxyglucose functional mapping in conscious animal models. However, higher anaesthetic concentrations or higher PCP challenge doses resulted in complete abolition of the positive response and in a widespread cortical deactivation (negative response). In the intermediate regime, we observed a dichotomic behaviour, with individual subjects showing one pattern or the other. These findings indicate a dose-dependent drug-anaesthetic interaction, with a complete reversal of the effects of PCP at higher challenge doses or HT concentrations.

Functional connectivity in the pharmacologically activated brain: resolving networks of correlated responses to d-amphetamine

We investigated the functional connectivity structure underlying the widespread relative cerebral blood volume (rCBV) response to d-amphetamine in the rat brain by systematically analyzing the intersubject correlations between the response amplitudes in 48 atlas-defined brain structures. A cluster analysis resolved three distinct networks of brain regions that exhibited closely coupled responses: one corresponding to primary dopamine projections from the midbrain to the striatum, a second consisting predominantly of forebrain cortical and basal ganglia regions that share a widespread correlation pattern resembling the univariate group response, and a third including structures in the periventricular dopamine system. These results suggest that different functional networks underlie the brain's response to d-amphetamine. This approach may provide important new insights regarding the central systems that underlie pharmacological action.

A multimodality investigation of cerebral hemodynamics and autoregulation in pharmacological MRI

Pharmacological MRI (phMRI) methods have been widely applied to assess the central hemodynamic response to pharmacological intervention as a surrogate for changes in the underlying neuronal activity. However, many psychoactive drugs can also affect cardiovascular parameters, including arterial blood pressure (BP). Abrupt changes in BP or the anesthetic agents used in preclinical phMRI may impair cerebral blood flow (CBF) autoregulation mechanisms, potentially introducing confounds in the phMRI response. Moreover, relative cerebral blood volume (rCBV), often measured in small-animal phMRI studies, may be sensitive to BP changes even in the presence of intact autoregulation. We applied laser Doppler flowmetry and MRI to measure changes in CBF and microvascular CBV induced by increasing doses of intravenous norepinephrine (NE) challenge in the halothane-anesthetized rat. NE is a potent vasopressor that does not cross the blood-brain barrier and mimics the rapid BP changes typically observed with acute drug challenges. We found that CBF autoregulation was maintained over a BP range of 60-120 mmHg. Under these conditions, no significant central rCBV responses were observed, suggesting that microvascular rCBV changes in response to abrupt changes in perfusion pressure are negligible within the autoregulatory range. Larger BP responses were accompanied by significant changes in both CBV and CBF that might confound the interpretation of phMRI results.

Pharmacological modulation of functional connectivity: the correlation structure underlying the phMRI response to d-amphetamine modified by selective dopamine D3 receptor antagonist SB277011A

Pharmacological MRI (phMRI) experiments utilise fMRI time series methods to map the central effect of pharmaceutical compounds. The typical univariate maps may, however, integrate the effects of several different neurotransmitter systems or underlying mechanisms. The results may thus be spatially and/or mechanistically nonspecific. Intersubject correlation analysis based on the phMRI response amplitude can more directly identify patterns of functional connectivity underlying the central effects of an acutely administered compound. In this article, we extend this approach to experiments where the effects of one compound in modulating the response to another are of interest. Specifically, we show a modulation of the correlation structure of a probe compound (d-amphetamine) by pretreatment with the selective dopamine D3 receptor antagonist SB277011A in the rat. The strongest modifications in the correlation patterns occurred in connection with the ventral tegmental area, the source of mesolimbic dopamine projections and a key substrate in the reward system.