Comparing functional connectivity matrices: A geometry-aware approach applied to participant identification

Separating group- and individual-level brain signatures in the newborn functional connectome: A deep learning approach

Recent studies indicate that differences in cognition among individuals may be partially attributed to unique brain wiring patterns. While functional connectivity (FC)-based fingerprinting has demonstrated high accuracy in identifying adults, early studies on neonates suggest that individualized FC signatures are absent. We posit that individual uniqueness is present in neonatal FC data and that conventional linear models fail to capture the rapid developmental trajectories characteristic of newborn brains. To explore this hypothesis, we employed a deep generative model, known as a variational autoencoder (VAE), leveraging two extensive public datasets: one comprising resting-state functional MRI (rs-fMRI) scans from 100 adults and the other from 464 neonates. VAE models trained on rs-fMRI from both adults and newborns produced superior age prediction performance (with r between predicted- and actual age ∼ 0.7) and individual identification accuracy (∼45 %) compared to models trained solely on adult or neonatal data. The VAE model also showed significantly higher individual identification accuracy than linear models (=10∼30 %). Importantly, the VAE differentiated connections reflecting age-related changes from those indicative of individual uniqueness, a distinction not possible with linear models. Moreover, we derived 20 latent variables, each corresponding to distinct patterns of cortical functional network (CFNs). These CFNs varied in their representation of brain maturation and individual signatures; notably, certain CFNs that failed to capture neurodevelopmental traits, in fact, exhibited individual signatures. CFNs associated with neonatal neurodevelopment predominantly encompassed unimodal regions such as visual and sensorimotor areas, whereas those linked to individual uniqueness spanned multimodal and transmodal brain regions. The VAE's capacity to extract features from rs-fMRI data beyond the capabilities of linear models positions it as a valuable tool for delineating cognitive traits inherent in rs-fMRI and exploring individualized imaging phenotypes.

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.

Tangent space functional reconfigurations in individuals at risk for alcohol use disorder

Human brain function dynamically adjusts to ever-changing stimuli from the external environment. Studies characterizing brain functional reconfiguration are nevertheless scarce. Here we present a principled mathematical framework to quantify brain functional reconfiguration when engaging and disengaging from a stop signal task (SST). We apply tangent space projection (a Riemannian geometry mapping technique) to transform functional connectomes (FCs) of 54 participants and quantify functional reconfiguration using the correlation distance of the resulting tangent-FCs. Our goal was to compare functional reconfigurations in individuals at risk for alcohol use disorder (AUD). We hypothesized that functional reconfigurations when transitioning to/from a task would be influenced by family history of alcohol use disorder (FHA) and other AUD risk factors. Multilinear regression models showed that engaging and disengaging functional reconfiguration were associated with FHA and recent drinking. When engaging in the SST after a rest condition, functional reconfiguration was negatively associated with recent drinking, while functional reconfiguration when disengaging from the SST was negatively associated with FHA. In both models, several other factors contributed to the functional reconfiguration. This study 1demonstrates that tangent-FCs can characterize task-induced functional reconfiguration, and that it is related to AUD risk.

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.

Extracting interpretable signatures of whole-brain dynamics through systematic comparison

The brain's complex distributed dynamics are typically quantified using a limited set of manually selected statistical properties, leaving the possibility that alternative dynamical properties may outperform those reported for a given application. Here, we address this limitation by systematically comparing diverse, interpretable features of both intra-regional activity and inter-regional functional coupling from resting-state functional magnetic resonance imaging (rs-fMRI) data, demonstrating our method using case-control comparisons of four neuropsychiatric disorders. Our findings generally support the use of linear time-series analysis techniques for rs-fMRI case-control analyses, while also identifying new ways to quantify informative dynamical fMRI structures. While simple statistical representations of fMRI dynamics performed surprisingly well (e.g., properties within a single brain region), combining intra-regional properties with inter-regional coupling generally improved performance, underscoring the distributed, multifaceted changes to fMRI dynamics in neuropsychiatric disorders. The comprehensive, data-driven method introduced here enables systematic identification and interpretation of quantitative dynamical signatures of multivariate time-series data, with applicability beyond neuroimaging to diverse scientific problems involving complex time-varying systems.

Decoding Autism: Uncovering patterns in brain connectivity through sparsity analysis with rs-fMRI data

BACKGROUND

In the realm of neuro-disorders, precise diagnosis and treatment rely heavily on objective imaging-based biomarker identification. This study employs a sparsity approach on resting-state fMRI to discern relevant brain region connectivity for predicting Autism.

NEW METHOD

The proposed methodology involves four key steps: (1) Utilizing three probabilistic brain atlases to extract functionally homogeneous brain regions from fMRI data. (2) Employing a hybrid approach of Graphical Lasso and Akaike Information Criteria to optimize sparse inverse covariance matrices for representing the brain functional connectivity. (3) Employing statistical techniques to scrutinize functional brain structures in Autism and Control subjects. (4) Implementing both autoencoder-based feature extraction and entire feature-based approach coupled with AI-based learning classifiers to predict Autism.

RESULTS

The ensemble classifier with the extracted feature set achieves a classification accuracy of 84.7% ± 0.3% using the MSDL atlas. Meanwhile, the 1D-CNN model, employing all features, exhibits superior classification accuracy of 88.6% ± 1.7% with the Smith 2009 (rsn70) atlas.

COMPARISON WITH EXISTING METHOD (S)

The proposed methodology outperforms the conventional correlation-based functional connectivity approach with a notably high prediction accuracy of more than 88%, whereas considering all direct and noisy indirect region-based functional connectivity, the traditional methods bound the prediction accuracy within 70% to 79%.

CONCLUSIONS

This study underscores the potential of sparsity-based FC analysis using rs-fMRI data as a prognostic biomarker for detecting Autism.

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.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Rapid fluctuations in functional connectivity of cortical networks encode spontaneous behavior

Experimental work across species has demonstrated that spontaneously generated behaviors are robustly coupled to variations in neural activity within the cerebral cortex. Functional magnetic resonance imaging data suggest that temporal correlations in cortical networks vary across distinct behavioral states, providing for the dynamic reorganization of patterned activity. However, these data generally lack the temporal resolution to establish links between cortical signals and the continuously varying fluctuations in spontaneous behavior observed in awake animals. Here, we used wide-field mesoscopic calcium imaging to monitor cortical dynamics in awake mice and developed an approach to quantify rapidly time-varying functional connectivity. We show that spontaneous behaviors are represented by fast changes in both the magnitude and correlational structure of cortical network activity. Combining mesoscopic imaging with simultaneous cellular-resolution two-photon microscopy demonstrated that correlations among neighboring neurons and between local and large-scale networks also encode behavior. Finally, the dynamic functional connectivity of mesoscale signals revealed subnetworks not predicted by traditional anatomical atlas-based parcellation of the cortex. These results provide new insights into how behavioral information is represented across the neocortex and demonstrate an analytical framework for investigating time-varying functional connectivity in neural networks.

Profiles of resting state functional connectivity in temporal lobe epilepsy associated with post-laser interstitial thermal therapy seizure outcomes and semiologies

Introduction

It is now understood that in focal epilepsy, impacted neural regions are not limited to the epileptogenic zone. As such, further investigation into the underlying functional connectivity (FC) patterns in those enduring Temporal Lobe Epilepsy (TLE) with Mesial Temporal Sclerosis (MTS) is imperative to understanding the intricacies of the disease.

Methods

The rsfMRIs of 17 healthy participants, 10 left-sided TLE-MTS patients with a pre-operative history of focal impaired awareness seizures (FIA), and 13 left-sided TLE-MTS patients with a pre-operative history of focal aware seizures (FA) were compared to determine the existence of distinct FC patterns with respect to seizure types. Similarly, the rsfMRIs of the above-mentioned healthy participants, 16 left-sided TLE-MTS individuals who were seizure-free (SF) 12 months postoperatively, and 16 left-sided TLE-MTS persons without seizure freedom (nSF) were interrogated. The ROI-to-ROI connectivity analysis included a total of 175 regions of interest (ROIs) and accounted for both age and duration of epileptic activity. Significant correlations were determined via two-sample t-tests and Bonferroni correction (α = 0.05).

Results

Comparisons of FA and FIA groups depicted significant correlations between the contralateral anterior cingulate gyrus, subgenual region, and the contralateral cerebellum, lobule III (p-value = 2.26e-4, mean z-score = -0.05 ± 0.28, T = -4.23). Comparisons of SF with nSF depicted two significantly paired-ROIs; the contralateral amygdala and the contralateral precuneus (p-value = 2.9e-5, mean z-score = -0.12 ± 0.19, T = 4.98), as well as the contralateral locus coeruleus and the ipsilateral intralaminar nucleus (p-value= 1.37e-4, mean z-score = 0.06 ± 0.17, T = -4.41).

Significance

FC analysis proves to be a lucrative modality for exploring unique signatures with respect to seizure types and postoperative outcomes. By furthering our understanding of the differences between epileptic phenotypes, we can achieve improvement in future treatment modalities not limited to targeting advancements.

Capturing functional connectomics using Riemannian partial least squares

For neurological disorders and diseases, functional and anatomical connectomes of the human brain can be used to better inform targeted interventions and treatment strategies. Functional magnetic resonance imaging (fMRI) is a non-invasive neuroimaging technique that captures spatio-temporal brain function through change in blood-oxygen-level-dependent (BOLD) signals over time. FMRI can be used to study the functional connectome through the functional connectivity matrix; that is, Pearson's correlation matrix between time series from the regions of interest of an fMRI image. One approach to analysing functional connectivity is using partial least squares (PLS), a multivariate regression technique designed for high-dimensional predictor data. However, analysing functional connectivity with PLS ignores a key property of the functional connectivity matrix; namely, these matrices are positive definite. To account for this, we introduce a generalisation of PLS to Riemannian manifolds, called R-PLS, and apply it to symmetric positive definite matrices with the affine invariant geometry. We apply R-PLS to two functional imaging datasets: COBRE, which investigates functional differences between schizophrenic patients and healthy controls, and; ABIDE, which compares people with autism spectrum disorder and neurotypical controls. Using the variable importance in the projection statistic on the results of R-PLS, we identify key functional connections in each dataset that are well represented in the literature. Given the generality of R-PLS, this method has the potential to investigate new functional connectomes in the brain, and with future application to structural data can open up further avenues of research in multi-modal imaging analysis.

Tangent functional connectomes uncover more unique phenotypic traits

Functional connectomes (FCs) containing pairwise estimations of functional couplings between pairs of brain regions are commonly represented by correlation matrices. As symmetric positive definite matrices, FCs can be transformed via tangent space projections, resulting into tangent-FCs. Tangent-FCs have led to more accurate models predicting brain conditions or aging. Motivated by the fact that tangent-FCs seem to be better biomarkers than FCs, we hypothesized that tangent-FCs have also a higher fingerprint. We explored the effects of six factors: fMRI condition, scan length, parcellation granularity, reference matrix, main-diagonal regularization, and distance metric. Our results showed that identification rates are systematically higher when using tangent-FCs across the "fingerprint gradient" (here including test-retest, monozygotic and dizygotic twins). Highest identification rates were achieved when minimally (0.01) regularizing FCs while performing tangent space projection using Riemann reference matrix and using correlation distance to compare the resulting tangent-FCs. Such configuration was validated in a second dataset (resting-state).

A joint subspace mapping between structural and functional brain connectomes

Understanding the connection between the brain's structural connectivity and its functional connectivity is of immense interest in computational neuroscience. Although some studies have suggested that whole brain functional connectivity is shaped by the underlying structure, the rule by which anatomy constraints brain dynamics remains an open question. In this work, we introduce a computational framework that identifies a joint subspace of eigenmodes for both functional and structural connectomes. We found that a small number of those eigenmodes are sufficient to reconstruct functional connectivity from the structural connectome, thus serving as low-dimensional basis function set. We then develop an algorithm that can estimate the functional eigen spectrum in this joint space from the structural eigen spectrum. By concurrently estimating the joint eigenmodes and the functional eigen spectrum, we can reconstruct a given subject's functional connectivity from their structural connectome. We perform elaborate experiments and demonstrate that the proposed algorithm for estimating functional connectivity from the structural connectome using joint space eigenmodes gives competitive performance as compared to the existing benchmark methods with better interpretability.

Cross-modality comparison between structural and metabolic networks in individual brain based on the Jensen-Shannon divergence method: a healthy Chinese population study

The study aimed to investigate the consistency and diversity between metabolic and structural brain networks at individual level constructed with divergence-based method in healthy Chinese population. The ¹⁸F-FDG PET and T1-weighted images of brain were collected from 209 healthy participants. The Jensen-Shannon divergence (JSD) was used to calculate metabolic or structural connectivities between any pair of brain regions and then individual brain networks were constructed. The global and regional topological properties of both networks were analyzed with graph theoretical analysis. Regional properties including nodal efficiency, degree, and betweenness centrality were used to define hub regions of networks. Cross-modality similarity of brain connectivity was analyzed with differential power (DP) analysis. The default mode network (DMN) had the largest number of brain connectivities with high DP values. The small-worldness indexes of metabolic and structural networks in all participants were greater than 1. The structural network showed higher assortativity and local efficiency than metabolic network, while hierarchy and global efficiency were greater in the metabolic network (all P < 0.001). Most of hubs in both networks were symmetrically spatial distributed in the regions of the DMN and subcortical nuclei including thalamus and amygdala, etc. The human brain presented small-world architecture both in perspective of individual metabolic and structural networks. There was a structural substrate that supported the brain to globally and efficiently integrate and process metabolic interaction across brain regions. The cross-modality cooperation or specialization in both networks might imply mechanisms of achieving higher-order brain functions.

Open-source statistical and data processing tools for wide-field optical imaging data in mice

Significance

Wide-field optical imaging (WOI) can produce concurrent hemodynamic and cell-specific calcium recordings across the entire cerebral cortex in animal models. There have been multiple studies using WOI to image mouse models with various environmental or genetic manipulations to understand various diseases. Despite the utility of pursuing mouse WOI alongside human functional magnetic resonance imaging (fMRI), and the multitude of analysis toolboxes in the fMRI literature, there is not an available open-source, user-friendly data processing and statistical analysis toolbox for WOI data.

Aim

To assemble a MATLAB toolbox for processing WOI data, as described and adapted to combine techniques from multiple WOI groups and fMRI.

Approach

We outline our MATLAB toolbox on GitHub with multiple data analysis packages and translate a commonly used statistical approach from the fMRI literature to the WOI data. To illustrate the utility of our MATLAB toolbox, we demonstrate the ability of the processing and analysis framework to detect a well-established deficit in a mouse model of stroke and plot activation areas during an electrical paw stimulus experiment.

Results

Our processing toolbox and statistical methods isolate a somatosensory-based deficit 3 days following photothrombotic stroke and cleanly localize sensory stimulus activations.

Conclusions

The toolbox presented here details an open-source, user-friendly compilation of WOI processing tools with statistical methods to apply to any biological question investigated with WOI techniques.

Revealing the spatiotemporal requirements for accurate subject identification with resting-state functional connectivity: a simultaneous fNIRS-fMRI study

SIGNIFICANCE

Brain fingerprinting refers to identifying participants based on their functional patterns. Despite its success with functional magnetic resonance imaging (fMRI), brain fingerprinting with functional near-infrared spectroscopy (fNIRS) still lacks adequate validation.

AIM

We investigated how fNIRS-specific acquisition features (limited spatial information and nonneural contributions) influence resting-state functional connectivity (rsFC) patterns at the intra-subject level and, therefore, brain fingerprinting.

APPROACH

We performed multiple simultaneous fNIRS and fMRI measurements in 29 healthy participants at rest. Data were preprocessed following the best practices, including the removal of motion artifacts and global physiology. The rsFC maps were extracted with the Pearson correlation coefficient. Brain fingerprinting was tested with pairwise metrics and a simple linear classifier.

RESULTS

Our results show that average classification accuracy with fNIRS ranges from 75% to 98%, depending on the number of runs and brain regions used for classification. Under the right conditions, brain fingerprinting with fNIRS is close to the 99.9% accuracy found with fMRI. Overall, the classification accuracy is more impacted by the number of runs and the spatial coverage than the choice of the classification algorithm.

CONCLUSIONS

This work provides evidence that brain fingerprinting with fNIRS is robust and reliable for extracting unique individual features at the intra-subject level once relevant spatiotemporal constraints are correctly employed.

Concussion Prone Scenarios: A Multi-Dimensional Exploration in Impact Directions, Brain Morphology, and Network Architectures Using Computational Models

While individual susceptibility to traumatic brain injury (TBI) has been speculated, past work does not provide an analysis considering how physical features of an individual's brain (e.g., brain size, shape), impact direction, and brain network features can holistically contribute to the risk of suffering a TBI from an impact. This work investigated each of these features simultaneously using computational modeling and analyses of simulated functional connectivity. Unlike the past studies that assess the severity of TBI based on the quantification of brain tissue damage (e.g., principal strain), we approached the brain as a complex network in which neuronal oscillations orchestrate to produce normal brain function (estimated by functional connectivity) and, to this end, both the anatomical damage location and its topological characteristics within the brain network contribute to the severity of brain function disruption and injury. To represent the variations in the population, we analyzed a publicly available database of brain imaging data and selected five distinct network architectures, seven different brain sizes, and three uniaxial head rotational conditions to study the consequences of 74 virtual impact scenarios. Results show impact direction produces the most significant change in connections across brain areas (structural connectome) and the functional coupling of activity across these brain areas (functional connectivity). Axial rotations were more injurious than those with sagittal and coronal rotations when the head kinematics were the same for each condition. When the impact direction was held constant, brain network architecture showed a significantly different vulnerability across axial and sagittal, but not coronal rotations. As expected, brain size significantly affected the expected change in structural and functional connectivity after impact. Together, these results provided groupings of predicted vulnerability to impact-a subgroup of male brain architectures exposed to axial impacts were most vulnerable, while a subgroup of female brain architectures was the most tolerant to the sagittal impacts studied. These findings lay essential groundwork for subject-specific analyses of concussion and provide invaluable guidance for designing personalized protection equipment.

Person-specific and precision neuroimaging: Current methods and future directions

Most neuroimaging studies of brain function analyze data in normalized space to identify regions of common activation across participants. These studies treat interindividual differences in brain organization as noise, but this approach can obscure important information about the brain's functional architecture. Recently, a number of studies have adopted a person-specific approach that aims to characterize these individual differences and explore their reliability and implications for behavior. A subset of these studies has taken a precision imaging approach that collects multiple hours of data from each participant to map brain function on a finer scale. In this review, we provide a broad overview of how person-specific and precision imaging techniques have used resting-state measures to examine individual differences in the brain's organization and their impact on behavior, followed by how task-based activity continues to add detail to these discoveries. We argue that person-specific and precision approaches demonstrate substantial promise in uncovering new details of the brain's functional organization and its relationship to behavior in many areas of cognitive neuroscience. We also discuss some current limitations in this new field and some new directions it may take.

Brainprints: identifying individuals from magnetoencephalograms

Magnetoencephalography (MEG) is used to study a wide variety of cognitive processes. Increasingly, researchers are adopting principles of open science and releasing their MEG data. While essential for reproducibility, sharing MEG data has unforeseen privacy risks. Individual differences may make a participant identifiable from their anonymized recordings. However, our ability to identify individuals based on these individual differences has not yet been assessed. Here, we propose interpretable MEG features to characterize individual difference. We term these features brainprints (brain fingerprints). We show through several datasets that brainprints accurately identify individuals across days, tasks, and even between MEG and Electroencephalography (EEG). Furthermore, we identify consistent brainprint components that are important for identification. We study the dependence of identifiability on the amount of data available. We also relate identifiability to the level of preprocessing and the experimental task. Our findings reveal specific aspects of individual variability in MEG. They also raise concerns about unregulated sharing of brain data, even if anonymized.

Improving Functional Connectome Fingerprinting with Degree-Normalization

Background: Functional connectivity quantifies the statistical dependencies between the activity of brain regions, measured using neuroimaging data such as functional magnetic resonance imaging (fMRI) blood-oxygenation-level dependent time series. The network representation of functional connectivity, called a functional connectome (FC), has been shown to contain an individual fingerprint allowing participants identification across consecutive testing sessions. Recently, researchers have focused on the extraction of these fingerprints, with potential applications in personalized medicine. Materials and Methods: In this study, we show that a mathematical operation denominated degree-normalization can improve the extraction of FC fingerprints. Degree-normalization has the effect of reducing the excessive influence of strongly connected brain areas in the whole-brain network. We adopt the differential identifiability framework and apply it to both original and degree-normalized FCs of 409 individuals from the Human Connectome Project, in resting-state and 7 fMRI tasks. Results: Our results indicate that degree-normalization systematically improves three fingerprinting metrics, namely differential identifiability, identification rate, and matching rate. Moreover, the results related to the matching rate metric suggest that individual fingerprints are embedded in a low-dimensional space. Discussion: The results suggest that low-dimensional functional fingerprints lie in part in weakly connected subnetworks of the brain and that degree-normalization helps uncovering them. This work introduces a simple mathematical operation that could lead to significant improvements in future FC fingerprinting studies.

An interdisciplinary computational model for predicting traumatic brain injury: Linking biomechanics and functional neural networks

The brain is a complex network consisting of neuron cell bodies in the gray matter and their axonal projections, forming the white matter tracts. These neurons are supported by an equally complex vascular network as well as glial cells. Traumatic brain injury (TBI) can lead to the disruption of the structural and functional brain networks due to disruption of both neuronal cell bodies in the gray matter as well as their projections and supporting cells. To explore how an impact can alter the function of brain networks, we integrated a finite element (FE) brain mechanics model with linked models of brain dynamics (Kuramoto oscillator) and vascular perfusion (Balloon-Windkessel) in this study. We used empirical resting-state functional magnetic resonance imaging (MRI) data to optimize the fit of our brain dynamics and perfusion models to clinical data. Results from the FE model were used to mimic injury in these optimized brain dynamics models: injury to the nodes (gray matter) led to a decrease in the nodal oscillation frequency, while damage to the edges (axonal connections/white matter) progressively decreased coupling among connected nodes. A total of 53 cases, including 33 non-injurious and 20 concussive head impacts experienced by professional American football players were simulated using this integrated model. We examined the correlation of injury outcomes with global measures of structural connectivity, neural dynamics, and functional connectivity of the brain networks when using different lesion methods. Results show that injurious head impacts cause significant alterations in global network topology regardless of lesion methods. Changes between the disrupted and healthy functional connectivity (measured by Pearson correlation) consistently correlated well with injury outcomes (AUC≥0.75), although the predictive performance is not significantly different (p>0.05) to that of traditional kinematic measures (angular acceleration). Intriguingly, our lesion model for gray matter damage predicted increases in global efficiency and clustering coefficient with increases in injury risk, while disrupting axonal connections led to lower network efficiency and clustering. When both injury mechanisms were combined into a single injury prediction model, the injury prediction performance depended on the thresholds used to determine neurodegeneration and mechanical tolerance for axonal injury. Together, these results point towards complex effects of mechanical trauma to the brain and provide a new framework for understanding brain injury at a causal mechanistic level and developing more effective diagnostic methods and therapeutic interventions.

Embedding Functional Brain Networks in Low Dimensional Spaces Using Manifold Learning Techniques

Background: fMRI data is inherently high-dimensional and difficult to visualize. A recent trend has been to find spaces of lower dimensionality where functional brain networks can be projected onto manifolds as individual data points, leading to new ways to analyze and interpret the data. Here, we investigate the potential of two powerful non-linear manifold learning techniques for functional brain networks representation: (1) T-stochastic neighbor embedding (t-SNE) and (2) Uniform Manifold Approximation Projection (UMAP) a recent breakthrough in manifold learning. Methods: fMRI data from the Human Connectome Project (HCP) and an independent study of aging were used to generate functional brain networks. We used fMRI data collected during resting state data and during a working memory task. The relative performance of t-SNE and UMAP were investigated by projecting the networks from each study onto 2D manifolds. The levels of discrimination between different tasks and the preservation of the topology were evaluated using different metrics. Results: Both methods effectively discriminated the resting state from the memory task in the embedding space. UMAP discriminated with a higher classification accuracy. However, t-SNE appeared to better preserve the topology of the high-dimensional space. When networks from the HCP and aging studies were combined, the resting state and memory networks in general aligned correctly. Discussion: Our results suggest that UMAP, a more recent development in manifold learning, is an excellent tool to visualize functional brain networks. Despite dramatic differences in data collection and protocols, networks from different studies aligned correctly in the embedding space.

Coordinated multivoxel coding beyond univariate effects is not likely to be observable in fMRI data

Simultaneous recording of activity across brain regions can contain additional information compared to regional recordings done in isolation. In particular, multivariate pattern analysis (MVPA) across voxels has been interpreted as evidence for distributed coding of cognitive or sensorimotor processes beyond what can be gleaned from a collection of univariate effects (UVE) using functional magnetic resonance imaging (fMRI). Here, we argue that regardless of patterns revealed, conventional MVPA is merely a decoding tool with increased sensitivity arising from considering a large number of 'weak classifiers' (i.e., single voxels) in higher dimensions. We propose instead that 'real' multivoxel coding should result in changes in higher-order statistics across voxels between conditions such as second-order multivariate effects (sMVE). Surprisingly, analysis of conditions with robust multivariate effects (MVE) revealed by MVPA failed to show significant sMVE in two species (humans and macaques). Further analysis showed that while both MVE and sMVE can be readily observed in the spiking activity of neuronal populations, the slow and nonlinear hemodynamic coupling and low spatial resolution of fMRI activations make the observation of higher-order statistics between voxels highly unlikely. These results reveal inherent limitations of fMRI signals for studying coordinated coding across voxels. Together, these findings suggest that care should be taken in interpreting significant MVPA results as representing anything beyond a collection of univariate effects.

Network alignment and similarity reveal atlas-based topological differences in structural connectomes

The interactions between different brain regions can be modeled as a graph, called connectome, whose nodes correspond to parcels from a predefined brain atlas. The edges of the graph encode the strength of the axonal connectivity between regions of the atlas that can be estimated via diffusion magnetic resonance imaging (MRI) tractography. Herein, we aim to provide a novel perspective on the problem of choosing a suitable atlas for structural connectivity studies by assessing how robustly an atlas captures the network topology across different subjects in a homogeneous cohort. We measure this robustness by assessing the alignability of the connectomes, namely the possibility to retrieve graph matchings that provide highly similar graphs. We introduce two novel concepts. First, the graph Jaccard index (GJI), a graph similarity measure based on the well-established Jaccard index between sets; the GJI exhibits natural mathematical properties that are not satisfied by previous approaches. Second, we devise WL-align, a new technique for aligning connectomes obtained by adapting the Weisfeiler-Leman (WL) graph-isomorphism test. We validated the GJI and WL-align on data from the Human Connectome Project database, inferring a strategy for choosing a suitable parcellation for structural connectivity studies. Code and data are publicly available.

An important part of our current understanding of the structure of the human brain relies on the concept of brain network, which is obtained by looking at how different brain regions are connected with each other. In this paper we present a strategy for choosing a suitable parcellation of the brain for structural connectivity studies by making use of the concepts of network alignment and similarity. To do so, we design a novel similarity measure between weighted networks called graph Jaccard index, and a new network alignment technique called WL-align. By assessing the possibility to retrieve graph matchings that provide highly similar graphs, we show that morphology- and structure-based atlases define brain networks that are more topologically robust across a wide range of resolutions.

Representation learning of resting state fMRI with variational autoencoder

Resting state functional magnetic resonance imaging (rsfMRI) data exhibits complex but structured patterns. However, the underlying origins are unclear and entangled in rsfMRI data. Here we establish a variational auto-encoder, as a generative model trainable with unsupervised learning, to disentangle the unknown sources of rsfMRI activity. After being trained with large data from the Human Connectome Project, the model has learned to represent and generate patterns of cortical activity and connectivity using latent variables. The latent representation and its trajectory represent the spatiotemporal characteristics of rsfMRI activity. The latent variables reflect the principal gradients of the latent trajectory and drive activity changes in cortical networks. Representational geometry captured as covariance or correlation between latent variables, rather than cortical connectivity, can be used as a more reliable feature to accurately identify subjects from a large group, even if only a short period of data is available in each subject. Our results demonstrate that VAE is a valuable addition to existing tools, particularly suited for unsupervised representation learning of resting state fMRI activity.

A simple permutation-based test of intermodal correspondence

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

The instability of functional connectomes across the first year of life

The uniqueness and stability of the adolescent and adult functional connectome has been demonstrated to be high (80-95 % identification) using connectome-based identification (ID) or "fingerprinting". However, it is unclear to what extent individuals exhibit similar distinctiveness and stability in infancy, a developmental period of rapid and unparalleled brain development. In this study, we examined connectome-based ID rates within and across the first year of life using a longitudinal infant dataset at 1.5 month and 9 months of age. We also calculated the test-retest reliability of individual connections across the first year of life using the intraclass correlation coefficient (ICC). Overall, we found substantially lower infant ID rates than have been reported in adult and adolescent populations. Within-session ID rates were moderate and significant (ID = 48.94-70.83 %). Between-session ID rates were very low and not significant, with task-to-task connectomes resulting in the highest between-session ID rate (ID = 26.6 %). Similarly, average edge-level test-retest reliability was higher within-session than between-session (mean within-session ICC = 0.17, mean between-session ICC = 0.10). These findings suggest a lack of uniqueness and stability in functional connectomes across the first year of life consistent with the unparalleled changes in brain functional organization during this critical period.

Which deep learning model can best explain object representations of within-category exemplars?

Deep neural network (DNN) models realize human-equivalent performance in tasks such as object recognition. Recent developments in the field have enabled testing the hierarchical similarity of object representation between the human brain and DNNs. However, the representational geometry of object exemplars within a single category using DNNs is unclear. In this study, we investigate which DNN model has the greatest ability to explain invariant within-category object representations by computing the similarity between representational geometries of visual features extracted at the high-level layers of different DNN models. We also test for the invariability of within-category object representations of these models by identifying object exemplars. Our results show that transfer learning models based on ResNet50 best explained both within-category object representation and object identification. These results suggest that the invariability of object representations in deep learning depends not on deepening the neural network but on building a better transfer learning model.

Interindividual variability of functional connectome in schizophrenia

Schizophrenia is a complex psychiatric disorder that displays an outstanding interindividual variability in clinical manifestation and neurobiological substrates. A better characterization and quantification of this heterogeneity could guide the search for both common abnormalities (linked to lower intersubject variability) and the presence of biological subtypes (leading to a greater heterogeneity across subjects). In the current study, we address interindividual variability in functional connectome by means of resting-state fMRI in a large sample of patients with schizophrenia and healthy controls. Among the different metrics of distance/dissimilarity used to assess variability, geodesic distance showed robust results to head motion. The main findings of the current study point to (i) a higher between subject heterogeneity in the functional connectome of patients, (ii) variable levels of heterogeneity throughout the cortex, with greater variability in frontoparietal and default mode networks, and lower variability in the salience network, and (iii) an association of whole-brain variability with levels of clinical symptom severity and with topological properties of brain networks, suggesting that the average functional connectome overrepresents those patients with lower functional integration and with more severe clinical symptoms. Moreover, after performing a graph theoretical analysis of brain networks, we found that patients with more severe clinical symptoms had decreased connectivity at both whole-brain level and within the salience network, and that patients with higher negative symptoms had large-scale functional integration deficits.

Construction of embedded fMRI resting-state functional connectivity networks using manifold learning

We construct embedded functional connectivity networks (FCN) from benchmark resting-state functional magnetic resonance imaging (rsfMRI) data acquired from patients with schizophrenia and healthy controls based on linear and nonlinear manifold learning algorithms, namely, Multidimensional Scaling, Isometric Feature Mapping, Diffusion Maps, Locally Linear Embedding and kernel PCA. Furthermore, based on key global graph-theoretic properties of the embedded FCN, we compare their classification potential using machine learning. We also assess the performance of two metrics that are widely used for the construction of FCN from fMRI, namely the Euclidean distance and the cross correlation metric. We show that diffusion maps with the cross correlation metric outperform the other combinations.

Structure-informed functional connectivity driven by identifiable and state-specific control regions

Describing how the brain anatomical wiring contributes to the emergence of coordinated neural activity underlying complex behavior remains challenging. Indeed, patterns of remote coactivations that adjust with the ongoing task-demand do not systematically match direct, static anatomical links. Here, we propose that observed coactivation patterns, known as functional connectivity (FC), can be explained by a controllable linear diffusion dynamics defined on the brain architecture. Our model, termed structure-informed FC, is based on the hypothesis that different sets of brain regions controlling the information flow on the anatomical wiring produce state-specific functional patterns. We thus introduce a principled framework for the identification of potential control centers in the brain. We find that well-defined, sparse, and robust sets of control regions, partially overlapping across several tasks and resting state, produce FC patterns comparable to empirical ones. Our findings suggest that controllability is a fundamental feature allowing the brain to reach different states.

Understanding how brain anatomy promotes particular patterns of coactivations among neural regions is a key challenge in neuroscience. This challenge can be addressed using network science and systems theory. Here, we propose that coactivations result from the diffusion of information through the network of anatomical links connecting brain regions, with certain regions controlling the dynamics. We translate this hypothesis into a model called structure-informed functional connectivity, and we introduce a framework for identifying control regions based on empirical data. We find that our model produces coactivation patterns comparable to empirical ones, and that distinct sets of control regions are associated with different functional states. These findings suggest that controllability is an important feature allowing the brain to reach different states.

Smooth graph learning for functional connectivity estimation

Functional connectivity (FC) estimated from functional magnetic resonance imaging (fMRI) signals is important in understanding neural representation and information processing in cortical networks. However, due to a lack of "ground truth" FC pattern, the reliability and robustness of FC estimates are usually examined in downstream FC analysis tasks, such as performing participant's identification (also known as "fingerprinting"). In this paper, we propose to learn FC via a smooth graph learning framework. In particular, we treat each time frame of the fMRI time series as a graph signal on an underlying functional brain graph, and estimate the smooth graph functional connectivity (SGFC) by learning the weighted graph adjacency matrix based on graph signal smoothness assumption. We demonstrate that our approach gives rise to a natural and sparse graph representation of FC from which reliable graph measures can be extracted. Reliability of SGFC is evaluated in the context of fingerprinting and compared to correlation FC (CFC). SGFC achieves higher fingerprinting accuracy across several different experiment settings; the improvement is even more significant when a shorter fMRI scanning length is used for FC estimation. In addition to being reliable, we also validate the cognitive relevance of SGFC by using it to predict fluid intelligence. Finally, in evaluating topological measures of the sparse graph, SGFC reveals a more small-world and modular structure compared to CFC. Together, our results suggest that the smooth graph learning framework produces a naturally sparse, reliable, and cognitive-relevant representation of functional connectivity.

Using distance on the Riemannian manifold to compare representations in brain and in models

Representational similarity analysis (RSA) summarizes activity patterns for a set of experimental conditions into a matrix composed of pairwise comparisons between activity patterns. Two examples of such matrices are the condition-by-condition inner product and correlation matrix. These representational matrices reside on the manifold of positive semidefinite matrices, called the Riemannian manifold. We hypothesize that representational similarities would be more accurately quantified by considering the underlying manifold of the representational matrices. Thus, we introduce the distance on the Riemannian manifold as a metric for comparing representations. Analyzing simulated and real fMRI data and considering a wide range of metrics, we show that the Riemannian distance is least susceptible to sampling bias, results in larger intra-subject reliability, and affords searchlight mapping with high sensitivity and specificity. Furthermore, we show that the Riemannian distance can be used for measuring multi-dimensional connectivity. This measure captures both univariate and multivariate connectivity and is also more sensitive to nonlinear regional interactions compared to the state-of-the-art measures. Applying our proposed metric to neural network representations of natural images, we demonstrate that it also possesses outstanding performance in quantifying similarity in models. Taken together, our results lend credence to the proposition that RSA should consider the manifold of the representational matrices to summarize response patterns in the brain and in models.

Feature selection framework for functional connectome fingerprinting

The ability to uniquely characterize individual subjects based on their functional connectome (FC) is a key requirement for progress toward precision psychiatry. FC fingerprinting is increasingly studied in the neuroimaging community for this purpose, where a variety of approaches have been developed for effective FC fingerprinting. Recent independent studies showed that fingerprinting accuracy suffers at large sample sizes and when coarser parcellations are used for computing the FC. Quantifying this problem and understanding the reasons these factors impact fingerprinting accuracy is crucial to develop more accurate fingerprinting methods for large sample sizes. Part of the challenge in fingerprinting is that FC captures both generic and subject‐specific information. A systematic approach for identifying subject‐specific FC information is crucial for making progress in addressing the fingerprinting problem. In this study, we addressed three gaps in our understanding of the FC fingerprinting problem. First, we studied the joint effect of sample size and parcellation granularity. Second, we explained the reason for reduced fingerprinting accuracy with increased sample size and reduced parcellation granularity. To this end, we used a clustering quality metric from the data mining community. Third, we developed a general feature selection framework for systematically identifying resting‐state functional connectivity (RSFC) elements that capture information to uniquely identify subjects. In sum, we evaluated six different approaches from this framework by quantifying both subject‐specific fingerprinting accuracy and the decrease in accuracy with an increase in sample size to identify which approach improved quality metrics the most.

Subject identification using edge-centric functional connectivity

Group-level studies do not capture individual differences in network organization, an important prerequisite for understanding neural substrates shaping behavior and for developing interventions in clinical conditions. Recent studies have employed 'fingerprinting' analyses on functional connectivity to identify subjects' idiosyncratic features. Here, we develop a complementary approach based on an edge-centric model of functional connectivity, which focuses on the co-fluctuations of edges. We first show whole-brain edge functional connectivity (eFC) to be a robust substrate that improves identifiability over nodal FC (nFC) across different datasets and parcellations. Next, we characterize subjects' identifiability at different spatial scales, from single nodes to the level of functional systems and clusters using k-means clustering. Across spatial scales, we find that heteromodal brain regions exhibit consistently greater identifiability than unimodal, sensorimotor, and limbic regions. Lastly, we show that identifiability can be further improved by reconstructing eFC using specific subsets of its principal components. In summary, our results highlight the utility of the edge-centric network model for capturing meaningful subject-specific features and sets the stage for future investigations into individual differences using edge-centric models.

Geodesic Distance on Optimally Regularized Functional Connectomes Uncovers Individual Fingerprints

Background: Functional connectomes (FCs) have been shown to provide a reproducible individual fingerprint, which has opened the possibility of personalized medicine for neuro/psychiatric disorders. Thus, developing accurate ways to compare FCs is essential to establish associations with behavior and/or cognition at the individual level. Methods: Canonically, FCs are compared using Pearson's correlation coefficient of the entire functional connectivity profiles. Recently, it has been proposed that the use of geodesic distance is a more accurate way of comparing FCs, one which reflects the underlying non-Euclidean geometry of the data. Computing geodesic distance requires FCs to be positive-definite and hence invertible matrices. As this requirement depends on the functional magnetic resonance imaging scanning length and the parcellation used, it is not always attainable and sometimes a regularization procedure is required. Results: In the present work, we show that regularization is not only an algebraic operation for making FCs invertible, but also that an optimal magnitude of regularization leads to systematically higher fingerprints. We also show evidence that optimal regularization is data set-dependent and varies as a function of condition, parcellation, scanning length, and the number of frames used to compute the FCs. Discussion: We demonstrate that a universally fixed regularization does not fully uncover the potential of geodesic distance on individual fingerprinting and indeed could severely diminish it. Thus, an optimal regularization must be estimated on each data set to uncover the most differentiable across-subject and reproducible within-subject geodesic distances between FCs. The resulting pairwise geodesic distances at the optimal regularization level constitute a very reliable quantification of differences between subjects.

A Deep Learning Approach to Predict Autism Spectrum Disorder Using Multisite Resting-State fMRI

Autism spectrum disorder (ASD) is a complex and degenerative neuro-developmental disorder. Most of the existing methods utilize functional magnetic resonance imaging (fMRI) to detect ASD with a very limited dataset which provides high accuracy but results in poor generalization. To overcome this limitation and to enhance the performance of the automated autism diagnosis model, in this paper, we propose an ASD detection model using functional connectivity features of resting-state fMRI data. Our proposed model utilizes two commonly used brain atlases, Craddock 200 (CC200) and Automated Anatomical Labelling (AAL), and two rarely used atlases Bootstrap Analysis of Stable Clusters (BASC) and Power. A deep neural network (DNN) classifier is used to perform the classification task. Simulation results indicate that the proposed model outperforms state-of-the-art methods in terms of accuracy. The mean accuracy of the proposed model was 88%, whereas the mean accuracy of the state-of-the-art methods ranged from 67% to 85%. The sensitivity, F1-score, and area under receiver operating characteristic curve (AUC) score of the proposed model were 90%, 87%, and 96%, respectively. Comparative analysis on various scoring strategies show the superiority of BASC atlas over other aforementioned atlases in classifying ASD and control.

Integrating new memories into the hippocampal network activity space

By investigating the topology of neuronal co-activity, we found that mnemonic information spans multiple operational axes in the mouse hippocampus network. High-activity principal cells form the core of each memory along a first axis, segregating spatial contexts and novelty. Low-activity cells join co-activity motifs across behavioral events and enable their crosstalk along two other axes. This reveals an organizational principle for continuous integration and interaction of hippocampal memories.

Identifiability of structural networks of nonlinear electronic oscillators

The interplay between structure and function is critical in the understanding of complex systems, their dynamics and their behavior. We investigated the interplay between structural and functional networks by means of the differential identifiability framework, which here quantifies the ability of identifying a particular network structure based on (1) the observation of its functional network and (2) the comparison with a prior observation under different initial conditions. We carried out an experiment consisting of the construction of [Formula: see text] different structural networks composed of [Formula: see text] nonlinear electronic circuits and studied the regions where network structures are identifiable. Specifically, we analyzed how differential identifiability is related to the coupling strength between dynamical units (modifying the level of synchronization) and what are the consequences of increasing the amount of noise existing in the functional networks. We observed that differential identifiability reaches its highest value for low to intermediate coupling strengths. Furthermore, it is possible to increase the identifiability parameter by including a principal component analysis in the comparison of functional networks, being especially beneficial for scenarios where noise reaches intermediate levels. Finally, we showed that the regime of the parameter space where differential identifiability is the highest is highly overlapped with the region where structural and functional networks correlate the most.

Typicality of functional connectivity robustly captures motion artifacts in rs-fMRI across datasets, atlases, and preprocessing pipelines

Functional connectivity analysis of resting-state fMRI data has recently become one of the most common approaches to characterizing individual brain function. It has been widely suggested that the functional connectivity matrix is a useful approximate representation of the brain's connectivity, potentially providing behaviorally or clinically relevant markers. However, functional connectivity estimates are known to be detrimentally affected by various artifacts, including those due to in-scanner head motion. Moreover, as individual functional connections generally covary only very weakly with head motion estimates, motion influence is difficult to quantify robustly, and prone to be neglected in practice. Although the use of individual estimates of head motion, or group-level correlation of motion and functional connectivity has been suggested, a sufficiently sensitive measure of individual functional connectivity quality has not yet been established. We propose a new intuitive summary index, Typicality of Functional Connectivity, to capture deviations from standard brain functional connectivity patterns. In a resting-state fMRI dataset of 245 healthy subjects, this measure was significantly correlated with individual head motion metrics. The results were further robustly reproduced across atlas granularity, preprocessing options, and other datasets, including 1,081 subjects from the Human Connectome Project. In principle, Typicality of Functional Connectivity should be sensitive also to other types of artifacts, processing errors, and possibly also brain pathology, allowing extensive use in data quality screening and quantification in functional connectivity studies as well as methodological investigations.

Functional connectivity of EEG is subject-specific, associated with phenotype, and different from fMRI

A variety of psychiatric, behavioral and cognitive phenotypes have been linked to brain ''functional connectivity'' -- the pattern of correlation observed between different brain regions. Most commonly assessed using functional magnetic resonance imaging (fMRI), here, we investigate the connectivity-phenotype associations with functional connectivity measured with electroencephalography (EEG), using phase-coupling. We analyzed data from the publicly available Healthy Brain Network Biobank. This database compiles a growing sample of children and adolescents, currently encompassing 1657 individuals. Among a variety of assessment instruments we focus on ten phenotypic and additional demographic measures that capture most of the variance in this sample. The largest effect sizes are found for age and sex for both fMRI and EEG. We replicate previous findings of an association of Intelligence Quotient (IQ) and Attention Deficit Hyperactivity Disorder (ADHD) with the pattern of fMRI functional connectivity. We also find an association with socioeconomic status, anxiety and the Child Behavior Checklist Score. For EEG we find a significant connectivity-phenotype relationship with IQ. The actual spatial patterns of functional connectivity are quite different between fMRI and source-space EEG. However, within EEG we observe clusters of functional connectivity that are consistent across frequency bands. Additionally we analyzed reproducibility of functional connectivity. We compare connectivity obtained with different tasks, including resting state, a video and a visual flicker task. For both EEG and fMRI the variation between tasks was smaller than the variability observed between subjects. We also found an increase of reliability with increasing frequency of the EEG, and increased sampling duration. We conclude that, while the patterns of functional connectivity are distinct between fMRI and phase-coupling of EEG, they are nonetheless similar in their robustness to the task, and similar in that idiosyncratic patterns of connectivity predict individual phenotypes.

Domino-like transient dynamics at seizure onset in epilepsy

The International League Against Epilepsy (ILAE) groups seizures into "focal", "generalized" and "unknown" based on whether the seizure onset is confined to a brain region in one hemisphere, arises in several brain region simultaneously, or is not known, respectively. This separation fails to account for the rich diversity of clinically and experimentally observed spatiotemporal patterns of seizure onset and even less so for the properties of the brain networks generating them. We consider three different patterns of domino-like seizure onset in Idiopathic Generalized Epilepsy (IGE) and present a novel approach to classification of seizures. To understand how these patterns are generated on networks requires understanding of the relationship between intrinsic node dynamics and coupling between nodes in the presence of noise, which currently is unknown. We investigate this interplay here in the framework of domino-like recruitment across a network. In particular, we use a phenomenological model of seizure onset with heterogeneous coupling and node properties, and show that in combination they generate a range of domino-like onset patterns observed in the IGE seizures. We further explore the individual contribution of heterogeneous node dynamics and coupling by interpreting in-vitro experimental data in which the speed of onset can be chemically modulated. This work contributes to a better understanding of possible drivers for the spatiotemporal patterns observed at seizure onset and may ultimately contribute to a more personalized approach to classification of seizure types in clinical practice.

Idiosynchrony: From shared responses to individual differences during naturalistic neuroimaging

Two ongoing movements in human cognitive neuroscience have researchers shifting focus from group-level inferences to characterizing single subjects, and complementing tightly controlled tasks with rich, dynamic paradigms such as movies and stories. Yet relatively little work combines these two, perhaps because traditional analysis approaches for naturalistic imaging data are geared toward detecting shared responses rather than between-subject variability. Here, we review recent work using naturalistic stimuli to study individual differences, and advance a framework for detecting structure in idiosyncratic patterns of brain activity, or "idiosynchrony". Specifically, we outline the emerging technique of inter-subject representational similarity analysis (IS-RSA), including its theoretical motivation and an empirical demonstration of how it recovers brain-behavior relationships during movie watching using data from the Human Connectome Project. We also consider how stimulus choice may affect the individual signal and discuss areas for future research. We argue that naturalistic neuroimaging paradigms have the potential to reveal meaningful individual differences above and beyond those observed during traditional tasks or at rest.

Uncovering differential identifiability in network properties of human brain functional connectomes

The identifiability framework (𝕀f) has been shown to improve differential identifiability (reliability across-sessions and -sites, and differentiability across-subjects) of functional connectomes for a variety of fMRI tasks. But having a robust single session/subject functional connectome is just the starting point to subsequently assess network properties for characterizing properties of integration, segregation, and communicability, among others. Naturally, one wonders whether uncovering identifiability at the connectome level also uncovers identifiability on the derived network properties. This also raises the question of where to apply the 𝕀f framework: on the connectivity data or directly on each network measurement? Our work answers these questions by exploring the differential identifiability profiles of network measures when 𝕀f is applied (a) on the functional connectomes, and (b) directly on derived network measurements. Results show that improving across-session reliability of functional connectomes (FCs) also improves reliability of derived network measures. We also find that, for specific network properties, application of 𝕀f directly on network properties is more effective. Finally, we discover that applying the framework, either way, increases task sensitivity of network properties. At a time when the neuroscientific community is focused on subject-level inferences, this framework is able to uncover FC fingerprints, which propagate to derived network properties.