Mapping morphological cortical networks with joint probability distributions from multiple morphological features.

• MCNs were constructed by applying Jensen-Shannon divergence to regional joint probability distributions derived from different morphological features. • MCNs showed efficient small-word organization and high test-retest reliability, aligned with cortical cytoarchitectonics and symmetry and axonal connectivity, and explained inter-individual variance in behavior and cognition. • MCNs built based on multiple features performed better than those derived from different single features in most aspects mentioned above. Morphological features sourced from structural magnetic resonance imaging can be used to infer human brain connectivity. Although integrating different morphological features may theoretically be beneficial for obtaining more precise morphological connectivity networks (MCNs), the empirical evidence to support this supposition is scarce. Moreover, the incorporation of different morphological features remains an open question. In this study, we proposed a method to construct cortical MCNs based on multiple morphological features. Specifically, we adopted a multi-dimensional kernel density estimation algorithm to fit regional joint probability distributions (PDs) from different combinations of four morphological features, and estimated inter-regional similarity in the joint PDs via Jensen-Shannon divergence. We evaluated the method by comparing the resultant MCNs with those built based on different single morphological features in terms of topological organization, test-retest reliability, biological plausibility, and behavioral and cognitive relevance. We found that, compared to MCNs built based on different single morphological features, MCNs derived from multiple morphological features displayed less segregated, but more integrated network architecture and different hubs, had higher test-retest reliability, encompassed larger proportions of inter-hemispheric edges and edges between brain regions within the same cytoarchitectonic class, and explained more inter-individual variance in behavior and cognition. These findings were largely reproducible when different brain atlases were used for cortical parcellation. Further analysis of macaque MCNs revealed weak, but significant correlations with axonal connectivity from tract-tracing, independent of the number of morphological features. Altogether, this paper proposes a new method for integrating different morphological features, which will be beneficial for constructing MCNs. [ABSTRACT FROM AUTHOR]