An automated pipeline for obtaining labeled ICA-templates corresponding to functional brain systems

The complexity of our actions and thinking is likely reflected in functional brain networks. Independent component analysis (ICA) is a popular data-driven method to compute group differences between such networks. A common way to investigate network differences is based on ICA maps which are generated from study-specific samples. However, this approach limits the generalizability and reproducibility of the results. Alternatively, network ICA templates can be used, but up to date, few such templates exist and are limited in terms of the functional systems they cover. Here, we propose a simple two-step procedure to obtain ICA-templates corresponding to functional brain systems of the researcher's choice: In step 1, the functional system of interest needs to be defined by means of a statistical parameter map (input), which one can generate with open-source software such as NeuroSynth or BrainMap. In step 2, that map is correlated to group-ICA maps provided by the Human Connectome Project (HCP), which is based on a large sample size and uses high quality and standardized acquisition procedures. The HCP-provided ICA-map with the highest correlation to the input map is then used as an ICA template representing the functional system of interest, for example, for subsequent analyses such as dual regression. We provide a toolbox to complete step 2 of the suggested procedure and demonstrate the usage of our pipeline by producing an ICA templates that corresponds to "motor function" and nine additional brain functional systems resulting in an ICA maps with excellent alignment with the gray matter/white matter boundaries of the brain. Our toolbox generates data in two different file formats: volumetric-based (NIFTI) and combined surface/volumetric files (CIFTI). Compared to 10 existing templates, our procedure output component maps with systematically stronger contribution of gray matter to the ICA z-values compared to white matter voxels in 9/10 cases by at least a factor of 2. The toolbox allows users to investigate functional networks of interest, which will enhance interpretability, reproducibility, and standardization of research investigating functional brain networks.