Functional brain mapping by blood oxygenation level-dependent contrast magnetic resonance imaging. A comparison of signal characteristics with a biophysical model

Human brain functional MRS reveals interplay of metabolites implicated in neurotransmission and neuroenergetics

While functional MRI (fMRI) localizes brain activation and deactivation, functional MRS (fMRS) provides insights into the underlying metabolic conditions. There is much interest in measuring task-induced and resting levels of metabolites implicated in neuroenergetics (e.g., lactate, glucose, or β-hydroxybutyrate (BHB)) and neurotransmission (e.g., γ-aminobutyric acid (GABA) or pooled glutamate and glutamine (Glx)). Ultra-high magnetic field (e.g., 7T) has boosted the fMRS quantification precision, reliability, and stability of spectroscopic observations using short echo-time (TE) ¹H-MRS techniques. While short TE ¹H-MRS lacks sensitivity and specificity for fMRS at lower magnetic fields (e.g., 3T or 4T), most of these metabolites can also be detected by J-difference editing (JDE) ¹H-MRS with longer TE to filter overlapping resonances. The ¹H-MRS studies show that JDE can detect GABA, Glx, lactate, and BHB at 3T, 4T and 7T. Most recently, it has also been demonstrated that JDE ¹H-MRS is capable of reliable detection of metabolic changes in different brain areas at various magnetic fields. Combining fMRS measurements with fMRI is important for understanding normal brain function, but also clinically relevant for mechanisms and/or biomarkers of neurological and neuropsychiatric disorders. We provide an up-to-date overview of fMRS research in the last three decades, both in terms of applications and technological advances. Overall the emerging fMRS techniques can be expected to contribute substantially to our understanding of metabolism for brain function and dysfunction.

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.

High-resolution relaxometry-based calibrated fMRI in murine brain: Metabolic differences between awake and anesthetized states

Functional magnetic resonance imaging (fMRI) techniques using the blood-oxygen level-dependent (BOLD) signal have shown great potential as clinical biomarkers of disease. Thus, using these techniques in preclinical rodent models is an urgent need. Calibrated fMRI is a promising technique that can provide high-resolution mapping of cerebral oxygen metabolism (CMR_O2). However, calibrated fMRI is difficult to use in rodent models for several reasons: rodents are anesthetized, stimulation-induced changes are small, and gas challenges induce noisy CMR_O2 predictions. We used, in mice, a relaxometry-based calibrated fMRI method which uses cerebral blood flow (CBF) and the BOLD-sensitive magnetic relaxation component, R₂', the same parameter derived in the deoxyhemoglobin-dilution model of calibrated fMRI. This method does not use any gas challenges, which we tested on mice in both awake and anesthetized states. As anesthesia induces a whole-brain change, our protocol allowed us to overcome the former limitations of rodent studies using calibrated fMRI. We revealed 1.5-2 times higher CMR_O2, dependent upon brain region, in the awake state versus the anesthetized state. Our results agree with alternative measurements of whole-brain CMR_O2 in the same mice and previous human anesthesia studies. The use of calibrated fMRI in rodents has much potential for preclinical fMRI.

From a Demand-Based to a Supply-Limited Framework of Brain Metabolism

What defines the rate of energy use by the brain, as well as per neurons of different sizes in different structures and animals, is one fundamental aspect of neuroscience for which much has been theorized, but very little data are available. The prevalent theories and models consider that energy supply from the vascular system to different brain regions is adjusted both dynamically and in the course of development and evolution to meet the demands of neuronal activity. In this perspective, we offer an alternative view: that regional rates of energy use might be mostly constrained by supply, given the properties of the brain capillary network, the highly stable rate of oxygen delivery to the whole brain under physiological conditions, and homeostatic constraints. We present evidence that these constraints, based on capillary density and tissue oxygen homeostasis, are similar between brain regions and mammalian species, suggesting they derive from fundamental biophysical limitations. The same constraints also determine the relationship between regional rates of brain oxygen supply and usage over the full physiological range of brain activity, from deep sleep to intense sensory stimulation, during which the apparent uncoupling of blood flow and oxygen use is still a predicted consequence of supply limitation. By carefully separating "energy cost" into energy supply and energy use, and doing away with the problematic concept of energetic "demands," our new framework should help shine a new light on the neurovascular bases of metabolic support of brain function and brain functional imaging. We speculate that the trade-offs between functional systems and even the limitation to a single attentional spot at a time might be consequences of a strongly supply-limited brain economy. We propose that a deeper understanding of brain energy supply constraints will provide a new evolutionary understanding of constraints on brain function due to energetics; offer new diagnostic insight to disturbances of brain metabolism; lead to clear, testable predictions on the scaling of brain metabolic cost and the evolution of brains of different sizes; and open new lines of investigation into the microvascular bases of progressive cognitive loss in normal aging as well as metabolic diseases.

Self-navigated prospective motion correction for 3D-EPI acquisition

PURPOSE

Although 3D EPI is more susceptible to motion artifacts than 2D EPI, it presents some benefits for functional MRI, including the absence of spin-history artifacts, greater potential for parallel imaging acceleration, and better functional sensitivity in high-resolution imaging. Here we present a self-navigated 3D-EPI sequence suitable for prospective motion-corrected functional MRI without additional hardware or pulses.

METHODS

For each volume acquisition, the first 24 of the 52 partitions being acquired are accumulated to a new feedback block that was added to the image reconstruction pipeline. After zero-filling the remaining partitions, the feedback block constructs a volumetric self-navigator (vSNav), co-registers it to the reference vSNav acquired during the first volume acquisition, and sends motion estimates to the sequence. The sequence then updates its FOV and acquires subsequent partitions with the adjusted FOV, until the next update is received. The sequence was validated without and with intentional motion in phantom and in vivo on a 3T Skyra.

RESULTS

For phantom scans, the FOV was updated 0.704 s after acquisition of the vSNav partitions, and for in vivo scans after 0.768 s. Both phantom and in vivo data demonstrated stable motion estimates in the absence of motion. For in vivo acquisitions, prospective head-pose estimates using the vSNav's and retrospective estimates with FLIRT (FMRIB's Linear Image Registration Tool) agreed to within 0.23 mm (< 10% of the slice thickness) and 0.14° in all directions.

CONCLUSION

Depending when motion occurs during a volume acquisition, the proposed method fully corrects the FOV and recovers image quality within one volume acquisition.

Association between tissue hypoxia, perfusion restrictions, and microvascular architecture alterations with lesion-induced impairment of neurovascular coupling

Functional magnetic resonance imaging (fMRI) has been mainly utilized for the preoperative localization of eloquent cortical areas. However, lesion-induced impairment of neurovascular coupling (NVC) in the lesion border zone may lead to false-negative fMRI results. The purpose of this study was to determine physiological factors impacting the NVC. Twenty patients suffering from brain lesions were preoperatively examined using multimodal neuroimaging including fMRI, magnetoencephalography (MEG) during language or sensorimotor tasks (depending on lesion location), and a novel physiologic MRI approach for the combined quantification of oxygen metabolism, perfusion state, and microvascular architecture. Congruence of brain activity patterns between fMRI and MEG were found in 13 patients. In contrast, we observed missing fMRI activity in perilesional cortex that demonstrated MEG activity in seven patients, which was interpreted as lesion-induced impairment of NVC. In these brain regions with impaired NVC, physiologic MRI revealed significant brain tissue hypoxia, as well as significantly decreased macro- and microvascular perfusion and microvascular architecture. We demonstrated that perilesional hypoxia with reduced vascular perfusion and architecture is associated with lesion-induced impairment of NVC. Our physiologic MRI approach is a clinically applicable method for preoperative risk assessment for the presence of false-negative fMRI results and may prevent severe postoperative functional deficits.

Advanced Imaging Techniques for Newly Diagnosed and Recurrent Gliomas

Management of gliomas following initial diagnosis requires thoughtful presurgical planning followed by regular imaging to monitor treatment response and survey for new tumor growth. Traditional MR imaging modalities such as T1 post-contrast and T2-weighted sequences have long been a staple of tumor diagnosis, surgical planning, and post-treatment surveillance. While these sequences remain integral in the management of gliomas, advances in imaging techniques have allowed for a more detailed characterization of tumor characteristics. Advanced MR sequences such as perfusion, diffusion, and susceptibility weighted imaging, as well as PET scans have emerged as valuable tools to inform clinical decision making and provide a non-invasive way to help distinguish between tumor recurrence and pseudoprogression. Furthermore, these advances in imaging have extended to the operating room and assist in making surgical resections safer. Nevertheless, surgery, chemotherapy, and radiation treatment continue to make the interpretation of MR changes difficult for glioma patients. As analytics and machine learning techniques improve, radiomics offers the potential to be more quantitative and personalized in the interpretation of imaging data for gliomas. In this review, we describe the role of these newer imaging modalities during the different stages of management for patients with gliomas, focusing on the pre-operative, post-operative, and surveillance periods. Finally, we discuss radiomics as a means of promoting personalized patient care in the future.

Personalized fMRI delineates functional regions preserved within brain tumors

Objective

Methods

Results

Interpretation

Reference-Free Adaptive Filtering of Extracellular Neural Signals Recording in Ultra-High Field Magnetic Resonance Imaging Scanners: Removal of Periodic Interferences

This paper focuses on the removal of periodic artifacts from neural signals recorded in rats in ultra-high field (UHF) MRI scanners, using a reference free adaptive feedforward method. Recording extracellular neural signals in the UHF environment is motivated by the desire to combine neural recording and UHF functional magnetic resonance imaging (fMRI) to better understand brain function. However, the neural signals are found to have extremely high noise artifacts of a periodic nature due to electromagnetic interference and due to small oscillatory motions. In particular, noise at 60 Hz and several harmonics of 60 Hz, sinusoidal noise from a pump, and low frequency breathing motion artifacts are observed. Due to significant overlap between the noise frequencies and the neural frequency region of interest, band pass filters cannot be effectively utilized in this application. Hence, this paper develops adaptive least squares feedforward cancellation filters to remove the periodic artifacts. The interference fundamental frequency is identified precisely using an implementation of k-means in an iterative approach. The paper includes significant animal data from rats recorded in an IACUC-approved procedure in 9.4T and 16.4T MRI machines. For breathing artifacts filtered from 4 rats, the mean signal cancellation values at the harmonic interference frequencies are 5.18, 12.97, and 20.87 dB/Hz for a sliding template subtraction, a single-stage impulse reference method, and the cascaded adaptive filtering approach respectively. For pump artifacts filtered from 2 chronically implanted rats, mean signal cancellation values are 2.85, 9.52 and 12.06 dB/Hz respectively. The experimental results show that periodic noise is very effectively removed by the developed cascaded adaptive least squares feedforward algorithm.

Remote ischemic postconditioning increased cerebral blood flow and oxygenation assessed by magnetic resonance imaging in newborn piglets after hypoxia-ischemia

BACKGROUND

We have previously investigated neurological outcomes following remote ischemic postconditioning (RIPC) in a newborn piglet model of hypoxic-ischemic encephalopathy. The aim of this study was to further investigate potential mechanisms of neuroprotection by comparing newborn piglets subjected to global hypoxia-ischemia (HI) treated with and without RIPC with regards to measures of cerebral blood flow and oxygenation assessed by functional magnetic resonance imaging.

MATERIALS AND METHODS

A total of 50 piglets were subjected to 45 min global HI and randomized to either no treatment or RIPC treatment. Magnetic resonance imaging was performed 72 h after the HI insult with perfusion-weighted (arterial spin labeling, ASL) and oxygenation-weighted (blood-oxygen-level-dependent, BOLD) sequences in the whole brain, basal ganglia, thalamus, and cortex. Four sham animals received anesthesia and mechanical ventilation only.

RESULTS

Piglets treated with RIPC had higher measures of cerebral blood flow in all regions of interest and the whole brain (mean difference: 2.6 ml/100 g/min, 95% CI: 0.1; 5.2) compared with the untreated controls. They also had higher BOLD values in the basal ganglia and the whole brain (mean difference: 4.2 T2*, 95% CI: 0.4; 7.9). Measures were similar between piglets treated with RIPC and sham animals.

CONCLUSION

Piglets treated with RIPC had higher measures of cerebral blood flow and oxygenation assessed by magnetic resonance imaging in the whole brain and several regions of interest compared with untreated controls 72 h after the HI insult. Whether this reflects a potential neuroprotective mechanism of RIPC requires further study.

Imaging human engrams using 7 Tesla magnetic resonance imaging

The investigation of the physical traces of memories (engrams) has made significant progress in the last decade due to optogenetics and fluorescent cell tagging applied in rodents. Engram cells were identified. The ablation of engram cells led to the loss of the associated memory, silent memories were reactivated, and artificial memories were implanted in the brain. Human engram research lags behind engram research in rodents due to methodological and ethical constraints. However, advances in multivariate analysis techniques of functional magnetic resonance imaging (fMRI) data and machine learning algorithms allowed the identification of stable engram patterns in humans. In addition, MRI scanners with an ultrahigh field strength of 7 Tesla (T) have left their prototype state and became more common around the world to assist human engram research. Although most engram research in humans is still being performed with a field strength of 3T, fMRI at 7T will push engram research. Here, we summarize the current state and findings of human engram research and discuss the advantages and disadvantages of applying 7 versus 3T fMRI to image human memory traces.

Cortical Surface-Informed Volumetric Spatial Smoothing of fMRI Data via Graph Signal Processing

Conventionally, as a preprocessing step, functional MRI (fMRI) data are spatially smoothed before further analysis, be it for activation mapping on task-based fMRI or functional connectivity analysis on resting-state fMRI data. When images are smoothed volumetrically, however, isotropic Gaussian kernels are generally used, which do not adapt to the underlying brain structure. Alternatively, cortical surface smoothing procedures provide the benefit of adapting the smoothing process to the underlying morphology, but require projecting volumetric data on to the surface. In this paper, leveraging principles from graph signal processing, we propose a volumetric spatial smoothing method that takes advantage of the gray-white and pial cortical surfaces, and as such, adapts the filtering process to the underlying morphological details at each point in the cortex.

Entropic regression with neurologically motivated applications

The ultimate goal of cognitive neuroscience is to understand the mechanistic neural processes underlying the functional organization of the brain. The key to this study is understanding the structure of both the structural and functional connectivity between anatomical regions. In this paper, we use an information theoretic approach, which defines direct information flow in terms of causation entropy, to improve upon the accuracy of the recovery of the true network structure over popularly used methods for this task such as correlation and least absolute shrinkage and selection operator regression. The method outlined above is tested on synthetic data, which is produced by following previous work in which a simple dynamical model of the brain is used, simulated on top of a real network of anatomical brain regions reconstructed from diffusion tensor imaging. We demonstrate the effectiveness of the method of AlMomani et al. [Chaos 30, 013107 (2020)] when applied to data simulated on the realistic diffusion tensor imaging network, as well as on randomly generated small-world and Erdös-Rényi networks.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

Simultaneous pure T2 and varying T2'-weighted BOLD fMRI using Echo Planar Time-resolved Imaging for mapping cortical-depth dependent responses

Spin-echo (SE) BOLD fMRI has high microvascular specificity, and thus provides a more reliable means to localize neural activity compared to conventional gradient-echo BOLD fMRI. However, the most common SE BOLD acquisition method, SE-EPI, is known to suffer from T₂′ contrast contamination with undesirable draining vein bias. To address this, in this study, we extended a recently developed distortion/blurring-free multi-shot EPI technique, Echo-Planar Time-resolved Imaging (EPTI), to cortical-depth dependent SE-fMRI at 7T to test whether it could provide purer SE BOLD contrast with minimal T₂′ contamination for improved neuronal specificity. From the same acquisition, the time-resolved feature of EPTI also provides a series of asymmetric SE (ASE) images with varying T₂′ weightings, and enables extraction of data equivalent to conventional SE EPI with different echo train lengths (ETLs). This allows us to systematically examine how T₂′-contribution affects different SE acquisition strategies using a single dataset. A low-rank spatiotemporal subspace reconstruction was implemented for the SE-EPTI acquisition, which incorporates corrections for both shot-to-shot phase variations and dynamic B₀ drifts. SE-EPTI was used in a visual task fMRI experiment to demonstrate that i) the pure SE image provided by EPTI results in the highest microvascular specificity; ii) the ASE EPTI series, with a graded introduction of T₂′ weightings at time points farther away from the pure SE, show a gradual sensitivity increase along with increasing draining vein bias; iii) the longer ETL seen in conventional SE EPI acquisitions will induce more draining vein bias. Consistent results were observed across multiple subjects, demonstrating the robustness of the proposed technique for SE-BOLD fMRI with high specificity.

Impaired cerebral vascular and metabolic responses to parametric N-back tasks in subjective cognitive decline

Previous studies reported abnormally increased and/or decreased blood oxygen level-dependent (BOLD) activations during functional tasks in subjective cognitive decline (SCD). The neurophysiological basis underlying these functional aberrations remains debated. This study aims to investigate vascular and metabolic responses and their dependence on cognitive processing loads during functional tasks in SCD. Twenty-one SCD and 18 control subjects performed parametric N-back working-memory tasks during MRI scans. Task-evoked percentage changes (denoted as δ) in cerebral blood volume (δCBV), cerebral blood flow (δCBF), BOLD signal (δBOLD) and cerebral metabolic rate of oxygen (δCMRO₂) were evaluated. In the frontal lobe, trends of decreased δCBV, δCBF and δCMRO₂ and increased δBOLD were observed in SCD compared with control subjects under lower loads, and these trends increased to significant differences under the 3-back load. δCBF was significantly correlated with δCMRO₂ in controls, but not in SCD subjects. As N-back loads increased, the differences between SCD and control subjects in δCBF and δCMRO₂ tended to enlarge. In the parietal lobe, no significant between-group difference was observed. Our findings suggested that impaired vascular and metabolic responses to functional tasks occurred in the frontal lobe of SCD, which contributed to unusual BOLD hyperactivation and was modulated by cognitive processing loads.

Topological Receptive Field Model for Human Retinotopic Mapping

The mapping between visual inputs on the retina and neuronal activations in the visual cortex, i.e., retinotopic map, is an essential topic in vision science and neuroscience. Human retinotopic maps can be revealed by analyzing the functional magnetic resonance imaging (fMRI) signal responses to designed visual stimuli in vivo. Neurophysiology studies summarized that visual areas are topological (i.e., nearby neurons have receptive fields at nearby locations in the image). However, conventional fMRI-based analyses frequently generate non-topological results because they process fMRI signals on a voxel-wise basis, without considering the neighbor relations on the surface. Here we propose a topological receptive field (tRF) model which imposes the topological condition when decoding retinotopic fMRI signals. More specifically, we parametrized the cortical surface to a unit disk, characterized the topological condition by tRF, and employed an efficient scheme to solve the tRF model. We tested our framework on both synthetic and human fMRI data. Experimental results showed that the tRF model could remove the topological violations, improve model explaining power, and generate biologically plausible retinotopic maps. The proposed framework is general and can be applied to other sensory maps.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

Advances in resting state fMRI acquisitions for functional connectomics

Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.

Effects of mild hypoxia on oxygen extraction fraction responses to brain stimulation

Characterizing the effect of limited oxygen availability on brain metabolism during brain activation is an essential step towards a better understanding of brain homeostasis and has obvious clinical implications. However, how the cerebral oxygen extraction fraction (OEF) depends on oxygen availability during brain activation remains unclear, which is mostly attributable to the scarcity and safety of measurement techniques. Recently, a magnetic resonance imaging (MRI) method that enables noninvasive and dynamic measurement of the OEF has been developed and confirmed to be applicable to functional MRI studies. Using this novel method, the present study investigated the motor-evoked OEF response in both normoxia (21% O₂) and hypoxia (12% O₂). Our results showed that OEF activation decreased in the brain areas involved in motor task execution. Decreases in the motor-evoked OEF response were greater under hypoxia (-21.7% ± 5.5%) than under normoxia (-11.8% ± 3.7%) and showed a substantial decrease as a function of arterial oxygen saturation. These findings suggest a different relationship between oxygen delivery and consumption during hypoxia compared to normoxia. This methodology may provide a new perspective on the effects of mild hypoxia on brain function.

Development of functional connectivity within and among the resting-state networks in anesthetized rhesus monkeys

OBJECTIVE

The age-related changes in the resting-state networks (RSNs) exhibited temporally specific patterns in humans, and humans and rhesus monkeys have similar RSNs. We hypothesized that the RSNs in rhesus monkeys experienced similar developmental patterns as humans.

METHODS

We acquired resting-state fMRI data from 62 rhesus monkeys, which were divided into childhood, adolescence, and early adulthood groups. Group independent component analysis (ICA) was used to identify monkey RSNs. We detected the between-group differences in the RSNs and static, dynamic, and effective functional network connections (FNCs) using one-way variance analysis (ANOVA) and post-hoc analysis.

RESULTS

Eight rhesus RSNs were identified, including cerebellum (CN), left and right lateral visual (LVN and RVN), posterior default mode (pDMN), visuospatial (VSN), frontal (FN), salience (SN), and sensorimotor networks (SMN). In internal connections, the CN, SN, FN, and SMN mainly matured in early adulthood. The static FNCs associated with FN, SN, pDMN primarily experienced fast descending slow ascending type (U-shaped) developmental patterns for maturation, and the dynamic FNCs related to pDMN (RVN, CN, and SMN) and SMN (CN) were mature in early adulthood. The effective FNC results showed that the pDMN and VSN (stimulated), SN (inhibited), and FN (first inhibited then stimulated) chiefly matured in early adulthood.

CONCLUSION

We identified eight monkey RSNs, which exhibited similar development patterns as humans. All the RSNs and FNCs in monkeys were not widely changed but fine-tuned. Our study clarified that the progressive synchronization, exploration, and regulation of cognitive RSNs within the pDMN, FN, SN, and VSN denoted potential maturation of the RSNs throughout development. We confirmed the development patterns of RSNs and FNCs would support the use of monkeys as a best animal model for human brain function.

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.

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

The Effect of Crystal Dependence on Brain Activity Related to the Perception of Pleasure Using fMRI

Background: Addiction is currently one of the problems of human society. Drug abuse is one of the most important issues in the field of addiction. Methamphetamine (crystal) is one of the drugs that has been abused in recent decades. Methods: In this case-control study, 10 individuals aged 20 to 40 years old with at least 2 years of experience of methamphetamine consumption without any history of drug use or other stimulants from clients and drug withdrawal centers in Tehran City, and 10 healthy volunteers were selected. Age, social status, and economic status of addicts were included in the fMRI apparatus, and 90 selected pleasurable, non-pleasurable, and neutral images (IAPS) were displayed by the projector through an event-related method. The playback time of each photo was 3 s, and after this process, the person outside the device, without the time limit selected the enjoyable and unpleasant images. Results: The results showed that there was no significant difference between the groups in terms of age, alcohol use, and smoking history (P < 0.05). There was no significant difference in terms of the age at first use between members of the methamphetamine-dependent group. Also, the methamphetamine-dependent group showed more brain activity in their pre-center and post-center gyrus than the normal (control) group. Conclusions: According to the results obtained in this study, in general, it can be concluded that there are some areas in the brain of addicts that are activated when watching pleasant photos, while these areas are not active in the brains of normal people.

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.

Impaired response of cerebral oxygen metabolism to visual stimulation in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disease caused by a CAG triplet repeat expansion in the Huntingtin gene. Metabolic and microvascular abnormalities in the brain may contribute to early physiological changes that subserve the functional impairments in HD. This study is intended to investigate potential abnormality in dynamic changes in cerebral blood volume (CBV) and cerebral blood flow (CBF), and cerebral metabolic rate of oxygen (CMRO2) in the brain in response to functional stimulation in premanifest and early manifest HD patients. A recently developed 3-D-TRiple-acquisition-after-Inversion-Preparation magnetic resonance imaging (MRI) approach was used to measure dynamic responses in CBV, CBF, and CMRO2 during visual stimulation in one single MRI scan. Experiments were conducted in 23 HD patients and 16 healthy controls. Decreased occipital cortex CMRO2 responses were observed in premanifest and early manifest HD patients compared to controls (P < 0.001), correlating with the CAG-Age Product scores in these patients (R2 = 0.4, P = 0.001). The results suggest the potential value of this reduced CMRO2 response during visual stimulation as a biomarker for HD and may illuminate the role of metabolic alterations in the pathophysiology of HD.

Changes in Temporal and Spatial Patterns of Intrinsic Brain Activity and Functional Connectivity in Upper-Limb Amputees: An fMRI Study

Background

Amputation in adults is a serious procedure or traumatic outcome, one that leads to a possible “remapping” of limb representations (somatotopy) in the motor and sensory cortex. The temporal and spatial extent underlying reorganization of somatotopy is unclear. The aim of this study was to better understand how local and global structural plasticity in sensory-motor cortical networks changes temporally and spatially after upper-limb amputation.

Methods

We studied 8 healthy nonamputee control subjects and 16 complete upper-limb amputees. Resting-state MRI (rs-fMRI) was used to measure local and large-scale relative differences (compared to controls) in both the amplitude of low-frequency fluctuations (ALFF) and degree of centrality (DC) at 2 months, 6 months, and 12 months after traumatic amputation.

Results

In amputees, rs-fMRI scans revealed differences in spatial patterns of ALFF and DC among brain regions over time. Significant relative increases in ALFF and DC were detected not only in the sensory and motor cortex but also in related cortical regions believed to be involved in cognition and motor planning. We observed changes in the magnitude of ALFFs in the pre- and postcentral gyrus and primary sensory cortex, as well as in the anterior cingulate, parahippocampal gyrus, and hippocampus, 2 months after the amputation. The regional distribution of increases/decreases in ALFFs and DC documented at 2-month postamputation was very different from those at 6 and 12-month postamputation.

Conclusion

Local and wide-spread changes in ALFFs in the sensorimotor cortex and cognitive-related brain regions after upper-limb amputation may imply dysfunction not only in sensory and motor function but also in areas responsible for sensorimotor integration and motor planning. These results suggest that cortical reorganization after upper extremity deafferentation is temporally and spatially more complicated than previously appreciated, affecting DC in widespread regions.

Brain microvasculature has a common topology with local differences in geometry that match metabolic load

The microvasculature underlies the supply networks that support neuronal activity within heterogeneous brain regions. What are common versus heterogeneous aspects of the connectivity, density, and orientation of capillary networks? To address this, we imaged, reconstructed, and analyzed the microvasculature connectome in whole adult mice brains with sub-micrometer resolution. Graph analysis revealed common network topology across the brain that leads to a shared structural robustness against the rarefaction of vessels. Geometrical analysis, based on anatomically accurate reconstructions, uncovered a scaling law that links length density, i.e., the length of vessel per volume, with tissue-to-vessel distances. We then derive a formula that connects regional differences in metabolism to differences in length density and, further, predicts a common value of maximum tissue oxygen tension across the brain. Last, the orientation of capillaries is weakly anisotropic with the exception of a few strongly anisotropic regions; this variation can impact the interpretation of fMRI data.

An Overview of ICA/BSS-Based Application to Alzheimer's Brain Signal Processing

Alzheimer's disease (AD) is by far the most common cause of dementia associated with aging. Early and accurate diagnosis of AD and ability to track progression of the disease is increasingly important as potential disease-modifying therapies move through clinical trials. With the advent of biomedical techniques, such as computerized tomography (CT), electroencephalography (EEG), magnetoencephalography (MEG), positron emission tomography (PET), magnetic resonance imaging (MRI), and functional magnetic resonance imaging (fMRI), large amounts of data from Alzheimer's patients have been acquired and processed from which AD-related information or "signals" can be assessed for AD diagnosis. It remains unknown how best to mine complex information from these brain signals to aid in early diagnosis of AD. An increasingly popular technique for processing brain signals is independent component analysis or blind source separation (ICA/BSS) that separates blindly observed signals into original signals that are as independent as possible. This overview focuses on ICA/BSS-based applications to AD brain signal processing.

Characterizing changes in network connectivity following chronic head trauma in special forces military personnel: a combined resting-fMRI and DTI study

BACKGROUND

Soldiers are exposed to significant repetitive head trauma, which may disrupt functional and structural brain connectivity patterns.

PURPOSE/HYPOTHESIS

Integrate resting-state functional MRI (rs-fMRI) and diffusion tensor imaging (DTI) to characterize changes in connectivity biomarkers within Canadian Special Operations Forces (CANSOF), hypothesizing that alterations in architectural organization of cortical hubs may follow chronic repetitive head trauma.

METHODS

Fifteen CANSOFs with a history of chronic exposure to sub-concussive head trauma and concussive injuries (1.9 ± 2.0 concussions (range: [0-6])), as well as an equal age-matched cohort of controls (CTLs) were recruited. BOLD-based rs-fMRI was combined with DTI to reconstruct functional and structural networks using independent component analyses and probabilistic tractography. Connectivity markers were computed based on the distance between functional seeds to assess for possible differences in injury susceptibility of short- and long-range connections.

RESULTS/DISCUSSION

Significant hyper- and hypo-connectivity differences in cortical connections were observed suggesting that chronic head trauma may predispose soldiers to changes in the functional organization of brain networks. Significant structural alterations in axonal fibers directly connecting disrupted functional nodes were specific to hyper-connected long-range connections, suggesting a potential relationship between axonal injury and increases in neural recruitment following repetitive head trauma from high-exposure military duties.

Imaging the transmembrane and transendothelial sodium gradients in gliomas

Under normal conditions, high sodium (Na⁺) in extracellular (Na⁺_e) and blood (Na⁺_b) compartments and low Na⁺ in intracellular milieu (Na⁺_i) produce strong transmembrane (ΔNa⁺_mem) and weak transendothelial (ΔNa⁺_end) gradients respectively, and these manifest the cell membrane potential (V_m) as well as blood–brain barrier (BBB) integrity. We developed a sodium (²³Na) magnetic resonance spectroscopic imaging (MRSI) method using an intravenously-administered paramagnetic polyanionic agent to measure ΔNa⁺_mem and ΔNa⁺_end. In vitro ²³Na-MRSI established that the ²³Na signal is intensely shifted by the agent compared to other biological factors (e.g., pH and temperature). In vivo ²³Na-MRSI showed Na⁺_i remained unshifted and Na⁺_b was more shifted than Na⁺_e, and these together revealed weakened ΔNa⁺_mem and enhanced ΔNa⁺_end in rat gliomas (vs. normal tissue). Compared to normal tissue, RG2 and U87 tumors maintained weakened ΔNa⁺_mem (i.e., depolarized V_m) implying an aggressive state for proliferation, whereas RG2 tumors displayed elevated ∆Na⁺_end suggesting altered BBB integrity. We anticipate that ²³Na-MRSI will allow biomedical explorations of perturbed Na⁺ homeostasis in vivo.

Early fMRI responses to somatosensory and optogenetic stimulation reflect neural information flow

fMRI has revolutionized how neuroscientists investigate human brain functions and networks. To further advance understanding of brain functions, identifying the direction of information flow, such as thalamocortical versus corticothalamic projections, is critical. Because the early hemodynamic response at microvessels near active neurons can be detected by ultrahigh field fMRI, we propose using the onset times of fMRI responses to discern the information flow. This approach was confirmed by observing the ultrahigh spatiotemporal resolution BOLD fMRI responses to bottom-up somatosensory stimulation and top-down optogenetic stimulation of the primary motor cortex in anesthetized mice. Because ultrahigh field MRI is increasingly available, ultrahigh spatiotemporal fMRI will significantly facilitate the investigation of functional circuits in humans.

Blood oxygenation level–dependent (BOLD) functional magnetic resonance imaging (fMRI) has been widely used to localize brain functions. To further advance understanding of brain functions, it is critical to understand the direction of information flow, such as thalamocortical versus corticothalamic projections. For this work, we performed ultrahigh spatiotemporal resolution fMRI at 15.2 T of the mouse somatosensory network during forepaw somatosensory stimulation and optogenetic stimulation of the primary motor cortex (M1). Somatosensory stimulation induced the earliest BOLD response in the ventral posterolateral nucleus (VPL), followed by the primary somatosensory cortex (S1) and then M1 and posterior thalamic nucleus. Optogenetic stimulation of excitatory neurons in M1 induced the earliest BOLD response in M1, followed by S1 and then VPL. Within S1, the middle cortical layers responded to somatosensory stimulation earlier than the upper or lower layers, whereas the upper cortical layers responded earlier than the other two layers to optogenetic stimulation in M1. The order of early BOLD responses was consistent with the canonical understanding of somatosensory network connections and cannot be explained by regional variabilities in the hemodynamic response functions measured using hypercapnic stimulation. Our data demonstrate that early BOLD responses reflect the information flow in the mouse somatosensory network, suggesting that high-field fMRI can be used for systems-level network analyses.

In-vivo imaging of targeting and modulation of depression-relevant circuitry by transcranial direct current stimulation: a randomized clinical trial

Recent clinical trials of transcranial direct current stimulation (tDCS) in depression have shown contrasting results. Consequently, we used in-vivo neuroimaging to confirm targeting and modulation of depression-relevant neural circuitry by tDCS. Depressed participants (N = 66, Baseline Hamilton Depression Rating Scale (HDRS) 17-item scores ≥14 and <24) were randomized into Active/Sham and High-definition (HD)/Conventional (Conv) tDCS groups using a double-blind, parallel design, and received tDCS individually targeted at the left dorsolateral prefrontal cortex (DLPFC). In accordance with Ampere's Law, tDCS currents were hypothesized to induce magnetic fields at the stimulation-target, measured in real-time using dual-echo echo-planar-imaging (DE-EPI) MRI. Additionally, the tDCS treatment trial (consisting of 12 daily 20-min sessions) was hypothesized to induce cerebral blood flow (CBF) changes post-treatment at the DLPFC target and in the reciprocally connected anterior cingulate cortex (ACC), measured using pseudo-continuous arterial spin labeling (pCASL) MRI. Significant tDCS current-induced magnetic fields were observed at the left DLPFC target for both active stimulation montages (Brodmann's area (BA) 46: p_HD = 0.048, Cohen's d_HD = 0.73; p_Conv = 0.018, d_Conv = 0.86; BA 9: p_HD = 0.011, d_HD = 0.92; p_Conv = 0.022, d_Conv = 0.83). Significant longitudinal CBF increases were observed (a) at the left DLPFC stimulation-target for both active montages (p_HD = 3.5E-3, d_HD = 0.98; p_Conv = 2.8E-3, d_Conv = 1.08), and (b) at ACC for the HD-montage only (p_HD = 2.4E-3, d_HD = 1.06; p_Conv = 0.075, d_Conv = 0.64). These results confirm that tDCS-treatment (a) engages the stimulation-target, and (b) modulates depression-relevant neural circuitry in depressed participants, with stronger network-modulations induced by the HD-montage. Although not primary outcomes, active HD-tDCS showed significant improvements of anhedonia relative to sham, though HDRS scores did not differ significantly between montages post-treatment.

Insights Into Cerebral Tissue-Specific Response to Respiratory Challenges at 7T: Evidence for Combined Blood Flow and CO2-Mediated Effects

Cerebrovascular reactivity (CVR) mapping is finding increasing clinical applications as a non-invasive probe for vascular health. Further analysis extracting temporal delay information from the CVR response provide additional insight that reflect arterial transit time, blood redistribution, and vascular response speed. Untangling these factors can help better understand the (patho)physiology and improve diagnosis/prognosis associated with vascular impairments. Here, we use hypercapnic (HC) and hyperoxic (HO) challenges to gather insight about factors driving temporal delays between gray-matter (GM) and white-matter (WM). Blood Oxygen Level Dependent (BOLD) datasets were acquired at 7T in nine healthy subjects throughout BLOCK- and RAMP-HC paradigms. In a subset of seven participants, a combined HC+HO block, referred as the “BOOST” protocol, was also acquired. Tissue-based differences in Rapid Interpolation at Progressive Time Delays (RIPTiDe) were compared across stimulus to explore dynamic (BLOCK-HC) versus progressive (RAMP-HC) changes in CO₂, as well as the effect of bolus arrival time on CVR delays (BLOCK-HC versus BOOST). While GM delays were similar between the BLOCK- (21.80 ± 10.17 s) and RAMP-HC (24.29 ± 14.64 s), longer WM lag times were observed during the RAMP-HC (42.66 ± 17.79 s), compared to the BLOCK-HC (34.15 ± 10.72 s), suggesting that the progressive stimulus may predispose WM vasculature to longer delays due to the smaller arterial content of CO₂ delivered to WM tissues, which in turn, decreases intravascular CO₂ gradients modulating CO₂ diffusion into WM tissues. This was supported by a maintained ∼10 s offset in GM (11.66 ± 9.54 s) versus WM (21.40 ± 11.17 s) BOOST-delays with respect to the BLOCK-HC, suggesting that the vasoactive effect of CO₂ remains constant and that shortening of BOOST delays was be driven by blood arrival reflected through the non-vasodilatory HO contrast. These findings support that differences in temporal and magnitude aspects of CVR between vascular networks reflect a component of CO₂ sensitivity, in addition to redistribution and steal blood flow effects. Moreover, these results emphasize that the addition of a BOOST paradigm may provide clinical insights into whether vascular diseases causing changes in CVR do so by way of severe blood flow redistribution effects, alterations in vascular properties associated with CO₂ diffusion, or changes in blood arrival time.

Functional connectivity of the default mode, dorsal attention and fronto-parietal executive control networks in glial tumor patients

PURPOSE

Resting state functional magnetic resonance imaging (rsfMRI) is an emerging tool to explore the functional connectivity of different brain regions. We aimed to assess the disruption of functional connectivity of the Default Mode Network (DMN), Dorsal Attention Network(DAN) and Fronto-Parietal Network (FPN) in patients with glial tumors.

METHODS

rsfMRI data acquired on 3T-MR of treatment-naive glioma patients prospectively recruited (2015-2019) and matched controls from the 1000 functional-connectomes-project were analyzed using the CONN functional toolbox. Seed-Based Connectivity Analysis (SBCA) and Independent Component Analysis (ICA, with 10 to 100 components) were performed to study reliably the three networks of interest.

RESULTS

35 patients with gliomas (17 WHO grade I-II, 18 grade III-IV) and 70 controls were included. Global increased DMN connectivity was consistently found with SBCA and ICA in patients compared to controls (Cluster1: Precuneus, height: p < 10^-6; Cluster2: subcallosum; height: p < 10^-5). However, an area of decreased connectivity was found in the posterior corpus callosum, particularly in high-grade gliomas (height: p < 10^-5). The DAN demonstrated small areas of increased connectivity in frontal and occipital regions (height: p < 10^-6). For the FPN, increased connectivity was noted in the precuneus, posterior cingulate gyrus, and frontal cortex. No difference in the connectivity of the networks of interest was demonstrated between low- and high-grade gliomas, as well as when stratified by their IDH1-R132H (isocitrate dehydrogenase) mutation status.

CONCLUSION

Altered functional connectivity is reliably found with SBCA and ICA in the DMN, DAN, and FPN in glioma patients, possibly explained by decreased connectivity between the cerebral hemispheres across the corpus callosum due to disruption of the connections.

Effects of phase regression on high-resolution functional MRI of the primary visual cortex

High-resolution functional MRI studies have become a powerful tool to non-invasively probe the sub-millimeter functional organization of the human cortex. Advances in MR hardware, imaging techniques and sophisticated post-processing methods have allowed high resolution fMRI to be used in both the clinical and academic neurosciences. However, consensus within the community regarding the use of gradient echo (GE) or spin echo (SE) based acquisition remains largely divided. On one hand, GE provides a high temporal signal-to-noise ratio (tSNR) technique sensitive to both the macro- and micro-vascular signal while SE based methods are more specific to microvasculature but suffer from lower tSNR and specific absorption rate limitations, especially at high field and with short repetition times. Fortunately, the phase of the GE-EPI signal is sensitive to vessel size and this provides a potential avenue to reduce the macrovascular weighting of the signal (phase regression, Menon 2002). In order to determine the efficacy of this technique at high-resolution, phase regression was applied to GE-EPI timeseries and compared to SE-EPI to determine if GE-EPI's specificity to the microvascular compartment improved. To do this, functional data was collected from seven subjects on a neuro-optimized 7 T system at 800 μm isotropic resolution with both GE-EPI and SE-EPI while observing an 8 Hz contrast reversing checkerboard. Phase data from the GE-EPI was used to create a microvasculature-weighted time series (GE-EPI-PR). Anatomical imaging (MP2RAGE) was also collected to allow for surface segmentation so that the functional results could be projected onto a surface. A multi-echo gradient echo sequence was collected and used to identify venous vasculature. The GE-EPI-PR surface activation maps showed a high qualitative similarity with SE-EPI and also produced laminar activity profiles similar to SE-EPI. When the GE-EPI and GE-EPI-PR distributions were compared to SE-EPI it was shown that GE-EPI-PR had similar distribution characteristics to SE-EPI (p < 0.05) across the top 60% of cortex. Furthermore, it was shown that GE-EPI-PR has a higher contrast-to-noise ratio (0.5 ± 0.2, mean ± std. dev. across layers) than SE-EPI (0.27 ± 0.07) demonstrating the technique has higher sensitivity than SE-EPI. Taken together this evidence suggests phase regression is a useful method in low SNR studies such as high-resolution fMRI.

Dependence of stimulus-induced rotary saturation on the direction of target oscillating magnetic fields: A phantom and simulation study

Several noninvasive techniques for the direct measurement of the neuronal activity using magnetic resonance imaging (MRI) have recently been reported. As a promising candidate, we focus on a spin-lock MRI sequence (i.e., stimulus-induced rotary saturation (SIRS)) directly measuring a tiny oscillating magnetic field. Previous phantom studies on SIRS have applied the target oscillating magnetic field parallel to the direction of the static magnetic field B₀. However, in practice, the neuromagnetic fields are not always aligned in the same direction as in such a condition. This study investigates the MR signal changes during SIRS when the target magnetic field direction is not the same as that of the B₀ field through both phantom experiments and Bloch simulations. The experimental results indicate that only the target magnetic field component along the B₀ field affects the signal change, indicating that SIRS has partial sensitivity, even if the target magnetic fields are tilted from the B₀ field. Furthermore, the simulation results show good agreements with the experimental results. These results clarify the sensitivity direction of SIRS-based fMRI and lead to the possibility that the direction of the generated neuromagnetic fields can be estimated, such that we can separate directional information from the other information contained in neuromagnetic fields (e.g., phase information).

Functional Interplay Between Posterior Parietal Cortex and Hippocampus During Detection of Memory Targets and Non-targets

Posterior parietal cortex is frequently activated during episodic memory retrieval but its role during retrieval and its interactions with the hippocampus are not yet clear. In this fMRI study, we investigated the neural bases of recognition memory when study repetitions and retrieval goals were manipulated. During encoding participants studied words either once or three times, and during retrieval they were rewarded more to detect either studied words or new words. We found that (1) dorsal parietal cortex (DPC) was more engaged during detection of items studied once compared to three times, whereas regions in the ventral parietal cortex (VPC) responded more to items studied multiple times; (2) DPC, within a network of brain regions functionally connected to the anterior hippocampus, responded more to items consistent with retrieval goals (associated with high reward); (3) VPC, within a network of brain regions functionally connected to the posterior hippocampus, responded more to items not aligned with retrieval goals (i.e., unexpected). These findings support the hypothesis that DPC and VPC regions contribute differentially to top-down vs. bottom-up attention to memory. Moreover, they reveal a dissociation in the functional profile of the anterior and posterior hippocampi.

Evaluation of a 16-channel transceiver loop + dipole antenna array for human head imaging at 10.5 tesla

We evaluated a 16-channel loop + dipole (LD) transceiver antenna array with improved specific absorption rate (SAR) efficiency for 10.5 Tesla (T) human head imaging apsplications. Three different array designs with equal inner dimensions were considered: an 8-channel dipole antenna, an 8-channel loop, and a 16-channel LD antenna arrays. Signal-to-noise ratio (SNR) and B1 + efficiency (in units of μT per √W) were simulated and measured in 10.5 T magnetic resonance imaging (MRI) experiments. For the safety validation, 10 g SAR and SAR efficiency (defined as the B1 + over √ (peak 10 g SAR)) were calculated through simulation. Finally, high resolution porcine brain images were acquired with the 16-channel LD antenna array, including a fast turbo-spin echo (TSE) sequence incorporating B1 shimming techniques. Both the simulation and experiments demonstrated that the combined 16-channel LD antenna array showed similar B1 + efficiency compared to the 8-channel dipole antenna and the 8-channel loop arrays in a circular polarized (CP) mode. In a central 2 mm × 2 mm region of the phantom, however, the 16-channel LD antenna array showed an improvement in peak 10 g SAR of 27.5 % and 32.5 % over the 8-channel dipole antenna and the 8-channel loop arrays, respectively. We conclude that the proposed 16-channel head LD antenna array design is capable of achieving ~7% higher SAR efficiency at 10.5 T compared to either the 8-channel loop-only or the 8-channel dipole-only antenna arrays of the same dimensions.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

Blood oxygenation level dependent magnetic resonance imaging and diffusion weighted MRI imaging for benign and malignant breast cancer discrimination

PURPOSE

The purpose of this study is to assess Blood oxygenation level dependent Magnetic Resonance Imaging (BOLD-MRI) and Diffusion Weighted Magnetic Resonance Imaging (DW-MRI) in the differentiation of benign and malignant breast lesions.

METHODS

Fifty-nine breast lesions (26 benign and 33 malignant lesions) pathologically proven in 59 patients were included in this retrospective study. As BOLD parameters were estimated basal signal S₀ and the relaxation rate R2*, diffusion and perfusion parameters were derived by DWI (pseudo-diffusion coefficient (Dp), perfusion fraction (fp) and tissue diffusivity (Dt)). Wilcoxon-Mann-Whitney U test and Receiver operating characteristic (ROC) analyses were calculated and area under ROC curve (AUC) was obtained. Moreover, pattern recognition approaches (linear discrimination analysis (LDA), support vector machine, k-nearest neighbours, decision tree) with least absolute shrinkage and selection operator (LASSO) method and leave one out cross validation approach were considered.

RESULTS

A significant discrimination was obtained by the standard deviation value of S0, as BOLD parameter, that reached an AUC of 0.76 with a sensitivity of 65%, a specificity of 85% and an accuracy of 76%. No significant discrimination was obtained considering diffusion and perfusion parameters. Considering LASSO results, the features to use as predictors were all extracted parameters except that the mean value of R2* and the best result was obtained by a LDA that obtained an AUC = 0.83, with a sensitivity of 88%, a specificity of 77% and an accuracy of 83%.

CONCLUSIONS

Good performance to discriminate benign and malignant lesions could be obtained using BOLD and DWI derived parameters with a LDA classification approach. However, these findings should be proven on larger and several dataset with different MR scanners.

Measurement of CMRO2 and its relationship with CBF in hypoxia with an extended calibrated BOLD method

Cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO₂) are physiological parameters that not only reflect brain health and disease but also jointly contribute to blood oxygen level-dependent (BOLD) signals. Nevertheless, unsolved issues remain concerning the CBF-CMRO₂ relationship in the working brain under various oxygen conditions. In particular, the CMRO₂ responses to functional tasks in hypoxia are less studied. We extended the calibrated BOLD model to incorporate CMRO₂ measurements in hypoxia. The extended model, which was cross-validated with a multicompartment BOLD model, considers the influences of the reduced arterial saturation level and increased baseline cerebral blood volume (CBV) and deoxyhemoglobin concentration on the changes of BOLD signals in hypoxia. By implementing a pulse sequence to simultaneously acquire the CBV-, CBF- and BOLD-weighted signals, we investigated the effects of mild hypoxia on the CBF and CMRO₂ responses to graded visual stimuli. Compared with normoxia, mild hypoxia caused significant alterations in both the amplitude and the trend of the CMRO₂ responses but did not impact the corresponding CBF responses. Our observations suggested that the flow-metabolism coupling strategies in the brain during mild hypoxia were different from those during normoxia.

High spatiotemporal vessel-specific hemodynamic mapping with multi-echo single-vessel fMRI

High-resolution fMRI enables noninvasive mapping of the hemodynamic responses from individual penetrating vessels in animal brains. Here, a 2D multi-echo single-vessel fMRI (MESV-fMRI) method has been developed to map the fMRI signal from arterioles and venules with a 100 ms sampling rate at multiple echo times (TE, 3-30 ms) and short acquisition windows (<1 ms). The T₂*-weighted signal shows the increased extravascular effect on venule voxels as a function of TE. In contrast, the arteriole voxels show an increased fMRI signal with earlier onset than venules voxels at the short TE (3 ms) with increased blood inflow and volume effects. MESV-fMRI enables vessel-specific T2* mapping and presents T2*-based fMRI time courses with higher contrast-to-noise ratios (CNRs) than the T2*-weighted fMRI signal at a given TE. The vessel-specific T2* mapping also allows semi-quantitative estimation of the oxygen saturation levels (Y) and their changes (ΔY) at a given blood volume fraction upon neuronal activation. The MESV-fMRI method enables vessel-specific T2* measurements with high spatiotemporal resolution for better modeling of the fMRI signal based on the hemodynamic parameters.

Venous cerebral blood volume mapping in the whole brain using venous-spin-labeled 3D turbo spin echo

PURPOSE

Venous cerebral blood volume (CBVv ) is a major contributor to BOLD contrast, and therefore is an important parameter for understanding the underlying mechanism. Here, we propose a velocity-selective venous spin labeling (VS-VSL)-prepared 3D turbo spin echo pulse sequence for whole-brain baseline CBVv mapping.

METHODS

Unlike previous CBVv measurement techniques that exploit the interrelationship between BOLD signals and CBVv , in the proposed VS-VSL technique both arterial blood and cerebrospinal fluid (CSF) signals were suppressed before the VS pulse train for exclusive labeling of venous blood, while a single-slab 3D turbo spin echo readout was used because of its relative immunity to magnetic field variations. Furthermore, two approximations were made to the VS-VSL signal model for simplified derivation of CBVv . In vivo studies were performed at 3T field strength in 8 healthy subjects. The performance of the proposed VS-VSL method in baseline CBVv estimation was first evaluated in comparison to the existing, hyperoxia-based method. Then, data were also acquired using VS-VSL under hypercapnic and hyperoxic gas breathing challenges for further validation of the technique.

RESULTS

The proposed technique yielded physiologically plausible baseline CBVv values, and when compared with the hyperoxia-based method, showed no statistical difference. Furthermore, data acquired using VS-VSL yielded average CBVv of 2.89%/1.78%, 3.71%/2.29%, and 2.88%/1.76% for baseline, hypercapnia, and hyperoxia, respectively, in gray/white matter regions. As expected, hyperoxia had negligible effect (P > .8), whereas hypercapnia demonstrated vasodilation (P << .01).

CONCLUSION

Upon further validation of the quantification model, the method is expected to have merit for 3D CBVv measurements across the entire brain.

Optimizing Visual Cortex Parameterization with Error-Tolerant Teichmüller Map in Retinotopic Mapping

The mapping between the visual input on the retina to the cortical surface, i.e., retinotopic mapping, is an important topic in vision science and neuroscience. Human retinotopic mapping can be revealed by analyzing cortex functional magnetic resonance imaging (fMRI) signals when the subject is under specific visual stimuli. Conventional methods process, smooth, and analyze the retinotopic mapping based on the parametrization of the (partial) cortical surface. However, the retinotopic maps generated by this approach frequently contradict neuropsychology results. To address this problem, we propose an integrated approach that parameterizes the cortical surface, such that the parametric coordinates linearly relates the visual coordinate. The proposed method helps the smoothing of noisy retinotopic maps and obtains neurophysiological insights in human vision systems. One key element of the approach is the Error-Tolerant Teichmüller Map, which uniforms the angle distortion and maximizes the alignments to self-contradicting landmarks. We validated our overall approach with synthetic and real retinotopic mapping datasets. The experimental results show the proposed approach is superior in accuracy and compatibility. Although we focus on retinotopic mapping, the proposed framework is general and can be applied to process other human sensory maps.

Use of perfusional CBCT imaging for intraprocedural evaluation of endovascular treatment in patients with diabetic foot: a concept paper

Diabetes mellitus (DM) is one of the most common metabolic diseases worldwide; its global burden has increased rapidly over the past decade, enough to be considered a public health emergency in many countries. Diabetic foot disease and, particularly diabetic foot ulceration, is the major complication of DM: through a skin damage of the foot, with a loss of epithelial tissue, it can deepen to muscles and bones and lead to the amputation of the lower limbs. Peripheral arterial disease (PAD) in patients with diabetes, manifests like a diffuse macroangiopathic multi-segmental involvement of the lower limb vessels, also connected to a damage of collateral circulation; it may also display characteristic microaneurysms and tortuosity in distal arteries. As validation method, Bold-MRI is used. The diabetic foot should be handled with a multidisciplinary team approach, as its management requires systemic and localized treatments, pain control, monitoring of cardiovascular risk factors and other comorbidities. CBCT is an emerging medical imaging technique with the original feature of divergent radiation, forming a cone, in contrast with the spiral slicing of conventional CT, and has become increasingly important in treatment planning and diagnosis: from small anatomical areas, such as implantology, to the world of interventional radiology, with a wide range of applications: as guidance for biopsies or ablation treatments. The aim of this project is to evaluate the usefulness of perfusion CBCT imaging, obtained during endovascular revascularization, for intraprocedural evaluation of endovascular treatment in patients with diabetic foot. (www.actabiomedica.it).

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.

Cerebrovascular Reactivity Assays Collateral Function in Carotid Stenosis

In patients with carotid artery stenosis (CAS), the risk of stroke, its severity, and response to revascularization are strongly related to the availability of collateral blood flow. Unfortunately, there is poor agreement between observers in assessing collateral flow using flow-based imaging. We used changes in blood-oxygen-level-dependent (BOLD) MRI as a surrogate of changes in regional cerebral blood flow in response to a hypercapnic stimulus [i.e., cerebrovascular reactivity (CVR)] as indicating flow reserve ipsilateral to CAS. We hypothesized that some patients with hemodynamically significant CAS develop functional collateral flow as indicated by normalization of ipsilateral CVR. We identified 55 patients in our CVR database with various degrees of CAS assessed by angiography and classed them as <50% stenosis, 50–69% stenosis, 70–90% stenosis, >90% stenosis, and full occlusion. CVR was measured as the change in BOLD signal in response to changes in end-tidal partial pressure of CO₂ (Δ BOLD/Δ PETCO₂) and normalized voxel-wise relative to the mean and standard deviation of the CVR in the corresponding voxels of an atlas of 46 healthy controls (CVR z scores). CVR and z scores were then averaged over gray matter (GM) and white matter (WM) on each side of the middle cerebral artery (MCA) territory. As hypothesized, CVR varied for each severity of CAS. Ipsilateral MCA territory CVR was less than normal in each class, including that with <50% stenosis (Student t-test, two-tailed; p = 0.0014 for GM and p = 0.030 for WM), with a trend of decreasing average CVR with increasing stenosis. Remarkably, the considerable individual variability in MCA CVR included some patients with normal CVR in each class – including that with complete occlusion. We conclude that, in general, CAS depresses downstream vascular reserve, but the extent of collateralization is highly variable and not predictable from the degree of stenosis, including both <50% stenosis and complete occlusion. CVR may be the more reliable marker for recruitable collateral blood flow than degree of CAS.

Dynamics of the cerebral blood flow response to brief neural activity in human visual cortex

The blood oxygen-level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal depends on an interplay of cerebral blood flow (CBF), oxygen metabolism, and cerebral blood volume. Despite wide usage of BOLD fMRI, it is not clear how these physiological components create the BOLD signal. Here, baseline CBF and its dynamics evoked by a brief stimulus (2 s) in human visual cortex were measured at 3T. We found a stereotypical CBF response: immediate increase, rising to a peak a few second after the stimulus, followed by a significant undershoot. The BOLD hemodynamic response function (HRF) was also measured in the same session. Strong correlations between HRF and CBF peak responses indicate that the flow responses evoked by neural activation in nearby gray matter drive the early HRF. Remarkably, peak CBF and HRF were also strongly modulated by baseline perfusion. The CBF undershoot was reliable and significantly correlated with the HRF undershoot. However, late-time dynamics of the HRF and CBF suggest that oxygen metabolism can also contribute to the HRF undershoot. Combined measurement of the CBF and HRF for brief neural activation is a useful tool to understand the temporal dynamics of neurovascular and neurometabolic coupling.