Simultaneous cortical, subcortical, and brainstem mapping of sensory activation.

Nonpainful tactile sensory stimuli are processed in the cortex, subcortex, and brainstem. Recent functional magnetic resonance imaging studies have highlighted the value of whole-brain, systems-level investigation for examining sensory processing. However, whole-brain functional magnetic resonance imaging studies are uncommon, in part due to challenges with signal to noise when studying the brainstem. Furthermore, differentiation of small sensory brainstem structures such as the cuneate and gracile nuclei necessitates high-resolution imaging. To address this gap in systems-level sensory investigation, we employed a whole-brain, multi-echo functional magnetic resonance imaging acquisition at 3T with multi-echo independent component analysis denoising and brainstem-specific modeling to enable detection of activation across the entire sensory system. In healthy participants, we examined patterns of activity in response to nonpainful brushing of the right hand, left hand, and right foot (n = 10 per location), and found the expected lateralization, with distinct cortical and subcortical responses for upper and lower limb stimulation. At the brainstem level, we differentiated the adjacent cuneate and gracile nuclei, corresponding to hand and foot stimulation respectively. Our findings demonstrate that simultaneous cortical, subcortical, and brainstem mapping at 3T could be a key tool to understand the sensory system in both healthy individuals and clinical cohorts with sensory deficits.

Toward precision brain health: accurate prediction of a cognitive index trajectory using neuroimaging metrics.

The goal of precision brain health is to accurately predict individuals' longitudinal patterns of brain change. We trained a machine learning model to predict changes in a cognitive index of brain health from neurophysiologic metrics. A total of 48 participants (ages 21-65) completed a sensorimotor task during 2 functional magnetic resonance imaging sessions 6 mo apart. Hemodynamic response functions (HRFs) were parameterized using traditional (amplitude, dispersion, latency) and novel (curvature, canonicality) metrics, serving as inputs to a neural network model that predicted gain on indices of brain health (cognitive factor scores) for each participant. The optimal neural network model successfully predicted substantial gain on the cognitive index of brain health with 90% accuracy (determined by 5-fold cross-validation) from 3 HRF parameters: amplitude change, dispersion change, and similarity to a canonical HRF shape at baseline. For individuals with canonical baseline HRFs, substantial gain in the index is overwhelmingly predicted by decreases in HRF amplitude. For individuals with non-canonical baseline HRFs, substantial gain in the index is predicted by congruent changes in both HRF amplitude and dispersion. Our results illustrate that neuroimaging measures can track cognitive indices in healthy states, and that machine learning approaches using novel metrics take important steps toward precision brain health.

fMRI, LFP, and anatomical evidence for hierarchical nociceptive routing pathway between somatosensory and insular cortices.

The directional organization of multiple nociceptive regions, particularly within obscure operculoinsular areas, underlying multidimensional pain processing remains elusive. This study aims to establish the fundamental organization between somatosensory and insular cortices in routing nociceptive information. By employing an integrated multimodal approach of high-field fMRI, intracranial electrophysiology, and transsynaptic viral tracing in rats, we observed a hierarchically organized connection of S1/S2 → posterior insula → anterior insula in routing nociceptive information. The directional nociceptive pathway determined by early fMRI responses was consistent with that examined by early evoked LFP, intrinsic effective connectivity, and anatomical projection, suggesting fMRI could provide a valuable facility to discern directional neural circuits in animals and humans non-invasively. Moreover, our knowledge of the nociceptive hierarchical organization of somatosensory and insular cortices and the interface role of the posterior insula may have implications for the development of targeted pain therapies.

Negative association between neurovascular coupling and cortical gray matter volume during the lifespan.

Recent studies have established the moment-to-moment turnover of the blood-oxygen-level-dependent signal (TBOLD) at resting state as a key measure of local cortical brain function. Here, we sought to extend that line of research by evaluating TBOLD in 70 cortical areas with respect to corresponding brain volume, age, and sex across the lifespan in 1,344 healthy participants including 633 from the Human Connectome Project (HCP)-Development cohort (294 males and 339 females, age range 8-21 yr) and 711 healthy participants from HCP-Aging cohort (316 males and 395 females, 36-90 yr old). In both groups, we found that 1) TBOLD increased with age, 2) volume decreased with age, and 3) TBOLD and volume were highly significantly negatively correlated, independent of age. The inverse association between TBOLD and volume was documented in nearly all 70 brain areas and for both sexes, with slightly stronger associations documented for males. The strong correspondence between TBOLD and volume across age and sex suggests a common influence such as chronic neuroinflammation contributing to reduced cortical volume and increased TBOLD across the lifespan.NEW & NOTEWORTHY We report a significant negative association between resting functional magnetic resonance imaging (fMRI) blood-oxygen-level-dependent (BOLD) signal turnover (TBOLD) and cortical gray matter volume across the lifespan, such that TBOLD increased whereas volume decreased. We attribute this association to a hypothesized chronic, low-grade neuroinflammation, probably induced by various neurotropic pathogens, including human herpes viruses known to be dormant in the brain in a latent state and reactivated by stress, fever, and various environmental exposures, such as ultraviolet light.

Reliability of variability and complexity measures for task and task-free BOLD fMRI.

Brain activity continuously fluctuates over time, even if the brain is in controlled (e.g., experimentally induced) states. Recent years have seen an increasing interest in understanding the complexity of these temporal variations, for example with respect to developmental changes in brain function or between-person differences in healthy and clinical populations. However, the psychometric reliability of brain signal variability and complexity measures-which is an important precondition for robust individual differences as well as longitudinal research-is not yet sufficiently studied. We examined reliability (split-half correlations) and test-retest correlations for task-free (resting-state) BOLD fMRI as well as split-half correlations for seven functional task data sets from the Human Connectome Project to evaluate their reliability. We observed good to excellent split-half reliability for temporal variability measures derived from rest and task fMRI activation time series (standard deviation, mean absolute successive difference, mean squared successive difference), and moderate test-retest correlations for the same variability measures under rest conditions. Brain signal complexity estimates (several entropy and dimensionality measures) showed moderate to good reliabilities under both, rest and task activation conditions. We calculated the same measures also for time-resolved (dynamic) functional connectivity time series and observed moderate to good reliabilities for variability measures, but poor reliabilities for complexity measures derived from functional connectivity time series. Global (i.e., mean across cortical regions) measures tended to show higher reliability than region-specific variability or complexity estimates. Larger subcortical regions showed similar reliability as cortical regions, but small regions showed lower reliability, especially for complexity measures. Lastly, we also show that reliability scores are only minorly dependent on differences in scan length and replicate our results across different parcellation and denoising strategies. These results suggest that the variability and complexity of BOLD activation time series are robust measures well-suited for individual differences research. Temporal variability of global functional connectivity over time provides an important novel approach to robustly quantifying the dynamics of brain function. PRACTITIONER POINTS: Variability and complexity measures of BOLD activation show good split-half reliability and moderate test-retest reliability. Measures of variability of global functional connectivity over time can robustly quantify neural dynamics. Length of fMRI data has only a minor effect on reliability.

Real-time quantification of in vivo cerebrospinal fluid velocity using the functional magnetic resonance imaging inflow effect.

In vivo estimation of cerebrospinal fluid (CSF) velocity is crucial for understanding the glymphatic system and its potential role in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease. Current cardiac or respiratory-gated approaches, such as 4D flow magnetic resonance imaging (MRI), cannot capture CSF movement in real time because of limited temporal resolution and, in addition, deteriorate in accuracy at low fluid velocities. Other techniques like real-time phase-contrast-MRI or time-spatial labeling inversion pulse are not limited by temporal averaging but have limited availability, even in research settings. This study aims to quantify the inflow effect of dynamic CSF motion on functional MRI (fMRI) for in vivo, real-time measurement of CSF flow velocity. We considered linear and nonlinear models of velocity waveforms and empirically fit them to fMRI data from a controlled flow experiment. To assess the utility of this methodology in human data, CSF flow velocities were computed from fMRI data acquired in eight healthy volunteers. Breath-holding regimens were used to amplify CSF flow oscillations. Our experimental flow study revealed that CSF velocity is nonlinearly related to inflow effect-mediated signal increase and well estimated using an extension of a previous nonlinear framework. Using this relationship, we recovered velocity from in vivo fMRI signal, demonstrating the potential of our approach for estimating CSF flow velocity in the human brain. This novel method could serve as an alternative approach to quantifying slow flow velocities in real time, such as CSF flow in the ventricular system, thereby providing valuable insights into the glymphatic system's function and its implications for neurological disorders.

Masked features of task states found in individual brain networks.

Completing complex tasks requires that we flexibly integrate information across brain areas. While studies have shown how functional networks are altered during different tasks, this work has generally focused on a cross-subject approach, emphasizing features that are common across people. Here we used extended sampling "precision" fMRI data to test the extent to which task states generalize across people or are individually specific. We trained classifiers to decode state using functional network data in single-person datasets across 5 diverse task states. Classifiers were then tested on either independent data from the same person or new individuals. Individualized classifiers were able to generalize to new participants. However, classification performance was significantly higher within a person, a pattern consistent across model types, people, tasks, feature subsets, and even for decoding very similar task conditions. Notably, these findings also replicated in a new independent dataset. These results suggest that individual-focused approaches can uncover robust features of brain states, including features obscured in cross-subject analyses. Individual-focused approaches have the potential to deepen our understanding of brain interactions during complex cognition.

Multivariate fMRI responses in superior temporal cortex predict visual contributions to, and individual differences in, the intelligibility of noisy speech.

Humans have the unique ability to decode the rapid stream of language elements that constitute speech, even when it is contaminated by noise. Two reliable observations about noisy speech perception are that seeing the face of the talker improves intelligibility and the existence of individual differences in the ability to perceive noisy speech. We introduce a multivariate BOLD fMRI measure that explains both observations. In two independent fMRI studies, clear and noisy speech was presented in visual, auditory and audiovisual formats to thirty-seven participants who rated intelligibility. An event-related design was used to sort noisy speech trials by their intelligibility. Individual-differences multidimensional scaling was applied to fMRI response patterns in superior temporal cortex and the dissimilarity between responses to clear speech and noisy (but intelligible) speech was measured. Neural dissimilarity was less for audiovisual speech than auditory-only speech, corresponding to the greater intelligibility of noisy audiovisual speech. Dissimilarity was less in participants with better noisy speech perception, corresponding to individual differences. These relationships held for both single word and entire sentence stimuli, suggesting that they were driven by intelligibility rather than the specific stimuli tested. A neural measure of perceptual intelligibility may aid in the development of strategies for helping those with impaired speech perception.

Hemodynamic timing in resting-state and breathing-task BOLD fMRI.

The blood flow response to a vasoactive stimulus demonstrates regional heterogeneity across both the healthy brain and in cerebrovascular pathology. The timing of a regional hemodynamic response is emerging as an important biomarker of cerebrovascular dysfunction, as well as a confound within fMRI analyses. Previous research demonstrated that hemodynamic timing is more robustly characterized when a larger systemic vascular response is evoked by a breathing challenge, compared to when only spontaneous fluctuations in vascular physiology are present (i.e., in resting-state data). However, it is not clear whether hemodynamic delays in these two conditions are physiologically interchangeable, and how methodological signal-to-noise factors may limit their agreement. To address this, we generated whole-brain maps of hemodynamic delays in nine healthy adults. We assessed the agreement of voxel-wise gray matter (GM) hemodynamic delays between two conditions: resting-state and breath-holding. We found that delay values demonstrated poor agreement when considering all GM voxels, but increasingly greater agreement when limiting analyses to voxels showing strong correlation with the GM mean time-series. Voxels showing the strongest agreement with the GM mean time-series were primarily located near large venous vessels, however these voxels explain some, but not all, of the observed agreement in timing. Increasing the degree of spatial smoothing of the fMRI data enhanced the correlation between individual voxel time-series and the GM mean time-series. These results suggest that signal-to-noise factors may be limiting the accuracy of voxel-wise timing estimates and hence their agreement between the two data segments. In conclusion, caution must be taken when using voxel-wise delay estimates from resting-state and breathing-task data interchangeably, and additional work is needed to evaluate their relative sensitivity and specificity to aspects of vascular physiology and pathology.

White matter BOLD signals at 7 Tesla reveal visual field maps in optic radiation and vertical occipital fasciculus.

There is growing evidence that blood-oxygen-level-dependent (BOLD) activity in the white matter (WM) can be detected by functional magnetic resonance imaging (fMRI). However, the functional relevance and significance of WM BOLD signals remain controversial. Here we investigated whether 7T BOLD fMRI can reveal fine-scale functional organizations of a WM bundle. Population receptive field (pRF) analyses of the 7T retinotopy dataset from the Human Connectome Project revealed clear contralateral retinotopic organizations of two visual WM bundles: the optic radiation (OR) and the vertical occipital fasciculus (VOF). The retinotopic maps of OR are highly consistent with post-mortem dissections and diffusion tractographies, while the VOF maps are compatible with the dorsal and ventral visual areas connected by the WM. Similar to the grey matter (GM) visual areas, both WM bundles show over-representations of the central visual field and increasing pRF size with eccentricity. Hemodynamic response functions of visual WM were slower and wider compared with those of GM areas. These findings clearly demonstrate that WM BOLD at 7 Tesla is closely coupled with neural activity related to axons, encoding highly specific information that can be used to characterize fine-scale functional organizations of a WM bundle.

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI.

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO2 (PETCO2) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable PETCO2 recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with PETCO2. We define a threshold to determine whether a PETCO2 recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient PETCO2 and 6 datasets with insufficient PETCO2. When PETCO2 quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from PETCO2 (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to PETCO2 with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when PETCO2 quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient PETCO2 recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

A selection and targeting framework of cortical locations for line-scanning fMRI.

Depth-resolved functional magnetic resonance imaging (fMRI) is an emerging field growing in popularity given the potential of separating signals from different computational processes in cerebral cortex. Conventional acquisition schemes suffer from low spatial and temporal resolutions. Line-scanning methods allow depth-resolved fMRI by sacrificing spatial coverage to sample blood oxygenated level-dependent (BOLD) responses at ultra-high temporal and spatial resolution. For neuroscience applications, it is critical to be able to place the line accurately to (1) sample the right neural population and (2) target that neural population with tailored stimuli or tasks. To this end, we devised a multi-session framework where a target cortical location is selected based on anatomical and functional properties. The line is then positioned according to this information in a separate second session, and we tailor the experiment to focus on the target location. Anatomically, the precision of the line placement was confirmed by projecting a nominal representation of the acquired line back onto the surface. Functional estimates of neural selectivities in the line, as quantified by a visual population-receptive field model, resembled the target selectivities well for most subjects. This functional precision was quantified in detail by estimating the distance between the visual field location of the targeted vertex and the location in visual cortex (V1) that most closely resembled the line-scanning estimates; this distance was on average ~5.5 mm. Given the dimensions of the line, differences in acquisition, session, and stimulus design, this validates that line-scanning can be used to probe local neural sensitivities across sessions. In summary, we present an accurate framework for line-scanning MRI; we believe such a framework is required to harness the full potential of line-scanning and maximize its utility. Furthermore, this approach bridges canonical fMRI experiments with electrophysiological experiments, which in turn allows novel avenues for studying human physiology non-invasively.

Rapid high-fidelity T 2 * mapping using single-shot overlapping-echo acquisition and deep learning reconstruction.

PURPOSE: To develop and evaluate a single-shot quantitative MRI technique called GRE-MOLED (gradient-echo multiple overlapping-echo detachment) for rapid T 2 * $$ {T}_2^{\ast } $$ mapping. METHODS: In GRE-MOLED, multiple echoes with different TEs are generated and captured in a single shot of the k-space through MOLED encoding and EPI readout. A deep neural network, trained by synthetic data, was employed for end-to-end parametric mapping from overlapping-echo signals. GRE-MOLED uses pure GRE acquisition with a single echo train to deliver T 2 * $$ {T}_2^{\ast } $$ maps less than 90 ms per slice. The self-registered B0 information modulated in image phase was utilized for distortion-corrected parametric mapping. The proposed method was evaluated in phantoms, healthy volunteers, and task-based FMRI experiments. RESULTS: The quantitative results of GRE-MOLED T 2 * $$ {T}_2^{\ast } $$ mapping demonstrated good agreement with those obtained from the multi-echo GRE method (Pearson's correlation coefficient = 0.991 and 0.973 for phantom and in vivo brains, respectively). High intrasubject repeatability (coefficient of variation <1.0%) were also achieved in scan-rescan test. Enabled by deep learning reconstruction, GRE-MOLED showed excellent robustness to geometric distortion, noise, and random subject motion. Compared to the conventional FMRI approach, GRE-MOLED also achieved a higher temporal SNR and BOLD sensitivity in task-based FMRI. CONCLUSION: GRE-MOLED is a new real-time technique for T 2 * $$ {T}_2^{\ast } $$ quantification with high efficiency and quality, and it has the potential to be a better quantitative BOLD detection method.

BOLD fMRI and DTI fiber tracking for preoperative mapping of eloquent cerebral regions in brain tumor patients: impact on surgical approach and outcome.

PURPOSE: Task-based BOLD fMRI and DTI-fiber tracking have become part of the routine presurgical work-up of brain tumor patients in many institutions. However, their potential impact on both surgical treatment and neurologic outcome remains unclear, in despite of the high costs and complex implementation. METHODS: We retrospectively investigated whether performing fMRI and DTI-ft preoperatively substantially impacted surgical planning and patient outcome in a series of brain tumor patients. We assessed (i) the quality of fMRI and DTI-ft results, by using a scale of 0-2 (0 = failed mapping; 1 = intermediate confidence; 2 = good confidence), (ii) whether functional planning substantially contributed to defining the surgical strategy to be undertaken (i.e., no surgery, biopsy, or resection, with or without ESM), the surgical entry point and extent of resection, and (iii) the incidence of neurological deficits post-operatively. RESULTS: Twenty-seven patients constituted the study population. The mean confidence rating was 1.9/2 for fMRI localization of the eloquent cortex and lateralization of the language function and 1.7/2 for DTI-ft results. Treatment strategy was altered in 33% (9/27) of cases. Surgical entry point was modified in 8% (2/25) of cases. The extent of resection was modified in 40% (10/25). One patient (1/25, 4%) developed one new functional deficit post-operatively. CONCLUSION: Functional MR mapping - which must not be considered an alternative to ESM - has a critical role preoperatively, potentially modifying treatment strategy or increasing the neurosurgeons' confidence in the surgical approach hypothesized based on conventional imaging.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex.

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2-mm voxels) in central gray matter to minimize variations caused by partial-volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across-session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long-term use of event-related functional magnetic resonance imaging protocols.

Multi-echo balanced SSFP with a sequential phase-encoding order for functional MR imaging at 7T.

PURPOSE: To develop a 2D multi-echo passband balanced SSFP (bSSFP) sequence using an echo-train readout with a sequential phase-encoding order (sequential multi-echo bSSFP), and evaluate its performance in fast functional brain imaging at 7 T. METHODS: As images of sequential multi-echo bSSFP exhibit multiple ghosts due to periodic k-space modulations, a GRAPPA-based reconstruction method was proposed to eliminate ghosting artifacts. MRI experiments were performed to compare the image quality of multi-echo bSSFP and conventional single-echo bSSFP. Submillimeter-resolution fMRI using a checkerboard visual stimulus was conducted to compare the activation characteristics of multi-echo bSSFP, conventional single-echo bSSFP and standard gradient-echo EPI (GE-EPI). RESULTS: A higher mean structural similarity index was found between images of single-echo bSSFP and multi-echo bSSFP with a shorter echo train length (ETL). Multi-echo bSSFP (ETL = 3) showed higher temporal SNR (tSNR) values than GRAPPA-accelerated single-echo bSSFP (R = 2). In submillimeter-resolution fMRI experiments, multi-echo bSSFP (ETL = 3) approached the imaging speed of GRAPPA-accelerated single-echo bSSFP (R = 2), but without tSNR penalty and reduced activation due to acceleration. The median t-value and the number of significantly activated voxels were comparable between GE-EPI and multi-echo bSSFP (ETL = 3) that provides virtually distortion-free functional images and inherits the activation patterns of conventional bSSFP. CONCLUSION: Sequential multi-echo bSSFP (ETL = 3) is suitable for fast fMRI with submillimeter in-plane resolution, and offers an option to accelerate bSSFP imaging without tSNR penalty like parallel imaging.

Metabolic and Hemodynamic Resting-State Connectivity of the Human Brain: A High-Temporal Resolution Simultaneous BOLD-fMRI and FDG-fPET Multimodality Study.

Simultaneous [18F]-fluorodeoxyglucose positron emission tomography functional magnetic resonance imaging (FDG-PET/fMRI) provides the capacity to image 2 sources of energetic dynamics in the brain-glucose metabolism and the hemodynamic response. fMRI connectivity has been enormously useful for characterizing interactions between distributed brain networks in humans. Metabolic connectivity based on static FDG-PET has been proposed as a biomarker for neurological disease, but FDG-sPET cannot be used to estimate subject-level measures of "connectivity," only across-subject "covariance." Here, we applied high-temporal resolution constant infusion functional positron emission tomography (fPET) to measure subject-level metabolic connectivity simultaneously with fMRI connectivity. fPET metabolic connectivity was characterized by frontoparietal connectivity within and between hemispheres. fPET metabolic connectivity showed moderate similarity with fMRI primarily in superior cortex and frontoparietal regions. Significantly, fPET metabolic connectivity showed little similarity with FDG-sPET metabolic covariance, indicating that metabolic brain connectivity is a nonergodic process whereby individual brain connectivity cannot be inferred from group-level metabolic covariance. Our results highlight the complementary strengths of fPET and fMRI in measuring the intrinsic connectivity of the brain and open up the opportunity for novel fundamental studies of human brain connectivity as well as multimodality biomarkers of neurological diseases.

Assessment of olfactory information in the human brain using 7-Tesla functional magnetic resonance imaging.

Olfaction could prove to be an early marker of neurodegenerative diseases, including Alzheimer's and Parkinson's diseases. To use olfaction for disease diagnosis, elucidating the standard olfactory functions in healthy humans is necessary. However, the olfactory function in the human brain is less frequently assessed because of methodological difficulties associated with olfactory-related cerebral areas. Using ultra-high fields (UHF), functional magnetic resonance imaging (fMRI) with high spatial resolution and sensitivity may allow for the measurement of activation in the cerebral areas. This study aimed to apply 7-Tesla fMRI to assess olfactory function in the human brain by exposing individuals to four different odorants for 8 s. We found that olfactory stimulation mainly activated the piriform and orbitofrontal cortex in addition to the amygdala. Among these regions, univariate fMRI analysis indicated that subjective odor intensity significantly correlated with the averaged fMRI signals in the piriform cortex but not with subjective hedonic tone in any region. In contrast, multivariate fMRI analysis showed that subjective hedonic tone could be discriminated from the fMRI response patterns in the posterior orbitofrontal cortex. Thus, the piriform cortex is mainly associated with subjective odor intensity, whereas the posterior orbitofrontal cortex are involved in the discrimination of the subjective hedonic tone of the odorant. UHF-fMRI may be useful for assessing olfactory function in the human brain.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function.

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2-3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO2 pressure: a breath-hold task to induce hypercapnia (CO2 increase) and a cued deep breathing task to induce hypocapnia (CO2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO2 by systematically shifting the CO2 regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

Resting-State functional networks of different topographic representations in the somatosensory cortex of macaque monkeys and humans.

Information processing in the brain is mediated through a complex functional network architecture whose comprising nodes integrate and segregate themselves on different timescales. To gain an understanding of the network function it is imperative to identify and understand the network structure with respect to the underlying anatomical connectivity and the topographic organization. Here we show that the previously described resting-state network for the somatosensory area 3b comprises of distinct networks that are characteristic for different topographic representations. Seed-based resting-state functional connectivity analysis in macaque monkeys and humans using BOLD-fMRI signals from the face, the hand and rest of the medial somatosensory representations of area 3b revealed different correlation patterns. Both monkeys and humans have many similarities in the connectivity networks, although the networks are more complex in humans with many more nodes. In both the species face area network has the highest ipsilateral and contralateral connectivity, which included areas 3b and 4, and ventral premotor area. The area 3b hand network included ipsilateral hand representation in area 4. The emergent functional network structures largely reflect the known anatomical connectivity. Our results show that different body part representations in area 3b have independent functional networks perhaps reflecting differences in the behavioral use of different body parts. The results also show that large cortical areas if considered together, do not give a complete and accurate picture of the network architecture.