Quasi-periodic patterns contribute to functional connectivity in the brain

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.

Creative tempo: Spatiotemporal dynamics of the default mode network in improvisational musicians

The intrinsic dynamics of human brain activity display a recurring pattern of anti-correlated activity between the default mode network (DMN), associated with internal processing and mentation, and task positive regions, associated with externally directed attention. In human functional magnetic resonance imaging (fMRI) data, this anti-correlated pattern is detectable on the infraslow timescale (<0.1 Hz) as a quasi-periodic pattern (QPP). While the DMN is implicated in creativity and musicality in traditional time-averaged functional connectivity studies, no one has yet explored how creative training may alter dynamic spatiotemporal patterns involving the DMN such as QPPs. In the present study, we compare the outputs of two QPP detection approaches, sliding window algorithm and complex principal components analysis (cPCA). We apply both methods to an existing dataset of musicians captured with resting state fMRI, grouped as either classical, improvisational, or minimally trained non-musicians. The original time-averaged functional connectivity (FC) analysis of this dataset used improvisation as a proxy for creative thinking and found that the DMN and visual networks (VIS) display higher connectivity in improvisational musicians. We expand upon this dataset's original study and find that QPP analysis detects convergent results at the group level with both methods. In improvisational musicians, dynamic functional correlation in the group-averaged QPP was found to be increased between the DMN-VIS and DMN-FPN for both the QPP algorithm and complex principal components analysis (cPCA) methods. Additionally, we found an unexpected increase in FC in the group-averaged QPP between the dorsal attention network and amygdala in improvisational musicians; this result was not reported in the original seed-based study of this dataset. The current study represents a novel application of two dynamic FC detection methods with results that replicate and expand upon previous seed-based FC findings. The results show the robustness of both the QPP phenomenon and its detection methods. This study also demonstrates the value of dynamic FC methods in reproducing seed-based findings and their promise in detecting group-wise or individual differences that may be missed by traditional seed-based resting state fMRI studies.

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.

QPPLab: A generally applicable software package for detecting, analyzing, and visualizing large-scale quasiperiodic spatiotemporal patterns (QPPs) of brain activity

One prominent brain dynamic process detected in functional neuroimaging data is large-scale quasi-periodic patterns (QPPs) which display spatiotemporal propagations along brain cortical gradients. QPP associates with the infraslow neural activity related to attention and arousal fluctuations and has been identified in both resting and task-evoked brains across various species. Several QPP detection and analysis tools were developed for distinct applications with study-specific parameter methods. This MATLAB package provides a simplified and user-friendly generally applicable toolbox for detecting, analyzing, and visualizing QPPs from fMRI timeseries of the brain. This paper describes the software functions and presents its ease of use on any brain datasets.

Metadata

[Table: see text].

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.

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.

Brain network integration, segregation and quasi-periodic activation and deactivation during tasks and rest

Previous studies have shown that a re-organization of the brain's functional connectome expressed in terms of network integration and segregation may play a pivotal role for brain function. However, it has been proven difficult to fully capture both processes independently in a single methodological framework. In this study, by starting from pair-wise assessments of instantaneous phase synchronization and community membership, we assemble spatiotemporally flexible networks that reflect changes in integration/segregation that occur at a spectrum of spatial as well as temporal scales. This is achieved by iteratively assembling smaller networks into larger units under the constraint that the smaller units should be internally integrated, i.e. belong to the same community. The assembled subnetworks can be partly overlapping and differ in size across time. Our results show that subnetwork integration and segregation occur simultaneously in the brain. During task performance, global changes in synchronization between networks arise that are tied to the underlying temporal design of the experiment. We show that a hallmark property of the dynamics of the brain's functional connectome is a presence of quasi-periodic patterns of network activation and deactivation, which during task performance becomes intertwined with the underlying temporal structure of the experimental paradigm. Additionally, we show that the degree of network integration throughout a n-back working memory task is correlated to performance.

Recent advances in psychoradiology

Psychiatry, as a field, lacks objective markers for diagnosis, progression, treatment planning, and prognosis, in part due to difficulties studying the brain in vivo, and diagnoses are based on self-reported symptoms and observation of patient behavior and cognition. Rapid advances in brain imaging techniques allow clinical investigators to noninvasively quantify brain features at the structural, functional, and molecular levels. Psychoradiology is an emerging discipline at the intersection of psychiatry and radiology. Psychoradiology applies medical imaging technologies to psychiatry and promises not only to improve insight into structural and functional brain abnormalities in patients with psychiatric disorders but also to have potential clinical utility. We searched for representative studies related to recent advances in psychoradiology through May 1, 2022, and conducted a selective review of 165 references, including 75 research articles. We summarize the novel dynamic imaging processing methods to model brain networks and present imaging genetics studies that reveal the relationship between various neuroimaging endophenotypes and genetic markers in psychiatric disorders. Furthermore, we survey recent advances in psychoradiology, with a focus on future psychiatric diagnostic approaches with dimensional analysis and a shift from group-level to individualized analysis. Finally, we examine the application of machine learning in psychoradiology studies and the potential of a novel option for brain stimulation treatment based on psychoradiological findings in precision medicine. Here, we provide a summary of recent advances in psychoradiology research, and we hope this review will help guide the practice of psychoradiology in the scientific and clinical fields.

An orderly sequence of autonomic and neural events at transient arousal changes

Resting-state functional magnetic resonance imaging (rsfMRI) allows the study of functional brain connectivity based on spatially structured variations in neuronal activity. Proper evaluation of connectivity requires removal of non-neural contributions to the fMRI signal, in particular hemodynamic changes associated with autonomic variability. Regression analysis based on autonomic indicator signals has been used for this purpose, but may be inadequate if neuronal and autonomic activities covary. To investigate this potential co-variation, we performed rsfMRI experiments while concurrently acquiring electroencephalography (EEG) and autonomic indicator signals, including heart rate, respiratory depth, and peripheral vascular tone. We identified a recurrent and systematic spatiotemporal pattern of fMRI (named as fMRI cascade), which features brief signal reductions in salience and default-mode networks and the thalamus, followed by a biphasic global change with a sensory-motor dominance. This fMRI cascade, which was mostly observed during eyes-closed condition, was accompanied by large EEG and autonomic changes indicative of arousal modulations. Importantly, the removal of the fMRI cascade dynamics from rsfMRI diminished its correlations with various signals. These results suggest that the rsfMRI correlations with various physiological and neural signals are not independent but arise, at least partly, from the fMRI cascades and associated neural and physiological changes at arousal modulations.

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.

Altered basal forebrain function during whole-brain network activity at pre- and early-plaque stages of Alzheimer's disease in TgF344-AD rats

BACKGROUND

Imbalanced synaptic transmission appears to be an early driver in Alzheimer's disease (AD) leading to brain network alterations. Early detection of altered synaptic transmission and insight into mechanisms causing early synaptic alterations would be valuable treatment strategies. This study aimed to investigate how whole-brain networks are influenced at pre- and early-plague stages of AD and if these manifestations are associated with concomitant cellular and synaptic deficits.  METHODS: To this end, we used an established AD rat model (TgF344-AD) and employed resting state functional MRI and quasi-periodic pattern (QPP) analysis, a method to detect recurrent spatiotemporal motifs of brain activity, in parallel with state-of-the-art immunohistochemistry in selected brain regions.

RESULTS

At the pre-plaque stage, QPPs in TgF344-AD rats showed decreased activity of the basal forebrain (BFB) and the default mode-like network. Histological analyses revealed increased astrocyte abundance restricted to the BFB, in the absence of amyloid plaques, tauopathy, and alterations in a number of cholinergic, gaba-ergic, and glutamatergic synapses. During the early-plaque stage, when mild amyloid-beta (Aβ) accumulation was observed in the cortex and hippocampus, QPPs in the TgF344-AD rats normalized suggesting the activation of compensatory mechanisms during this early disease progression period. Interestingly, astrogliosis observed in the BFB at the pre-plaque stage was absent at the early-plaque stage. Moreover, altered excitatory/inhibitory balance was observed in cortical regions belonging to the default mode-like network. In wild-type rats, at both time points, peak activity in the BFB preceded peak activity in other brain regions-indicating its modulatory role during QPPs. However, this pattern was eliminated in TgF344-AD suggesting that alterations in BFB-directed neuromodulation have a pronounced impact in network function in AD.

CONCLUSIONS

This study demonstrates the value of rsfMRI and advanced network analysis methods to detect early alterations in BFB function in AD, which could aid early diagnosis and intervention in AD. Restoring the global synaptic transmission, possibly by modulating astrogliosis in the BFB, might be a promising therapeutic strategy to restore brain network function and delay the onset of symptoms in AD.

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.

A parsimonious description of global functional brain organization in three spatiotemporal patterns

Resting-state functional magnetic resonance imaging (MRI) has yielded seemingly disparate insights into large-scale organization of the human brain. The brain's large-scale organization can be divided into two broad categories: zero-lag representations of functional connectivity structure and time-lag representations of traveling wave or propagation structure. In this study, we sought to unify observed phenomena across these two categories in the form of three low-frequency spatiotemporal patterns composed of a mixture of standing and traveling wave dynamics. We showed that a range of empirical phenomena, including functional connectivity gradients, the task-positive/task-negative anti-correlation pattern, the global signal, time-lag propagation patterns, the quasiperiodic pattern and the functional connectome network structure, are manifestations of these three spatiotemporal patterns. These patterns account for much of the global spatial structure that underlies functional connectivity analyses and unifies phenomena in resting-state functional MRI previously thought distinct.

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.

Latency structure of BOLD signals within white matter in resting-state fMRI

PURPOSE

Previous studies have demonstrated that BOLD signals in gray matter in resting-state functional MRI (RSfMRI) have variable time lags, representing apparent propagations of fMRI BOLD signals in gray matter. We complemented existing findings and explored the corresponding variations of signal latencies in white matter.

METHODS

We used data from the Brain Genomics Superstruct Project, consisting of 1412 subjects (both sexes included) and divided the dataset into ten equal groups to study both the patterns and reproducibility of latency estimates within white matter. We constructed latency matrices by computing cross-covariances between voxel pairs. We also applied a clustering analysis to identify functional networks within white matter, based on which latency analysis was also performed to investigate lead/lag relationship at network level. A dataset consisting of various sensory states (eyes closed, eyes open and eyes open with fixation) was also included to examine the relationship between latency structure and different states.

RESULTS

Projections of voxel latencies from the latency matrices were highly correlated (average Pearson correlation coefficient = 0.89) across the subgroups, confirming the reproducibility and structure of signal lags in white matter. Analysis of latencies within and between networks revealed a similar pattern of inter- and intra-network communication to that reported for gray matter. Moreover, a dominant direction, from inferior to superior regions, of BOLD signal propagation was revealed by higher resolution clustering. The variations of lag structure within white matter are associated with different sensory states.

CONCLUSIONS

These findings provide additional insight into the character and roles of white matter BOLD signals in brain functions.

Comparison of Resting-State Functional MRI Methods for Characterizing Brain Dynamics

Resting-state functional MRI (fMRI) exhibits time-varying patterns of functional connectivity. Several different analysis approaches have been developed for examining these resting-state dynamics including sliding window connectivity (SWC), phase synchrony (PS), co-activation pattern (CAP), and quasi-periodic patterns (QPP). Each of these approaches can be used to generate patterns of activity or inter-areal coordination which vary across time. The individual frames can then be clustered to produce temporal groupings commonly referred to as “brain states.” Several recent publications have investigated brain state alterations in clinical populations, typically using a single method for quantifying frame-wise functional connectivity. This study directly compares the results of k-means clustering in conjunction with three of these resting-state dynamics methods (SWC, CAP, and PS) and quantifies the brain state dynamics across several metrics using high resolution data from the human connectome project. Additionally, these three dynamics methods are compared by examining how the brain state characterizations vary during the repeated sequences of brain states identified by a fourth dynamic analysis method, QPP. The results indicate that the SWC, PS, and CAP methods differ in the clusters and trajectories they produce. A clear illustration of these differences is given by how each one results in a very different clustering profile for the 24s sequences explicitly identified by the QPP algorithm. PS clustering is sensitive to QPPs with the mid-point of most QPP sequences grouped into the same single cluster. CAPs are also highly sensitive to QPPs, separating each phase of the QPP sequences into different sets of clusters. SWC (60s window) is less sensitive to QPPs. While the QPPs are slightly more likely to occur during specific SWC clusters, the SWC clustering does not vary during the 24s QPP sequences, the goal of this work is to improve both the practical and theoretical understanding of different resting-state dynamics methods, thereby enabling investigators to better conceptualize and implement these tools for characterizing functional brain networks.

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.

Spatiotemporal trajectories in resting-state FMRI revealed by convolutional variational autoencoder

Recent resting-state fMRI studies have shown that brain activity exhibits temporal variations in functional connectivity by using various approaches including sliding window correlation, co-activation patterns, independent component analysis, quasi-periodic patterns, and hidden Markov models. These methods often model the brain activity as a discretized hopping among several brain states that are defined by the spatial configurations of network activity. However, the discretized states are merely a simplification of what is likely to be a continuous process, where each network evolves over time following its unique path. To model these characteristic spatiotemporal trajectories, we trained a variational autoencoder using rs-fMRI data and evaluated the spatiotemporal features of the latent variables obtained from the trained networks. Our results suggest that there are a relatively small number of approximately orthogonal whole-brain spatiotemporal patterns that capture the most prominent features of rs-fMRI data, which can serve as the building blocks to construct all possible spatiotemporal dynamics in resting state fMRI. These spatiotemporal patterns provide insight into how activity flows across the brain in concordance with known network structures and functional connectivity gradients.

Single-neuron firing cascades underlie global spontaneous brain events

The resting brain consumes enormous energy and shows highly organized spontaneous activity. To investigate how this activity is manifest among single neurons, we analyzed spiking discharges of ∼10,000 isolated cells recorded from multiple cortical and subcortical regions of the mouse brain during immobile rest. We found that firing of a significant proportion (∼70%) of neurons conformed to a ubiquitous, temporally sequenced cascade of spiking that was synchronized with global events and elapsed over timescales of 5 to 10 s. Across the brain, two intermixed populations of neurons supported orthogonal cascades. The relative phases of these cascades determined, at each moment, the response magnitude evoked by an external visual stimulus. Furthermore, the spiking of individual neurons embedded in these cascades was time locked to physiological indicators of arousal, including local field potential power, pupil diameter, and hippocampal ripples. These findings demonstrate that the large-scale coordination of low-frequency spontaneous activity, which is commonly observed in brain imaging and linked to arousal, sensory processing, and memory, is underpinned by sequential, large-scale temporal cascades of neuronal spiking across the brain.

Recent developments in representations of the connectome

Research into the human connectome (i.e., all connections in the human brain) with the use of resting state functional MRI has rapidly increased in popularity in recent years, especially with the growing availability of large-scale neuroimaging datasets. The goal of this review article is to describe innovations in functional connectome representations that have come about in the past 8 years, since the 2013 NeuroImage special issue on 'Mapping the Connectome'. In the period, research has shifted from group-level brain parcellations towards the characterization of the individualized connectome and of relationships between individual connectomic differences and behavioral/clinical variation. Achieving subject-specific accuracy in parcel boundaries while retaining cross-subject correspondence is challenging, and a variety of different approaches are being developed to meet this challenge, including improved alignment, improved noise reduction, and robust group-to-subject mapping approaches. Beyond the interest in the individualized connectome, new representations of the data are being studied to complement the traditional parcellated connectome representation (i.e., pairwise connections between distinct brain regions), such as methods that capture overlapping and smoothly varying patterns of connectivity ('gradients'). These different connectome representations offer complimentary insights into the inherent functional organization of the brain, but challenges for functional connectome research remain. Interpretability will be improved by future research towards gaining insights into the neural mechanisms underlying connectome observations obtained from functional MRI. Validation studies comparing different connectome representations are also needed to build consensus and confidence to proceed with clinical trials that may produce meaningful clinical translation of connectome insights.

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.

Brain Activity Fluctuations Propagate as Waves Traversing the Cortical Hierarchy

The brain exhibits highly organized patterns of spontaneous activity as measured by resting-state functional magnetic resonance imaging (fMRI) fluctuations that are being widely used to assess the brain's functional connectivity. Some evidence suggests that spatiotemporally coherent waves are a core feature of spontaneous activity that shapes functional connectivity, although this has been difficult to establish using fMRI given the temporal constraints of the hemodynamic signal. Here, we investigated the structure of spontaneous waves in human fMRI and monkey electrocorticography. In both species, we found clear, repeatable, and directionally constrained activity waves coursed along a spatial axis approximately representing cortical hierarchical organization. These cortical propagations were closely associated with activity changes in distinct subcortical structures, particularly those related to arousal regulation, and modulated across different states of vigilance. The findings demonstrate a neural origin of spatiotemporal fMRI wave propagation at rest and link it to the principal gradient of resting-state fMRI connectivity.

Global waves synchronize the brain's functional systems with fluctuating arousal

We propose and empirically support a parsimonious account of intrinsic, brain-wide spatiotemporal organization arising from traveling waves linked to arousal. We hypothesize that these waves are the predominant physiological process reflected in spontaneous functional magnetic resonance imaging (fMRI) signal fluctuations. The correlation structure ("functional connectivity") of these fluctuations recapitulates the large-scale functional organization of the brain. However, a unifying physiological account of this structure has so far been lacking. Here, using fMRI in humans, we show that ongoing arousal fluctuations are associated with global waves of activity that slowly propagate in parallel throughout the neocortex, thalamus, striatum, and cerebellum. We show that these waves can parsimoniously account for many features of spontaneous fMRI signal fluctuations, including topographically organized functional connectivity. Last, we demonstrate similar, cortex-wide propagation of neural activity measured with electrocorticography in macaques. These findings suggest that traveling waves spatiotemporally pattern brain-wide excitability in relation to arousal.

Propagating patterns of intrinsic activity along macroscale gradients coordinate functional connections across the whole brain

The intrinsic activity of the human brain, observed with resting-state fMRI (rsfMRI) and functional connectivity, exhibits macroscale spatial organization such as functional networks and gradients. Dynamic analysis techniques have shown that functional connectivity is a mere summary of time-varying patterns with distinct spatial and temporal characteristics. A better understanding of these patterns might provide insight into aspects of the brain's intrinsic activity that cannot be inferred by functional connectivity or the spatial maps derived from it, such as functional networks and gradients. Here, we describe three spatiotemporal patterns of coordinated activity across the whole brain obtained by averaging similar ~20-second-long segments of rsfMRI timeseries. In each of these patterns, activity propagates along a particular macroscale functional gradient, simultaneously across the cerebral cortex and in most other brain regions. In some regions, like the thalamus, the propagation suggests previously-undescribed gradients. The coordinated activity across areas is consistent with known tract-based connections, and nuanced differences in the timing of peak activity between regions point to plausible driving mechanisms. The magnitude of correlation within and particularly between functional networks is remarkably diminished when these patterns are regressed from the rsfMRI timeseries, a quantitative demonstration of their significant role in functional connectivity. Taken together, our results suggest that a few recurring patterns of propagating intrinsic activity along macroscale gradients give rise to and coordinate functional connections across the whole brain.

Origin of the Time Lag Phenomenon and the Global Signal in Resting-State fMRI

The global mean signal of resting-state fMRI (rs-fMRI) shows a characteristic spatiotemporal pattern that is closely related to the pattern of vascular perfusion. Although being increasingly adopted in the mapping of the flow of neural activity, the mechanism that gives rise to the BOLD signal time lag remains controversial. In the present study, we compared the time lag of the global mean signal with those of the local network components obtained by applying temporal independent component analysis to the resting-state fMRI data, as well as by using simultaneous wide-field visual stimulation, and demonstrated that the time lag patterns are highly similar across all types of data. These results suggest that the time lag of the rs-fMRI signal reflects the local variance of the hemodynamic responses rather than the arrival or transit time of the stimulus, whether the trigger is neuronal or non-neuronal in origin as long as it is mediated by local hemodynamic responses. Examinations of the internal carotid artery signal further confirmed that the arterial signal is tightly inversely coupled with the global mean signal in accordance with previous studies, presumably reflecting the blood flow or blood pressure changes that are occurring almost simultaneously in the internal carotid artery and the cerebral pial/capillary arteries, within the low-frequency component in human rs-fMRI.

Atypical Inter-Network Deactivation Associated With the Posterior Default-Mode Network in Autism Spectrum Disorder

Previous studies have suggested that atypical deactivation of functional brain networks contributes to the complex cognitive and behavioral profile associated with autism spectrum disorder (ASD). However, these studies have not considered the temporal dynamics of deactivation mechanisms between the networks. In this study, we examined (a) mutual deactivation and (b) mutual activation-deactivation (i.e., anticorrelated) time-lag patterns between resting-state networks (RSNs) in young adults with ASD (n = 20) and controls (n = 20) by applying the recently defined dynamic lag analysis (DLA) method, which measures time-lag variations peak-by-peak between the networks. In order to achieve temporally accurate lag patterns, the brain imaging data was acquired with a fast functional magnetic resonance imaging (fMRI) sequence (TR = 100 ms). Group-level independent component analysis was used to identify 16 RSNs for the DLA. We found altered mutual deactivation timings in ASD in (a) three of the deactivated and (b) two of the transiently anticorrelated (activated-deactivated) RSN pairs, which survived the strict threshold for significance of surrogate data. Of the significant RSN pairs, 80% included the posterior default-mode network (DMN). We propose that temporally altered deactivation mechanisms, including timings and directionality, between the posterior DMN and RSNs mediating processing of socially relevant information may contribute to the ASD phenotype. LAY SUMMARY: To understand autistic traits on a neural level, we examined temporal fluctuations in information flow between brain regions in young adults with autism spectrum disorder (ASD) and controls. We used a fast neuroimaging procedure to investigate deactivation mechanisms between brain regions. We found that timings and directionality of communication between certain brain regions were temporally altered in ASD, suggesting atypical deactivation mechanisms associated with the posterior default-mode network.

Dynamic Brain Connectivity in Resting State Functional MR Imaging

Dynamic functional connectivity adds another dimension to resting-state functional MR imaging analysis. In recent years, dynamic functional connectivity has been increasingly used in resting-state functional MR imaging, and several studies have demonstrated that dynamic functional connectivity patterns correlate with different physiologic and pathologic brain states. In fact, evidence suggests that dynamic functional connectivity is a more sensitive marker than static functional connectivity; therefore, it might be a promising tool to add to clinical functional neuroimaging. This article provides a broad overview of dynamic functional connectivity and reviews its general principles, techniques, and potential clinical applications.

Resting Brain Fluctuations Are Intrinsically Coupled to Visual Response Dynamics

How do intrinsic brain dynamics interact with processing of external sensory stimuli? We sought new insights using functional magnetic resonance imaging to track spatiotemporal activity patterns at the whole brain level in lightly anesthetized mice, during both resting conditions and visual stimulation trials. Our results provide evidence that quasiperiodic patterns (QPPs) are the most prominent component of mouse resting brain dynamics. These QPPs captured the temporal alignment of anticorrelation between the default mode (DMN)- and task-positive (TPN)-like networks, with global brain fluctuations, and activity in neuromodulatory nuclei of the reticular formation. Specifically, the phase of QPPs prior to stimulation could significantly stratify subsequent visual response magnitude, suggesting QPPs relate to brain state fluctuations. This is the first observation in mice that dynamics of the DMN- and TPN-like networks, and particularly their anticorrelation, capture a brain state dynamic that affects sensory processing. Interestingly, QPPs also displayed transient onset response properties during visual stimulation, which covaried with deactivations in the reticular formation. We conclude that QPPs appear to capture a brain state fluctuation that may be orchestrated through neuromodulation. Our findings provide new frontiers to understand the neural processes that shape functional brain states and modulate sensory input processing.

15 Years MR-encephalography

Objective

This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere ‘Gedankenexperiment’ in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions.

Methods

The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and—most recently—the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality.

Applications

Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility.

Discussion

MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

Decoding Analysis of Alpha Oscillation Networks on Maintaining Driver Alertness

The countermeasure of driver fatigue is valuable for reducing the risk of accidents caused by vigilance failure during prolonged driving. Listening to the radio (RADIO) has been proven to be a relatively effective “in-car” countermeasure. However, the connectivity analysis, which can be used to investigate its alerting effect, is subject to the issue of signal mixing. In this study, we propose a novel framework based on clustering and entropy to improve the performance of the connectivity analysis to reveal the effect of RADIO to maintain driver alertness. Regardless of reducing signal mixing, we introduce clustering algorithm to classify the functional connections with their nodes into different categories to mine the effective information of the alerting effect. Differential entropy (DE) is employed to measure the information content in different brain regions after clustering. Compared with the Louvain-based community detection method, the proposed method shows more superior ability to present RADIO effectin confused functional connection matrices. Our experimental results reveal that the active connection clusters distinguished by the proposed method gradually move from frontal region to parieto-occipital regionwith the progress of fatigue, consistent with the alpha energy changes in these two brain areas. The active class of the clusters in parieto-occipital region significantly decreases and the most active clusters remain in the frontal region when RADIO is taken. The estimation results of DE confirm the significant change (p < 0.05) of information content due to the cluster movements. Hence, preventing the movement of the active clusters from frontal region to parieto-occipital region may correlate with maintaining driver alertness. The revelation of alerting effect is helpful for the targeted upgrade of fatigue countermeasures.

Aberrant static and dynamic functional patterns of frontoparietal control network in antipsychotic-naïve first-episode psychosis subjects

Psychotic disorders are disabling clinical syndromes characterized by widespread alterations in cortical information processing. Disruption of frontoparietal network (FPN) connectivity has emerged as a common footprint across the psychosis spectrum. Our goal was to characterize the static and dynamic resting‐state functional connectivity (FC) of the FPN in antipsychotic‐naïve first‐episode psychosis (FEP) subjects. We compared the static FC of the FPN in 40 FEP and 40 healthy control (HC) subjects, matched on age, sex, and socioeconomic status. To study the dynamic FC, we measured quasiperiodic patterns (QPPs) that consist of infraslow spatioemporal patterns embedded in the blood oxygen level‐dependent signal that repeats over time, exhibiting alternation of high and low activity. Relative to HC, we found functional hypoconnectivity between the right middle frontal gyrus and the right middle temporal gyrus, as well as the left inferior temporal gyrus and the left inferior parietal gyrus in FEP (p < .05, false discovery rate corrected). The correlation of the QPP with all functional scans was significantly stronger for FEP compared to HC, suggesting a greater impact of the QPPs to intrinsic brain activity in psychotic population. Regressing the QPP from the functional scans erased all significant group differences in static FC, suggesting that abnormal connectivity in FEP could result from altered QPP. Our study supports that alterations of cortical information processing are not a function of psychotic chronicity or antipsychotic medication exposure and may be regarded as trait specific. In addition, static connectivity abnormality may be partly related to altered brain network temporal dynamics.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

Low Frequency Phase-locking of Brain Signals Contribute to Efficient Face Recognition

Low frequency phase synchronization is an essential mechanism of information communication among brain regions. In the infra-slow frequency range (<0.1 Hz), inter-regional phase lag is of importance for brain function (e.g., anti-phase between the default mode network and task positive network). However, the role of phase lag in cognitive processing remains unclear. Based on the frequency tagging experimental paradigm and functional magnetic resonance imaging (fMRI) technique, we investigated inter-regional phase lag and phase coherence using a face recognition task (n = 30, 15 males/15 females). Phase coherence within the face processing system was significantly increased during task state, highlighting the importance of regular inter-regional phase relationship for face recognition. Moreover, results showed decreased phase lag within the core and extended face areas (face processing system) and increased phase lag between the face processing system and frontoparietal network, indicating a reorganization of inter-regional relationships of the two systems. Inter-regional phase lag was modulated by the task at ascending and descending phases of the fMRI signal, suggesting a phase-dependent inter-regional relationship. Furthermore, phase lags between visual cortex and amygdala and between visual cortex and motor area were positively related to reaction time, indicating better task performance depends on both rapid emotional detection pathway and visual-motor pathway. Overall, inter-regional phase synchronization in the infra-slow frequency range is of important for effective information communication and cognitive performance.

Global vs. Network-Specific Regulations as the Source of Intrinsic Coactivations in Resting-State Networks

Spontaneous neural activities are endowed with specific patterning characterized by synchronizations within functionally relevant distant regions that are termed as resting-state networks (RSNs). Although the mechanisms that organize the large-scale neural systems are still largely unknown, recent studies have proposed a hypothesis that network-specific coactivations indeed emerge as the result of globally propagating neural activities with specific paths of transmission. However, the extent to which such a centralized global regulation, rather than network-specific control, contributes to the RSN synchronization remains unknown. In the present study, we investigated the contribution from each mechanism by directly identifying the global as well as local component of resting-state functional MRI (fMRI) data provided by human connectome project, using temporal independent component analysis (ICA). Based on the spatial distribution pattern, each ICA component was classified as global or local. Time lag mapping of each IC revealed several paths of global or semi-global propagations that are partially overlapping yet spatially distinct to each other. Consistent with previous studies, the time lag of global oscillation, although being less spatially homogenous than what was assumed to be, contributed to the RSN synchronization. However, an equivalent contribution was also shown on the part of the more locally confined activities that are independent to each other. While allowing the view that network-specific coactivation occurs as part of the sequences of global neural activities, these results further confirm an equally important role of the network-specific regulation for its coactivation, regardless of whether vascular artifacts contaminate the global component in fMRI measures.