Integrative Bayesian Analysis of Brain Functional Networks Incorporating Anatomical Knowledge

TL;DR

This paper introduces a hierarchical Bayesian Gaussian graphical model that integrates anatomical knowledge to improve the accuracy and reproducibility of brain functional network estimation from multimodal imaging data.

Contribution

The paper presents a novel Bayesian modeling approach that incorporates anatomical pathways into functional network estimation, robust to inaccuracies in anatomical data.

Findings

Enhanced accuracy in functional connectivity estimation.

Greater reproducibility of network metrics across sessions.

Effective identification of structure-supported functional connections.

Abstract

Recently, there has been increased interest in fusing multimodal imaging to better understand brain organization. Specifically, accounting for knowledge of anatomical pathways connecting brain regions should lead to desirable outcomes such as increased accuracy in functional brain network estimates and greater reproducibility of topological features across scanning sessions. Despite the clear merits, major challenges persist in integrative analyses including an incomplete understanding of the structure-function relationship and inaccuracies in mapping anatomical structures due to deficiencies in existing imaging technology. Clearly advanced network modeling tools are needed to appropriately incorporate anatomical structure in constructing brain functional networks. We propose a hierarchical Bayesian Gaussian graphical modeling approach that estimates the functional networks via sparse precision matrices whose degree of edge-specific shrinkage is informed by anatomical structure and an independent baseline component. The approach flexibly identifies functional connections supported by structural connectivity knowledge. This enables robust brain network estimation even in the presence of mis-specified anatomical knowledge, while accommodating heterogeneity in the structure-function relationship. We implement the approach via an efficient optimization algorithm yielding maximum a posteriori estimates. Extensive numerical studies reveal the clear advantages of our approach over competing methods in accurately estimating brain functional connectivity, even when the anatomical knowledge is mis-specified. An application of the approach to the Philadelphia Neurodevelopmental Cohort (PNC) study reveals gender based connectivity differences across multiple age groups, and higher reproducibility in the estimation of network metrics compared to alternative methods.