The Virtual Child Brain: Modeling Neuromaturational Trajectories

Introduction Understanding the mechanisms of human neuromaturation constitutes one of the fundamental questions of neuroscience. While it is well described that large-scale brain maturation is initiated within sensorimotor brain regions and progresses to associative cortex, the underlying developmental neurobiology remains to be fully characterized. Animal models have indicated that cortical inhibitory upregulation might be a driver of neurodevelopment. To investigate the hypothesis that cortical inhibitory upregulation plays a similar role in human neuromaturation, we developed a The Virtual Brain (TVB) based computational model (TVB-Child) to explore potential mechanisms of human neurodevelopment. Material and method We created neurodevelopmental dynamic brain network models capturing neurobiological maturation by using the large-scale brain simulator TVB and fitting brain network models to developmental functional MRI (fMRI) from the Human Connectome Project-Development (HCP-D) data set with 640 subjects with an age range of 6-21 years. Age-dependent trajectories in the fMRI data set were first analyzed by combined group-ICA/Dual Regression extracting subject-specific resting-state networks (RSN). Maturational topographical and topological redistribution of these networks were analyzed by linear and non-linear regression of RSN size and degree and strength centrality. Brain network models were fitted to the fMRI functional connectivity obtained from the HCP-D data set. Hypothesizing that cortical inhibition is a driver of neuromaturation, we analyzed spatiotemporal inhibition parameter gradients in the dynamic brain network model for the hypothesized significant correlations with fMRI RSN maturational trajectories. Results While during development frontoparietal (FP) and default mode network (DMN) grew and exhibited an increase in both degree and strength centrality, becoming dominant network hubs, the attention network underwent network pruning with a decrease in size and node degree. The primary sensory network changed little. For the fitted brain network models, we obtained a high degree of reproduction with correlation coefficients between empirical and simulated functional connectivities ranging between 0.80 and 0.95. Values of the feed forward inhibition model parameter wijFFI representing the strength of regional feedforward inhibitory input exhibited the most significant increase with age within the FP and DMN networks. A less pronounced, but significant, age-dependent increase of the inhibitory parameter values were seen in attention networks and no change within primary sensory networks. Conclusion Our study shows that high order (FP, DMN), attention and primary sensory networks exhibit distinct topographical and topological maturation trajectories. Moreover, brain network modeling revealed RSN-specific age-dependent inhibition trajectories, indicating that the model is able to reproduce and thus support candidate mechanisms of neurodevelopment.