Model-based physiological noise removal in fast fMRI

Advances in functional and structural imaging of the brainstem: implications for disease

PURPOSE OF REVIEW

The brainstem's complex anatomy and relatively small size means that structural and functional assessment of this structure is done less frequently compared to other brain areas. However, recent years have seen substantial progress in brainstem imaging, enabling more detailed investigations into its structure and function, as well as its role in neuropathology.

RECENT FINDINGS

Advancements in ultrahigh field MRI technology have allowed for unprecedented spatial resolution in brainstem imaging, facilitating the new creation of detailed brainstem-specific atlases. Methodological improvements have significantly enhanced the accuracy of physiological (cardiac and respiratory) noise correction within brainstem imaging studies. These technological and methodological advancements have allowed for in-depth analyses of the brainstem's anatomy, including quantitative assessments and examinations of structural connectivity within both gray and white matter. Furthermore, functional studies, including assessments of activation patterns and functional connectivity, have revealed the brainstem's roles in both specialized functions and broader neural integration. Notably, these investigations have identified alterations in brainstem structure and function associated with various neurological disorders.

SUMMARY

The aforementioned developments have allowed for a greater appreciation of the importance of the brainstem in the wider context of neuroscience and clinical neurology.

Comparing data-driven physiological denoising approaches for resting-state fMRI: implications for the study of aging

Introduction

Physiological nuisance contributions by cardiac and respiratory signals have a significant impact on resting-state fMRI data quality. As these physiological signals are often not recorded, data-driven denoising methods are commonly used to estimate and remove physiological noise from fMRI data. To investigate the efficacy of these denoising methods, one of the first steps is to accurately capture the cardiac and respiratory signals, which requires acquiring fMRI data with high temporal resolution.

Methods

In this study, we used such high-temporal resolution fMRI data to evaluate the effectiveness of several data-driven denoising methods, including global-signal regression (GSR), white matter and cerebrospinal fluid regression (WM-CSF), anatomical (aCompCor) and temporal CompCor (tCompCor), ICA-AROMA. Our analysis focused on the consequence of changes in low-frequency, cardiac and respiratory signal power, as well as age-related differences in terms of functional connectivity (fcMRI).

Results

Our results confirm that the ICA-AROMA and GSR removed the most physiological noise but also more low-frequency signals. These methods are also associated with substantially lower age-related fcMRI differences. On the other hand, aCompCor and tCompCor appear to be better at removing high-frequency physiological signals but not low-frequency signal power. These methods are also associated with relatively higher age-related fcMRI differences, whether driven by neuronal signal or residual artifact. These results were reproduced in data downsampled to represent conventional fMRI sampling frequency. Lastly, methods differ in performance depending on the age group.

Discussion

While this study cautions direct comparisons of fcMRI results based on different denoising methods in the study of aging, it also enhances the understanding of different denoising methods in broader fcMRI applications.

Distributed fMRI dynamics predict distinct EEG rhythms across sleep and wakefulness

The brain exhibits rich oscillatory dynamics that vary across tasks and states, such as the EEG oscillations that define sleep. These oscillations play critical roles in cognition and arousal, but the brainwide mechanisms underlying them are not yet described. Using simultaneous EEG and fast fMRI in subjects drifting between sleep and wakefulness, we developed a machine learning approach to investigate which brainwide fMRI dynamics predict alpha (8-12 Hz) and delta (1-4 Hz) rhythms. We predicted moment-by-moment EEG power from fMRI activity in held-out subjects, and found that information about alpha power was represented by a remarkably small set of regions, segregated in two distinct networks linked to arousal and visual systems. Conversely, delta rhythms were diffusely represented on a large spatial scale across the cortex. These results identify distributed networks that predict delta and alpha rhythms, and establish a computational framework for investigating fMRI brainwide dynamics underlying EEG oscillations.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here, we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

EEG-LLAMAS: A low-latency neurofeedback platform for artifact reduction in EEG-fMRI

Simultaneous EEG-fMRI is a powerful multimodal technique for imaging the brain, but its use in neurofeedback experiments has been limited by EEG noise caused by the MRI environment. Neurofeedback studies typically require analysis of EEG in real time, but EEG acquired inside the scanner is heavily contaminated with ballistocardiogram (BCG) artifact, a high-amplitude artifact locked to the cardiac cycle. Although techniques for removing BCG artifacts do exist, they are either not suited to real-time, low-latency applications, such as neurofeedback, or have limited efficacy. We propose and validate a new open-source artifact removal software called EEG-LLAMAS (Low Latency Artifact Mitigation Acquisition Software), which adapts and advances existing artifact removal techniques for low-latency experiments. We first used simulations to validate LLAMAS in data with known ground truth. We found that LLAMAS performed better than the best publicly-available real-time BCG removal technique, optimal basis sets (OBS), in terms of its ability to recover EEG waveforms, power spectra, and slow wave phase. To determine whether LLAMAS would be effective in practice, we then used it to conduct real-time EEG-fMRI recordings in healthy adults, using a steady state visual evoked potential (SSVEP) task. We found that LLAMAS was able to recover the SSVEP in real time, and recovered the power spectra collected outside the scanner better than OBS. We also measured the latency of LLAMAS during live recordings, and found that it introduced a lag of less than 50 ms on average. The low latency of LLAMAS, coupled with its improved artifact reduction, can thus be effectively used for EEG-fMRI neurofeedback. A limitation of the method is its use of a reference layer, a piece of EEG equipment which is not commercially available, but can be assembled in-house. This platform enables closed-loop experiments which previously would have been prohibitively difficult, such as those that target short-duration EEG events, and is shared openly with the neuroscience community.

Less is more: balancing noise reduction and data retention in fMRI with data-driven scrubbing

Functional MRI (fMRI) data may be contaminated by artifacts arising from a myriad of sources, including subject head motion, respiration, heartbeat, scanner drift, and thermal noise. These artifacts cause deviations from common distributional assumptions, introduce spatial and temporal outliers, and reduce the signal-to-noise ratio of the data-all of which can have negative consequences for the accuracy and power of downstream statistical analysis. Scrubbing is a technique for excluding fMRI volumes thought to be contaminated by artifacts and generally comes in two flavors. Motion scrubbing based on subject head motion-derived measures is popular but suffers from a number of drawbacks, among them the need to choose a threshold, a lack of generalizability to multiband acquisitions, and high rates of censoring of individual volumes and entire subjects. Alternatively, data-driven scrubbing methods like DVARS are based on observed noise in the processed fMRI timeseries and may avoid some of these issues. Here we propose "projection scrubbing", a novel data-driven scrubbing method based on a statistical outlier detection framework and strategic dimension reduction, including independent component analysis (ICA), to isolate artifactual variation. We undertake a comprehensive comparison of motion scrubbing with data-driven projection scrubbing and DVARS. We argue that an appropriate metric for the success of scrubbing is maximal data retention subject to reasonable performance on typical benchmarks such as the validity, reliability, and identifiability of functional connectivity. We find that stringent motion scrubbing yields worsened validity, worsened reliability, and produced small improvements to fingerprinting. Meanwhile, data-driven scrubbing methods tend to yield greater improvements to fingerprinting while not generally worsening validity or reliability. Importantly, however, data-driven scrubbing excludes a fraction of the number of volumes or entire sessions compared to motion scrubbing. The ability of data-driven fMRI scrubbing to improve data retention without negatively impacting the quality of downstream analysis has major implications for sample sizes in population neuroscience research.

Improving robustness of 3D multi-shot EPI by structured low-rank reconstruction of segmented CAIPI sampling for fMRI at 7T

Three-dimensional (3D) encoding methods are increasingly being explored as alternatives to two-dimensional (2D) multi-slice acquisitions in fMRI, particularly in cases where high isotropic resolution is needed. 3D multi-shot EPI acquisition, as the workhorse of 3D fMRI imaging, is susceptible to physiological fluctuations which can induce inter-shot phase variations, and thus reducing the achievable tSNR, negating some of the benefit of 3D encoding. This issue can be particularly problematic at ultra-high fields like 7T, which have more severe off-resonance effects. In this work, we aim to improve the temporal stability of 3D multi-shot EPI at 7T by improving its robustness to inter-shot phase variations. We presented a 3D segmented CAIPI sampling trajectory ("seg-CAIPI") and an improved reconstruction method based on Hankel structured low-rank matrix recovery. Simulation and in-vivo results demonstrate that the combination of the seg-CAIPI sampling scheme and the proposed structured low-rank reconstruction is a promising way to effectively reduce the unwanted temporal variance induced by inter-shot physiological fluctuations, and thus improve the robustness of 3D multi-shot EPI for fMRI.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, as differences can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR)

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

WHOCARES: WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions

Objective. To spatio-temporally resolve cardiac signals in functional magnetic resonance imaging (fMRI) time-series of the human brain using neither external physiological measurements nor ad hoc modelling assumptions.Approach. Cardiac pulsation is a physiological confound of fMRI time-series that introduces spurious signal fluctuations in proximity to blood vessels. fMRI alone is not sufficiently fast to resolve cardiac pulsation. Depending on the ratio between the instantaneous heart-rate and the acquisition sampling frequency (1/TR, with TR being the repetition time), the cardiac signal may alias into the frequency band of neural activation so that its removal through spectral filtering techniques is generally not possible. In this paper, we show that it is feasible to temporally and spatially resolve cardiac signals throughout the brain even when cardiac aliasing occurs by combining fMRI hyper-sampling with simultaneous multislice (SMS) imaging. The technique, which we name WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions (WHOCARES), was developed on 695 healthy subjects selected from the Human Connectome Project and its performance validated against the RETROICOR, HAPPY and the pulse oxymeter signal regression methods.Main results.WHOCARES is capable of retrieving voxel-wise cardiac signal regressors. This is achieved without employing external physiological recordings nor through ad hoc modelling assumptions. The performance of WHOCARES was, on average, superior to RETROICOR, HAPPY and the pulse oxymeter regression methods.Significance.WHOCARES holds basis for the reliable mapping of cardiac activity in fMRI time-series. WHOCARES can be employed for the retrospective removal of cardiac noise in publicly available fMRI datasets where physiological recordings are not available. WHOCARES is freely available athttps://github.com/gferrazzi/WHOCARES.

The influence of the subcortex and brain stem on overeating: How advances in functional neuroimaging can be applied to expand neurobiological models to beyond the cortex

Functional neuroimaging has become a widely used tool in obesity and eating disorder research to explore the alterations in neurobiology that underlie overeating and binge eating behaviors. Current and traditional neurobiological models underscore the importance of impairments in brain systems supporting reward, cognitive control, attention, and emotion regulation as primary drivers for overeating. Due to the technical limitations of standard field strength functional magnetic resonance imaging (fMRI) scanners, human neuroimaging research to date has focused largely on cortical and basal ganglia effects on appetitive behaviors. The present review draws on animal and human research to highlight how neural signaling encoding energy regulation, reward-learning, and habit formation converge on hypothalamic, brainstem, thalamic, and striatal regions to contribute to overeating in humans. We also consider the role of regions such as the mediodorsal thalamus, ventral striatum, lateral hypothalamus and locus coeruleus in supporting habit formation, inhibitory control of food craving, and attentional biases. Through these discussions, we present proposals on how the neurobiology underlying these processes could be examined using functional neuroimaging and highlight how ultra-high field 7-Tesla (7 T) fMRI may be leveraged to elucidate the potential functional alterations in subcortical networks. Focus is given to how interactions of these regions with peripheral endocannabinoids and neuropeptides, such as orexin, could be explored. Technical and methodological aspects regarding the use of ultra-high field 7 T fMRI to study eating behaviors are also reviewed.

Effects of the Prolong Life With Nine Turn Method (Yan Nian Jiu Zhuan) Qigong on Brain Functional Changes in Patients With Chronic Fatigue Syndrome in Terms of Fatigue and Quality of Life

Background

Chronic fatigue syndrome (CFS) is characterized by persistent fatigue, which often leads to physical and psychological damage. The Prolong Life with Nine Turn method (PLWNT) Qigong is considered as one of the complementary treatments for improving symptoms in patients with CFS. In this study, we used functional magnetic resonance imaging (fMRI) to explore the effects of PLWNT intervention on the subjects with CFS.

Methods

Thirty four CFS patients were randomly divided into PLWNT group and cognitive behavioral therapy (CBT) group. Both groups were taught by a highly qualified professor at the Shanghai University of Traditional Chinese Medicine once a week and were supervised online during the remaining 6 days at home, over 12 consecutive weeks. We calculated the regional rs-fMRI index amplitude of low-frequency fluctuations (ALFF) for all subjects. To study the changes of the brain network, we used the brain regions with significant differences in ALFF as the regions of interest for whole-brain functional connectivity (FC) analysis. The Multi-dimensional Fatigue Inventory 20 (MFI-20) and Short Form 36-item Health Survey (SF-36) were used for clinical symptom assessment to explore the possible correlation between the rs-fMRI indicators and clinical variations.

Results

The ALFF values of the right superior frontal gyrus (SFG), and left median cingulate gyrus (DCG) were increased, whereas those of the left middle occipital gyrus (OG), right middle OG and left middle temporal gyrus (MTG) were decreased in PLWNT group. The FC values between the DCG and middle temporal gyrus (MTG), and those between the left OG and the right OG were enhanced. In addition, the SF-36 were positively with the left OG (r = 0.524), SFG (r = 0.517), and DCG (r = 0.533), MFI-20 were negatively with the SFG (r = −0.542) and DCG (r = −0.578). These results were all corrected by FWE (voxel level p < 0.001, cluster level p < 0.05).

Conclusion

CFS patients have abnormal regional spontaneous neuronal activity and abnormal functional connections between regions after PLWNT intervention. PLWNT can relieve the fatigue symptoms of CFS patients and improve their quality of life. The study was registered in the American Clinical Trial Registry (12/04/2018). Registration Number is NCT03496961.

Functional MRI in major depressive disorder: A review of findings, limitations, and future prospects

Objective diagnosis and prognosis in major depressive disorder (MDD) remains a challenge due to the absence of biomarkers based on physiological parameters or medical tests. Numerous studies have been conducted to identify functional magnetic resonance imaging‐based biomarkers of depression that either objectively differentiate patients with depression from healthy subjects, predict personalized treatment outcome, or characterize biological subtypes of depression. While there are some findings of consistent functional biomarkers, there is still lack of robust data acquisition and analysis methodology. According to current findings, primarily, the anterior cingulate cortex, prefrontal cortex, and default mode network play a crucial role in MDD. Yet, there are also less consistent results and the involvement of other regions or networks remains ambiguous. We further discuss image acquisition, processing, and analysis limitations that might underlie these inconsistencies. Finally, the current review aims to address and discuss possible remedies and future opportunities that could improve the search for consistent functional imaging biomarkers of depression. Novel acquisition techniques, such as multiband and multiecho imaging, and neural network‐based cleaning approaches can enhance the signal quality in limbic and frontal regions. More comprehensive analyses, such as directed or dynamic functional features or the identification of biological depression subtypes, can improve objective diagnosis or treatment outcome prediction and mitigate the heterogeneity of MDD. Overall, these improvements in functional MRI imaging techniques, processing, and analysis could advance the search for biomarkers and ultimately aid patients with MDD and their treatment course.

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.

Frequency-specific coactivation patterns in resting-state and their alterations in schizophrenia: An fMRI study

The resting‐state human brain is a dynamic system that shows frequency‐dependent characteristics. Recent studies demonstrate that coactivation pattern (CAP) analysis can identify recurring brain states with similar coactivation configurations. However, it is unclear whether and how CAPs depend on the frequency bands. The current study investigated the spatial and temporal characteristics of CAPs in the four frequency sub‐bands from slow‐5 (0.01–0.027 Hz), slow‐4 (0.027–0.073 Hz), slow‐3 (0.073–0.198 Hz), to slow‐2 (0.198–0.25 Hz), in addition to the typical low‐frequency range (0.01–0.08 Hz). In the healthy subjects, six CAP states were obtained at each frequency band in line with our prior study. Similar spatial patterns with the typical range were observed in slow‐5, 4, and 3, but not in slow‐2. While the frequency increased, all CAP states displayed shorter persistence, which caused more between‐state transitions. Specifically, from slow‐5 to slow‐4, the coactivation not only changed significantly in distributed cortical networks, but also increased in the basal ganglia as well as the amygdala. Schizophrenia patients showed significant alteration in the persistence of CAPs of slow‐5. Using leave‐one‐pair‐out, hold‐out and resampling validations, the highest classification accuracy (84%) was achieved by slow‐4 among different frequency bands. In conclusion, our findings provide novel information about spatial and temporal characteristics of CAP states at different frequency bands, which contributes to a better understanding of the frequency aspect of biomarkers for schizophrenia and other disorders.

Deep Attentive Spatio-Temporal Feature Learning for Automatic Resting-State fMRI Denoising

Resting-state functional magnetic resonance imaging (rs-fMRI) is a non-invasive functional neuroimaging modality that has been widely used to investigate functional connectomes in the brain. Since noise and artifacts generated by non-neuronal physiological activities are predominant in raw rs-fMRI data, effective noise removal is one of the most important preprocessing steps prior to any subsequent analysis. For rs-fMRI denoising, a common trend is to decompose rs-fMRI data into multiple components and then regress out noise-related components. Therefore, various machine learning techniques have been used in such analyses with predefined procedures and manually engineered features. However, the lack of a universal definition of a noise-related source or artifact complicates manual feature engineering. Manual feature selection can result in the failure to capture unknown types of noise. Furthermore, the possibility that the hand-crafted features will only work for the broader population (e.g., healthy adults) but not for "outliers" (e.g., infants or subjects that belong to a disease cohort) is quite high. In practice, we have limited knowledge of which features should be extracted; thus, multi-classifier assembly must be implemented to improve performance, although this process is quite time-consuming. However, in real rs-fMRI applications, fast and accurate automatic identification of noise-related components on different datasets is critical. To solve this problem, we propose a novel, automatic, and end-to-end deep learning framework dedicated to noise-related component identification via a faster and more effective multi-layer feature extraction strategy that learns deeply embedded spatio-temporal features of the components. In this study, we achieved remarkable performance on various rs-fMRI datasets, including multiple adult rs-fMRI datasets from different rs-fMRI studies and an infant rs-fMRI dataset, which is quite heterogeneous and differs from that of adults. Our proposed framework also dramatically increases the noise detection speed owing to its inherent ability for deep learning (< 1s for single-component classification). It can be easily integrated into any preprocessing pipeline, even those that do not use standard procedures but depend on alternative toolboxes.

Cardiovascular Pulsatility Increases in Visual Cortex Before Blood Oxygen Level Dependent Response During Stimulus

The physiological pulsations that drive tissue fluid homeostasis are not well characterized during brain activation. Therefore, we used fast magnetic resonance encephalography (MREG) fMRI to measure full band (0–5 Hz) blood oxygen level-dependent (BOLD_FB) signals during a dynamic visual task in 23 subjects. This revealed brain activity in the very low frequency (BOLD_VLF) as well as in cardiac and respiratory bands. The cardiovascular hemodynamic envelope (CHe) signal correlated significantly with the visual BOLD_VLF response, considered as an independent signal source in the V1-V2 visual cortices. The CHe preceded the canonical BOLD_VLF response by an average of 1.3 (± 2.2) s. Physiologically, the observed CHe signal could mark increased regional cardiovascular pulsatility following vasodilation.

Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI

Despite extensive studies detecting laminar functional magnetic resonance imaging (fMRI) signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across cortical regions remain to be fully investigated in both evoked and resting conditions at different brain states. Here, we developed a multislice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 μm) and temporal resolution (100 ms) in anesthetized rats. Across different trials, we detected either laminar-specific positive or negative blood-oxygen-level-dependent (BOLD) responses in the surrounding cortical region adjacent to the most activated cortex under the evoked condition. Specifically, in contrast to typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for trials with negative BOLD in adjacent regions, which indicated brain state-dependent laminar-fMRI responses based on corticocortical interaction. Also, in resting-state (rs-) fMRI study, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from caudal to rostral slices. In summary, our study provided a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity in evoked and resting-state MS-LS.

Brain tumor detection and classification using machine learning: a comprehensive survey

Brain tumor occurs owing to uncontrolled and rapid growth of cells. If not treated at an initial phase, it may lead to death. Despite many significant efforts and promising outcomes in this domain, accurate segmentation and classification remain a challenging task. A major challenge for brain tumor detection arises from the variations in tumor location, shape, and size. The objective of this survey is to deliver a comprehensive literature on brain tumor detection through magnetic resonance imaging to help the researchers. This survey covered the anatomy of brain tumors, publicly available datasets, enhancement techniques, segmentation, feature extraction, classification, and deep learning, transfer learning and quantum machine learning for brain tumors analysis. Finally, this survey provides all important literature for the detection of brain tumors with their advantages, limitations, developments, and future trends.

Imaging the temporal dynamics of brain states with highly sampled fMRI

The spontaneous dynamics of the brain modulate its function from moment to moment, shaping neural computation and cognition. Functional MRI (fMRI), while classically used as a tool for spatial localization, is increasingly being used to identify the temporal dynamics of brain activity. fMRI analyses focused on the temporal domain have revealed important new information about the dynamics underlying states such as arousal, attention, and sleep. Dense temporal sampling – either by using fast fMRI acquisition, or multiple repeated scan sessions within individuals – can further enrich the information present in these studies. This review focuses on recent developments in using fMRI to identify dynamics across brain states, particularly vigilance and sleep states, and the potential for highly temporally sampled fMRI to answer these questions.

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.

Statistical power or more precise insights into neuro-temporal dynamics? Assessing the benefits of rapid temporal sampling in fMRI

Functional magnetic resonance imaging (fMRI), a non-invasive and widely used human neuroimaging method, is most known for its spatial precision. However, there is a growing interest in its temporal sensitivity. This is despite the temporal blurring of neuronal events by the blood oxygen level dependent (BOLD) signal, the peak of which lags neuronal firing by 4-6 seconds. Given this, the goal of this review is to answer a seemingly simple question - "What are the benefits of increased temporal sampling for fMRI?". To answer this, we have combined fMRI data collected at multiple temporal scales, from 323 to 1000 milliseconds, with a review of both historical and contemporary temporal literature. After a brief discussion of technological developments that have rekindled interest in temporal research, we next consider the potential statistical and methodological benefits. Most importantly, we explore how fast fMRI can uncover previously unobserved neuro-temporal dynamics - effects that are entirely missed when sampling at conventional 1 to 2 second rates. With the intrinsic link between space and time in fMRI, this temporal renaissance also delivers improvements in spatial precision. Far from producing only statistical gains, the array of benefits suggest that the continued temporal work is worth the effort.

Mapping dependencies of BOLD signal change to end-tidal CO2: linear and nonlinear modeling, and effect of physiological noise correction

BACKGROUND

Disentangling physiological noise and signal of interest is a major issue when evaluating BOLD-signal changes in response to breath holding. Currently-adopted approaches for retrospective noise correction are general-purpose, and have non-negligible effects in studies on hypercapnic challenges.

NEW METHOD

We provide a novel approach to the analysis of specific and non-specific BOLD-signal changes related to end-tidal CO₂ (P_ETCO₂) in breath-hold fMRI studies. Multiple-order nonlinear predictors for P_ETCO₂ model a region-dependent nonlinear input-output relationship hypothesized in literature and possibly playing a crucial role in disentangling noise. We explore Retrospective Image-based Correction (RETROICOR) effects on the estimated BOLD response, applying our analysis both with and without RETROICOR and analyzing the linear and non-linear correlation between P_ETCO₂ and RETROICOR regressors.

RESULTS

The RETROICOR model of noise related to respiratory activity correlated with P_ETCO₂ both linearly and non-linearly. The correction affected the shape of the estimated BOLD response to hypercapnia but allowed to discard spurious activity in ventricles and white matter. Activation clusters were best detected using non-linear components in the BOLD response model.

COMPARISON WITH EXISTING METHOD

We evaluated the side-effects of standard physiological noise correction procedure, tailoring our analysis on challenging understudied brainstem and subcortical regions. Our novel approach allowed to characterize delays and non-linearities in BOLD response.

CONCLUSIONS

RETROICOR successfully avoided false positives, still broadly affecting the estimated non-linear BOLD responses. Non-linearities in the model better explained CO₂-related BOLD signal fluctuations. The necessity to modify the standard procedure for physiological-noise correction in breath-hold studies was addressed, stating its crucial importance.

Evaluating the sensitivity of functional connectivity measures to motion artifact in resting-state fMRI data

Functional connectivity (FC) networks are typically inferred from resting-state fMRI data using the Pearson correlation between BOLD time series from pairs of brain regions. However, alternative methods of estimating functional connectivity have not been systematically tested for their sensitivity or robustness to head motion artifact. Here, we evaluate the sensitivity of eight different functional connectivity measures to motion artifact using resting-state data from the Human Connectome Project. We report that FC estimated using full correlation has a relatively high residual distance-dependent relationship with motion compared to partial correlation, coherence, and information theory-based measures, even after implementing rigorous methods for motion artifact mitigation. This disadvantage of full correlation, however, may be offset by higher test-retest reliability, fingerprinting accuracy, and system identifiability. FC estimated by partial correlation offers the best of both worlds, with low sensitivity to motion artifact and intermediate system identifiability, with the caveat of low test-retest reliability and fingerprinting accuracy. We highlight spatial differences in the sub-networks affected by motion with different FC metrics. Further, we report that intra-network edges in the default mode and retrosplenial temporal sub-networks are highly correlated with motion in all FC methods. Our findings indicate that the method of estimating functional connectivity is an important consideration in resting-state fMRI studies and must be chosen carefully based on the parameters of the study.

Reconstruction of respiratory variation signals from fMRI data

Functional MRI signals can be heavily influenced by systemic physiological processes in addition to local neural activity. For example, widespread hemodynamic fluctuations across the brain have been found to correlate with natural, low-frequency variations in the depth and rate of breathing over time. Acquiring peripheral measures of respiration during fMRI scanning not only allows for modeling such effects in fMRI analysis, but also provides valuable information for interrogating brain-body physiology. However, physiological recordings are frequently unavailable or have insufficient quality. Here, we propose a computational technique for reconstructing continuous low-frequency respiration volume (RV) fluctuations from fMRI data alone. We evaluate the performance of this approach across different fMRI preprocessing strategies. Further, we demonstrate that the predicted RV signals can account for similar patterns of temporal variation in resting-state fMRI data compared to measured RV fluctuations. These findings indicate that fluctuations in respiration volume can be extracted from fMRI alone, in the common scenario of missing or corrupted respiration recordings. The results have implications for enriching a large volume of existing fMRI datasets through retrospective addition of respiratory variations information.

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.