Laminar differences in functional oxygen metabolism in monkey visual cortex measured with calibrated fMRI.

Blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) of cortical layers relies on the hemodynamic response and is biased toward large veins on the cortical surface. Functional changes in the cerebral metabolic rate of oxygen (ΔCMRO2) may reflect neural cortical function better than BOLD fMRI, but it is unknown whether the calibrated BOLD model for functional CMRO2 measurement remains valid at high resolution. Here, we measure laminar ΔCMRO2 elicited by visual stimulation in macaque primary visual cortex (V1) and find that ΔCMRO2 peaks in the middle of the cortex, in agreement with autoradiographic measures of metabolism. ΔCMRO2 values in gray matter are similar as found previously. Reductions in CMRO2 are associated with veins at the cortical surface, suggesting that techniques for vein removal may improve the accuracy of the model at very high resolution. However, our results show feasibility of laminar ΔCMRO2 measurement, providing a physiologically meaningful metric of laminar functional metabolism.

Validation of a patient-specific one-dimensional model of the systemic arterial tree.

The aim of this study is to develop and validate a patient-specific distributed model of the systemic arterial tree. This model is built using geometric and hemodynamic data measured on a specific person and validated with noninvasive measurements of flow and pressure on the same person, providing thus a patient-specific model and validation. The systemic arterial tree geometry was obtained from MR angiographic measurements. A nonlinear viscoelastic constitutive law for the arterial wall is considered. Arterial wall distensibility is based on literature data and adapted to match the wave propagation velocity of the main arteries of the specific subject, which were estimated by pressure waves traveling time. The intimal shear stress is modeled using the Witzig-Womersley theory. Blood pressure is measured using applanation tonometry and flow rate using transcranial ultrasound and phase-contrast-MRI. The model predicts pressure and flow waveforms in good qualitative and quantitative agreement with the in vivo measurements, in terms of wave shape and specific wave features. Comparison with a generic one-dimensional model shows that the patient-specific model better predicts pressure and flow at specific arterial sites. These results obtained let us conclude that a patient-specific one-dimensional model of the arterial tree is able to predict well pressure and flow waveforms in the main systemic circulation, whereas this is not always the case for a generic one-dimensional model.

fMRI at High Spatial Resolution: Implications for BOLD-Models.

As high-resolution functional magnetic resonance imaging (fMRI) and fMRI of cortical layers become more widely used, the question how well high-resolution fMRI signals reflect the underlying neural processing, and how to interpret laminar fMRI data becomes more and more relevant. High-resolution fMRI has shown laminar differences in cerebral blood flow (CBF), volume (CBV), and neurovascular coupling. Features and processes that were previously lumped into a single voxel become spatially distinct at high resolution. These features can be vascular compartments such as veins, arteries, and capillaries, or cortical layers and columns, which can have differences in metabolism. Mesoscopic models of the blood oxygenation level dependent (BOLD) response therefore need to be expanded, for instance, to incorporate laminar differences in the coupling between neural activity, metabolism and the hemodynamic response. Here we discuss biological and methodological factors that affect the modeling and interpretation of high-resolution fMRI data. We also illustrate with examples from neuropharmacology and the negative BOLD response how combining BOLD with CBF- and CBV-based fMRI methods can provide additional information about neurovascular coupling, and can aid modeling and interpretation of high-resolution fMRI.

Haemodynamics in the mouse aortic arch computed from MRI-derived velocities at the aortic root.

Mice are widely used to investigate atherogenesis, which is known to be influenced by stresses related to blood flow. However, numerical characterization of the haemodynamic environment in the commonly studied aortic arch has hitherto been based on idealizations of inflow into the aorta. Our purpose in this work was to numerically characterize the haemodynamic environment in the mouse aortic arch using measured inflow velocities, and to relate the resulting shear stress patterns to known locations of high- and low-lesion prevalence. Blood flow velocities were measured in the aortic root of C57/BL6 mice using phase-contrast MRI. Arterial geometries were obtained by micro-CT of corrosion casts. These data were used to compute blood flow and wall shear stress (WSS) patterns in the arch. WSS profiles computed using realistic and idealized aortic root velocities differed significantly. An unexpected finding was that average WSS in the high-lesion-probability region on the inner wall was actually higher than the WSS in the low-probability region on the outer wall. Future studies of mouse aortic arch haemodynamics should avoid the use of idealized inflow velocity profiles. Lesion formation does not seem to uniquely associate with low or oscillating WSS in this segment, suggesting that other factors may also play a role in lesion localization.