Transient neuronal coactivations embedded in globally propagating waves underlie resting-state functional connectivity

Dynamic resting state fMRI analysis in mice reveals a set of Quasi-Periodic Patterns and illustrates their relationship with the global signal

Time-resolved 'dynamic' over whole-period 'static' analysis of low frequency (LF) blood-oxygen level dependent (BOLD) fluctuations provides many additional insights into the macroscale organization and dynamics of neural activity. Although there has been considerable advancement in the development of mouse resting state fMRI (rsfMRI), very little remains known about its dynamic repertoire. Here, we report for the first time the detection of a set of recurring spatiotemporal Quasi-Periodic Patterns (QPPs) in mice, which show spatial similarity with known resting state networks. Furthermore, we establish a close relationship between several of these patterns and the global signal. We acquired high temporal rsfMRI scans under conditions of low (LA) and high (HA) medetomidine-isoflurane anesthesia. We then employed the algorithm developed by Majeed et al. (2011), previously applied in rats and humans, which detects and averages recurring spatiotemporal patterns in the LF BOLD signal. One type of observed patterns in mice was highly similar to those originally observed in rats, displaying propagation from lateral to medial cortical regions, which suggestively pertain to a mouse Task-Positive like network (TPN) and Default Mode like network (DMN). Other QPPs showed more widespread or striatal involvement and were no longer detected after global signal regression (GSR). This was further supported by diminished detection of subcortical dynamics after GSR, with cortical dynamics predominating. Observed QPPs were both qualitatively and quantitatively determined to be consistent across both anesthesia conditions, with GSR producing the same outcome. Under LA, QPPs were consistently detected at both group and single subject level. Under HA, consistency and pattern occurrence rate decreased, whilst cortical contribution to the patterns diminished. These findings confirm the robustness of QPPs across species and demonstrate a new approach to study mouse LF BOLD spatiotemporal dynamics and mechanisms underlying functional connectivity. The observed impact of GSR on QPPs might help better comprehend its controversial role in conventional resting state studies. Finally, consistent detection of QPPs at single subject level under LA promises a step forward towards more reliable mouse rsfMRI and further confirms the importance of selecting an optimal anesthesia regime.

Neurogenetic profiles delineate large-scale connectivity dynamics of the human brain

Experimental and modeling work of neural activity has described recurrent and attractor dynamic patterns in cerebral microcircuits. However, it is still poorly understood whether similar dynamic principles exist or can be generalizable to the large-scale level. Here, we applied dynamic graph theory-based analyses to evaluate the dynamic streams of whole-brain functional connectivity over time across cognitive states. Dynamic connectivity in local networks is located in attentional areas during tasks and primary sensory areas during rest states, and dynamic connectivity in distributed networks converges in the default mode network (DMN) in both task and rest states. Importantly, we find that distinctive dynamic connectivity patterns are spatially associated with Allen Human Brain Atlas genetic transcription levels of synaptic long-term potentiation and long-term depression-related genes. Our findings support the neurobiological basis of large-scale attractor-like dynamics in the heteromodal cortex within the DMN, irrespective of cognitive state.

Decreased BOLD fluctuations in lateral temporal cortices of premature born adults

Lasting volume reductions in subcortical and temporal-insular cortices after premature birth suggest altered ongoing activity in these areas. We hypothesized altered fluctuations in ongoing neural excitability and activity, as measured by slowly fluctuating blood oxygenation of resting-state functional MRI (rs-fMRI), in premature born adults, with altered fluctuations being linked with underlying brain volume reductions. To investigate this hypothesis, 94 very preterm/very low birth weight (VP/VLBW) and 92 full-term born young adults underwent structural and rs-fMRI data acquisition with voxel-based morphometry and amplitude of low-frequency fluctuations (ALFF) as main outcome measure. In VP/VLBW adults, ALFF was reduced in lateral temporal cortices, and this reduction was positively associated with lower birth weight. Regions of reduced ALFF overlapped with reduced brain volume. On the one hand, ALFF reduction remained after controlling for volume loss, supporting the functional nature of ALFF reductions. On the other hand, ALFF decreases were positively associated with underlying brain volume loss, indicating a relation between structural and functional changes. Furthermore, within the VP/VLBW group, reduced ALFF was associated with reduced IQ, indicating the behavioral relevance of ALFF decreases in temporal cortices. These results demonstrate long-term impact of premature birth on ongoing BOLD fluctuations in lateral temporal cortices, which are linked with brain volume reductions. Data suggest permanently reduced fluctuations in ongoing neural excitability and activity in structurally altered lateral temporal cortices after premature birth.

Longitudinal observations using simultaneous fMRI, multiple channel electrophysiology recording, and chemical microiontophoresis in the rat brain

BACKGROUND

fMRI blood oxygenation level-dependent (BOLD) signal has been widely used as a surrogate for neural activity. However, interpreting differences in BOLD fMRI based on underlying neuronal activity remains a challenge. Concurrent rsMRI data collection and electrophysiological recording in combination with microiontophoretically injected modulatory chemicals allows for improved understanding of the relationship between resting state BOLD and neuronal activity.

NEW METHODS

Simultaneous fMRI, multi-channel intracortical electrophysiology and focal pharmacological manipulation data to be acquired longitudinally in rats for up to 2 months. Our artifact replacing technique is optimized for combined LFP and rsMRI data collection.

RESULTS

Intracortical implantation of a multichannel microelectrode array resulted in minimal distortion and signal loss in fMRI images inside a 9.4T MRI scanner. rsMRI-induced electrophysiology artifacts were replaced using an in-house developed algorithm. Microinjection of AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) enhanced dopaminergic neuronal activity in the ventral tegmental area (VTA) and altered LFP signal and fMRI functional connectivity in the striatum.

COMPARISONS WITH EXISTING METHOD(S)

Nanomanufacturing advances permit the production of MRI-compatible microelectrode arrays (with 16 or more channels), extending research beyond conventional methods limited to fewer channels. Our method permits longitudinal data collection of LFP and rsMRI and our algorithm effectively detects and replaces fMRI-induced electrophysiological noise, permitting rsMRI data collection concomitant with LFP recordings.

CONCLUSIONS

Our model consists of longitudinal concurrent fMRI and multichannel intracortical electrophysiological recording during microinjection of pharmacological agents to modulate neural activity in the rat brain. We used commercial micro-electrodes and recording system and can be readily generalized to other labs.

Theta and Alpha Oscillations Are Traveling Waves in the Human Neocortex

Human cognition requires the coordination of neural activity across widespread brain networks. Here, we describe a new mechanism for large-scale coordination in the human brain: traveling waves of theta and alpha oscillations. Examining direct brain recordings from neurosurgical patients performing a memory task, we found contiguous clusters of cortex in individual patients with oscillations at specific frequencies within 2 to 15 Hz. These oscillatory clusters displayed spatial phase gradients, indicating that they formed traveling waves that propagated at ∼0.25-0.75 m/s. Traveling waves were relevant behaviorally because their propagation correlated with task events and was more consistent when subjects performed the task well. Human traveling theta and alpha waves can be modeled by a network of coupled oscillators because the direction of wave propagation correlated with the spatial orientation of local frequency gradients. Our findings suggest that oscillations support brain connectivity by organizing neural processes across space and time.

Short timescale abnormalities in the states of spontaneous synchrony in the functional neural networks in Alzheimer's disease

Alzheimer's disease (AD) is a prevalent neurodegenerative condition that can lead to severe cognitive and functional deterioration. Functional magnetic resonance imaging (fMRI) revealed abnormalities in AD in intrinsic synchronization between spatially separate regions in the so-called default mode network (DMN) of the brain. To understand the relationship between this disruption in large-scale synchrony and the cognitive impairment in AD, it is critical to determine whether and how the deficit in the low frequency hemodynamic fluctuations recorded by fMRI translates to much faster timescales of memory and other cognitive processes. The present study employed magnetoencephalography (MEG) and a Hidden Markov Model (HMM) approach to estimate spontaneous synchrony variations in the functional neural networks with high temporal resolution. In a group of cognitively healthy (CH) older adults, we found transient (mean duration of 150-250 ms) network activity states, which were visited in a rapid succession, and were characterized by spatially coordinated changes in the amplitude of source-localized electrophysiological oscillations. The inferred states were similar to those previously observed in younger participants using MEG, and the estimated cortical source distributions of the state-specific activity resembled the classic functional neural networks, such as the DMN. In patients with AD, inferred network states were different from those of the CH group in short-scale timing and oscillatory features. The state of increased oscillatory amplitudes in the regions overlapping the DMN was visited less often in AD and for shorter periods of time, suggesting that spontaneous synchronization in this network was less likely and less stable in the patients. During the visits to this state, in some DMN nodes, the amplitude change in the higher-frequency (8-30 Hz) oscillations was less robust in the AD than CH group. These findings highlight relevance of studying short-scale temporal evolution of spontaneous activity in functional neural networks to understanding the AD pathophysiology. Capacity of flexible intrinsic synchronization in the DMN may be crucial for memory and other higher cognitive functions. Our analysis yielded metrics that quantify distinct features of the neural synchrony disorder in AD and may offer sensitive indicators of the neural network health for future investigations.

Time-Resolved Resting-State Functional Magnetic Resonance Imaging Analysis: Current Status, Challenges, and New Directions

Time-resolved analysis of resting-state functional magnetic resonance imaging (rs-fMRI) data allows researchers to extract more information about brain function than traditional functional connectivity analysis, yet a number of challenges in data analysis and interpretation remain. This article briefly summarizes common methods for time-resolved analysis and presents some of the pressing issues and opportunities in the field. From there, the discussion moves to interpretation of the network dynamics observed with rs-fMRI and the role that rs-fMRI can play in elucidating the large-scale organization of brain activity.

Spontaneous Infra-slow Brain Activity Has Unique Spatiotemporal Dynamics and Laminar Structure

Systems-level organization in spontaneous infra-slow (<0.1Hz) brain activity, measured using blood oxygen signals in fMRI and optical imaging, has become a major theme in the study of neural function in both humans and animal models. Yet the neurophysiological basis of infra-slow activity (ISA) remains unresolved. In particular, is ISA a distinct physiological process, or is it a low-frequency analog of faster neural activity? Here, using whole-cortex calcium/hemoglobin imaging in mice, we show that ISA in each of these modalities travels through the cortex along stereotypical spatiotemporal trajectories that are state dependent (wake versus anesthesia) and distinct from trajectories in delta (1-4 Hz) activity. Moreover, mouse laminar electrophysiology reveals that ISA travels through specific cortical layers and is organized into unique cross-laminar temporal dynamics that are different from higher frequency local field potential activity. These findings suggest that ISA is a distinct neurophysiological process that is reflected in fMRI blood oxygen signals.

Fiber-optic implant for simultaneous fluorescence-based calcium recordings and BOLD fMRI in mice

Despite the growing popularity of blood oxygen level-dependent (BOLD) functional MRI (fMRI), understanding of its underlying principles is still limited. This protocol describes a technique for simultaneous measurement of neural activity using fluorescent calcium indicators together with the corresponding hemodynamic BOLD fMRI response in the mouse brain. Our early work using small-molecule fluorophores in rats gave encouraging results but was limited to acute measurements using synthetic dyes. Our latest procedure combines fMRI with optical detection of cell-type-specific virally delivered GCaMP6, a genetically encoded calcium indicator (GECI). GCaMP6 fluorescence, which increases upon calcium binding, is collected by a chronically implanted optical fiber, allowing longitudinal studies in mice. The chronic implant, placed horizontally on the skull, has an angulated tip that reflects light into the brain and is connected via fiber optics to a remote optical setup. The technique allows access to the neocortex and does not require adaptations of commercial MRI hardware. The hybrid approach permits fiber-optic calcium recordings with simultaneous artifact-free BOLD fMRI with full brain coverage and 1-s temporal resolution using standard gradient-echo echo-planar imaging (GE-EPI) sequences. The method provides robust, cell-type-specific readouts to link neural activity to BOLD signals, as emonstrated for task-free ('resting-state') conditions and in response to hind-paw stimulation. These results highlight the power of fiber photometry combined with fMRI, which we aim to further advance in this protocol. The approach can be easily adapted to study other molecular processes using suitable fluorescent indicators.

Distinctive Correspondence Between Separable Visual Attention Functions and Intrinsic Brain Networks

Separable visual attention functions are assumed to rely on distinct but interacting neural mechanisms. Bundesen's “theory of visual attention” (TVA) allows the mathematical estimation of independent parameters that characterize individuals' visual attentional capacity (i.e., visual processing speed and visual short-term memory storage capacity) and selectivity functions (i.e., top-down control and spatial laterality). However, it is unclear whether these parameters distinctively map onto different brain networks obtained from intrinsic functional connectivity, which organizes slowly fluctuating ongoing brain activity. In our study, 31 demographically homogeneous healthy young participants performed whole- and partial-report tasks and underwent resting-state functional magnetic resonance imaging (rs-fMRI). Report accuracy was modeled using TVA to estimate, individually, the four TVA parameters. Networks encompassing cortical areas relevant for visual attention were derived from independent component analysis of rs-fMRI data: visual, executive control, right and left frontoparietal, and ventral and dorsal attention networks. Two TVA parameters were mapped on particular functional networks. First, participants with higher (vs. lower) visual processing speed showed lower functional connectivity within the ventral attention network. Second, participants with more (vs. less) efficient top-down control showed higher functional connectivity within the dorsal attention network and lower functional connectivity within the visual network. Additionally, higher performance was associated with higher functional connectivity between networks: specifically, between the ventral attention and right frontoparietal networks for visual processing speed, and between the visual and executive control networks for top-down control. The higher inter-network functional connectivity was related to lower intra-network connectivity. These results demonstrate that separable visual attention parameters that are assumed to constitute relatively stable traits correspond distinctly to the functional connectivity both within and between particular functional networks. This implies that individual differences in basic attention functions are represented by differences in the coherence of slowly fluctuating brain activity.

Ridding fMRI data of motion-related influences: Removal of signals with distinct spatial and physical bases in multiecho data

"Functional connectivity" techniques are commonplace tools for studying brain organization. A critical element of these analyses is to distinguish variance due to neurobiological signals from variance due to nonneurobiological signals. Multiecho fMRI techniques are a promising means for making such distinctions based on signal decay properties. Here, we report that multiecho fMRI techniques enable excellent removal of certain kinds of artifactual variance, namely, spatially focal artifacts due to motion. By removing these artifacts, multiecho techniques reveal frequent, large-amplitude blood oxygen level-dependent (BOLD) signal changes present across all gray matter that are also linked to motion. These whole-brain BOLD signals could reflect widespread neural processes or other processes, such as alterations in blood partial pressure of carbon dioxide (pCO₂) due to ventilation changes. By acquiring multiecho data while monitoring breathing, we demonstrate that whole-brain BOLD signals in the resting state are often caused by changes in breathing that co-occur with head motion. These widespread respiratory fMRI signals cannot be isolated from neurobiological signals by multiecho techniques because they occur via the same BOLD mechanism. Respiratory signals must therefore be removed by some other technique to isolate neurobiological covariance in fMRI time series. Several methods for removing global artifacts are demonstrated and compared, and were found to yield fMRI time series essentially free of motion-related influences. These results identify two kinds of motion-associated fMRI variance, with different physical mechanisms and spatial profiles, each of which strongly and differentially influences functional connectivity patterns. Distance-dependent patterns in covariance are nearly entirely attributable to non-BOLD artifacts.

Stereotypical modulations in dynamic functional connectivity explained by changes in BOLD variance

Spontaneous activity measured in human subject under the absence of any task exhibits complex patterns of correlation that largely correspond to large-scale functional topographies obtained with a wide variety of cognitive and perceptual tasks. These "resting state networks" (RSNs) fluctuate over time, forming and dissolving on the scale of seconds to minutes. While these fluctuations, most prominently those of the default mode network, have been linked to cognitive function, it remains unclear whether they result from random noise or whether they index a nonstationary process which could be described as state switching. In this study, we use a sliding windows-approach to relate temporal dynamics of RSNs to global modulations in correlation and BOLD variance. We compare empirical data, phase-randomized surrogate data, and data simulated with a stationary model. We find that RSN time courses exhibit a large amount of coactivation in all three cases, and that the modulations in their activity are closely linked to global dynamics of the underlying BOLD signal. We find that many properties of the observed fluctuations in FC and BOLD, including their ranges and their correlations amongst each other, are explained by fluctuations around the average FC structure. However, we also report some interesting characteristics that clearly support nonstationary features in the data. In particular, we find that the brain spends more time in the troughs of modulations than can be expected from stationary dynamics.

Identifying and characterizing systematic temporally-lagged BOLD artifacts

Residual noise in the BOLD signal remains problematic for fMRI - particularly for techniques such as functional connectivity, where findings can be spuriously influenced by noise sources that can covary with individual differences. Many such potential noise sources - for instance, motion and respiration - can have a temporally lagged effect on the BOLD signal. Thus, here we present a tool for assessing residual lagged structure in the BOLD signal that is associated with nuisance signals, using a construction similar to a peri-event time histogram. Using this method, we find that framewise displacements - both large and very small - were followed by structured, prolonged, and global changes in the BOLD signal that depend on the magnitude of the preceding displacement and extend for tens of seconds. This residual lagged BOLD structure was consistent across datasets, and independently predicted considerable variance in the global cortical signal (as much as 30-40% in some subjects). Mean functional connectivity estimates varied similarly as a function of displacements occurring many seconds in the past, even after strict censoring. Similar patterns of residual lagged BOLD structure were apparent following respiratory fluctuations (which covaried with framewise displacements), implicating respiration as one likely mechanism underlying the displacement-linked structure observed. Global signal regression largely attenuates this artifactual structure. These findings suggest the need for caution in interpreting results of individual difference studies where noise sources might covary with the individual differences of interest, and highlight the need for further development of preprocessing techniques for mitigating such structure in a more nuanced and targeted manner.

Macroscale variation in resting-state neuronal activity and connectivity assessed by simultaneous calcium imaging, hemodynamic imaging and electrophysiology

Functional imaging of spontaneous activity continues to play an important role in the field of connectomics. The most common imaging signal used for these experiments is the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal, but how this signal relates to spontaneous neuronal activity remains incompletely understood. Genetically encoded calcium indicators represent a promising tool to study this problem, as they can provide brain-wide measurements of neuronal activity compared to point measurements afforded by electrophysiological recordings. However, the relationship between the calcium signal and neurophysiological parameters at the mesoscopic scale requires further systematic characterization. Therefore, we collected simultaneous resting-state measurements of electrophysiology, along with calcium and hemodynamic imaging, in lightly anesthetized mice to investigate two aims. First, we examined the relationship between each imaging signal and the simultaneously recorded electrophysiological signal in a single brain region, finding that both signals are better correlated with multi-unit activity compared to local field potentials, with the calcium signal possessing greater signal-to-noise ratio and regional specificity. Second, we used the resting-state imaging data to model the relationship between the calcium and hemodynamic signals across the brain. We found that this relationship varied across brain regions in a way that is consistent across animals, with delays increasing by600 ms towards posterior cortical regions. Furthermore, while overall functional connectivity (FC) measured by the hemodynamic signal is significantly correlated with FC measured by calcium, the two estimates were found to be significantly different. We hypothesize that these differences arise at least in part from the observed regional variation in the hemodynamic response. In total, this work highlights some of the caveats needed in interpreting hemodynamic-based measurements of FC, as well as the need for improved modeling methods to reduce this potential source of bias.

How networks communicate: propagation patterns in spontaneous brain activity

Initially regarded as 'noise', spontaneous (intrinsic) activity accounts for a large portion of the brain's metabolic cost. Moreover, it is now widely known that infra-slow (less than 0.1 Hz) spontaneous activity, measured using resting state functional magnetic resonance imaging of the blood oxygen level-dependent (BOLD) signal, is correlated within functionally defined resting state networks (RSNs). However, despite these advances, the temporal organization of spontaneous BOLD fluctuations has remained elusive. By studying temporal lags in the resting state BOLD signal, we have recently shown that spontaneous BOLD fluctuations consist of remarkably reproducible patterns of whole brain propagation. Embedded in these propagation patterns are unidirectional 'motifs' which, in turn, give rise to RSNs. Additionally, propagation patterns are markedly altered as a function of state, whether physiological or pathological. Understanding such propagation patterns will likely yield deeper insights into the role of spontaneous activity in brain function in health and disease.This article is part of the themed issue 'Interpreting blood oxygen level-dependent: a dialogue between cognitive and cellular neuroscience'.

Making connections in the brain

Simultaneous measurements of neuronal activity and fMRI signals in the rat brain have shed new light on the origins of resting-state fMRI connectivity networks.

Functional connectivity structure of cortical calcium dynamics in anesthetized and awake mice

The interplay between hemodynamic-based markers of cortical activity (e.g. fMRI and optical intrinsic signal imaging), which are an indirect and relatively slow report of neural activity, and underlying synaptic electrical and metabolic activity through neurovascular coupling is a topic of ongoing research and debate. As application of resting state functional connectivity measures is extended further into topics such as brain development, aging and disease, the importance of understanding the fundamental physiological basis for functional connectivity will grow. Here we extend functional connectivity analysis from hemodynamic- to calcium-based imaging. Transgenic mice (n = 7) expressing a fluorescent calcium indicator (GCaMP6) driven by the Thy1 promoter in glutamatergic neurons were imaged transcranially in both anesthetized (using ketamine/xylazine) and awake states. Sequential LED illumination (λ = 454, 523, 595, 640nm) enabled concurrent imaging of both GCaMP6 fluorescence emission (corrected for hemoglobin absorption) and hemodynamics. Functional connectivity network maps were constructed for infraslow (0.009–0.08Hz), intermediate (0.08–0.4Hz), and high (0.4–4.0Hz) frequency bands. At infraslow and intermediate frequencies, commonly used in BOLD fMRI and fcOIS studies of functional connectivity and implicated in neurovascular coupling mechanisms, GCaMP6 and HbO₂ functional connectivity structures were in high agreement, both qualitatively and also quantitatively through a measure of spatial similarity. The spontaneous dynamics of both contrasts had the highest correlation when the GCaMP6 signal was delayed with a ~0.6–1.5s temporal offset. Within the higher-frequency delta band, sensitive to slow wave sleep oscillations in non-REM sleep and anesthesia, we evaluate the speed with which the connectivity analysis stabilized and found that the functional connectivity maps captured putative network structure within time window lengths as short as 30 seconds. Homotopic GCaMP6 functional connectivity maps at 0.4–4.0Hz in the anesthetized states show a striking correlated and anti-correlated structure along the anterior to posterior axis. This structure is potentially explained in part by observed propagation of delta-band activity from frontal somatomotor regions to visuoparietal areas. During awake imaging, this spatio-temporal quality is altered, and a more complex and detailed functional connectivity structure is observed. The combined calcium/hemoglobin imaging technique described here will enable the dissociation of changes in ionic and hemodynamic functional structure and neurovascular coupling and provide a framework for subsequent studies of neurological disease such as stroke.

Cortex-wide BOLD fMRI activity reflects locally-recorded slow oscillation-associated calcium waves

Spontaneous slow oscillation-associated slow wave activity represents an internally generated state which is characterized by alternations of network quiescence and stereotypical episodes of neuronal activity - slow wave events. However, it remains unclear which macroscopic signal is related to these active periods of the slow wave rhythm. We used optic fiber-based calcium recordings of local neural populations in cortex and thalamus to detect neurophysiologically defined slow calcium waves in isoflurane anesthetized rats. The individual slow wave events were used for an event-related analysis of simultaneously acquired whole-brain BOLD fMRI. We identified BOLD responses directly related to onsets of slow calcium waves, revealing a cortex-wide BOLD correlate: the entire cortex was engaged in this specific type of slow wave activity. These findings demonstrate a direct relation of defined neurophysiological events to a specific BOLD activity pattern and were confirmed for ongoing slow wave activity by independent component and seed-based analyses.

When a person is in a deep non-dreaming sleep, neurons in their brain alternate slowly between periods of silence and periods of activity. This gives rise to low-frequency brain rhythms called slow waves, which are thought to help stabilize memories. Slow wave activity can be detected on multiple scales, from the pattern of electrical impulses sent by an individual neuron to the collective activity of the brain’s entire outer layer, the cortex. But does slow wave activity in an individual group of neurons in the cortex affect the activity of the rest of the brain?

To find out, Schwalm, Schmid, Wachsmuth et al. took advantage of the fact that slow waves also occur under general anesthesia, and placed anesthetized rats inside miniature whole-brain scanners. A small region of cortex in each rat had been injected with a dye that fluoresces whenever the neurons in that region are active. An optical fiber was lowered into the rat’s brain to transmit the fluorescence signals to a computer. Monitoring these signals while the animals lay inside the scanner revealed that slow-wave activity in any one group of cortical neurons was accompanied by slow-wave activity across the cortex as a whole. This relationship was seen only for slow waves, and not for other brain rhythms.

Slow waves seem to occur in all species of animal with a backbone, and in both healthy and diseased brains. While it is not possible to inject fluorescent dyes into the human brain, it is possible to monitor neuronal activity using electrodes. Comparing local electrode recordings with measures of whole-brain activity from scanners could thus allow similar experiments to be performed in people. There is growing evidence – from animal models and from studies of patients – that slow waves may be altered in Alzheimer’s disease. Further work is required to determine whether detecting these changes could help diagnose disease at earlier stages, and whether reversing them may have therapeutic potential.

Phase-synchronization-based parcellation of resting state fMRI signals reveals topographically organized clusters in early visual cortex

Resting-state fMRI is widely used to study brain function and connectivity. However, interpreting patterns of resting state (RS) fMRI activity remains challenging as they may arise from different neuronal mechanisms than those triggered by exogenous events. Currently, this limits the use of RS-fMRI for understanding cortical function in health and disease. Here, we examine the phase synchronization (PS) properties of blood-oxygen level dependent (BOLD) signals obtained during visual field mapping (VFM) and RS with 7T fMRI. This data-driven approach exploits spatiotemporal covariations in the phase of BOLD recordings to establish the presence of clusters of synchronized activity. We find that, in both VFM and RS data, selecting the most synchronized neighboring recording sites identifies spatially localized PS clusters that follow the topographic organization of the visual cortex. However, in activity obtained during VFM, PS is spatially more extensive than in RS activity, likely reflecting stimulus-driven interactions between local responses. Nevertheless, the similarity of the PS clusters obtained for RS and stimulus-driven fMRI suggest that they share a common neuroanatomical origin. Our finding justifies and facilitates direct comparison of RS and stimulus-evoked activity.

Functional Segregation and Development of Mouse Higher Visual Areas

Recent studies suggest that higher visual areas (HVAs) in the mouse visual cortex are segregated anatomically into two visual streams, likely analogous to the ventral and dorsal streams in primates. However, HVAs in mice have yet to be characterized functionally. Moreover, it is unknown when the functional segregation of HVAs occurs during development. Here, we investigated spatiotemporal selectivity of HVAs and their development using wide-field calcium imaging. We found that lateral HVAs in the anatomical ventral stream shared similar spatiotemporal selectivity, whereas the spatiotemporal selectivity of anterior and medial HVAs in the anatomical dorsal stream was not uniform and these areas were segregated functionally into multiple groups. This functional segregation of HVAs developed and reached an adult-like pattern ∼10 d after eye opening (EO). These results suggest, not only the functional segregation of ventral and dorsal streams, but also the presence of multiple substreams in the dorsal stream, and indicate that the functional segregation of visual streams occurs gradually after EO.SIGNIFICANCE STATEMENT Investigation of the spatiotemporal selectivity of nine higher visual areas (HVAs) in adult and developing mice revealed that lateral HVAs belonging to the putative ventral stream shared similar spatiotemporal selectivity, whereas the spatiotemporal selectivity of anterior and medial HVAs belonging to the putative dorsal stream was not uniform and these areas were segregated functionally into multiple groups. These results suggest the presence of multiple substreams within the putative dorsal stream for visuospatial processing. Furthermore, we found that initially immature functional segregation among HVAs developed to an adult-like pattern ∼10 d after eye opening. These results provide a foundation for using mouse HVAs as a model to understand parallel processing and its developmental mechanism.

Low-frequency hippocampal-cortical activity drives brain-wide resting-state functional MRI connectivity

The hippocampus, including the dorsal dentate gyrus (dDG), and cortex engage in bidirectional communication. We propose that low-frequency activity in hippocampal-cortical pathways contributes to brain-wide resting-state connectivity to integrate sensory information. Using optogenetic stimulation and brain-wide fMRI and resting-state fMRI (rsfMRI), we determined the large-scale effects of spatiotemporal-specific downstream propagation of hippocampal activity. Low-frequency (1 Hz), but not high-frequency (40 Hz), stimulation of dDG excitatory neurons evoked robust cortical and subcortical brain-wide fMRI responses. More importantly, it enhanced interhemispheric rsfMRI connectivity in various cortices and hippocampus. Subsequent local field potential recordings revealed an increase in slow oscillations in dorsal hippocampus and visual cortex, interhemispheric visual cortical connectivity, and hippocampal-cortical connectivity. Meanwhile, pharmacological inactivation of dDG neurons decreased interhemispheric rsfMRI connectivity. Functionally, visually evoked fMRI responses in visual regions also increased during and after low-frequency dDG stimulation. Together, our results indicate that low-frequency activity robustly propagates in the dorsal hippocampal-cortical pathway, drives interhemispheric cortical rsfMRI connectivity, and mediates visual processing.

Mesoscale Mapping of Mouse Cortex Reveals Frequency-Dependent Cycling between Distinct Macroscale Functional Modules

Connectivity mapping based on resting-state activity in mice has revealed functional motifs of correlated activity. However, the rules by which motifs organize into larger functional modules that lead to hemisphere wide spatial-temporal activity sequences is not clear. We explore cortical activity parcellation in head-fixed, quiet awake GCaMP6 mice from both sexes by using mesoscopic calcium imaging. Spectral decomposition of spontaneous cortical activity revealed the presence of two dominant frequency modes (<1 and ∼3 Hz), each of them associated with a unique spatial signature of cortical macro-parcellation not predicted by classical cytoarchitectonic definitions of cortical areas. Based on assessment of 0.1-1 Hz activity, we define two macro-organizing principles: the first being a rotating polymodal-association pinwheel structure around which activity flows sequentially from visual to barrel then to hindlimb somatosensory; the second principle is correlated activity symmetry planes that exist on many levels within a single domain such as intrahemispheric reflections of sensory and motor cortices. In contrast, higher frequency activity >1 Hz yielded two larger clusters of coactivated areas with an enlarged default mode network-like posterior region. We suggest that the apparent constrained structure for intra-areal cortical activity flow could be exploited in future efforts to normalize activity in diseases of the nervous system.SIGNIFICANCE STATEMENT Increasingly, functional connectivity mapping of spontaneous activity is being used to reveal the organization of the brain. However, because the brain operates across multiple space and time domains a more detailed understanding of this organization is necessary. We used in vivo wide-field calcium imaging of the indicator GCaMP6 in head-fixed, awake mice to characterize the organization of spontaneous cortical activity at different spatiotemporal scales. Correlation analysis defines the presence of two to three superclusters of activity that span traditionally defined functional territories and were frequency dependent. This work helps define the rules for how different cortical areas interact in time and space. We provide a framework necessary for future studies that explore functional reorganization of brain circuits in disease models.

A Comprehensive Analysis of the Correlations between Resting-State Oscillations in Multiple-Frequency Bands and Big Five Traits

Recently, the association between human personality traits and resting-state brain activity has gained interest in neuroimaging studies. However, it remains unclear if Big Five personality traits are represented in frequency bands (~0.25 Hz) of resting-state functional magnetic resonance imaging (fMRI) activity. Based on earlier neurophysiological studies, we investigated the correlation between the five personality traits assessed by the NEO Five-Factor Inventory (NEO-FFI), and the fractional amplitude of low-frequency fluctuation (fALFF) at four distinct frequency bands (slow-5 (0.01–0.027 Hz), slow-4 (0.027–0.073 Hz), slow-3 (0.073–0.198 Hz) and slow-2 (0.198–0.25 Hz)). We enrolled 835 young subjects and calculated the correlations of resting-state fMRI signals using a multiple regression analysis. We found a significant and consistent correlation between fALFF and the personality trait of extraversion at all frequency bands. Furthermore, significant correlations were detected in distinct brain regions for each frequency band. This finding supports the frequency-specific spatial representations of personality traits as previously suggested. In conclusion, our data highlight an association between human personality traits and fALFF at four distinct frequency bands.

Reconsidering Tonotopic Maps in the Auditory Cortex and Lemniscal Auditory Thalamus in Mice

The auditory thalamus and auditory cortex (AC) are pivotal structures in the central auditory system. However, the thalamocortical mechanisms of processing sounds are largely unknown. Investigation of this process benefits greatly from the use of mice because the mouse is a powerful animal model in which various experimental techniques, especially genetic tools, can be applied. However, the use of mice has been limited in auditory research, and thus even basic anatomical knowledge of the mouse central auditory system has not been sufficiently collected. Recently, optical imaging combined with morphological analyses has enabled the elucidation of detailed anatomical properties of the mouse auditory system. These techniques have uncovered fine AC maps with multiple frequency-organized regions, each of which receives point-to-point thalamocortical projections from different origins inside the lemniscal auditory thalamus, the ventral division of the medial geniculate body (MGv). This precise anatomy now provides a platform for physiological research. In this mini review article, we summarize these recent achievements that will facilitate physiological investigations in the mouse auditory system.

The global signal in fMRI: Nuisance or Information?

The global signal is widely used as a regressor or normalization factor for removing the effects of global variations in the analysis of functional magnetic resonance imaging (fMRI) studies. However, there is considerable controversy over its use because of the potential bias that can be introduced when it is applied to the analysis of both task-related and resting-state fMRI studies. In this paper we take a closer look at the global signal, examining in detail the various sources that can contribute to the signal. For the most part, the global signal has been treated as a nuisance term, but there is growing evidence that it may also contain valuable information. We also examine the various ways that the global signal has been used in the analysis of fMRI data, including global signal regression, global signal subtraction, and global signal normalization. Furthermore, we describe new ways for understanding the effects of global signal regression and its relation to the other approaches.

Brain Rhythms of Pain

Pain is an integrative phenomenon that results from dynamic interactions between sensory and contextual (i.e., cognitive, emotional, and motivational) processes. In the brain the experience of pain is associated with neuronal oscillations and synchrony at different frequencies. However, an overarching framework for the significance of oscillations for pain remains lacking. Recent concepts relate oscillations at different frequencies to the routing of information flow in the brain and the signaling of predictions and prediction errors. The application of these concepts to pain promises insights into how flexible routing of information flow coordinates diverse processes that merge into the experience of pain. Such insights might have implications for the understanding and treatment of chronic pain.

Combined spatiotemporal ICA (stICA) for continuous and dynamic lag structure analysis of MREG data

This study investigated lag structure in the resting-state fMRI by applying a novel independent component (ICA) method to magnetic resonance encephalography (MREG) data. Briefly, the spatial ICA (sICA) was used for defining the frontal and back nodes of the default mode network (DMN), and the temporal ICA (tICA), which is enabled by the high temporal resolution of MREG (TR=100ms), was used to separate both neuronal and physiological components of these two spatial map regions. Subsequently, lag structure was investigated between the frontal (DMNvmpf) and posterior (DMNpcc) DMN nodes using both conventional method with all-time points and a sliding-window approach. A rigorous noise exclusion criterion was applied for tICs to remove physiological pulsations, motion and system artefacts. All the de-noised tICs were used to calculate the null-distributions both for expected lag variability over time and over subjects. Lag analysis was done for the three highest correlating denoised tICA pairs. Mean time lag of 0.6s (± 0.5 std) and mean absolute correlation of 0.69 (± 0.08) between the highest correlating tICA pairs of DMN nodes was observed throughout the whole analyzed period. In dynamic 2min window analysis, there was large variability over subjects as ranging between 1-10sec. Directionality varied between these highly correlating sources an average 28.8% of the possible number of direction changes. The null models show highly consistent correlation and lag structure between DMN nodes both in continuous and dynamic analysis. The mean time lag of a null-model over time between all denoised DMN nodes was 0.0s and, thus the probability of having either DMNpcc or DMNvmpf as a preceding component is near equal. All the lag values of highest correlating tICA pairs over subjects lie within the standard deviation range of a null-model in whole time window analysis, supporting the earlier findings that there is a consistent temporal lag structure across groups of individuals. However, in dynamic analysis, there are lag values exceeding the threshold of significance of a null-model meaning that there might be biologically meaningful variation in this measure. Taken together the variability in lag and the presence of high activity peaks during strong connectivity indicate that individual avalanches may play an important role in defining dynamic independence in resting state connectivity within networks.

Degradation in intrinsic connectivity networks across the Alzheimer's disease spectrum

Introduction

Changes in intrinsic functional connectivity (iFC) have been reported at various stages of the Alzheimer's disease (AD) spectrum. We aimed to investigate such alterations over a variety of large-scale intrinsic brain networks (iBNs) across the spectrum of amyloid β positivity and uncover their relation to cognitive impairment.

Methods

Eight iBNs were defined from resting-state functional magnetic resonance imaging data. In amyloid β–positive healthy subjects, prodromal, and AD patients (N = 70), within-network iFC (intra-iFC) and between-network iFC (inter-iFC) were correlated with scores of cognitive impairment.

Results

Across all iBNs, a general degradation in intra-iFC along the scale of cognitive impairment severity was found. Only subtle changes in inter-iFC were identified.

Discussion

Across the AD spectrum, changes in iFC that are strongly related to cognitive impairment occur within an extensive variety of networks.

Sources and implications of whole-brain fMRI signals in humans

Whole-brain fMRI signals are a subject of intense interest: variance in the global fMRI signal (the spatial mean of all signals in the brain) indexes subject arousal, and psychiatric conditions such as schizophrenia and autism have been characterized by differences in the global fMRI signal. Further, vigorous debates exist on whether global signals ought to be removed from fMRI data. However, surprisingly little research has focused on the empirical properties of whole-brain fMRI signals. Here we map the spatial and temporal properties of the global signal, individually, in 1000+ fMRI scans. Variance in the global fMRI signal is strongly linked to head motion, to hardware artifacts, and to respiratory patterns and their attendant physiologic changes. Many techniques used to prepare fMRI data for analysis fail to remove these uninteresting kinds of global signal fluctuations. Thus, many studies include, at the time of analysis, prominent global effects of yawns, breathing changes, and head motion, among other signals. Such artifacts will mimic dynamic neural activity and will spuriously alter signal covariance throughout the brain. Methods capable of isolating and removing global artifactual variance while preserving putative "neural" variance are needed; this paper adopts no position on the topic of global signal regression.