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

Evolution of aberrant brain-wide spatiotemporal dynamics of resting-state networks in a Huntington's disease mouse model

BACKGROUND

Huntington's disease (HD) is marked by irreversible loss of neuronal function for which currently no availability for disease-modifying treatment exists. Advances in the understanding of disease progression can aid biomarker development, which in turn can accelerate therapeutic discovery.

METHODS

We characterised the progression of altered dynamics of whole-brain network states in the zQ175DN mouse model of HD using a dynamic functional connectivity (FC) approach to resting-state fMRI and identified quasi-periodic patterns (QPPs) of brain activity constituting the most prominent resting-state networks.

RESULTS

The occurrence of the normative QPPs, as observed in healthy controls, was reduced in the HD model as the phenotype progressed. This uncovered progressive cessation of synchronous brain activity with phenotypic progression, which is not observed with the conventional static FC approaches. To better understand the potential underlying cause of the observed changes in these brain states, we further assessed how mutant huntingtin (mHTT) protein deposition affects astrocytes and pericytes - one of the most important effectors of neurovascular coupling, along phenotypic progression. Increased cell-type dependent mHTT deposition was observed at the age of onset of motor anomalies, in the caudate putamen, somatosensory and motor cortex, regions that are prominently involved in HD pathology as seen in humans.

CONCLUSION

Our findings provide meaningful insights into the development and progression of altered functional brain dynamics in this HD model and open new avenues in assessing the dynamics of whole brain states, through QPPs, in clinical HD research.

HIGHLIGHTS

Hyperactivity in the LCN-linked regions within short QPPs observed before motor impairment onset. DMLN QPP presents a progressive decrease in DMLN activity and occurrence along HD-like phenotype development. Breakdown of the LCN DMLN state flux at motor onset leads to a subsequent absence of the LCN DMLN QPP at an advanced HD-like stage.

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.

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.

Early altered directionality of resting brain network state transitions in the TgF344-AD rat model of Alzheimer's disease

Introduction

Alzheimer's disease (AD) is a progressive neurodegenerative disease resulting in memory loss and cognitive decline. Synaptic dysfunction is an early hallmark of the disease whose effects on whole-brain functional architecture can be identified using resting-state functional MRI (rsfMRI). Insights into mechanisms of early, whole-brain network alterations can help our understanding of the functional impact of AD's pathophysiology.

Methods

Here, we obtained rsfMRI data in the TgF344-AD rat model at the pre- and early-plaque stages. This model recapitulates the major pathological and behavioral hallmarks of AD. We used co-activation pattern (CAP) analysis to investigate if and how the dynamic organization of intrinsic brain functional networks states, undetectable by earlier methods, is altered at these early stages.

Results

We identified and characterized six intrinsic brain states as CAPs, their spatial and temporal features, and the transitions between the different states. At the pre-plaque stage, the TgF344-AD rats showed reduced co-activation of hub regions in the CAPs corresponding to the default mode-like and lateral cortical network. Default mode-like network activity segregated into two distinct brain states, with one state characterized by high co-activation of the basal forebrain. This basal forebrain co-activation was reduced in TgF344-AD animals mainly at the pre-plaque stage. Brain state transition probabilities were altered at the pre-plaque stage between states involving the default mode-like network, lateral cortical network, and basal forebrain regions. Additionally, while the directionality preference in the network-state transitions observed in the wild-type animals at the pre-plaque stage had diminished at the early-plaque stage, TgF344-AD animals continued to show directionality preference at both stages.

Discussion

Our study enhances the understanding of intrinsic brain state dynamics and how they are impacted at the early stages of AD, providing a nuanced characterization of the early, functional impact of the disease's neurodegenerative process.

Global spatiotemporal synchronizing structures of spontaneous neural activities in different cell types

Increasing evidence has revealed the large-scale nonstationary synchronizations as traveling waves in spontaneous neural activity. However, the interplay of various cell types in fine-tuning these spatiotemporal patters remains unclear. Here, we performed comprehensive exploration of spatiotemporal synchronizing structures across different cell types, states (awake, anesthesia, motion) and developmental axis in male mice. We found traveling waves in glutamatergic neurons exhibited greater variety than those in GABAergic neurons. Moreover, the synchronizing structures of GABAergic neurons converged toward those of glutamatergic neurons during development, but the evolution of waves exhibited varying timelines for different sub-type interneurons. Functional connectivity arises from both standing and traveling waves, and negative connections can be elucidated by the spatial propagation of waves. In addition, some traveling waves were correlated with the spatial distribution of gene expression. Our findings offer further insights into the neural underpinnings of traveling waves, functional connectivity, and resting-state networks, with cell-type specificity and developmental perspectives.

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.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Abnormal Spatial and Temporal Overlap of Time-Varying Brain Functional Networks in Patients with Schizophrenia

Schizophrenia (SZ) is a complex psychiatric disorder with unclear etiology and pathological features. Neuroscientists are increasingly proposing that schizophrenia is an abnormality in the dynamic organization of brain networks. Previous studies have found that the dynamic brain networks of people with SZ are abnormal in both space and time. However, little is known about the interactions and overlaps between hubs of the brain underlying spatiotemporal dynamics. In this study, we aimed to investigate different patterns of spatial and temporal overlap of hubs between SZ patients and healthy individuals. Specifically, we obtained resting-state functional magnetic resonance imaging data from the public dataset for 43 SZ patients and 49 healthy individuals. We derived a representation of time-varying functional connectivity using the Jackknife Correlation (JC) method. We employed the Betweenness Centrality (BC) method to identify the hubs of the brain's functional connectivity network. We then applied measures of temporal overlap, spatial overlap, and hierarchical clustering to investigate differences in the organization of brain hubs between SZ patients and healthy controls. Our findings suggest significant differences between SZ patients and healthy controls at the whole-brain and subnetwork levels. Furthermore, spatial overlap and hierarchical clustering analysis showed that quasi-periodic patterns were disrupted in SZ patients. Analyses of temporal overlap revealed abnormal pairwise engagement preferences in the hubs of SZ patients. These results provide new insights into the dynamic characteristics of the network organization of the SZ brain.

Mapping and comparing fMRI connectivity networks across species

Technical advances in neuroimaging, notably in fMRI, have allowed distributed patterns of functional connectivity to be mapped in the human brain with increasing spatiotemporal resolution. Recent years have seen a growing interest in extending this approach to rodents and non-human primates to understand the mechanism of fMRI connectivity and complement human investigations of the functional connectome. Here, we discuss current challenges and opportunities of fMRI connectivity mapping across species. We underscore the critical importance of physiologically decoding neuroimaging measures of brain (dys)connectivity via multiscale mechanistic investigations in animals. We next highlight a set of general principles governing the organization of mammalian connectivity networks across species. These include the presence of evolutionarily conserved network systems, a dominant cortical axis of functional connectivity, and a common repertoire of topographically conserved fMRI spatiotemporal modes. We finally describe emerging approaches allowing comparisons and extrapolations of fMRI connectivity findings across species. As neuroscientists gain access to increasingly sophisticated perturbational, computational and recording tools, cross-species fMRI offers novel opportunities to investigate the large-scale organization of the mammalian brain in health and disease.

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.

Longitudinal investigation of changes in resting-state co-activation patterns and their predictive ability in the zQ175 DN mouse model of Huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder caused by expanded (≥ 40) glutamine-encoding CAG repeats in the huntingtin gene, which leads to dysfunction and death of predominantly striatal and cortical neurons. While the genetic profile and clinical signs and symptoms of the disease are better known, changes in the functional architecture of the brain, especially before the clinical expression becomes apparent, are not fully and consistently characterized. In this study, we sought to uncover functional changes in the brain in the heterozygous (HET) zQ175 delta-neo (DN) mouse model at 3, 6, and 10 months of age, using resting-state functional magnetic resonance imaging (RS-fMRI). This mouse model shows molecular, cellular and circuitry alterations that worsen through age. Motor function disturbances are manifested in this model at 6 and 10 months of age. Specifically, we investigated, longitudinally, changes in co-activation patterns (CAPs) that are the transient states of brain activity constituting the resting-state networks (RSNs). Most robust changes in the temporal properties of CAPs occurred at the 10-months time point; the durations of two anti-correlated CAPs, characterized by simultaneous co-activation of default-mode like network (DMLN) and co-deactivation of lateral-cortical network (LCN) and vice-versa, were reduced in the zQ175 DN HET animals compared to the wild-type mice. Changes in the spatial properties, measured in terms of activation levels of different brain regions, during CAPs were found at all three ages and became progressively more pronounced at 6-, and 10 months of age. We then assessed the cross-validated predictive power of CAP metrics to distinguish HET animals from controls. Spatial properties of CAPs performed significantly better than the chance level at all three ages with 80% classification accuracy at 6 and 10 months of age.

Influence of vigilance-related arousal on brain dynamics: Potentials of new approaches

Growing research has focused on how mesoscopic activity in the brain develops over time and space. Recent influential studies using functional imaging have characterized brain dynamics in terms of the spread of activation across the brain following a unimodal to transmodal axis. In parallel, a number of studies have assessed changes of brain connectivity in terms of vigilance-linked arousal. Here I offer a view on how these two lines of research can lead to a deeper understanding of how arousal shapes the brain's dynamic behavior. This knowledge could have great impact on the investigation of mental disease.

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.

Functional network antagonism and consciousness

Spontaneous brain activity changes across states of consciousness. A particular consciousness-mediated configuration is the anticorrelations between the default mode network and other brain regions. What this antagonistic organization implies about consciousness to date remains inconclusive. In this Perspective Article, we propose that anticorrelations are the physiological expression of the concept of segregation, namely the brain's capacity to show selectivity in the way areas will be functionally connected. We postulate that this effect is mediated by the process of neural inhibition, by regulating global and local inhibitory activity. While recognizing that this effect can also result from other mechanisms, neural inhibition helps the understanding of how network metastability is affected after disrupting local and global neural balance. In combination with relevant theories of consciousness, we suggest that anticorrelations are a physiological prior that can work as a marker of preserved consciousness. We predict that if the brain is not in a state to host anticorrelations, then most likely the individual does not entertain subjective experience. We believe that this link between anticorrelations and the underlying physiology will help not only to comprehend how consciousness happens, but also conceptualize effective interventions for treating consciousness disorders in which anticorrelations seem particularly affected.

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.

Age-related change in task-evoked amygdala-prefrontal circuitry: A multiverse approach with an accelerated longitudinal cohort aged 4-22 years

The amygdala and its connections with medial prefrontal cortex (mPFC) play central roles in the development of emotional processes. While several studies have suggested that this circuitry exhibits functional changes across the first two decades of life, findings have been mixed - perhaps resulting from differences in analytic choices across studies. Here we used multiverse analyses to examine the robustness of task-based amygdala-mPFC function findings to analytic choices within the context of an accelerated longitudinal design (4-22 years-old; N = 98; 183 scans; 1-3 scans/participant). Participants recruited from the greater Los Angeles area completed an event-related emotional face (fear, neutral) task. Parallel analyses varying in preprocessing and modeling choices found that age-related change estimates for amygdala reactivity were more robust than task-evoked amygdala-mPFC functional connectivity to varied analytical choices. Specification curves indicated evidence for age-related decreases in amygdala reactivity to faces, though within-participant changes in amygdala reactivity could not be differentiated from between-participant differences. In contrast, amygdala-mPFC functional connectivity results varied across methods much more, and evidence for age-related change in amygdala-mPFC connectivity was not consistent. Generalized psychophysiological interaction (gPPI) measurements of connectivity were especially sensitive to whether a deconvolution step was applied. Our findings demonstrate the importance of assessing the robustness of findings to analysis choices, although the age-related changes in our current work cannot be overinterpreted given low test-retest reliability. Together, these findings highlight both the challenges in estimating developmental change in longitudinal cohorts and the value of multiverse approaches in developmental neuroimaging for assessing robustness of results.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

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.

Monitoring Neuronal Network Disturbances of Brain Diseases: A Preclinical MRI Approach in the Rodent Brain

Functional and structural neuronal networks, as recorded by resting-state functional MRI and diffusion MRI-based tractography, gain increasing attention as data driven whole brain imaging methods not limited to the foci of the primary pathology or the known key affected regions but permitting to characterize the entire network response of the brain after disease or injury. Their connectome contents thus provide information on distal brain areas, directly or indirectly affected by and interacting with the primary pathological event or affected regions. From such information, a better understanding of the dynamics of disease progression is expected. Furthermore, observation of the brain's spontaneous or treatment-induced improvement will contribute to unravel the underlying mechanisms of plasticity and recovery across the whole-brain networks. In the present review, we discuss the values of functional and structural network information derived from systematic and controlled experimentation using clinically relevant animal models. We focus on rodent models of the cerebral diseases with high impact on social burdens, namely, neurodegeneration, and stroke.

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.

An isotropic EPI database and analytical pipelines for rat brain resting-state fMRI

Resting-state functional magnetic resonance imaging (fMRI) has drastically expanded the scope of brain research by advancing our knowledge about the topologies, dynamics, and interspecies translatability of functional brain networks. Several databases have been developed and shared in accordance with recent key initiatives in the rodent fMRI community to enhance the transparency, reproducibility, and interpretability of data acquired at various sites. Despite these pioneering efforts, one notable challenge preventing efficient standardization in the field is the customary choice of anisotropic echo planar imaging (EPI) schemes with limited spatial coverage. Imaging with anisotropic resolution and/or reduced brain coverage has significant shortcomings including reduced registration accuracy and increased deviation in brain feature detection. Here we proposed a high-spatial-resolution (0.4 mm), isotropic, whole-brain EPI protocol for the rat brain using a horizontal slicing scheme that can maintain a functionally relevant repetition time (TR), avoid high gradient duty cycles, and offer unequivocal whole-brain coverage. Using this protocol, we acquired resting-state EPI fMRI data from 87 healthy rats under the widely used dexmedetomidine sedation supplemented with low-dose isoflurane on a 9.4 T MRI system. We developed an EPI template that closely approximates the Paxinos and Watson's rat brain coordinate system and demonstrated its ability to improve the accuracy of group-level approaches and streamline fMRI data pre-processing. Using this database, we employed a multi-scale dictionary-learning approach to identify reliable spatiotemporal features representing rat brain intrinsic activity. Subsequently, we performed k-means clustering on those features to obtain spatially discrete, functional regions of interest (ROIs). Using Euclidean-based hierarchical clustering and modularity-based partitioning, we identified the topological organizations of the rat brain. Additionally, the identified group-level FC network appeared robust across strains and sexes. The "triple-network" commonly adapted in human fMRI were resembled in the rat brain. Through this work, we disseminate raw and pre-processed isotropic EPI data, a rat brain EPI template, as well as identified functional ROIs and networks in standardized rat brain coordinates. We also make our analytical pipelines and scripts publicly available, with the hope of facilitating rat brain resting-state fMRI study standardization.

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.

Emerging imaging methods to study whole-brain function in rodent models

In the past decade, the idea that single populations of neurons support cognition and behavior has gradually given way to the realization that connectivity matters and that complex behavior results from interactions between remote yet anatomically connected areas that form specialized networks. In parallel, innovation in brain imaging techniques has led to the availability of a broad set of imaging tools to characterize the functional organization of complex networks. However, each of these tools poses significant technical challenges and faces limitations, which require careful consideration of their underlying anatomical, physiological, and physical specificity. In this review, we focus on emerging methods for measuring spontaneous or evoked activity in the brain. We discuss methods that can measure large-scale brain activity (directly or indirectly) with a relatively high temporal resolution, from milliseconds to seconds. We further focus on methods designed for studying the mammalian brain in preclinical models, specifically in mice and rats. This field has seen a great deal of innovation in recent years, facilitated by concomitant innovation in gene-editing techniques and the possibility of more invasive recordings. This review aims to give an overview of currently available preclinical imaging methods and an outlook on future developments. This information is suitable for educational purposes and for assisting scientists in choosing the appropriate method for their own research question.

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.

Characteristics of fMRI responses to visual stimulation in anesthetized vs. awake mice

The functional characteristics of the mouse visual system have not previously been well explored using fMRI. In this research, we examined 9.4 T BOLD fMRI responses to visual stimuli of varying pulse durations (1 - 50 ms) and temporal frequencies (1 - 10 Hz) under ketamine and xylazine anesthesia, and compared fMRI responses of anesthetized and awake mice. Under anesthesia, significant positive BOLD responses were detected bilaterally in the major structures of the visual pathways, including the dorsal lateral geniculate nuclei, superior colliculus, lateral posterior nucleus of thalamus, primary visual area, and higher-order visual area. BOLD responses increased slightly with pulse duration, were maximal at 3 - 5 Hz stimulation, and significantly decreased at 10 Hz, which were all consistent with previous neurophysiological findings. When the mice were awake, the BOLD fMRI response was faster in all active regions and stronger in the subcortical areas compared with the anesthesia condition. In the V1, the BOLD response was biphasic for 5 Hz stimulation and negative for 10 Hz stimulation under wakefulness, whereas prolonged positive BOLD responses were observed at both frequencies under anesthesia. Unexpected activation was detected in the extrastriate postrhinal area and non-visual subiculum complex under anesthesia, but not under wakefulness. Widespread positive BOLD activity under anesthesia likely results from the disinhibition and sensitization of excitatory neurons induced by ketamine. Overall, fMRI can be a viable tool for mapping brain-wide functional networks.

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.

Prevalent and sex-biased breathing patterns modify functional connectivity MRI in young adults

Resting state functional connectivity magnetic resonance imaging (fMRI) is a tool for investigating human brain organization. Here we identify, visually and algorithmically, two prevalent influences on fMRI signals during 440 h of resting state scans in 440 healthy young adults, both caused by deviations from normal breathing which we term deep breaths and bursts. The two respiratory patterns have distinct influences on fMRI signals and signal covariance, distinct timescales, distinct cardiovascular correlates, and distinct tendencies to manifest by sex. Deep breaths are not sex-biased. Bursts, which are serial taperings of respiratory depth typically spanning minutes at a time, are more common in males. Bursts share features of chemoreflex-driven clinical breathing patterns that also occur primarily in males, with notable neurological, psychiatric, medical, and lifespan associations. These results identify common breathing patterns in healthy young adults with distinct influences on functional connectivity and an ability to differentially influence resting state fMRI studies.

Functional connectivity measured from fMRI data is widely used in neuroscience. Here the authors report an association between two types of breathing signature and obtained BOLD data, and associated sex differences.

Cholinergic Modulation of the Default Mode Like Network in Rats

The discovery of the default mode network (DMN), a large-scale brain network that is suppressed during attention-demanding tasks, had major impact in neuroscience. This network exhibits an antagonistic relationship with attention-related networks. A better understanding of the processes underlying modulation of DMN is imperative, as this network is compromised in several neurological diseases. Cholinergic neuromodulation is one of the major regulatory networks for attention, and studies suggest a role in regulation of the DMN. In this study, we unilaterally activated the right basal forebrain cholinergic neurons and observed decreased right intra-hemispheric and interhemispheric FC in the default mode like network (DMLN). Our findings provide critical insights into the interplay between cholinergic neuromodulation and DMLN, demonstrate that differential effects can be exerted between the two hemispheres by unilateral stimulation, and open windows for further studies involving directed modulations of DMN in treatments for diseases demonstrating compromised DMN activity.

Bring the Noise: Reconceptualizing Spontaneous Neural Activity

Definitions of what constitutes the 'signal of interest' in neuroscience can be controversial, due in part to continuously evolving notions regarding the significance of spontaneous neural activity. This review highlights how the challenge of separating brain signal from noise has led to new conceptualizations of brain functional organization at both the micro- and macroscopic level. Recent debates in the functional neuroimaging community surrounding artifact removal processes have revived earlier discussions surrounding how to appropriately isolate and measure neuronal signals against a background of noise from various sources. Insights from electrophysiological studies and computational modeling can inform current theory and data analytic practices in human functional neuroimaging, given that signal and noise may be inextricably linked in the brain.

Brain states govern the spatio-temporal dynamics of resting-state functional connectivity

Previously, using simultaneous resting-state functional magnetic resonance imaging (fMRI) and photometry-based neuronal calcium recordings in the anesthetized rat, we identified blood oxygenation level-dependent (BOLD) responses directly related to slow calcium waves, revealing a cortex-wide and spatially organized correlate of locally recorded neuronal activity (Schwalm et al., 2017). Here, using the same techniques, we investigate two distinct cortical activity states: persistent activity, in which compartmentalized network dynamics were observed; and slow wave activity, dominated by a cortex-wide BOLD component, suggesting a strong functional coupling of inter-cortical activity. During slow wave activity, we find a correlation between the occurring slow wave events and the strength of functional connectivity between different cortical areas. These findings suggest that down-up transitions of neuronal excitability can drive cortex-wide functional connectivity. This study provides further evidence that changes in functional connectivity are dependent on the brain's current state, directly linked to the generation of slow waves.

The (Un)Conscious Mouse as a Model for Human Brain Functions: Key Principles of Anesthesia and Their Impact on Translational Neuroimaging

In recent years, technical and procedural advances have brought functional magnetic resonance imaging (fMRI) to the field of murine neuroscience. Due to its unique capacity to measure functional activity non-invasively, across the entire brain, fMRI allows for the direct comparison of large-scale murine and human brain functions. This opens an avenue for bidirectional translational strategies to address fundamental questions ranging from neurological disorders to the nature of consciousness. The key challenges of murine fMRI are: (1) to generate and maintain functional brain states that approximate those of calm and relaxed human volunteers, while (2) preserving neurovascular coupling and physiological baseline conditions. Low-dose anesthetic protocols are commonly applied in murine functional brain studies to prevent stress and facilitate a calm and relaxed condition among animals. Yet, current mono-anesthesia has been shown to impair neural transmission and hemodynamic integrity. By linking the current state of murine electrophysiology, Ca2+ imaging and fMRI of anesthetic effects to findings from human studies, this systematic review proposes general principles to design, apply and monitor anesthetic protocols in a more sophisticated way. The further development of balanced multimodal anesthesia, combining two or more drugs with complementary modes of action helps to shape and maintain specific brain states and relevant aspects of murine physiology. Functional connectivity and its dynamic repertoire as assessed by fMRI can be used to make inferences about cortical states and provide additional information about whole-brain functional dynamics. Based on this, a simple and comprehensive functional neurosignature pattern can be determined for use in defining brain states and anesthetic depth in rest and in response to stimuli. Such a signature can be evaluated and shared between labs to indicate the brain state of a mouse during experiments, an important step toward translating findings across species.

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.

fMRI study of olfactory processing in mice under three anesthesia protocols: Insight into the effect of ketamine on olfactory processing

Functional magnetic resonance imaging (fMRI) is a valuable tool for studying neural activations in the central nervous system of animals due to its wide spatial coverage and non-invasive nature. However, the advantages of fMRI have not been fully realized in functional studies in mice, especially in the olfactory system, possibly due to the lack of suitable anesthesia protocols with spontaneous breathing. Since mice are widely used in biomedical research, it is desirable to evaluate different anesthesia protocols for olfactory fMRI studies in mice. Dexmedetomidine (DEX) as a sedative/anesthetic has been introduced to fMRI studies in mice, but it has a limited anesthesia duration. To extend the anesthesia duration, DEX has been combined with a low dose of isoflurane (ISO) or ketamine (KET) in previous functional studies in mice. In this report, olfactory fMRI studies were performed under three anesthesia protocols (DEX alone, DEX/ISO, and DEX/KET) in three different groups of mice. Isoamyl-acetate was used as an odorant, and the odorant-induced neural activations were measured by blood oxygenation-level dependent (BOLD) fMRI. BOLD fMRI responses were observed in the olfactory bulb (OB), anterior olfactory nuclei (AON), and piriform cortex (Pir). Interestingly, BOLD fMRI activations were also observed in the prefrontal cortical region (PFC), which are most likely caused by the draining vein effect. The response in the OB showed no adaptation to either repeated odor stimulations or continuous odor exposure, but the response in the Pir showed adaptation during the continuous odor exposure. The data also shows that ISO suppresses the olfactory response in the OB and AON, while KET enhances the olfactory response in the Pir. Thus, DEX/KET should be an attractive anesthesia for olfactory fMRI in mice.

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.

Common functional networks in the mouse brain revealed by multi-centre resting-state fMRI analysis

Preclinical applications of resting-state functional magnetic resonance imaging (rsfMRI) offer the possibility to non-invasively probe whole-brain network dynamics and to investigate the determinants of altered network signatures observed in human studies. Mouse rsfMRI has been increasingly adopted by numerous laboratories worldwide. Here we describe a multi-centre comparison of 17 mouse rsfMRI datasets via a common image processing and analysis pipeline. Despite prominent cross-laboratory differences in equipment and imaging procedures, we report the reproducible identification of several large-scale resting-state networks (RSN), including a mouse default-mode network, in the majority of datasets. A combination of factors was associated with enhanced reproducibility in functional connectivity parameter estimation, including animal handling procedures and equipment performance. RSN spatial specificity was enhanced in datasets acquired at higher field strength, with cryoprobes, in ventilated animals, and under medetomidine-isoflurane combination sedation. Our work describes a set of representative RSNs in the mouse brain and highlights key experimental parameters that can critically guide the design and analysis of future rodent rsfMRI investigations.

Propagations of spontaneous brain activity in awake rats

Slow propagations of spontaneous brain activity have been reported in multiple species. However, systematical investigation of the organization of such brain activity is still lacking. In this study, we analyzed propagations of spontaneous brain activity using a reference library of characteristic resting-state functional connectivity (RSFC) patterns in awake rodents. We found that transitions through multiple distinct RSFC patterns were reproducible not only in transition sequences but also in transition time delays. In addition, the organization of these transitions and their spatiotemporal dynamic patterns were revealed using a graphical model. We further identified prominent brain regions involved in these transitions. These results provide a comprehensive framework of brainwide propagations of spontaneous activity in awake rats. This study also offers a new tool to study the spatiotemporal dynamics of activity in the resting brain.

Quasi-periodic patterns contribute to functional connectivity in the brain

Functional connectivity is widely used to study the coordination of activity between brain regions over time. Functional connectivity in the default mode and task positive networks is particularly important for normal brain function. However, the processes that give rise to functional connectivity in the brain are not fully understood. It has been postulated that low-frequency neural activity plays a key role in establishing the functional architecture of the brain. Quasi-periodic patterns (QPPs) are a reliably observable form of low-frequency neural activity that involve the default mode and task positive networks. Here, QPPs from resting-state and working memory task-performing individuals were acquired. The spatiotemporal pattern, strength, and frequency of the QPPs between the two groups were compared and the contribution of QPPs to functional connectivity in the brain was measured. In task-performing individuals, the spatiotemporal pattern of the QPP changes, particularly in task-relevant regions, and the QPP tends to occur with greater strength and frequency. Differences in the QPPs between the two groups could partially account for the variance in functional connectivity between resting-state and task-performing individuals. The QPPs contribute strongly to connectivity in the default mode and task positive networks and to the strength of anti-correlation seen between the two networks. Many of the connections affected by QPPs are also disrupted during several neurological disorders. These findings contribute to understanding the dynamic neural processes that give rise to functional connectivity in the brain and how they may be disrupted during disease.

Dynamic properties of simulated brain network models and empirical resting-state data

Brain network models (BNMs) have become a promising theoretical framework for simulating signals that are representative of whole-brain activity such as resting-state fMRI. However, it has been difficult to compare the complex brain activity obtained from simulations to empirical data. Previous studies have used simple metrics to characterize coordination between regions such as functional connectivity. We extend this by applying various different dynamic analysis tools that are currently used to understand empirical resting-state fMRI (rs-fMRI) to the simulated data. We show that certain properties correspond to the structural connectivity input that is shared between the models, and certain dynamic properties relate more to the mathematical description of the brain network model. We conclude that the dynamic properties that explicitly examine patterns of signal as a function of time rather than spatial coordination between different brain regions in the rs-fMRI signal seem to provide the largest contrasts between different BNMs and the unknown empirical dynamical system. Our results will be useful in constraining and developing more realistic simulations of whole-brain activity.

Functional Brain Connectivity Revealed by Sparse Coding of Large-Scale Local Field Potential Dynamics

Exploration of brain dynamics patterns has attracted increasing attention due to its fundamental significance in understanding the working mechanism of the brain. However, due to the lack of effective modeling methods, how the simultaneously recorded LFP can inform us about the brain dynamics remains a general challenge. In this paper, we propose a novel sparse coding based method to investigate brain dynamics of freely-behaving mice from the perspective of functional connectivity, using super-long local field potential (LFP) recordings from 13 distinct regions of the mouse brain. Compared with surrogate datasets, six and four reproducible common functional connectivities were discovered to represent the space of brain dynamics in the frequency bands of alpha and theta respectively. Modeled by a finite state machine, temporal transition framework of functional connectivities was inferred for each frequency band, and evident preference was discovered. Our results offer a novel perspective for analyzing neural recording data at such high temporal resolution and recording length, as common functional connectivities and their transition framework discovered in this work reveal the nature of the brain dynamics in freely behaving mice.

Quasi-Periodic Patterns of Neural Activity improve Classification of Alzheimer's Disease in Mice

Resting state (rs)fMRI allows measurement of brain functional connectivity and has identified default mode (DMN) and task positive (TPN) network disruptions as promising biomarkers for Alzheimer's disease (AD). Quasi-periodic patterns (QPPs) of neural activity describe recurring spatiotemporal patterns that display DMN with TPN anti-correlation. We reasoned that QPPs could provide new insights into AD network dysfunction and improve disease diagnosis. We therefore used rsfMRI to investigate QPPs in old TG2576 mice, a model of amyloidosis, and age-matched controls. Multiple QPPs were determined and compared across groups. Using linear regression, we removed their contribution from the functional scans and assessed how they reflected functional connectivity. Lastly, we used elastic net regression to determine if QPPs improved disease classification. We present three prominent findings: (1) Compared to controls, TG2576 mice were marked by opposing neural dynamics in which DMN areas were anti-correlated and displayed diminished anti-correlation with the TPN. (2) QPPs reflected lowered DMN functional connectivity in TG2576 mice and revealed significantly decreased DMN-TPN anti-correlations. (3) QPP-derived measures significantly improved classification compared to conventional functional connectivity measures. Altogether, our findings provide insight into the neural dynamics of aberrant network connectivity in AD and indicate that QPPs might serve as a translational diagnostic tool.