Mapping vascular network architecture in primate brain using ferumoxytol-weighted laminar MRI

Curation statements for this article:
  • Curated by eLife

    eLife logo

    eLife assessment

    This study presents valuable findings on the relative cerebral blood volume of non-human primates that move us closer to uncovering the functional and architectonic principles that govern the interplay between neuronal and vascular networks. The evidence of areal variations is solid, but that of vessel counting and laminar analysis is incomplete. The lack of a direct comparison of their approach against better-established MRI-based methods for measuring hemodynamics and vascular structure weakens the evidence provided in the current paper version. The work will be of interest to NHP imaging scientists.

This article has been Reviewed by the following groups

Read the full article See related articles

Abstract

Mapping the vascular organization of the brain is of great importance across various domains of basic neuroimaging research, diagnostic radiology, and neurology. However, the intricate task of precisely mapping vasculature across brain regions and cortical layers presents formidable challenges, resulting in a limited understanding of neurometabolic factors influencing the brain’s microvasculature. Addressing this gap, our study investigates whole-brain vascular volume using ferumoxytol-weighted laminar-resolution multi-echo gradient-echo imaging in macaque monkeys. We validate the results with published data for vascular densities and compare them with cytoarchitecture, neuron and synaptic densities. The ferumoxytol-induced change in transverse relaxation rate (ΔR 2 *), an indirect proxy measure of cerebral blood volume (CBV), was mapped onto twelve equivolumetric laminar cortical surfaces. Our findings reveal that CBV varies 3-fold across the brain, with the highest vascular volume observed in the inferior colliculus and lowest in the corpus callosum. In the cerebral cortex, CBV is notably high in early primary sensory areas and low in association areas responsible for higher cognitive functions. Classification of CBV into distinct groups unveils extensive replication of translaminar vascular network motifs, suggesting distinct computational energy supply requirements in areas with varying cytoarchitecture types. Regionally, baseline R 2 * and CBV exhibit positive correlations with neuron density and negative correlations with receptor densities. Adjusting image resolution based on the critical sampling frequency of penetrating cortical vessels allows us to delineate approximately 30% of the arterial-venous vessels. Collectively, these results mark significant methodological and conceptual advancements, contributing to the refinement of cerebrovascular MRI. Furthermore, our study establishes a linkage between neurometabolic factors and the vascular network architecture in the primate brain.

Highlights

Cortical layer vascular mapping using ferumoxytol-weighted R 2 * MRI

Vascular volume is high in primary sensory areas and low in association areas

Correlation between R 2 * and vascular volume with neuron and receptor densities

Vascularization co-varies with densities of specific interneuron types

Article activity feed

  1. eLife assessment

    This study presents valuable findings on the relative cerebral blood volume of non-human primates that move us closer to uncovering the functional and architectonic principles that govern the interplay between neuronal and vascular networks. The evidence of areal variations is solid, but that of vessel counting and laminar analysis is incomplete. The lack of a direct comparison of their approach against better-established MRI-based methods for measuring hemodynamics and vascular structure weakens the evidence provided in the current paper version. The work will be of interest to NHP imaging scientists.

  2. Reviewer #1 (Public review):

    Summary:

    Audio et al. measured cerebral blood volume (CBV) across cortical areas and layers using high-resolution MRI with contrast agents in non-human primates. While the non-invasive CBV MRI methodology is often used to enhance fMRI sensitivity in NHPs, its application for baseline CBV measurement is rare due to the complexities of susceptibility contrast mechanisms. The authors determined the number of large vessels and the areal and laminar variations of CBV in NHP and compared those with various other metrics.

    Strengths:

    Non-invasive mapping of relative cerebral blood volume is novel for non-human primates. A key finding was the observation of variations in CBV across regions; primary sensory cortices had high CBV, whereas other higher areas had low CBV. The measured CBV values correlated with previously reported neuronal and receptor densities.

    Weaknesses:

    A weakness of this manuscript is that the quantification of CBV with postprocessing approaches to remove susceptibility effects from pial and penetrating vessels, as well as orientation dependency, is not fully validated, especially on a laminar scale. Further specific comments follow.

    (1) Baseline CBV indices were determined using contrast agent-enhanced MRI (deltaR2*). Although this approach is suitable for areal comparisons, its application on a laminar scale has not been validated in the literature or in this study. By comparing with histological vascular information of V1, the authors attempted to validate their approach. However, the generalization of their method is questionable. The main issue is whether the large vessel contribution is minimized by processing approaches properly in various cortical areas (such as clusters 1-3 in Figure 5). It would be beneficial to compare deltaR2* with deltaR2 induced by contrast agents in a few selected slices, as deltaR2 is supposed to be sensitive to microvessels, not macrovessels. Please discuss this issue.

    (2) High-resolution MRI with a critical sampling frequency estimated from previous studies (Weber 2008, Zheng 1991) was performed to separate penetrating vessels, which is considered one of the major advancements in this study. However, this approach is still insufficient to accurately identify the number of vessels due to the blooming effects of susceptibility and insufficient spatial resolution. There was no detailed description of the detection criteria. More importantly, the number of observable penetrating vessels is dependent on imaging parameters and the dose of the contrast agent. If imaging slices were obtained in parallel to the cortex with higher in-plane resolution, it would likely improve the detection of penetrating vessels. Using higher-field MRI would further enhance the detection of penetrating vessels. Therefore, the reported value is only applicable to the experimental and processing conditions used in this study. Detailed selection criteria should be mentioned, and all potential pitfalls should be discussed.

    (3) Attempts to obtain pial vascular structures were made (Figure 2). As mentioned in this manuscript, the blooming effect of susceptibility contrasts is problematic. In the MRI community, T1-based Gd contrast agents have been used for mapping large vasculature, which is a better approach for obtaining pial vascular structures. Alternatively, computer tomography with a blood contrast agent can be used for mapping blood vasculature noninvasively. This issue should be discussed.

    (4) Since baseline R2* is related to baseline R2, vascular volume, iron content, and susceptibility gradients, it is difficult to correlate it with physiological parameters. Baseline R2* is also sensitive to imaging parameters; higher spatial resolution tends to result in lower R2* values (closer to the R2 value). Therefore, baseline R2* findings need to be emphasized.

    (5) CBV-weighted deltaR2* is correlated with various other metrics (cytoarchitectural parcellation, myelin/receptor density, cortical thickness, CO, cell-type specificity, etc.). While testing the correlation between deltaR2* and these other metrics may be acceptable as an exploratory analysis, it is challenging for readers to discern a causal relationship between them. A critical question is whether CBV-weighted deltaR2* can provide insights into other metrics in diseased or abnormal brain states. If this is the case, then high-resolution deltaR2* will be useful. Please comment on this possibility.

    (6) There is no discussion about the deltaR2* difference across subcortical areas (Figure 1). This finding is intriguing and warrants a thorough discussion in the context of the cortical findings.

    (7) Figure 3 is missing. Several statements in the manuscript require statistics (e.g., bimodality in Figure 2D, Figure 3F).

  3. Reviewer #2 (Public review):

    Summary:

    This manuscript presents a new approach for non-invasive, MRI-based measurements of cerebral blood volume (CBV). Here, the authors use ferumoxytol, a high-contrast agent, and apply specific sequences to infer CBV. The authors then move to statistically compare measured regional CBV with the known distribution of different types of neurons, markers of metabolic load, and others. While the presented methodology captures an estimated 30% of the vasculature, the authors corroborated previous findings regarding the lack of vascular compartmentalization around functional neuronal units in the primary visual cortex.

    Strengths:

    Non-invasive methodology geared to map vascular properties in vivo.

    Implementation of a highly sensitive approach for measuring blood volume.

    Ability to map vascular structural and functional vascular metrics to other types of published data.

    Weaknesses:

    The key issue here is the underlying assumption about the appropriate spatial sampling frequency needed to capture the architecture of the brain vasculature. Namely, ~7 penetrating vessels / mm2 as derived from Weber et al 2008 (Cer Cor). The cited work begins by characterizing the spacing of penetrating arteries and ascending veins using a vascular cast of 7 monkeys (Macaca mulatta, same as in the current paper). The ~7 penetrating vessels / mm2 are computed by dividing the total number of identified vessels by the area imaged. The problem here is that all measurements were made in a "non-volumetric" manner and only in V1. Extrapolating from here to the entire brain seems like an over-assumption, particularly given the region-dependent heterogeneity that the current paper reports.