Search

Abstract Across development, experience has a strong impact on the way we think and adapt. School experience affects academic and social-emotional outcomes, yet whether differences in pedagogical experience modulate underlying brain network development is still unknown. In this study, we compared the brain network dynamics of students with different pedagogical backgrounds. Specifically, we characterized the diversity and stability of brain activity at rest by combining both resting-state fMRI and diffusion-weighted structural imaging data of 87 4–18 years old students experiencing either the Montessori pedagogy (i.e., student-led, trial-and-error pedagogy) or the traditional pedagogy (i.e., teacher-led, test-based pedagogy). Our results revealed spatiotemporal brain dynamics differences between students as a function of schooling experience at the whole-brain level. Students from Montessori schools showed overall higher functional integration (higher system diversity) and neural stability (lower spatiotemporal diversity) compared to traditionally schooled students. Higher integration was explained mainly through the cerebellar (CBL) functional network. In contrast, higher temporal stability was observed in the ventral attention, dorsal attention, somatomotor, frontoparietal, and CBL functional networks. This study suggests a form of experience-dependent dynamic functional connectivity plasticity, in learning-related networks.

Abstract While research has unveiled and quantified brain markers of abnormal neurodevelopment, clinicians still work with qualitative metrics for MRI brain investigation. The purpose of the current article is to bridge the knowledge gap between case-control cohort studies and individual patient care. Here, we provide a unique dataset of seventy-three 3-to-17 years-old healthy subjects acquired with a 6-minute MRI protocol encompassing T1 and T2 relaxation quantitative sequence that can be readily implemented in the clinical setting; MP2RAGE for T1 mapping and the prototype sequence GRAPPATINI for T2 mapping. White matter and grey matter volumes were automatically quantified. We further provide normative developmental curves based on these two imaging sequences; T1, T2 and volume normative ranges with respect to age were computed, for each ROI of a pediatric brain atlas. This open-source dataset provides normative values allowing to position individual patients acquired with the same protocol on the brain maturation curve and as such provides potentially useful quantitative biomarkers facilitating precise and personalized care.

Brain oscillations are produced by the coordinated activity of large groups of neurons and different rhythms are thought to reflect different modes of information processing. These modes, in turn, are known to occur at different spatial scales. Nevertheless, how these rhythms support different spatial modes of information processing at the brain scale is not yet fully understood. Here we use ''Joint Time-Vertex Spectral Analysis'' to characterize the joint spectral content of brain activity both in time (temporal frequencies) and in space over the connectivity graph (spatial connectome harmonics). This method allows us to characterize the relationship between spatially localized or distributed neural processes on one side and their respective temporal frequency bands in source-reconstructed M/EEG signals. We explore this approach on two different datasets, an auditory steady-state response (ASSR) and a visual grating task. Our results suggest that different information processing mechanisms are carried out at different frequency bands: while spatially distributed activity (which may also be interpreted as integration) specifically occurs at low temporal frequencies (alpha and theta) and low graph spatial frequencies, localized electrical activity (i.e., segregation) is observed at high temporal frequencies (high and low gamma) over restricted high spatial graph frequencies. Crucially, the estimated contribution of the distributed and localized neural activity predicts performance in a behavioral task, demonstrating the neurophysiological relevance of the joint time-vertex spectral representation.

Abstract How behavior arises from brain physiology has been one central topic of investigation in neuroscience. Considering the recent interest in predicting behavior from brain imaging using open datasets, there is the need for a principled approach to the categorization of behavioral variables. However, this is not trivial, as the definitions of psychological constructs and their relationships—their ontology—are not always clear. Here, we propose to use exploratory factor analysis (EFA) as a data-driven approach to find robust and interpretable domains of behavior in the Human Connectome Project (HCP) dataset. Additionally, we explore the clustering of behavioral variables using consensus clustering. We find that four and five factors offer the best description of the data, a result corroborated by the consensus clustering. In the four-factor solution, factors for Mental Health, Cognition, Processing Speed, and Substance Use arise. With five factors, Mental Health splits into Well-Being and Internalizing. Clustering results show a similar pattern, with clusters for Cognition, Processing Speed, Positive Affect, Negative Affect, and Substance Use. The factor structure is replicated in an independent dataset using confirmatory factor analysis (CFA). We discuss how the content of the factors fits with previous conceptualizations of general behavioral domains.

Measurement(s) brain connectivity measurement Technology Type(s) Fiber tracking

Abstract Accurate characterization of in utero human brain maturation is critical as it involves complex and interconnected structural and functional processes that may influence health later in life. Magnetic resonance imaging is a powerful tool to investigate equivocal neurological patterns during fetal development. However, the number of acquisitions of satisfactory quality available in this cohort of sensitive subjects remains scarce, thus hindering the validation of advanced image processing techniques. Numerical phantoms can mitigate these limitations by providing a controlled environment with a known ground truth. In this work, we present FaBiAN, an open-source Fetal Brain magnetic resonance Acquisition Numerical phantom that simulates clinical T2-weighted fast spin echo sequences of the fetal brain. This unique tool is based on a general, flexible and realistic setup that includes stochastic fetal movements, thus providing images of the fetal brain throughout maturation comparable to clinical acquisitions. We demonstrate its value to evaluate the robustness and optimize the accuracy of an algorithm for super-resolution fetal brain magnetic resonance imaging from simulated motion-corrupted 2D low-resolution series compared to a synthetic high-resolution reference volume. We also show that the images generated can complement clinical datasets to support data-intensive deep learning methods for fetal brain tissue segmentation.

Measurement(s) visual evoked potentials • brain activity measurement • diffusion weighted images Technology Type(s) high-density electroencephalography • magnetic resonance imaging Factor Type(s) visual stimulus Sample Characteristic - Organism Homo sapiens Sample Characteristic - Environment laboratory environment Machine-accessible metadata file describing the reported data: https://doi.org/10.6084/m9.figshare.17074685

AbstractThe dynamic repertoire of functional brain networks is constrained by the underlying topology of structural connections. Despite this intrinsic relationship between structural connectivity (SC) and functional connectivity (FC), integrative and multimodal approaches to combine the two remain limited. Here, we propose a new adaptive filter for estimating dynamic and directed FC using structural connectivity information as priors. We tested the filter in rat epicranial recordings and human event-related EEG data, using SC priors from a meta-analysis of tracer studies and diffusion tensor imaging metrics, respectively. We show that, particularly under conditions of low signal-to-noise ratio, SC priors can help to refine estimates of directed FC, promoting sparse functional networks that combine information from structure and function. In addition, the proposed filter provides intrinsic protection against SC-related false negatives, as well as robustness against false positives, representing a valuable new tool for multimodal imaging in the context of dynamic and directed FC analysis.

Aim: Pathological states of recovery after coma as a result of a severe brain injury are marked with changes in structural connectivity of the brain. This study aimed to identify a topological correlation between white matter integrity and the level of functional and cognitive impairment in patients recovering after coma. Methods: Structural connectomes were computed based on fractional anisotropy maps from 40 patients using a probabilistic human connectome atlas. We used a network based statistics approach to identify potential brain networks associated with a more favorable outcome, assessed with clinical neurobehavioral scores at the patient’s discharge from the acute neurorehabilitation unit. Results: We identified a subnetwork whose strength of connectivity correlated with a more favorable outcome as measured with the Disability Rating Scale (network based statistics: t >3.5, P =.010). The subnetwork predominated in the left hemisphere and included the thalamic nuclei, putamen, precentral and postcentral gyri, and medial parietal regions. Spearman correlation between the mean fractional anisotropy value of the subnetwork and the score was ρ = −0.60 (P 3.5, P =.033; Spearman’s ρ = 0.58, P

The implementation of adequate quality assessment (QA) and quality control (QC) protocols within the magnetic resonance imaging (MRI) research workflow is resource- and time-consuming and even more so is their execution. As a result, QA/QC practices highly vary across laboratories and “MRI schools”, ranging from highly specialized knowledge spots to environments where QA/QC is considered overly onerous and costly despite evidence showing that below-standard data increase the false positive and false negative rates of the final results. Here, we demonstrate a protocol based on the visual assessment of images one-by-one with reports generated by MRIQC and fMRIPrep, for the QC of data in functional (blood-oxygen dependent-level; BOLD) MRI analyses. We particularize the proposed, open-ended scope of application to whole-brain voxel-wise analyses of BOLD to correspondingly enumerate and define the exclusion criteria applied at the QC checkpoints. We apply our protocol on a composite dataset (n = 181 subjects) drawn from open fMRI studies, resulting in the exclusion of 97% of the data (176 subjects). This high exclusion rate was expected because subjects were selected to showcase artifacts. We describe the artifacts and defects more commonly found in the dataset that justified exclusion. We moreover release all the materials we generated in this assessment and document all the QC decisions with the expectation of contributing to the standardization of these procedures and engaging in the discussion of QA/QC by the community.

Abstract Although shared behavioral and neural mechanisms between working memory (WM) and motor sequence learning (MSL) have been suggested, the additive and interactive effects of training have not been studied. This study aimed at investigating changes in brain functional connectivity (FC) induced by sequential (WM + MSL and MSL + WM) and combined (WM × MSL) training programs. 54 healthy subjects (27 women; mean age: 30.2 ± 8.6 years) allocated to three training groups underwent twenty-four 40-min training sessions over 6 weeks and four cognitive assessments including functional MRI. A double-baseline approach was applied to account for practice effects. Test performances were compared using linear mixed-effects models and t-tests. Resting state fMRI data were analysed using FSL. Processing speed, verbal WM and manual dexterity increased following training in all groups. MSL + WM training led to additive effects in processing speed and verbal WM. Increased FC was found after training in a network including the right angular gyrus, left superior temporal sulcus, right superior parietal gyrus, bilateral middle temporal gyri and left precentral gyrus. No difference in FC was found between double baselines. Results indicate distinct patterns of resting state FC modulation related to sequential and combined WM and MSL training suggesting a relevance of the order of training performance. These observations could provide new insight for the planning of effective training/rehabilitation.

AbstractThe wiring of the brain is organized around a putative unimodal-transmodal hierarchy. Here we investigate how this intrinsic hierarchical organization of the brain shapes the transmission of information among regions. The hierarchical positioning of individual regions was quantified by applying diffusion map embedding to resting-state functional MRI networks. Structural networks were reconstructed from diffusion spectrum imaging and topological shortest paths among all brain regions were computed. Sequences of nodes encountered along a path were then labeled by their hierarchical position, tracing out path motifs. We find that the cortical hierarchy guides communication in the network. Specifically, nodes are more likely to forward signals to nodes closer in the hierarchy and cover a range of unimodal and transmodal regions, potentially enriching or diversifying signals en route. We also find evidence of systematic detours, particularly in attention networks, where communication is rerouted. Altogether, the present work highlights how the cortical hierarchy shapes signal exchange and imparts behaviorally relevant communication patterns in brain networks.

AbstractRecently, EEG recording techniques and source analysis have improved, making it feasible to tap into fast network dynamics. Yet, analyzing whole-cortex EEG signals in source space is not standard, partly because EEG suffers from volume conduction: Functional connectivity (FC) reflecting genuine functional relationships is impossible to disentangle from spurious FC introduced by volume conduction. Here, we investigate the relationship between white matter structural connectivity (SC) and large-scale network structure encoded in EEG-FC. We start by confirming that FC (power envelope correlations) is predicted by SC beyond the impact of Euclidean distance, in line with the assumption that SC mediates genuine FC. We then use information from white matter structural connectivity in order to smooth the EEG signal in the space spanned by graphs derived from SC. Thereby, FC between nearby, structurally connected brain regions increases while FC between nonconnected regions remains unchanged, resulting in an increase in genuine, SC-mediated FC. We analyze the induced changes in FC, assessing the resemblance between EEG-FC and volume-conduction- free fMRI-FC, and find that smoothing increases resemblance in terms of overall correlation and community structure. This result suggests that our method boosts genuine FC, an outcome that is of interest for many EEG network neuroscience questions.

The functional organization of neural processes is constrained by the brain's intrinsic structural connectivity, i.e., the connectome. Here, we explore how structural connectivity can improve the representation of brain activity signals and their dynamics. Using a multi-modal imaging dataset (electroencephalography, structural MRI, and diffusion MRI), we represent electrical brain activity at the cortical surface as a time-varying composition of harmonic modes of structural connectivity. These harmonic modes are known as connectome harmonics. Here we describe brain activity signal as a time-varying combination of connectome harmonics. We term this description as the connectome spectrum of the signal. We found that: first, the brain activity signal is represented more compactly by the connectome spectrum than by the traditional area-based representation; second, the connectome spectrum characterizes fast brain dynamics in terms of signal broadcasting profile, revealing different temporal regimes of integration and segregation that are consistent across participants. And last, the connectome spectrum characterizes fast brain dynamics with fewer degrees of freedom than area-based signal representations. Specifically, we show that a smaller number of dimensions capture the differences between low-level and high-level visual processing in the connectome spectrum. Also, we demonstrate that connectome harmonics capture more sensitively the topological properties of brain activity. In summary, this work provides statistical, functional, and topological evidence indicating that the description of brain activity in terms of structural connectivity fosters a more comprehensive understanding of large-scale dynamic neural functioning.

Signaling in brain networks unfolds over multiple topological scales. Areas may exchange information over local circuits, encompassing direct neighbours and areas with similar functions, or over global circuits, encompassing distant neighbours with dissimilar functions. Here we study how the organization of cortico-cortical networks mediate localized and global communication by parametrically tuning the range at which signals are transmitted on the white matter connectome. We show that brain regions vary in their preferred communication scale. By investigating the propensity for brain areas to communicate with their neighbors across multiple scales, we naturally reveal their functional diversity: unimodal regions show preference for local communication and multimodal regions show preferences for global communication. We show that these preferences manifest as region- and scale-specific structure-function coupling. Namely, the functional connectivity of unimodal regions emerges from monosynaptic communication in small-scale circuits, while the functional connectivity of transmodal regions emerges from polysynaptic communication in large-scale circuits. Altogether, the present findings reveal that communication preferences are highly heterogeneous across the cortex, shaping regional differences in structure-function coupling.

We present the comparison of two-dimensional (2D) fetal brain biometry on magnetic resonance (MR) images using orthogonal 2D T2-weighted sequences (T2WSs) vs. one 3D super-resolution (SR) reconstructed volume and evaluation of the level of confidence and concordance between an experienced pediatric radiologist (obs1) and a junior radiologist (obs2). Twenty-five normal fetal brain MRI scans (18–34 weeks of gestation) including orthogonal 3-mm-thick T2WSs were analyzed retrospectively. One 3D SR volume was reconstructed per subject based on multiple series of T2WSs. The two observers performed 11 2D biometric measurements (specifying their level of confidence) on T2WS and SR volumes. Measurements were compared using the paired Wilcoxon rank sum test between observers for each dataset (T2WS and SR) and between T2WS and SR for each observer. Bland–Altman plots were used to assess the agreement between each pair of measurements. Measurements were made with low confidence in three subjects by obs1 and in 11 subjects by obs2 (mostly concerning the length of the corpus callosum on T2WS). Inter-rater intra-dataset comparisons showed no significant difference (p > 0.05), except for brain axial biparietal diameter (BIP) on T2WS and for brain and skull coronal BIP and coronal transverse cerebellar diameter (DTC) on SR. None of them remained significant after correction for multiple comparisons. Inter-dataset intra-rater comparisons showed statistical differences in brain axial and coronal BIP for both observers, skull coronal BIP for obs1, and axial and coronal DTC for obs2. After correction for multiple comparisons, only axial brain BIP remained significantly different, but differences were small (2.95 ± 1.73 mm). SR allows similar fetal brain biometry as compared to using the conventional T2WS while improving the level of confidence in the measurements and using a single reconstructed volume.

Summary: The human brain consists of specialized areas that flexibly interact to form a multitude of functional networks. Complementary to this notion of modular organization, brain function has been shown to vary along a smooth continuum across the whole cortex. We demonstrate a mathematical framework that accounts for both of these perspectives: harmonic modes. We calculate the harmonic modes of the brain’s functional connectivity graph, called “functional harmonics,” revealing a multi-dimensional, frequency-ordered set of basis functions. Functional harmonics link characteristics of cortical organization across several spatial scales, capturing aspects of intra-areal organizational features (retinotopy, somatotopy), delineating brain areas, and explaining macroscopic functional networks as well as global cortical gradients. Furthermore, we show how the activity patterns elicited by seven different tasks are reconstructed from a very small subset of functional harmonics. Our results suggest that the principle of harmonicity, ubiquitous in nature, also underlies functional cortical organization in the human brain.

Large-scale structural brain networks encode white matter connectivity patterns among distributed brain areas. These connection patterns are believed to support cognitive processes and, when compromised, can lead to neurocognitive deficits and maladaptive behavior. A powerful approach for studying the organizing principles of brain networks is to construct group-representative networks from multisubject cohorts. Doing so amplifies signal to noise ratios and provides a clearer picture of brain network organization. Here, we show that current approaches for generating sparse group-representative networks overestimate the proportion of short-range connections present in a network and, as a result, fail to match subject-level networks along a wide range of network statistics. We present an alternative approach that preserves the connection-length distribution of individual subjects. We have used this method in previous papers to generate group-representative networks, though to date its performance has not been appropriately benchmarked and compared against other methods. As a result of this simple modification, the networks generated using this approach successfully recapitulate subject-level properties, outperforming similar approaches by better preserving features that promote integrative brain function rather than segregative. The method developed here holds promise for future studies investigating basic organizational principles and features of large-scale structural brain networks. Sparse structural connectivity data from many subjects can be succinctly represented using appropriate averaging procedures. We show, however, that several popular procedures for doing so generate group-averaged networks with statistics that are dissimilar from the subject-level networks they are intended to represent. These dissimilarities, we argue, arise from the over- and underexpression of short-range and long-distance connections, respectively, in the group-averaged matrix. We present a distance-dependent thresholding procedure that preserves connection length distributions and consequently better matches subject-level networks and their statistics. These findings inform data-driven exploratory analyses of connectomes.

An accurate evaluation and detection of awareness after a severe brain injury is crucial to a patient’s diagnosis, therapy, and end-of-life decisions. Misdiagnosis is frequent as behavior-based assessments often overlook subtle signs of consciousness. This study aimed to identify brain MRI characteristics of patients with residual consciousness after a severe brain injury and to develop a simple MRI-based scoring system according to the findings.We retrieved data from 128 patients and split them into a development or validation set. Structural brain MRIs were qualitatively assessed for lesions in 18 brain regions. We used logistic regression and support vector machine algorithms to first identify the most relevant brain regions predicting a patient’s outcome in the development set. We next built a diagnostic MRI-based score and estimated its optimal diagnostic cut-off point. The classifiers were then tested on the validation set and their performance compared using the receiver operating characteristic curve.Relevant brain regions predicting negative outcome highly overlapped between both classifiers and included the left mesencephalon, right basal ganglia, right thalamus, right parietal cortex, and left frontal cortex. The support vector machine classifier showed higher accuracy (0.93, 95% CI: 0.81–0.96) and specificity (0.97, 95% CI: 0.85–1) than logistic regression (accuracy: 0.87, 95% CI: 0.73 – 0.95; specificity: 0.90, 95% CI: 0.75–0.97), but equal sensitivity (0.67, 95% CI: 0.24–0.94 and 0.22–0.96, respectively) for distinguishing patients with and without residual consciousness.The novel MRI-based score assessing brain lesions in patients with disorders of consciousness accurately detects patients with residual consciousness. It could complement valuably behavioral evaluation as it is time-efficient and requires only conventional MRI.

We present an approach for tracking fast spatiotemporal cortical dynamics in which we combine white matter connectivity data with source-projected electroencephalographic (EEG) data. We employ the mathematical framework of graph signal processing in order to derive the Fourier modes of the brain structural connectivity graph, or “network harmonics”. These network harmonics are naturally ordered by smoothness. Smoothness in this context can be understood as the amount of variation along the cortex, leading to a multi-scale representation of brain connectivity. We demonstrate that network harmonics provide a sparse representation of the EEG signal, where, at certain times, the smoothest 15 network harmonics capture 90% of the signal power. This suggests that network harmonics are functionally meaningful, which we demonstrate by using them as a basis for the functional EEG data recorded from a face detection task. There, only 13 network harmonics are sufficient to track the large-scale cortical activity during the processing of the stimuli with a 50 ​ms resolution, reproducing well-known activity in the fusiform face area as well as revealing co-activation patterns in somatosensory/motor and frontal cortices that an unconstrained ROI-by-ROI analysis fails to capture. The proposed approach is simple and fast, provides a means of integration of multimodal datasets, and is tied to a theoretical framework in mathematics and physics. Thus, network harmonics point towards promising research directions both theoretically - for example in exploring the relationship between structure and function in the brain - and practically - for example for network tracking in different tasks and groups of individuals, such as patients.

A brief session of rightward prismatic adaptation (R-PA) has been shown to alleviate neglect symptoms in patients with right hemispheric damage, very likely by switching hemispheric dominance of the ventral attentional network (VAN) from the right to the left and by changing task-related activity within the dorsal attentional network (DAN). We have investigated this very rapid change in functional organisation with a network approach by comparing resting-state connectivity before and after a brief exposure i) to R-PA (14 normal subjects; experimental condition) or ii) to plain glasses (12 normal subjects; control condition). A whole brain analysis (comprising 129 regions of interest) highlighted R-PA-induced changes within a bilateral, fronto-temporal network, which consisted of 13 nodes and 11 edges; all edges involved one of 4 frontal nodes, which were part of VAN. The analysis of network characteristics within VAN and DAN revealed a R-PA-induced decrease in connectivity strength between nodes and a decrease in local efficiency within VAN but not within DAN. These results indicate that the resting-state connectivity configuration of VAN is modulated by R-PA, possibly by decreasing its modularity.

Schizophrenia, as a psychiatric disorder, has recognized brain alterations both at the structural and at the functional magnetic resonance imaging level. The developing field of connectomics has attracted much attention as it allows researchers to take advantage of powerful tools of network analysis in order to study structural and functional connectivity abnormalities in schizophrenia. Many methods have been proposed to identify biomarkers in schizophrenia, focusing mainly on improving the classification performance or performing statistical comparisons between groups. However, the stability of biomarkers selection has been for long overlooked in the connectomics field. In this study, we follow a machine learning approach where the identification of biomarkers is addressed as a feature selection problem for a classification task. We perform a recursive feature elimination and support vector machines (RFE-SVM) approach to identify the most meaningful biomarkers from the structural, functional, and multi-modal connectomes of healthy controls and patients. Furthermore, the stability of the retrieved biomarkers is assessed across different subsamplings of the dataset, allowing us to identify the affected core of the pathology. Considering our technique altogether, it demonstrates a principled way to achieve both accurate and stable biomarkers while highlighting the importance of multi-modal approaches to brain pathology as they tend to reveal complementary information.

Objective: Epilepsy diagnosis can be difficult in the absence of interictal epileptic discharges (IED) on scalp EEG. We used high-density EEG to measure connectivity in large‐scale functional networks of patients with focal epilepsy (Temporal and Extratemporal Lobe Epilepsy, TLE and ETLE) and tested for network alterations during resting wakefulness without IEDs, compared to healthy controls. We measured global efficiency as a marker of integration within networks. Methods: We analysed 49 adult patients with focal epilepsy and 16 healthy subjects who underwent high-density-EEG and structural MRI. We estimated cortical activity using electric source analysis in 82 atlas-based cortical regions based on the individual MRI. We applied directed connectivity analysis (Partial Directed Coherence) on these sources and performed graph analysis: we computed the Global Efficiency on the whole brain and on each resting state network. We tested these features in different group of patients. Results: Compared to controls, efficiency was increased in both TLE and ETLE (p

The function of a brain region is determined by the network it is embedded in. Here the authors implement the word2vec algorithm for connectomes generating a vector embedding of the connectivity structure for each node allowing inference about functional relationships between brain regions.

Abstract The last two decades of network science have discovered stunning similarities in the topological characteristics of real life networks (many biological, social, transportation and organizational networks) on a strong empirical basis. However our knowledge about the operational paths used in these networks is very limited, which prohibits the proper understanding of the principles of their functioning. Today, the most widely adopted hypothesis about the structure of the operational paths is the shortest path assumption. Here we present a striking result that the paths in various networks are significantly stretched compared to their shortest counterparts. Stretch distributions are also found to be extremely similar. This phenomenon is empirically confirmed on four networks from diverse areas of life. We also identify the high-level path selection rules nature seems to use when picking its paths.

Communication of signals among nodes in a complex network poses fundamental problems of efficiency and cost. Routing of messages along shortest paths requires global information about the topology, while spreading by diffusion, which operates according to local topological features, is informationally "cheap" but inefficient. We introduce a stochastic model for network communication that combines local and global information about the network topology to generate biased random walks on the network. The model generates a continuous spectrum of dynamics that converge onto shortest-path and random-walk (diffusion) communication processes at the limiting extremes. We implement the model on two cohorts of human connectome networks and investigate the effects of varying the global information bias on the network's communication cost. We identify routing strategies that approach a (highly efficient) shortest-path communication process with a relatively small global information bias on the system's dynamics. Moreover, we show that the cost of routing messages from and to hub nodes varies as a function of the global information bias driving the system's dynamics. Finally, we implement the model to identify individual subject differences from a communication dynamics point of view. The present framework departs from the classical shortest paths vs. diffusion dichotomy, unifying both models under a single family of dynamical processes that differ by the extent to which global information about the network topology influences the routing patterns of neural signals traversing the network.

Studies have shown scholastic, creative, and social benefits of Montessori education, benefits that were hypothesized to result from better executive functioning on the part of those so educated. As these previous studies have not reported consistent outcomes supporting this idea, we therefore evaluated scholastic development in a cross-sectional study of kindergarten and elementary school-age students, with an emphasis on the three core executive measures of cognitive flexibility, working memory update, and selective attention (inhibition). Two hundred and one (201) children underwent a complete assessment: half of the participants were from Montessori settings, while the other half were controls from traditional schools. The results confirmed that Montessori participants outperformed peers from traditional schools both in academic outcomes and in creativity skills across age groups and in self-reported well-being at school at kindergarten age. No differences were found in global executive functions, except working memory. Moreover, a multiple mediations model revealed a significant impact of creative skills on academic outcomes influenced by the school experience. These results shed light on the possibly overestimated contribution of executive functions as the main contributor to scholastic success of Montessori students and call for further investigation. Here, we propose that Montessori school-age children benefit instead from a more balanced development stemming from self-directed creative execution.

Chromosome 22q11.2 deletion syndrome (22q11DS) is a genetic disease known to lead to cerebral structural alterations, which we study using the framework of the macroscopic white-matter connectome. We create weighted connectomes of 44 patients with 22q11DS and 44 healthy controls using diffusion tensor magnetic resonance imaging, and perform a weighted graph theoretical analysis. After confirming global network integration deficits in 22q11DS (previously identified using binary connectomes), we identify the spatial distribution of regions responsible for global deficits. Next, we further characterize the dysconnectivity of the deficient regions in terms of sub-network properties, and investigate their relevance with respect to clinical profiles. We define the subset of regions with decreased nodal integration (evaluated using the closeness centrality measure) as the affected core (A-core) of the 22q11DS structural connectome. A-core regions are broadly bilaterally symmetric and consist of numerous network hubs — chiefly parietal and frontal cortical, as well as subcortical regions. Using a simulated lesion approach, we demonstrate that these core regions and their connections are particularly important to efficient network communication. Moreover, these regions are generally densely connected, but less so in 22q11DS. These specific disturbances are associated to a rerouting of shortest network paths that circumvent the A-core in 22q11DS, “de-centralizing” the network. Finally, the efficiency and mean connectivity strength of an orbito-frontal/cingulate circuit, included in the affected regions, correlate negatively with the extent of negative symptoms in 22q11DS patients, revealing the clinical relevance of present findings. The identified A-core overlaps numerous regions previously identified as affected in 22q11DS as well as in schizophrenia, which approximately 30–40% of 22q11DS patients develop.

Error-monitoring is a crucial cognitive process that enables us to adapt to the constantly changing environment. The anterior cingulate cortex (ACC) plays a vital role in error-monitoring, and its prolonged maturation suggests that it can be influenced by experience-dependent plasticity. To explore this possibility, we collected morphometric magnetic resonance imaging (MRI) measures of the ACC and error-related response-locked event-related potentials (ERPs) in 26 schoolchildren, aged 6-12 years, enrolled in either a Montessori or a traditional curriculum in Switzerland. We show that the caudal ACC undergoes significant morphometric changes during this developmental age range that seem related to error detection ERP activity. Furthermore, we observed differences in source localization activity related to error detection within the caudal ACC between Montessori and traditionally-schooled children, indicating a potential difference in the development of error-monitoring in these groups. Our study provides preliminary evidence for a potential window of opportunity to influence error-monitoring during development and calls for more work in that direction.

OBJECTIVES:We studied in a clinical setting the age dependent T1 relaxation time as a marker of normal late brain maturation and compared it to conventional techniques, namely the apparent diffusion coefficient (ADC). MATERIALS AND METHODS:Forty-two healthy subjects ranging from ages 1 year to 20 years were included in our study. T1 brain maps in which the intensity of each pixel corresponded to T1 relaxation times were generated based on MR imaging data acquired using a MP2RAGE sequence. During the same session, diffusion tensor imaging data was collected. T1 relaxation times and ADC in white matter and grey matter were measured in seven clinically relevant regions of interest and were correlated to subjects' age. RESULTS:In the basal ganglia, there was a small, yet significant, decrease in T1 relaxation time (-0.45 ≤R≤-0.59, p

The maturation of the cortical gray matter (GM) and white matter (WM) are described as sequential processes following multiple, but distinct rules. However, neither the mechanisms driving brain maturation processes, nor the relationship between GM and WM maturation are well understood. Here we use connectomics and two MRI measures reflecting maturation related changes in cerebral microstructure, namely the Apparent Diffusion Coefficient (ADC) and the T1 relaxation time (T1), to study brain development. We report that the advancement of GM and WM maturation are inter-related and depend on the underlying brain connectivity architecture. Particularly, GM regions and their incident WM connections show corresponding maturation levels, which is also observed for GM regions connected through a WM tract. Based on these observations, we propose a simple computational model supporting a key role for the connectome in propagating maturation signals sequentially from external stimuli, through primary sensory structures to higher order functional cortices.

How a stimulus or a task alters the spontaneous dynamics of the brain remains a fundamental open question in neuroscience. One of the most robust hallmarks of task/stimulus-driven brain dynamics is the decrease of variability with respect to the spontaneous level, an effect seen across multiple experimental conditions and in brain signals observed at different spatiotemporal scales. Recently, it was observed that the trial-to-trial variability and temporal variance of functional magnetic resonance imaging (fMRI) signals decrease in the task-driven activity. Here we examined the dynamics of a large-scale model of the human cortex to provide a mechanistic understanding of these observations. The model allows computing the statistics of synaptic activity in the spontaneous condition and in putative tasks determined by external inputs to a given subset of brain regions. We demonstrated that external inputs decrease the variance, increase the covariances, and decrease the autocovariance of synaptic activity as a consequence of single node and large-scale network dynamics. Altogether, these changes in network statistics imply a reduction of entropy, meaning that the spontaneous synaptic activity outlines a larger multidimensional activity space than does the task-driven activity. We tested this model's prediction on fMRI signals from healthy humans acquired during rest and task conditions and found a significant decrease of entropy in the stimulus-driven activity. Altogether, our study proposes a mechanism for increasing the information capacity of brain networks by enlarging the volume of possible activity configurations at rest and reliably settling into a confined stimulus-driven state to allow better transmission of stimulus-related information.

Spatial patterns of coherent activity across different brain areas have been identified during the resting-state fluctuations of the brain. However, recent studies indicate that resting-state activity is not stationary, but shows complex temporal dynamics. We were interested in the spatiotemporal dynamics of the phase interactions among resting-state fMRI BOLD signals from human subjects. We found that the global phase synchrony of the BOLD signals evolves on a characteristic ultra-slow (

Functional imaging with intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI) is demonstrated. Images were acquired at 3 Tesla using a standard Stejskal-Tanner diffusion-weighted echo-planar imaging sequence with multiple b-values. Cerebro-spinal fluid signal, which is highly incoherent, was suppressed with an inversion recovery preparation pulse. IVIM microvascular perfusion parameters were calculated according to a two-compartment (vascular and non-vascular) diffusion model. The results obtained in 8 healthy human volunteers during visual stimulation are presented. The IVIM blood flow related parameter fD* increased 170% during stimulation in the visual cortex, and 70% in the underlying white matter.

PURPOSE: To suppress the noise, by sacrificing some of the signal homogeneity for numerical stability, in uniform T1 weighted (T1w) images obtained with the magnetization prepared 2 rapid gradient echoes sequence (MP2RAGE) and to compare the clinical utility of these robust T1w images against the uniform T1w images. MATERIALS AND METHODS: 8 healthy subjects (29.0 ± 4.1 years; 6 Male), who provided written consent, underwent two scan sessions within a 24 hour period on a 7T head-only scanner. The uniform and robust T1w image volumes were calculated inline on the scanner. Two experienced radiologists qualitatively rated the images for: general image quality; 7T specific artefacts; and, local structure definition. Voxel-based and volume-based morphometry packages were used to compare the segmentation quality between the uniform and robust images. Statistical differences were evaluated by using a positive sided Wilcoxon rank test. RESULTS: The robust image suppresses background noise inside and outside the skull. The inhomogeneity introduced was ranked as mild. The robust image was significantly ranked higher than the uniform image for both observers (observer 1/2, p-value = 0.0006/0.0004). In particular, an improved delineation of the pituitary gland, cerebellar lobes was observed in the robust versus uniform T1w image. The reproducibility of the segmentation results between repeat scans improved (p-value = 0.0004) from an average volumetric difference across structures of ≈ 6.6% to ≈ 2.4% for the uniform image and robust T1w image respectively. CONCLUSIONS: The robust T1w image enables MP2RAGE to produce, clinically familiar T1w images, in addition to T1 maps, which can be readily used in uniform morphometry packages.

Measurement of microvascular perfusion with Intravoxel Incoherent Motion (IVIM) MRI is gaining interest. Yet, the physiological influences on the IVIM perfusion parameters ("pseudo-diffusion" coefficient D*, perfusion fraction f, and flow related parameter fD*) remain insufficiently characterized. In this article, we hypothesize that D* and fD*, which depend on blood speed, should vary during the cardiac cycle. We extended the IVIM model to include time dependence of D* = D*(t), and demonstrate in the healthy human brain that both parameters D* and fD* are significantly larger during systole than diastole, while the diffusion coefficient D and f do not vary significantly. The results non-invasively demonstrate the pulsatility of the brain's microvasculature.

Diffusion MRI has evolved towards an important clinical diagnostic and research tool. Though clinical routine is using mainly diffusion weighted and tensor imaging approaches, Q-ball imaging and diffusion spectrum imaging techniques have become more widely available. They are frequently used in research-oriented investigations in particular those aiming at measuring brain network connectivity. In this work, we aim at assessing the dependency of connectivity measurements on various diffusion encoding schemes in combination with appropriate data modeling. We process and compare the structural connection matrices computed from several diffusion encoding schemes, including diffusion tensor imaging, q-ball imaging and high angular resolution schemes, such as diffusion spectrum imaging with a publically available processing pipeline for data reconstruction, tracking and visualization of diffusion MR imaging. The results indicate that the high angular resolution schemes maximize the number of obtained connections when applying identical processing strategies to the different diffusion schemes. Compared to the conventional diffusion tensor imaging, the added connectivity is mainly found for pathways in the 50-100mm range, corresponding to neighboring association fibers and long-range associative, striatal and commissural fiber pathways. The analysis of the major associative fiber tracts of the brain reveals striking differences between the applied diffusion schemes. More complex data modeling techniques (beyond tensor model) are recommended 1) if the tracts of interest run through large fiber crossings such as the centrum semi-ovale, or 2) if non-dominant fiber populations, e.g. the neighboring association fibers are the subject of investigation. An important finding of the study is that since the ground truth sensitivity and specificity is not known, the comparability between results arising from different strategies in data reconstruction and/or tracking becomes implausible to understand.

Researchers working in the field of global connectivity analysis using diffusion magnetic resonance imaging (MRI) can count on a wide selection of software packages for processing their data, with methods ranging from the reconstruction of the local intra-voxel axonal structure to the estimation of the trajectories of the underlying fibre tracts. However, each package is generally task-specific and uses its own conventions and file formats. In this article we present the Connectome Mapper, a software pipeline aimed at helping researchers through the tedious process of organising, processing and analysing diffusion MRI data to perform global brain connectivity analyses. Our pipeline is written in Python and is freely available as open-source at www.cmtk.org.

We study an adaptive statistical approach to analyze brain networks represented by brain connection matrices of interregional connectivity (connectomes). Our approach is at a middle level between a global analysis and single connections analysis by considering subnetworks of the global brain network. These subnetworks represent either the inter-connectivity between two brain anatomical regions or by the intra-connectivity within the same brain anatomical region. An appropriate summary statistic, that characterizes a meaningful feature of the subnetwork, is evaluated. Based on this summary statistic, a statistical test is performed to derive the corresponding p-value. The reformulation of the problem in this way reduces the number of statistical tests in an orderly fashion based on our understanding of the problem. Considering the global testing problem, the p-values are corrected to control the rate of false discoveries. Finally, the procedure is followed by a local investigation within the significant subnetworks. We contrast this strategy with the one based on the individual measures in terms of power. We show that this strategy has a great potential, in particular in cases where the subnetworks are well defined and the summary statistics are properly chosen. As an application example, we compare structural brain connection matrices of two groups of subjects with a 22q11.2 deletion syndrome, distinguished by their IQ scores.

Lesions of anatomical brain networks result in functional disturbances of brain systems and behavior which depend sensitively, often unpredictably, on the lesion site. The availability of whole-brain maps of structural connections within the human cerebrum and our increased understanding of the physiology and large-scale dynamics of cortical networks allow us to investigate the functional consequences of focal brain lesions in a computational model. We simulate the dynamic effects of lesions placed in different regions of the cerebral cortex by recording changes in the pattern of endogenous ("resting-state") neural activity. We find that lesions produce specific patterns of altered functional connectivity among distant regions of cortex, often affecting both cortical hemispheres. The magnitude of these dynamic effects depends on the lesion location and is partly predicted by structural network properties of the lesion site. In the model, lesions along the cortical midline and in the vicinity of the temporo-parietal junction result in large and widely distributed changes in functional connectivity, while lesions of primary sensory or motor regions remain more localized. The model suggests that dynamic lesion effects can be predicted on the basis of specific network measures of structural brain networks and that these effects may be related to known behavioral and cognitive consequences of brain lesions.

Structurally segregated and functionally specialized regions of the human cerebral cortex are interconnected by a dense network of cortico-cortical axonal pathways. By using diffusion spectrum imaging, we noninvasively mapped these pathways within and across cortical hemispheres in individual human participants. An analysis of the resulting large-scale structural brain networks reveals a structural core within posterior medial and parietal cerebral cortex, as well as several distinct temporal and frontal modules. Brain regions within the structural core share high degree, strength, and betweenness centrality, and they constitute connector hubs that link all major structural modules. The structural core contains brain regions that form the posterior components of the human default network. Looking both within and outside of core regions, we observed a substantial correspondence between structural connectivity and resting-state functional connectivity measured in the same participants. The spatial and topological centrality of the core within cortex suggests an important role in functional integration.

BackgroundSince the emergence of diffusion tensor imaging, a lot of work has been done to better understand the properties of diffusion MRI tractography. However, the validation of the reconstructed fiber connections remains problematic in many respects. For example, it is difficult to assess whether a connection is the result of the diffusion coherence contrast itself or the simple result of other uncontrolled parameters like for example: noise, brain geometry and algorithmic characteristics.Methodology/principal findingsIn this work, we propose a method to estimate the respective contributions of diffusion coherence versus other effects to a tractography result by comparing data sets with and without diffusion coherence contrast. We use this methodology to assign a confidence level to every gray matter to gray matter connection and add this new information directly in the connectivity matrix.Conclusions/significanceOur results demonstrate that whereas we can have a strong confidence in mid- and long-range connections obtained by a tractography experiment, it is difficult to distinguish between short connections traced due to diffusion coherence contrast from those produced by chance due to the other uncontrolled factors of the tractography methodology.

Understanding the large-scale structural network formed by neurons is a major challenge in system neuroscience. A detailed connectivity map covering the entire brain would therefore be of great value. Based on diffusion MRI, we propose an efficient methodology to generate large, comprehensive and individual white matter connectional datasets of the living or dead, human or animal brain. This non-invasive tool enables us to study the basic and potentially complex network properties of the entire brain. For two human subjects we find that their individual brain networks have an exponential node degree distribution and that their global organization is in the form of a small world.

The development of error monitoring is central to learning and academic achievement. However, few studies exist on the neural correlates of children's error monitoring, and no studies have examined its susceptibility to educational influences. Pedagogical methods differ on how they teach children to learn from errors. Here, 32 students (aged 8-12 years) from high-quality Swiss traditional or Montessori schools performed a math task with feedback during fMRI. Although the groups' accuracies were similar, Montessori students skipped fewer trials, responded faster and showed more neural activity in right parietal and frontal regions involved in math processing. While traditionally-schooled students showed greater functional connectivity between the ACC, involved in error monitoring, and hippocampus following correct trials, Montessori students showed greater functional connectivity between the ACC and frontal regions following incorrect trials. The findings suggest that pedagogical experience influences the development of error monitoring and its neural correlates, with implications for neurodevelopment and education.