Search

Recent advances in imaging and tracing technology provide increasingly detailed reconstructions of brain connectomes. Concomitant analytic advances enable rigorous identification and quantification of functionally important features of brain network architecture. Null models are a flexible tool to statistically benchmark the presence or magnitude of features of interest, by selectively preserving specific architectural properties of brain networks while systematically randomizing others. Here we describe the logic, implementation and interpretation of null models of connectomes. We introduce randomization and generative approaches to constructing null networks, and outline a taxonomy of network methods for statistical inference. We highlight the spectrum of null models - from liberal models that control few network properties, to conservative models that recapitulate multiple properties of empirical networks - that allow us to operationalize and test detailed hypotheses about the structure and function of brain networks. We review emerging scenarios for the application of null models in network neuroscience, including for spatially embedded networks, annotated networks and correlation-derived networks. Finally, we consider the limits of null models, as well as outstanding questions for the field.

Contains fulltext : 248878.pdf (Publisher’s version ) (Open Access) Current neuroimaging acquisition and processing approaches tend to be optimised for quality rather than speed. However, rapid acquisition and processing of neuroimaging data can lead to novel neuroimaging paradigms, such as adaptive acquisition, where rapidly processed data is used to inform subsequent image acquisition steps. Here we first evaluate the impact of several processing steps on the processing time and quality of registration of manually labelled T(1) -weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T(1) -FLAIR, T(2) , T(2) *, T(2) -FLAIR, DWI and ADC contrasts, acquired in ~1 min), as well as to slower, more standard single-contrast T(1) -weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix T(1) -FLAIR and single-contrast T(1) -weighted scans, using correlations between voxels and regions of interest across participants, measures of within- and between-participant identifiability as well as regional structural covariance networks. Furthermore, we explore the use of EPImix for the rapid construction of morphometric similarity networks. Finally, we quantify the reliability of EPImix-derived data using test-retest scans of 10 participants. Our results demonstrate that quantitative information can be derived from a neuroimaging scan acquired and processed within minutes, which could further be used to implement adaptive multimodal imaging and tailor neuroimaging examinations to individual patients.

Significance How does the human brain change during adolescence? We found 2 distinct modes of change in functional connectivity between brain regions, “conservative” and “disruptive,” measured using functional MRI (fMRI) in healthy young people (14 to 26 y old). Conservative regions, often specialized for basic sensory and motor functions, were strongly connected at age 14 before strengthening more by age 26, whereas disruptive regions that were activated by complex tasks comprised both connections that were weak at age 14 but strengthened by age 26 and connections that were strong at age 14 but weakened by age 26. Disruptive maturation of fMRI connectivity between cortex and subcortex could represent metabolically costly remodeling that underpins development of adult faculties., Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: “conservative” and “disruptive.” Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman’s correlation between edgewise baseline FC (at 14 y, FC14) and adolescent change in FC (ΔFC14−26), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas.

The study of complex systems deals with emergent behavior that arises as a result of nonlinear spatiotemporal interactions between a large number of components both within the system, as well as between the system and its environment. There is a strong case to be made that neural systems as well as their emergent behavior and disorders can be studied within the framework of complexity science. In particular, the field of neuroimaging has begun to apply both theoretical and experimental procedures originating in complexity science—usually in parallel with traditional methodologies. Here, we illustrate the basic properties that characterize complex systems and evaluate how they relate to what we have learned about brain structure and function from neuroimaging experiments. We then argue in favor of adopting a complex systems-based methodology in the study of neuroimaging, alongside appropriate experimental paradigms, and with minimal influences from noncomplex system approaches. Our exposition includes a review of the fundamental mathematical concepts, combined with practical examples and a compilation of results from the literature.

Current neuroimaging acquisition and processing approaches tend to be optimised for quality rather than speed. However, rapid acquisition and processing of neuroimaging data can lead to novel neuroimaging paradigms, such as adaptive acquisition, where rapidly processed data is used to inform subsequent image acquisition steps. Here we first evaluate the impact of several processing steps on the processing time and quality of registration of manually labelled T1-weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T1-FLAIR, T2, T2*, T2-FLAIR, DWI & ADC contrasts, acquired in ∼1 minute), as well as to slower, more standard single-contrast T1-weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix and single-contrast T1-weighted scans, using correlations between voxels and regions of interest across participants, measures of within- and between-participant identifiability as well as regional structural covariance networks. Furthermore, we explore the use of EPImix for the rapid construction of morphometric similarity networks. Finally, we quantify the reliability of EPImix-derived data using test-retest scans of 10 participants. Our results demonstrate that quantitative information can be derived from a neuroimaging scan acquired and processed within minutes, which could further be used to implement adaptive multimodal imaging and tailor neuroimaging examinations to individual patients.Abstract FigureGraphical abstract.

For most neuroimaging questions the range of possible analytic choices makes it unclear how to evaluate conclusions from any single analytic method. One possible way to address this issue is to evaluate all possible analyses using a multiverse approach, however, this can be computationally challenging and sequential analyses on the same data can compromise predictive power. Here, we establish how active learning on a low-dimensional space capturing the inter-relationships between pipelines can efficiently approximate the full spectrum of analyses. This approach balances the benefits of a multiverse analysis without incurring the cost on computational and predictive power. We illustrate this approach with two functional MRI datasets (predicting brain age and autism diagnosis) demonstrating how a multiverse of analyses can be efficiently navigated and mapped out using active learning. Furthermore, our presented approach not only identifies the subset of analysis techniques that are best able to predict age or classify individuals with autism spectrum disorder and healthy controls, but it also allows the relationships between analyses to be quantified.

For most neuroimaging questions the huge range of possible analytic choices leads to the possibility that conclusions from any single analytic approach may be misleading. Examples of possible choices include the motion regression approach used and smoothing and threshold factors applied during the processing pipeline. Although it is possible to perform a multiverse analysis that evaluates all possible analytic choices, this can be computationally challenging and repeated sequential analyses on the same data can compromise inferential and predictive power. Here, we establish how active learning on a low-dimensional space that captures the inter-relationships between analysis approaches can be used to efficiently approximate the whole multiverse of analyses. This approach balances the benefits of a multiverse analysis without the accompanying cost to statistical power, computational power and the integrity of inferences. We illustrate this approach with a functional MRI dataset of functional connectivity across adolescence, demonstrating how a multiverse of graph theoretic and simple pre-processing steps can be efficiently navigated using active learning. Our study shows how this approach can identify the subset of analysis techniques (i.e., pipelines) which are best able to predict participants’ ages, as well as allowing the performance of different approaches to be quantified.

The study of complex systems deals with emergent behaviour that arises as a result of nonlinear spatio-temporal interactions between a large number of components both within the system, as well as between the system and its environment. There is a strong case to be made that neural systems as well as their emergent behaviour and disorders, can be studied within the framework of complexity science. In particular, the field of neuroimaging has begun to apply both theoretical and experimental procedures originating in complexity science – usually in parallel with traditional methodologies. Here, we demonstrate that the use of such traditional models may distort the outcomes of neuroimaging experiments – hence affecting their interpretability and raising questions about their reliability.Therefore, we argue in favor of adopting a complex systems-based methodology in the study of neuroimaging, alongside appropriate experimental paradigms, and with minimal influences from non-complex systems approaches. Our exposition includes a review of the fundamental mathematical concepts, combined with practical examples and a compilation of results from the literature.

We hypothesized that there are sexual differences in human brain network development underlying the female > male divergence in adolescent depression. We tested for sex differences in parameters of brain network development (accelerated longitudinal fMRI, N=298 healthy adolescents, each scanned 1 to 3 times). Sexually divergent development of functional connectivity was located in default mode network (DMN), limbic cortex, and subcortical nuclei. Females had a more “disruptive” pattern of development, where weak functional connectivity at age 14 became stronger during adolescence. This fMRI-derived map of divergent adolescent development was co-located with (i) a map of functional dysconnectivity associated with adult major depressive disorder (MDD); and (ii) an adult brain gene expression pattern enriched for genes on the X chromosome, neurodevelopmental genes, and risk genes for MDD. Sexual divergence in disruptive development of DMN, limbic and subcortical functional networks is potentially relevant to the increased risk of depression in adolescent females.

Significance Goal-directed behavior is impaired in disorders of compulsivity. Here, we characterize the developmental trajectory of model-based control and show a progressive strengthening from adolescence to early adulthood. We found that the presence of compulsivity traits impacts on this trajectory as well as on the degree of remodeling in functional connectivity within frontostriatal circuits. These findings have implications for understanding the interplay between compulsivity, the developmental trajectory of model-based planning, and functional connectivity in frontostriatal circuits., A characteristic of adaptive behavior is its goal-directed nature. An ability to act in a goal-directed manner is progressively refined during development, but this refinement can be impacted by the emergence of psychiatric disorders. Disorders of compulsivity have been framed computationally as a deficit in model-based control, and have been linked also to abnormal frontostriatal connectivity. However, the developmental trajectory of model-based control, including an interplay between its maturation and an emergence of compulsivity, has not been characterized. Availing of a large sample of healthy adolescents (n = 569) aged 14 to 24 y, we show behaviorally that over the course of adolescence there is a within-person increase in model-based control, and this is more pronounced in younger participants. Using a bivariate latent change score model, we provide evidence that the presence of higher compulsivity traits is associated with an atypical profile of this developmental maturation in model-based control. Resting-state fMRI data from a subset of the behaviorally assessed subjects (n = 230) revealed that compulsivity is associated with a less pronounced change of within-subject developmental remodeling of functional connectivity, specifically between the striatum and a frontoparietal network. Thus, in an otherwise clinically healthy population sample, in early development, individual differences in compulsivity are linked to the developmental trajectory of model-based control and a remodeling of frontostriatal connectivity.

Adolescent changes in human brain function are not entirely understood. Here we used multi-echo functional magnetic resonance imaging (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in N=298 healthy adolescents. Participants were aged 14-26 years and were scanned on two or more occasions at least 6 months apart. We found two distinct modes of age-related change in FC: “conservative” and “disruptive”. Conservative development was characteristic of primary cortex, which was strongly connected at 14 years and became even more connected in the period 14-26 years. Disruptive development was characteristic of association cortex, hippocampus and amygdala, which were not strongly connected at 14 years but became more strongly connected during adolescence. We defined the maturational index (MI) as the signed coefficient of the linear relationship between baseline FC (at 14 years,FC14) and adolescent change in FC (∆FC14−26). Disruptive systems (with negative MI) were functionally specialised for social cognition and autobiographical memory and were significantly co-located with prior maps of aerobic glycolysis (AG), AG-related gene expression, post-natal expansion of cortical surface area, and adolescent shrinkage of cortical depth. We conclude that human brain organization is disrupted during adolescence by the emergence of strong functional connectivity of subcortical nuclei and association cortical areas, representing metabolically expensive re-modelling of synaptic connectivity between brain regions that were not strongly connected in childhood. We suggest that this re-modelling process may support emergence of social skills and self-awareness during healthy human adolescence.

Neurodevelopmental disorders are highly heritable and associated with spatially-selective disruptions of brain anatomy. The logic that translates genetic risks into spatially patterned brain vulnerabilities remains unclear but is a fundamental question in disease pathogenesis. Here, we approach this question by integrating (i)in vivoneuroimaging data from patient subgroups with known causal genomic copy number variations (CNVs), and (ii) bulk and single-cell gene expression data from healthy cortex. First, for each of six different CNV disorders, we show that spatial patterns of cortical anatomy change in youth are correlated with spatial patterns of expression for CNV region genes in bulk cortical tissue from typically-developing adults. Next, by transforming normative bulk-tissue cortical expression data into cell-type expression maps, we further link each disorder’s anatomical change map to specific cell classes and specific CNV-region genes that these cells express. Finally, we establish convergent validity of this “transcriptional vulnerability model” by inter-relating patient neuroimaging data with measures of altered gene expression in both brain and blood-derived patient tissue. Our work clarifies general biological principles that govern the mapping of genetic risks onto regional brain disruption in neurodevelopmental disorders. We present new methods that can harness these principles to screen for potential cellular and molecular determinants of disease from readily available patient neuroimaging data.

Seminal human brain histology work has demonstrated developmental waves of myelination. Here, using a micro-structural magnetic resonance imaging (MRI) marker sensitive to myelin, we studied fine-grained age differences to deduce waves of growth, stability, and decline of cortical myelination over the life-cycle. In 484 participants, aged 8-85 years, we fitted smooth growth curves to T1- to T2-weighted ratio in each of 360 regions from one of 7 cytoarchitectonic classes. From the first derivatives of these generally inverted-U trajectories, we defined three milestones: the age at peak growth; the age at onset of a stable plateau; and the age at the onset of decline. Age at peak growth had a bimodal distribution comprising an early (pre-pubertal) wave of primary sensory and motor cortices and a later (post-pubertal) wave of association, insular and limbic cortices. Most regions reached stability in the 30s but there was a second wave reaching stability in the 50s. Age at onset of decline was also bimodal: in some right hemisphere regions, the curve declined from the 60s, but in other left hemisphere regions, there was no significant decline from the stable plateau. Network analysis of the micro-structural connectome revealed that late developing (second wave) regions had significantly higher degree and other measures of centrality. These results are consistent with regionally heterogeneous waves of intracortical myelinogenesis and age-related demyelination, which may be relevant to the onset of neuropsychiatric disorders, and dementia, and the staged acquisition and decline of motor and cognitive skills over the course of the life-cycle.

How does human brain structure mature during adolescence? We used MRI to measure cortical thickness and intracortical myelination in 297 population volunteers aged 14-24 y old. We found and replicated that association cortical areas were thicker and less myelinated than primary cortical areas at 14 y. However, association cortex had faster rates of shrinkage and myelination over the course of adolescence. Age-related increases in cortical myelination were maximized approximately at the internal layer of projection neurons. Adolescent cortical myelination and shrinkage were coupled and specifically associated with a dorsoventrally patterned gene expression profile enriched for synaptic, oligodendroglial- and schizophrenia-related genes. Topologically efficient and biologically expensive hubs of the brain anatomical network had greater rates of shrinkage/myelination and were associated with overexpression of the same transcriptional profile as cortical consolidation. We conclude that normative human brain maturation involves a genetically patterned process of consolidating anatomical network hubs. We argue that developmental variation of this consolidation process may be relevant both to normal cognitive and behavioral changes and the high incidence of schizophrenia during human brain adolescence.

Functional connectomes are commonly analysed as sparse graphs, constructed by thresholding cross-correlations between regional neurophysiological signals. Thresholding generally retains the strongest edges (correlations), either by retaining edges surpassing a given absolute weight, or by constraining the edge density. The latter (more widely used) method risks inclusion of false positive edges at high edge densities and exclusion of true positive edges at low edge densities. Here we apply new wavelet-based methods, which enable construction of probabilistically-thresholded graphs controlled for type I error, to a dataset of resting-state fMRI scans of 56 patients with schizophrenia and 71 healthy controls. By thresholding connectomes to fixed edge-specific P value, we found that functional connectomes of patients with schizophrenia were more dysconnected than those of healthy controls, exhibiting a lower edge density and a higher number of (dis)connected components. Furthermore, many participants’ connectomes could not be built up to the fixed edge densities commonly studied in the literature (~5-30%), while controlling for type I error. Additionally, we showed that the topological randomisation previously reported in the schizophrenia literature is likely attributable to “non-significant” edges added when thresholding connectomes to fixed density based on correlation. Finally, by explicitly comparing connectomes thresholded by increasing P value and decreasing correlation, we showed that probabilistically thresholded connectomes show decreased randomness and increased consistency across participants. Our results have implications for future analysis of functional connectivity using graph theory, especially within datasets exhibiting heterogenous distributions of edge weights (correlations), between groups or across participants.

Macroscopic cortical networks are important for cognitive function, but it remains challenging to construct anatomically plausible individual structural connectomes from human neuroimaging. We introduce a new technique for cortical network mapping based on inter-regional similarity of multiple morphometric parameters measured using multimodal MRI. In three cohorts (two human, one macaque), we find that the resulting morphometric similarity networks (MSNs) have a complex topological organization comprising modules and high-degree hubs. Human MSN modules recapitulate known cortical cytoarchitectonic divisions, and greater inter-regional morphometric similarity was associated with stronger inter-regional co-expression of genes enriched for neuronal terms. Comparing macaque MSNs with tract-tracing data confirmed that morphometric similarity was related to axonal connectivity. Finally, variation in the degree of human MSN nodes accounted for about 40% of between-subject variability in IQ. Morphometric similarity mapping provides a novel, robust, and biologically plausible approach to understanding how human cortical networks underpin individual differences in psychological functions.

Complex network topology is characteristic of many biological systems, including anatomical and functional brain networks (connectomes). Here, we first constructed a structural covariance network from MRI measures of cortical thickness on 296 healthy volunteers, aged 14–24 years. Next, we designed a new algorithm for matching sample locations from the Allen Brain Atlas to the nodes of the SCN. Subsequently we used this to define, transcriptomic brain networks by estimating gene co-expression between pairs of cortical regions. Finally, we explored the hypothesis that transcriptional networks and structural MRI connectomes are coupled. A transcriptional brain network (TBN) and a structural covariance network (SCN) were correlated across connection weights and showed qualitatively similar complex topological properties: assortativity, small-worldness, modularity, and a rich-club. In both networks, the weight of an edge was inversely related to the anatomical (Euclidean) distance between regions. There were differences between networks in degree and distance distributions: the transcriptional network had a less fat-tailed degree distribution and a less positively skewed distance distribution than the SCN. However, cortical areas connected to each other within modules of the SCN had significantly higher levels of whole genome co-expression than expected by chance. Nodes connected in the SCN had especially high levels of expression and co-expression of a human supragranular enriched (HSE) gene set that has been specifically located to supragranular layers of human cerebral cortex and is known to be important for large-scale, long-distance cortico-cortical connectivity. This coupling of brain transcriptome and connectome topologies was largely but not entirely accounted for by the common constraint of physical distance on both networks., Highlights • Transcriptomic Brain Network (TBN) is defined as inter-regional gene co-expression. • TBN has complex topological properties partially overlapped with structural networks. • Structural modules have higher gene co-expression than expected by chance. • Human Supragranular genes are highly expressed and coexpressed in structural hubs.

Complex network topology is characteristic of many biological systems, including anatomical and functional brain networks (connectomes). Here, we first constructed a structural covariance network (SCN) from MRI measures of cortical thickness on 296 healthy volunteers, aged 14-24 years. Next, we designed a new algorithm for matching sample locations from the Allen Brain Atlas to the nodes of the SCN. Subsequently we use this to define, transcriptomic brain networks (TBN) by estimating gene co-expression between pairs of cortical regions. Finally, we explore the hypothesis that TBN and the SCN are coupled.TBN and SCN were correlated across connection weights and showed qualitatively similar complex topological properties. There were differences between networks in degree and distance distributions. However, cortical areas connected to each other within modules of the SCN network had significantly higher levels of whole genome co-expression than expected by chance.Nodes connected in the SCN had significantly higher levels of expression and co-expression of a Human Supragranular Enriched (HSE) gene set that are known to be important for large-scale cortico-cortical connectivity. This coupling of brain transcriptome and connectome topologies was largely but not completely related to the common constraint of physical distance on both networks.

SummaryMacroscopic cortical networks are important for cognitive function, but it remains challenging to construct anatomically plausible individual structural connectomes from human neuroimaging. We introduce a new technique for cortical network mapping, based on inter-regional similarity of multiple morphometric parameters measured using multimodal MRI. In three cohorts (two human, one macaque), we find that the resulting morphometric similarity networks (MSNs) have a complex topological organisation comprising modules and high-degree hubs. Human MSN modules recapitulate known cortical cytoarchitectonic divisions, and greater inter-regional morphometric similarity was associated with stronger inter-regional co-expression of genes enriched for neuronal terms. Comparing macaque MSNs to tract-tracing data confirmed that morphometric similarity was related to axonal connectivity. Finally, variation in the degree of human MSN nodes accounted for about 40% of between-subject variability in IQ. Morphometric similarity mapping provides a novel, robust and biologically plausible approach to understanding how human cortical networks underpin individual differences in psychological functions.

Motivated by prior data on local cortical shrinkage and intracortical myelination, we predicted age-related changes in topological organization of cortical structural networks during adolescence. We estimated structural correlation from magnetic resonance imaging measures of cortical thickness at 308 regions in a sample of N = 297 healthy participants, aged 14–24 years. We used a novel sliding-window analysis to measure age-related changes in network attributes globally, locally and in the context of several community partitions of the network. We found that the strength of structural correlation generally decreased as a function of age. Association cortical regions demonstrated a sharp decrease in nodal degree (hubness) from 14 years, reaching a minimum at approximately 19 years, and then levelling off or even slightly increasing until 24 years. Greater and more prolonged age-related changes in degree of cortical regions within the brain network were associated with faster rates of adolescent cortical myelination and shrinkage. The brain regions that demonstrated the greatest age-related changes were concentrated within prefrontal modules. We conclude that human adolescence is associated with biologically plausible changes in structural imaging markers of brain network organization, consistent with the concept of tuning or consolidating anatomical connectivity between frontal cortex and the rest of the connectome.

At the macroscopic scale, the human brain can be described as a complex network of white matter tracts integrating grey matter assemblies — the human connectome. The structure of the connectome, which is often described using graph theoretic approaches, can be used to model macroscopic brain function at low computational cost. Here, we use the Kuramoto model of coupled oscillators with time-delays, calibrated with respect to empirical functional MRI data, to study the relation between the structure of the connectome and two aspects of functional brain dynamics — synchrony, a measure of general coherence, and metastability, a measure of dynamical flexibility. Specifically, we investigate the relationship between the local structure of the connectome, quantified using graph theory, and the synchrony and metastability of the model's dynamics. By removing individual nodes and all of their connections from the model, we study the effect of lesions on both global and local dynamics. Of the nine nodal graph-theoretical properties tested, two were able to predict effects of node lesion on the global dynamics. The removal of nodes with high eigenvector centrality leads to decreases in global synchrony and increases in global metastability, as does the removal of hub nodes joining topologically segregated network modules. At the level of local dynamics in the neighbourhood of the lesioned node, structural properties of the lesioned nodes hold more predictive power, as five nodal graph theoretical measures are related to changes in local dynamics following node lesions. We discuss these results in the context of empirical studies of stroke and functional brain dynamics.