The analysis of time‐varying activity and connectivity patterns (i.e., the chronnectome) using resting‐state magnetic resonance imaging has become an important part of ongoing neuroscience discussions. The majority of previous work has focused on variations of temporal coupling among fixed spatial nodes or transition of the dominant activity/connectivity pattern over time. Here, we introduce an approach to capture spatial dynamics within functional domains (FDs), as well as temporal dynamics within and between FDs. The approach models the brain as a hierarchical functional architecture with different levels of granularity, where lower levels have higher functional homogeneity and less dynamic behavior and higher levels have less homogeneity and more dynamic behavior. First, a high‐order spatial independent component analysis is used to approximate functional units. A functional unit is a pattern of regions with very similar functional activity over time. Next, functional units are used to construct FDs. Finally, functional modules (FMs) are calculated from FDs, providing an overall view of brain dynamics. Results highlight the spatial fluidity within FDs, including a broad spectrum of changes in regional associations, from strong coupling to complete decoupling. Moreover, FMs capture the dynamic interplay between FDs. Patients with schizophrenia show transient reductions in functional activity and state connectivity across several FDs, particularly the subcortical domain. Activity and connectivity differences convey unique information in many cases (e.g., the default mode) highlighting their complementarity information. The proposed hierarchical model to capture FD spatiotemporal variations provides new insight into the macroscale chronnectome and identifies changes hidden from existing approaches.
The brain is highly dynamic, reorganizing its activity at different interacting spatial and temporal scales, including variation within and between brain networks. The chronnectome is a model of the brain in which nodal activity and connectivity patterns change in fundamental and recurring ways over time. Most literature assumes fixed spatial nodes/networks, ignoring the possibility that spatial nodes/networks may vary in time. Here, we introduce an approach to calculate a spatially fluid chronnectome (called the spatial chronnectome for clarity), which focuses on the variations of networks coupling at the voxel level, and identify a novel set of spatially dynamic features. Results reveal transient spatially fluid interactions between intra‐ and internetwork relationships in which brain networks transiently merge and separate, emphasizing dynamic segregation and integration. Brain networks also exhibit distinct spatial patterns with unique temporal characteristics, potentially explaining a broad spectrum of inconsistencies in previous studies that assumed static networks. Moreover, we show anticorrelative connections to brain networks are transient as opposed to constant across the entire scan. Preliminary assessments using a multi‐site dataset reveal the ability of the approach to obtain new information and nuanced alterations that remain undetected during static analysis. Patients with schizophrenia (SZ) display transient decreases in voxel‐wise network coupling within visual and auditory networks, and higher intradomain coupling variability. In summary, the spatial chronnectome represents a new direction of research enabling the study of functional networks which are transient at the voxel level, and the identification of mechanisms for within‐ and between‐subject spatial variability.
more » « less- NSF-PAR ID:
- 10461178
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Human Brain Mapping
- Volume:
- 40
- Issue:
- 10
- ISSN:
- 1065-9471
- Page Range / eLocation ID:
- p. 3058-3077
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract There is ample evidence of atypical functional connectivity (FC) in autism spectrum disorders (ASDs). However, transient relationships between neural networks cannot be captured by conventional static FC analyses. Dynamic FC (dFC) approaches have been used to identify repeating, transient connectivity patterns (“states”), revealing spatiotemporal network properties not observable in static FC. Recent studies have found atypical dFC in ASDs, but questions remain about the nature of group differences in transient connectivity, and the degree to which states persist or change over time. This study aimed to: (a) describe and relate static and dynamic FC in typical development and ASDs, (b) describe group differences in transient states and compare them with static FC patterns, and (c) examine temporal stability and flexibility between identified states. Resting‐state functional magnetic resonance imaging (fMRI) data were collected from 62 ASD and 57 typically developing (TD) children and adolescents. Whole‐brain, data‐driven regions of interest were derived from group independent component analysis. Sliding window analysis and k‐means clustering were used to explore dFC and identify transient states. Across all regions, static overconnnectivity and increased variability over time in ASDs predominated. Furthermore, significant patterns of group differences emerged in two transient states that were not observed in the static FC matrix, with group differences in one state primarily involving sensory and motor networks, and in the other involving higher‐order cognition networks. Default mode network segregation was significantly reduced in ASDs in both states. Results highlight that dynamic approaches may reveal more nuanced transient patterns of atypical FC in ASDs.
-
Analysis of time-evolving data is crucial to understand the functioning of dynamic systems such as the brain. For instance, analysis of functional magnetic resonance imaging (fMRI) data collected during a task may reveal spatial regions of interest, and how they evolve during the task. However, capturing underlying spatial patterns as well as their change in time is challenging. The traditional approach in fMRI data analysis is to assume that underlying spatial regions of interest are static. In this article, using fractional amplitude of low-frequency fluctuations (fALFF) as an effective way to summarize the variability in fMRI data collected during a task, we arrange time-evolving fMRI data as a subjects by voxels by time windows tensor, and analyze the tensor using a tensor factorization-based approach called a PARAFAC2 model to reveal spatial dynamics. The PARAFAC2 model jointly analyzes data from multiple time windows revealing subject-mode patterns, evolving spatial regions (also referred to as networks) and temporal patterns. We compare the PARAFAC2 model with matrix factorization-based approaches relying on independent components, namely, joint independent component analysis (ICA) and independent vector analysis (IVA), commonly used in neuroimaging data analysis. We assess the performance of the methods in terms of capturing evolving networks through extensive numerical experiments demonstrating their modeling assumptions. In particular, we show that (i) PARAFAC2 provides a compact representation in all modes, i.e., subjects, time , and voxels , revealing temporal patterns as well as evolving spatial networks, (ii) joint ICA is as effective as PARAFAC2 in terms of revealing evolving networks but does not reveal temporal patterns, (iii) IVA's performance depends on sample size, data distribution and covariance structure of underlying networks. When these assumptions are satisfied, IVA is as accurate as the other methods, (iv) when subject-mode patterns differ from one time window to another, IVA is the most accurate. Furthermore, we analyze real fMRI data collected during a sensory motor task, and demonstrate that a component indicating statistically significant group difference between patients with schizophrenia and healthy controls is captured, which includes primary and secondary motor regions, cerebellum, and temporal lobe, revealing a meaningful spatial map and its temporal change.more » « less
-
De Lellis, Pietro (Ed.)Many temporal networks exhibit multiple system states, such as weekday and weekend patterns in social contact networks. The detection of such distinct states in temporal network data has recently been studied as it helps reveal underlying dynamical processes. A commonly used method is network aggregation over a time window, which aggregates a subsequence of multiple network snapshots into one static network. This method, however, necessarily discards temporal dynamics within the time window. Here we propose a new method for detecting dynamic states in temporal networks using connection series (i.e., time series of connection status) between nodes. Our method consists of the construction of connection series tensors over nonoverlapping time windows, similarity measurement between these tensors, and community detection in the similarity network of those time windows. Experiments with empirical temporal network data demonstrated that our method outperformed the conventional approach using simple network aggregation in revealing interpretable system states. In addition, our method allows users to analyze hierarchical temporal structures and to uncover dynamic states at different spatial/temporal resolutions.more » « less
-
Jonathan R. Whitlock (Ed.)more » « less
Introduction Understanding the neural code has been one of the central aims of neuroscience research for decades. Spikes are commonly referred to as the units of information transfer, but multi-unit activity (MUA) recordings are routinely analyzed in aggregate forms such as binned spike counts, peri-stimulus time histograms, firing rates, or population codes. Various forms of averaging also occur in the brain, from the spatial averaging of spikes within dendritic trees to their temporal averaging through synaptic dynamics. However, how these forms of averaging are related to each other or to the spatial and temporal units of information representation within the neural code has remained poorly understood.
Materials and methods In this work we developed NeuroPixelHD, a symbolic hyperdimensional model of MUA, and used it to decode the spatial location and identity of static images shown to
n = 9 mice in the Allen Institute Visual Coding—NeuroPixels dataset from large-scale MUA recordings. We parametrically varied the spatial and temporal resolutions of the MUA data provided to the model, and compared its resulting decoding accuracy.Results For almost all subjects, we found 125ms temporal resolution to maximize decoding accuracy for both the spatial location of Gabor patches (81 classes for patches presented over a 9×9 grid) as well as the identity of natural images (118 classes corresponding to 118 images) across the whole brain. This optimal temporal resolution nevertheless varied greatly between different regions, followed a sensory-associate hierarchy, and was significantly modulated by the central frequency of theta-band oscillations across different regions. Spatially, the optimal resolution was at either of two mesoscale levels for almost all mice: the area level, where the spiking activity of all neurons within each brain area are combined, and the population level, where neuronal spikes within each area are combined across fast spiking (putatively inhibitory) and regular spiking (putatively excitatory) neurons, respectively. We also observed an expected interplay between optimal spatial and temporal resolutions, whereby increasing the amount of averaging across one dimension (space or time) decreases the amount of averaging that is optimal across the other dimension, and vice versa.
Discussion Our findings corroborate existing empirical practices of spatiotemporal binning and averaging in MUA data analysis, and provide a rigorous computational framework for optimizing the level of such aggregations. Our findings can also synthesize these empirical practices with existing knowledge of the various sources of biological averaging in the brain into a new theory of neural information processing in which the
unit of information varies dynamically based on neuronal signal and noise correlations across space and time.
