Quantitative theory for the transverse relaxation time of blood water

Quantitative BOLD (qBOLD) imaging of oxygen metabolism and blood oxygenation in the human body: A scoping review

PURPOSE

There are many approaches to the quantitative BOLD (qBOLD) technique described in the literature, differing in pulse sequences, MRI parameters and data processing. Thus, in this review, we summarized the acquisition methods, approaches used for oxygenation quantification and clinical populations investigated.

METHODS

Three databases were systematically searched (Medline, Embase, and Web of Science) for published research that used qBOLD methods for quantification of oxygen metabolism. Data extraction and synthesis were performed by one author and reviewed by a second author.

RESULTS

A total of 93 relevant papers were identified. Acquisition strategies were summarized, and oxygenation parameters were found to have been investigated in many pathologies such as steno-occlusive diseases, stroke, glioma, and multiple sclerosis disease.

CONCLUSION

A summary of qBOLD approaches for oxygenation measurements and applications could help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

Imaging arterial and venous vessels using Iron Dextran enhanced multi-echo 3D gradient echo MRI at 7T

Iron Dextran is a widely used iron oxide compound to treat iron-deficiency anemia patients in the clinic. Similar to other iron oxide compounds such as Ferumoxytol, it can also be used off-label as an intravascular magnetic resonance imaging (MRI) contrast agent due to its strong iron-induced T2 and T2* shortening effects. In this study, we seek to evaluate the feasibility of using Iron Dextran enhanced multi-echo susceptibility weighted imaging (SWI) MRI at 7T to image arterial and venous blood vessels in the human brain. Phantom experiments were performed to measure the r2* relaxivity for Iron Dextran in blood, based on which the SWI sequence was optimized. Pre- and post-infusion MR images were acquired in human subjects from which maps of arteries and veins were extracted. The post-contrast SWI images showed enhanced susceptibility difference between blood and the surrounding tissue in both arteries and veins. Our results showed that the proposed Iron Dextran enhanced multi-echo SWI approach allowed the visualization of blood vessels with diameters down to ~100 μm, including small blood vessels supplying and draining small brain structures such as the hippocampus. We conclude that Iron Dextran can be an alternative iron-based MRI contrast agent for blood vessel imaging in the human brain.

Reducing femoral flow artefacts in radial magnetic resonance fingerprinting of the prostate using region-optimised virtual coils

High acceleration factors in radial magnetic resonance fingerprinting (MRF) of the prostate lead to strong streak-like artefacts from flow in the femoral blood vessels, possibly concealing important anatomical information. Region-optimised virtual (ROVir) coils is a beamforming-based framework to create virtual coils that maximise signal in a region of interest while minimising signal in a region of interference. In this study, the potential of removing femoral flow streak artefacts in prostate MRF using ROVir coils is demonstrated in silico and in vivo. The ROVir framework was applied to radial MRF k-space data in an automated pipeline designed to maximise prostate signal while minimising signal from the femoral vessels. The method was tested in 15 asymptomatic volunteers at 3 T. The presence of streaks was visually assessed and measurements of whole prostate T1, T2 and signal-to-noise ratio (SNR) with and without streak correction were examined. In addition, a purpose-built simulation framework in which blood flow through the femoral vessels can be turned on and off was used to quantitatively evaluate ROVir's ability to suppress streaks in radial prostate MRF. In vivo it was shown that removing selected ROVir coils visibly reduces streak-like artefacts from the femoral blood flow, without increasing the reconstruction time. On average, 80% of the prostate SNR was retained. A similar reduction of streaks was also observed in silico, while the quantitative accuracy of T1 and T2 mapping was retained. In conclusion, ROVir coils efficiently suppress streaking artefacts from blood flow in radial MRF of the prostate, thereby improving the visual clarity of the images, without significant sacrifices to acquisition time, reconstruction time and accuracy of quantitative values. This is expected to help enable T1 and T2 mapping of prostate cancer in clinically viable times, aiding differentiation between prostate cancer from noncancer and healthy prostate tissue.

Effect of inhaled oxygen level on dynamic glucose-enhanced MRI in mouse brain

PURPOSE

To investigate the effect of inhaled oxygen level on dynamic glucose enhanced (DGE) MRI in mouse brain tissue and CSF at 3 T.

METHODS

DGE data of brain tissue and CSF from mice under normoxia or hyperoxia were acquired in independent and interleaved experiments using on-resonance variable delay multi-pulse (onVDMP) MRI. A bolus of 0.15 mL filtered 50% D-glucose was injected through the tail vein over 1 min during DGE acquisition. MRS was acquired before and after DGE experiments to confirm the presence of D-glucose.

RESULTS

A significantly higher DGE effect under normoxia than under hyperoxia was observed in brain tissue (p = 0.0001 and p = 0.0002 for independent and interleaved experiments, respectively), but not in CSF (p > 0.3). This difference is attributed to the increased baseline MR tissue signal under hyperoxia induced by a shortened T1 and an increased BOLD effect. When switching from hyperoxia to normoxia without glucose injection, a signal change of ˜3.0% was found in brain tissue and a signal change of ˜1.5% was found in CSF.

CONCLUSIONS

DGE signal was significantly lower under hyperoxia than that under normoxia in brain tissue, but not in CSF. The reason is that DGE effect size of brain tissue is affected by the baseline signal, which could be influenced by T1 change and BOLD effect. Therefore, DGE experiments in which the oxygenation level is changed from baseline need to be interpreted carefully.

Recent Advances in MR Imaging-based Quantification of Brain Oxygen Metabolism

The metabolic rate of oxygen (MRO2) is fundamental to tissue metabolism. Determination of MRO2 demands knowledge of the arterio-venous difference in hemoglobin-bound oxygen concentration, typically expressed as oxygen extraction fraction (OEF), and blood flow rate (BFR). MRI is uniquely suited for measurement of both these quantities, yielding MRO2 in absolute physiologic units of µmol O2 min-1/100 g tissue. Two approaches are discussed, both relying on hemoglobin magnetism. Emphasis will be on cerebral oxygen metabolism expressed in terms of the cerebral MRO2 (CMRO2), but translation of the relevant technologies to other organs, including kidney and placenta will be touched upon as well. The first class of methods exploits the blood's bulk magnetic susceptibility, which can be derived from field maps. The second is based on measurement of blood water T2, which is modulated by diffusion and exchange in the local-induced fields within and surrounding erythrocytes. Some whole-organ methods achieve temporal resolution adequate to permit time-series studies of brain energetics, for instance, during sleep in the scanner with concurrent electroencephalogram (EEG) sleep stage monitoring. Conversely, trading temporal for spatial resolution has led to techniques for spatially resolved approaches based on quantitative blood oxygen level dependent (BOLD) or calibrated BOLD models, allowing regional assessment of vascular-metabolic parameters, both also exploiting deoxyhemoglobin paramagnetism like their whole-organ counterparts.

Unveiling metabolic patterns in dementia: Insights from high-resolution quantitative blood-oxygenation-level-dependent MRI

BACKGROUND

Oxygen extraction fraction (OEF) and deoxyhemoglobin (DoHb) levels reflect variations in cerebral oxygen metabolism in demented patients.

PURPOSE

Delineating the metabolic profiles evident throughout different phases of dementia necessitates an integrated analysis of OEF and DoHb levels. This is enabled by leveraging high-resolution quantitative blood oxygenation level dependent (qBOLD) analysis of magnitude images obtained from a multi-echo gradient-echo MRI (mGRE) scan performed on a 3.0 Tesla scanner.

METHODS

Achieving superior spatial resolution in qBOLD necessitates the utilization of an mGRE scan with only four echoes, which in turn limits the number of measurements compared to the parameters within the qBOLD model. Consequently, it becomes imperative to discard non-essential parameters to facilitate further analysis. This process entails transforming the qBOLD model into a format suitable for fitting the log-magnitude difference (L-MDif) profiles of the four echo magnitudes present in each brain voxel. In order to bolster spatial specificity, the log-difference qBOLD model undergoes refinement into a representative form, termed as r-qBOLD, particularly when applied to class-averaged L-MDif signals derived through k-means clustering of L-MDif signals from all brain voxels into a predetermined number of clusters. The agreement between parameters estimated using r-qBOLD for different cluster sizes is validated using Bland-Altman analysis, and the model's goodness-of-fit is evaluated using aχ 2 ${\chi ^2}$ -test. Retrospective MRI data of Alzheimer's disease (AD), mild cognitive impairment (MCI), and non-demented patients without neuropathological disorders, pacemakers, other implants, or psychiatric disorders, who completed a minimum of three visits prior to MRI enrolment, are utilized for the study.

RESULTS

Utilizing a cohort comprising 30 demented patients aged 65-83 years in stages 4-6 representing mild, moderate, and severe stages according to the clinical dementia rating (CDR), matched with an age-matched non-demented control group of 18 individuals, we conducted joint observations of OEF and DoHb levels estimated using r-qBOLD. The observations elucidate metabolic signatures in dementia based on OEF and DoHb levels in each voxel. Our principal findings highlight the significance of spatial patterns of metabolic profiles (metabolic patterns) within two distinct regimes: OEF levels exceeding the normal range (S1-regime), and OEF levels below the normal range (S2-regime). The S1-regime, accompanied by low DoHb levels, predominantly manifests in fronto-parietal and perivascular regions with increase in dementia severity. Conversely, the S2-regime, accompanied by low DoHb levels, is observed in medial temporal (MTL) regions. Other regions with abnormal metabolic patterns included the orbitofrontal cortex (OFC), medial-orbital prefrontal cortex (MOPFC), hypothalamus, ventro-medial prefrontal cortex (VMPFC), and retrosplenial cortex (RSP). Dysfunction in the OFC and MOPFC indicated cognitive and emotional impairment, while hypothalamic involvement potentially indicated preclinical dementia. Reduced metabolic activity in the RSP suggested early-stage AD related functional abnormalities.

CONCLUSIONS

Integrated analysis of OEF and DoHb levels using r-qBOLD reveals distinct metabolic signatures across dementia phases, highlighting regions susceptible to neuronal loss, vascular involvement, and preclinical indicators.

On the optimization of 3D inflow-based vascular-space-occupancy (iVASO) MRI for the quantification of arterial cerebral blood volume (CBVa)

PURPOSE

The inflow-based vascular-space-occupancy (iVASO) MRI was originally developed in a single-slice mode to measure arterial cerebral blood volume (CBVa). When vascular crushers are applied in iVASO, the signals can be sensitized predominantly to small pial arteries and arterioles. The purpose of this study is to perform a systematic optimization and evaluation of a 3D iVASO sequence on both 3 T and 7 T for the quantification of CBVa values in the human brain.

METHODS

Three sets of experiments were performed in three separate cohorts. (1) 3D iVASO MRI protocols were compared to single-slice iVASO, and the reproducibility of whole-brain 3D iVASO MRI was evaluated. (2) The effects from different vascular crushers in iVASO were assessed. (3) 3D iVASO MRI results were evaluated in arterial and venous blood vessels identified using ultrasmall-superparamagnetic-iron-oxides-enhanced MRI to validate its arterial origin.

RESULTS

3D iVASO scans showed signal-to-noise ratio (SNR) and CBVa measures consistent with single-slice iVASO with reasonable intrasubject reproducibility. Among the iVASO scans performed with different vascular crushers, the whole-brain 3D iVASO scan with a motion-sensitized-driven-equilibrium preparation with two binomial refocusing pulses and an effective TE of 50 ms showed the best suppression of macrovascular signals, with a relatively low specific absorption rate. When no vascular crusher was applied, the CBVa maps from 3D iVASO scans showed large CBVa values in arterial vessels but well-suppressed signals in venous vessels.

CONCLUSION

A whole-brain 3D iVASO MRI scan was optimized for CBVa measurement in the human brain. When only microvascular signals are desired, a motion-sensitized-driven-equilibrium-based vascular crusher with binomial refocusing pulses can be applied in 3D iVASO.

Quantification of whole-organ individual and bilateral renal metabolic rate of oxygen

PURPOSE

Renal metabolic rate of oxygen (rMRO2 ) is a potentially important biomarker of kidney function. The key parameters for rMRO2 quantification include blood flow rate (BFR) and venous oxygen saturation (SvO2 ) in a draining vessel. Previous approaches to quantify renal metabolism have focused on the single organ. Here, both kidneys are considered as one unit to quantify bilateral rMRO2 . A pulse sequence to facilitate bilateral rMRO2 quantification is introduced.

METHODS

To quantify bilateral rMRO2 , measurements of BFR and SvO2 are made along the inferior vena cava (IVC) at suprarenal and infrarenal locations. From the continuity equation, these four parameters can be related to derive an expression for bilateral rMRO2 . The recently reported K-MOTIVE pulse sequence was implemented at four locations: left kidney, right kidney, suprarenal IVC, and infrarenal IVC. A dual-band variant of K-MOTIVE (db-K-MOTIVE) was developed by incorporating simultaneous-multi-slice imaging principles. The sequence simultaneously measures BFR and SvO2 at suprarenal and infrarenal locations in a single pass of 21 s, yielding bilateral rMRO2 .

RESULTS

SvO2 and BFR are higher in suprarenal versus infrarenal IVC, and the renal veins are highly oxygenated (SvO2  >90%). Bilateral rMRO2 quantified in 10 healthy subjects (8 M, 30 ± 8 y) was found to be 291 ± 247 and 349 ± 300 (μmolO2 /min)/100 g, derived from K-MOTIVE and db-K-MOTIVE, respectively. In comparison, total rMRO2 from combining left and right was 329 ± 273 (μmolO2 /min)/100 g.

CONCLUSION

The present work demonstrates that bilateral rMRO2 quantification is feasible with fair reproducibility and physiological plausibility. The indirect method is a promising approach to compute bilateral rMRO2 when individual rMRO2 quantification is difficult.

MRI-based quantification of whole-organ renal metabolic rate of oxygen

During the early stages of diabetes, kidney oxygen utilization increases. The mismatch between oxygen demand and supply contributes to tissue hypoxia, a key driver of chronic kidney disease. Thus, whole-organ renal metabolic rate of oxygen (rMRO2 ) is a potentially valuable biomarker of kidney function. The key parameters required to determine rMRO2 include the renal blood flow rate (RBF) in the feeding artery and oxygen saturation in the draining renal vein (SvO2 ). However, there is currently no noninvasive method to quantify rMRO2 in absolute physiologic units. Here, a new MRI pulse sequence, Kidney Metabolism of Oxygen via T2 and Interleaved Velocity Encoding (K-MOTIVE), is described, along with evaluation of its performance in the human kidney in vivo. K-MOTIVE interleaves a phase-contrast module before a background-suppressed T2 -prepared balanced steady-state-free-precession (bSSFP) readout to measure RBF and SvO2 in a single breath-hold period of 22 s, yielding rMRO2 via Fick's principle. Variants of K-MOTIVE to evaluate alternative bSSFP readout strategies were studied. Kidney mass was manually determined from multislice gradient recalled echo images. Healthy subjects were recruited to quantify rMRO2 of the left kidney at 3-T field strength (N = 15). Assessments of repeat reproducibility and comparisons with individual measurements of RBF and SvO2 were performed, and the method's sensitivity was evaluated with a high-protein meal challenge (N = 8). K-MOTIVE yielded the following metabolic parameters: T2  = 157 ± 19 ms; SvO2  = 92% ± 6%; RBF = 400 ± 110 mL/min; and rMRO2  = 114 ± 117(μmol O2 /min)/100 g tissue. Reproducibility studies of T2 and RBF (parameters directly measured by K-MOTIVE) resulted in coefficients of variation less than 10% and intraclass correlation coefficients more than 0.75. The high-protein meal elicited an increase in rMRO2 , which was corroborated by serum biomarkers. The K-MOTIVE sequence measures SvO2 and RBF, the parameters necessary to quantify whole-organ rMRO2 , in a single breath-hold. The present work demonstrates that rMRO2 quantification is feasible with good reproducibility. rMRO2 is a potentially valuable physiological biomarker.

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

Ion influence on surface water dynamics and proton exchange at protein surfaces - A unified model for transverse and longitudinal NMR relaxation dispersion

In all biologically relevant media, proteins interact in the presence of surrounding ions, and such interactions are water-mediated. Water molecules play a crucial role in the restructuring of proteins in solution and indeed in their biological activity. Surface water dynamics and proton exchange at protein surfaces is investigated here using NMR relaxometry, for two well-known globular proteins, lysozyme and bovine serum albumin, with particular attention to the role of surface ions. We present a unified model of surface water dynamics and proton exchange, accounting simultaneously for the observed longitudinal and transverse relaxation rates. The most notable effect of salt (0.1 M) concerns the slow surface water dynamics, related to rare water molecules embedded in energy wells on the protein surface. This response is protein-specific. On the other hand, the proton exchange time between labile protein-protons and water-protons at the protein surface seems to be very similar for the two proteins and is insensitive to the addition of salts at the concentration studied.

In-cell NMR: Why and how?

NMR spectroscopy has been applied to cells and tissues analysis since its beginnings, as early as 1950. We have attempted to gather here in a didactic fashion the broad diversity of data and ideas that emerged from NMR investigations on living cells. Covering a large proportion of the periodic table, NMR spectroscopy permits scrutiny of a great variety of atomic nuclei in all living organisms non-invasively. It has thus provided quantitative information on cellular atoms and their chemical environment, dynamics, or interactions. We will show that NMR studies have generated valuable knowledge on a vast array of cellular molecules and events, from water, salts, metabolites, cell walls, proteins, nucleic acids, drugs and drug targets, to pH, redox equilibria and chemical reactions. The characterization of such a multitude of objects at the atomic scale has thus shaped our mental representation of cellular life at multiple levels, together with major techniques like mass-spectrometry or microscopies. NMR studies on cells has accompanied the developments of MRI and metabolomics, and various subfields have flourished, coined with appealing names: fluxomics, foodomics, MRI and MRS (i.e. imaging and localized spectroscopy of living tissues, respectively), whole-cell NMR, on-cell ligand-based NMR, systems NMR, cellular structural biology, in-cell NMR… All these have not grown separately, but rather by reinforcing each other like a braided trunk. Hence, we try here to provide an analytical account of a large ensemble of intricately linked approaches, whose integration has been and will be key to their success. We present extensive overviews, firstly on the various types of information provided by NMR in a cellular environment (the "why", oriented towards a broad readership), and secondly on the employed NMR techniques and setups (the "how", where we discuss the past, current and future methods). Each subsection is constructed as a historical anthology, showing how the intrinsic properties of NMR spectroscopy and its developments structured the accessible knowledge on cellular phenomena. Using this systematic approach, we sought i) to make this review accessible to the broadest audience and ii) to highlight some early techniques that may find renewed interest. Finally, we present a brief discussion on what may be potential and desirable developments in the context of integrative studies in biology.

T2 -oximetry-based cerebral venous oxygenation mapping using Fourier-transform-based velocity-selective pulse trains

Purpose

To develop a T₂‐oximetry method for quantitative mapping of cerebral venous oxygenation fraction (Y_v) using Fourier‐transform–based velocity‐selective (FT‐VS) pulse trains.

Methods

The venous isolation preparation was achieved by using an FT‐VS inversion plus a nonselective inversion (NSI) pulse to null the arterial blood signal while minimally affected capillary blood flows out into the venular vasculature during the outflow time (TO), and then applying an Fourier transform based velocity selective saturation (FT‐VSS) pulse to suppress the tissue signal. A multi‐echo readout was employed to obtain venous T₂ (T_2,v) efficiently with the last echo used to detect the residual CSF signal and correct its contamination in the fitting. Here we compared the performance of this FT‐VS–based venous isolation preparations with a traditional velocity‐selective saturation (VSS)–based approach (quantitative imaging of extraction of oxygen and tissue consumption [QUIXOTIC]) with different cutoff velocities for Y_v mapping on 6 healthy volunteers at 3 Tesla.

Results

The FT‐VS–based methods yielded higher venous blood signal and temporal SNR with less CSF contamination than the velocity‐selective saturation–based results. The averaged Y_v values across the whole slice measured in different experiments were close to the global Y_v measured from the individual internal jugular vein.

Conclusion

The feasibility of the FT‐VS–based Y_v estimation was demonstrated on healthy volunteers. The obtained high venous signal as well as the mitigation of CSF contamination led to a good agreement between the T_2,v and Y_v measured in the proposed method with the values in the literature.

Click here for author‐reader discussions

Cerebral oxygen extraction fraction MRI: Techniques and applications

The human brain constitutes 2% of the body's total mass but uses 20% of the oxygen. The rate of the brain's oxygen utilization can be derived from a knowledge of cerebral blood flow and the oxygen extraction fraction (OEF). Therefore, OEF is a key physiological parameter of the brain's function and metabolism. OEF has been suggested to be a useful biomarker in a number of brain diseases. With recent advances in MRI techniques, several MRI-based methods have been developed to measure OEF in the human brain. These MRI OEF techniques are based on the T₂ of blood, the blood signal phase, the magnetic susceptibility of blood-containing voxels, the effect of deoxyhemoglobin on signal behavior in extravascular tissue, and the calibration of the BOLD signal using gas inhalation. Compared to ¹⁵ O PET, which is considered the "gold standard" for OEF measurement, MRI-based techniques are non-invasive, radiation-free, and are more widely available. This article provides a review of these emerging MRI-based OEF techniques. We first briefly introduce the role of OEF in brain oxygen homeostasis. We then review the methodological aspects of different categories of MRI OEF techniques, including their signal mechanisms, acquisition methods, and data analyses. The strengths and limitations of the techniques are discussed. Finally, we review key applications of these techniques in physiological and pathological conditions.

Dependency of R2 and R2 * relaxation on Gd-DTPA concentration in arterial blood: Influence of hematocrit and magnetic field strength

Dynamic susceptibility contrast (DSC) MRI is clinically used to measure brain perfusion by monitoring the dynamic passage of a bolus of contrast agent through the brain. For quantitative analysis of the DSC images, the arterial input function is required. It is known that the original assumption of a linear relation between the R₂^(*) relaxation and the arterial contrast agent concentration is invalid, although the exact relation is as of yet unknown. Studying this relation in vitro is time-consuming, because of the widespread variations in field strengths, MRI sequences, contrast agents, and physiological conditions. This study aims to simulate the R₂^(*) versus contrast concentration relation under varying physiological and technical conditions using an adapted version of an open-source simulation tool. The approach was validated with previously acquired data in human whole blood at 1.5 T by means of a gradient-echo sequence (proof-of-concept). Subsequently, the impact of hematocrit, field strength, and oxygen saturation on this relation was studied for both gradient-echo and spin-echo sequences. The results show that for both gradient-echo and spin-echo sequences, the relaxivity increases with hematocrit and field strength, while the hematocrit dependency was nonlinear for both types of MRI sequences. By contrast, oxygen saturation has only a minor effect. In conclusion, the simulation setup has proven to be an efficient method to rapidly calibrate and estimate the relation between R₂^(*) and gadolinium concentration in whole blood. This knowledge will be useful in future clinical work to more accurately retrieve quantitative information on brain perfusion.

Dynamic contrast-enhanced CEST MRI using a low molecular weight dextran

Natural and synthetic sugars have great potential for developing highly biocompatible and translatable chemical exchange saturation transfer (CEST) MRI contrast agents. In this study, we aimed to develop the smallest clinically available form of dextran, Dex1 (molecular weight, MW ~ 1 kDa), as a new CEST agent. We first characterized the CEST properties of Dex1 in vitro at 11.7 T and showed that the Dex1 had a detectable CEST signal at ~1.2 ppm, attributed to hydroxyl protons. In vivo CEST MRI studies were then carried out on C57BL6 mice bearing orthotopic GL261 brain tumors (n = 5) using a Bruker BioSpec 11.7 T MRI scanner. Both steady-state full Z-spectral images and single offset (1.2 ppm) dynamic dextran-enhanced (DDE) images were acquired before and after the intravenous injection of Dex1 (2 g/kg). The steady-state Z-spectral analysis showed a significantly higher CEST contrast enhancement in the tumor than in contralateral brain (∆MTRasym1.2 ppm  = 0.010 ± 0.006 versus 0.002 ± 0.008, P = 0.0069) at 20 min after the injection of Dex1. Pharmacokinetic analyses of DDE were performed using the area under the curve (AUC) in the first 10 min after Dex1 injection, revealing a significantly higher uptake of Dex1 in the tumor than in brain tissue for tumor-bearing mice (AUC[0-10 min] = 21.9 ± 4.2 versus 5.3 ± 6.4%·min, P = 0.0294). In contrast, no Dex1 uptake was foundling in the brains of non-tumor-bearing mice (AUC[0-10 min] = -1.59 ± 2.43%·min). Importantly, the CEST MRI findings were consistent with the measurements obtained using DCE MRI and fluorescence microscopy, demonstrating the potential of Dex1 as a highly translatable CEST MRI contrast agent for assessing tumor hemodynamics.

Oxygen saturation-dependent effects on blood transverse relaxation at low fields

Objective

Blood oxygenation can be measured using magnetic resonance using the paramagnetic effect of deoxy-haemoglobin, which decreases the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 relaxation time of blood. This \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 contrast has been well characterised at the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{B}_{{0}}$$\end{document}B0 fields used in MRI (1.5 T and above). However, few studies have characterised this effect at lower magnetic fields. Here, the feasibility of blood oximetry at low field based on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 changes that are within a physiological relevant range is explored. This study could be used for specifying requirements for construction of a monitoring device based on low field permanent magnet systems.

Methods

A continuous flow circuit was used to control parameters such as oxygen saturation and temperature in a sample of blood. It flowed through a variable field magnet, where CPMG experiments were performed to measure its \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2. In addition, the oxygen saturation was monitored by an optical sensor for comparison with the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 changes.

Results

These results show that at low \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{B}_{{0}}$$\end{document}B0 fields, the change in blood \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 due to oxygenation is small, but still detectable. The data measured at low fields are also in agreement with theoretical models for the oxy-deoxy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 effect.

Conclusion

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{T}_{2}$$\end{document}T2 changes in blood due to oxygenation were observed at fields as low as 0.1 T. These results suggest that low field NMR relaxometry devices around 0.3 T could be designed to detect changes in blood oxygenation.

Quantitative validation of MRI mapping of cerebral venous oxygenation with direct blood sampling: A graded-O2 study in piglets

PURPOSE

To validate two neonatal cerebral venous oxygenation (Y_v ) MRI techniques, T₂ relaxation under phase contrast (TRUPC) and accelerated TRUPC (aTRUPC) MRI, with oxygenation measured with direct blood sampling.

METHODS

In vivo experiments were performed on seven healthy newborn piglets. For each piglet, a catheter was placed in the superior sagittal sinus to obtain venous blood samples for blood gas oximetry measurement as a gold standard. During the MRI experiment, three to five venous oxygenation levels were achieved in each piglet by varying inhaled O₂ content and breathing rate. Under each condition, Y_v values of the superior sagittal sinus measured by TRUPC, aTRUPC, and blood gas oximetry were obtained. The Y_v quantification in TRUPC and aTRUPC used a standard bovine blood calibration model. The aTRUPC scan was repeated twice to assess its reproducibility. Agreements among TRUPC Y_v , aTRUPC Y_v , and blood gas oximetry were evaluated by intraclass correlation coefficient (ICC) and paired Student's t-test.

RESULTS

The mean hematocrit was 23.6 ± 6.5% among the piglets. Across all measurements, Y_v values were 51.9 ± 21.3%, 54.1 ± 18.8%, and 53.7 ± 19.2% for blood gas oximetry, TRUPC and aTRUPC, respectively, showing no significant difference between any two methods (P > .3). There were good correlations between TRUPC and blood gas Y_v (ICC = 0.801; P < .0001), between aTRUPC and blood gas Y_v (ICC = 0.809; P < .0001), and between aTRUPC and TRUPC Y_v (ICC = 0.887; P < .0001). The coefficient of variation of aTRUPC Y_v was 8.1 ± 9.9%.

CONCLUSION

The values of Y_v measured by TRUPC and aTRUPC were in good agreement with blood gas oximetry. These findings suggest that TRUPC and aTRUPC can provide accurate quantifications of Y_v in major cerebral veins.

Brain metabolism in tau and amyloid mouse models of Alzheimer's disease: An MRI study

Alzheimer's disease (AD) is the leading cause of cognitive impairment and dementia in elderly individuals. According to the current biomarker framework for "unbiased descriptive classification", biomarkers of neurodegeneration, "N", constitute a critical component in the tri-category "A/T/N" system. Current biomarkers of neurodegeneration suffer from potential drawbacks such as requiring invasive lumbar puncture, involving ionizing radiation, or representing a late, irreversible marker. Recent human studies have suggested that reduced brain oxygen metabolism may be a new functional marker of neurodegeneration in AD, but the heterogeneity and the presence of mixed pathology in human patients did not allow a full understanding of the role of oxygen extraction and metabolism in AD. In this report, global brain oxygen metabolism and related physiological parameters were studied in two AD mouse models with relatively pure pathology, using advanced MRI techniques including T2 -relaxation-under-spin-tagging (TRUST) and phase contrast (PC) MRI. Additionally, regional cerebral blood flow (CBF) was determined with pseudocontinuous arterial spin labeling. Reduced global oxygen extraction fraction (by -18.7%, p = 0.008), unit-mass cerebral metabolic rate of oxygen (CMRO2 ) (by -17.4%, p = 0.04) and total CMRO2 (by -30.8%, p < 0.001) were observed in Tau4RΔK mice-referred to as the tau AD model-which manifested pronounced neurodegeneration, as measured by diminished brain volume (by -15.2%, p < 0.001). Global and regional CBF in these mice were not different from those of wild-type mice (p > 0.05), suggesting normal vascular function. By contrast, in B6;SJL-Tg [APPSWE]2576Kha (APP) mice-referred to as the amyloid AD model-no brain volume reduction, as well as relatively intact brain oxygen extraction and metabolism, were found (p > 0.05). Consistent with the imaging data, behavioral measures of walking distance were impaired in Tau4RΔK mice (p = 0.004), but not in APP mice (p = 0.88). Collectively, these findings support the hypothesis that noninvasive MRI measurement of brain oxygen metabolism may be a promising biomarker of neurodegeneration in AD.

Hyperpolarized MRI, functional MRI, MR spectroscopy and CEST to provide metabolic information in vivo

Access to metabolic information in vivo using magnetic resonance (MR) technologies has generally been the niche of MR spectroscopy (MRS) and spectroscopic imaging (MRSI). Metabolic fluxes can be studied using the infusion of substrates labeled with magnetic isotopes, with the use of hyperpolarization especially powerful. Unfortunately, these promising methods are not yet accepted clinically, where fast, simple, and reliable measurement and diagnosis are key. Recent advances in functional MRI and chemical exchange saturation transfer (CEST) MRI allow the use of water imaging to study oxygen metabolism and tissue metabolite levels. These, together with the use of novel data analysis approaches such as machine learning for all of these metabolic MR approaches, are increasing the likelihood of their clinical translation.

Single and Combined Neuroimaging Techniques for Alzheimer's Disease Detection

Alzheimer's disease (AD) consists of the gradual process of decreasing volume and quality of neuron connection in the brain, which consists of gradual synaptic integrity and loss of cognitive functions. In recent years, there has been significant attention in AD classification and early detection with machine learning algorithms. There are different neuroimaging techniques for capturing data and using it for the classification task. Input data as images will help machine learning models to detect different biomarkers for AD classification. This marker has a more critical role for AD detection than other diseases because beta-amyloid can extract complex structures with some metal ions. Most researchers have focused on using 3D and 4D convolutional neural networks for AD classification due to reasonable amounts of data. Also, combination neuroimaging techniques like functional magnetic resonance imaging and positron emission tomography for AD detection have recently gathered much attention. However, gathering a combination of data can be expensive, complex, and tedious. For time consumption reasons, most patients prefer to throw one of the neuroimaging techniques. So, in this review article, we have surveyed different research studies with various neuroimaging techniques and ML methods to see the effect of using combined data as input. The result has shown that the use of the combination method would increase the accuracy of AD detection. Also, according to the sensitivity metrics from different machine learning methods, MRI and fMRI showed promising results.

Three-dimensional whole-brain mapping of cerebral blood volume and venous cerebral blood volume using Fourier transform-based velocity-selective pulse trains

PURPOSE

To develop 3D MRI methods for cerebral blood volume (CBV) and venous cerebral blood volume (vCBV) estimation with whole-brain coverage using Fourier transform-based velocity-selective (FT-VS) pulse trains.

METHODS

For CBV measurement, FT-VS saturation pulse trains were used to suppress static tissue, whereas CSF contamination was corrected voxel-by-voxel using a multi-readout acquisition and a fast CSF T₂ scan. The vCBV mapping was achieved by inserting an arterial-nulling module that included a FT-VS inversion pulse train. Using these methods, CBV and vCBV maps were obtained on 6 healthy volunteers at 3 T.

RESULTS

The mean CBV and vCBV values in gray matter and white matter in different areas of the brain showed high correlation (r = 0.95 and P < .0001). The averaged CBV and vCBV values of the whole brain were 5.4 ± 0.6 mL/100 g and 2.5 ± 0.3 mL/100 g in gray matter, and 2.6 ± 0.5 mL/100 g and 1.5 ± 0.2 mL/100 g in white matter, respectively, comparable to the literature.

CONCLUSION

The feasibility of FT-VS-based CBV and vCBV estimation was demonstrated for 3D acquisition with large spatial coverage.

Intravascular BOLD signal characterization of balanced SSFP experiments in human blood at high to ultrahigh fields

PURPOSE

To investigate the intravascular contribution to the overall balanced SSFP (bSSFP) BOLD effect in human blood at high to ultrahigh field strengths (3 T, 9.4 T, and 14.1 T).

METHODS

Venous blood prepared at two different oxygenation levels (deoxygenated: Y ≈ 71%, oxygenated: Y ≈ 94%) was measured with phase-cycled bSSFP for varying TRs/flip angles at 3 T, 9.4 T, and 14.1 T. The oxygen sensitivity was analyzed by intrinsic MIRACLE (motion-insensitive rapid configuration relaxometry)-R₂ estimation and passband signal differences. The intravascular BOLD-related signal change was extracted from the measured data for microvasculature and macrovasculature, and compared with the extravascular contribution obtained by Monte Carlo simulations.

RESULTS

The MIRACLE-R₂ values showed a characteristic increase with longer TRs in deoxygenated blood, corroborating that SE-R₂ data cannot be used to assess the intravascular bSSFP BOLD effect. Passband bSSFP signal differences measured at optimal flip angles of 30° at 3 T and 20° at 9.4 T/14.1 T revealed considerable relative intravascular contributions of 95%/70% at 3 T, 74%/43% at 9.4 T, 66%/46% at 14.1 T for TR = 5 ms, and 90%/65% at 3 T, 36%/27% at 9.4 T, 13%/15% at 14.1 T for TR = 10 ms in macrovascular/microvascular regimes.

CONCLUSION

The results indicate that intravascular effects have to be considered to better understand the origin of bSSFP BOLD contrast in functional MRI experiments, especially at short TRs. The MIRACLE-R₂ method demonstrated the ability to quantify the apparent decrease in R₂ due to rapid RF refocusing.

Acute-stage MRI cerebral oxygen consumption biomarkers predict 24-hour neurological outcome in a rat cardiac arrest model

Brain injury following cardiac arrest (CA) is thought to be caused by a sudden loss of blood flow resulting in disruption in oxygen delivery, neural function and metabolism. However, temporal trajectories of the brain's physiology in the first few hours following CA have not been fully characterized. Furthermore, the extent to which these early measures can predict future neurological outcomes has not been determined. The present study sought to perform dynamic measurements of cerebral blood flow (CBF), oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO₂ ) with MRI in the first 3 hours following the return of spontaneous circulation (ROSC) in a rat CA model. It was found that CBF, OEF and CMRO₂ all revealed a time-dependent increase during the first 3 hours after the ROSC. Furthermore, the temporal trajectories of CBF and CMRO₂ , but not OEF, were different across rats and related to neurologic outcomes at a later time (24 hours after the ROSC) (P < .001). Rats who manifested better outcomes revealed faster increases in CBF and CMRO₂ during the acute stage. When investigating physiological parameters measured at a single time point, CBF (ρ = 0.82, P = .004) and CMRO₂ (ρ = 0.80, P = .006) measured at ~ 3 hours post-ROSC were positively associated with neurologic outcome scores at 24 hours. These findings shed light on brain physiological changes following CA, and suggest that MRI measures of brain perfusion and metabolism may provide a potential biomarker to guide post-CA management.

Blood-brain barrier permeability measurement by biexponentially modeling whole-brain arterial spin labeling data with multiple T2 -weightings

Blood-brain barrier (BBB) permeability assessment remains of ongoing interest in clinical practice and research. Transitions between intravascular (IV) and extravascular (EV) gray matter (GM) compartments may provide information regarding the microstructural status of the BBB. Due to different transverse relaxation times (T₂ ) of water protons in vessels and GM, it is possible to determine the compartment in which these protons are located. This work presents and investigates the feasibility of a simplified analytical approach for compartmentalizing the proportions of magnetically marked water protons into IV and EV GM components by biexponentially modeling T₂ -weighted arterial spin labeling (ASL) data. Numerous model assumptions were used to stabilize the fit and achieve in vivo applicability. Particularly, transverse relaxation times of IV and EV water protons were determined from the analysis of two supporting T₂ -weighted ASL measurements, utilizing a monoexponential signal model. This stabilized a two-parameter biexponential fit of ASL data with T₂ preparation (PLD = 0.9/1.2/1.5/1.8 s, TE_T2Prep  = 0/30/40/60/80/120/160 ms), which thereby robustly provided estimates of the IV and EV compartment fractions. Experiments were conducted with three healthy volunteers in a 3 T scanner. Averaged over all subjects, the labeled water protons inherit T_2,IV  = 200 ± 18 ms initially and adapt T_2,EV  = 91 ± 2 ms with a longer retention time in cerebral structures. Accordingly, the EVlocated ASL signal fraction rises with increasing PLD from 0.31 ± 0.11 at the shortest PLD of 0.9 s to 0.73 ± 0.02 at the longest PLD of 1.8s. These results indicate a transition of the water protons from IV to EV space. The findings support the potential of biexponential modeling for compartmentalizing ASL spin fractions between IV and EV space. The novel integration of monoexponential parameter estimates stabilizes the two-compartment model fit, suggesting that this technique is suitable for robustly estimating the BBB permeability in vivo.

Oxygen extraction fraction mapping with multi-parametric quantitative BOLD MRI: Reduced transverse relaxation bias using 3D-GraSE imaging

Magnetic resonance imaging (MRI)-based quantification of the blood-oxygenation-level-dependent (BOLD) effect allows oxygen extraction fraction (OEF) mapping. The multi-parametric quantitative BOLD (mq-BOLD) technique facilitates relative OEF (rOEF) measurements with whole brain coverage in clinically applicable scan times. Mq-BOLD requires three separate scans of cerebral blood volume and transverse relaxation rates measured by gradient-echo (1/T₂*) and spin-echo (1/T₂). Although the current method is of clinical merit in patients with stroke, glioma and internal carotid artery stenosis (ICAS), there are relaxation measurement artefacts that impede the sensitivity of mq-BOLD and artificially elevate reported rOEF values.

We posited that T₂-related biases caused by slice refocusing imperfections during rapid 2D-GraSE (Gradient and Spin Echo) imaging can be reduced by applying 3D-GraSE imaging sequences, because the latter requires no slice selective pulses. The removal of T₂-related biases would decrease overestimated rOEF values measured by mq-BOLD. We characterized effects of T₂-related bias in mq-BOLD by comparing the initially employed 2D-GraSE and two proposed 3D-GraSE sequences to multiple single spin-echo reference measurements, both in vitro and in vivo. A phantom and 25 participants, including young and elderly healthy controls as well as ICAS-patients, were scanned. We additionally proposed a procedure to reliably identify and exclude artefact affected voxels. In the phantom, 3D-GraSE derived T₂ values had 57% lower deviation from the reference. For in vivo scans, the formerly overestimated rOEF was reduced by −27% (p < 0.001). We obtained rOEF = 0.51, which is much closer to literature values from positron emission tomography (PET) measurements. Furthermore, increased sensitivity to a focal rOEF elevation in an ICAS-patient was demonstrated.

In summary, the application of 3D-GraSE improves the mq-BOLD-based rOEF quantification while maintaining clinically feasible scan times. Thus, mq-BOLD with non-slice selective T₂ imaging is highly promising to improve clinical diagnostics of cerebrovascular diseases such as ICAS.

Age-Related Alterations in Brain Perfusion, Venous Oxygenation, and Oxygen Metabolic Rate of Mice: A 17-Month Longitudinal MRI Study

Background: Characterization of physiological parameters of the aging brain, such as perfusion and brain metabolism, is important for understanding brain function and diseases. Aging studies on human brain have mostly been based on the cross-sectional design, while the few longitudinal studies used relatively short follow-up time compared to the lifespan. Objectives: To determine the longitudinal time courses of cerebral physiological parameters across the adult lifespan in mice. Methods: The present work examined longitudinal changes in cerebral blood flow (CBF), cerebral venous oxygenation (Yv), and cerebral metabolic rate of oxygen (CMRO2) using MRI in healthy C57BL/6 mice from 3 to 20 months of age. Each mouse received 16 imaging sessions at an ~1-month interval. Results: Significant increases with age were observed in CBF (p = 0.017) and CMRO2 (p < 0.001). Meanwhile, Yv revealed a significant decrease (p = 0.002) with a non-linear pattern (p = 0.013). The rate of change was 0.87, 2.26, and -0.24% per month for CBF, CMRO2, and Yv, respectively. On the other hand, systemic parameters such as heart rate did not show a significant age dependence (p = 0.47). No white-matter-hyperintensities (WMH) were observed on the T2-weighted image at any age of the mice. Conclusion: With age, the mouse brain revealed an increase in oxygen consumption. This observation is consistent with previous findings in humans using a cross-sectional design and suggests a degradation of the brain's energy production or utilization machinery. Cerebral perfusion remains relatively intact in aged mice, at least until 20 months of age, consistent with the absence of WMH in mice.

Morphology of perivascular spaces and enclosed blood vessels in young to middle-aged healthy adults at 7T: Dependences on age, brain region, and breathing gas

Perivascular spaces (PVSs) are fluid-filled spaces surrounding penetrating blood vessels in the brain and are an integral pathway of the glymphatic system. A PVS and the enclosed blood vessel are commonly visualized as a single vessel-like complex (denoted as PVSV) in high-resolution MRI images. Quantitative characterization of the PVSV morphology in MRI images in healthy subjects may serve as a reference for detecting disease related PVS and/or blood vessel alterations in patients with brain diseases.

To this end, we evaluated the age dependences, spatial heterogeneities, and dynamic properties of PVSV morphological features in 45 healthy subjects (21–55 years old), using an ultra-high-resolution three-dimensional transverse relaxation time weighted MRI sequence (0.41 × 0.41 × 0.4 mm³) at 7T. Quantitative PVSV parameters, including apparent diameter, count, volume fraction (VF), and relative contrast to noise ratio (rCNR) were calculated in the white matter and subcortical structures. Dynamic changes were induced by carbogen breathing which are known to induce vasodilation and increase the blood oxygenation level in the brain.

PVSV count and VF significantly increased with age in basal ganglia (BG), so did rCNR in BG, midbrain, and white matter (WM). Apparent PVSV diameter also showed a positive association with age in the three brain regions, although it did not reach statistical significance. The PVSV VF and count showed large inter-subject variations, with coefficients of variation ranging from 0.17 to 0.74 after regressing out age and gender effects. Both apparent diameter and VF exhibited significant spatial heterogeneity, which cannot be explained solely by radio-frequency field inhomogeneities. Carbogen breathing significantly increased VF in BG and WM, and rCNR in thalamus, BG, and WM compared to air breathing.

Our results are consistent with gradual dilation of PVSs with age in healthy adults. The PVSV morphology exhibited spatial heterogeneity and large inter-subject variations and changed during carbogen breathing compared to air breathing.