Parameter estimation for connectome generative models: Accuracy, reliability, and a fast parameter fitting method.

Generative models of the human connectome enable in silico generation of brain networks based on probabilistic wiring rules. These wiring rules are governed by a small number of parameters that are typically fitted to individual connectomes and quantify the extent to which geometry and topology shape the generative process. A significant shortcoming of generative modeling in large cohort studies is that parameter estimation is computationally burdensome, and the accuracy and reliability of current estimation methods remain untested. Here, we propose a fast, reliable, and accurate parameter estimation method for connectome generative models that is scalable to large sample sizes. Our method achieves improved estimation accuracy and reliability and reduces computational cost by orders of magnitude, compared to established methods. We demonstrate an inherent tradeoff between accuracy, reliability, and computational expense in parameter estimation and provide recommendations for leveraging this tradeoff. To enable power analyses in future studies, we empirically approximate the minimum sample size required to detect between-group differences in generative model parameters. While we focus on the classic two-parameter generative model based on connection length and the topological matching index, our method can be generalized to other growth-based generative models. Our work provides a statistical and practical guide to parameter estimation for connectome generative models.

On the intersection between data quality and dynamical modelling of large-scale fMRI signals.

Large-scale dynamics of the brain are routinely modelled using systems of nonlinear dynamical equations that describe the evolution of population-level activity, with distinct neural populations often coupled according to an empirically measured structural connectivity matrix. This modelling approach has been used to generate insights into the neural underpinnings of spontaneous brain dynamics, as recorded with techniques such as resting state functional MRI (fMRI). In fMRI, researchers have many degrees of freedom in the way that they can process the data and recent evidence indicates that the choice of pre-processing steps can have a major effect on empirical estimates of functional connectivity. However, the potential influence of such variations on modelling results are seldom considered. Here we show, using three popular whole-brain dynamical models, that different choices during fMRI preprocessing can dramatically affect model fits and interpretations of findings. Critically, we show that the ability of these models to accurately capture patterns in fMRI dynamics is mostly driven by the degree to which they fit global signals rather than interesting sources of coordinated neural dynamics. We show that widespread deflections can arise from simple global synchronisation. We introduce a simple two-parameter model that captures these fluctuations and performs just as well as more complex, multi-parameter biophysical models. From our combined analyses of data and simulations, we describe benchmarks to evaluate model fit and validity. Although most models are not resilient to denoising, we show that relaxing the approximation of homogeneous neural populations by more explicitly modelling inter-regional effective connectivity can improve model accuracy at the expense of increased model complexity. Our results suggest that many complex biophysical models may be fitting relatively trivial properties of the data, and underscore a need for tighter integration between data quality assurance and model development.

The efficacy of different preprocessing steps in reducing motion-related confounds in diffusion MRI connectomics.

Head motion is a major confounding factor in neuroimaging studies. While numerous studies have investigated how motion impacts estimates of functional connectivity, the effects of motion on structural connectivity measured using diffusion MRI have not received the same level of attention, despite the fact that, like functional MRI, diffusion MRI relies on elaborate preprocessing pipelines that require multiple choices at each step. Here, we report a comprehensive analysis of how these choices influence motion-related contamination of structural connectivity estimates. Using a healthy adult sample (N = 294), we evaluated 240 different preprocessing pipelines, devised using plausible combinations of different choices related to explicit head motion correction, tractography propagation algorithms, track seeding methods, track termination constraints, quantitative metrics derived for each connectome edge, and parcellations. We found that an approach to motion correction that includes outlier replacement and within-slice volume correction led to a dramatic reduction in cross-subject correlations between head motion and structural connectivity strength, and that motion contamination is more severe when quantifying connectivity strength using mean tract fractional anisotropy rather than streamline count. We also show that the choice of preprocessing strategy can significantly influence subsequent inferences about network organization, with the location of network hubs varying considerably depending on the specific preprocessing steps applied. Our findings indicate that the impact of motion on structural connectivity can be successfully mitigated using recent motion-correction algorithms that include outlier replacement and within-slice motion correction.

Individual differences in haemoglobin concentration influence bold fMRI functional connectivity and its correlation with cognition.

Resting-state connectivity measures the temporal coherence of the spontaneous neural activity of spatially distinct regions, and is commonly measured using BOLD-fMRI. The BOLD response follows neuronal activity, when changes in the relative concentration of oxygenated and deoxygenated haemoglobin cause fluctuations in the MRI T2* signal. Since the BOLD signal detects changes in relative concentrations of oxy/deoxy-haemoglobin, individual differences in haemoglobin levels may influence the BOLD signal-to-noise ratio in a manner independent of the degree of neural activity. In this study, we examined whether group differences in haemoglobin may confound measures of functional connectivity. We investigated whether relationships between measures of functional connectivity and cognitive performance could be influenced by individual variability in haemoglobin. Finally, we mapped the neuroanatomical distribution of the influence of haemoglobin on functional connectivity to determine where group differences in functional connectivity are manifest. In a cohort of 518 healthy elderly subjects (259 men), each sex group was median-split into two groups with high and low haemoglobin concentration. Significant differences were obtained in functional connectivity between the high and low haemoglobin groups for both men and women (Cohen's d 0.17 and 0.03 for men and women respectively). The haemoglobin connectome in males showed a widespread systematic increase in functional connectivity correlation values, whilst the female connectome showed predominantly parietal and subcortical increases and temporo-parietal decreases. Despite the haemoglobin groups having no differences in cognitive measures, significant differences in the linear relationships between cognitive performance and functional connectivity were obtained for all 5 cognitive tests in males, and 4 out of 5 tests in females. Our findings confirm that individual variability in haemoglobin levels that give rise to group differences are an important confounding variable in BOLD-fMRI-based studies of functional connectivity. Controlling for haemoglobin variability as a potentially confounding variable is crucial to ensure the reproducibility of human brain connectome studies, especially in studies that compare groups of individuals, compare sexes, or examine connectivity-cognition relationships.

The anticipation and outcome phases of reward and loss processing: A neuroimaging meta-analysis of the monetary incentive delay task.

The processing of rewards and losses are crucial to everyday functioning. Considerable interest has been attached to investigating the anticipation and outcome phases of reward and loss processing, but results to date have been inconsistent. It is unclear if anticipation and outcome of a reward or loss recruit similar or distinct brain regions. In particular, while the striatum has widely been found to be active when anticipating a reward, whether it activates in response to the anticipation of losses as well remains ambiguous. Furthermore, concerning the orbitofrontal/ventromedial prefrontal regions, activation is often observed during reward receipt. However, it is unclear if this area is active during reward anticipation as well. We ran an Activation Likelihood Estimation meta-analysis of 50 fMRI studies, which used the Monetary Incentive Delay Task (MIDT), to identify which brain regions are implicated in the anticipation of rewards, anticipation of losses, and the receipt of reward. Anticipating rewards and losses recruits overlapping areas including the striatum, insula, amygdala and thalamus, suggesting that a generalised neural system initiates motivational processes independent of valence. The orbitofrontal/ventromedial prefrontal regions were recruited only during the reward outcome, likely representing the value of the reward received. Our findings help to clarify the neural substrates of the different phases of reward and loss processing, and advance neurobiological models of these processes.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization.

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Erythrocyte membrane fatty acid concentrations and myelin integrity in young people at ultra-high risk of psychosis.

Decreased white matter (WM) integrity and disturbance in fatty acid composition have been reported in individuals at ultra-high risk of psychosis (UHR). The current study is the first to investigate both WM integrity and erythrocyte membrane polyunsaturated fatty acid (PUFA) levels as potential risk biomarkers for persistent UHR status, and global functioning in UHR individuals. Forty UHR individuals were analysed at baseline for erythrocyte membrane PUFA concentrates. Tract-based spatial statistics (TBSS) was used to analyse fractional anisotropy (FA) and diffusivity measures. Measures of global functioning and psychiatric symptoms were evaluated at baseline and at 12-months. Fatty acids and WM indices did not predict functional outcomes at baseline or 12-months. Significant differences were found in FA between UHR remitters and non-remitters (individuals who no longer met UHR criteria versus those who continued to meet criteria at 12-months). Docosahexaenoic acid (DHA) was found to be a significant predictor of UHR status at 12-months, as was the interaction between the sum of ώ-3 and whole brain FA, and the interaction between the right anterior limb of the internal capsule and the sum of ώ-3. The results confirm that certain fatty acids have a unique relationship with WM integrity in UHR individuals.

A phylogenetically-conserved axis of thalamocortical connectivity in the human brain.

The thalamus enables key sensory, motor, emotive, and cognitive processes via connections to the cortex. These projection patterns are traditionally considered to originate from discrete thalamic nuclei, however recent work showing gradients of molecular and connectivity features in the thalamus suggests the organisation of thalamocortical connections occurs along a continuous dimension. By performing a joint decomposition of densely sampled gene expression and non-invasive diffusion tractography in the adult human thalamus, we define a principal axis of genetic and connectomic variation along a medial-lateral thalamic gradient. Projections along this axis correspond to an anterior-posterior cortical pattern and are aligned with electrophysiological properties of the cortex. The medial-lateral axis demonstrates phylogenetic conservation, reflects transitions in neuronal subtypes, and shows associations with neurodevelopment and common brain disorders. This study provides evidence for a supra-nuclear axis of thalamocortical organisation characterised by a graded transition in molecular properties and anatomical connectivity.

Can hubs of the human connectome be identified consistently with diffusion MRI?

Recent years have seen a surge in the use of diffusion MRI to map connectomes in humans, paralleled by a similar increase in processing and analysis choices. Yet these different steps and their effects are rarely compared systematically. Here, in a healthy young adult population (n = 294), we characterized the impact of a range of analysis pipelines on one widely studied property of the human connectome: its degree distribution. We evaluated the effects of 40 pipelines (comparing common choices of parcellation, streamline seeding, tractography algorithm, and streamline propagation constraint) and 44 group-representative connectome reconstruction schemes on highly connected hub regions. We found that hub location is highly variable between pipelines. The choice of parcellation has a major influence on hub architecture, and hub connectivity is highly correlated with regional surface area in most of the assessed pipelines (ρ > 0.70 in 69% of the pipelines), particularly when using weighted networks. Overall, our results demonstrate the need for prudent decision-making when processing diffusion MRI data, and for carefully considering how different processing choices can influence connectome organization.

Network-Based Spreading of Gray Matter Changes Across Different Stages of Psychosis.

Importance: Psychotic illness is associated with anatomically distributed gray matter reductions that can worsen with illness progression, but the mechanisms underlying the specific spatial patterning of these changes is unknown. Objective: To test the hypothesis that brain network architecture constrains cross-sectional and longitudinal gray matter alterations across different stages of psychotic illness and to identify whether certain brain regions act as putative epicenters from which volume loss spreads. Design, Settings, and Participants: This case-control study included 534 individuals from 4 cohorts, spanning early and late stages of psychotic illness. Early-stage cohorts included patients with antipsychotic-naive first-episode psychosis (n = 59) and a group of patients receiving medications within 3 years of psychosis onset (n = 121). Late-stage cohorts comprised 2 independent samples of people with established schizophrenia (n = 136). Each patient group had a corresponding matched control group (n = 218). A sample of healthy adults (n = 356) was used to derive representative structural and functional brain networks for modeling of network-based spreading processes. Longitudinal illness-related and antipsychotic-related gray matter changes over 3 and 12 months were examined using a triple-blind randomized placebo-control magnetic resonance imaging study of the antipsychotic-naive patients. All data were collected between April 29, 2008, and January 15, 2020, and analyses were performed between March 1, 2021, and January 14, 2023. Main Outcomes and Measures: Coordinated deformation models were used to estimate the extent of gray matter volume (GMV) change in each of 332 parcellated areas by the volume changes observed in areas to which they were structurally or functionally coupled. To identify putative epicenters of volume loss, a network diffusion model was used to simulate the spread of pathology from different seed regions. Correlations between estimated and empirical spatial patterns of GMV alterations were used to quantify model performance. Results: Of 534 included individuals, 354 (66.3%) were men, and the mean (SD) age was 28.4 (7.4) years. In both early and late stages of illness, spatial patterns of cross-sectional volume differences between patients and controls were more accurately estimated by coordinated deformation models constrained by structural, rather than functional, network architecture (r range, >0.46 to <0.57; P < .01). The same model also robustly estimated longitudinal volume changes related to illness (r ≥ 0.52; P < .001) and antipsychotic exposure (r ≥ 0.50; P < .004). Network diffusion modeling consistently identified, across all 4 data sets, the anterior hippocampus as a putative epicenter of pathological spread in psychosis. Epicenters of longitudinal GMV loss were apparent in posterior cortex early in the illness and shifted to the prefrontal cortex with illness progression. Conclusion and Relevance: These findings highlight a central role for white matter fibers as conduits for the spread of pathology across different stages of psychotic illness, mirroring findings reported in neurodegenerative conditions. The structural connectome thus represents a fundamental constraint on brain changes in psychosis, regardless of whether these changes are caused by illness or medication. Moreover, the anterior hippocampus represents a putative epicenter of early brain pathology from which dysfunction may spread to affect connected areas.

Early and late development of hub connectivity in the human brain.

Human brain networks undergo pronounced changes during development. The emergence of highly connected hub regions that can support integrated brain function is central to this maturational process, with these areas undergoing a particularly protracted period of development that extends into adulthood. The location of cortical network hubs emerges early but connections to and from hubs continue to strengthen throughout childhood and adolescence. Patterns of functional coupling in cortical association hubs are immature and incomplete at birth, but gradually strengthen during development. Early establishment of hub connectivity may provide a stable substrate that is refined by changes in tissue organization and microstructure, resulting in the emergence of complex functional dynamics by adulthood.

Modeling spatial, developmental, physiological, and topological constraints on human brain connectivity.

The complex connectivity of nervous systems is thought to have been shaped by competitive selection pressures to minimize wiring costs and support adaptive function. Accordingly, recent modeling work indicates that stochastic processes, shaped by putative trade-offs between the cost and value of each connection, can successfully reproduce many topological properties of macroscale human connectomes measured with diffusion magnetic resonance imaging. Here, we derive a new formalism that more accurately captures the competing pressures of wiring cost minimization and topological complexity. We further show that model performance can be improved by accounting for developmental changes in brain geometry and associated wiring costs, and by using interregional transcriptional or microstructural similarity rather than topological wiring rules. However, all models struggled to capture topographical (i.e., spatial) network properties. Our findings highlight an important role for genetics in shaping macroscale brain connectivity and indicate that stochastic models offer an incomplete account of connectome organization.

Investigating longitudinal changes to frontal cortico-striatal tracts in Huntington's disease: the IMAGE-HD study.

The striatum is the principal site of disease pathology in Huntington's disease and contains neural connections to numerous cortical brain regions. Studies examining abnormalities to neural connections find that white matter integrity is compromised in HD; however, further regional, and longitudinal investigation is required. This paper is the first longitudinal investigation into region-based white-matter integrity changes in Huntington's Disease. The aim of this study was to better understand how disease progression impacts white matter tracts connecting the striatum to the prefrontal and motor cortical regions in HD. We used existing neuroimaging data from IMAGE-HD, comprised of 25 pre-symptomatic, 27 symptomatic, and 25 healthy controls at three separate time points (baseline, 18-months, 30-months). Fractional anisotropy, axial diffusivity and radial diffusivity were derived as measures of white matter microstructure. The anatomical regions of interest were identified using the Desikan-Killiany brain atlas. A Group by Time repeated measures ANCOVA was conducted for each tract of interest and for each measure. We found significantly lower fractional anisotropy and significantly higher radial diffusivity in the symptomatic group, compared to both the pre-symptomatic group and controls (the latter two groups did not differ from each other), in the rostral middle frontal and superior frontal tracts; as well as significantly higher axial diffusivity in the rostral middle tracts only. We did not find a Group by Time interaction for any of the white matter integrity measures. These findings demonstrate that whilst the microstructure of white matter tracts, extending from the striatum to these regions of interest, are compromised during the symptomatic stages of Huntington's disease, 36-month follow-up did not show progressive changes in these measures. Additionally, no correlations were found between clinical measures and tractography changes, indicating further investigations into the relationship between tractography changes and clinical symptoms in Huntington's disease are required.

Associations of cardiorespiratory fitness and exercise with brain white matter in healthy adults: A systematic review and meta-analysis.

Magnetic resonance imaging (MRI) studies have revealed positive associations between brain structure and physical activity, cardiorespiratory fitness, and exercise (referred to here as PACE). While a considerable body of research has investigated the effects of PACE on grey matter, much less is known about effects on white matter (WM). Hence, we conducted a systematic review of peer-reviewed literature published prior to 5th July 2021 using online databases (PubMed and Scopus) and PRISMA guidelines to synthesise what is currently known about the relationship between PACE and WM in healthy adults. A total of 60 studies met inclusion criteria and were included in the review. Heterogeneity across studies was calculated using Qochran's q test, and publication bias was assessed for each meta-analysis using Begg and Mazumdar rank correlation test. A meta-regression was also conducted to explore factors contributing to any observed heterogeneity. Overall, we observed evidence of positive associations between PACE and global WM volume (effect size (Hedges's g) = 0.137, p < 0.001), global WM anomalies (effect size = 0.182, p < 0.001), and local microstructure integrity (i.e., corpus callosum: effect size = 0.345, p < 0.001, and anterior limb of internal capsule: effect size = 0.198, p < 0.001). These findings suggest that higher levels of PACE are associated with improved global WM volume and local integrity. We appraise the quality of evidence, and discuss the implications of these findings for the preservation of WM across the lifespan. We conclude by providing recommendations for future research in order to advance our understanding of the specific PACE parameters and neurobiological mechanisms underlying these effects.

Dynamical consequences of regional heterogeneity in the brain's transcriptional landscape.

Brain regions vary in their molecular and cellular composition, but how this heterogeneity shapes neuronal dynamics is unclear. Here, we investigate the dynamical consequences of regional heterogeneity using a biophysical model of whole-brain functional magnetic resonance imaging (MRI) dynamics in humans. We show that models in which transcriptional variations in excitatory and inhibitory receptor (E:I) gene expression constrain regional heterogeneity more accurately reproduce the spatiotemporal structure of empirical functional connectivity estimates than do models constrained by global gene expression profiles or MRI-derived estimates of myeloarchitecture. We further show that regional transcriptional heterogeneity is essential for yielding both ignition-like dynamics, which are thought to support conscious processing, and a wide variance of regional-activity time scales, which supports a broad dynamical range. We thus identify a key role for E:I heterogeneity in generating complex neuronal dynamics and demonstrate the viability of using transcriptomic data to constrain models of large-scale brain function.

Genetic influences on hub connectivity of the human connectome.

Brain network hubs are both highly connected and highly inter-connected, forming a critical communication backbone for coherent neural dynamics. The mechanisms driving this organization are poorly understood. Using diffusion-weighted magnetic resonance imaging in twins, we identify a major role for genes, showing that they preferentially influence connectivity strength between network hubs of the human connectome. Using transcriptomic atlas data, we show that connected hubs demonstrate tight coupling of transcriptional activity related to metabolic and cytoarchitectonic similarity. Finally, comparing over thirteen generative models of network growth, we show that purely stochastic processes cannot explain the precise wiring patterns of hubs, and that model performance can be improved by incorporating genetic constraints. Our findings indicate that genes play a strong and preferential role in shaping the functionally valuable, metabolically costly connections between connectome hubs.

Timescales of spontaneous fMRI fluctuations relate to structural connectivity in the brain.

Intrinsic timescales of activity fluctuations vary hierarchically across the brain. This variation reflects a broad gradient of functional specialization in information storage and processing, with integrative association areas displaying slower timescales that are thought to reflect longer temporal processing windows. The organization of timescales is associated with cognitive function, distinctive between individuals, and disrupted in disease, but we do not yet understand how the temporal properties of activity dynamics are shaped by the brain's underlying structural connectivity network. Using resting-state fMRI and diffusion MRI data from 100 healthy individuals from the Human Connectome Project, here we show that the timescale of resting-state fMRI dynamics increases with structural connectivity strength, matching recent results in the mouse brain. Our results hold at the level of individuals, are robust to parcellation schemes, and are conserved across a range of different timescale- related statistics. We establish a comprehensive BOLD dynamical signature of structural connectivity strength by comparing over 6,000 time series features, highlighting a range of new temporal features for characterizing BOLD dynamics, including measures of stationarity and symbolic motif frequencies. Our findings indicate a conserved property of mouse and human brain organization in which a brain region's spontaneous activity fluctuations are closely related to their surrounding structural scaffold.

The development of brain network hubs.

Some brain regions have a central role in supporting integrated brain function, marking them as network hubs. Given the functional importance of hubs, it is natural to ask how they emerge during development and to consider how they shape the function of the maturing brain. Here, we review evidence examining how brain network hubs, both in structural and functional connectivity networks, develop over the prenatal, neonate, childhood, and adolescent periods. The available evidence suggests that structural hubs of the brain arise in the prenatal period and show a consistent spatial topography through development, but undergo a protracted period of consolidation that extends into late adolescence. In contrast, the hubs of brain functional networks show a more variable topography, being predominantly located in primary cortical areas in early development, before moving to association areas by late childhood. These findings suggest that while the basic anatomical infrastructure of hubs may be established early, the functional viability and integrative capacity of these areas undergoes extensive postnatal maturation. Not all findings are consistent with this view however. We consider methodological factors that might drive these inconsistencies, and which should be addressed to promote a more rigorous investigation of brain network development.

Consistency and differences between centrality measures across distinct classes of networks.

The roles of different nodes within a network are often understood through centrality analysis, which aims to quantify the capacity of a node to influence, or be influenced by, other nodes via its connection topology. Many different centrality measures have been proposed, but the degree to which they offer unique information, and whether it is advantageous to use multiple centrality measures to define node roles, is unclear. Here we calculate correlations between 17 different centrality measures across 212 diverse real-world networks, examine how these correlations relate to variations in network density and global topology, and investigate whether nodes can be clustered into distinct classes according to their centrality profiles. We find that centrality measures are generally positively correlated to each other, the strength of these correlations varies across networks, and network modularity plays a key role in driving these cross-network variations. Data-driven clustering of nodes based on centrality profiles can distinguish different roles, including topological cores of highly central nodes and peripheries of less central nodes. Our findings illustrate how network topology shapes the pattern of correlations between centrality measures and demonstrate how a comparative approach to network centrality can inform the interpretation of nodal roles in complex networks.