Reliability of energy landscape analysis of resting-state functional MRI data.

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

State-transition dynamics of resting-state functional magnetic resonance imaging data: model comparison and test-to-retest analysis.

Electroencephalogram (EEG) microstate analysis entails finding dynamics of quasi-stable and generally recurrent discrete states in multichannel EEG time series data and relating properties of the estimated state-transition dynamics to observables such as cognition and behavior. While microstate analysis has been widely employed to analyze EEG data, its use remains less prevalent in functional magnetic resonance imaging (fMRI) data, largely due to the slower timescale of such data. In the present study, we extend various data clustering methods used in EEG microstate analysis to resting-state fMRI data from healthy humans to extract their state-transition dynamics. We show that the quality of clustering is on par with that for various microstate analyses of EEG data. We then develop a method for examining test-retest reliability of the discrete-state transition dynamics between fMRI sessions and show that the within-participant test-retest reliability is higher than between-participant test-retest reliability for different indices of state-transition dynamics, different networks, and different data sets. This result suggests that state-transition dynamics analysis of fMRI data could discriminate between different individuals and is a promising tool for performing fingerprinting analysis of individuals.

Quantifying the contribution of subject and group factors in brain activation.

Research in neuroscience often assumes universal neural mechanisms, but increasing evidence points toward sizeable individual differences in brain activations. What remains unclear is the extent of the idiosyncrasy and whether different types of analyses are associated with different levels of idiosyncrasy. Here we develop a new method for addressing these questions. The method consists of computing the within-subject reliability and subject-to-group similarity of brain activations and submitting these values to a computational model that quantifies the relative strength of group- and subject-level factors. We apply this method to a perceptual decision-making task (n = 50) and find that activations related to task, reaction time, and confidence are influenced equally strongly by group- and subject-level factors. Both group- and subject-level factors are dwarfed by a noise factor, though higher levels of smoothing increases their contributions relative to noise. Overall, our method allows for the quantification of group- and subject-level factors of brain activations and thus provides a more detailed understanding of the idiosyncrasy levels in brain activations.

Single trial variability in neural activity during a working memory task reveals multiple distinct information processing sequences.

Successful encoding, maintenance, and retrieval of information stored in working memory requires persistent coordination of activity among multiple brain regions. It is generally assumed that the pattern of such coordinated activity remains consistent for a given task. Thus, to separate this task-relevant signal from noise, multiple trials of the same task are completed, and the neural response is averaged across trials to generate an event-related potential (ERP). However, from trial to trial, the neuronal activity recorded with electroencephalogram (EEG) is actually spatially and temporally diverse, conflicting with the assumption of a single pattern of activity for a given task. Here, we show that variability in neuronal activity among single time-locked trials arises from the presence of multiple forms of stimulus dependent synchronized activity (i.e., distinct ERPs). We develop a data-driven classification method based on community detection to identify three discrete spatio-temporal clusters, or subtypes, of trials with different patterns of activation that are further associated with differences in decision-making processes. These results demonstrate that differences in the patterns of neural activity during working memory tasks represent fluctuations in the engagement of distinct brain networks and cognitive processes, suggesting that the brain can choose from multiple mechanisms to perform a given task.

Synchronization and Clustering in Complex Quadratic Networks.

Synchronization and clustering are well studied in the context of networks of oscillators, such as neuronal networks. However, this relationship is notoriously difficult to approach mathematically in natural, complex networks. Here, we aim to understand it in a canonical framework, using complex quadratic node dynamics, coupled in networks that we call complex quadratic networks (CQNs). We review previously defined extensions of the Mandelbrot and Julia sets for networks, focusing on the behavior of the node-wise projections of these sets and on describing the phenomena of node clustering and synchronization. One aspect of our work consists of exploring ties between a network's connectivity and its ensemble dynamics by identifying mechanisms that lead to clusters of nodes exhibiting identical or different Mandelbrot sets. Based on our preliminary analytical results (obtained primarily in two-dimensional networks), we propose that clustering is strongly determined by the network connectivity patterns, with the geometry of these clusters further controlled by the connection weights. Here, we first explore this relationship further, using examples of synthetic networks, increasing in size (from 3, to 5, to 20 nodes). We then illustrate the potential practical implications of synchronization in an existing set of whole brain, tractography-based networks obtained from 197 human subjects using diffusion tensor imaging. Understanding the similarities to how these concepts apply to CQNs contributes to our understanding of universal principles in dynamic networks and may help extend theoretical results to natural, complex systems.

Interneuron Desynchronization Precedes Seizures in a Mouse Model of Dravet Syndrome.

Recurrent seizures, which define epilepsy, are transient abnormalities in the electrical activity of the brain. The mechanistic basis of seizure initiation, and the contribution of defined neuronal subtypes to seizure pathophysiology, remains poorly understood. We performed in vivo two-photon calcium imaging in neocortex during temperature-induced seizures in male and female Dravet syndrome (Scn1a+/-) mice, a neurodevelopmental disorder with prominent temperature-sensitive epilepsy. Mean activity of both putative principal cells and parvalbumin-positive interneurons (PV-INs) was higher in Scn1a+/- relative to wild-type controls during quiet wakefulness at baseline and at elevated core body temperature. However, wild-type PV-INs showed a progressive synchronization in response to temperature elevation that was absent in PV-INs from Scn1a+/- mice. Hence, PV-IN activity remains intact interictally in Scn1a+/- mice, yet exhibits decreased synchrony immediately before seizure onset. We suggest that impaired PV-IN synchronization may contribute to the transition to the ictal state during temperature-induced seizures in Dravet syndrome.SIGNIFICANCE STATEMENT Epilepsy is a common neurological disorder defined by recurrent, unprovoked seizures. However, basic mechanisms of seizure initiation and propagation remain poorly understood. We performed in vivo two-photon calcium imaging in an experimental model of Dravet syndrome (Scn1a+/- mice)-a severe neurodevelopmental disorder defined by temperature-sensitive, treatment-resistant epilepsy-and record activity of putative excitatory neurons and parvalbumin-positive GABAergic neocortical interneurons (PV-INs) during naturalistic seizures induced by increased core body temperature. PV-IN activity was higher in Scn1a+/- relative to wild-type controls during quiet wakefulness. However, wild-type PV-INs showed progressive synchronization in response to temperature elevation that was absent in PV-INs from Scn1a+/- mice before seizure onset. Hence, impaired PV-IN synchronization may contribute to transition to seizure in Dravet syndrome.

Behavior can be decoded across the cortex when individual differences are considered.

Group-level analyses have typically associated behavioral signatures with a constrained set of brain areas. Here we show that two behavioral metrics - reaction time (RT) and confidence - can be decoded across the cortex when each individual is considered separately. Subjects (N=50) completed a perceptual decision-making task with confidence. We built models decoding trial-level RT and confidence separately for each subject using the activation patterns in one brain area at a time after splitting the entire cortex into 200 regions of interest (ROIs). At the group level, we replicated previous results by showing that both RT and confidence could be decoded from a small number of ROIs (12.0% and 3.5%, respectively). Critically, at the level of the individual, both RT and confidence could be decoded from most brain regions even after Bonferroni correction (90.0% and 72.5%, respectively). Surprisingly, we observed that many brain regions exhibited opposite brain-behavior relationships across individuals, such that, for example, higher activations predicted fast RTs in some subjects but slow RTs in others. These results were further replicated in a second dataset. Lastly, we developed a simple test to determine the robustness of decoding performance, which showed that several hundred trials per subject are required for robust decoding. These results show that behavioral signatures can be decoded from a much broader range of cortical areas than previously recognized and suggest the need to study the brain-behavior relationship at both the group and the individual level.

Multiple brain activation patterns for the same task.

Meaningful variation in internal states that impacts cognition and behavior remains challenging to discover and characterize. Here we leveraged trial-to-trial fluctuations in the brain-wide signal recorded using functional MRI to test if distinct sets of brain regions are activated on different trials when accomplishing the same task. Across three different perceptual decision-making experiments, we estimated the brain activations for each trial. We then clustered the trials based on their similarity using modularity-maximization, a data-driven classification method. In each experiment, we found multiple distinct but stable subtypes of trials, suggesting that the same task can be accomplished in the presence of widely varying brain activation patterns. Surprisingly, in all experiments, one of the subtypes exhibited strong activation in the default mode network, which is typically thought to decrease in activity during tasks that require externally focused attention. The remaining subtypes were characterized by activations in different task-positive areas. The default mode network subtype was characterized by behavioral signatures that were similar to the other subtypes exhibiting activation with task-positive regions. Finally, in a fourth experiment, we tested whether multiple activation patterns would also appear for a qualitatively different, working memory task. We again found multiple subtypes of trials with differential activation in frontoparietal control, dorsal attention, and ventral attention networks. Overall, these findings demonstrate that the same cognitive tasks are accomplished through multiple brain activation patterns.

Brain signatures indexing variation in internal processing during perceptual decision-making.

Brain activity is highly variable even while performing the same cognitive task with consequences for performance. Discovering, characterizing, and linking variability in brain activity to internal processes has primarily relied on experimentally inducing changes (e.g., via attention manipulation) to identify neuronal and behavioral consequences or studying spontaneous changes in ongoing brain dynamics. However, changes in internal processing could arise from many factors, such as variation in strategy or arousal, that are independent of experimental conditions. Here we utilize a data-driven clustering method based on modularity-maximation to identify consistent spatial-temporal EEG activity patterns across individual trials and relate this activity to behavioral performance. Subjects (N = 25) performed a motion direction discrimination task with six interleaved levels of motion coherence. Modularity-maximization based clustering identified two discrete spatial-temporal clusters, or subtypes, of trials with different patterns of brain activity. Surprisingly, even though Subtype 1 occurred more frequently with lower motion coherence, it was nonetheless associated with faster response times. Computational modeling suggests that Subtype 1 was characterized by a lower threshold for reaching a decision. These results highlight trial-to-trial variability in decision processes usually masked to experimenters and provide a method for identifying endogenous brain state variability relevant to cognition and behavior.

Brain signatures indexing variation in internal processing during perceptual decision-making.

Brain activity is highly variable during a task. Discovering, characterizing, and linking variability in brain activity to internal processes has primarily relied on experimental manipulations. However, changes in internal processing could arise from many factors independent of experimental conditions. Here we utilize a data-driven clustering method based on modularity-maximation to identify consistent spatial-temporal EEG activity patterns across individual trials. Subjects (N = 25) performed a motion discrimination task with six interleaved levels of coherence. Clustering identified two discrete subtypes of trials with different patterns of activity. Surprisingly, Subtype 1 occurred more frequently in trials with lower motion coherence but was associated with faster response times. Computational modeling suggests that Subtype 1 was characterized by a lower threshold for reaching a decision. These results highlight across-trial variability in decision processes traditionally hidden to experimenters and provide a method for identifying endogenous brain state variability relevant to cognition and behavior.

Reliability of energy landscape analysis of resting-state functional MRI data.

Energy landscape analysis is a data-driven method to analyze multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e., within-participant reliability) than across different sets of sessions from different participants (i.e., between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Within-subject reproducibility varies in multi-modal, longitudinal brain networks.

Network neuroscience provides important insights into brain function by analyzing complex networks constructed from diffusion Magnetic Resonance Imaging (dMRI), functional MRI (fMRI) and Electro/Magnetoencephalography (E/MEG) data. However, in order to ensure that results are reproducible, we need a better understanding of within- and between-subject variability over long periods of time. Here, we analyze a longitudinal, 8 session, multi-modal (dMRI, and simultaneous EEG-fMRI), and multiple task imaging data set. We first confirm that across all modalities, within-subject reproducibility is higher than between-subject reproducibility. We see high variability in the reproducibility of individual connections, but observe that in EEG-derived networks, during both rest and task, alpha-band connectivity is consistently more reproducible than connectivity in other frequency bands. Structural networks show a higher reliability than functional networks across network statistics, but synchronizability and eigenvector centrality are consistently less reliable than other network measures across all modalities. Finally, we find that structural dMRI networks outperform functional networks in their ability to identify individuals using a fingerprinting analysis. Our results highlight that functional networks likely reflect state-dependent variability not present in structural networks, and that the type of analysis should depend on whether or not one wants to take into account state-dependent fluctuations in connectivity.

Differential Patterns of Change in Brain Connectivity Resulting from Severe Traumatic Brain Injury.

Background: Traumatic brain injury (TBI) damages white matter tracts, disrupting brain network structure and communication. There exists a wide heterogeneity in the pattern of structural damage associated with injury, as well as a large heterogeneity in behavioral outcomes. However, little is known about the relationship between changes in network connectivity and clinical outcomes. Materials and Methods: We utilize the rat lateral fluid-percussion injury model of severe TBI to study differences in brain connectivity in 8 animals that received the insult and 11 animals that received only a craniectomy. Diffusion tensor imaging is performed 5 weeks after the injury and network theory is used to investigate changes in white matter connectivity. Results: We find that (1) global network measures are not able to distinguish between healthy and injured animals; (2) injury induced alterations predominantly exist in a subset of connections (subnetworks) distributed throughout the brain; and (3) injured animals can be divided into subgroups based on changes in network motifs-measures of local structural connectivity. In addition, alterations in predicted functional connectivity indicate that the subgroups have different propensities to synchronize brain activity, which could relate to the heterogeneity of clinical outcomes. Discussion: These results suggest that network measures can be used to quantify progressive changes in brain connectivity due to injury and differentiate among subpopulations with similar injuries, but different pathological trajectories.

Correlation of Histomorphometric Changes with Diffusion Tensor Imaging for Evaluation of Blast-Induced Auditory Neurodegeneration in Chinchilla.

In the present study, we have evaluated the blast-induced auditory neurodegeneration in chinchilla by correlating the histomorphometric changes with diffusion tensor imaging. The chinchillas were exposed to single unilateral blast-overpressure (BOP) at ∼172dB peak sound pressure level (SPL) and the pathological changes were compared at 1 week and 1 month after BOP. The functional integrity of the auditory system was assessed by auditory brainstem response (ABR) and distortion product otoacoustic emissions (DPOAE). The axonal integrity was assessed using diffusion tensor imaging at regions of interests (ROIs) of the central auditory neuraxis (CAN) including the cochlear nucleus (CN), inferior colliculus (IC), and auditory cortex (AC). Post-BOP, cyto-architecture metrics such as viable cells, degenerating neurons, and apoptotic cells were quantified at the CAN ROIs using light microscopic studies using cresyl fast violet, hematoxylin and eosin, and modified Crossmon's trichrome stains. We observed mean ABR threshold shifts of 30- and 10-dB SPL at 1 week and 1 month after BOP, respectively. A similar pattern was observed in DPAOE amplitudes shift. In the CAN ROIs, diffusion tensor imaging studies showed a decreased axial diffusivity in CN 1 month after BOP and a decreased mean diffusivity and radial diffusivity at 1 week after BOP. However, morphometric measures such as decreased viable cells and increased degenerating neurons and apoptotic cells were observed at CN, IC, and AC. Specifically, increased degenerating neurons and reduced viable cells were high on the ipsilateral side when compared with the contralateral side. These results indicate that a single blast significantly damages structural and functional integrity at all levels of CAN ROIs.

Personalized brain network models for assessing structure-function relationships.

Many recent efforts in computational modeling of macro-scale brain dynamics have begun to take a data-driven approach by incorporating structural and/or functional information derived from subject data. Here, we discuss recent work using personalized brain network models to study structure-function relationships in human brains. We describe the steps necessary to build such models and show how this computational approach can provide previously unobtainable information through the ability to perform virtual experiments. Finally, we present examples of how personalized brain network models can be used to gain insight into the effects of local stimulation and improve surgical outcomes in epilepsy.