A 3D Brain Geometry Toolkit for Multisite Neuroimaging Analysis

Compared to traditional gross volumetrics, surface-based models provide greater spatial precision for understanding brain alterations related to developmental, neurological, and psychiatric disorders. Large-scale brain initiatives are combining data from around the world to discover and improve illness-related brain markers. Here, we present a toolkit for 3D brain geometry analysis aimed at addressing key challenges facing large-scale neuroimaging studies. Our framework incorporates scalable methods for multisite data integration, site-specific confound correction, accelerated statistical modeling, interpretable machine learning, and interactive results visualization. The toolkit was tested on data from 21 independently collected study samples participating in the ENIGMA Bipolar Disorder Working Group (N = 3,373). Compared to traditional volume features, we show how subcortical shape measures can be combined across study sites to capture spatially complex differences between diagnostic groups and associations with common treatments. Statistical modeling was accelerated using the Fast and Efficient Mixed-Effects Algorithm (FEMA) and achieved a 16-fold reduction in computation time compared to traditional approaches. Machine learning models showed shape features may provide greater predictive performance over traditional volumes for both diagnostic and treatment prediction tasks, with interpretable weight maps providing insights into the local features driving model performance.