High-accuracy individual identification using a "thin slice" of the functional connectome

Using independent component analysis to extract a cross-modality and individual-specific brain baseline pattern

The ongoing brain activity serves as a baseline that supports both internal and external cognitive processes. However, its precise nature remains unclear. Considering that people display various patterns of brain activity even when engaging in the same task, it is reasonable to believe that individuals possess their unique brain baseline pattern. Using spatial independent component analysis on a large sample of fMRI data from the Human Connectome Project (HCP), we found an individual-specific component which can be consistently extracted from either resting-state or different task states and is reliable over months. Compared to functional connectome fingerprinting, it is much more stable across different fMRI modalities. Its stability is closely related to high explained variance and is minimally influenced by factors such as noise, scan duration, and scan interval. We propose that this component underlying the ongoing activity represents an individual-specific baseline pattern of brain activity.

Multimodal gradients of basal forebrain connectivity across the neocortex

Cortical cholinergic projections originate from subregions of the basal forebrain (BF). To examine its organization in humans, we computed multimodal gradients of BF connectivity by combining 7 T diffusion and resting state functional MRI. Moving from anteromedial to posterolateral BF, we observe reduced tethering between structural and functional connectivity gradients, with the lowest tethering in the nucleus basalis of Meynert. In the neocortex, this gradient is expressed by progressively reduced tethering from unimodal sensory to transmodal cortex, with the lowest tethering in the midcingulo-insular network, and is also spatially correlated with the molecular concentration of VAChT, measured by [18F]fluoroethoxy-benzovesamicol (FEOBV) PET. In mice, viral tracing of BF cholinergic projections and [18F]FEOBV PET confirm a gradient of axonal arborization. Altogether, our findings reveal that BF cholinergic neurons vary in their branch complexity, with certain subpopulations exhibiting greater modularity and others greater diffusivity in the functional integration with their cortical targets.

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.

Integrative, segregative, and degenerate harmonics of the structural connectome

Unifying integration and segregation in the brain has been a fundamental puzzle in neuroscience ever since the conception of the "binding problem." Here, we introduce a framework that places integration and segregation within a continuum based on a fundamental property of the brain-its structural connectivity graph Laplacian harmonics and a new feature we term the gap-spectrum. This framework organizes harmonics into three regimes-integrative, segregative, and degenerate-that together account for various group-level properties. Integrative and segregative harmonics occupy the ends of the continuum, and they share properties such as reproducibility across individuals, stability to perturbation, and involve "bottom-up" sensory networks. Degenerate harmonics are in the middle of the continuum, and they are subject-specific, flexible, and involve "top-down" networks. The proposed framework accommodates inter-subject variation, sensitivity to changes, and structure-function coupling in ways that offer promising avenues for studying cognition and consciousness in the brain.

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.

Multimodal gradients of basal forebrain connectivity across the neocortex

The cholinergic innervation of the cortex originates almost entirely from populations of neurons in the basal forebrain (BF). Structurally, the ascending BF cholinergic projections are highly branched, with individual cells targeting multiple different cortical regions. However, it is not known whether the structural organization of basal forebrain projections reflects their functional integration with the cortex. We therefore used high-resolution 7T diffusion and resting state functional MRI in humans to examine multimodal gradients of BF cholinergic connectivity with the cortex. Moving from anteromedial to posterolateral BF, we observed reduced tethering between structural and functional connectivity gradients, with the most pronounced dissimilarity localized in the nucleus basalis of Meynert (NbM). The cortical expression of this structure-function gradient revealed progressively weaker tethering moving from unimodal to transmodal cortex, with the lowest tethering in midcingulo-insular cortex. We used human [ 18 F] fluoroethoxy-benzovesamicol (FEOBV) PET to demonstrate that cortical areas with higher concentrations of cholinergic innervation tend to exhibit lower tethering between BF structural and functional connectivity, suggesting a pattern of increasingly diffuse axonal arborization. Anterograde viral tracing of cholinergic projections and [ 18 F] FEOBV PET in mice confirmed a gradient of axonal arborization across individual BF cholinergic neurons. Like humans, cholinergic neurons with the highest arborization project to cingulo-insular areas of the mouse isocortex. Altogether, our findings reveal that BF cholinergic neurons vary in their branch complexity, with certain subpopulations exhibiting greater modularity and others greater diffusivity in the functional integration of their cortical targets.

Single-subject cortical morphological brain networks: Phenotypic associations and neurobiological substrates

Although single-subject morphological brain networks provide an important way for human connectome studies, their roles and origins are poorly understood. Combining cross-sectional and repeated structural magnetic resonance imaging scans from adults, children and twins with behavioral and cognitive measures and brain-wide transcriptomic, cytoarchitectonic and chemoarchitectonic data, this study examined phenotypic associations and neurobiological substrates of single-subject morphological brain networks. We found that single-subject morphological brain networks explained inter-individual variance and predicted individual outcomes in Motor and Cognition domains, and distinguished individuals from each other. The performance can be further improved by integrating different morphological indices for network construction. Low-moderate heritability was observed for single-subject morphological brain networks with the highest heritability for sulcal depth-derived networks and higher heritability for inter-module connections. Furthermore, differential roles of genetic, cytoarchitectonic and chemoarchitectonic factors were observed for single-subject morphological brain networks. Cortical thickness-derived networks were related to the three factors with contributions from genes enriched in membrane and transport related functions, genes preferentially located in supragranular and granular layers, overall thickness in the molecular layer and thickness of wall in the infragranular layers, and metabotropic glutamate receptor 5 and dopamine transporter; fractal dimension-, gyrification index- and sulcal depth-derived networks were only associated with the chemoarchitectonic factor with contributions from different sets of neurotransmitter receptors. Most results were reproducible across different parcellation schemes and datasets. Altogether, this study demonstrates phenotypic associations and neurobiological substrates of single-subject morphological brain networks, which provide intermediate endophenotypes to link molecular and cellular architecture and behavior and cognition.

Myelination and excitation-inhibition balance synergistically shape structure-function coupling across the human cortex

Recent work has demonstrated that the relationship between structural and functional connectivity varies regionally across the human brain, with reduced coupling emerging along the sensory-association cortical hierarchy. The biological underpinnings driving this expression, however, remain largely unknown. Here, we postulate that intracortical myelination and excitation-inhibition (EI) balance mediate the heterogeneous expression of structure-function coupling (SFC) and its temporal variance across the cortical hierarchy. We employ atlas- and voxel-based connectivity approaches to analyze neuroimaging data acquired from two groups of healthy participants. Our findings are consistent across six complementary processing pipelines: 1) SFC and its temporal variance respectively decrease and increase across the unimodal-transmodal and granular-agranular gradients; 2) increased myelination and lower EI-ratio are associated with more rigid SFC and restricted moment-to-moment SFC fluctuations; 3) a gradual shift from EI-ratio to myelination as the principal predictor of SFC occurs when traversing from granular to agranular cortical regions. Collectively, our work delivers a framework to conceptualize structure-function relationships in the human brain, paving the way for an improved understanding of how demyelination and/or EI-imbalances induce reorganization in brain disorders.

Few temporally distributed brain connectivity states predict human cognitive abilities

Human functional brain connectivity can be temporally decomposed into states of high and low cofluctuation, defined as coactivation of brain regions over time. Rare states of particularly high cofluctuation have been shown to reflect fundamentals of intrinsic functional network architecture and to be highly subject-specific. However, it is unclear whether such network-defining states also contribute to individual variations in cognitive abilities - which strongly rely on the interactions among distributed brain regions. By introducing CMEP, a new eigenvector-based prediction framework, we show that as few as 16 temporally separated time frames (< 1.5% of 10min resting-state fMRI) can significantly predict individual differences in intelligence (N = 263, p < .001). Against previous expectations, individual's network-defining time frames of particularly high cofluctuation do not predict intelligence. Multiple functional brain networks contribute to the prediction, and all results replicate in an independent sample (N = 831). Our results suggest that although fundamentals of person-specific functional connectomes can be derived from few time frames of highest connectivity, temporally distributed information is necessary to extract information about cognitive abilities. This information is not restricted to specific connectivity states, like network-defining high-cofluctuation states, but rather reflected across the entire length of the brain connectivity time series.

Uncovering individual differences in fine-scale dynamics of functional connectivity

Functional connectivity (FC) profiles contain subject-specific features that are conserved across time and have potential to capture brain-behavior relationships. Most prior work has focused on spatial features (nodes and systems) of these FC fingerprints, computed over entire imaging sessions. We propose a method for temporally filtering FC, which allows selecting specific moments in time while also maintaining the spatial pattern of node-based activity. To this end, we leverage a recently proposed decomposition of FC into edge time series (eTS). We systematically analyze functional magnetic resonance imaging frames to define features that enhance identifiability across multiple fingerprinting metrics, similarity metrics, and data sets. Results show that these metrics characteristically vary with eTS cofluctuation amplitude, similarity of frames within a run, transition velocity, and expression of functional systems. We further show that data-driven optimization of features that maximize fingerprinting metrics isolates multiple spatial patterns of system expression at specific moments in time. Selecting just 10% of the data can yield stronger fingerprints than are obtained from the full data set. Our findings support the idea that FC fingerprints are differentially expressed across time and suggest that multiple distinct fingerprints can be identified when spatial and temporal characteristics are considered simultaneously.

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.

Functional connectomes become more longitudinally self-stable, but not more distinct from others, across early childhood

Functional connectomes, as measured with functional magnetic resonance imaging (fMRI), are highly individualized, and evidence suggests this individualization may increase across childhood. A connectome can become more individualized either by increasing self-stability or decreasing between-subject-similarity. Here we used a longitudinal early childhood dataset to investigate age associations with connectome self-stability, between-subject-similarity, and developmental individualization, defined as an individual's self-stability across a 12-month interval relative to their between-subject-similarity. fMRI data were collected during an 18-minute passive viewing scan from 73 typically developing children aged 4-7 years, at baseline and 12-month follow-up. We found that young children had highly individualized connectomes, with sufficient self-stability across 12-months for 98% identification accuracy. Linear models showed a significant relationship between age and developmental individualization across the whole brain and in most networks. This association appeared to be largely driven by an increase in self-stability with age, with only weak evidence for relationships between age and similarity across participants. Together our findings suggest that children's connectomes become more individualized across early childhood, and that this effect is driven by increasing self-stability rather than decreasing between-subject-similarity.

Sharing individualised template MRI data for MEG source reconstruction: A solution for open data while keeping subject confidentiality

The increasing requirements for adoption of FAIR data management and sharing original research data from neuroimaging studies can be at odds with protecting the anonymity of the research participants due to the person-identifiable anatomical features in the data. We propose a solution to this dilemma for anatomical MRIs used in MEG source analysis. In MEG analysis, the channel-level data is reconstructed to the source-level using models derived from anatomical MRIs. Sharing data, therefore, requires sharing the anatomical MRI to replicate the analysis. The suggested solution is to replace the individual anatomical MRIs with individualised warped templates that can be used to carry out the MEG source analysis and that provide sufficient geometrical similarity to the original participants' MRIs. First, we demonstrate how the individualised template warping can be implemented with one of the leading open-source neuroimaging analysis toolboxes. Second, we compare results from four different MEG source reconstruction methods performed with an individualised warped template to those using the participant's original MRI. While the source reconstruction results are not numerically identical, there is a high similarity between the results for single dipole fits, dynamic imaging of coherent sources beamforming, and atlas-based virtual channel beamforming. There is a moderate similarity between minimum-norm estimates, as anticipated due to this method being anatomically constrained and dependent on the exact morphological features of the cortical sheet. We also compared the morphological features of the warped template to those of the original MRI. These showed a high similarity in grey matter volume and surface area, but a low similarity in the average cortical thickness and the mean folding index within cortical parcels. Taken together, this demonstrates that the results obtained by MEG source reconstruction can be preserved with the warped templates, whereas the anatomical and morphological fingerprint is sufficiently altered to protect the anonymity of research participants. In cases where participants consent to sharing anatomical MRI data, it remains preferable to share the original defaced data with an appropriate data use agreement. In cases where participants did not consent to share their MRIs, the individualised warped MRI template offers a good compromise in sharing data for reuse while retaining anonymity for research participants.

Benchmarking common preprocessing strategies in early childhood functional connectivity and intersubject correlation fMRI

Preprocessing choices present a particular challenge for researchers working with functional magnetic resonance imaging (fMRI) data from young children. Steps which have been shown to be important for mitigating head motion, such as censoring and global signal regression (GSR), remain controversial, and benchmarking studies comparing preprocessing pipelines have been conducted using resting data from older participants who tend to move less than young children. Here, we conducted benchmarking of fMRI preprocessing steps in a population with high head-motion, children aged 4-8 years, leveraging a unique longitudinal, passive viewing fMRI dataset. We systematically investigated combinations of global signal regression (GSR), volume censoring, and ICA-AROMA. Pipelines were compared using previously established metrics of noise removal as well as metrics sensitive to recovery of individual differences (i.e., connectome fingerprinting), and stimulus-evoked responses (i.e., intersubject correlations; ISC). We found that: 1) the most efficacious pipeline for both noise removal and information recovery included censoring, GSR, bandpass filtering, and head motion parameter (HMP) regression, 2) ICA-AROMA performed similarly to HMP regression and did not obviate the need for censoring, 3) GSR had a minimal impact on connectome fingerprinting but improved ISC, and 4) the strictest censoring approaches reduced motion correlated edges but negatively impacted identifiability.

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.

Existence of Functional Connectome Fingerprint during Infancy and Its Stability over Months

The functional connectome fingerprint is a cluster of individualized brain functional connectivity patterns that are capable of distinguishing one individual from others. Although its existence has been demonstrated in adolescents and adults, whether such individualized patterns exist during infancy is barely investigated despite its importance in identifying the origin of the intrinsic connectome patterns that potentially mirror distinct behavioral phenotypes. To fill this knowledge gap, capitalizing on a longitudinal high-resolution structural and resting-state functional MRI dataset with 104 human infants (53 females) with 806 longitudinal scans (age, 16-876 d) and infant-specific functional parcellation maps, we observe that the brain functional connectome fingerprint may exist since infancy and keeps stable over months during early brain development. Specifically, we achieve an ∼78% individual identification rate by using ∼5% selected functional connections, compared with the best identification rate of 60% without connection selection. The frontoparietal networks recognized as the most contributive networks in adult functional connectome fingerprinting retain their superiority in infants despite being widely acknowledged as rapidly developing systems during childhood. The existence and stability of the functional connectome fingerprint are further validated on adjacent age groups. Moreover, we show that the infant frontoparietal networks can reach similar accuracy in predicting individual early learning composite scores as the whole-brain connectome, again resembling the observations in adults and highlighting the relevance of functional connectome fingerprint to cognitive performance. For the first time, these results suggest that each individual may retain a unique and stable marker of functional connectome during early brain development.SIGNIFICANCE STATEMENT Functional connectome fingerprinting during infancy featuring rapid brain development remains almost uninvestigated even though it is essential for understanding the early individual-level intrinsic pattern of functional organization and its relationship with distinct behavioral phenotypes. With an infant-tailored functional connection selection and validation strategy, we strive to provide the delineation of the infant functional connectome fingerprint by examining its existence, stability, and relationship with early cognitive performance. We observe that the brain functional connectome fingerprint may exist since early infancy and remains stable over months during the first 2 years. The identified key contributive functional connections and networks for fingerprinting are also verified to be highly predictive for cognitive score prediction, which reveals the association between infant connectome fingerprint and cognitive performance.

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.

Inferring the heritability of large-scale functional networks with a multivariate ACE modeling approach

Recent evidence suggests that the human functional connectome is stable at different timescales and is unique. These characteristics posit the functional connectome not only as an individual marker but also as a powerful discriminatory measure characterized by high intersubject variability. Among distinct sources of intersubject variability, the long-term sources include functional patterns that emerge from genetic factors. Here, we sought to investigate the contribution of additive genetic factors to the variability of functional networks by determining the heritability of the connectivity strength in a multivariate fashion. First, we reproduced and extended the connectome fingerprinting analysis to the identification of twin pairs. Then, we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. Twin pairs were identified above chance level using connectome fingerprinting, with monozygotic twin identification accuracy equal to 57.2% on average for whole-brain connectome. Additionally, we found that a visual (0.37), the medial frontal (0.31), and the motor (0.30) functional networks were the most influenced by additive genetic factors. Our findings suggest that genetic factors not only partially determine intersubject variability of the functional connectome, such that twins can be identified using connectome fingerprinting, but also differentially influence connectivity strength in large-scale functional networks.

The functional connectome is a unique representation of the functional organization of the human brain. As such, it has been extensively used as an individual marker, a “fingerprint,” because of its high intersubject variability. Here, we sought to investigate the influence of genetic factors on intersubject variability of functional networks. Therefore, we extended the connectome fingerprinting analysis to the identification of twin pairs, and we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. We found that genetic factors not only partially determine intersubject variability of the functional connectome, such that monozygotic twin identification accuracy achieved 57.2% on average using whole-brain connectome in the fingerprinting analysis, but also differentially influence connectivity strength in large-scale functional networks.

Beyond fingerprinting: Choosing predictive connectomes over reliable connectomes

Recent years have seen a surge of research on variability in functional brain connectivity within and between individuals, with encouraging progress toward understanding the consequences of this variability for cognition and behavior. At the same time, well-founded concerns over rigor and reproducibility in psychology and neuroscience have led many to question whether functional connectivity is sufficiently reliable, and call for methods to improve its reliability. The thesis of this opinion piece is that when studying variability in functional connectivity-both across individuals and within individuals over time-we should use behavior prediction as our benchmark rather than optimize reliability for its own sake. We discuss theoretical and empirical evidence to compel this perspective, both when the goal is to study stable, trait-level differences between people, as well as when the goal is to study state-related changes within individuals. We hope that this piece will be useful to the neuroimaging community as we continue efforts to characterize inter- and intra-subject variability in brain function and build predictive models with an eye toward eventual real-world applications.

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.

Dynamic expression of brain functional systems disclosed by fine-scale analysis of edge time series

Functional connectivity (FC) describes the statistical dependence between neuronal populations or brain regions in resting-state fMRI studies and is commonly estimated as the Pearson correlation of time courses. Clustering or community detection reveals densely coupled sets of regions constituting resting-state networks or functional systems. These systems manifest most clearly when FC is sampled over longer epochs but appear to fluctuate on shorter timescales. Here, we propose a new approach to reveal temporal fluctuations in neuronal time series. Unwrapping FC signal correlations yields pairwise co-fluctuation time series, one for each node pair or edge, and allows tracking of fine-scale dynamics across the network. Co-fluctuations partition the network, at each time step, into exactly two communities. Sampled over time, the overlay of these bipartitions, a binary decomposition of the original time series, very closely approximates functional connectivity. Bipartitions exhibit characteristic spatiotemporal patterns that are reproducible across participants and imaging runs, capture individual differences, and disclose fine-scale temporal expression of functional systems. Our findings document that functional systems appear transiently and intermittently, and that FC results from the overlay of many variable instances of system expression. Potential applications of this decomposition of functional connectivity into a set of binary patterns are discussed.

Constructing Compact Signatures for Individual Fingerprinting of Brain Connectomes

Recent neuroimaging studies have shown that functional connectomes are unique to individuals, i.e., two distinct fMRIs taken over different sessions of the same subject are more similar in terms of their connectomes than those from two different subjects. In this study, we present new results that identify specific parts of resting state and task-specific connectomes that are responsible for the unique signatures. We show that a very small part of the connectome can be used to derive features for discriminating between individuals. A network of these features is shown to achieve excellent training and test accuracy in matching imaging datasets. We show that these features are statistically significant, robust to perturbations, invariant across populations, and are localized to a small number of structural regions of the brain. Furthermore, we show that for task-specific connectomes, the regions identified by our method are consistent with their known functional characterization. We present a new matrix sampling technique to derive computationally efficient and accurate methods for identifying the discriminating sub-connectome and support all of our claims using state-of-the-art statistical tests and computational techniques.

Jumping over baselines with new methods to predict activation maps from resting-state fMRI

Cognitive fMRI research primarily relies on task-averaged responses over many subjects to describe general principles of brain function. Nonetheless, there exists a large variability between subjects that is also reflected in spontaneous brain activity as measured by resting state fMRI (rsfMRI). Leveraging this fact, several recent studies have therefore aimed at predicting task activation from rsfMRI using various machine learning methods within a growing literature on 'connectome fingerprinting'. In reviewing these results, we found lack of an evaluation against robust baselines that reliably supports a novelty of predictions for this task. On closer examination to reported methods, we found most underperform against trivial baseline model performances based on massive group averaging when whole-cortex prediction is considered. Here we present a modification to published methods that remedies this problem to large extent. Our proposed modification is based on a single-vertex approach that replaces commonly used brain parcellations. We further provide a summary of this model evaluation by characterizing empirical properties of where prediction for this task appears possible, explaining why some predictions largely fail for certain targets. Finally, with these empirical observations we investigate whether individual prediction scores explain individual behavioral differences in a task.

High-amplitude cofluctuations in cortical activity drive functional connectivity

Resting-state functional connectivity is used throughout neuroscience to study brain organization and to generate biomarkers of development, disease, and cognition. The processes that give rise to correlated activity are, however, poorly understood. Here we decompose resting-state functional connectivity using a temporal unwrapping procedure to assess the contributions of moment-to-moment activity cofluctuations to the overall connectivity pattern. This approach temporally resolves functional connectivity at a timescale of single frames, which enables us to make direct comparisons of cofluctuations of network organization with fluctuations in the blood oxygen level-dependent (BOLD) time series. We show that surprisingly, only a small fraction of frames exhibiting the strongest cofluctuation amplitude are required to explain a significant fraction of variance in the overall pattern of connection weights as well as the network's modular structure. These frames coincide with frames of high BOLD activity amplitude, corresponding to activity patterns that are remarkably consistent across individuals and identify fluctuations in default mode and control network activity as the primary driver of resting-state functional connectivity. Finally, we demonstrate that cofluctuation amplitude synchronizes across subjects during movie watching and that high-amplitude frames carry detailed information about individual subjects (whereas low-amplitude frames carry little). Our approach reveals fine-scale temporal structure of resting-state functional connectivity and discloses that frame-wise contributions vary across time. These observations illuminate the relation of brain activity to functional connectivity and open a number of directions for future research.

Disentangled Intensive Triplet Autoencoder for Infant Functional Connectome Fingerprinting

Functional connectome "fingerprint" is a highly characterized brain pattern that distinguishes one individual from others. Although its existence has been demonstrated in adults, an unanswered but fundamental question is whether such individualized pattern emerges since infancy. This problem is barely investigated despites its importance in identifying the origin of the intrinsic connectome patterns that mirror distinct behavioral phenotypes. However, addressing this knowledge gap is challenging because the conventional methods are only applicable to developed brains with subtle longitudinal changes and typically fail on the dramatically developing infant brains. To tackle this challenge, we invent a novel model, namely, disentangled intensive triplet autoencoder (DI-TAE). First, we introduce the triplet autoencoder to embed the original connectivity into a latent space with higher discriminative capability among infant individuals. Then, a disentanglement strategy is proposed to separate the latent variables into identity-code, age-code, and noise-code, which not only restrains the interference from age-related developmental variance, but also captures the identity-related invariance. Next, a cross-reconstruction loss and an intensive triplet loss are designed to guarantee the effectiveness of the disentanglement and enhance the inter-subject dissimilarity for better discrimination. Finally, a variance-guided bootstrap aggregating is developed for DI-TAE to further improve the performance of identification. DI-TAE is validated on three longitudinal resting-state fMRI datasets with 394 infant scans aged 16 to 874 days. Our proposed model outperforms other state-of-the-art methods by increasing the identification rate by more than 50%, and for the first time suggests the plausible existence of brain functional connectome "fingerprint" since early infancy.

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.

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.

Accurate prediction of individual subject identity and task, but not autism diagnosis, from functional connectomes

Despite enthusiasm about the potential for using fMRI-based functional connectomes in the development of biomarkers for autism spectrum disorder (ASD), the literature is full of negative findings-failures to distinguish ASD functional connectomes from those of typically developing controls (TD)-and positive findings that are inconsistent across studies. Here, we report on a new study designed to either better differentiate ASD from TD functional connectomes-or, alternatively, to refine our understanding of the factors underlying the current state of affairs. We scanned individuals with ASD and controls both at rest and while watching videos with social content. Using multiband fMRI across repeat sessions, we improved both data quantity and scanning duration by collecting up to 2 hr of data per individual. This is about 50 times the typical number of temporal samples per individual in ASD fcMRI studies. We obtained functional connectomes that were discriminable, allowing for near-perfect individual identification regardless of diagnosis, and equally reliable in both groups. However, contrary to what one might expect, we did not consistently or robustly observe in the ASD group either reductions in similarity to TD functional connectivity (FC) patterns or shared atypical FC patterns. Accordingly, FC-based predictions of diagnosis group achieved accuracy levels around chance. However, using the same approaches to predict scan type (rest vs. video) achieved near-perfect accuracy. Our findings suggest that neither the limitations of resting state as a "task," data resolution, data quantity, or scan duration can be considered solely responsible for failures to differentiate ASD from TD functional connectomes.

Functional Connectivity Fingerprints at Rest Are Similar across Youths and Adults and Vary with Genetic Similarity

Distinguishing individuals from brain connectivity, and studying the genetic influences on that identification across different ages, improves our basic understanding of functional brain network organization. We applied support vector machine classifiers to two datasets of twins (adult, pediatric) and two datasets of repeat-scan individuals (adult, pediatric). Classifiers were trained on resting state functional connectivity magnetic resonance imaging (rs-fcMRI) data and used to predict individuals and co-twin pairs from independent data. The classifiers successfully identified individuals from a previous scan with 100% accuracy, even when scans were separated by months. In twin samples, classifier accuracy decreased as genetic similarity decreased. Our results demonstrate that classification is stable within individuals, similar within families, and contains similar representations of functional connections over a few decades of life. Moreover, the degree to which these patterns of connections predict siblings' data varied by genetic relatedness, suggesting that genetic influences on rs-fcMRI connectivity are established early in life.

Nonreplication of functional connectivity differences in autism spectrum disorder across multiple sites and denoising strategies

A rapidly growing number of studies on autism spectrum disorder (ASD) have used resting‐state fMRI to identify alterations of functional connectivity, with the hope of identifying clinical biomarkers or underlying neural mechanisms. However, results have been largely inconsistent across studies, and there remains a pressing need to determine the primary factors influencing replicability. Here, we used resting‐state fMRI data from the Autism Brain Imaging Data Exchange to investigate two potential factors: denoising strategy and data site (which differ in terms of sample, data acquisition, etc.). We examined the similarity of both group‐averaged functional connectomes and group‐level differences (ASD vs. control) across 33 denoising pipelines and four independently‐acquired datasets. The group‐averaged connectomes were highly consistent across pipelines (r = 0.92 ± 0.06) and sites (r = 0.88 ± 0.02). However, the group differences, while still consistent within site across pipelines (r = 0.76 ± 0.12), were highly inconsistent across sites regardless of choice of denoising strategies (r = 0.07 ± 0.04), suggesting lack of replication may be strongly influenced by site and/or cohort differences. Across‐site similarity remained low even when considering the data at a large‐scale network level or when considering only the most significant edges. We further show through an extensive literature survey that the parameters chosen in the current study (i.e., sample size, age range, preprocessing methods) are quite representative of the published literature. These results highlight the importance of examining replicability in future studies of ASD, and, more generally, call for extra caution when interpreting alterations in functional connectivity across groups of individuals.