Development of neonatal connectome dynamics and its prediction for cognitive and language outcomes at age 2.

The functional brain connectome is highly dynamic over time. However, how brain connectome dynamics evolves during the third trimester of pregnancy and is associated with later cognitive growth remains unknown. Here, we use resting-state functional Magnetic Resonance Imaging (MRI) data from 39 newborns aged 32 to 42 postmenstrual weeks to investigate the maturation process of connectome dynamics and its role in predicting neurocognitive outcomes at 2 years of age. Neonatal brain dynamics is assessed using a multilayer network model. Network dynamics decreases globally but increases in both modularity and diversity with development. Regionally, module switching decreases with development primarily in the lateral precentral gyrus, medial temporal lobe, and subcortical areas, with a higher growth rate in primary regions than in association regions. Support vector regression reveals that neonatal connectome dynamics is predictive of individual cognitive and language abilities at 2  years of age. Our findings highlight network-level neural substrates underlying early cognitive development.

Different computational relations in language are captured by distinct brain systems.

A critical way for humans to acquire information is through language, yet whether and how language experience drives specific neural semantic representations is still poorly understood. We considered statistical properties captured by 3 different computational principles of language (simple co-occurrence, network-(graph)-topological relations, and neural-network-vector-embedding relations) and tested the extent to which they can explain the neural patterns of semantic representations, measured by 2 functional magnetic resonance imaging experiments that shared common semantic processes. Distinct graph-topological word relations, and not simple co-occurrence or neural-network-vector-embedding relations, had unique explanatory power for the neural patterns in the anterior temporal lobe (capturing graph-common-neighbors), inferior frontal gyrus, and posterior middle/inferior temporal gyrus (capturing graph-shortest-path). These results were relatively specific to language: they were not explained by sensory-motor similarities and the same computational relations of visual objects (based on visual image database) showed effects in the visual cortex in the picture naming experiment. That is, different topological properties within language and the same topological computations (common-neighbors) for language and visual inputs are captured by different brain regions. These findings reveal the specific neural semantic representations along graph-topological properties of language, highlighting the information type-specific and statistical property-specific manner of semantic representations in the human brain.

Association of aerobic glycolysis with the structural connectome reveals a benefit-risk balancing mechanism in the human brain.

Aerobic glycolysis (AG), that is, the nonoxidative metabolism of glucose, contributes significantly to anabolic pathways, rapid energy generation, task-induced activity, and neuroprotection; yet high AG is also associated with pathological hallmarks such as amyloid-ß deposition. An important yet unresolved question is whether and how the metabolic benefits and risks of brain AG is structurally shaped by connectome wiring. Using positron emission tomography and magnetic resonance imaging techniques as well as computational models, we investigate the relationship between brain AG and the macroscopic connectome. Specifically, we propose a weighted regional distance-dependent model to estimate the total axonal projection length of a brain node. This model has been validated in a macaque connectome derived from tract-tracing data and shows a high correspondence between experimental and estimated axonal lengths. When applying this model to the human connectome, we find significant associations between the estimated total axonal projection length and AG across brain nodes, with higher levels primarily located in the default-mode and prefrontal regions. Moreover, brain AG significantly mediates the relationship between the structural and functional connectomes. Using a wiring optimization model, we find that the estimated total axonal projection length in these high-AG regions exhibits a high extent of wiring optimization. If these high-AG regions are randomly rewired, their total axonal length and vulnerability risk would substantially increase. Together, our results suggest that high-AG regions have expensive but still optimized wiring cost to fulfill metabolic requirements and simultaneously reduce vulnerability risk, thus revealing a benefit-risk balancing mechanism in the human brain.

3D-MASNet: 3D mixed-scale asymmetric convolutional segmentation network for 6-month-old infant brain MR images.

Precise segmentation of infant brain magnetic resonance (MR) images into gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) are essential for studying neuroanatomical hallmarks of early brain development. However, for 6-month-old infants, the extremely low-intensity contrast caused by inherent myelination hinders accurate tissue segmentation. Existing convolutional neural networks (CNNs) based segmentation models for this task generally employ single-scale symmetric convolutions, which are inefficient for encoding the isointense tissue boundaries in baby brain images. Here, we propose a 3D mixed-scale asymmetric convolutional segmentation network (3D-MASNet) framework for brain MR images of 6-month-old infants. We replaced the traditional convolutional layer of an existing to-be-trained network with a 3D mixed-scale convolution block consisting of asymmetric kernels (MixACB) during the training phase and then equivalently converted it into the original network. Five canonical CNN segmentation models were evaluated using both T1- and T2-weighted images of 23 6-month-old infants from iSeg-2019 datasets, which contained manual labels as ground truth. MixACB significantly enhanced the average accuracy of all five models and obtained the most considerable improvement in the fully convolutional network model (CC-3D-FCN) and the highest performance in the Dense U-Net model. This approach further obtained Dice coefficient accuracies of 0.931, 0.912, and 0.961 in GM, WM, and CSF, respectively, ranking first among 30 teams on the validation dataset of the iSeg-2019 Grand Challenge. Thus, the proposed 3D-MASNet can improve the accuracy of existing CNNs-based segmentation models as a plug-and-play solution that offers a promising technique for future infant brain MRI studies.

Progressive Stabilization of Brain Network Dynamics during Childhood and Adolescence.

Functional brain networks require dynamic reconfiguration to support flexible cognitive function. However, the developmental principles shaping brain network dynamics remain poorly understood. Here, we report the longitudinal development of large-scale brain network dynamics during childhood and adolescence, and its connection with gene expression profiles. Using a multilayer network model, we show the temporally varying modular architecture of child brain networks, with higher network switching primarily in the association cortex and lower switching in the primary regions. This topographical profile exhibits progressive maturation, which manifests as reduced modular dynamics, particularly in the transmodal (e.g., default-mode and frontoparietal) and sensorimotor regions. These developmental refinements mediate age-related enhancements of global network segregation and are linked with the expression profiles of genes associated with the enrichment of ion transport and nucleobase-containing compound transport. These results highlight a progressive stabilization of brain dynamics, which expand our understanding of the neural mechanisms that underlie cognitive development.

Structural insight into the individual variability architecture of the functional brain connectome.

Human cognition and behaviors depend upon the brain's functional connectomes, which vary remarkably across individuals. However, whether and how the functional connectome individual variability architecture is structurally constrained remains largely unknown. Using tractography- and morphometry-based network models, we observed the spatial convergence of structural and functional connectome individual variability, with higher variability in heteromodal association regions and lower variability in primary regions. We demonstrated that functional variability is significantly predicted by a unifying structural variability pattern and that this prediction follows a primary-to-heteromodal hierarchical axis, with higher accuracy in primary regions and lower accuracy in heteromodal regions. We further decomposed group-level connectome variability patterns into individual unique contributions and uncovered the structural-functional correspondence that is associated with individual cognitive traits. These results advance our understanding of the structural basis of individual functional variability and suggest the importance of integrating multimodal connectome signatures for individual differences in cognition and behaviors.

Individual Uniqueness in the Neonatal Functional Connectome.

The functional connectome is highly distinctive in adults and adolescents, underlying individual differences in cognition and behavior. However, it remains unknown whether the individual uniqueness of the functional connectome is present in neonates, who are far from mature. Here, we utilized the multiband resting-state functional magnetic resonance imaging data of 40 healthy neonates from the Developing Human Connectome Project and a split-half analysis approach to characterize the uniqueness of the functional connectome in the neonatal brain. Through functional connectome-based individual identification analysis, we found that all the neonates were correctly identified, with the most discriminative regions predominantly confined to the higher-order cortices (e.g., prefrontal and parietal regions). The connectivities with the highest contributions to individual uniqueness were primarily located between different functional systems, and the short- (0-30 mm) and middle-range (30-60 mm) connectivities were more distinctive than the long-range (>60 mm) connectivities. Interestingly, we found that functional data with a scanning length longer than 3.5 min were able to capture the individual uniqueness in the functional connectome. Our results highlight that individual uniqueness is present in the functional connectome of neonates and provide insights into the brain mechanisms underlying individual differences in cognition and behavior later in life.

Development of the default-mode network during childhood and adolescence: A longitudinal resting-state fMRI study.

The default-mode network (DMN) is a set of functionally connected regions that play crucial roles in internal cognitive processing. Previous resting-state fMRI studies have demonstrated that the intrinsic functional organization of the DMN undergoes remarkable reconfigurations during childhood and adolescence. However, these studies have mainly focused on cross-sectional designs with small sample sizes, limiting the consistency and interpretations of the findings. Here, we used a large sample of longitudinal resting-state fMRI data comprising 305 typically developing children (6-12 years of age at baseline, 491 scans in total) and graph theoretical approaches to delineate the developmental trajectories of the functional architecture of the DMN. For each child, the DMN was constructed according to a prior parcellation with 32 brain nodes. We showed that the overall connectivity increased in strength from childhood to adolescence and became spatially similar to that in the young adult group (N = 61, 18-28 years of age). These increases were primarily located in the midline structures. Global and local network efficiency in the DMN also increased with age, indicating an enhanced capability in parallel information communication within the brain system. Based on the divergent developmental rates of nodal centrality, we identified three subclusters within the DMN, with the fastest rates in the cluster mainly comprising the anterior medial prefrontal cortex and posterior cingulate cortex. Together, our findings highlight the developmental patterns of the functional architecture in the DMN from childhood to adolescence, which has implications for the understanding of network mechanisms underlying the cognitive development of individuals.

The spatial organization of the chronnectome associates with cortical hierarchy and transcriptional profiles in the human brain.

The chronnectome of the human brain represents dynamic connectivity patterns of brain networks among interacting regions, but its organization principle and related transcriptional signatures remain unclear. Using task-free fMRI data from the Human Connectome Project (681 participants) and microarray-based gene expression data from the Allen Institute for Brain Science (1791 brain tissue samples from six donors), we conduct a transcriptome-chronnectome association study to investigate the spatial configurations of dynamic brain networks and their linkages with transcriptional profiles. We first classify the dynamic brain networks into four categories of nodes according to their time-varying characteristics in global connectivity and modular switching: the primary sensorimotor regions with large global variations, the paralimbic/limbic regions with frequent modular switching, the frontoparietal cortex with both high global and modular dynamics, and the sensorimotor association cortex with limited dynamics. Such a spatial layout reflects the cortical functional hierarchy, microarchitecture, and primary connectivity gradient spanning from primary to transmodal areas, and the cognitive spectrum from perception to abstract processing. Importantly, the partial least squares regression analysis reveals that the transcriptional profiles could explain 28% of the variation in this spatial layout of network dynamics. The top-related genes in the transcriptional profiles are enriched for potassium ion channel complex and activity and mitochondrial part of the cellular component. These findings highlight the hierarchically spatial arrangement of dynamic brain networks and their coupling with the variation in transcriptional signatures, which provides indispensable implications for the organizational principle and cellular and molecular functions of spontaneous network dynamics.

Development and Emergence of Individual Variability in the Functional Connectivity Architecture of the Preterm Human Brain.

Individual variability in human brain networks underlies individual differences in cognition and behaviors. However, researchers have not conclusively determined when individual variability patterns of the brain networks emerge and how they develop in the early phase. Here, we employed resting-state functional MRI data and whole-brain functional connectivity analyses in 40 neonates aged around 31-42 postmenstrual weeks to characterize the spatial distribution and development modes of individual variability in the functional network architecture. We observed lower individual variability in primary sensorimotor and visual areas and higher variability in association regions at the third trimester, and these patterns are generally similar to those of adult brains. Different functional systems showed dramatic differences in the development of individual variability, with significant decreases in the sensorimotor network; decreasing trends in the visual, subcortical, and dorsal and ventral attention networks, and limited change in the default mode, frontoparietal and limbic networks. The patterns of individual variability were negatively correlated with the short- to middle-range connection strength/number and this distance constraint was significantly strengthened throughout development. Our findings highlight the development and emergence of individual variability in the functional architecture of the prenatal brain, which may lay network foundations for individual behavioral differences later in life.

Unbiased age-specific structural brain atlases for Chinese pediatric population.

In magnetic resonance (MR) imaging studies of child brain development, structural brain atlases usually serve as important references for the pediatric population, in which individual images are spatially normalized into a common or standard stereotactic space. However, the popular existing pediatric brain atlases (e.g., National Institutes of Health pediatric atlases, NIH-PD) are mostly based on MR images obtained from Caucasian populations and thus are not ideal for the characterization of the brains of Chinese children due to neuroanatomical differences related to genetic and environmental factors. Here, we use an unbiased template construction algorithm to create a set of age-specific Chinese pediatric (CHN-PD) atlases based on high-quality T1-and T2-weighted MR images from 328 cognitively normal Chinese children aged 6-12 years. The CHN-PD brain atlases include asymmetric and symmetric templates, sex-specific templates and tissue probability templates, and contain multiple age-specific templates at one-year intervals. A direct comparison of the CHN-PD and NIH-PD atlases reveals dramatic anatomical differences mainly in the bilateral frontal and parietal regions. After applying the CHN-PD and NIH-PD atlases to two independent Chinese pediatric datasets (N = 114 and N = 71), we find that the CHN-PD atlases result in significantly higher accuracy than the NIH-PD atlases in both predicting "brain age" and guiding brain tissue segmentation. These results suggest that the CHN-PD brain atlases are necessary for studies of the typical and atypical development of the Chinese pediatric population. These CHN-PD atlases have been released on the Neuroimaging Informatics Tools and Resources Clearinghouse (NITRC) website (https://www.nitrc.org/projects/chn-pd).

Chronnectome fingerprinting: Identifying individuals and predicting higher cognitive functions using dynamic brain connectivity patterns.

The human brain is a large, interacting dynamic network, and its architecture of coupling among brain regions varies across time (termed the "chronnectome"). However, very little is known about whether and how the dynamic properties of the chronnectome can characterize individual uniqueness, such as identifying individuals as a "fingerprint" of the brain. Here, we employed multiband resting-state functional magnetic resonance imaging data from the Human Connectome Project (N = 105) and a sliding time-window dynamic network analysis approach to systematically examine individual time-varying properties of the chronnectome. We revealed stable and remarkable individual variability in three dynamic characteristics of brain connectivity (i.e., strength, stability, and variability), which was mainly distributed in three higher order cognitive systems (i.e., default mode, dorsal attention, and fronto-parietal) and in two primary systems (i.e., visual and sensorimotor). Intriguingly, the spatial patterns of these dynamic characteristics of brain connectivity could successfully identify individuals with high accuracy and could further significantly predict individual higher cognitive performance (e.g., fluid intelligence and executive function), which was primarily contributed by the higher order cognitive systems. Together, our findings highlight that the chronnectome captures inherent functional dynamics of individual brain networks and provides implications for individualized characterization of health and disease.

Topological analyses of functional connectomics: A crucial role of global signal removal, brain parcellation, and null models.

Recently, functional connectome studies based on resting-state functional magnetic resonance imaging (R-fMRI) and graph theory have greatly advanced our understanding of the topological principles of healthy and diseased brains. However, how different strategies for R-fMRI data preprocessing and for connectome analyses jointly affect topological characterization and contrastive research of brain networks remains to be elucidated. Here, we used two R-fMRI data sets, a healthy young adult data set and an Alzheimer's disease (AD) patient data set, and up to 42 analysis strategies to comprehensively investigate the joint influence of three key factors (global signal regression, regional parcellation schemes, and null network models) on the topological analysis and contrastive research of whole-brain functional networks. At the global level, we first found that these three factors affected not only the quantitative values but also the individual variability profile in small-world related metrics and modularity, wherein global signal regression exhibited the predominant influence. Moreover, strategies without global signal regression and with topological randomization null model enhanced the sensitivity of the detection of differences between AD and control groups in small-worldness and modularity. At the nodal level, strategies of global signal regression dominantly influenced the spatial distribution of both hubs and between-group differences in terms of nodal degree centrality. Together, we highlight the remarkable joint influence of global signal regression, regional parcellation schemes and null network models on functional connectome analyses in both health and diseases, which may provide guidance for the choice of analysis strategies in future functional network studies.

PAGANI Toolkit: Parallel graph-theoretical analysis package for brain network big data.

The recent collection of unprecedented quantities of neuroimaging data with high spatial resolution has led to brain network big data. However, a toolkit for fast and scalable computational solutions is still lacking. Here, we developed the PArallel Graph-theoretical ANalysIs (PAGANI) Toolkit based on a hybrid central processing unit-graphics processing unit (CPU-GPU) framework with a graphical user interface to facilitate the mapping and characterization of high-resolution brain networks. Specifically, the toolkit provides flexible parameters for users to customize computations of graph metrics in brain network analyses. As an empirical example, the PAGANI Toolkit was applied to individual voxel-based brain networks with ∼200,000 nodes that were derived from a resting-state fMRI dataset of 624 healthy young adults from the Human Connectome Project. Using a personal computer, this toolbox completed all computations in ∼27 h for one subject, which is markedly less than the 118 h required with a single-thread implementation. The voxel-based functional brain networks exhibited prominent small-world characteristics and densely connected hubs, which were mainly located in the medial and lateral fronto-parietal cortices. Moreover, the female group had significantly higher modularity and nodal betweenness centrality mainly in the medial/lateral fronto-parietal and occipital cortices than the male group. Significant correlations between the intelligence quotient and nodal metrics were also observed in several frontal regions. Collectively, the PAGANI Toolkit shows high computational performance and good scalability for analyzing connectome big data and provides a friendly interface without the complicated configuration of computing environments, thereby facilitating high-resolution connectomics research in health and disease.

Intrinsic Brain Hub Connectivity Underlies Individual Differences in Spatial Working Memory.

Spatial working memory (SWM) is an important component of working memory and plays an essential role in driving high-level cognitive abilities. Recent studies have demonstrated that individual SWM is associated with global brain communication. However, whether specific network nodal connectivity, such as brain hub connectivity, is involved in individual SWM performances remains largely unknown. Here, we collected resting-state fMRI (R-fMRI) data from a large group of 130 young healthy participants and evaluated their SWM performances. A voxel-wise whole-brain network analysis approach was employed to study the relationship between the nodal functional connectivity strength (FCS) and the SWM behavioral scores and to further estimate the participation of brain hubs in individual SWM. We showed significant associations between nodal FCS and SWM performance primarily in the default mode, visual, dorsal attention, and fronto-parietal systems. Moreover, over 41% of these nodal regions were identified as brain network hubs, and these hubs' FCS values contributed to 57% of the variance of the individual SWM performances that all SWM-related regions could explain. Collectively, our findings highlight the cognitive significance of the brain network hubs in SWM, which furthers our understanding of how intrinsic brain network architectures underlie individual differences in SWM processing.

Early Development of Functional Network Segregation Revealed by Connectomic Analysis of the Preterm Human Brain.

Human brain functional networks are topologically organized with nontrivial connectivity characteristics such as small-worldness and densely linked hubs to support highly segregated and integrated information processing. However, how they emerge and change at very early developmental phases remains poorly understood. Here, we used resting-state functional MRI and voxel-based graph theory analysis to systematically investigate the topological organization of whole-brain networks in 40 infants aged around 31 to 42 postmenstrual weeks. The functional connectivity strength and heterogeneity increased significantly in primary motor, somatosensory, visual, and auditory regions, but much less in high-order default-mode and executive-control regions. The hub and rich-club structures in primary regions were already present at around 31 postmenstrual weeks and exhibited remarkable expansions with age, accompanied by increased local clustering and shortest path length, indicating a transition from a relatively random to a more organized configuration. Moreover, multivariate pattern analysis using support vector regression revealed that individual brain maturity of preterm babies could be predicted by the network connectivity patterns. Collectively, we highlighted a gradually enhanced functional network segregation manner in the third trimester, which is primarily driven by the rapid increases of functional connectivity of the primary regions, providing crucial insights into the topological development patterns prior to birth.

Individual differences and time-varying features of modular brain architecture.

Recent studies have suggested that human brain functional networks are topologically organized into functionally specialized but inter-connected modules to facilitate efficient information processing and highly flexible cognitive function. However, these studies have mainly focused on group-level network modularity analyses using "static" functional connectivity approaches. How these extraordinary modular brain structures vary across individuals and spontaneously reconfigure over time remain largely unknown. Here, we employed multiband resting-state functional MRI data (N=105) from the Human Connectome Project and a graph-based modularity analysis to systematically investigate individual variability and dynamic properties in modular brain networks. We showed that the modular structures of brain networks dramatically vary across individuals, with higher modular variability primarily in the association cortex (e.g., fronto-parietal and attention systems) and lower variability in the primary systems. Moreover, brain regions spontaneously changed their module affiliations on a temporal scale of seconds, which cannot be simply attributable to head motion and sampling error. Interestingly, the spatial pattern of intra-subject dynamic modular variability largely overlapped with that of inter-subject modular variability, both of which were highly reproducible across repeated scanning sessions. Finally, the regions with remarkable individual/temporal modular variability were closely associated with network connectors and the number of cognitive components, suggesting a potential contribution to information integration and flexible cognitive function. Collectively, our findings highlight individual modular variability and the notable dynamic characteristics in large-scale brain networks, which enhance our understanding of the neural substrates underlying individual differences in a variety of cognition and behaviors.

Identifying topological motif patterns of human brain functional networks.

Recent imaging connectome studies demonstrated that the human functional brain network follows an efficient small-world topology with cohesive functional modules and highly connected hubs. However, the functional motif patterns that represent the underlying information flow remain largely unknown. Here, we investigated motif patterns within directed human functional brain networks, which were derived from resting-state functional magnetic resonance imaging data with controlled confounding hemodynamic latencies. We found several significantly recurring motifs within the network, including the two-node reciprocal motif and five classes of three-node motifs. These recurring motifs were distributed in distinct patterns to support intra- and inter-module functional connectivity, which also promoted integration and segregation in network organization. Moreover, the significant participation of several functional hubs in the recurring motifs exhibited their critical role in global integration. Collectively, our findings highlight the basic architecture governing brain network organization and provide insight into the information flow mechanism underlying intrinsic brain activities. Hum Brain Mapp 38:2734-2750, 2017. © 2017 Wiley Periodicals, Inc.

Long-period rhythmic synchronous firing in a scale-free network.

Stimulus information is encoded in the spatial-temporal structures of external inputs to the neural system. The ability to extract the temporal information of inputs is fundamental to brain function. It has been found that the neural system can memorize temporal intervals of visual inputs in the order of seconds. Here we investigate whether the intrinsic dynamics of a large-size neural circuit alone can achieve this goal. The network models we consider have scale-free topology and the property that hub neurons are difficult to be activated. The latter is implemented by either including abundant electrical synapses between neurons or considering chemical synapses whose efficacy decreases with the connectivity of the postsynaptic neuron. We find that hub neurons trigger synchronous firing across the network, loops formed by low-degree neurons determine the rhythm of synchronous firing, and the hardness of exciting hub neurons avoids epileptic firing of the network. Our model successfully reproduces the experimentally observed rhythmic synchronous firing with long periods and supports the notion that the neural system can process temporal information through the dynamics of local circuits in a distributed way.

Development of segregation and integration of functional connectomes during the first 1,000 days.

The first 1,000 days of human life lay the foundation for brain development and later cognitive growth. However, the developmental rules of the functional connectome during this critical period remain unclear. Using high-resolution, longitudinal, task-free functional magnetic resonance imaging data from 930 scans of 665 infants aged 28 postmenstrual weeks to 3 years, we report the early maturational process of connectome segregation and integration. We show the dominant development of local connections alongside a few global connections, the shift of brain hubs from primary regions to high-order association cortices, the developmental divergence of network segregation and integration along the anterior-posterior axis, the prediction of neurocognitive outcomes, and their associations with gene expression signatures of microstructural development and neuronal metabolic pathways. These findings advance our understanding of the principles of connectome remodeling during early life and its neurobiological underpinnings and have implications for studying typical and atypical development.