Optimization of graph construction can significantly increase the power of structural brain network studies

Modeling the impact of MRI acquisition bias on structural connectomes: Harmonizing structural connectomes

One way to increase the statistical power and generalizability of neuroimaging studies is to collect data at multiple sites or merge multiple cohorts. However, this usually comes with site-related biases due to the heterogeneity of scanners and acquisition parameters, negatively impacting sensitivity. Brain structural connectomes are not an exception: Being derived from T1-weighted and diffusion-weighted magnetic resonance images, structural connectivity is impacted by differences in imaging protocol. Beyond minimizing acquisition parameter differences, removing bias with postprocessing is essential. In this work we create, from the exhaustive Human Connectome Project Young Adult dataset, a resampled dataset of different b-values and spatial resolutions, modeling a cohort scanned across multiple sites. After demonstrating the statistical impact of acquisition parameters on connectivity, we propose a linear regression with explicit modeling of b-value and spatial resolution, and validate its performance on separate datasets. We show that b-value and spatial resolution affect connectivity in different ways and that acquisition bias can be reduced using a linear regression informed by the acquisition parameters while retaining interindividual differences and hence boosting fingerprinting performance. We also demonstrate the generative potential of our model, and its generalization capability in an independent dataset reflective of typical acquisition practices in clinical settings.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

Similarity between structural and proxy estimates of brain connectivity

Functional magnetic resonance and diffusion weighted imaging have so far made a major contribution to delineation of the brain connectome at the macroscale. While functional connectivity (FC) was shown to be related to structural connectivity (SC) to a certain degree, their spatial overlap is unknown. Even less clear are relations of SC with estimates of connectivity from inter-subject covariance of regional F18-fluorodeoxyglucose uptake (FDGcov) and grey matter volume (GMVcov). Here, we asked to what extent SC underlies three proxy estimates of brain connectivity: FC, FDGcov and GMVcov. Simultaneous PET/MR acquisitions were performed in 56 healthy middle-aged individuals. Similarity between four networks was assessed using Spearman correlation and convergence ratio (CR), a measure of spatial overlap. Spearman correlation coefficient was 0.27 for SC-FC, 0.40 for SC-FDGcov, and 0.15 for SC-GMVcov. Mean CRs were 51% for SC-FC, 48% for SC-FDGcov, and 37% for SC-GMVcov. These results proved to be reproducible and robust against image processing steps. In sum, we found a relevant similarity of SC with FC and FDGcov, while GMVcov consistently showed the weakest similarity. These findings indicate that white matter tracts underlie FDGcov to a similar degree as FC, supporting FDGcov as estimate of functional brain connectivity.

Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Normative values of the topological metrics of the structural connectome: A multi-site reproducibility study across the Italian Neuroscience network

PURPOSE

The use of topological metrics to derive quantitative descriptors from structural connectomes is receiving increasing attention but deserves specific studies to investigate their reproducibility and variability in the clinical context. This work exploits the harmonization of diffusion-weighted acquisition for neuroimaging data performed by the Italian Neuroscience and Neurorehabilitation Network initiative to obtain normative values of topological metrics and to investigate their reproducibility and variability across centers.

METHODS

Different topological metrics, at global and local level, were calculated on multishell diffusion-weighted data acquired at high-field (e.g. 3 T) Magnetic Resonance Imaging scanners in 13 different centers, following the harmonization of the acquisition protocol, on young and healthy adults. A "traveling brains" dataset acquired on a subgroup of subjects at 3 different centers was also analyzed as reference data. All data were processed following a common processing pipeline that includes data pre-processing, tractography, generation of structural connectomes and calculation of graph-based metrics. The results were evaluated both with statistical analysis of variability and consistency among sites with the traveling brains range. In addition, inter-site reproducibility was assessed in terms of intra-class correlation variability.

RESULTS

The results show an inter-center and inter-subject variability of <10%, except for "clustering coefficient" (variability of 30%). Statistical analysis identifies significant differences among sites, as expected given the wide range of scanners' hardware.

CONCLUSIONS

The results show low variability of connectivity topological metrics across sites running a harmonised protocol.

Strong intercorrelations among global graph-theoretic indices of structural connectivity in the human brain

Graph-theoretic metrics derived from neuroimaging data have been heralded as powerful tools for uncovering neural mechanisms of psychological traits, psychiatric disorders, and neurodegenerative diseases. In N = 8,185 human structural connectomes from UK Biobank, we examined the extent to which 11 commonly-used global graph-theoretic metrics index distinct versus overlapping information with respect to interindividual differences in brain organization. Using unthresholded, FA-weighted networks we found that all metrics other than Participation Coefficient were highly intercorrelated, both with each other (mean |r| = 0.788) and with a topologically-naïve summary index of brain structure (mean edge weight; mean |r| = 0.873). In a series of sensitivity analyses, we found that overlap between metrics is influenced by the sparseness of the network and the magnitude of variation in edge weights. Simulation analyses representing a range of population network structures indicated that individual differences in global graph metrics may be intrinsically difficult to separate from mean edge weight. In particular, Closeness, Characteristic Path Length, Global Efficiency, Clustering Coefficient, and Small Worldness were nearly perfectly collinear with one another (mean |r| = 0.939) and with mean edge weight (mean |r| = 0.952) across all observed and simulated conditions. Global graph-theoretic measures are valuable for their ability to distill a high-dimensional system of neural connections into summary indices of brain organization, but they may be of more limited utility when the goal is to index separable components of interindividual variation in specific properties of the human structural connectome.

Genetic risk for schizophrenia is associated with increased proportion of indirect connections in brain networks revealed by a semi-metric analysis: evidence from population sample stratified for polygenic risk

Research studies based on tractography have revealed a prominent reduction of asymmetry in some key white-matter tracts in schizophrenia (SCZ). However, we know little about the influence of common genetic risk factors for SCZ on the efficiency of routing on structural brain networks (SBNs). Here, we use a novel recall-by-genotype approach, where we sample young adults from a population-based cohort (ALSPAC:N genotyped = 8,365) based on their burden of common SCZ risk alleles as defined by polygenic risk score (PRS). We compared 181 individuals at extremes of low (N = 91) or high (N = 90) SCZ-PRS under a robust diffusion MRI-based graph theoretical SBN framework. We applied a semi-metric analysis revealing higher SMR values for the high SCZ-PRS group compared with the low SCZ-PRS group in the left hemisphere. Furthermore, a hemispheric asymmetry index showed a higher leftward preponderance of indirect connections for the high SCZ-PRS group compared with the low SCZ-PRS group (PFDR < 0.05). These findings might indicate less efficient structural connectivity in the higher genetic risk group. This is the first study in a population-based sample that reveals differences in the efficiency of SBNs associated with common genetic risk variants for SCZ.

Increased structural connectivity in high schizotypy

The link between brain structural connectivity and schizotypy was explored in two healthy participant cohorts, collected at two different neuroimaging centres, comprising 140 and 115 participants, respectively. The participants completed the Schizotypal Personality Questionnaire (SPQ), through which their schizotypy scores were calculated. Diffusion-MRI data were used to perform tractography and to generate the structural brain networks of the participants. The edges of the networks were weighted with the inverse radial diffusivity. Graph theoretical metrics of the default mode, sensorimotor, visual, and auditory subnetworks were derived and their correlation coefficients with the schizotypy scores were calculated. To the best of our knowledge, this is the first time that graph theoretical measures of structural brain networks are investigated in relation to schizotypy. A positive correlation was found between the schizotypy score and the mean node degree and mean clustering coefficient of the sensorimotor and the default mode subnetworks. The nodes driving these correlations were the right postcentral gyrus, the left paracentral lobule, the right superior frontal gyrus, the left parahippocampal gyrus, and the bilateral precuneus, that is, nodes that exhibit compromised functional connectivity in schizophrenia. Implications for schizophrenia and schizotypy are discussed.

Assessing the Repeatability of Multi-Frequency Multi-Layer Brain Network Topologies Across Alternative Researcher's Choice Paths

There is a growing interest in the neuroscience community on the advantages of multilayer functional brain networks. Researchers usually treated different frequencies separately at distinct functional brain networks. However, there is strong evidence that these networks share complementary information while their interdependencies could reveal novel findings. For this purpose, neuroscientists adopt multilayer networks, which can be described mathematically as an extension of trivial single-layer networks. Multilayer networks have become popular in neuroscience due to their advantage to integrate different sources of information. Here, Ι will focus on the multi-frequency multilayer functional connectivity analysis on resting-state fMRI (rs-fMRI) recordings. However, constructing a multilayer network depends on selecting multiple pre-processing steps that can affect the final network topology. Here, I analyzed the rs-fMRI dataset from a single human performing scanning over a period of 18 months (84 scans in total), and the rs-fMRI dataset containing 25 subjects with 3 repeat scans. I focused on assessing the reproducibility of multi-frequency multilayer topologies exploring the effect of two filtering methods for extracting frequencies from BOLD activity, three connectivity estimators, with or without a topological filtering scheme, and two spatial scales. Finally, I untangled specific combinations of researchers' choices that yield consistently brain networks with repeatable topologies, giving me the chance to recommend best practices over consistent topologies.

The impact of genetic risk for Alzheimer's disease on the structural brain networks of young adults

Introduction

We investigated the structural brain networks of 562 young adults in relation to polygenic risk for Alzheimer's disease, using magnetic resonance imaging (MRI) and genotype data from the Avon Longitudinal Study of Parents and Children.

Methods

Diffusion MRI data were used to perform whole-brain tractography and generate structural brain networks for the whole-brain connectome, and for the default mode, limbic and visual subnetworks. The mean clustering coefficient, mean betweenness centrality, characteristic path length, global efficiency and mean nodal strength were calculated for these networks, for each participant. The connectivity of the rich-club, feeder and local connections was also calculated. Polygenic risk scores (PRS), estimating each participant's genetic risk, were calculated at genome-wide level and for nine specific disease pathways. Correlations were calculated between the PRS and (a) the graph theoretical metrics of the structural networks and (b) the rich-club, feeder and local connectivity of the whole-brain networks.

Results

In the visual subnetwork, the mean nodal strength was negatively correlated with the genome-wide PRS (r = -0.19, p = 1.4 × 10-3), the mean betweenness centrality was positively correlated with the plasma lipoprotein particle assembly PRS (r = 0.16, p = 5.5 × 10-3), and the mean clustering coefficient was negatively correlated with the tau-protein binding PRS (r = -0.16, p = 0.016). In the default mode network, the mean nodal strength was negatively correlated with the genome-wide PRS (r = -0.14, p = 0.044). The rich-club and feeder connectivities were negatively correlated with the genome-wide PRS (r = -0.16, p = 0.035; r = -0.15, p = 0.036).

Discussion

We identified small reductions in brain connectivity in young adults at risk of developing Alzheimer's disease in later life.

Improved sensitivity and precision in multicentre diffusion MRI network analysis using thresholding and harmonization

PURPOSE

To investigate if network thresholding and raw data harmonization improve consistency of diffusion MRI (dMRI)-based brain networks while also increasing precision and sensitivity to detect disease effects in multicentre datasets.

METHODS

Brain networks were reconstructed from dMRI of five samples with cerebral small vessel disease (SVD; 629 patients, 166 controls), as a clinically relevant exemplar condition for studies on network integrity. We evaluated consistency of network architecture in age-matched controls, by calculating cross-site differences in connection probability and fractional anisotropy (FA). Subsequently we evaluated precision and sensitivity to disease effects by identifying connections with low FA in sporadic SVD patients relative to controls, using more severely affected patients with a pure form of genetically defined SVD as reference.

RESULTS

In controls, thresholding and harmonization improved consistency of network architecture, minimizing cross-site differences in connection probability and FA. In patients relative to controls, thresholding improved precision to detect disrupted connections by removing false positive connections (precision, before: 0.09-0.19; after: 0.38-0.70). Before harmonization, sensitivity was low within individual sites, with few connections surviving multiple testing correction (k = 0-25 connections). Harmonization and pooling improved sensitivity (k = 38), while also achieving higher precision when combined with thresholding (0.97).

CONCLUSION

We demonstrated that network consistency, precision and sensitivity to detect disease effects in SVD are improved by thresholding and harmonization. We recommend introducing these techniques to leverage large existing multicentre datasets to better understand the impact of disease on brain networks.

Accuracy and reliability of diffusion imaging models

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927-1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3's Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

Impact of thresholding on the consistency and sensitivity of diffusion MRI-based brain networks in patients with cerebral small vessel disease

INTRODUCTION

Thresholding of low-weight connections of diffusion MRI-based brain networks has been proposed to remove false-positive connections. It has been previously established that this yields more reproducible scan-rescan network architecture in healthy subjects. In patients with brain disease, network measures are applied to assess inter-individual variation and changes over time. Our aim was to investigate whether thresholding also achieves improved consistency in network architecture in patients, while maintaining sensitivity to disease effects for these applications.

METHODS

We applied fixed-density and absolute thresholding on brain networks in patients with cerebral small vessel disease (SVD, n = 86; ≈24 months follow-up), as a clinically relevant exemplar condition. In parallel, we applied the same methods in healthy young subjects (n = 44; scan-rescan interval ≈4 months) as a frame of reference. Consistency of network architecture was assessed with dice similarity of edges and intraclass correlation coefficient (ICC) of edge-weights and hub-scores. Sensitivity to disease effects in patients was assessed by evaluating interindividual variation, changes over time, and differences between those with high and low white matter hyperintensity burden, using correlation analyses and mixed ANOVA.

RESULTS

Compared to unthresholded networks, both thresholding methods generated more consistent architecture over time in patients (unthresholded: dice = .70; ICC: .70-.78; thresholded: dice = .77; ICC: .73-.83). However, absolute thresholding created fragmented nodes. Similar observations were made in the reference group. Regarding sensitivity to disease effects in patients, fixed-density thresholds that were optimal in terms of consistency (densities: .10-.30) preserved interindividual variation in global efficiency and node strength as well as the sensitivity to detect effects of time and group. Absolute thresholding produced larger fluctuations of interindividual variation.

CONCLUSIONS

Our results indicate that thresholding of low-weight connections, particularly when using fixed-density thresholding, results in more consistent network architecture in patients with longer rescan intervals, while preserving sensitivity to disease effects.

Regional Brain Fusion: Graph Convolutional Network for Alzheimer's Disease Prediction and Analysis

Alzheimer's disease (AD) has raised extensive concern in healthcare and academia as one of the most prevalent health threats to the elderly. Due to the irreversible nature of AD, early and accurate diagnoses are significant for effective prevention and treatment. However, diverse clinical symptoms and limited neuroimaging accuracy make diagnoses challenging. In this article, we built a brain network for each subject, which assembles several commonly used neuroimaging data simply and reasonably, including structural magnetic resonance imaging (MRI), diffusion-weighted imaging (DWI), and amyloid positron emission tomography (PET). Based on some existing research results, we applied statistical methods to analyze (i) the distinct affinity of AD burden on each brain region, (ii) the topological lateralization between left and right hemispheric sub-networks, and (iii) the asymmetry of the AD attacks on the left and right hemispheres. In the light of advances in graph convolutional networks for graph classifications and summarized characteristics of brain networks and AD pathologies, we proposed a regional brain fusion-graph convolutional network (RBF-GCN), which is constructed with an RBF framework mainly, including three sub-modules, namely, hemispheric network generation module, multichannel GCN module, and feature fusion module. In the multichannel GCN module, the improved GCN by our proposed adaptive native node attribute (ANNA) unit embeds within each channel independently. We not only fully verified the effectiveness of the RBF framework and ANNA unit but also achieved competitive results in multiple sets of AD stages' classification tasks using hundreds of experiments over the ADNI clinical dataset.

Dementia Risk Factors Modify Hubs but Leave Other Connectivity Measures Unchanged in Asymptomatic Individuals: A Graph Theoretical Analysis

Background: Alzheimer's disease (AD) is the most common form of dementia with genetic and environmental risk contributing to its development. Graph theoretical analyses of brain networks constructed from structural and functional magnetic resonance imaging (MRI) measurements have identified connectivity changes in AD and individuals with mild cognitive impairment. However, brain connectivity in asymptomatic individuals at risk of AD remains poorly understood. Methods: We analyzed diffusion-weighted MRI data from 161 asymptomatic individuals (38-71 years) from the Cardiff Ageing and Risk of Dementia Study (CARDS). We calculated white matter tracts and constructed whole-brain, default mode network (DMN) and visual structural brain networks that incorporate multiple structural metrics as edge weights. We then calculated the relationship of three AD risk factors, namely Apolipoprotein-E ɛ4 (APOE4) genotype, family history of dementia (FH), and central obesity (Waist-Hip-Ratio [WHR]), on graph theoretical measures and hubs. Results: We observed no risk-related differences in clustering coefficients, characteristic path lengths, eccentricity, diameter, and radius across the whole-brain, DMN or visual system. However, a hub in the right paracentral lobule was present in all the high-risk groups (FH, APOE4, obese), but absent in low-risk groups (no FH, APOE4-ve, healthy WHR). Discussion: We identified no risk-related effects on graph theoretical metrics in the structural brain networks of cognitively healthy individuals. However, high risk was associated with a hub in the right paracentral lobule, a medial fronto-parietal cortical area with motor and sensory functions. This finding is consistent with accumulating evidence for right parietal cortex contributions in AD. If this phenotype is shown to predict symptom development in longitudinal studies, it could be used as an early biomarker of AD. Impact statement Alzheimer's disease (AD) is a common form of dementia that to date has no cure. Identifying early biomarkers will aid the discovery and development of treatments that may slow AD progression in the future. In this article, we report that asymptomatic individuals at heightened risk of dementia due to their family history, Apolipoprotein-E ɛ4 genotype, and central adiposity have a hub in the right paracentral lobule, which is absent in low-risk groups. If this phenotype were to predict the development of symptoms in a longitudinal study of the same cohort, it could provide an early biomarker of disease progression.

Structural Brain Network Reproducibility: Influence of Different Diffusion Acquisition and Tractography Reconstruction Schemes on Graph Metrics

Background: Graph metrics of structural brain networks demonstrate to be a powerful tool for investigating brain topology at a large scale. However, the variability of the results related to applying different magnetic resonance acquisition schemes and tractography reconstruction techniques is not fully characterized. Materials and Methods: The present work aims to evaluate the influence of different combinations of diffusion acquisition schemes (single and multishell), diffusion models (tensor and spherical deconvolution), and tractography reconstruction approaches (deterministic and probabilistic) on the reproducibility of graph metrics derived from structural connectome on test/retest (TRT) data released by the Human Connectome Project. From each implemented experimental setup, both global and local graph metrics were evaluated and their reproducibility was estimated by the intraclass correlation coefficient (ICC). Moreover, the percentage relative standard deviation (pRSD) from the ICC values of local graph metrics was calculated to quantify how much the reproducibility varied across nodes within each experimental setup. Results: The presented results show that different combinations of diffusion acquisition schemes, diffusion models, and tractography algorithms can strongly affect the reproducibility of global and local graph metrics. The combination of constrained spherical deconvolution (CSD) and deterministic tractography gave generally high reproducibility (ICCs >0.75) and lowest pRSD for the considered graph metrics, meanwhile probabilistic CSD with a high b-value returned the highest reproducibility. Notably, the introduction of streamline selection filters on CSD can substantially affect the reproducibility. Discussion: This work demonstrates that the TRT reproducibility of graph metrics is generally high but can vary substantially with different combinations of acquisition and reconstruction schemes. Impact statement This work demonstrates the influence of different diffusion acquisition schemes, diffusion models, and tractography reconstruction approaches on the reproducibility of graph metrics derived from structural connectome. The presented findings impact on the choice of both acquisition protocol and processing pipeline for topological analyses to produce reproducible measurements for brain network studies.

Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition

Characterization of brain networks by diffusion MRI (dMRI) has rapidly evolved, and there are ongoing movements toward data sharing and multi-center studies. To extract meaningful information from multi-center data, methods to correct for the bias caused by scanner differences, that is, harmonization, are urgently needed. In this work, we report the cross-scanner differences in structural network analyses using data from nine traveling subjects (four males and five females, 21-49 years-old) who underwent scanning using four 3T scanners (public database available from the Brain/MINDS Beyond Human Brain MRI project (http://mriportal.umin.jp/)). The reliability and reproducibility were compared to those of data from another set of four subjects (all males, 29-42 years-old) who underwent scan-rescan (interval, 105-147 days) with the same scanner as well as scan-rescan data from the Human Connectome Project database. The results demonstrated that the reliability of the edge weights and graph theory metrics was lower for data including different scanners, compared to the scan-rescan with the same scanner. Besides, systematic differences between scanners were observed, indicating the risk of bias in comparing networks obtained from different scanners directly. We further demonstrate that it is feasible to reduce inter-scanner variabilities while preserving the inter-subject differences among healthy individuals by modeling the scanner effects at the level of network matrices, when traveling-subject data are available for calibration between scanners. The present data and results are expected to serve as a basis for developing and evaluating novel harmonization methods.

Dyconnmap: Dynamic connectome mapping-A neuroimaging python module

Despite recent progress in the analysis of neuroimaging data sets, our comprehension of the main mechanisms and principles which govern human brain cognition and function remains incomplete. Network neuroscience makes substantial efforts to manipulate these challenges and provide real answers. For the last decade, researchers have been modelling brain structure and function via a graph or network that comprises brain regions that are either anatomically connected via tracts or functionally via a more extensive repertoire of functional associations. Network neuroscience is a relatively new multidisciplinary scientific avenue of the study of complex systems by pursuing novel ways to analyze, map, store and model the essential elements and their interactions in complex neurobiological systems, particularly the human brain, the most complex system in nature. Due to a rapid expansion of neuroimaging data sets' size and complexity, it is essential to propose and adopt new empirical tools to track dynamic patterns between neurons and brain areas and create comprehensive maps. In recent years, there is a rapid growth of scientific interest in moving functional neuroimaging analysis beyond simplified group or time‐averaged approaches and sophisticated algorithms that can capture the time‐varying properties of functional connectivity. We describe algorithms and network metrics that can capture the dynamic evolution of functional connectivity under this perspective. We adopt the word ‘chronnectome’ (integration of the Greek word ‘Chronos’, which means time, and connectome) to describe this specific branch of network neuroscience that explores how mutually informed brain activity correlates across time and brain space in a functional way. We also describe how good temporal mining of temporally evolved dynamic functional networks could give rise to the detection of specific brain states over which our brain evolved. This characteristic supports our complex human mind. The temporal evolution of these brain states and well‐known network metrics could give rise to new analytic trends. Functional brain networks could also increase the multi‐faced nature of the dynamic networks revealing complementary information. Finally, we describe a python module (https://github.com/makism/dyconnmap) which accompanies this article and contains a collection of dynamic complex network analytics and measures and demonstrates its great promise for the study of a healthy subject's repeated fMRI scans.

The impact of graph construction scheme and community detection algorithm on the repeatability of community and hub identification in structural brain networks

A critical question in network neuroscience is how nodes cluster together to form communities, to form the mesoscale organisation of the brain. Various algorithms have been proposed for identifying such communities, each identifying different communities within the same network. Here, (using test–retest data from the Human Connectome Project), the repeatability of thirty‐three community detection algorithms, each paired with seven different graph construction schemes were assessed. Repeatability of community partition depended heavily on both the community detection algorithm and graph construction scheme. Hard community detection algorithms (in which each node is assigned to only one community) outperformed soft ones (in which each node can belong to more than one community). The highest repeatability was observed for the fast multi‐scale community detection algorithm paired with a graph construction scheme that combines nine white matter metrics. This pair also gave the highest similarity between representative group community affiliation and individual community affiliation. Connector hubs had higher repeatability than provincial hubs. Our results provide a workflow for repeatable identification of structural brain networks communities, based on the optimal pairing of community detection algorithm and graph construction scheme.

Predicting MEG resting-state functional connectivity from microstructural information

Understanding how human brain microstructure influences functional connectivity is an important endeavor. In this work, magnetic resonance imaging data from 90 healthy participants were used to calculate structural connectivity matrices using the streamline count, fractional anisotropy, radial diffusivity, and a myelin measure (derived from multicomponent relaxometry) to assign connection strength. Unweighted binarized structural connectivity matrices were also constructed. Magnetoencephalography resting-state data from those participants were used to calculate functional connectivity matrices, via correlations of the Hilbert envelopes of beamformer time series in the delta, theta, alpha, and beta frequency bands. Nonnegative matrix factorization was performed to identify the components of the functional connectivity. Shortest path length and search-information analyses of the structural connectomes were used to predict functional connectivity patterns for each participant. The microstructure-informed algorithms predicted the components of the functional connectivity more accurately than they predicted the total functional connectivity. This provides a methodology to understand functional mechanisms better. The shortest path length algorithm exhibited the highest prediction accuracy. Of the weights of the structural connectivity matrices, the streamline count and the myelin measure gave the most accurate predictions, while the fractional anisotropy performed poorly. Overall, different structural metrics paint very different pictures of the structural connectome and its relationship to functional connectivity.

Mapping Structural Connectivity Using Diffusion MRI: Challenges and Opportunities

Diffusion MRI-based tractography is the most commonly-used technique when inferring the structural brain connectome, i.e., the comprehensive map of the connections in the brain. The utility of graph theory-a powerful mathematical approach for modeling complex network systems-for analyzing tractography-based connectomes brings important opportunities to interrogate connectome data, providing novel insights into the connectivity patterns and topological characteristics of brain structural networks. When applying this framework, however, there are challenges, particularly regarding methodological and biological plausibility. This article describes the challenges surrounding quantitative tractography and potential solutions. In addition, challenges related to the calculation of global network metrics based on graph theory are discussed.Evidence Level: 5Technical Efficacy: Stage 1.

Tractography density affects whole-brain structural architecture and resting-state dynamical modeling

Dynamical modeling of the resting-state brain dynamics essentially relies on the empirical neuroimaging data utilized for the model derivation and validation. There is however still no standardized data processing for magnetic resonance imaging pipelines and the structural and functional connectomes involved in the models. In this study, we thus address how the parameters of diffusion-weighted data processing for structural connectivity (SC) can influence the validation results of the whole-brain mathematical models informed by SC. For this, we introduce a set of simulation conditions including the varying number of total streamlines of the whole-brain tractography (WBT) used for extraction of SC, cortical parcellations based on functional and anatomical brain properties and distinct model fitting modalities. The main objective of this study is to explore how the quality of the model validation can vary across the considered simulation conditions. We observed that the graph-theoretical network properties of structural connectome can be affected by varying tractography density and strongly relate to the model performance. We also found that the optimal number of the total streamlines of WBT can vary for different brain atlases. Consequently, we suggest a way how to improve the model performance based on the network properties and the optimal parameter configurations from multiple WBT conditions. Furthermore, the population of subjects can be stratified into subgroups with divergent behaviors induced by the varying WBT density such that different recommendations can be made with respect to the data processing for individual subjects and brain parcellations.

The R1-weighted connectome: complementing brain networks with a myelin-sensitive measure

Myelin plays a crucial role in how well information travels between brain regions. Complementing the structural connectome, obtained with diffusion MRI tractography, with a myelin-sensitive measure could result in a more complete model of structural brain connectivity and give better insight into white-matter myeloarchitecture. In this work we weight the connectome by the longitudinal relaxation rate (R1), a measure sensitive to myelin, and then we assess its added value by comparing it with connectomes weighted by the number of streamlines (NOS). Our analysis reveals differences between the two connectomes both in the distribution of their weights and the modular organization. Additionally, the rank-based analysis shows that R1 can be used to separate transmodal regions (responsible for higher-order functions) from unimodal regions (responsible for low-order functions). Overall, the R1-weighted connectome provides a different perspective on structural connectivity taking into account white matter myeloarchitecture.

In the present work, we show that by using a myelin-sensitive measure we can complement the diffusion MRI-based connectivity and provide a different picture of the brain organization. We show that the R1-weighted average distribution does not follow the same trend as the number of streamlines strength distribution, and the two connectomes exhibit different modular organization. We also show that unimodal cortical regions tend to be connected by more streamlines, but the connections exhibit a lower R1-weighted average, while the transmodal regions have higher R1-weighted average but fewer streamlines. This could imply that the unimodal regions require more connections with lower myelination, whereas the transmodal regions rely on connections with higher myelination.

Combining network topology and information theory to construct representative brain networks

Network neuroscience employs graph theory to investigate the human brain as a complex network, and derive generalizable insights about the brain's network properties. However, graph-theoretical results obtained from network construction pipelines that produce idiosyncratic networks may not generalize when alternative pipelines are employed. This issue is especially pressing because a wide variety of network construction pipelines have been employed in the human network neuroscience literature, making comparisons between studies problematic. Here, we investigate how to produce networks that are maximally representative of the broader set of brain networks obtained from the same neuroimaging data. We do so by minimizing an information-theoretic measure of divergence between network topologies, known as the portrait divergence. Based on functional and diffusion MRI data from the Human Connectome Project, we consider anatomical, functional, and multimodal parcellations at three different scales, and 48 distinct ways of defining network edges. We show that the highest representativeness can be obtained by using parcellations in the order of 200 regions and filtering functional networks based on efficiency-cost optimization-though suitable alternatives are also highlighted. Overall, we identify specific node definition and thresholding procedures that neuroscientists can follow in order to derive representative networks from their human neuroimaging data.

Network diffusion modeling predicts neurodegeneration in traumatic brain injury

Objective

Traumatic brain injury (TBI) is a heterogeneous disease with multiple neurological deficits that evolve over time. It is also associated with an increased incidence of neurodegenerative diseases. Accordingly, clinicians need better tools to predict a patient’s long‐term prognosis.

Methods

Diffusion‐weighted and anatomical MRI data were collected from 17 adolescents (mean age = 15y8mo) with moderate‐to‐severe TBI and 19 healthy controls. Using a network diffusion model (NDM), we examined the effect of progressive deafferentation and gray matter thinning in young TBI patients. Moreover, using a novel automated inference method, we identified several injury epicenters in order to determine the neural degenerative patterns in each TBI patient.

Results

We were able to identify the subject‐specific patterns of degeneration in each patient. In particular, the hippocampus, temporal cortices, and striatum were frequently found to be the epicenters of degeneration across the TBI patients. Orthogonal transformation of the predicted degeneration, using principal component analysis, identified distinct spatial components in the temporal–hippocampal network and the cortico‐striatal network, confirming the vulnerability of these networks to injury. The NDM model, best predictive of the degeneration, was significantly correlated with time since injury, indicating that NDM can potentially capture the pathological progression in the chronic phase of TBI.

Interpretation

These findings suggest that network spread may help explain patterns of distant gray matter thinning, which would be consistent with Wallerian degeneration of the white matter connections (i.e., “diaschisis”) from diffuse axonal injuries and multifocal contusive injuries, and the neurodegenerative patterns of abnormal protein aggregation and transmission, which are hallmarks of brain changes in TBI. NDM approaches could provide highly subject‐specific biomarkers relevant for disease monitoring and personalized therapies in TBI.