Temporal complexity of fMRI is reproducible and correlates with higher order cognition

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

BACKGROUND

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state functional magnetic resonance imaging (rs-fMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

METHODS

rs-fMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron emission tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

RESULTS

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta-power. Exploratory analyses revealed a close statistical relationship between LEN and positive symptom severity in patients.

CONCLUSION

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

Complexity organization of resting-state functional-MRI networks

Entropy measures are increasingly being used to analyze the structure of neural activity observed by functional magnetic resonance imaging (fMRI), with resting-state networks (RSNs) being of interest for their reproducible descriptions of the brain's functional architecture. Temporal correlations have shown a dichotomy among these networks: those that engage with the environment, known as extrinsic, which include the visual and sensorimotor networks; and those associated with executive control and self-referencing, known as intrinsic, which include the default mode network and the frontoparietal control network. While these inter-voxel temporal correlations enable the assessment of synchrony among the components of individual networks, entropic measures introduce an intra-voxel assessment that quantifies signal features encoded within each blood oxygen level-dependent (BOLD) time series. As a result, this framework offers insights into comprehending the representation and processing of information within fMRI signals. Multiscale entropy (MSE) has been proposed as a useful measure for characterizing the entropy of neural activity across different temporal scales. This measure of temporal entropy in BOLD data is dependent on the length of the time series; thus, high-quality data with fine-grained temporal resolution and a sufficient number of time frames is needed to improve entropy precision. We apply MSE to the Midnight Scan Club, a highly sampled and well-characterized publicly available dataset, to analyze the entropy distribution of RSNs and evaluate its ability to distinguish between different functional networks. Entropy profiles are compared across temporal scales and RSNs. Our results have shown that the spatial distribution of entropy at infra-slow frequencies (0.005-0.1 Hz) reproduces known parcellations of RSNs. We found a complexity hierarchy between intrinsic and extrinsic RSNs, with intrinsic networks robustly exhibiting higher entropy than extrinsic networks. Finally, we found new evidence that the topography of entropy in the posterior cerebellum exhibits high levels of entropy comparable to that of intrinsic RSNs.

The heritability and structural correlates of resting-state fMRI complexity

The complexity of fMRI signals quantifies temporal dynamics of spontaneous neural activity, which has been increasingly recognized as providing important insights into cognitive functions and psychiatric disorders. However, its heritability and structural underpinnings are not well understood. Here, we utilize multi-scale sample entropy to extract resting-state fMRI complexity in a large healthy adult sample from the Human Connectome Project. We show that fMRI complexity at multiple time scales is heritable in broad brain regions. Heritability estimates are modest and regionally variable. We relate fMRI complexity to brain structure including surface area, cortical myelination, cortical thickness, subcortical volumes, and total brain volume. We find that surface area is negatively correlated with fine-scale complexity and positively correlated with coarse-scale complexity in most cortical regions, especially the association cortex. Most of these correlations are related to common genetic and environmental effects. We also find positive correlations between cortical myelination and fMRI complexity at fine scales and negative correlations at coarse scales in the prefrontal cortex, lateral temporal lobe, precuneus, lateral parietal cortex, and cingulate cortex, with these correlations mainly attributed to common environmental effects. We detect few significant associations between fMRI complexity and cortical thickness. Despite the non-significant association with total brain volume, fMRI complexity exhibits significant correlations with subcortical volumes in the hippocampus, cerebellum, putamen, and pallidum at certain scales. Collectively, our work establishes the genetic basis and structural correlates of resting-state fMRI complexity across multiple scales, supporting its potential application as an endophenotype for psychiatric disorders.

Reliability of variability and complexity measures for task and task-free BOLD fMRI

Brain activity continuously fluctuates over time, even if the brain is in controlled (e.g., experimentally induced) states. Recent years have seen an increasing interest in understanding the complexity of these temporal variations, for example with respect to developmental changes in brain function or between-person differences in healthy and clinical populations. However, the psychometric reliability of brain signal variability and complexity measures-which is an important precondition for robust individual differences as well as longitudinal research-is not yet sufficiently studied. We examined reliability (split-half correlations) and test-retest correlations for task-free (resting-state) BOLD fMRI as well as split-half correlations for seven functional task data sets from the Human Connectome Project to evaluate their reliability. We observed good to excellent split-half reliability for temporal variability measures derived from rest and task fMRI activation time series (standard deviation, mean absolute successive difference, mean squared successive difference), and moderate test-retest correlations for the same variability measures under rest conditions. Brain signal complexity estimates (several entropy and dimensionality measures) showed moderate to good reliabilities under both, rest and task activation conditions. We calculated the same measures also for time-resolved (dynamic) functional connectivity time series and observed moderate to good reliabilities for variability measures, but poor reliabilities for complexity measures derived from functional connectivity time series. Global (i.e., mean across cortical regions) measures tended to show higher reliability than region-specific variability or complexity estimates. Larger subcortical regions showed similar reliability as cortical regions, but small regions showed lower reliability, especially for complexity measures. Lastly, we also show that reliability scores are only minorly dependent on differences in scan length and replicate our results across different parcellation and denoising strategies. These results suggest that the variability and complexity of BOLD activation time series are robust measures well-suited for individual differences research. Temporal variability of global functional connectivity over time provides an important novel approach to robustly quantifying the dynamics of brain function. PRACTITIONER POINTS: Variability and complexity measures of BOLD activation show good split-half reliability and moderate test-retest reliability. Measures of variability of global functional connectivity over time can robustly quantify neural dynamics. Length of fMRI data has only a minor effect on reliability.

Individual characteristics outperform resting-state fMRI for the prediction of behavioral phenotypes

In this study, we aimed to compare imaging-based features of brain function, measured by resting-state fMRI (rsfMRI), with individual characteristics such as age, gender, and total intracranial volume to predict behavioral measures. We developed a machine learning framework based on rsfMRI features in a dataset of 20,000 healthy individuals from the UK Biobank, focusing on temporal complexity and functional connectivity measures. Our analysis across four behavioral phenotypes revealed that both temporal complexity and functional connectivity measures provide comparable predictive performance. However, individual characteristics consistently outperformed rsfMRI features in predictive accuracy, particularly in analyses involving smaller sample sizes. Integrating rsfMRI features with demographic data sometimes enhanced predictive outcomes. The efficacy of different predictive modeling techniques and the choice of brain parcellation atlas were also examined, showing no significant influence on the results. To summarize, while individual characteristics are superior to rsfMRI in predicting behavioral phenotypes, rsfMRI still conveys additional predictive value in the context of machine learning, such as investigating the role of specific brain regions in behavioral phenotypes.

Sex differences of signal complexity at resting-state functional magnetic resonance imaging and their associations with the estrogen-signaling pathway in the brain

Sex differences in the brain have been widely reported and may hold the key to elucidating sex differences in many medical conditions and drug response. However, the molecular correlates of these sex differences in structural and functional brain measures in the human brain remain unclear. Herein, we used sample entropy (SampEn) to quantify the signal complexity of resting-state functional magnetic resonance imaging (rsfMRI) in a large neuroimaging cohort (N = 1,642). The frontoparietal control network and the cingulo-opercular network had high signal complexity while the cerebellar and sensory motor networks had low signal complexity in both men and women. Compared with those in male brains, we found greater signal complexity in all functional brain networks in female brains with the default mode network exhibiting the largest sex difference. Using the gene expression data in brain tissues, we identified genes that were significantly associated with sex differences in brain signal complexity. The significant genes were enriched in the gene sets that were differentially expressed between the brain cortex and other tissues, the estrogen-signaling pathway, and the biological function of neural plasticity. In particular, the G-protein-coupled estrogen receptor 1 gene in the estrogen-signaling pathway was expressed more in brain regions with greater sex differences in SampEn. In conclusion, greater complexity in female brains may reflect the interactions between sex hormone fluctuations and neuromodulation of estrogen in women.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-09954-y.

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

Background

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state fMRI (rfMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

Methods

rfMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron-Emission Tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

Results

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta power. Exploratory analyses revealed a close statistical relationship between LEN and positive PSD symptoms.

Conclusion

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

CRediT Authorship Contribution Statement

Fabian Hirsch: Conceptualization, Methodology, Software, Formal analysis, Writing - Original Draft, Writing - Review & Editing, Visualization; Ângelo Bumanglag: Methodology, Software, Formal analysis, Writing - Review & Editing; Yifei Zhang: Methodology, Software, Formal analysis, Writing - Review & Editing; Afra Wohlschlaeger: Methodology, Writing - Review & Editing, Supervision, Project administration.

Information decomposition and the informational architecture of the brain

To explain how the brain orchestrates information-processing for cognition, we must understand information itself. Importantly, information is not a monolithic entity. Information decomposition techniques provide a way to split information into its constituent elements: unique, redundant, and synergistic information. We review how disentangling synergistic and redundant interactions is redefining our understanding of integrative brain function and its neural organisation. To explain how the brain navigates the trade-offs between redundancy and synergy, we review converging evidence integrating the structural, molecular, and functional underpinnings of synergy and redundancy; their roles in cognition and computation; and how they might arise over evolution and development. Overall, disentangling synergistic and redundant information provides a guiding principle for understanding the informational architecture of the brain and cognition.

Multivariate multiscale entropy (mMSE) as a tool for understanding the resting-state EEG signal dynamics: the spatial distribution and sex/gender-related differences

BACKGROUND

The study aimed to determine how the resting-state EEG (rsEEG) complexity changes both over time and space (channels). The complexity of rsEEG and its sex/gender differences were examined using the multivariate Multiscale Entropy (mMSE) in 95 healthy adults. Following the probability maps (Giacometti et al. in J Neurosci Methods 229:84-96, 2014), channel sets have been identified that correspond to the functional networks. For each channel set the area under curve (AUC), which represents the total complexity, MaxSlope-the maximum complexity change of the EEG signal at thefine scales (1:4 timescales), and AvgEnt-to the average entropy level at coarse-grained scales (9:12 timescales), respectively, were extracted. To check dynamic changes between the entropy level at the fine and coarse-grained scales, the difference in mMSE between the #9 and #4 timescale (DiffEnt) was also calculated.

RESULTS

We found the highest AUC for the channel sets corresponding to the somatomotor (SMN), dorsolateral network (DAN) and default mode (DMN) whereas the visual network (VN), limbic (LN), and frontoparietal (FPN) network showed the lowest AUC. The largest MaxSlope were in the SMN, DMN, ventral attention network (VAN), LN and FPN, and the smallest in the VN. The SMN and DAN were characterized by the highest and the LN, FPN, and VN by the lowest AvgEnt. The most stable entropy were for the DAN and VN while the LN showed the greatest drop of entropy at the coarse scales. Women, compared to men, showed higher MaxSlope and DiffEnt but lower AvgEnt in all channel sets.

CONCLUSIONS

Novel results of the present study are: (1) an identification of the mMSE features that capture entropy at the fine and coarse timescales in the channel sets corresponding to the main resting-state networks; (2) the sex/gender differences in these features.

Relation of resting brain signal variability to cognitive and socioemotional measures in an adult lifespan sample

Temporal variability of the fMRI-derived blood-oxygen-level-dependent (BOLD) signal during cognitive tasks shows important associations with individual differences in age and performance. Less is known about relations between spontaneous BOLD variability measured at rest and relatively stable cognitive measures, such as IQ or socioemotional function. Here, we examined associations among resting BOLD variability, cognitive/socioemotional scores from the NIH Toolbox and optimal time of day for alertness (chronotype) in a sample of 157 adults from 20 to 86 years of age. To investigate individual differences in these associations independently of age, we regressed age out from both behavioral and BOLD variability scores. We hypothesized that greater BOLD variability would be related to higher fluid cognition scores, more positive scores on socioemotional scales and a morningness chronotype. Consistent with this idea, we found positive correlations between resting BOLD variability, positive socioemotional scores (e.g. self-efficacy) and morning chronotype, as well as negative correlations between variability and negative emotional scores (e.g. loneliness). Unexpectedly, we found negative correlations between BOLD variability and fluid cognition. These results suggest that greater resting brain signal variability facilitates optimal socioemotional function and characterizes those with morning-type circadian rhythms, but individuals with greater fluid cognition may be more likely to show less temporal variability in spontaneous measures of BOLD activity.

Accounting for temporal variability in functional magnetic resonance imaging improves prediction of intelligence

Neuroimaging-based prediction methods for intelligence have seen a rapid development. Among different neuroimaging modalities, prediction using functional connectivity (FC) has shown great promise. Most literature has focused on prediction using static FC, with limited investigations on the merits of such analysis compared to prediction using dynamic FC or region-level functional magnetic resonance imaging (fMRI) times series that encode temporal variability. To account for the temporal dynamics in fMRI, we propose a bi-directional long short-term memory (bi-LSTM) approach that incorporates feature selection mechanism. The proposed pipeline is implemented via an efficient algorithm and applied for predicting intelligence using region-level time series and dynamic FC. We compare the prediction performance using different fMRI features acquired from the Adolescent Brain Cognitive Development (ABCD) study involving nearly 7000 individuals. Our detailed analysis illustrates the consistently inferior performance of static FC compared to region-level time series or dynamic FC for single and combined rest and task fMRI experiments. The joint analysis of task and rest fMRI leads to improved intelligence prediction under all models compared to using fMRI from only one experiment. In addition, the proposed bi-LSTM pipeline based on region-level time series identifies several shared and differential important brain regions across fMRI experiments that drive intelligence prediction. A test-retest analysis of the selected regions shows strong reliability across cross-validation folds. Given the large sample size of ABCD study, our results provide strong evidence that superior prediction of intelligence can be achieved by accounting for temporal variations in fMRI.

Hierarchical fluctuation shapes a dynamic flow linked to states of consciousness

Consciousness arises from the spatiotemporal neural dynamics, however, its relationship with neural flexibility and regional specialization remains elusive. We identified a consciousness-related signature marked by shifting spontaneous fluctuations along a unimodal-transmodal cortical axis. This simple signature is sensitive to altered states of consciousness in single individuals, exhibiting abnormal elevation under psychedelics and in psychosis. The hierarchical dynamic reflects brain state changes in global integration and connectome diversity under task-free conditions. Quasi-periodic pattern detection revealed that hierarchical heterogeneity as spatiotemporally propagating waves linking to arousal. A similar pattern can be observed in macaque electrocorticography. Furthermore, the spatial distribution of principal cortical gradient preferentially recapitulated the genetic transcription levels of the histaminergic system and that of the functional connectome mapping of the tuberomammillary nucleus, which promotes wakefulness. Combining behavioral, neuroimaging, electrophysiological, and transcriptomic evidence, we propose that global consciousness is supported by efficient hierarchical processing constrained along a low-dimensional macroscale gradient.

Recommendations for validating hierarchical clustering in consumer sensory projects

Choosing the proper hierarchical clustering algorithm and number of clusters is always a key question in consumer sensory projects. In many cases, researchers do not publish any reason why it was chosen a given distance measure and linkage rule along with cluster numbers. The reason behind this could be that different cluster validation and comparison techniques give contradictory results in most cases. A complex evaluation to define the proper clustering might be time-consuming and tedious. The paper introduces the clustering of three sensory data sets using different distance metrics and linkage rules for different numbers of clusters. The results of the validation methods deviate, suggesting that clustering depends heavily on the data set in question. Although Euclidean distance, Ward's method seems a safe choice, testing, and validation of different clustering combinations is strongly suggested.

From local properties to brain-wide organization: A review of intraregional temporal features in functional magnetic resonance imaging data

Based on the fluctuations ensembled over neighbouring neurons, blood oxygen level-dependent (BOLD) signal is a mesoscale measurement of brain signals. Intraregional temporal features (IRTFs) of BOLD signal, extracted from regional neural activities, are utilized to investigate how the brain functions in local brain areas. This literature highlights four types of IRTFs and their representative calculations including variability in the temporal domain, variability in the frequency domain, entropy, and intrinsic neural timescales, which are tightly related to cognitions. In the brain-wide spatial organization, these brain features generally organized into two spatial hierarchies, reflecting structural constraints of regional dynamics and hierarchical functional processing workflow in brain. Meanwhile, the spatial organization gives rise to the link between neuronal properties and cognitive performance. Disrupted or unbalanced spatial conditions of IRTFs emerge with suboptimal cognitive states, which improved our understanding of the aging process and/or neuropathology of brain disease. This review concludes that IRTFs are important properties of the brain functional system and IRTFs should be considered in a brain-wide manner.

Temporal Stability of the Dynamic Resting-State Functional Brain Network: Current Measures, Clinical Research Progress, and Future Perspectives

Based on functional magnetic resonance imaging and multilayer dynamic network model, the brain network's quantified temporal stability has shown potential in predicting altered brain functions. This manuscript aims to summarize current knowledge, clinical research progress, and future perspectives on brain network's temporal stability. There are a variety of widely used measures of temporal stability such as the variance/standard deviation of dynamic functional connectivity strengths, the temporal variability, the flexibility (switching rate), and the temporal clustering coefficient, while there is no consensus to date which measure is the best. The temporal stability of brain networks may be associated with several factors such as sex, age, cognitive functions, head motion, circadian rhythm, and data preprocessing/analyzing strategies, which should be considered in clinical studies. Multiple common psychiatric disorders such as schizophrenia, major depressive disorder, and bipolar disorder have been found to be related to altered temporal stability, especially during the resting state; generally, both excessively decreased and increased temporal stabilities were thought to reflect disorder-related brain dysfunctions. However, the measures of temporal stability are still far from applications in clinical diagnoses for neuropsychiatric disorders partly because of the divergent results. Further studies with larger samples and in transdiagnostic (including schizoaffective disorder) subjects are warranted.

The complexity of the stream of consciousness

Typical consciousness can be defined as an individual-specific stream of experiences. Modern consciousness research on dynamic functional connectivity uses clustering techniques to create common bases on which to compare different individuals. We propose an alternative approach by combining modern theories of consciousness and insights arising from phenomenology and dynamical systems theory. This approach enables a representation of an individual's connectivity dynamics in an intrinsically-defined, individual-specific landscape. Given the wealth of evidence relating functional connectivity to experiential states, we assume this landscape is a proxy measure of an individual's stream of consciousness. By investigating the properties of this landscape in individuals in different states of consciousness, we show that consciousness is associated with short term transitions that are less predictable, quicker, but, on average, more constant. We also show that temporally-specific connectivity states are less easily describable by network patterns that are distant in time, suggesting a richer space of possible states. We show that the cortex, cerebellum and subcortex all display consciousness-relevant dynamics and discuss the implication of our results in forming a point of contact between dynamical systems interpretations and phenomenology.

On the Spatial Distribution of Temporal Complexity in Resting State and Task Functional MRI

Measuring the temporal complexity of functional MRI (fMRI) time series is one approach to assess how brain activity changes over time. In fact, hemodynamic response of the brain is known to exhibit critical behaviour at the edge between order and disorder. In this study, we aimed to revisit the spatial distribution of temporal complexity in resting state and task fMRI of 100 unrelated subjects from the Human Connectome Project (HCP). First, we compared two common choices of complexity measures, i.e., Hurst exponent and multiscale entropy, and observed a high spatial similarity between them. Second, we considered four tasks in the HCP dataset (Language, Motor, Social, and Working Memory) and found high task-specific complexity, even when the task design was regressed out. For the significance thresholding of brain complexity maps, we used a statistical framework based on graph signal processing that incorporates the structural connectome to develop the null distributions of fMRI complexity. The results suggest that the frontoparietal, dorsal attention, visual, and default mode networks represent stronger complex behaviour than the rest of the brain, irrespective of the task engagement. In sum, the findings support the hypothesis of fMRI temporal complexity as a marker of cognition.

Network dynamics scale with levels of awareness

Small world topologies are thought to provide a valuable insight into human brain organisation and consciousness. However, functional magnetic resonance imaging studies in consciousness have not yielded consistent results. Given the importance of dynamics for both consciousness and cognition, here we investigate how the diversity of small world dynamics (quantified by sample entropy; dSW-E¹) scales with decreasing levels of awareness (i.e., sedation and disorders of consciousness). Paying particular attention to result reproducibility, we show that dSW-E is a consistent predictor of levels of awareness even when controlling for the underlying functional connectivity dynamics. We find that dSW-E of subcortical and cortical areas are predictive, with the former showing higher and more robust effect sizes across analyses. We find that the network dynamics of intermodular communication in the cerebellum also have unique predictive power for levels of awareness. Consequently, we propose that the dynamic reorganisation of the functional information architecture, in particular of the subcortex, is a characteristic that emerges with awareness and has explanatory power beyond that of the complexity of dynamic functional connectivity.

Monofractal analysis of functional magnetic resonance imaging: An introductory review

The following review will aid readers in providing an overview of scale-free dynamics and monofractal analysis, as well as its applications and potential in functional magnetic resonance imaging (fMRI) neuroscience and clinical research. Like natural phenomena such as the growth of a tree or crashing ocean waves, the brain expresses scale-invariant, or fractal, patterns in neural signals that can be measured. While neural phenomena may represent both monofractal and multifractal processes and can be quantified with many different interrelated parameters, this review will focus on monofractal analysis using the Hurst exponent (H). Monofractal analysis of fMRI data is an advanced analysis technique that measures the complexity of brain signaling by quantifying its degree of scale-invariance. As such, the H value of the blood oxygenation level-dependent (BOLD) signal specifies how the degree of correlation in the signal may mediate brain functions. This review presents a brief overview of the theory of fMRI monofractal analysis followed by notable findings in the field. Through highlighting the advantages and challenges of the technique, the article provides insight into how to best conduct fMRI fractal analysis and properly interpret the findings with physiological relevance. Furthermore, we identify the future directions necessary for its progression towards impactful functional neuroscience discoveries and widespread clinical use. Ultimately, this presenting review aims to build a foundation of knowledge among readers to facilitate greater understanding, discussion, and use of this unique yet powerful imaging analysis technique.

Avalanche criticality in individuals, fluid intelligence, and working memory

The critical brain hypothesis suggests that efficient neural computation can be achieved through critical brain dynamics. However, the relationship between human cognitive performance and scale‐free brain dynamics remains unclear. In this study, we investigated the whole‐brain avalanche activity and its individual variability in the human resting‐state functional magnetic resonance imaging (fMRI) data. We showed that though the group‐level analysis was inaccurate because of individual variability, the subject wise scale‐free avalanche activity was significantly associated with maximal synchronization entropy of their brain activity. Meanwhile, the complexity of functional connectivity, as well as structure–function coupling, is maximized in subjects with maximal synchronization entropy. We also observed order–disorder phase transitions in resting‐state brain dynamics and found that there were longer times spent in the subcritical regime. These results imply that large‐scale brain dynamics favor the slightly subcritical regime of phase transition. Finally, we showed evidence that the neural dynamics of human participants with higher fluid intelligence and working memory scores are closer to criticality. We identified brain regions whose critical dynamics showed significant positive correlations with fluid intelligence performance and found that these regions were located in the prefrontal cortex and inferior parietal cortex, which were believed to be important nodes of brain networks underlying human intelligence. Our results reveal the possible role that avalanche criticality plays in cognitive performance and provide a simple method to identify the critical point and map cortical states on a spectrum of neural dynamics, ranging from subcriticality to supercriticality.

Intrinsic dynamics and topography of sensory input systems

The brain is continuously bombarded by external stimuli, which are processed in different input systems. The intrinsic features of these sensory input systems remain yet unclear. Investigating topography and dynamics of input systems is the goal of our study in order to better understand the intrinsic features that shape their neural processing. Using a functional magnetic resonance imaging dataset, we measured neural topography and dynamics of the input systems during rest and task states. Neural dynamics were probed by scale-free activity, measured with the power-law exponent (PLE), as well as by order/disorder as measured with sample entropy (SampEn). Our main findings during both rest and task states are: 1) differences in neural dynamics (PLE, SampEn) between regions within each of the three sensory input systems 2) differences in topography and dynamics among the three input systems; 3) PLE and SampEn correlate and, as demonstrated in simulation, show non-linear relationship in the critical range of PLE; 4) scale-free activity during rest mediates the transition of SampEn from rest to task as probed in a mediation model. We conclude that the sensory input systems are characterized by their intrinsic topographic and dynamic organization which, through scale-free activity, modulates their input processing.

Interindividual Signatures of fMRI Temporal Fluctuations

The complexity and variability of human brain activity, such as quantified from Functional Magnetic Resonance Imaging (fMRI) time series, have been widely studied as potential markers of healthy and pathological states. However, the extent to which fMRI temporal features exhibit stable markers of inter-individual differences in brain function across healthy young adults is currently an open question. In this study, we draw upon two widely used time-series measures-a nonlinear complexity measure (sample entropy; SampEn) and a spectral measure of low-frequency content (fALFF)-to capture dynamic properties of resting-state fMRI in a large sample of young adults from the Human Connectome Project. We observe that these two measures are closely related, and that both generate reproducible patterns across brain regions over four different fMRI runs, with intra-class correlations of up to 0.8. Moreover, we find that both metrics can uniquely differentiate subjects with high identification rates (ca. 89%). Canonical correlation analysis revealed a significant relationship between multivariate brain temporal features and behavioral measures. Overall, these findings suggest that regional profiles of fMRI temporal characteristics may provide stable markers of individual differences, and motivate future studies to further probe relationships between fMRI time series metrics and behavior.