Quantifying the Microvascular Origin of BOLD-fMRI from First Principles with Two-Photon Microscopy and an Oxygen-Sensitive Nanoprobe

The blood oxygenation level-dependent (BOLD) contrast is widely used in functional magnetic resonance imaging (fMRI) studies aimed at investigating neuronal activity. However, the BOLD signal reflects changes in blood volume and oxygenation rather than neuronal activity per se. Therefore, understanding the transformation of microscopic vascular behavior into macroscopic BOLD signals is at the foundation of physiologically informed noninvasive neuroimaging. Here, we use oxygen-sensitive two-photon microscopy to measure the BOLD-relevant microvascular physiology occurring within a typical rodent fMRI voxel and predict the BOLD signal from first principles using those measurements. The predictive power of the approach is illustrated by quantifying variations in the BOLD signal induced by the morphological folding of the human cortex. This framework is then used to quantify the contribution of individual vascular compartments and other factors to the BOLD signal for different magnet strengths and pulse sequences.National Institutes of Health (U.S.) (Grant P41RR14075)National Institutes of Health (U.S.) (Grant R01NS067050)National Institutes of Health (U.S.) (Grant R01NS057198)National Institutes of Health (U.S.) (Grant R01EB000790)American Heart Association (Grant 11SDG7600037)Advanced Multimodal NeuroImaging Training Program (R90DA023427

Three-Dimensional Iron Oxide Nanoparticle-Based Contrast-Enhanced Magnetic Resonance Imaging for Characterization of Cerebral Arteriogenesis in the Mouse Neocortex

Purpose: Subsurface blood vessels in the cerebral cortex have been identified as a bottleneck in cerebral perfusion with the potential for collateral remodeling. However, valid techniques for non-invasive, longitudinal characterization of neocortical microvessels are still lacking. In this study, we validated contrast-enhanced magnetic resonance imaging (CE-MRI) for in vivo characterization of vascular changes in a model of spontaneous collateral outgrowth following chronic cerebral hypoperfusion. Methods: C57BL/6J mice were randomly assigned to unilateral internal carotid artery occlusion or sham surgery and after 21 days, CE-MRI based on T2*-weighted imaging was performed using ultra-small superparamagnetic iron oxide nanoparticles to obtain subtraction angiographies and steady-state cerebral blood volume (ss-CBV) maps. First pass dynamic susceptibility contrast MRI (DSC-MRI) was performed for internal validation of ss-CBV. Further validation at the histological level was provided by ex vivo serial two-photon tomography (STP). Results: Qualitatively, an increase in vessel density was observed on CE-MRI subtraction angiographies following occlusion; however, a quantitative vessel tracing analysis was prone to errors in our model. Measurements of ss-CBV reliably identified an increase in cortical vasculature, validated by DSC-MRI and STP. Conclusion: Iron oxide nanoparticle-based ss-CBV serves as a robust, non-invasive imaging surrogate marker for neocortical vessels, with the potential to reduce and refine preclinical models targeting the development and outgrowth of cerebral collateralization

A simulation study investigating potential diffusion-based MRI signatures of microstrokes

ABSTRACT: Recent studies suggested that cerebrovascular micro-occlusions, i.e. microstokes, could lead to ischemic tissue infarctions and cognitive deficits. Due to their small size, identifying measurable biomarkers of these microvascular lesions remains a major challenge. This work aims to simulate potential MRI signatures combining arterial spin labeling (ASL) and multi-directional diffusion-weighted imaging (DWI). Driving our hypothesis are 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 optical coherence tomography (OCT) angiograms, acquired in the barrel cortex of mice (n = 5) before and after inducing targeted photothrombosis, were analyzed. Computational vascular graphs combined with a 3D Monte-Carlo simulator were used to characterize the magnetic resonance (MR) response, encompassing the effects of magnetic field perturbations caused by deoxyhemoglobin, and the advection and diffusion of the nuclear spins. We quantified the minimal intravoxel signal loss ratio when applying multiple gradient directions, at varying sequence parameters with and without ASL. With ASL, our 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 measured at week 4. Without ASL, no reliable signal change was found. We found that higher ratios, and accordingly improved significance, were achieved at lower magnetic field strengths (e.g., B0 = 3T) and shorter echo time TE (< 16 ms). Our simulations suggest that microstrokes might be characterized through ASL-DWI sequence, providing necessary insights for posterior experimental validations, and ultimately, future translational trials

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

Disambiguating the role of blood flow and global signal with partial information decomposition

Global signal (GS) is an ubiquitous construct in resting state functional magnetic resonance imaging (rs-fMRI), associated to nuisance, but containing by definition most of the neuronal signal. Global signal regression (GSR) effectively removes the impact of physiological noise and other artifacts, but at the same time it alters correlational patterns in unpredicted ways. Performing GSR taking into account the underlying physiology (mainly the blood arrival time) has been proven to be beneficial. From these observations we aimed to: 1) characterize the effect of GSR on network-level functional connectivity in a large dataset; 2) assess the complementary role of global signal and vessels; and 3) use the framework of partial information decomposition to further look into the joint dynamics of the global signal and vessels, and their respective influence on the dynamics of cortical areas. We observe that GSR affects intrinsic connectivity networks in the connectome in a non-uniform way. Furthermore, by estimating the predictive information of blood flow and the global signal using partial information decomposition, we observe that both signals are present in different amounts across intrinsic connectivity networks. Simulations showed that differences in blood arrival time can largely explain this phenomenon, while using hemodynamic and calcium mouse recordings we were able to confirm the presence of vascular effects, as calcium recordings lack hemodynamic information. With these results we confirm network-specific effects of GSR and the importance of taking blood flow into account for improving de-noising methods. Additionally, and beyond the mere issue of data denoising, we quantify the diverse and complementary effect of global and vessel BOLD signals on the dynamics of cortical areas

Current Understanding of the Anatomy, Physiology, and Magnetic Resonance Imaging of Neurofluids: Update From the 2022 "ISMRM Imaging Neurofluids Study group" Workshop in Rome

Neurofluids is a term introduced to define all fluids in the brain and spine such as blood, cerebrospinal fluid, and interstitial fluid. Neuroscientists in the past millennium have steadily identified the several different fluid environments in the brain and spine that interact in a synchronized harmonious manner to assure a healthy microenvironment required for optimal neuroglial function. Neuroanatomists and biochemists have provided an incredible wealth of evidence revealing the anatomy of perivascular spaces, meninges and glia and their role in drainage of neuronal waste products. Human studies have been limited due to the restricted availability of noninvasive imaging modalities that can provide a high spatiotemporal depiction of the brain neurofluids. Therefore, animal studies have been key in advancing our knowledge of the temporal and spatial dynamics of fluids, for example, by injecting tracers with different molecular weights. Such studies have sparked interest to identify possible disruptions to neurofluids dynamics in human diseases such as small vessel disease, cerebral amyloid angiopathy, and dementia. However, key differences between rodent and human physiology should be considered when extrapolating these findings to understand the human brain. An increasing armamentarium of noninvasive MRI techniques is being built to identify markers of altered drainage pathways. During the three-day workshop organized by the International Society of Magnetic Resonance in Medicine that was held in Rome in September 2022, several of these concepts were discussed by a distinguished international faculty to lay the basis of what is known and where we still lack evidence. We envision that in the next decade, MRI will allow imaging of the physiology of neurofluid dynamics and drainage pathways in the human brain to identify true pathological processes underlying disease and to discover new avenues for early diagnoses and treatments including drug delivery. Evidence level: 1. Technical Efficacy: Stage 3

Biophysical Simulation of the Functional Magnetic Resonance Signal Formation in Realistic Neurovascular Networks

The human brain is one of the most complex living systems. The scientific study of the brain’s anatomy and neurophysiology are fundamental to understand the basic principles of mental processes such as cognition and behavior. For this reason, in the endeavor to investigate the underlying neural mechanisms that drive these processes, the neuroscientific research has employed all the available technological resources and methodologies from different fields like anatomy, histology, electrophysiology, neurobiology, etc. Likewise, great advances have been provided by neuroimaging techniques such as PET and MRI, in order to comprehend the neural activity and the metabolic reactions that occur in the central nervous system. In particular, functional MRI provides an indirect measurement of the neural activity throughout local hemodynamic changes, thus related to the neurovascular coupling, as a response to a particular task-evoked stimulus. This MR signal behavior modulated by oxygenation level changes is better known as the BOLD signal. Although progress has been done in order to understand the BOLD signal change under well-defined nonrealistic vascular geometries, on the other hand, realistic neurovascular networks might give valuable information to resolve the influence on the BOLD signal evolution from a particular vascular tissue and specific hemodynamic responses. In order to extend the analysis of the BOLD signal change obtained by randomly oriented cylinders and spheres, throughout this thesis, the geometrical features of a realistic neurovascular network as well as the biophysical effects related to the hemodynamic response and thermal motion were investigated by means of Monte Carlo simulations in pursuance to resolve the functional MR signal formation. In the Introduction of this thesis, I made a small recapitulation on the MR physics and spin dynamics; magnetic susceptibility and thermal motion as crucial modulators of the BOLD signal behavior. In addition, a summary of the problem and the aims of the project. Therefore, I described the importance of the use of the Monte Carlo method to calculate the MR signal under nonrealistic vascular models. I summarized the seminal analytical and numerical results that provide important insights to characterize the main parameters that influence the MR signal formation. Finally, I described the importance of the use of realistic neurovascular structures in order to disentangle the specific tissue contribution and the direct impact on the BOLD signal change

Biodistribution, biocompatibility and targeted accumulation of magnetic nanoporous silica nanoparticles as drug carrier in orthopedics

Background: In orthopedics, the treatment of implant-associated infections represents a high challenge. Especially, potent antibacterial effects at implant surfaces can only be achieved by the use of high doses of antibiotics, and still often fail. Drug-loaded magnetic nanoparticles are very promising for local selective therapy, enabling lower systemic antibiotic doses and reducing adverse side effects. The idea of the following study was the local accumulation of such nanoparticles by an externally applied magnetic field combined with a magnetizable implant. The examination of the biodistribution of the nanoparticles, their effective accumulation at the implant and possible adverse side effects were the focus. In a BALB/c mouse model (n = 50) ferritic steel 1.4521 and Ti90Al6V4 (control) implants were inserted subcutaneously at the hindlimbs. Afterwards, magnetic nanoporous silica nanoparticles (MNPSNPs), modified with rhodamine B isothiocyanate and polyethylene glycol-silane (PEG), were administered intravenously. Directly/1/7/21/42 day(s) after subsequent application of a magnetic field gradient produced by an electromagnet, the nanoparticle biodistribution was evaluated by smear samples, histology and multiphoton microscopy of organs. Additionally, a pathohistological examination was performed. Accumulation on and around implants was evaluated by droplet samples and histology. Results: Clinical and histological examinations showed no MNPSNP-associated changes in mice at all investigated time points. Although PEGylated, MNPSNPs were mainly trapped in lung, liver, and spleen. Over time, they showed two distributional patterns: early significant drops in blood, lung, and kidney and slow decreases in liver and spleen. The accumulation of MNPSNPs on the magnetizable implant and in its area was very low with no significant differences towards the control. Conclusion: Despite massive nanoparticle capture by the mononuclear phagocyte system, no significant pathomorphological alterations were found in affected organs. This shows good biocompatibility of MNPSNPs after intravenous administration. The organ uptake led to insufficient availability of MNPSNPs in the implant region. For that reason, among others, the nanoparticles did not achieve targeted accumulation in the desired way, manifesting future research need. However, with different conditions and dimensions in humans and further modifications of the nanoparticles, this principle should enable reaching magnetizable implant surfaces at any time in any body region for a therapeutic reason. © 2020 The Author(s)

Biophysical modeling of hemodynamic-based neuroimaging techniques

Thesis (Ph. D. in Medical Engineering and Medical Physics)--Harvard-MIT Program in Health Sciences and Technology, 2013.Cataloged from PDF version of thesis.Includes bibliographical references (pages 163-182).Two different hemodynamic-based neuroimaging techniques were studied in this work. Near-Infrared Spectroscopy (NIRS) is a promising technique to measure cerebral hemodynamics in a clinical setting due to its potential for continuous monitoring. However, the presence of strong systemic interference in the signal significantly limits our ability to recover the hemodynamic response without averaging tens of trials. Developing a new methodology to clean the NIRS signal from systemic interference and isolate the cortical signal would therefore significantly increase our ability to recover the hemodynamic response opening the door for clinical NIRS studies such as epilepsy. Toward this goal, a new method based on multi-distance measurements and state-space modeling was developed and further optimized to remove systemic physiological oscillations contaminating the NIRS signal. Furthermore, the cortical and pial contributions to the NIRS signal were quantified using a new multimodal regression analysis. Functional Magnetic Resonance Imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) response has become the method of choice for exploring brain function, and yet the physiological basis of this technique is still poorly understood. Despite the effort, a detailed and validated model relating the signal measured to the physiological changes occurring in the cortical tissue is still lacking. Modeling the BOLD signal is challenging because of the difficulty to take into account the complex morphology of the cortical microvasculature, the distribution of oxygen in those microvessels and its dynamics during neuronal activation. Here, we overcome this difficulty by performing Monte Carlo simulations over real microvascular networks and oxygen distributions measured in vivo on rodents, at rest and during forepaw stimulation, using two-photon microscopy. Our model reveals for the first time the specific contribution of individual vascular compartment to the BOLD signal, for different field strengths and different cortical orientations. Our model makes a new prediction: the amplitude of the BOLD signal produced by a given physiological change during neuronal activation depends on the spatial orientation of the cortical region in the MRI scanner. This occurs because veins are preferentially oriented either perpendicular or parallel to the cortical surface in the gray matter.by Louis Gagnon.Ph.D.in Medical Engineering and Medical Physic

BOLD signal physiology: Models and applications

The BOLD contrast mechanism has a complex relationship with functional brain activity, oxygen metabolism, and neurovascular factors. Accurate interpretation of the BOLD signal for neuroscience and clinical applications necessitates a clear understanding of the sources of BOLD contrast and its relationship to underlying physiology. This review describes the physiological components that contribute to the BOLD signal and the steady-state calibrated BOLD models that enable quantification of functional changes with a separate challenge paradigm. The principles derived from these biophysical models are then used to interpret BOLD measurements in different neurological disorders in the presence of confounding vascular factors related to disease