Dominance of layer-specific microvessel dilation in contrast-enhanced high-resolution fMRI: Comparison between hemodynamic spread and vascular architecture with CLARITY

Spatiotemporal Clustering of Functional Ultrasound Signals at the Single-Voxel Level

Functional ultrasound (fUS) imaging is a well-established neuroimaging technology that offers high spatiotemporal resolution and a large field of view. Typical strategies for analyzing fUS data comprise either region-based averaging, typically based on reference atlases, or correlation with experimental events. Nevertheless, these methodologies possess several inherent limitations, including a restricted utilization of the spatial dimension and a pronounced bias influenced by preconceived notions about the recorded activity. In this study, we put forth single-voxel clustering as a third method to address these issues. A comparison was conducted between the three strategies on a typical dataset comprising visually evoked activity in the superior colliculus in awake mice. The application of single-voxel clustering yielded the generation of detailed activity maps, which revealed a consistent layout of activity and a clear separation between hemodynamic responses. This method is best considered as a complement to region-based averaging and correlation. It has direct applicability to challenging contexts, such as paradigm-free analysis on behaving subjects and brain decoding.

Mesoscale Brain Mapping: Bridging Scales and Modalities in Neuroimaging - A Symposium Review

Advances in the spatiotemporal resolution and field-of-view of neuroimaging tools are driving mesoscale studies for translational neuroscience. On October 10, 2023, the Center for Mesoscale Mapping (CMM) at the Massachusetts General Hospital (MGH) Athinoula A. Martinos Center for Biomedical Imaging and the Massachusetts Institute of Technology (MIT) Health Sciences Technology based Neuroimaging Training Program (NTP) hosted a symposium exploring the state-of-the-art in this rapidly growing area of research. "Mesoscale Brain Mapping: Bridging Scales and Modalities in Neuroimaging" brought together researchers who use a broad range of imaging techniques to study brain structure and function at the convergence of the microscopic and macroscopic scales. The day-long event centered on areas in which the CMM has established expertise, including the development of emerging technologies and their application to clinical translational needs and basic neuroscience questions. The in-person symposium welcomed more than 150 attendees, including 57 faculty members, 61 postdoctoral fellows, 35 students, and four industry professionals, who represented institutions at the local, regional, and international levels. The symposium also served the training goals of both the CMM and the NTP. The event content, organization, and format were planned collaboratively by the faculty and trainees. Many CMM faculty presented or participated in a panel discussion, thus contributing to the dissemination of both the technologies they have developed under the auspices of the CMM and the findings they have obtained using those technologies. NTP trainees who benefited from the symposium included those who helped to organize the symposium and/or presented posters and gave "flash" oral presentations. In addition to gaining experience from presenting their work, they had opportunities throughout the day to engage in one-on-one discussions with visiting scientists and other faculty, potentially opening the door to future collaborations. The symposium presentations provided a deep exploration of the many technological advances enabling progress in structural and functional mesoscale brain imaging. Finally, students worked closely with the presenting faculty to develop this report summarizing the content of the symposium and putting it in the broader context of the current state of the field to share with the scientific community. We note that the references cited here include conference abstracts corresponding to the symposium poster presentations.

Functional ultrasound imaging and neuronal activity: how accurate is the spatiotemporal match?

Over the last decade, functional ultrasound (fUS) has risen as a critical tool in functional neuroimaging, leveraging hemodynamic changes to infer neural activity indirectly. Recent studies have established a strong correlation between neural spike rates (SR) and functional ultrasound signals. However, understanding their spatial distribution and variability across different brain areas is required to thoroughly interpret fUS signals. In this regard, we conducted simultaneous fUS imaging and Neuropixels recordings during stimulus-evoked activity in awake mice within three regions the visual pathway. Our findings indicate that the temporal dynamics of fUS and SR signals are linearly correlated, though the correlation coefficients vary among visual regions. Conversely, the spatial correlation between the two signals remains consistent across all regions with a spread of approximately 300 micrometers. Finally, we introduce a model that integrates the spatial and temporal components of the fUS signal, allowing for a more accurate interpretation of fUS images.

Odor-evoked layer-specific fMRI activities in the awake mouse olfactory bulb

Awake rodent fMRI is increasingly common over the use of anesthesia since it permits behavioral paradigms and does not confound normal brain function or neurovascular coupling. It is well established that adequate acclimation to the loud fMRI environment and head fixation reduces stress in the rodents and allows for whole brain imaging with little contamination from motion. However, it is unknown whether high-resolution fMRI with increased susceptibility to motion and lower sensitivity can measure small, but spatially discrete, activations in awake mice. To examine this, we used contrast-enhanced cerebral blood volume-weighted (CBVw) fMRI in the mouse olfactory bulb for its enhanced sensitivity and neural specificity. We determined that activation patterns in the glomerular layer to four different odors were spatially distinct and were consistent with previously established histological patterns. In addition, odor-evoked laminar activations were greatest in superficial layers that decreased with laminar depth, similar to previous observations. Interestingly, the fMRI response strengths in the granule cell layer were greater in awake mice than our previous anesthetized rat studies, suggesting that feedback neural activities were intact with wakefulness. We finally determined that fMRI signal changes to repeated odor exposure (i.e., olfactory adaptation) attenuated relatively more in the feedback granule cell layer compared to the input glomerular layer, which is consistent with prior observations. We, therefore, conclude that high-resolution CBVw fMRI can measure odor-specific activation patterns and distinguish changes in laminar activity of head and body restrained awake mice.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Dissection of brain-wide resting-state and functional somatosensory circuits by fMRI with optogenetic silencing

To further advance functional MRI (fMRI)-based brain science, it is critical to dissect fMRI activity at the circuit level. To achieve this goal, we combined brain-wide fMRI with neuronal silencing in well-defined regions. Since focal inactivation suppresses excitatory output to downstream pathways, intact input and suppressed output circuits can be separated. Highly specific cerebral blood volume-weighted fMRI was performed with optogenetic stimulation of local GABAergic neurons in mouse somatosensory regions. Brain-wide spontaneous somatosensory networks were found mostly in ipsilateral cortical and subcortical areas, which differed from the bilateral homotopic connections commonly observed in resting-state fMRI data. The evoked fMRI responses to somatosensory stimulation in regions of the somatosensory network were successfully dissected, allowing the relative contributions of spinothalamic (ST), thalamocortical (TC), corticothalamic (CT), corticocortical (CC) inputs, and local intracortical circuits to be determined. The ventral posterior thalamic nucleus receives ST inputs, while the posterior medial thalamic nucleus receives CT inputs from the primary somatosensory cortex (S1) with TC inputs. The secondary somatosensory cortex (S2) receives mostly direct CC inputs from S1 and a few TC inputs from the ventral posterolateral nucleus. The TC and CC input layers in cortical regions were identified by laminar-specific fMRI responses with a full width at half maximum of <150 µm. Long-range synaptic inputs in cortical areas were amplified approximately twofold by local intracortical circuits, which is consistent with electrophysiological recordings. Overall, whole-brain fMRI with optogenetic inactivation revealed brain-wide, population-based, long-range circuits, which could complement data typically collected in conventional microscopic functional circuit studies.

Layer-dependent functional connectivity methods

Recent methodological advances in fMRI contrast and readout strategies have allowed researchers to approach the mesoscopic spatial regime of cortical layers. This has revolutionized the ability to map cortical information processing within and across brain systems. However, until recently, most layer-fMRI studies have been confined to primary cortices using basic block-design tasks and macro-vascular-contaminated sequence contrasts. To become an established method for user-friendly applicability in neuroscience practice, layer-fMRI acquisition and analysis methods need to be extended to more flexible connectivity-based experiment designs; they must be able to capture subtle changes in brain networks of higher-order cognitive areas, and they should not be spatially biased with unwanted vein signals. In this article, we review the most pressing challenges of layer-dependent fMRI for large-scale neuroscientific applicability and describe recently developed acquisition methodologies that can resolve them. In doing so, we review technical tradeoffs and capabilities of modern MR-sequence approaches to achieve measurements that are free of locally unspecific vein signal, with whole-brain coverage, sub-second sampling, high resolutions, and with a combination of those capabilities. The presented approaches provide whole-brain layer-dependent connectivity data that open a new window to investigate brain network connections. We exemplify this by reviewing a number of candidate tools for connectivity analyses that will allow future studies to address new questions in network neuroscience. The considered network analysis tools include: hierarchy mapping, directional connectomics, source-specific connectivity mapping, and network sub-compartmentalization. We conclude: Whole-brain layer-fMRI without large-vessel contamination is applicable for human neuroscience and opens the door to investigate biological mechanisms behind any number of psychological and psychiatric phenomena, such as selective attention, hallucinations and delusions, and even conscious perception.

Huntingtin silencing delays onset and slows progression of Huntington's disease: a biomarker study

Huntington's disease is a dominantly inherited, fatal neurodegenerative disorder caused by a CAG expansion in the huntingtin (HTT) gene, coding for pathological mutant HTT protein (mHTT). Because of its gain-of-function mechanism and monogenic aetiology, strategies to lower HTT are being actively investigated as disease-modifying therapies. Most approaches are currently targeted at the manifest stage, where clinical outcomes are used to evaluate the effectiveness of therapy. However, as almost 50% of striatal volume has been lost at the time of onset of clinical manifest, it would be preferable to begin therapy in the premanifest period. An unmet challenge is how to evaluate therapeutic efficacy before the presence of clinical symptoms as outcome measures. To address this, we aim to develop non-invasive sensitive biomarkers that provide insight into therapeutic efficacy in the premanifest stage of Huntington's disease. In this study, we mapped the temporal trajectories of arteriolar cerebral blood volumes (CBVa) using inflow-based vascular-space-occupancy (iVASO) MRI in the heterozygous zQ175 mice, a full-length mHTT expressing and slowly progressing model with a premanifest period as in human Huntington's disease. Significantly elevated CBVa was evident in premanifest zQ175 mice prior to motor deficits and striatal atrophy, recapitulating altered CBVa in human premanifest Huntington's disease. CRISPR/Cas9-mediated non-allele-specific HTT silencing in striatal neurons restored altered CBVa in premanifest zQ175 mice, delayed onset of striatal atrophy, and slowed the progression of motor phenotype and brain pathology. This study-for the first time-shows that a non-invasive functional MRI measure detects therapeutic efficacy in the premanifest stage and demonstrates long-term benefits of a non-allele-selective HTT silencing treatment introduced in the premanifest Huntington's disease.

Tortuosity Index Calculations in Retinal Images: Some Criticalities Arising from Commonly Used Approaches

A growing body of research in retinal imaging is recently considering vascular tortuosity measures or indexes, with definitions and methods mostly derived from cardiovascular research. However, retinal microvasculature has its own peculiarities that must be considered in order to produce reliable measurements. This study analyzed and compared various derived metrics (e.g., TI, TI_avg, TI*CV) across four existing computational workflows. Specifically, the implementation of the models on two critical OCT images highlighted main pitfalls of the methods, which may fail in reliably differentiating a highly tortuous image from a normal one. A tentative, encouraging approach to mitigate the issue on the same OCT exemplificative images is described in the paper, based on the suggested index TI*CV.

Postsynaptic activity of inhibitory neurons evokes hemodynamic fMRI responses

Functional MRI responses are localized to the synaptic sites of evoked inhibitory neurons, but it is unknown whether, or by what mechanisms, these neurons initiate functional hyperemia. Here, the neuronal origins of these hemodynamic responses were investigated by fMRI or local field potential and blood flow measurements during topical application of pharmacological agents when GABAergic granule cells in the rat olfactory bulb were synaptically targeted. First, to examine if postsynaptic activation of these inhibitory neurons was required for neurovascular coupling, we applied an NMDA receptor antagonist during cerebral blood volume-weighted fMRI acquisition and found that responses below the drug application site (up to ~1.5 mm) significantly decreased within ~30 min. Similarly, large decreases in granule cell postsynaptic activities and blood flow responses were observed when AMPA or NMDA receptor antagonists were applied. Second, inhibition of nitric oxide synthase preferentially decreased the initial, fast component of the blood flow response, while inhibitors of astrocyte-specific glutamate transporters and vasoactive intestinal peptide receptors did not decrease blood flow responses. Third, inhibition of GABA release with a presynaptic GABAB receptor agonist caused less reduction of neuronal and blood flow responses compared to the postsynaptic glutamate receptor antagonists. In conclusion, local hyperemia by synaptically-evoked inhibitory neurons was primarily driven by their postsynaptic activities, possibly through NMDA receptor-dependent calcium signaling that was not wholly dependent on nitric oxide.

Time-dependent spatial specificity of high-resolution fMRI: insights into mesoscopic neurovascular coupling

High-resolution functional magnetic resonance imaging (fMRI) is becoming increasingly popular because of the growing availability of ultra-high magnetic fields which are capable of improving sensitivity and spatial resolution. However, it is debatable whether increased spatial resolutions for haemodynamic-based techniques, like fMRI, can accurately detect the true location of neuronal activity. We have addressed this issue in functional columns and layers of animals with haemoglobin-based optical imaging and different fMRI contrasts, such as blood oxygenation level-dependent, cerebral blood flow and cerebral blood volume fMRI. In this review, we describe empirical evidence primarily from our own studies on how well these fMRI signals are spatially specific to the neuronally active site and discuss insights into neurovascular coupling at the mesoscale. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Ultra-slow Oscillations in fMRI and Resting-State Connectivity: Neuronal and Vascular Contributions and Technical Confounds

Ultra-slow, ∼0.1-Hz variations in the oxygenation level of brain blood are widely used as an fMRI-based surrogate of "resting-state" neuronal activity. The temporal correlations among these fluctuations across the brain are interpreted as "functional connections" for maps and neurological diagnostics. Ultra-slow variations in oxygenation follow a cascade. First, they closely track changes in arteriole diameter. Second, interpretable functional connections arise when the ultra-slow changes in amplitude of γ-band neuronal oscillations, which are shared across even far-flung but synaptically connected brain regions, entrain the ∼0.1-Hz vasomotor oscillation in diameter of local arterioles. Significant confounds to estimates of functional connectivity arise from residual vasomotor activity as well as arteriole dynamics driven by self-generated movements and subcortical common modulatory inputs. Last, methodological limitations of fMRI can lead to spurious functional connections. The neuronal generator of ultra-slow variations in γ-band amplitude, including that associated with self-generated movements, remains an open issue.

Mapping optogenetically-driven single-vessel fMRI with concurrent neuronal calcium recordings in the rat hippocampus

Extensive in vivo imaging studies investigate the hippocampal neural network function, mainly focusing on the dorsal CA1 region given its optical accessibility. Multi-modality fMRI with simultaneous hippocampal electrophysiological recording reveal broad cortical correlation patterns, but the detailed spatial hippocampal functional map remains lacking given the limited fMRI resolution. In particular, hemodynamic responses linked to specific neural activity are unclear at the single-vessel level across hippocampal vasculature, which hinders the deciphering of the hippocampal malfunction in animal models and the translation to critical neurovascular coupling (NVC) patterns for human fMRI. We simultaneously acquired optogenetically-driven neuronal Ca2+ signals with single-vessel blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-fMRI from individual venules and arterioles. Distinct spatiotemporal patterns of hippocampal hemodynamic responses were correlated to optogenetically evoked and spreading depression-like calcium events. The calcium event-related single-vessel hemodynamic modeling revealed significantly reduced NVC efficiency upon spreading depression-like (SDL) events, providing a direct measure of the NVC function at various hippocampal states.

Magnetic Resonance Imaging technology-bridging the gap between noninvasive human imaging and optical microscopy

Technological advances in Magnetic Resonance Imaging (MRI) have provided substantial gains in the sensitivity and specificity of functional neuroimaging. Mounting evidence demonstrates that the hemodynamic changes utilized in functional MRI can be far more spatially and thus neuronally specific than previously believed. This has motivated a push toward novel, high-resolution MR imaging strategies that can match this biological resolution limit while recording from the entire human brain. Although sensitivity increases are a necessary component, new MR encoding technologies are required to convert improved sensitivity into higher resolution. These new sampling strategies improve image acquisition efficiency and enable increased image encoding in the time-frame needed to follow hemodynamic changes associated with brain activation.

Ultra-Slow Single-Vessel BOLD and CBV-Based fMRI Spatiotemporal Dynamics and Their Correlation with Neuronal Intracellular Calcium Signals

Functional MRI has been used to map brain activity and functional connectivity based on the strength and temporal coherence of neurovascular-coupled hemodynamic signals. Here, single-vessel fMRI reveals vessel-specific correlation patterns in both rodents and humans. In anesthetized rats, fluctuations in the vessel-specific fMRI signal are correlated with the intracellular calcium signal measured in neighboring neurons. Further, the blood-oxygen-level-dependent (BOLD) signal from individual venules and the cerebral-blood-volume signal from individual arterioles show correlations at ultra-slow (<0.1 Hz), anesthetic-modulated rhythms. These data support a model that links neuronal activity to intrinsic oscillations in the cerebral vasculature, with a spatial correlation length of ∼2 mm for arterioles. In complementary data from awake human subjects, the BOLD signal is spatially correlated among sulcus veins and specified intracortical veins of the visual cortex at similar ultra-slow rhythms. These data support the use of fMRI to resolve functional connectivity at the level of single vessels.