Search

Association fibers connect the cortical regions and experience rapid development involving myelination and axonal growth during infancy. Yet, the spatiotemporal patterns of microstructural changes along these tracts, as well as the developmental interaction between the white matter (WM) tracts and the cortical gray matter (cGM) connected to them, are mostly unknown during infancy. In this study, we performed a diffusion MRI-based tractography and microstructure study in a cohort of 89 healthy preterm-born infants with gestational age at birth between 28.1∼36.4 weeks and postmenstrual age at scan between 39.9∼59.9 weeks. Results revealed that several C-shaped fibers, such as the arcuate fasciculus, cingulum, and uncinate fasciculus, demonstrated symmetrical along-tract profiles; and the horizontally oriented running fibers, including the inferior fronto-occipital fasciculus and the inferior longitudinal fasciculus, demonstrated an anterior-posterior developmental gradient. This study characterized the along-tract profiles using fixel-based analysis and revealed that the fiber cross-section (FC) of all five association fibers demonstrated a fluctuating increase with age, while the fiber density (FD) monotonically increase with age. NODDI was utilized to analyze the microstructural development of cGM and indicated cGM connected to the anterior end of the association fibers developed faster than that of the posterior end during 0-5 months. Notably, a mediation analysis was used to explore the relation between the development of WM and associated cGM, and demonstrated a partial mediation effect of FD in WM on the development of intracellular volume (ICV) in cGM and a full mediation effect of ICV on the growth of FD in most fibers, suggesting a predominant mediation of cGM on the WM development. Furthermore, for assessing whether those results were biased by prematurity, we compared preterm- and term-born neonates with matched scan age, gender, and multiple births from the developing human connectome project (dHCP) dataset to assess the effect of preterm-birth, and the results indicated a similar developmental pattern of the association fibers and their attached cGM. These findings presented a comprehensive picture of the major association fibers during early infancy and deciphered the developmental interaction between WM and cGM in this period.

The 'core' regions of the default mode network (DMN) - the medial prefrontal cortex (MPFC), the posterior cingulate cortex (PCC), and inferior parietal lobules (IPL) - show consistent engagement across mental states that involve self-oriented processing. Precisely how these regions interact in support of such processes remains an important unanswered question. In the current functional magnetic resonance imaging (fMRI) study, we examined dynamic interactions of the 'core-self' DMN regions during two forms of self-referential cognition: direct self-appraisal (thinking about oneself) and reflected self-appraisal (thinking about oneself from a third-person perspective). One-hundred and eleven participants completed our dual self-appraisal task during fMRI, and general linear models were used to characterize common and distinct neural responses to these conditions. Informed by these results, we then applied dynamic causal modelling to examine causal interactions among the 'core-self' regions, and how they were specifically modulated under the influence of direct and reflected self-appraisal. As a primary observation, this network modelling revealed a distinct inhibitory influence of the left IPL on the PCC during reflected compared to direct self-appraisal, which was accompanied by evidence of greater activation in both regions during the reflected self-appraisal condition. We suggest that the greater engagement of posterior DMN regions during reflected self-appraisal is a function of the higher-order processing needed for this form of self-appraisal, with the left IPL supporting abstract self-related processes including episodic memory retrieval and shifts of perspective. Overall, we show that core DMN regions interact in functionally unique ways in support of self-referential processes, even when these processes are inter-related. Further characterization of DMN functional interactions across self-related mental states is likely to inform a deeper understanding of how this brain network orchestrates the self.

Each human brain has a unique functional synchronisation pattern (functional connectome) analogous to a fingerprint that underpins brain functions and related behaviours. Here we examine functional connectome (whole-brain and 13 networks) maturation by measuring its uniqueness in adolescents who underwent brain scans longitudinally from 12 years of age every four months. The uniqueness of a functional connectome is defined as its ratio of self-similarity (from the same subject at a different time point) to the maximal similarity-to-others (from a given subject and any others at a different time point). We found that the unique whole brain connectome exists in 12 years old adolescents, with 92% individuals having a whole brain uniqueness value greater than one. The cingulo-opercular network (CON; a long-acting 'brain control network' configuring information processing) demonstrated marginal uniqueness in early adolescence with 56% of individuals showing uniqueness greater than one (i.e., more similar to her/his own CON four months later than those from any other subjects) and this increased longitudinally. Notably, the low uniqueness of the CON correlates (β = -18.6, FDR-Q 0.001) with K10 levels at the subsequent time point. This association suggests that the individualisation of CON network is related to psychological distress levels. Our findings highlight the potential of longitudinal neuroimaging to capture mental health problems in young people who are undergoing profound neuroplasticity and environment sensitivity period.

Structural and functional brain networks are modular. Canonical functional systems, such as the default mode network, are well-known modules of the human brain and have been implicated in a large number of cognitive, behavioral and clinical processes. However, modules delineated in structural brain networks inferred from tractography generally do not recapitulate canonical functional systems. Neuroimaging evidence suggests that functional connectivity between regions in the same systems is not always underpinned by anatomical connections. As such, direct structural connectivity alone would be insufficient to characterize the functional modular organization of the brain. Here, we demonstrate that augmenting structural brain networks with models of indirect (polysynaptic) communication unveils a modular network architecture that more closely resembles the brain’s established functional systems. We find that diffusion models of polysynaptic connectivity, particularly communicability, narrow the gap between the modular organization of structural and functional brain networks by 20–60%, whereas routing models based on single efficient paths do not improve mesoscopic structure-function correspondence. This suggests that functional modules emerge from the constraints imposed by local network structure that facilitates diffusive neural communication. Our work establishes the importance of modeling polysynaptic communication to understand the structural basis of functional systems.

Mapping the human connectome and understanding its relationship to brain function holds tremendous clinical potential. The connectome has two fundamental components: the nodes and the sconnections between them. While much attention has been given to deriving atlases and measuring the connections between nodes, there have been no studies examining the networks within nodes. Here we demonstrate that each node contains significant connectivity information, that varies systematically across task-induced states and subjects, such that measures based on these variations can be used to classify tasks and identify subjects. The results are not specific for any particular atlas but hold across different atlas resolutions. To date, studies examining changes in connectivity have focused on edge changes and assumed there is no useful information within nodes. Our findings illustrate that for typical atlases, within-node changes can be significant and may account for a substantial fraction of the variance currently attributed to edge changes .

Control processes allow us to constrain the retrieval of semantic information from long-term memory so that it is appropriate for the task or context. Control demands are influenced by the strength of the target information itself and by the circumstances in which it is retrieved, with more control needed when relatively weak aspects of knowledge are required and after the sustained retrieval of related concepts. To investigate the neurocognitive basis of individual differences in these aspects of semantic control, we used resting-state fMRI to characterise the intrinsic connectivity of left ventrolateral prefrontal cortex (VLPFC), implicated in controlled retrieval, and examined associations on a paced serial semantic task, in which participants were asked to detect category members amongst distractors. This task manipulated both the strength of target associations and the requirement to sustain retrieval within a narrow semantic category over time. We found that individuals with stronger connectivity between VLPFC and medial prefrontal cortex within the default mode network (DMN) showed better retrieval of strong associations (which are thought to be recalled more automatically). Stronger connectivity between the same VLPFC seed and another DMN region in medial parietal cortex was associated with larger declines in retrieval over the course of the category. In contrast, participants with stronger connectivity between VLPFC and cognitive control regions within the ventral attention network (VAN) had better controlled retrieval of weak associations and were better able to sustain their comprehension throughout the category. These effects overlapped in left insular cortex within the VAN, indicating that a common pattern of connectivity is associated with different aspects of controlled semantic retrieval induced by both the structure of long-term knowledge and the sustained retrieval of related information.

The pathological mechanism of attention deficit hyperactivity disorder (ADHD) is incompletely specified, which leads to difficulty in precise diagnosis. Functional magnetic resonance imaging (fMRI) has emerged as a common neuroimaging technique for studying the brain functional connectome. Most existing methods that have either ignored or simply utilized graph structure, do not fully leverage the potentially important topological information which may be useful in characterizing brain disorders. There is a crucial need for designing novel and efficient approaches which can capture such information. To this end, we propose a new dynamic graph convolutional network (dGCN), which is trained with sparse brain regional connections from dynamically calculated graph features. We also develop a novel convolutional readout layer to improve graph representation. Our extensive experimental analysis demonstrates significantly improved performance of dGCN for ADHD diagnosis compared with existing machine learning and deep learning methods. Visualizations of the salient regions of interest (ROIs) and connectivity based on informative features learned by our model show that the identified functional abnormalities mainly involve brain regions in temporal pole, gyrus rectus, and cerebellar gyri from temporal lobe, frontal lobe, and cerebellum, respectively. A positive correlation was further observed between the identified connectomic abnormalities and ADHD symptom severity. The proposed dGCN model shows great promise in providing a functional network-based precision diagnosis of ADHD and is also broadly applicable to brain connectome-based study of mental disorders.

Studies assessing relationships between neural and cognitive changes in healthy aging have shown that a variety of aspects of brain structure and function explain a significant portion of the variability in cognitive outcomes throughout adulthood. Many studies assessing relationships between brain function and cognition have utilized time-averaged, or static functional connectivity methods to explore ways in which brain network organization may contribute to aspects of cognitive aging. However, recent studies in this field have suggested that time-varying, or dynamic measures of functional connectivity, which assess changes in functional connectivity over the course of a scan session, may play a stronger role in explaining cognitive outcomes in healthy young adults. Further, both static and dynamic functional connectivity studies suggest that there may be differences in patterns of brain-cognition relationships as a function of whether or not the participant is performing a task during the scan. Thus, the goals of the present study were threefold: (1) assess whether neural flexibility during both resting as well as task-based scans is related to participant age and cognitive performance in a lifespan aging sample, (2) determine whether neural flexibility moderates relationships between age and cognitive performance, and (3) explore differences in neural flexibility between rest and task. Participants in the study were 386 healthy adults between the ages of 20–80 who provided resting state and/or task-based (Matrix Reasoning) functional magnetic resonance imaging (fMRI) scan data as part of their participation in two ongoing studies of cognitive aging. Neural flexibility measures from both resting and task-based scans reflected the number of times each node changed network assignment, and were averaged both across the whole brain (global neural flexibility) as well as within ten somatosensory/cognitive networks. Results showed that neural flexibility was not related to participant age, and that task-based global neural flexibility, as well as task-based neural flexibility in several networks, tended to be negatively related to reaction times during the Matrix Reasoning task, however these effects did not survive strict multiple comparisons correction. Resting state neural flexibility was not significantly related to either participant age or cognitive performance. Additionally, no neural flexibility measures significantly moderated relationships between participant age and cognitive outcomes. Further, neural flexibility differed as a function of scan type, with resting state neural flexibility being significantly greater than task-based neural flexibility. Thus, neural flexibility measures computed during a cognitive task may be more meaningfully related to cognitive performance across the adult lifespan then resting state measures of neural flexibility.

The brain consumes the most energy per relative mass amongst the organs in the human body. Theoretical and empirical studies have shown that behavioral processes are relatively inexpensive metabolically, and that most energy goes to maintaining the status quo, i.e., the balance of cell membranes' resting potentials and subthreshold spontaneous activity. Spontaneous activity fluctuates across brain regions in a correlated fashion that defines multi-scale hierarchical networks called resting-state networks (RSNs). Different regions of the brain display different metabolic consumption, but the relationship between regional brain metabolism and RSNs is still under investigation. Here, we examine the variability of glucose metabolism across brain regions, measured with the relative standard uptake value (SUVR) using

Dynamic Causal Modeling (DCM) is a Bayesian framework for inferring on hidden (latent) neuronal states, based on measurements of brain activity. Since its introduction in 2003 for functional magnetic resonance imaging data, DCM has been extended to electrophysiological data, and several variants have been developed. Their biophysically motivated formulations make these models promising candidates for providing a mechanistic understanding of human brain dynamics, both in health and disease. However, due to their complexity and reliance on concepts from several fields, fully understanding the mathematical and conceptual basis behind certain variants of DCM can be challenging. At the same time, a solid theoretical knowledge of the models is crucial to avoid pitfalls in the application of these models and interpretation of their results. In this paper, we focus on one of the most advanced formulations of DCM, i.e. conductance-based DCM for cross-spectral densities, whose components are described across multiple technical papers. The aim of the present article is to provide an accessible exposition of the mathematical background, together with an illustration of the model's behavior. To this end, we include step-by-step derivations of the model equations, point to important aspects in the software implementation of those models, and use simulations to provide an intuitive understanding of the type of responses that can be generated and the role that specific parameters play in the model. Furthermore, all code utilized for our simulations is made publicly available alongside the manuscript to allow readers an easy hands-on experience with conductance-based DCM., NeuroImage, 245, ISSN:1053-8119, ISSN:1095-9572

The ability to extract cues from faces is fundamental for social animals, including humans. An individual's profile of functional connectivity across a face network can be shaped by common organizing principles, stable individual traits, and time-varying mental states. In the present study, we used data obtained with functional magnetic resonance imaging in two cohorts, IMAGEN (N = 534) and ALSPAC (N = 465), to investigate - both at group and individual levels - the consistency of the regional profile of functional connectivity across populations (IMAGEN, ALSPAC) and time (Visits 1 to 3 in IMAGEN; age 14 to 22 years). At the group level, we found a robust canonical profile of connectivity both across populations and time. At the individual level, connectivity profiles deviated from the canonical profile, and the magnitude of this deviation related to the presence of psychopathology. These findings suggest that the brain processes faces in a highly stereotypical manner, and that the deviations from this normative pattern may be related to the risk of mental illness.

Prior studies have shown that the left posterior superior temporal sulcus (pSTS) and left temporo-parietal junction (TPJ) both contribute to phonological short-term memory, speech perception and speech production. Here, by conducting a within-subjects multi-factorial fMRI study, we dissociate the response profiles of these regions and a third region – the anterior ascending terminal branch of the left superior temporal sulcus (atSTS), which lies dorsal to pSTS and ventral to TPJ. First, we show that each region was more activated by (i) 1-back matching on visually presented verbal stimuli (words or pseudowords) compared to 1-back matching on visually presented non-verbal stimuli (pictures of objects or non-objects), and (ii) overt speech production than 1-back matching, across 8 types of stimuli (visually presented words, pseudowords, objects and non-objects and aurally presented words, pseudowords, object sounds and meaningless hums). The response properties of the three regions dissociated within the auditory modality. In left TPJ, activation was higher for auditory stimuli that were non-verbal (sounds of objects or meaningless hums) compared to verbal (words and pseudowords), irrespective of task (speech production or 1-back matching). In left pSTS, activation was higher for non-semantic stimuli (pseudowords and hums) than semantic stimuli (words and object sounds) on the dorsal pSTS surface (dpSTS), irrespective of task. In left atSTS, activation was not sensitive to either semantic or verbal content. The contrasting response properties of left TPJ, dpSTS and atSTS was cross-validated in an independent sample of 59 participants, using region-by-condition interactions. We also show that each region participates in non-overlapping networks of frontal, parietal and cerebellar regions. Our results challenge previous claims about functional specialisation in the left posterior superior temporal lobe and motivate future studies to determine the timing and directionality of information flow in the brain networks involved in speech perception and production.

The search for an ‘ideal’ approach to investigate the functional connections in the human brain is an ongoing challenge for the neuroscience community. While resting-state functional magnetic resonance imaging (fMRI) has been widely used to study individual functional connectivity patterns, recent work has highlighted the benefits of collecting functional connectivity data while participants are exposed to naturalistic stimuli, such as watching a movie or listening to a story. For example, functional connectivity data collected during movie-watching were shown to predict cognitive and emotional scores more accurately than resting-state-derived functional connectivity. We have previously reported a tight link between resting-state functional connectivity and task-derived neural activity, such that the former successfully predicts the latter. In the current work we use data from the Human Connectome Project to demonstrate that naturalistic-stimulus-derived functional connectivity predicts task-induced brain activation maps more accurately than resting-state-derived functional connectivity. We then show that activation maps predicted using naturalistic stimuli are better predictors of individual intelligence scores than activation maps predicted using resting-state. We additionally examine the influence of naturalistic-stimulus type on prediction accuracy. Our findings emphasize the potential of naturalistic stimuli as a promising alternative to resting-state fMRI for connectome-based predictive modelling of individual brain activity and cognitive traits.

Both cortical and subcortical regions can be functionally organized into networks. Regions of the basal ganglia are extensively interconnected with the cortex via reciprocal connections that relay and modulate cortical function. Here we employ an edge-centric approach, which computes co-fluctuations among region pairs in a network to investigate the role and interaction of subcortical regions with cortical systems. By clustering edges into communities, we show that cortical systems and subcortical regions couple via multiple edge communities, with hippocampus and amygdala having a distinct pattern from striatum and thalamus. We show that the edge community structure of cortical networks is highly similar to one obtained from cortical nodes when the subcortex is present in the network. Additionally, we show that the edge community profile of both cortical and subcortical nodes can be estimates solely from cortico-subcortical interactions. Finally, we used a motif analysis focusing on edge community triads where a subcortical region coupled to two cortical regions and found that two community triads where one community couples the subcortex to the cortex were overrepresented. In summary, our results show organized coupling of the subcortex to the cortex that may play a role in cortical organization of primary sensorimotor/attention and heteromodal systems and puts forth the motif analysis of edge community triads as a promising method for investigation of communication patterns in networks.

Introduction Resting-state oscillatory activity has been extensively studied across a wide array of disorders. Establishing which spectrally- and spatially-specific oscillatory components exhibit test-retest reliability is essential to move the field forward. While studies have shown short-term reliability of MEG resting-state activity, no studies have examined test-retest reliability across an extended period of time to establish the stability of these signals, which is critical for reproducibility. Methods We examined 18 healthy adults age 23 - 61 who completed three visits across three years. For each visit, participants completed both a resting state MEG and structural MRI scan. MEG data were source imaged, and the cortical power in canonical frequency bands (delta, theta, alpha, beta, low gamma, high gamma) was computed. Intra-class correlation coefficients (ICC) were then calculated across the cortex for each frequency band. Results Over three years, power in the alpha and beta bands displayed the highest reliability estimates, while gamma showed the lowest estimates of three-year reliability. Spatially, delta, alpha, and beta all showed the highest degrees of reliability in the parietal cortex. Interestingly, the peak signal for each of these frequency bands was located outside of the parietal cortex, suggesting that reliability estimates were not solely dependent on the signal-to-noise ratio. Conclusion Oscillatory resting-state power in parietal delta, posterior beta, and alpha across most of the cortex are reliable across three years and future MEEG studies may focus on these measures for the development of specific markers.

Current views on the neural network subserving reading and its deficits in dyslexia rely largely on evidence derived from functional neuroimaging studies. However, understanding the structural organization of reading and its aberrations in dyslexia requires a hodological approach, studies of which have not provided consistent findings. Here, we adopted a whole brain hodological approach and investigated relationships between structural white matter connectivity and reading skills and phonological processing in a cross-sectional study of 44 adults using individual local connectome matrix from diffusion MRI data. Moreover, we performed quantitative anisotropy aided differential tractography to uncover structural white matter anomalies in dyslexia (23 dyslexics and 21 matched controls) and their correlation to reading-related skills. The connectometry analyses indicated that reading skills and phonological processing were both associated with corpus callosum (tapetum), forceps major and minor, as well as cerebellum bilaterally. Furthermore, the left dorsal and right thalamic pathways were associated with phonological processing. Differential tractography analyses revealed structural white matter anomalies in dyslexics in the left ventral route and bilaterally in the dorsal route compared to the controls. Connectivity deficits were also observed in the corpus callosum, forceps major, vertical occipital fasciculus and corticostriatal and thalamic pathways. Altered structural connectivity in the observed differential tractography results correlated with poor reading skills and phonological processing. Using a hodological approach, the current study provides novel evidence for the extent of the reading-related connectome and its aberrations in dyslexia. The results conform current functional neuroanatomical models of reading and developmental dyslexia but provide novel network-level and tract-level evidence on structural connectivity anomalies in dyslexia, including the vertical occipital fasciculus.

Functional connectivity (FC) networks are typically inferred from resting-state fMRI data using the Pearson correlation between BOLD time series from pairs of brain regions. However, alternative methods of estimating functional connectivity have not been systematically tested for their sensitivity or robustness to head motion artifact. Here, we evaluate the sensitivity of eight different functional connectivity measures to motion artifact using resting-state data from the Human Connectome Project. We report that FC estimated using full correlation has a relatively high residual distance-dependent relationship with motion compared to partial correlation, coherence, and information theory-based measures, even after implementing rigorous methods for motion artifact mitigation. This disadvantage of full correlation, however, may be offset by higher test-retest reliability, fingerprinting accuracy, and system identifiability. FC estimated by partial correlation offers the best of both worlds, with low sensitivity to motion artifact and intermediate system identifiability, with the caveat of low test-retest reliability and fingerprinting accuracy. We highlight spatial differences in the sub-networks affected by motion with different FC metrics. Further, we report that intra-network edges in the default mode and retrosplenial temporal sub-networks are highly correlated with motion in all FC methods. Our findings indicate that the method of estimating functional connectivity is an important consideration in resting-state fMRI studies and must be chosen carefully based on the parameters of the study.

Research into the human connectome (i.e., all connections in the human brain) with the use of resting state functional MRI has rapidly increased in popularity in recent years, especially with the growing availability of large-scale neuroimaging datasets. The goal of this review article is to describe innovations in functional connectome representations that have come about in the past 8 years, since the 2013 NeuroImage special issue on ‘Mapping the Connectome’. In the period, research has shifted from group-level brain parcellations towards the characterization of the individualized connectome and of relationships between individual connectomic differences and behavioral/clinical variation. Achieving subject-specific accuracy in parcel boundaries while retaining cross-subject correspondence is challenging, and a variety of different approaches are being developed to meet this challenge, including improved alignment, improved noise reduction, and robust group-to-subject mapping approaches. Beyond the interest in the individualized connectome, new representations of the data are being studied to complement the traditional parcellated connectome representation (i.e., pairwise connections between distinct brain regions), such as methods that capture overlapping and smoothly varying patterns of connectivity (‘gradients’). These different connectome representations offer complimentary insights into the inherent functional organization of the brain, but challenges for functional connectome research remain. Interpretability will be improved by future research towards gaining insights into the neural mechanisms underlying connectome observations obtained from functional MRI. Validation studies comparing different connectome representations are also needed to build consensus and confidence to proceed with clinical trials that may produce meaningful clinical translation of connectome insights.

Highlights • Applied deep learning to sex classification in UK BioBank fMRI connectomes. • Deep learning classifies sex better in resting-state than in task fMRI. • Algorithm to balance out multiple confounds from an fMRI dataset. • Adapted two deep learning visualization methods to fMRI connectome classification. • Analyzed role of three brain a priori networks in sex classification., Classification of whole-brain functional connectivity MRI data with convolutional neural networks (CNNs) has shown promise, but the complexity of these models impedes understanding of which aspects of brain activity contribute to classification. While visualization techniques have been developed to interpret CNNs, bias inherent in the method of encoding abstract input data, as well as the natural variance of deep learning models, detract from the accuracy of these techniques. We introduce a stochastic encoding method in an ensemble of CNNs to classify functional connectomes by sex. We applied our method to resting-state and task data from the UK BioBank, using two visualization techniques to measure the salience of three brain networks involved in task- and resting-states, and their interaction. To regress confounding factors such as head motion, age, and intracranial volume, we introduced a multivariate balancing algorithm to ensure equal distributions of such covariates between classes in our data. We achieved a final AUROC of 0.8459. We found that resting-state data classifies more accurately than task data, with the inner salience network playing the most important role of the three networks overall in classification of resting-state data and connections to the central executive network in task data.

Parent-child book reading is important for fostering the development of various lifelong cognitive and social abilities in young children. Despite numerous reports describing the effects of familiarity on shared reading for children, the exact neural basis of the functional network architecture remains unclear. We conducted Magnet-Encephalographic (MEG) experiments using graph theory to elucidate the role of familiarity in shared reading in a child's brain network and to measure the connectivity dynamics of a child while Listening to Storybook Reading (LSBR), which represents the daily activity of shared book reading between the child and caregiver. The LSBR task was performed with normally developing preschool- and school-age children (N = 15) under two conditions: reading by their own mother (familiar condition) vs. an experimenter (unfamiliar condition). We used the phase lag index (PLI), which captures synchronization of MEG signals, to estimate functional connectivity. For the whole brain network topology, an undirected weighted graph was produced using 68 brain regions as nodes and interregional PLI values as edges for five frequency bands. Behavioral data (i.e., the degree of attention and facial expressions) were evaluated from video images of the child's face during the two conditions. Our results showed enhanced widespread functional connectivity in the alpha band during the mother condition. In the mother condition, the whole brain network in the alpha band exhibited topographically high local segregation with high global integration, indicating an increased small-world property. Results of the behavioral analysis revealed that children were more attentive and showed more positive facial expressions in the mother condition than in the experimenter condition. Behavioral data were significantly correlated with graph metrics in the mother condition but not in the experimenter condition. In this study, we identified the neural correlates of a familiarity effect in children's brain connectivity dynamics during LSBR. Furthermore, these familiarity-related brain dynamics were closely linked to the child's behavior. Graph theory applied to MEG data may provide useful insight into the familiarity-related child brain response in a naturalistic setting and its relevance to child attitudes.

Ongoing brain function is largely determined by the underlying wiring of the brain, but the specific rules governing this relationship remain unknown. Emerging literature has suggested that functional interactions between brain regions emerge from the structural connections through mono- as well as polysynaptic mechanisms. Here, we propose a novel approach based on diffusion maps and Riemannian optimization to emulate this dynamic mechanism in the form of random walks on the structural connectome and predict functional interactions as a weighted combination of these random walks. Our proposed approach was evaluated in two different cohorts of healthy adults (Human Connectome Project, HCP; Microstructure-Informed Connectomics, MICs). Our approach outperformed existing approaches and showed that performance plateaus approximately around the third random walk. At macroscale, we found that the largest number of walks was required in nodes of the default mode and frontoparietal networks, underscoring an increasing relevance of polysynaptic communication mechanisms in transmodal cortical networks compared to primary and unimodal systems.

The adult brain is organized into distinct functional networks, forming the basis of information processing and determining individual differences in behavior. Is this network organization genetically determined and present at birth? And what is the individual variability in this organization in neonates? Here, we use unsupervised learning to uncover intrinsic functional brain organization using resting-state connectivity from a large cohort of neonates (Developing Human Connectome Project). We identified a set of symmetric, hierarchical, and replicable networks: sensorimotor, visual, default mode, ventral attention, and high-level vision. We quantified individual variability across neonates, and found the most individual variability in the ventral attention networks. Crucially, the variability of these networks was not driven by SNR differences or differences from adult networks (Yeo et al., 2011). Finally, differential gene expression provided a potential explanation for the emergence of these distinct networks and identified potential genes of interest for future developmental and individual variability research. Overall, we found neonatal connectomes (even at the voxel-level) can reveal broad individual-specific information processing units. The presence of individual differences in neonates and the framework for personalized parcellations demonstrated here has the potential to improve prediction of behavior and future outcomes from neonatal and infant brain data.

Dynamic resting state functional connectivity (RSFC) characterizes fluctuations that occur over time in functional brain networks. Existing methods to extract dynamic RSFCs, such as sliding-window and clustering methods that are inherently non-adaptive, have various limitations such as high-dimensionality, an inability to reconstruct brain signals, insufficiency of data for reliable estimation, insensitivity to rapid changes in dynamics, and a lack of generalizability across multiply functional imaging modalities. To overcome these deficiencies, we develop a novel and unifying time-varying dynamic network (TVDN) framework for examining dynamic resting state functional connectivity. TVDN includes a generative model that describes the relation between a low-dimensional dynamic RSFC and the brain signals, and an inference algorithm that automatically and adaptively learns the low-dimensional manifold of dynamic RSFC and detects dynamic state transitions in data. TVDN is applicable to multiple modalities of functional neuroimaging such as fMRI and MEG/EEG. The estimated low-dimensional dynamic RSFCs manifold directly links to the frequency content of brain signals. Hence we can evaluate TVDN performance by examining whether learnt features can reconstruct observed brain signals. We conduct comprehensive simulations to evaluate TVDN under hypothetical settings. We then demonstrate the application of TVDN with real fMRI and MEG data, and compare the results with existing benchmarks. Results demonstrate that TVDN is able to correctly capture the dynamics of brain activity and more robustly detect brain state switching both in resting state fMRI and MEG data.

Functional magnetic resonance imaging (fMRI) has been widely employed to study stroke pathophysiology. In particular, analyses of fMRI signals at rest were directed at quantifying the impact of stroke on spatial features of brain networks. However, brain networks have intrinsic time features that were, so far, disregarded in these analyses. In consequence, standard fMRI analysis failed to capture temporal imbalance resulting from stroke lesions, hence restricting their ability to reveal the interdependent pathological changes in structural and temporal network features following stroke. Here, we longitudinally analyzed hemodynamic-informed transient activity in a large cohort of stroke patients (n = 103) to assess spatial and temporal changes of brain networks after stroke. Metrics extracted from the hemodynamic-informed transient activity were replicable within- and between-individuals in healthy participants, hence supporting their robustness and their clinical applicability. While large-scale spatial patterns of brain networks were preserved after stroke, their durations were altered, with stroke subjects exhibiting a varied pattern of longer and shorter network activations compared to healthy individuals. Specifically, patients showed a longer duration in the lateral precentral gyrus and anterior cingulum, and a shorter duration in the occipital lobe and in the cerebellum. These temporal alterations were associated with white matter damage in projection and association pathways. Furthermore, they were tied to deficits in specific behavioral domains as restoration of healthy brain dynamics paralleled recovery of cognitive functions (attention, language and spatial memory), but was not significantly correlated to motor recovery. These findings underscore the critical importance of network temporal properties in dissecting the pathophysiology of brain changes after stroke, thus shedding new light on the clinical potential of time-resolved methods for fMRI analysis.

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, “tractometry” can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named “myelin streamline decomposition” (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

Resting functional MRI studies of the infant brain are increasingly becoming an important tool in developmental neuroscience. Whereas the test-retest reliability of functional connectivity (FC) measures derived from resting fMRI data have been characterized in the adult and child brain, similar assessments have not been conducted in infants. In this study, we examined the intra-session test-retest reliability of FC measures from 119 infant brain MRI scans from four neurodevelopmental studies. We investigated edge-level and subject-level reliability within one MRI session (between and within runs) measured by the Intraclass correlation coefficient (ICC). First, using an atlas-based approach, we examined whole-brain connectivity as well as connectivity within two common resting fMRI networks – the default mode network (DMN) and the sensorimotor network (SMN). Second, we examined the influence of run duration, study site, and scanning manufacturer (e.g., Philips and General Electric) on ICCs. Lastly, we tested spatial similarity using the Jaccard Index from networks derived from independent component analysis (ICA). Consistent with resting fMRI studies from adults, our findings indicated poor edge-level reliability (ICC = 0.14 - 0.18), but moderate-to-good subject-level intra-session reliability for whole-brain, DMN, and SMN connectivity (ICC = 0.40 - 0.78). We also found significant effects of run duration, site, and scanning manufacturer on reliability estimates. Some ICA-derived networks showed strong spatial reproducibility (e.g., DMN, SMN, and Visual Network), and were labelled based on their spatial similarity to analogous networks measured in adults. These networks were reproducibly found across different study studies. However, other ICA-networks (e.g. Executive Control Network) did not show strong spatial reproducibility, suggesting that the reliability and/or maturational course of functional connectivity may vary by network. In sum, our findings suggest that developmental scientist may be on safe ground examining the functional organization of some major neural networks (e.g. DMN and SMN), but judicious interpretation of functional connectivity is essential to its ongoing success.HighlightsInfant functional connectivity (FC) shows poor edge-level reliability (ICCs)However, subject-level infant FC estimates show good-to-excellent ICCsSpatial reproducibility is better for some resting networks (DMN, SMN) than others (ECN)Reliability estimates differ across study site and MRI scannerConclusion - Infant FC can be a reliable measurement, but judicious use is needed

Infants’ socio-cognitive ability develops dramatically during the first year of life. From the perspective of ontogeny, the early development of social behavior allows for parent-child attachment, which in turn enhances survival. Thus, it is theorized that the development of social behavior, driven by social brain networks, forms the core of developmental acquisitions during this period. Further, understanding the maturation within the neural networks during social development is crucial to obtain a better grasp of the development of social developmental disorders. Therefore, we performed a longitudinal study in 854 infants measured at around 5 and 10 months to map the development of functional networks in the brain when infants were processing social and non-social videos. Using EEG, we focused on the frequency bands most commonly connected to social behavior: theta and alpha. We found that alpha networks remained relatively stable over the first year of life and showed no selectivity for social versus non-social stimuli, theta networks, showed strong global reconfigurations. The development of the theta networks progressed from a parietal occipital network in early infancy to a frontoparietal network towards the end of the first year of life. This reconfiguration coincided with selectivity for social versus non-social stimuli, with infants approaching the end of their first year of life showing increased synchronicity of theta communication when watching social videos versus non-social videos. Our findings provide strong evidence for the involvement of a frontoparietal theta network in the development of the social brain.

Humans use different spatial reference frames (allocentric or egocentric) to navigate successfully toward their destination in different spatial scale spaces (environmental or vista). However, it remains unclear how the brain represents different spatial scales and different spatial reference frames. Thus, we conducted an activation likelihood estimation (ALE) meta-analysis of 47 fMRI articles involving human spatial navigation. We found that both the environmental and vista spaces activated the parahippocampal place area (PPA), retrosplenial complex (RSC), and occipital place area in the right hemisphere. The environmental space showed stronger activation than the vista space in the occipital and frontal regions. No brain region exhibited stronger activation for the vista than the environmental space. The allocentric and egocentric reference frames activated the bilateral PPA and right RSC. The allocentric frame showed more stronger activations than the egocentric frame in the right culmen, left middle frontal gyrus, and precuneus. No brain region displayed stronger activation for the egocentric than the allocentric navigation. Our findings suggest that navigation in different spatial scale spaces can evoke specific and common brain regions, and that the brain regions representing spatial reference frames are not absolutely separated.

Brain connectome fingerprinting is rapidly rising as a novel influential field in brain network analysis. Yet, it is still unclear whether connectivity fingerprints could be effectively used for mapping and predicting disease progression from human brain data. We hypothesize that dysregulation of brain activity in disease would reflect in worse subject identification. We propose a novel framework, Clinical Connectome Fingerprinting, to detect individual connectome features from clinical populations. We show that “clinical fingerprints” can map individual variations between elderly healthy subjects and patients with mild cognitive impairment in functional connectomes extracted from magnetoencephalography data. We find that identifiability is reduced in patients as compared to controls, and show that these connectivity features are predictive of the individual Mini-Mental State Examination (MMSE) score in patients. We hope that the proposed methodology can help in bridging the gap between connectivity features and biomarkers of brain dysfunction in large-scale brain networks.

Highlights • We propose a novel end-to-end neural network that employs resting-state and task-based functional MRI (fMRI) datasets, obtained one month after trauma exposure, to predict PTSD one, six and 14-months after the exposure. • The method utilizes connectivity maps extracted from pairs of brain regions which are subsequently updated by applying the algorithmic technique of pairwise attention. • The proposed deep learning method predicts PTSD status, PTSD symptom clusters and survival analysis within the prospective design. We demonstrate a significant improvement in performance on all the datasets and experiments in comparison to other relevant analytical techniques. • Pairwise association analysis reveals several significant functional connectivity patterns, in line with previous PTSD neuroimaging literature., Early intervention following exposure to a traumatic life event could change the clinical path from the development of post traumatic stress disorder (PTSD) to recovery, hence the interest in early detection and underlying biological mechanisms involved in the development of post traumatic sequelae. We introduce a novel end-to-end neural network that employs resting-state and task-based functional MRI (fMRI) datasets, obtained one month after trauma exposure, to predict PTSD symptoms at one-, six- and fourteen-months after the exposure. FMRI data, as well as PTSD status and symptoms, were collected from adults at risk for PTSD development, after admission to emergency room following a traumatic event. Our computational method utilized a per-region encoder to extract brain regions embedding, which were subsequently updated by applying the algorithmic technique of pairwise attention. The affinities obtained between each pair of regions were combined to create a pairwise co-activation map used to perform multi-label classification. The results demonstrate that the novel method’s performance in predicting PTSD symptoms, in a prospective manner, outperforms previous analytical techniques reported in the fMRI literature, all trained on the same dataset. We further show a high predictive ability for predicting PTSD symptom clusters and PTSD persistence. To the best of our knowledge, this is the first deep learning method applied on fMRI data with respect to prospective clinical outcomes, to predict PTSD status, severity and symptom clusters. Future work could further delineate the mechanisms that underlie such a prediction, and potentially improve single patient characterization.

Temporal concatenation group ICA (TC-GICA) is a widely used data-driven method to extract common functional brain networks among individuals. TC-GICA concatenates the time series of individual fMRI data and applies dimension reduction and ICA algorithms to decompose the data into group-level components. The default mode network (DMN) estimated using TC-GICA at relatively high model orders (i.e., large numbers of components) is split into multiple components. The split DMNs are topographically different from those estimated using other methods (e.g., seed-based correlation, clustering, graph theoretical analysis, and other ICA methods like gRAICAR and IVA-GL) and are inconsistent with the existing knowledge of DMN. We hypothesize that the ''DMN-splitting'' phenomenon reflects the impact of inter-individual variability in data, which is propagated into the ICA decomposition via the data-concatenation step of TC-GICA. By systematically manipulating the amount of variability involved in the temporal concatenation in both simulated and several realistic datasets, we observed that as more variability was involved, the estimated DMN became less similar to the averaged functional connectivity (FC) pattern obtained using seed-based correlation analysis. The performance of the DMN estimation in TC-GICA also exhibited remarkable dependence on the model order settings. Further analyses revealed that the “DMN-splitting” in TC-GICA could be reproduced when involving large variability in the data-concatenation and performing ICA at high model orders. These results were replicated across multiple datasets and various software implementations. When applying ICA approaches that avoid temporal concatenation, such as gRAICAR and IVA-GL, to the same datasets, the estimated group-level DMN was more consistent with the seed-based FC pattern and was more robust to various model order settings. This study calls for caution when applying TC-GICA to datasets expected to have large inter-individual variability, such as pooling different experimental groups of subjects.

Anesthetics are known to disrupt neural interactions in cortical and subcortical brain circuits. While the effect of anesthetic drugs on consciousness is reversible, the neural mechanism mediating induction and recovery may be different. Insight into these distinct mechanisms can be gained from a systematic comparison of neural dynamics during slow induction of and emergence from anesthesia. To this end, we used functional magnetic resonance imaging (fMRI) data obtained in healthy volunteers before, during, and after the administration of propofol at incrementally adjusted target concentrations. We analyzed functional connectivity of corticocortical and subcorticocortical networks and the temporal autocorrelation of fMRI signal as an index of neural processing timescales. We found that en route to unconsciousness, temporal autocorrelation across the entire brain gradually increased, whereas functional connectivity gradually decreased. In contrast, regaining consciousness was associated with an abrupt restoration of cortical but not subcortical temporal autocorrelation and an abrupt boost of subcorticocortical functional connectivity. Pharmacokinetic effects could not account for the difference in neural dynamics between induction and emergence. We conclude that the induction and recovery phases of anesthesia follow asymmetric neural dynamics. A rapid increase in the speed of cortical neural processing and subcorticocortical neural interactions may be a mechanism that reboots consciousness.

Brain signal variability has been consistently linked to functional integration; however, whether this coupling is associated with cognitive functions and/or psychiatric diseases has not been clarified. Using multiple multimodality datasets, including resting-state functional magnetic resonance imaging (rsfMRI) data from the Human Connectome Project (HCP: N = 927) and a Beijing sample (N = 416) and cerebral blood flow (CBF) and rsfMRI data from a Hangzhou sample (N = 29), we found that, compared with the existing variability measure (i.e., SDBOLD), the mean-scaled (standardized) fractional standard deviation of the BOLD signal (mfSDBOLD) maintained very high test-retest reliability, showed greater cross-site reliability and was less affected by head motion. We also found strong reproducible couplings between the mfSDBOLD and functional integration measured by the degree centrality (DC), both cross-voxel and cross-subject, which were robust to scanning and preprocessing parameters. Moreover, both mfSDBOLD and DC were correlated with CBF, suggesting a common physiological basis for both measures. Critically, the degree of coupling between mfSDBOLD and long-range DC was positively correlated with individuals’ cognitive total composite scores. Brain regions with greater mismatches between mfSDBOLD and long-range DC were more vulnerable to brain diseases. Our results suggest that BOLD signal variability could serve as a meaningful index of local function that underlies functional integration in the human brain and that a strong coupling between BOLD signal variability and functional integration may serve as a hallmark of balanced brain networks that are associated with optimal brain functions.

Functional connectivity (FC) and resting-state network (RSN) analyses using functional magnetic resonance imaging (fMRI) have evolved into a growing field of research and have provided useful biomarkers for the assessment of brain function in neurological disorders. However, the underlying mechanisms of the blood oxygen level-dependant (BOLD) signal are not fully resolved due to its inherent complexity. In contrast, [18F]fluorodeoxyglucose positron emission tomography ([18F]FDG-PET) has been shown to provide a more direct measure of local synaptic activity and may have additional value for the readout and interpretation of brain connectivity. We performed an RSN analysis from simultaneously acquired PET/fMRI data on a single-subject level to directly compare fMRI and [18F]FDG-PET-derived networks during the resting state. Simultaneous [18F]FDG-PET/fMRI scans were performed in 30 rats. Pairwise correlation analysis, as well as independent component analysis (ICA), were used to compare the readouts of both methods. We identified three RSNs with a high degree of similarity between PET and fMRI-derived readouts: the default-mode-like network (DMN), the basal ganglia network and the cerebellar-midbrain network. Overall, [18F]FDG connectivity indicated increased integration between different, often distant, brain areas compared to the results indicated by the more segregated fMRI-derived FC. Additionally, several networks exclusive to either modality were observed using ICA. These networks included mainly bilateral cortical networks of a limited spatial extent for fMRI and more spatially widespread networks for [18F]FDG-PET, often involving several subcortical areas. This is the first study using simultaneous PET/fMRI to report RSNs subject-wise from dynamic [18F]FDG tracer delivery and BOLD fluctuations with both independent component analysis (ICA) and pairwise correlation analysis in small animals. Our findings support previous studies, which show a close link between local synaptic glucose consumption and BOLD-fMRI-derived FC. However, several brain regions were exclusively attributed to either [18F]FDG or BOLD-derived networks underlining the complementarity of this hybrid imaging approach, which may contribute to the understanding of brain functional organization and could be of interest for future clinical applications.

Transcranial alternating current stimulation (tACS), applied to two brain sites with different phase lags, has been shown to modulate stimulation-outlasting functional EEG connectivity between the targeted regions. Given the lack of knowledge on mechanisms of tACS aftereffects, it is difficult to further enhance effect sizes and reduce variability in experiments. In this computational study, we tested if spike-timing-dependent plasticity (STDP) can explain stimulation-outlasting connectivity modulation by dual-site tACS and explored the effects of tACS parameter choices. Two populations of spiking neurons were coupled with synapses subject to STDP, and results were validated via a re-analysis of EEG data. Our simulations showed stimulation-outlasting connectivity changes between in- and anti-phase tACS, dependent on both tACS frequency and synaptic conduction delays. Importantly, both a simple network entraining to a wide range of tACS frequencies as well as a more realistic network that spontaneously oscillated at alpha frequency predicted that the largest effects would occur for short conduction delays between the stimulated regions. This finding agreed with experimental EEG connectivity modulation by 10 Hz tACS, showing a clear negative correlation of tACS effects with estimated conduction delays between regions. In conclusion, STDP can explain connectivity aftereffects of dual-site tACS. However, not all combinations of tACS frequency and application sites are expected to effectively modulate connectivity via STDP. We therefore suggest using appropriate computational models and/or EEG analysis for planning and interpretation of dual-site tACS studies relying on aftereffects.

Visual shape completion recovers object shape, size, and number from spatially segregated edges. Despite being extensively investigated, the process's underlying brain regions, networks, and functional connections are still not well understood. To shed light on the topic, we scanned (fMRI) healthy adults during rest and during a task in which they discriminated pac-man configurations that formed or failed to form completed shapes (illusory and fragmented condition, respectively). Task activation differences (illusory-fragmented), resting-state functional connectivity, and multivariate pattern differences were identified on the cortical surface using 360 predefined parcels and 12 functional networks composed of such parcels. Brain activity flow mapping (ActFlow) was used to evaluate the likely involvement of resting-state connections for shape completion. We identified 36 differentially-active parcels including a posterior temporal region, PH, whose activity was consistent across 95% of observers. Significant task regions primarily occupied the secondary visual network but also incorporated the frontoparietal, dorsal attention, default mode, and cingulo-opercular networks. Each parcel's task activation difference could be modeled via its resting-state connections with the remaining parcels (r=.62, p

Accurate epileptogenic zone (EZ) or seizure onset zone (SOZ) localization is crucial for epilepsy surgery optimization. Previous animal and human studies on epilepsy have reported that changes in blood oxygen level-dependent (BOLD) signals induced by epileptic events could be used as diagnostic markers for EZ or SOZ localization. Simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI) recording is gaining interest as a non-invasive tool for preoperative epilepsy evaluation. However, EEG-fMRI studies have reported inconsistent and ambiguous findings. Therefore, it remains unclear whether BOLD responses can be used for accurate EZ or SOZ localization. In this study, we used simultaneous EEG-fMRI recording in a rat model of 4-aminopyridine-induced acute focal seizures to assess the spatial concordance between individual BOLD responses and the SOZ. This was to determine the optimal use of simultaneous EEG-fMRI recording in the SOZ localization. We observed a high spatial consistency between BOLD responses and the SOZ. Further, dynamic BOLD responses were consistent with the regions where the seizures were propagated. These results suggested that simultaneous EEG-fMRI recording could be used as a noninvasive clinical diagnostic technique for localizing the EZ or SOZ and could be an effective tool for mapping epileptic networks.