Parcellation of fMRI datasets with ICA and PLS--a data driven approach

One Size Does Not Fit All: Methodological Considerations for Brain-Based Predictive Modeling in Psychiatry

Psychiatric illnesses are heterogeneous in nature. No illness manifests in the same way across individuals, and no two patients with a shared diagnosis exhibit identical symptom profiles. Over the last several decades, group-level analyses of in vivo neuroimaging data have led to fundamental advances in our understanding of the neurobiology of psychiatric illnesses. More recently, access to computational resources and large, publicly available datasets alongside the rise of predictive modeling and precision medicine approaches have facilitated the study of psychiatric illnesses at an individual level. Data-driven machine learning analyses can be applied to identify disease-relevant biological subtypes, predict individual symptom profiles, and recommend personalized therapeutic interventions. However, when developing these predictive models, methodological choices must be carefully considered to ensure accurate, robust, and interpretable results. Choices pertaining to algorithms, neuroimaging modalities and states, data transformation, phenotypes, parcellations, sample sizes, and populations we are specifically studying can influence model performance. Here, we review applications of neuroimaging-based machine learning models to study psychiatric illnesses and discuss the effects of different methodological choices on model performance. An understanding of these effects is crucial for the proper implementation of predictive models in psychiatry and will facilitate more accurate diagnoses, prognoses, and therapeutics.

Evaluation of node-inhomogeneity effects on the functional brain network properties using an anatomy-constrained hierarchical brain parcellation

To investigate functional brain networks, many graph-theoretical studies have defined nodes in a graph using an anatomical atlas with about a hundred partitions. Although use of anatomical node definition is popular due to its convenience, functional inhomogeneity within each node may lead to bias or systematic errors in the graph analysis. The current study was aimed to show functional inhomogeneity of a node defined by an anatomical atlas and to show its effects on the graph topology. For this purpose, we compared functional connectivity defined using 138 resting state fMRI data among 90 cerebral nodes from the automated anatomical labeling (AAL), which is an anatomical atlas, and among 372 cerebral nodes defined using a functional connectivity-based atlas as a ground truth, which was obtained using anatomy-constrained hierarchical modularity optimization algorithm (AHMO) that we proposed to evaluate the graph properties for anatomically defined nodes. We found that functional inhomogeneity in the anatomical parcellation induced significant biases in estimating both functional connectivity and graph-theoretical network properties. We also found very high linearity in major global network properties and nodal strength at all brain regions between anatomical atlas and functional atlas with reasonable network-forming thresholds for graph construction. However, some nodal properties such as betweenness centrality did not show significant linearity in some regions. The current study suggests that the use of anatomical atlas may be biased due to its inhomogeneity, but may generally be used in most neuroimaging studies when a single atlas is used for analysis.

Functional connectivity-based identification of subdivisions of the basal ganglia and thalamus using multilevel independent component analysis of resting state fMRI

This study aimed to identify subunits of the basal ganglia and thalamus and to investigate the functional connectivity among these anatomically segregated subdivisions and the cerebral cortex in healthy subjects. For this purpose, we introduced multilevel independent component analysis (ICA) of the resting-state functional magnetic resonance imaging (fMRI). After applying ICA to the whole brain gray matter, we applied second-level ICA restrictively to the basal ganglia and the thalamus area to identify discrete functional subunits of those regions. As a result, the basal ganglia and the thalamus were parcelled into 31 functional subdivisions according to their temporal activity patterns. The extracted parcels showed functional network connectivity between hemispheres, between subdivisions of the basal ganglia and thalamus, and between the extracted subdivisions and cerebral functional components. Grossly, these findings correspond to cortico-striato-thalamo-cortical circuits in the brain. This study also showed the utility of multilevel ICA of resting state fMRI in brain network research.