More Homogeneous Capillary Flow and Oxygenation in Deeper Cortical Layers Correlate with Increased Oxygen Extraction

Our understanding of how capillary blood flow and oxygen distribute across cortical layers to meet the local metabolic demand is incomplete. We addressed this question by using two-photon imaging of resting-state microvascular oxygen partial pressure (PO2) and flow in the whisker barrel cortex in awake mice. Our measurements in layers I-V show that the capillary red-blood-cell flux and oxygenation heterogeneity, and the intracapillary resistance to oxygen delivery, all decrease with depth, reaching a minimum around layer IV, while the depth-dependent oxygen extraction fraction is increased in layer IV, where oxygen demand is presumably the highest. Our findings suggest that more homogeneous distribution of the physiological observables relevant to oxygen transport to tissue is an important part of the microvascular network adaptation to local brain metabolism. These results will inform the biophysical models of layer-specific cerebral oxygen delivery and consumption and improve our understanding of the diseases that affect cerebral microcirculation

Seizure epicenter depth and translaminar field potential synchrony underlie complex variations in tissue oxygenation during ictal initiation

Whether functional hyperemia during epileptic activity is adequate to meet the heightened metabolic demand of such events is controversial. Whereas some studies have demonstrated hyperoxia during ictal onsets, other work has reported transient hypoxic episodes that are spatially dependent on local surface microvasculature. Crucially, how laminar differences in ictal evolution can affect subsequent cerebrovascular responses has not been thus far investigated, and is likely significant in view of possible laminar-dependent neurovascular mechanisms and angioarchitecture. We addressed this open question using a novel multi-modal methodology enabling concurrent measurement of cortical tissue oxygenation, blood flow and hemoglobin concentration, alongside laminar recordings of neural activity, in a urethane anesthetized rat model of recurrent seizures induced by 4-aminopyridine. We reveal there to be a close relationship between seizure epicenter depth, translaminar LFP synchrony and tissue oxygenation during the early stages of recurrent seizures, whereby deep layer seizures are associated with decreased cross laminar synchrony and prolonged periods of hypoxia, and middle layer seizures are accompanied by increased cross-laminar synchrony and hyperoxia. Through comparison with functional activation by somatosensory stimulation and graded hypercapnia, we show that these seizure-related cerebrovascular responses occur in the presence of conserved neural-hemodynamic and blood flow-volume coupling. Our data provide new insights into the laminar dependency of seizure-related neurovascular responses, which may reconcile inconsistent observations of seizure-related hypoxia in the literature, and highlight a potential layer-dependent vulnerability that may contribute to the harmful effects of clinical recurrent seizures. The relevance of our findings to perfusion-related functional neuroimaging techniques in epilepsy are also discussed

OCT methods for capillary velocimetry

To date, two main categories of OCT techniques have been described for imaging hemodynamics: Doppler OCT and OCT angiography. Doppler OCT can measure axial velocity profiles and flow in arteries and veins, while OCT angiography can determine vascular morphology, tone, and presence or absence of red blood cell (RBC) perfusion. However, neither method can quantify RBC velocity in capillaries, where RBC flow is typically transverse to the probe beam and single-file. Here, we describe new methods that potentially address these limitations. Firstly, we describe a complex-valued OCT signal in terms of a static scattering component, dynamic scattering component, and noise. Secondly, we propose that the time scale of random fluctuations in the dynamic scattering component are related to red blood cell velocity. Analysis was performed along the slow axis of repeated B-scans to parallelize measurements. We correlate our purported velocity measurements against two-photon microscopy measurements of RBC velocity, and investigate changes during hypercapnia. Finally, we image the ischemic stroke penumbra during distal middle cerebral artery occlusion (dMCAO), where OCT velocimetry methods provide additional insight that is not afforded by either Doppler OCT or OCT angiography

Structural and functional brain imaging using extended-focus optical coherence tomography and microscopy

Neuroimaging techniques aim at revealing the anatomy and functional organisation of cerebral structures. Over the past decades, functional magnetic resonance imaging (fMRI) has revolutionized our understanding of human cerebral physiology through its ability to probe neural activity throughout the entire brain in a non-invasive fashion. Nevertheless, despite recent technological improvements, the spatial resolution of fMRI remains limited to a few hundreds of microns, restricting its use to macroscopic studies. Microscopic imaging solutions have been proposed to circumvent this limitation and have enabled revealing the existence of various cerebral structures, such as neuronal and vascular networks and their contribution to information processing and blood flow regulation within the brain. Optical imaging has proven, through its increased resolution and available contrast mechanisms, to be a valuable complement to fMRI for cellular-scale imaging. In this context, we demonstrate here the capabilities of an extension of optical coherence tomography, termed extended-focus optical coherence tomography (xf-OCT), in imaging cerebral structure and function at high resolution and very high acquisitions rates. Optical coherence tomography is an interferometric imaging technique using a low-coherence illumination source to provide fast, three-dimensional imaging of the back-scattering of tissue and cells. By multiplexing the interferometric ranging over several spectral channels, Fourier-domain OCT performs depth-resolved imaging at very high acquisition rates and high sensitivity. Increasing the lateral resolution of optical systems typically reduces the available depth-of-field and thus hampers this depth multiplexing advantage of OCT. Extended-focus systems aim at alleviating this trade-off between imaging depth and lateral resolution through the use of diffraction-less beams such as Bessel beams, providing high resolution imaging over large depths. The xf-OCT system therefore combines fast acquisition rates and high resolution, both characteristics required to image and study the structure and function of microscopic constituents of cerebral tissue. In this work, we performed functional brain imaging using the ability of xf-OCT to obtain quantita- tive measurements of blood flow in the brain. Changes in blood velocity evoked by neuronal activation were monitored and maps of hemodynamic activity were generated by adapting tools routinely used in fMRI to xf-OCT imaging. Additionally, three novel xf-OCT instruments are presented, wherein the advantages of different spectral ranges are exploited to reach specific imaging parameters. The increased contrast and resolution afforded by an illumination in the visible spectral range was used in two extended-focus optical coherence microscopy (xf-OCM) implementations for subcellular imaging of ex-vivo brain slices and cellular imaging of neurons, capillaries and myelinated axons in the superficial cortex in-vivo. Subsequently, an xf-OCT system is presented, operating in the infrared spectral range, wherein the reduced scattering enabled imaging the smallest capillaries deep in the murine cortex in-vivo

Functional Connectivity of the Rodent Brain Using Optical Imaging

RÉSUMÉ L'objectif de cette thèse de doctorat est d’appliquer la connectivité fonctionnelle dans une variété de modèles animaux, à l’aide de plusieurs techniques d’imagerie optique. Le cerveau, même au repos, montre une activité métabolique élevée : la corrélation des fluctuations spontanées lentes permet d’identifier des régions cérébrales distantes mais connectées; d’où le terme connectivité fonctionnelle. Les changements dans l’activité spontanée peuvent donner un aperçu des processus neuronaux qui comprennent la majorité de l’activité métabolique du cerveau, et constituent en conséquent une vaste source de changements reliés aux maladies. L’hémodynamique du cerveau peut être modifiée grâce à des affections neurovasculaires et avoir un effet sur l’activité au repos. Cette thèse vise la compréhension des changements de connectivité fonctionnelle induits par des maladies, à l’aide de l’imagerie optique fonctionnelle. Les techniques d’imagerie explorées dans les deux premières contributions de cette thèse sont l’Imagerie Optique Intrinsèque et l’Imagerie par Granularité Laser. Ensemble, elles peuvent estimer les changements de consommation d'oxygène, étroitement liés à l’activité neuronale. Ces techniques possèdent des résolutions temporelles et spatiales adéquates et bien adaptées pour imager la convexité du cortex cérébral. Dans le dernier article, une modalité d’imagerie en profondeur, la Tomographie Photoacoustique a été utilisée chez le rat nouveau-né. La Tomographie par Cohérence Optique et la Tomographie Laminaire Optique font également partie de la gamme des techniques d’imagerie développées et appliquées dans d’autres collaborations. La première partie des résultats mesure les changements de connectivité fonctionnelle dans un modèle d’activité épileptiforme aiguë chez le rongeur. Il y a des augmentations ainsi que des diminutions entre les corrélations homologues, avec une faible dépendance aux crises épileptiques. Ces changements suggèrent un découplage potentiel entre les paramètres hémodynamiques dans les réseaux au repos, en soulignant l’importance d’investiguer les réseaux épileptiques à l’aide de plusieurs mesures hémodynamiques indépendantes. La deuxième partie des travaux étudie un nouveau modèle de rigidité artérielle chez la souris : la calcification unilatérale de la carotide droite. L’analyse de connectivité basé sur les régions d’intérêt montre une tendance décroissante de corrélation homologue dans les cortex moteur et cingulum. L’analyse de graphes montre une randomisation des réseaux corticaux, ce qui suggère une perte de connectivité; plus spécifiquement, dans le cortex moteur ipsilateral à la carotide----------ABSTRACT The aim of this thesis is to apply functional connectivity in a variety of animal models, using several optical imaging modalities. Even at rest, the brain shows high metabolic activity: the correlation in slow spontaneous fluctuations identifies remotely connected areas of the brain; hence the term “functional connectivity”. Ongoing changes in spontaneous activity may provide insight into the neural processing that takes most of the brain metabolic activity, and so may provide a vast source of disease related changes. Brain hemodynamics may be modified during disease and affect resting-state activity. The thesis aims to better understand these changes in functional connectivity due to disease, using functional optical imaging. The optical imaging techniques explored in the first two contributions of this thesis are Optical Imaging of Intrinsic Signals and Laser Speckle Contrast Imaging, together they can estimate the metabolic rate of oxygen consumption, that closely parallels neural activity. They both have adequate spatial and temporal resolution and are well adapted to image the convexity of the mouse cortex. In the last article, a depth-sensitive modality called photoacoustic tomography was used in the newborn rat. Optical coherence tomography and laminar optical tomography were also part of the array of imaging techniques developed and applied in other collaborations. The first article of this work shows the changes in functional connectivity in an acute murine model of epileptiform activity. Homologous correlations are both increased and decreased with a small dependence on seizure duration. These changes suggest a potential decoupling between the hemodynamic parameters in resting-state networks, underlining the importance to investigate epileptic networks with several independent hemodynamic measures. The second study examines a novel murine model of arterial stiffness: the unilateral calcification of the right carotid. Seed-based connectivity analysis showed a decreasing trend of homologous correlation in the motor and cingulate cortices. Graph analyses showed a randomization of the cortex functional networks, suggesting a loss of connectivity, more specifically in the motor cortex ipsilateral to the treated carotid; however these changes are not reflected in differentiated metabolic estimates. Confounds remain due to the fact that carotid rigidification gives rise to neural decline in the hippocampus as well as unilateral alteration of vascular pulsatility; howeve

Cerebrovascular Dysfunction and Degeneration in Alzheimer’s Disease Pathophysiology

Alzheimer’s disease (AD) is a terminal illness and the most common form of dementia, which disproportionately affects the aged population. The pathophysiology of AD is characterized by neurodegeneration that slowly progresses, affecting regions of the brain that are involved in learning, memory, language, and executive function. In patients with the disease, early symptoms include non-disruptive forgetfulness that evolves into the inability to form new memories and ultimately the loss of autonomy at late stages. Histopathological hallmarks in the brain from patients with AD is the presence of amyloid-β (Aβ)-plaques and neurofibrillary tangles (NFT) deposited in the parenchyma. Since the discovery of these hallmarks, the majority of AD research has disproportionately focused on Aβ -plaques and NFT. Although the etiology of AD remains unknown, considerable advances have been made describing the cellular, molecular, and genetic contributions to the disease. Aging is the important risk factor for the development of AD, many other factors that increase the risk of developing AD later in life are vascular in nature. The function of the cardiovascular system is known to decline during healthy aging, and the same is true for the cerebrovasculature. Empirical evidence has demonstrated a decline cerebrovascular function in AD that exceeds the decline that occurs in healthy aging. Cerebrovascular dysfunction is the major contributor to the development of hypoperfusion and hypometabolism in patients diagnosed with AD. Cerebral amyloid angiopathy (CAA) is a neuropathological condition defined by the abnormal accumulation of Aβ on the walls of the cerebrovasculature. CAA occurs in as many as 90% of patients with AD and is implicated in the weakening of the walls of cerebral blood vessels. The occurrence of microhemorrhages, aneurysms, and microinfarctions are pathological manifestations associated with weakened walls of cerebral blood vessels in the brains of patients with confirmed AD. Noteworthy, cerebrovascular dysfunction, hypoperfusion, and hypometabolism occur before the onset of Aβ-plaque and NFT deposition in the brain of patients and animal models with AD. These findings provide a compelling basis that suggest a prominent role of dysfunctional cerebrovasculature in the etiology and for the progression of AD. Although the overwhelming evidence that implicates cerebrovascular dysfunction in AD, a thorough account of the changes that occur to the cerebrovasculature nor the mechanisms that drive these changes during the development and progression of AD has not been previously reported. The overarching goal(s) of this work are to; (1) provide a thorough description of the changes that occur to the cerebrovasculature during age and the progression of AD; (2) describe the mechanisms involved in cerebrovascular damage in AD; and (3) characterize the degeneration that results from cerebrovascular hypoperfusion. These overarching goals were achieved by completing five separate studies. Described in study 1, we investigated the effects of hypoxia on astrocytic mitochondria by assessing mitochondrial fission-fusion dynamics, reactive oxygen species production, synthesis of ATP, and mitophagy. Overall, we found a drastic mitochondrial network change that is triggered by metabolic crisis during hypoxia; these changes are followed by mitochondrial degradation and retraction of astrocytic extensions during reoxygenation. In study 2, we provide a novel model for the gradual development of cerebrovascular hypoperfusion in mice. Cerebrovascular hypoperfusion developed over 34-days by inserting an ameroid constrictor ring and microcoil bilaterally around the external carotid arteries. We investigated the neurodegenerative effects of hypoperfusion in mice by assessing both gray and white matter pathology. Histopathological analyses of the brain revealed neuronal and axonal degeneration as well as necrotic lesions. The most severely affected regions were located in the hippocampus and corpus callosum. Described in study 3, we performed a series of experiments to investigate the effects of Aβ on cerebrovascular endothelial cells. In this study, we focused on characterizing the changes to mitochondrial oxidative phosphorylation, superoxide production, mitochondrial calcium, ATP synthesis, and endothelial cell death. These results describe a mechanism for mitochondrial degeneration caused by the production of mitochondrial superoxide, which was driven by increased mitochondrial Ca2+ uptake. We found that persistent superoxide production injures mitochondria and disrupts electron transport in cerebrovascular endothelial cells. In study 4, we developed a method to evaluate the cerebrovasculature of the whole-brain and constructed analyses to assess the angioarchitecture. We used vascular corrosion casting method to replicate the cerebrovasculature in adult mice and used MicroCT to acquire volumetric imaging data of the cerebrovascular network at a resolution required to investigate the microvasculature. Our analyses of the cerebrovasculature evaluated the morphology, topology, and organization of the angioarchitecture. With these developments, we investigated the effects of age and progression of disease on the cerebrovasculature in wild type mice and the triple transgenic mouse model of AD. Study 5 provides data describing degenerative changes to the microvascular network that progress with age in the triple transgenic mouse model of AD. These changes to the microvasculature occurred early, before the onset of Aβ-plaque deposition and NFT development. Overall, this body of work provides evidence of an early cerebrovascular disruption in the etiology of AD that progresses with age. Aβ mediates early cerebrovascular damage through direct interaction with vascular endothelial cells. Microvascular degeneration can lead to hypoperfusion which damages both gray and white matter. Hypoperfusion-associated hypoxia may mediate parenchymal damage by disrupting mitochondrial fission-fusion dynamics and enhancing mitophagy. These data provide a basis for the development of novel therapeutic strategies that target the changes to the cerebrovasculature for the treatment of AD. These observations may substantiate a prophylactic strategy for the treatment of AD by preventing the initial factors that lead to compromised cerebrovasculature

Pericyte remodeling is deficient in the aged brain and contributes to impaired capillary flow and structure

Deterioration of brain capillary flow and architecture is a hallmark of aging and dementia. It remains unclear how loss of brain pericytes in these conditions contributes to capillary dysfunction. Here, we conduct cause-and-effect studies by optically ablating pericytes in adult and aged mice in vivo. Focal pericyte loss induces capillary dilation without blood-brain barrier disruption. These abnormal dilations are exacerbated in the aged brain, and result in increased flow heterogeneity in capillary networks. A subset of affected capillaries experience reduced perfusion due to flow steal. Some capillaries stall in flow and regress, leading to loss of capillary connectivity. Remodeling of neighboring pericytes restores endothelial coverage and vascular tone within days. Pericyte remodeling is slower in the aged brain, resulting in regions of persistent capillary dilation. These findings link pericyte loss to disruption of capillary flow and structure. They also identify pericyte remodeling as a therapeutic target to preserve capillary flow dynamics

Anatomical Modeling of Cerebral Microvascular Structures: Application to Identify Biomarkers of Microstrokes

Les réseaux microvasculaires corticaux sont responsables du transport de l’oxygène et des substrats énergétiques vers les neurones. Ces réseaux réagissent dynamiquement aux demandes énergétiques lors d’une activation neuronale par le biais du couplage neurovasculaire. Afin d’élucider le rôle de la composante microvasculaire dans ce processus de couplage, l’utilisation de la modélisation in-formatique pourrait se révéler un élément clé. Cependant, la manque de méthodologies de calcul appropriées et entièrement automatisées pour modéliser et caractériser les réseaux microvasculaires reste l’un des principaux obstacles. Le développement d’une solution entièrement automatisée est donc important pour des explorations plus avancées, notamment pour quantifier l’impact des mal-formations vasculaires associées à de nombreuses maladies cérébrovasculaires. Une observation courante dans l’ensemble des troubles neurovasculaires est la formation de micro-blocages vascu-laires cérébraux (mAVC) dans les artérioles pénétrantes de la surface piale. De récents travaux ont démontré l’impact de ces événements microscopiques sur la fonction cérébrale. Par conséquent, il est d’une importance vitale de développer une approche non invasive et comparative pour identifier leur présence dans un cadre clinique. Dans cette thèse,un pipeline de traitement entièrement automatisé est proposé pour aborder le prob-lème de la modélisation anatomique microvasculaire. La méthode de modélisation consiste en un réseau de neurones entièrement convolutif pour segmenter les capillaires sanguins, un générateur de modèle de surface 3D et un algorithme de contraction de la géométrie pour produire des mod-èles graphiques vasculaires ne comportant pas de connections multiples. Une amélioration de ce pipeline est développée plus tard pour alléger l’exigence de maillage lors de la phase de représen-tation graphique. Un nouveau schéma permettant de générer un modèle de graphe est développé avec des exigences d’entrée assouplies et permettant de retenir les informations sur les rayons des vaisseaux. Il est inspiré de graphes géométriques déformants construits en respectant les morpholo-gies vasculaires au lieu de maillages de surface. Un mécanisme pour supprimer la structure initiale du graphe à chaque exécution est implémenté avec un critère de convergence pour arrêter le pro-cessus. Une phase de raffinement est introduite pour obtenir des modèles vasculaires finaux. La modélisation informatique développée est ensuite appliquée pour simuler les signatures IRM po-tentielles de mAVC, combinant le marquage de spin artériel (ASL) et l’imagerie multidirectionnelle pondérée en diffusion (DWI). L’hypothèse est basée sur des observations récentes démontrant une réorientation radiale de la microvascularisation dans la périphérie du mAVC lors de la récupéra-tion chez la souris. Des lits capillaires synthétiques, orientés aléatoirement et radialement, et des angiogrammes de tomographie par cohérence optique (OCT), acquis dans le cortex de souris (n = 5) avant et après l’induction d’une photothrombose ciblée, sont analysés. Les graphes vasculaires informatiques sont exploités dans un simulateur 3D Monte-Carlo pour caractériser la réponse par résonance magnétique (MR), tout en considérant les effets des perturbations du champ magnétique causées par la désoxyhémoglobine, et l’advection et la diffusion des spins nucléaires. Le pipeline graphique proposé est validé sur des angiographies synthétiques et réelles acquises avec différentes modalités d’imagerie. Comparé à d’autres méthodes effectuées dans le milieu de la recherche, les expériences indiquent que le schéma proposé produit des taux d’erreur géométriques et topologiques amoindris sur divers angiogrammes. L’évaluation confirme également l’efficacité de la méthode proposée en fournissant des modèles représentatifs qui capturent tous les aspects anatomiques des structures vasculaires. Ensuite, afin de trouver des signatures de mAVC basées sur le signal IRM, la modélisation vasculaire proposée est exploitée pour quantifier le rapport de perte de signal intravoxel minimal lors de l’application de plusieurs directions de gradient, à des paramètres de séquence variables avec et sans ASL. Avec l’ASL, les résultats démontrent une dif-férence significative (p <0,05) entre le signal calculé avant et 3 semaines après la photothrombose. La puissance statistique a encore augmenté (p <0,005) en utilisant des angiogrammes capturés à la semaine suivante. Sans ASL, aucun changement de signal significatif n’est trouvé. Des rapports plus élevés sont obtenus à des intensités de champ magnétique plus faibles (par exemple, B0 = 3) et une lecture TE plus courte (<16 ms). Cette étude suggère que les mAVC pourraient être carac-térisés par des séquences ASL-DWI, et fournirait les informations nécessaires pour les validations expérimentales postérieures et les futurs essais comparatifs.----------ABSTRACT Cortical microvascular networks are responsible for carrying the necessary oxygen and energy substrates to our neurons. These networks react to the dynamic energy demands during neuronal activation through the process of neurovascular coupling. A key element in elucidating the role of the microvascular component in the brain is through computational modeling. However, the lack of fully-automated computational frameworks to model and characterize these microvascular net-works remains one of the main obstacles. Developing a fully-automated solution is thus substantial for further explorations, especially to quantify the impact of cerebrovascular malformations associ-ated with many cerebrovascular diseases. A common pathogenic outcome in a set of neurovascular disorders is the formation of microstrokes, i.e., micro occlusions in penetrating arterioles descend-ing from the pial surface. Recent experiments have demonstrated the impact of these microscopic events on brain function. Hence, it is of vital importance to develop a non-invasive and translatable approach to identify their presence in a clinical setting. In this thesis, a fully automatic processing pipeline to address the problem of microvascular anatom-ical modeling is proposed. The modeling scheme consists of a fully-convolutional neural network to segment microvessels, a 3D surface model generator and a geometry contraction algorithm to produce vascular graphical models with a single connected component. An improvement on this pipeline is developed later to alleviate the requirement of water-tight surface meshes as inputs to the graphing phase. The novel graphing scheme works with relaxed input requirements and intrin-sically captures vessel radii information, based on deforming geometric graphs constructed within vascular boundaries instead of surface meshes. A mechanism to decimate the initial graph struc-ture at each run is formulated with a convergence criterion to stop the process. A refinement phase is introduced to obtain final vascular models. The developed computational modeling is then ap-plied to simulate potential MRI signatures of microstrokes, combining arterial spin labeling (ASL) and multi-directional diffusion-weighted imaging (DWI). The hypothesis is driven based on recent observations demonstrating a radial reorientation of microvasculature around the micro-infarction locus during recovery in mice. Synthetic capillary beds, randomly- and radially oriented, and op-tical coherence tomography (OCT) angiograms, acquired in the barrel cortex of mice (n=5) before and after inducing targeted photothrombosis, are analyzed. The computational vascular graphs are exploited within a 3D Monte-Carlo simulator to characterize the magnetic resonance (MR) re-sponse, encompassing the effects of magnetic field perturbations caused by deoxyhemoglobin, and the advection and diffusion of the nuclear spins. The proposed graphing pipeline is validated on both synthetic and real angiograms acquired with different imaging modalities. Compared to other efficient and state-of-the-art graphing schemes, the experiments indicate that the proposed scheme produces the lowest geometric and topological error rates on various angiograms. The evaluation also confirms the efficiency of the proposed scheme in providing representative models that capture all anatomical aspects of vascular struc-tures. Next, searching for MRI-based signatures of microstokes, the proposed vascular modeling is exploited to quantify the minimal intravoxel signal loss ratio when applying multiple gradient di-rections, at varying sequence parameters with and without ASL. With ASL, the results demonstrate a significant difference (p<0.05) between the signal-ratios computed at baseline and 3 weeks after photothrombosis. The statistical power further increased (p<0.005) using angiograms captured at week 4. Without ASL, no reliable signal change is found. Higher ratios with improved significance are achieved at low magnetic field strengths (e.g., at 3 Tesla) and shorter readout TE (<16 ms). This study suggests that microstrokes might be characterized through ASL-DWI sequences, and provides necessary insights for posterior experimental validations, and ultimately, future transla-tional trials

Depth-Dependent Physiological Modulators of the BOLD Response in the Human Motor Cortex

This dissertation proposes a set of methods for improving spatial localization of cerebral metabolic changes using functional magnetic resonance imaging (fMRI). Blood oxygen level dependent (BOLD) fMRI estabilished itself as the most frequently used technique for mapping brain activity in humans. It is non-invasive and allows to obtain information about brain oxygenation changes in a few minutes. It was discovered in 1990 and, since then, it contributed enormously to the developments in neuroscientific research. Nevertheless, the BOLD contrast suffers from inherent limitations. This comes from the fact that the observed response is the result of a complex interplay between cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen consumption (CMRO2) and has a strong dependency on baseline blood volume and oxygenation. Therefore, the observed response is mislocalized from the site where the metabolic activity takes place and it is subject to high variability across experiments due to normal brain physiology. Since the peak of BOLD changes can be as much as 4 mm apart from the site of metabolic changes, the problem of spatial mislocalization is particularly constraining at submillimeter resolution. Three methods are proposed in this work in order to overcome this limitation and make data more comparable. The first method involves a modification of an estabilished model for calibration of BOLD responses (the dilution model), in order to render it applicable at higher resolutions. The second method proposes a model-free scaling of the BOLD response, based on spatial normalization by a purely vascular response pattern. The third method takes into account the hypothesis that the cortical vasculature could act as a low-pass filter for BOLD fluctuations as the blood is carried downstream, and investigates differences in frequency composition of cortical laminae. All methods are described and tested on a depth-dependent scale in the human motor cortex

Cortical Layer-Dependent Hemodynamic Regulation Investigated by Functional Magnetic Resonance Imaging

Functional magnetic resonance imaging (fMRI) is currently one of the most widely used non-invasive neuroimaging modalities for mapping brain activation. Techniques such as blood oxygenation level dependent (BOLD) fMRI or cerebral blood volume (CBV)-weighted fMRI are based on the assumption that hemodynamic responses are tightly regulated by neural activity. However, the relationship between fMRI responses and neural activity is still unclear. To investigate this relationship, the unique properties of temporal frequency tuning of primary visual cortex neurons was used as a model since it can be used to separate the neural input and output activities of this area. During moving grating stimuli of 1, 2, 10 and 20 Hz temporal frequencies, two fMRI studies, areal and laminar studies, were conducted with different spatial resolution in a 9.4-T Varian spectrometer. In areal studies, BOLD fMRI was able to detect the difference in tuning properties between area 17 (A17), area 18 (A18) and lateral geniculate nucleus. In A17, the BOLD tuning curve seemed to reflect the local field potential (LFP) low frequency band (<12 Hz) rather than spiking activity and LFP gamma band (25-90 Hz). In laminar studies, a high spatial resolution protocol was adopted to resolve the different cortical layers in A17. In addition to BOLD fMRI, CBV-weighted fMRI was performed to eliminate the contamination from the superficial draining veins. These results showed that BOLD and CBV tuning curves do not reflect the underlying spiking activity or the LFP activity at infragranular layers (the bottom layer of three cortical layers). This implies that the hemodynamic response may not be regulated on a laminar level. Therefore, caution should be taken when interpreting BOLD responses as the sole indicator of different aspects of neural activity in areal and laminar scales