Group-Level Ranking-Based Hubness Analysis of Human Brain Connectome Reveals Significant Interhemispheric Asymmetry and Intraparcel Heterogeneities

Objective: Graph theory applications are commonly used in connectomics research to better understand connectivity architecture and characterize its role in cognition, behavior and disease conditions. One of the numerous open questions in the field is how to represent inter-individual differences with graph theoretical methods to make inferences for the population. Here, we proposed and tested a simple intuitive method that is based on finding the correlation between the rank-ordering of nodes within each connectome with respect to a given metric to quantify the differences/similarities between different connectomes. Methods: We used the diffusion imaging data of the entire HCP-1065 dataset of the Human Connectome Project (HCP) (n = 1,065 subjects). A customized cortical subparcellation of HCP-MMP atlas (360 parcels) (yielding a total of 1,598 ROIs) was used to generate connectivity matrices. Six graph measures including degree, strength, coreness, betweenness, closeness, and an overall “hubness” measure combining all five were studied. Group-level ranking-based aggregation method (“measure-then-aggregate”) was used to investigate network properties on population level. Results: Measure-then-aggregate technique was shown to represent population better than commonly used aggregate-then-measure technique (overall r(s): 0.7 vs 0.5). Hubness measure was shown to highly correlate with all five graph measures (r(s): 0.88–0.99). Minimum sample size required for optimal representation of population was found to be 50 to 100 subjects. Network analysis revealed a widely distributed set of cortical hubs on both hemispheres. Although highly-connected hub clusters had similar distribution between two hemispheres, average ranking values of homologous parcels of two hemispheres were significantly different in 71% of all cortical parcels on group-level. Conclusion: In this study, we provided experimental evidence for the robustness, limits and applicability of a novel group-level ranking-based hubness analysis technique. Graph-based analysis of large HCP dataset using this new technique revealed striking hemispheric asymmetry and intraparcel heterogeneities in the structural connectivity of the human brain.