Infraslow LFP correlates to resting-state fMRI BOLD signals

Networks extracted from nonlinear fMRI connectivity exhibit unique spatial variation and enhanced sensitivity to differences between individuals with schizophrenia and controls

Schizophrenia is a chronic brain disorder associated with widespread alterations in functional brain connectivity. Although data-driven approaches such as independent component analysis are often used to study how schizophrenia impacts linearly connected networks, alterations within the underlying nonlinear functional connectivity structure remain largely unknown. Here we report the analysis of networks from explicitly nonlinear functional magnetic resonance imaging connectivity in a case-control dataset. We found systematic spatial variation, with higher nonlinear weight within core regions, suggesting that linear analyses underestimate functional connectivity within network centers. We also found that a unique nonlinear network incorporating default-mode, cingulo-opercular and central executive regions exhibits hypoconnectivity in schizophrenia, indicating that typically hidden connectivity patterns may reflect inefficient network integration in psychosis. Moreover, nonlinear networks including those previously implicated in auditory, linguistic and self-referential cognition exhibit heightened statistical sensitivity to schizophrenia diagnosis, collectively underscoring the potential of our methodology to resolve complex brain phenomena and transform clinical connectivity analysis.

Neuronal basis of high frequency fMRI fluctuation: direct evidence from simultaneous recording

Resting-state functional magnetic resonance imaging (RS-fMRI) has been extensively utilized for noninvasive investigation of human brain activity. While studies employing simultaneous recordings of fMRI and electrophysiology have established a connection between the low-frequency fluctuation (< 0.1 Hz) observed in RS-fMRI and the local field potential (LFP), it remains unclear whether the RS-fMRI signal exhibits frequency-dependent modulation, which is a well-documented phenomenon in LFP. The present study concurrently recorded resting-state functional magnetic resonance imaging (RS-fMRI) and local field potentials (LFP) in the striatum of 8 rats before and after a pharmacological manipulation. We observed a highly similar frequency-dependent pattern of amplitude changes in both RS-fMRI and LFP following the manipulation, specifically an increase in high-frequency band amplitudes accompanied by a decrease in low-frequency band amplitudes. These findings provide direct evidence that the enhanced high-frequency fluctuations and reduced low-frequency fluctuations observed in RS-fMRI may reflect heightened neuronal activity.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by 'electrophysiology-invisible' signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Variation in the distribution of large-scale spatiotemporal patterns of activity across brain states

A few large-scale spatiotemporal patterns of brain activity (quasiperiodic patterns or QPPs) account for most of the spatial structure observed in resting state functional magnetic resonance imaging (rs-fMRI). The QPPs capture well-known features such as the evolution of the global signal and the alternating dominance of the default mode and task positive networks. These widespread patterns of activity have plausible ties to neuromodulatory input that mediates changes in nonlocalized processes, including arousal and attention. To determine whether QPPs exhibit variations across brain conditions, the relative magnitude and distribution of the three strongest QPPs were examined in two scenarios. First, in data from the Human Connectome Project, the relative incidence and magnitude of the QPPs was examined over the course of the scan, under the hypothesis that increasing drowsiness would shift the expression of the QPPs over time. Second, using rs-fMRI in rats obtained with a novel approach that minimizes noise, the relative incidence and magnitude of the QPPs was examined under three different anesthetic conditions expected to create distinct types of brain activity. The results indicate that both the distribution of QPPs and their magnitude changes with brain state, evidence of the sensitivity of these large-scale patterns to widespread changes linked to alterations in brain conditions.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by "electrophysiology-invisible" signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by "electrophysiology-invisible" signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Brain states as wave-like motifs

There is ample evidence of wave-like activity in the brain at multiple scales and levels. This emerging literature supports the broader adoption of a wave perspective of brain activity. Specifically, a brain state can be described as a set of recurring, sequential patterns of propagating brain activity, namely a wave. We examine a collective body of experimental work investigating wave-like properties. Based on these works, we consider brain states as waves using a scale-agnostic framework across time and space. Emphasis is placed on the sequentiality and periodicity associated with brain activity. We conclude by discussing the implications, prospects, and experimental opportunities of this framework.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and the underlying patterns of neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the amygdala of two male macaques. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5 Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-modal approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

The Effects of Locus Coeruleus Optogenetic Stimulation on Global Spatiotemporal Patterns in Rats

Whole-brain intrinsic activity as detected by resting-state fMRI can be summarized by three primary spatiotemporal patterns. These patterns have been shown to change with different brain states, especially arousal. The noradrenergic locus coeruleus (LC) is a key node in arousal circuits and has extensive projections throughout the brain, giving it neuromodulatory influence over the coordinated activity of structurally separated regions. In this study, we used optogenetic-fMRI in rats to investigate the impact of LC stimulation on the global signal and three primary spatiotemporal patterns. We report small, spatially specific changes in global signal distribution as a result of tonic LC stimulation, as well as regional changes in spatiotemporal patterns of activity at 5 Hz tonic and 15 Hz phasic stimulation. We also found that LC stimulation had little to no effect on the spatiotemporal patterns detected by complex principal component analysis. These results show that the effects of LC activity on the BOLD signal in rats may be small and regionally concentrated, as opposed to widespread and globally acting.

Variation in the Distribution of Large-scale Spatiotemporal Patterns of Activity Across Brain States

A few large-scale spatiotemporal patterns of brain activity (quasiperiodic patterns or QPPs) account for most of the spatial structure observed in resting state functional magnetic resonance imaging (rs-fMRI). The QPPs capture well-known features such as the evolution of the global signal and the alternating dominance of the default mode and task positive networks. These widespread patterns of activity have plausible ties to neuromodulatory input that mediates changes in nonlocalized processes, including arousal and attention. To determine whether QPPs exhibit variations across brain conditions, the relative magnitude and distribution of the three strongest QPPs were examined in two scenarios. First, in data from the Human Connectome Project, the relative incidence and magnitude of the QPPs was examined over the course of the scan, under the hypothesis that increasing drowsiness would shift the expression of the QPPs over time. Second, using rs-fMRI in rats obtained with a novel approach that minimizes noise, the relative incidence and magnitude of the QPPs was examined under three different anesthetic conditions expected to create distinct types of brain activity. The results indicate that both the distribution of QPPs and their magnitude changes with brain state, evidence of the sensitivity of these large-scale patterns to widespread changes linked to alterations in brain conditions.

Chronic full-band recordings with graphene microtransistors as neural interfaces for discrimination of brain states

Brain states such as sleep, anesthesia, wakefulness, or coma are characterized by specific patterns of cortical activity dynamics, from local circuits to full-brain emergent properties. We previously demonstrated that full-spectrum signals, including the infraslow component (DC, direct current-coupled), can be recorded acutely in multiple sites using flexible arrays of graphene solution-gated field-effect transistors (gSGFETs). Here, we performed chronic implantation of 16-channel gSGFET arrays over the rat cerebral cortex and recorded full-band neuronal activity with two objectives: (1) to test the long-term stability of implanted devices; and (2) to investigate full-band activity during the transition across different levels of anesthesia. First, we demonstrate it is possible to record full-band signals with stability, fidelity, and spatiotemporal resolution for up to 5.5 months using chronic epicortical gSGFET implants. Second, brain states generated by progressive variation of levels of anesthesia could be identified as traditionally using the high-pass filtered (AC, alternating current-coupled) spectrogram: from synchronous slow oscillations in deep anesthesia through to asynchronous activity in the awake state. However, the DC signal introduced a highly significant improvement for brain-state discrimination: the DC band provided an almost linear information prediction of the depth of anesthesia, with about 85% precision, using a trained algorithm. This prediction rose to about 95% precision when the full-band (AC + DC) spectrogram was taken into account. We conclude that recording infraslow activity using gSGFET interfaces is superior for the identification of brain states, and further supports the preclinical and clinical use of graphene neural interfaces for long-term recordings of cortical activity.

Longitudinal monitoring of the mouse brain reveals heterogenous network trajectories during aging

The human aging brain is characterized by changes in network efficiency that are currently best captured through longitudinal resting-state functional MRI (rs-fMRI). These studies however are challenging due to the long human lifespan. Here we show that the mouse animal model with a much shorter lifespan allows us to follow the functional network organization over most of the animal's adult lifetime. We used a longitudinal study of the functional connectivity of different brain regions with rs-fMRI under anesthesia. Our analysis uncovers network modules similar to those reported in younger mice and in humans (i.e., prefrontal/default mode network (DMN), somatomotor and somatosensory networks). Statistical analysis reveals different patterns of network reorganization during aging. Female mice showed a pattern akin to human aging, with de-differentiation of the connectome, mainly due to increases in connectivity of the prefrontal/DMN cortical networks to other modules. Our male cohorts revealed heterogenous aging patterns with only one group confirming the de- differentiation, while the majority showed an increase in connectivity of the somatomotor cortex to the Nucleus accumbens. In summary, in line with human work, our analysis in mice supports the concept of de-differentiation in the aging mammalian brain and reveals additional trajectories in aging mice networks.

Time-varying functional connectivity predicts fluctuations in sustained attention in a serial tapping task

The mechanisms for how large-scale brain networks contribute to sustained attention are unknown. Attention fluctuates from moment to moment, and this continuous change is consistent with dynamic changes in functional connectivity between brain networks involved in the internal and external allocation of attention. In this study, we investigated how brain network activity varied across different levels of attentional focus (i.e., "zones"). Participants performed a finger-tapping task, and guided by previous research, in-the-zone performance or state was identified by low reaction time variability and out-of-the-zone as the inverse. In-the-zone sessions tended to occur earlier in the session than out-of-the-zone blocks. This is unsurprising given the way attention fluctuates over time. Employing a novel method of time-varying functional connectivity, called the quasi-periodic pattern analysis (i.e., reliable, network-level low-frequency fluctuations), we found that the activity between the default mode network (DMN) and task positive network (TPN) is significantly more anti-correlated during in-the-zone states versus out-of-the-zone states. Furthermore, it is the frontoparietal control network (FPCN) switch that differentiates the two zone states. Activity in the dorsal attention network (DAN) and DMN were desynchronized across both zone states. During out-of-the-zone periods, FPCN synchronized with DMN, while during in-the-zone periods, FPCN switched to synchronized with DAN. In contrast, the ventral attention network (VAN) synchronized more closely with DMN during in-the-zone periods compared with out-of-the-zone periods. These findings demonstrate that time-varying functional connectivity of low frequency fluctuations across different brain networks varies with fluctuations in sustained attention or other processes that change over time.

Spatial and Spectral Components of the BOLD Global Signal in Rat Resting-State Functional MRI

PURPOSE

In resting-state fMRI (rs-fMRI), the global signal average captures widespread fluctuations related to unwanted sources of variance such as motion and respiration, as well as widespread neural activity; however, relative contributions of neural and non-neural sources to the global signal remain poorly understood. This study sought to tackle this problem through the comparison of the BOLD global signal to an adjacent non-brain tissue signal, where neural activity was absent, from the same rs-fMRI scan obtained from anesthetized rats. In this dataset, motion was minimal and ventilation was phase-locked to image acquisition to minimize respiratory fluctuations. Data were acquired using three different anesthetics: isoflurane, dexmedetomidine, and a combination of dexmedetomidine and light isoflurane.

METHODS

A power spectral density estimate, a voxel-wise spatial correlation via Pearson's correlation, and a co-activation pattern analysis were performed using the global signal and the non-brain tissue signal. Functional connectivity was calculated using Pearson's linear correlation on default mode network (DMN) regions.

RESULTS

We report differences in the spectral composition of the two signals and show spatial selectivity within DMN structures that show an increased correlation to the global signal and decreased intra-network connectivity after global signal regression. All of the observed differences between the global signal and the non-brain tissue signal were maintained across anesthetics.

CONCLUSION

These results show that the global signal is distinct from the noise contained in the tissue signal, as support for a neural contribution. This study provides a unique perspective to the contents of the global signal and their origins.

Networks extracted from nonlinear fMRI connectivity exhibit unique spatial variation and enhanced sensitivity to differences between individuals with schizophrenia and controls

Functional magnetic resonance imaging (fMRI) studies often estimate brain intrinsic connectivity networks (ICNs) from temporal relationships between hemodynamic signals using approaches such as independent component analysis (ICA). While ICNs are thought to represent functional sources that play important roles in various psychological phenomena, current approaches have been tailored to identify ICNs that mainly reflect linear statistical relationships. However, the elements comprising neural systems often exhibit remarkably complex nonlinear interactions that may be involved in cognitive operations and altered in psychiatric conditions such as schizophrenia. Consequently, there is a need to develop methods capable of effectively capturing ICNs from measures that are sensitive to nonlinear relationships. Here, we advance a novel approach to estimate ICNs from explicitly nonlinear whole-brain functional connectivity (ENL-wFC) by transforming resting-state fMRI (rsfMRI) data into the connectivity domain, allowing us to capture unique information from distance correlation patterns that would be missed by linear whole-brain functional connectivity (LIN-wFC) analysis. Our findings provide evidence that ICNs commonly extracted from linear (LIN) relationships are also reflected in explicitly nonlinear (ENL) connectivity patterns. ENL ICN estimates exhibit higher reliability and stability, highlighting our approach's ability to effectively quantify ICNs from rsfMRI data. Additionally, we observed a consistent spatial gradient pattern between LIN and ENL ICNs with higher ENL weight in core ICN regions, suggesting that ICN function may be subserved by nonlinear processes concentrated within network centers. We also found that a uniquely identified ENL ICN distinguished individuals with schizophrenia from healthy controls while a uniquely identified LIN ICN did not, emphasizing the valuable complementary information that can be gained by incorporating measures that are sensitive to nonlinearity in future analyses. Moreover, the ENL estimates of ICNs associated with auditory, linguistic, sensorimotor, and self-referential processes exhibit heightened sensitivity towards differentiating between individuals with schizophrenia and controls compared to LIN counterparts, demonstrating the translational value of our approach and of the ENL estimates of ICNs that are frequently reported as disrupted in schizophrenia. In summary, our findings underscore the tremendous potential of connectivity domain ICA and nonlinear information in resolving complex brain phenomena and revolutionizing the landscape of clinical FC analysis.

Neuromodulation of striatal D1 cells shapes BOLD fluctuations in anatomically connected thalamic and cortical regions

Understanding how the brain's macroscale dynamics are shaped by underlying microscale mechanisms is a key problem in neuroscience. In animal models, we can now investigate this relationship in unprecedented detail by directly manipulating cellular-level properties while measuring the whole-brain response using resting-state fMRI. Here, we focused on understanding how blood-oxygen-level-dependent (BOLD) dynamics, measured within a structurally well-defined striato-thalamo-cortical circuit in mice, are shaped by chemogenetically exciting or inhibiting D1 medium spiny neurons (MSNs) of the right dorsomedial caudate putamen (CPdm). We characterize changes in both the BOLD dynamics of individual cortical and subcortical brain areas, and patterns of inter-regional coupling (functional connectivity) between pairs of areas. Using a classification approach based on a large and diverse set of time-series properties, we found that CPdm neuromodulation alters BOLD dynamics within thalamic subregions that project back to dorsomedial striatum. In the cortex, changes in local dynamics were strongest in unimodal regions (which process information from a single sensory modality) and weakened along a hierarchical gradient towards transmodal regions. In contrast, a decrease in functional connectivity was observed only for cortico-striatal connections after D1 excitation. Our results show that targeted cellular-level manipulations affect local BOLD dynamics at the macroscale, such as by making BOLD dynamics more predictable over time by increasing its self-correlation structure. This contributes to ongoing attempts to understand the influence of structure-function relationships in shaping inter-regional communication at subcortical and cortical levels.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

The interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14~62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

Interacting spiral wave patterns underlie complex brain dynamics and are related to cognitive processing

The large-scale activity of the human brain exhibits rich and complex patterns, but the spatiotemporal dynamics of these patterns and their functional roles in cognition remain unclear. Here by characterizing moment-by-moment fluctuations of human cortical functional magnetic resonance imaging signals, we show that spiral-like, rotational wave patterns (brain spirals) are widespread during both resting and cognitive task states. These brain spirals propagate across the cortex while rotating around their phase singularity centres, giving rise to spatiotemporal activity dynamics with non-stationary features. The properties of these brain spirals, such as their rotational directions and locations, are task relevant and can be used to classify different cognitive tasks. We also demonstrate that multiple, interacting brain spirals are involved in coordinating the correlated activations and de-activations of distributed functional regions; this mechanism enables flexible reconfiguration of task-driven activity flow between bottom-up and top-down directions during cognitive processing. Our findings suggest that brain spirals organize complex spatiotemporal dynamics of the human brain and have functional correlates to cognitive processing.

The interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14∼62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

MRI texture-based radiomics analysis for the identification of altered functional networks in alcoholic patients and animal models

Alcohol use disorder (AUD) is a complex condition representing a leading risk factor for death, disease and disability. Its high prevalence and severe health consequences make necessary a better understanding of the brain network alterations to improve diagnosis and treatment. The purpose of this study was to evaluate the potential of resting-state fMRI 3D texture features as a novel source of biomarkers to identify AUD brain network alterations following a radiomics approach. A longitudinal study was conducted in Marchigian Sardinian alcohol-preferring msP rats (N = 36) who underwent resting-state functional and structural MRI before and after 30 days of alcohol or water consumption. A cross-sectional human study was also conducted among 33 healthy controls and 35 AUD patients. The preprocessed functional data corresponding to control and alcohol conditions were used to perform a probabilistic independent component analysis, identifying seven independent components as resting-state networks. Forty-three radiomic features extracted from each network were compared using a Wilcoxon signed-rank test with Holm correction to identify the network most affected by alcohol consumption. Features extracted from this network were then used in the machine learning process, evaluating two feature selection methods and six predictive models within a nested cross-validation structure. The classification was evaluated by computing the area under the ROC curve. Images were quantized using different numbers of gray-levels to test their influence on the results. The influence of ageing, data preprocessing, and brain iron accumulation were also analyzed. The methodology was validated using structural scans. The striatal network in alcohol-exposed msP rats presented the most significant number of altered features. The radiomics approach supported this result achieving good classification performance in animals (AUC = 0.915 ± 0.100, with 12 features) and humans (AUC = 0.724 ± 0.117, with 9 features) using a random forest model. Using the structural scans, high accuracy was achieved with a multilayer perceptron in both species (animals: AUC > 0.95 with 2 features, humans: AUC > 0.82 with 18 features). The best results were obtained using a feature selection method based on the p-value. The proposed radiomics approach is able to identify AUD patients and alcohol-exposed rats with good accuracy, employing a subset of 3D features extracted from fMRI. Furthermore, it can help identify relevant networks in drug addiction.

Intrinsic macroscale oscillatory modes driving long range functional connectivity in female rat brains detected by ultrafast fMRI

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals correlate across distant brain areas, shaping functionally relevant intrinsic networks. However, the generative mechanism of fMRI signal correlations, and in particular the link with locally-detected ultra-slow oscillations, are not fully understood. To investigate this link, we record ultrafast ultrahigh field fMRI signals (9.4 Tesla, temporal resolution = 38 milliseconds) from female rats across three anesthesia conditions. Power at frequencies extending up to 0.3 Hz is detected consistently across rat brains and is modulated by anesthesia level. Principal component analysis reveals a repertoire of modes, in which transient oscillations organize with fixed phase relationships across distinct cortical and subcortical structures. Oscillatory modes are found to vary between conditions, resonating at faster frequencies under medetomidine sedation and reducing both in number, frequency, and duration with the addition of isoflurane. Peaking in power within clear anatomical boundaries, these oscillatory modes point to an emergent systemic property. This work provides additional insight into the origin of oscillations detected in fMRI and the organizing principles underpinning spontaneous long-range functional connectivity.

Exploring electroencephalographic infraslow neurofeedback treatment for chronic low back pain: a double-blinded safety and feasibility randomized placebo-controlled trial

Chronic low back pain (CLBP) is a disabling condition worldwide. In CLBP, neuroimaging studies demonstrate abnormal activities in cortical areas responsible for pain modulation, emotional, and sensory components of pain experience [i.e., pregenual and dorsal anterior cingulate cortex (pgACC, dACC), and somatosensory cortex (SSC), respectively]. This pilot study, conducted in a university setting, evaluated the feasibility, safety, and acceptability of a novel electroencephalography-based infraslow-neurofeedback (EEG ISF-NF) technique for retraining activities in pgACC, dACC and SSC and explored its effects on pain and disability. Participants with CLBP (n = 60), recruited between July'20 to March'21, received 12 sessions of either: ISF-NF targeting pgACC, dACC + SSC, a ratio of pgACC*2/dACC + SSC, or Placebo-NF. Descriptive statistics demonstrated that ISF-NF training is feasible [recruitment rate (7 participants/month), dropouts (25%; 20-27%), and adherence (80%; 73-88%)], safe (no adverse events reported), and was moderate to highly acceptable [Mean ± SD: 7.8 ± 2.0 (pgACC), 7.5 ± 2.7 (dACC + SCC), 8.2 ± 1.9 (Ratio), and 7.7 ± 1.5 (Placebo)]. ISF-NF targeting pgACC demonstrated the most favourable clinical outcomes, with a higher proportion of participants exhibiting a clinically meaningful reduction in pain severity [53%; MD (95% CI): - 1.9 (- 2.7, - 1.0)], interference [80%; MD (95% CI): - 2.3 (- 3.5, - 1.2)], and disability [73%; MD (95% CI): - 4.5 (- 6.1, - 2.9)] at 1-month follow-up. ISF-NF training is a feasible, safe, and an acceptable treatment approach for CLBP.

The impact of vasomotion on analysis of rodent fMRI data

Introduction

Small animal fMRI is an essential part of translational research in the cognitive neurosciences. Due to small dimensions and animal physiology preclinical fMRI is prone to artifacts that may lead to misinterpretation of the data. To reach unbiased translational conclusions, it is, therefore, crucial to identify potential sources of experimental noise and to develop correction methods for contributions that cannot be avoided such as physiological noise. Aim of this study was to assess origin and prevalence of hemodynamic oscillations (HDO) in preclinical fMRI in rat, as well as their impact on data analysis.

Methods

Following the development of algorithms for HDO detection and suppression, HDO prevalence in fMRI measurements was investigated for different anesthetic regimens, comprising isoflurane and medetomidine, and for both gradient echo and spin echo fMRI sequences. In addition to assessing the effect of vasodilation on HDO, it was studied if HDO have a direct neuronal correlate using local field potential (LFP) recordings. Finally, the impact of HDO on analysis of fMRI data was assessed, studying both the impact on calculation of activation maps as well as the impact on brain network analysis. Overall, 303 fMRI measurements and 32 LFP recordings were performed in 71 rats.

Results

In total, 62% of the fMRI measurements showed HDO with a frequency of (0.20 ± 0.02) Hz. This frequent occurrence indicated that HDO cannot be generally neglected in fMRI experiments. Using the developed algorithms, HDO were detected with a specificity of 95%, and removed efficiently from the signal time courses. HDO occurred brain-wide under vasoconstrictive conditions in both small and large blood vessels. Vasodilation immediately interrupted HDO, which, however, returned within 1 h under vasoconstrictive conditions. No direct neuronal correlate of HDO was observed in LFP recordings. HDO significantly impacted analysis of fMRI data, leading to altered cluster sizes and F-values for activated voxels, as well as altered brain networks, when comparing data with and without HDO.

Discussion

We therefore conclude that HDO are caused by vasomotion under certain anesthetic conditions and should be corrected during fMRI data analysis to avoid bias.

Triple network activation causes tinnitus in patients with sudden sensorineural hearing loss: A model-based volume-entropy analysis

Tinnitus can be defined as the conscious perception of phantom sounds in the absence of corresponding external auditory signals. Tinnitus can develop in the setting of sudden sensorineural hearing loss (SSNHL), but the underlying mechanism is largely unknown. Using electroencephalography, we investigated differences in afferent node capacity between 15 SSNHL patients without tinnitus (NT) and 30 SSNHL patients with tinnitus (T). Where the T group showed increased afferent node capacity in regions constituting a "triple brain network" [default mode network (DMN), central executive network (CEN), and salience network (SN)], the NT group showed increased information flow in regions implicated in temporal auditory processing and noise-canceling pathways. Our results demonstrate that when all components of the triple network are activated due to sudden-onset auditory deprivation, tinnitus ensues. By contrast, auditory processing-associated and tinnitus-suppressing networks are highly activated in the NT group, to overcome the activation of the triple network and effectively suppress the generation of tinnitus.

Infraslow closed-loop brain training for anxiety and depression (ISAD): a protocol for a randomized, double-blind, sham-controlled pilot trial in adult females with internalizing disorders

BACKGROUND

The core intrinsic connectivity networks (core-ICNs), encompassing the default-mode network (DMN), salience network (SN) and central executive network (CEN), have been shown to be dysfunctional in individuals with internalizing disorders (IDs, e.g. major depressive disorder, MDD; generalized anxiety disorder, GAD; social anxiety disorder, SOC). As such, source-localized, closed-loop brain training of electrophysiological signals, also known as standardized low-resolution electromagnetic tomography (sLORETA) neurofeedback (NFB), targeting key cortical nodes within these networks has the potential to reduce symptoms associated with IDs and restore normal core ICN function. We intend to conduct a randomized, double-blind (participant and assessor), sham-controlled, parallel-group (3-arm) trial of sLORETA infraslow (<0.1 Hz) fluctuation neurofeedback (sLORETA ISF-NFB) 3 times per week over 4 weeks in participants (n=60) with IDs. Our primary objectives will be to examine patient-reported outcomes (PROs) and neurophysiological measures to (1) compare the potential effects of sham ISF-NFB to either genuine 1-region ISF-NFB or genuine 2-region ISF-NFB, and (2) assess for potential associations between changes in PRO scores and modifications of electroencephalographic (EEG) activity/connectivity within/between the trained regions of interest (ROIs). As part of an exploratory analysis, we will investigate the effects of additional training sessions and the potential for the potentiation of the effects over time.

METHODS

We will randomly assign participants who meet the criteria for MDD, GAD, and/or SOC per the MINI (Mini International Neuropsychiatric Interview for DSM-5) to one of three groups: (1) 12 sessions of posterior cingulate cortex (PCC) ISF-NFB up-training (n=15), (2) 12 sessions of concurrent PCC ISF up-training and dorsal anterior cingulate cortex (dACC) ISF-NFB down-training (n=15), or (3) 6 sessions of yoked-sham training followed by 6 sessions genuine ISF-NFB (n=30). Transdiagnostic PROs (Hospital Anxiety and Depression Scale, HADS; Inventory of Depression and Anxiety Symptoms - Second Version, IDAS-II; Multidimensional Emotional Disorder Inventory, MEDI; Intolerance of Uncertainty Scale - Short Form, IUS-12; Repetitive Thinking Questionnaire, RTQ-10) as well as resting-state neurophysiological measures (full-band EEG and ECG) will be collected from all subjects during two baseline sessions (approximately 1 week apart) then at post 6 sessions, post 12 sessions, and follow-up (1 month later). We will employ Bayesian methods in R and advanced source-localisation software (i.e. exact low-resolution brain electromagnetic tomography; eLORETA) in our analysis.

DISCUSSION

This protocol will outline the rationale and research methodology for a clinical pilot trial of sLORETA ISF-NFB targeting key nodes within the core-ICNs in a female ID population with the primary aims being to assess its potential efficacy via transdiagnostic PROs and relevant neurophysiological measures.

TRIAL REGISTRATION

Our study was prospectively registered with the Australia New Zealand Clinical Trials Registry (ANZCTR; Trial ID: ACTRN12619001428156). Registered on October 15, 2019.

Individual slow wave events give rise to macroscopic fMRI signatures and drive the strength of the BOLD signal in human resting-state EEG-fMRI recordings

The slow wave state is a general state of quiescence interrupted by sudden bursts of activity or so-called slow wave events (SWEs). Recently, the relationship between SWEs and blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals was assessed in rodent models which revealed cortex-wide BOLD activation. However, it remains unclear which macroscopic signature corresponds to these specific neurophysiological events in the human brain. Therefore, we analyzed simultaneous electroencephalographic (EEG)-fMRI data during human non-REM sleep. SWEs individually detected in the EEG data were used as predictors in event-related fMRI analyses to examine the relationship between SWEs and fMRI signals. For all 10 subjects we identified significant changes in BOLD activity associated with SWEs covering substantial parts of the gray matter. As demonstrated in rodents, we observed a direct relation of a neurophysiological event to specific BOLD activation patterns. We found a correlation between the number of SWEs and the spatial extent of these BOLD activation patterns and discovered that the amplitude of the BOLD response strongly depends on the SWE amplitude. As altered SWE propagation has recently been found in neuropsychiatric diseases, it is critical to reveal the brain's physiological slow wave state networks to potentially establish early imaging biomarkers for various diseases long before disease onset.

Receptor-informed network control theory links LSD and psilocybin to a flattening of the brain's control energy landscape

Psychedelics including lysergic acid diethylamide (LSD) and psilocybin temporarily alter subjective experience through their neurochemical effects. Serotonin 2a (5-HT2a) receptor agonism by these compounds is associated with more diverse (entropic) brain activity. We postulate that this increase in entropy may arise in part from a flattening of the brain's control energy landscape, which can be observed using network control theory to quantify the energy required to transition between recurrent brain states. Using brain states derived from existing functional magnetic resonance imaging (fMRI) datasets, we show that LSD and psilocybin reduce control energy required for brain state transitions compared to placebo. Furthermore, across individuals, reduction in control energy correlates with more frequent state transitions and increased entropy of brain state dynamics. Through network control analysis that incorporates the spatial distribution of 5-HT2a receptors (obtained from publicly available positron emission tomography (PET) data under non-drug conditions), we demonstrate an association between the 5-HT2a receptor and reduced control energy. Our findings provide evidence that 5-HT2a receptor agonist compounds allow for more facile state transitions and more temporally diverse brain activity. More broadly, we demonstrate that receptor-informed network control theory can model the impact of neuropharmacological manipulation on brain activity dynamics.

A systematic data-driven approach to analyze sensor-level EEG connectivity: Identifying robust phase-synchronized network components using surface Laplacian with spectral-spatial PCA

Although conventional averaging across predefined frequency bands reduces the complexity of EEG functional connectivity (FC), it obscures the identification of resting-state brain networks (RSN) and impedes accurate estimation of FC reliability. Extending prior work, we combined scalp current source density (CSD; spherical spline surface Laplacian) and spectral-spatial PCA to identify FC components. Phase-based FC was estimated via debiased-weighted phase-locking index from CSD-transformed resting EEGs (71 sensors, 8 min, eyes open/closed, 35 healthy adults, 1-week retest). Spectral PCA extracted six robust alpha and theta components (86.6% variance). Subsequent spatial PCA for each spectral component revealed seven robust regionally focused (posterior, central, and frontal) and long-range (posterior-anterior) alpha components (peaks at 8, 10, and 13 Hz) and a midfrontal theta (6 Hz) component, accounting for 37.0% of FC variance. These spatial FC components were consistent with well-known networks (e.g., default mode, visual, and sensorimotor), and four were sensitive to eyes open/closed conditions. Most FC components had good-to-excellent internal consistency (odd/even epochs, eyes open/closed) and test-retest reliability (ICCs ≥ .8). Moreover, the FC component structure was generally present in subsamples (session × odd/even epoch, or smaller subgroups [n = 7-10]), as indicated by high similarity of component loadings across PCA solutions. Apart from systematically reducing FC dimensionality, our approach avoids arbitrary thresholds and allows quantification of meaningful and reliable network components that may prove to be of high relevance for basic and clinical research applications.

Dose-response relationship between the variables of unilateral optogenetic stimulation and transcallosal evoked responses in rat motor cortex

Efficient interhemispheric integration of neural activity between left and right primary motor cortex (M1) is critical for inter-limb motor control. We employed optogenetic stimulation to establish a framework for probing transcallosal M1–M1 interactions in rats. We performed optogenetic stimulation of excitatory neurons in right M1 of male Sprague-Dawley rats. We recorded the transcallosal evoked potential in contralateral left M1 via chronically implanted electrodes. Recordings were performed under anesthesia combination of dexmedetomidine and a low concentration of isoflurane. We systematically varied the stimulation intensity and duration to characterize the relationship between stimulation parameters in right M1 and the characteristics of the evoked intracortical potentials in left M1. Optogenetic stimulation of right M1 consistently evoked a transcallosal response in left M1 with a consistent negative peak (N1) that sometimes was preceded by a smaller positive peak (P1). Higher stimulation intensity or longer stimulation duration gradually increased N1 amplitude and reduced N1 variability across trials. A combination of stimulation intensities of 5–10 mW with stimulus durations of 1–10 ms were generally sufficient to elicit a robust transcallosal response in most animal, with our optic fiber setup. Optogenetically stimulated excitatory neurons in M1 can reliably evoke a transcallosal response in anesthetized rats. Characterizing the relationship between “stimulation dose” and “response magnitude” (i.e., the gain function) of transcallosal M1-to-M1 excitatory connections can be used to optimize the variables of optogenetic stimulation and ensure stimulation efficacy.

Digital reconstruction of infraslow activity in human intracranial ictal recordings using a deconvolution-based inverse filter

Infraslow activity (ISA) is a biomarker that has recently become of interest in the characterization of seizure recordings. Recent data from a small number of studies have suggested that the epileptogenic zone may be identified by the presence of ISA. Investigation of low frequency activity in clinical seizure recordings, however, has been hampered by technical limitations. EEG systems necessarily include a high-pass filter early in the measurement chain to remove large artifactual drifts that can saturate recording elements such as the amplifier. This filter, unfortunately, attenuates legitimately seizure-related low frequencies, making ISA difficult to study in clinical EEG recordings. In this study, we present a deconvolution-based digital inverse filter that allows recovery of attenuated low frequency activity in intracranial recordings of temporal lobe epilepsy patients. First, we show that the unit impulse response (UIR) of an EEG system can be characterized by differentiation of the system's step response. As proof of method, we present several examples that show that the low frequency component of a high-pass filtered signal can be restored by deconvolution with the UIR. We then demonstrate that this method can be applied to biologically relevant signals including clinical EEG recordings obtained from seizure patients. Finally, we discuss how this method can be applied to study ISA to identify and assess the seizure onset zone.

Selective blockade of rat brain T-type calcium channels provides insights on neurophysiological basis of arousal dependent resting state functional magnetic resonance imaging signals

A number of studies point to slow (0.1–2 Hz) brain rhythms as the basis for the resting-state functional magnetic resonance imaging (rsfMRI) signal. Slow waves exist in the absence of stimulation, propagate across the cortex, and are strongly modulated by vigilance similar to large portions of the rsfMRI signal. However, it is not clear if slow rhythms serve as the basis of all neural activity reflected in rsfMRI signals, or just the vigilance-dependent components. The rsfMRI data exhibit quasi-periodic patterns (QPPs) that appear to increase in strength with decreasing vigilance and propagate across the brain similar to slow rhythms. These QPPs can complicate the estimation of functional connectivity (FC) via rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel-time courses of functionally connected brain regions. In this study, we examined the relationship between cortical slow rhythms and the rsfMRI signal, using a well-established pharmacological model of slow wave suppression. Suppression of cortical slow rhythms led to significant reduction in the amplitude of QPPs but increased rsfMRI measures of intrinsic FC in rats. The results suggest that cortical slow rhythms serve as the basis of only the vigilance-dependent components (e.g., QPPs) of rsfMRI signals. Further attenuation of these non-specific signals enhances delineation of brain functional networks.

The enigma of infra-slow fluctuations in the human EEG

Spontaneous Infra-Slow Fluctuations (ISFs) of the human EEG (EEG-ISFs) were discovered 60 years ago when appropriate amplifiers for their recordings were designed. To avoid skin-related artifacts the recording of EEG-ISFs required puncturing the skin under the electrode. In the beginning of the 21st century the interest in EEG-ISFs was renewed with the appearance of commercially available DC-coupled amplified and by observation of ISFs of the blood oxygen level–dependent (BOLD) functional magnetic resonance imaging signal at a similar frequency. The independent components of irregular EEG-ISFs were shown to correlate with BOLD signals which in turn were driven by changes in arousal level measured by galvanic skin response (GSR), pupil size and HRV. There is no consensus regarding the temporal organization of EEG-ISFs: some studies emphasize the absence of peaks on EEG-ISFs spectra, some studies report prominent oscillations with frequency around 0.1 or 0.02 Hz, while some emphasize multiple discrete infraslow oscillations. No studies used parameters of EEG-ISFs as neuromarkers to discriminate psychiatric patients from healthy controls. Finally, a set of working hypotheses is suggested that must be tested in future research to solve the enigma of EEG-ISFs.

Source localized infraslow neurofeedback training in people with chronic painful knee osteoarthritis: A randomized, double-blind, sham-controlled feasibility clinical trial

Persistent pain is a key symptom in people living with knee osteoarthritis (KOA). Infra-slow Neurofeedback (ISF-NF) training is a recent development focusing on modulating cortical slow-wave activity to improve pain outcomes. A parallel, two-armed double-blinded, randomized sham-controlled, feasibility clinical trial aimed to determine the feasibility and safety of a novel electroencephalography-based infraslow fluctuation neurofeedback (EEG ISF-NF) training in people with KOA and determine the variability of clinical outcomes and EEG changes following NF training. Eligible participants attended nine 30-min ISF-NF training sessions involving three cortical regions linked to pain. Feasibility measures were monitored during the trial period. Pain and functional outcomes were measured at baseline, post-intervention, and follow-up after 2 weeks. Resting-state EEG was recorded at baseline and immediate post-intervention. Participants were middle-aged (61.7 ± 7.6 years), New Zealand European (90.5%), and mostly females (62%) with an average knee pain duration of 4 ± 3.4 years. The study achieved a retention rate of 91%, with 20/22 participants completing all the sessions. Participants rated high levels of acceptance and “moderate to high levels of perceived effectiveness of the training.” No serious adverse events were reported during the trial. Mean difference (95% CI) for clinical pain and function measures are as follows for pain severity [active: 0.89 ± 1.7 (−0.27 to 2.0); sham: 0.98 ± 1.1 (0.22–1.7)], pain interference [active: 0.75 ± 2.3 (−0.82 to 2.3); Sham: 0.89 ± 2.1 (−0.60 to 2.4)], pain unpleasantness [active: 2.6 ± 3.7 (0.17–5.1); sham: 2.8 ± 3 (0.62–5.0)] and physical function [active: 6.2 ± 13 (−2.6 to 15); sham: 1.6 ± 12 (−6.8 to 10)]. EEG sources demonstrated frequency-specific neuronal activity, functional connectivity, and ISF ratio changes following NF training. The findings of the study indicated that the ISF-NF training is a feasible, safe, and acceptable intervention for pain management in people with KOA, with high levels of perceived effectiveness. The study also reports the variability in clinical, brain activity, and connectivity changes following training.

Spatiotemporal patterns of spontaneous brain activity: a mini-review

The brain exists in a state of constant activity in the absence of any external sensory input. The spatiotemporal patterns of this spontaneous brain activity have been studied using various recording and imaging techniques. This has enabled considerable progress to be made in elucidating the cellular and network mechanisms that are involved in the observed spatiotemporal dynamics. This mini-review outlines different spatiotemporal dynamic patterns that have been identified in four commonly used modalities: electrophysiological recordings, optical imaging, functional magnetic resonance imaging, and electroencephalography. Signal sources for each modality, possible sources of the observed dynamics, and future directions are also discussed.

Functional Connectivity Dynamics Altered of the Resting Brain in Subjective Cognitive Decline

Background

Subjective cognitive decline (SCD) appears in the preclinical stage of the Alzheimer's disease continuum. In this stage, dynamic features are more sensitive than static features to reflect early subtle changes in functional brain connectivity. Therefore, we studied local and extended dynamic connectivity of the resting brain of people with SCD to determine their intrinsic brain changes.

Methods

We enrolled cognitively normal older adults from the communities and divided them into SCD and normal control (NC) groups. We used mean dynamic amplitude of low-frequency fluctuation (mdALFF) to evaluate region of interest (ROI)-wise local dynamic connectivity of resting-state functional MRI. The dynamic functional connectivity (dFC) between ROIs was tested by whole-brain-based statistics.

Results

When comparing SCD (N = 40) with NC (N = 45), mdALFFmean decreased at right inferior parietal lobule (IPL) of the frontoparietal network (FPN). Still, it increased at the right middle temporal gyrus (MTG) of the ventral attention network (VAN) and right calcarine of the visual network (VIS). Also, the mdALFFvar (variance) increased at the left superior temporal gyrus of AUD, right MTG of VAN, right globus pallidum of the cingulo-opercular network (CON), and right lingual gyrus of VIS. Furthermore, mdALFFmean at right IPL of FPN are correlated negatively with subjective complaints and positively with objective cognitive performance. In the dFC seeded from the ROIs with local mdALFF group differences, SCD showed a generally lower dFCmean and higher dFCvar (variance) to other regions of the brain. These weakened and unstable functional connectivity appeared among FPN, CON, the default mode network, and the salience network, the large-scale networks of the triple network model for organizing neural resource allocations.

Conclusion

The local dynamic connectivity of SCD decreased in brain regions of cognitive executive control. Meanwhile, compensatory visual efforts and bottom-up attention rose. Mixed decrease and compensatory increase of dynamics of intrinsic brain activity suggest the transitional nature of SCD. The FPN local dynamics balance subjective and objective cognition and maintain cognitive preservation in preclinical dementia. Aberrant triple network model features the dFC alternations of SCD. Finally, the right lateralization phenomenon emerged early in the dementia continuum and affected local dynamic connectivity.

High-definition transcranial infraslow pink noise stimulation for chronic low back pain: protocol for a pilot, safety and feasibility randomised placebo-controlled trial

INTRODUCTION

Chronic low back pain (CLBP) is a common disabling health condition. Current treatments demonstrate modest effects, warranting newer therapies. Brain imaging demonstrates altered electrical activities in cortical areas responsible for pain modulation, emotional and sensory components of pain experience. Treatments targeting to change electrical activities of these key brain regions may produce clinical benefits. This pilot study aims to (1) evaluate feasibility, safety and acceptability of a novel neuromodulation technique, high-definition transcranial infraslow pink noise stimulation (HD-tIPNS), in people with CLBP, (2) explore the trend of effect of HD-tIPNS on pain and function, and (3) derive treatment estimates to support sample size calculation for a fully powered trial should trends of effectiveness be present.

METHODS AND ANALYSIS

A pilot, triple-blinded randomised two-arm placebo-controlled parallel trial. Participants (n=40) with CLBP will be randomised to either sham stimulation or HD-tIPNS (targeting somatosensory cortex and dorsal and pregenual anterior cingulate cortex). Primary outcomes include feasibility and safety measures, and clinical outcomes of pain (Brief Pain Inventory) and disability (Roland-Morris disability questionnaire). Secondary measures include clinical, psychological, quantitative sensory testing and electroencephalography collected at baseline, immediately postintervention, and at 1-week, 1-month and 3 months postintervention. All data will be analysed descriptively. A nested qualitative study will assess participants perceptions about acceptability of intervention and analysed thematically.

ETHICS AND DISSEMINATION

Ethical approval has been obtained from Health and Disability Ethics Committee (Ref:20/NTB/67). Findings will be reported to regulatory and funding bodies, presented at conferences, and published in a scientific journal.

TRIAL REGISTRATION NUMBER

ACTRN12620000505909p.

Brain Network Allostasis after Chronic Alcohol Drinking Is Characterized by Functional Dedifferentiation and Narrowing

Alcohol use disorder (AUD) causes complex alterations in the brain that are poorly understood. The heterogeneity of drinking patterns and the high incidence of comorbid factors compromise mechanistic investigations in AUD patients. Here we used male Marchigian Sardinian alcohol-preferring (msP) rats, a well established animal model of chronic alcohol drinking, and a combination of longitudinal resting-state fMRI and manganese-enhanced MRI to provide objective measurements of brain connectivity and activity, respectively. We found that 1 month of chronic alcohol drinking changed the correlation between resting-state networks. The change was not homogeneous, resulting in the reorganization of pairwise interactions and a shift in the equilibrium of functional connections. We identified two fundamentally different forms of network reorganization. First is functional dedifferentiation, which is defined as a regional increase in neuronal activity and overall correlation, with a concomitant decrease in preferential connectivity between specific networks. Through this mechanism, occipital cortical areas lost their specific interaction with sensory-insular cortex, striatal, and sensorimotor networks. Second is functional narrowing, which is defined as an increase in neuronal activity and preferential connectivity between specific brain networks. Functional narrowing strengthened the interaction between striatal and prefrontocortical networks, involving the anterior insular, cingulate, orbitofrontal, prelimbic, and infralimbic cortices. Importantly, these two types of alterations persisted after alcohol discontinuation, suggesting that dedifferentiation and functional narrowing rendered persistent network states. Our results support the idea that chronic alcohol drinking, albeit at moderate intoxicating levels, induces an allostatic change in the brain functional connectivity that propagates into early abstinence.SIGNIFICANCE STATEMENT Excessive consumption of alcohol is positioned among the top five risk factors for disease and disability. Despite this priority, the transformations that the nervous system undergoes from an alcohol-naive state to a pathologic alcohol drinking are not well understood. In our study, we use an animal model with proven translational validity to study this transformation longitudinally. The results show that shortly after chronic alcohol consumption there is an increase in redundant activity shared by brain structures, and the specific communication shrinks to a set of pathways. This functional dedifferentiation and narrowing are not reversed immediately after alcohol withdrawal but persist during early abstinence. We causally link chronic alcohol drinking with an early and abstinence-persistent retuning of the functional equilibrium of the brain.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

Brain Infraslow Activity Correlates With Arousal Levels

The functional role of the brain’s infraslow activity (ISA, 0.01–0.1 Hz) in human behavior has yet to be elucidated. To date, it has been shown that the brain’s ISA correlates with behavioral performance; task performance is more likely to increase when executed at a specific ISA phase. However, it is unclear how the ISA correlates behavioral performance. We hypothesized that the ISA phase correlation of behavioral performance is mediated by arousal. Our data analysis results showed that the electroencephalogram (EEG) ISA phase was correlated with the galvanic skin response (GSR) amplitude, a measure of the arousal level. Furthermore, subjects whose EEG ISA phase correlated with the GSR amplitude more strongly also showed greater EEG ISA modulation during meditation, which implies an intimate relationship between brain ISA and arousal. These results may help improve understanding of the functional role of the brain’s ISA.

Increased fMRI connectivity upon chemogenetic inhibition of the mouse prefrontal cortex

While shaped and constrained by axonal connections, fMRI-based functional connectivity reorganizes in response to varying interareal input or pathological perturbations. However, the causal contribution of regional brain activity to whole-brain fMRI network organization remains unclear. Here we combine neural manipulations, resting-state fMRI and in vivo electrophysiology to probe how inactivation of a cortical node causally affects brain-wide fMRI coupling in the mouse. We find that chronic inhibition of the medial prefrontal cortex (PFC) via overexpression of a potassium channel increases fMRI connectivity between the inhibited area and its direct thalamo-cortical targets. Acute chemogenetic inhibition of the PFC produces analogous patterns of fMRI overconnectivity. Using in vivo electrophysiology, we find that chemogenetic inhibition of the PFC enhances low frequency (0.1–4 Hz) oscillatory power via suppression of neural firing not phase-locked to slow rhythms, resulting in increased slow and δ band coherence between areas that exhibit fMRI overconnectivity. These results provide causal evidence that cortical inactivation can counterintuitively increase fMRI connectivity via enhanced, less-localized slow oscillatory processes.

Pathological perturbation affects whole brain network activity. Here the authors show in mice that cortical inactivation unexpectedly results in increased fMRI connectivity between the manipulated regions and its direct axonal targets.

Gaining insight into the neural basis of resting-state fMRI signal

The blood oxygenation level-dependent (BOLD)-based resting-state functional magnetic resonance imaging (rsfMRI) has been widely used as a non-invasive tool to map brain-wide connectivity architecture. However, the neural basis underpinning the resting-state BOLD signal remains elusive. In this study, we combined simultaneous calcium-based fiber photometry with rsfMRI in awake animals to examine the relationship of the BOLD signal and spiking activity at the resting state. We observed robust couplings between calcium and BOLD signals in the dorsal hippocampus as well as other distributed areas in the default mode network (DMN), suggesting that the calcium measurement can reliably predict the rsfMRI signal. In addition, using the calcium signal recorded as the ground truth, we assessed the impacts of different rsfMRI data preprocessing pipelines on functional connectivity mapping. Overall, our results provide important evidence suggesting that spiking activity measured by the calcium signal plays a key role in the neural mechanism of resting-state BOLD signal.

Relationships between correlated spikes, oxygen and LFP in the resting-state primate

Resting-state functional MRI (rsfMRI) provides a view of human brain organization based on correlation patterns of blood oxygen level dependent (BOLD) signals recorded across the whole brain. The neural basis of resting-state BOLD fluctuations and their correlation remains poorly understood. We simultaneously recorded oxygen level, spikes, and local field potential (LFP) at multiple sites in awake, resting monkeys. Following a spike, the average local oxygen and LFP voltage responses each resemble a task-driven BOLD response, with LFP preceding oxygen by 0.5 s. Between sites, features of the long-range correlation patterns of oxygen, LFP, and spikes are similar to features seen in rsfMRI. Most of the variance shared between sites lies in the infraslow frequency band (0.01-0.1 Hz) and in the infraslow envelope of higher-frequency bands (e.g. gamma LFP). While gamma LFP and infraslow LFP are both strong correlates of local oxygen, infraslow LFP explains significantly more of the variance shared between correlated oxygen signals than any other electrophysiological signal. Together these findings are consistent with a causal relationship between infraslow LFP and long-range oxygen correlations in the resting state.

MRI magnetic compatible electrical neural interface: From materials to application

Neural electrical interfaces are important tools for local neural stimulation and recording, which potentially have wide application in the diagnosis and treatment of neural diseases, as well as in the transmission of neural activity for brain-computer interface (BCI) systems. At the same time, magnetic resonance imaging (MRI) is one of the effective and non-invasive techniques for recording whole-brain signals, providing details of brain structures and also activation pattern maps. Simultaneous recording of extracellular neural signals and MRI combines two expressions of the same neural activity and is believed to be of great importance for the understanding of brain function. However, this combination makes requests on the magnetic and electronic performance of neural interface devices. MRI-compatibility refers here to a technological approach to simultaneous MRI and electrode recording or stimulation without artifacts in imaging. Trade-offs between materials magnetic susceptibility selection and electrical function should be considered. Herein, prominent trends in selecting materials of suitable magnetic properties are analyzed and material design, function and application of neural interfaces are outlined together with the remaining challenge to fabricate MRI-compatible neural interface.

State-related neural influences on fMRI connectivity estimation

The spatiotemporal structure of functional magnetic resonance imaging (fMRI) signals has provided a valuable window into the network underpinnings of human brain function and dysfunction. Although some cross-regional temporal correlation patterns (functional connectivity; FC) exhibit a high degree of stability across individuals and species, there is growing acknowledgment that measures of FC can exhibit marked changes over a range of temporal scales. Further, FC can covary with experimental task demands and ongoing neural processes linked to arousal, consciousness and perception, cognitive and affective state, and brain-body interactions. The increased recognition that such interrelated neural processes modulate FC measurements has raised both challenges and new opportunities in using FC to investigate brain function. Here, we review recent advances in the quantification of neural effects that shape fMRI FC and discuss the broad implications of these findings in the design and analysis of fMRI studies. We also discuss how a more complete understanding of the neural factors that shape FC measurements can resolve apparent inconsistencies in the literature and lead to more interpretable conclusions from fMRI studies.

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.

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.