BrainVoyager QX v2.8
Functionally Informed Cortex-Based Alignment
While cortex-based alignment (CBA) achieves high corresondance of macro-anatomical landmarks (sulci, gyri) across subjects, functional areas still exhibit residual variability. This residual variability contributes to a decrease in statistical power of any multi-subject whole-cortex fMRI analysis. In order to address this issue, functionally informed CBA (fCBA) uses homologue functional information to drive alignment in addition to the minimization of the macroanatomical variability across subjects. By using proper weighting of anatomical and functional forces, fCBA creates a group aligned space with strong spatial correspondence of both structure and function while maintaining a high level of anatomical fidelity. In case that homologue ROIs are available from functional localizer runs for all included subjects, fCBA allows to perform a single main analysis that integrates whole-cortex and ROI-based analysis as compared to the usual approach that involves running two separate analyses, i.e. one analysis for aligned ROIs and one separate whole-cortex analysis. In both cases, a separate definition of homologue regions based on a separate localizer is needed requiring a second (fCBA) or third (CBA, Talairach) analysis.
It is important to note that fCBA should not be used when the goal of a study is the assessment of functional variability of homologue regions across subjects (Frost & Goebel, 2012). FCBA is recommended in cases where one wants to optimize group-level statistical analyses in the same way as one usually would use ROI-aligned multi-subject analysis (Frost & Goebel, submitted).
It has been recently demonstrated that, after curvature driven cortex based alignment (CBA), macroantomical variability between subjects is substantially reduced and indeed almost completely removed for major sulci. Furthermore as a consequence of increasing correspondence of anatomy, a concomitant increase in the overlap of a set of investigated functional areas across subjects was observed (Frost and Goebel, 2012). The gain in overlap after cortex based alignment was not uniform across all function areas, suggesting that some functional areas are strongly related to macro-anatomy e.g. the frontal eye fields (Paus, 1996) whereas others are not as tightly bound e.g. the fusiform face area (FFA, but see Weiner & Grillspector (2010) and discussion in Frost & Goebel (2012)).
In order to overcome the problem that there remain differences in the spatial location of functional areas after registration to a standard template, and even after advanced alignment schemes, a region of interest (ROI) approach is often utilised. Restricting analysis to a few circumscribed, independently localised regions allows to carry out group comparisons of functionally homologous regions across subjects with high statistical power and substantially reduces the correction needed for multiple comparisons as is needed for the many thousands of voxels of a whole brain analysis. However this approach ignores activity in areas outside the ROIs which would be revealed by a whole-brain analysis. While one could perform both kinds of analyses separately (as is indeed reported in many published studies), an integrated approach combining ROI and whole brain analysis with maximal correspondence across subjects would be desirable. Here we propose an integrated alignment technique that incorporates the spatial specificity and statistical power of a region of interest analysis, with the benefits of a whole brain analysis (Frost & Goebel, submitted). To indicate that the whole-cortex alignment incorporates functional information from localised ROIs, we refer to the developed procedure as functionally informed cortex based alignment (fCBA). Due to the fact that a purely macro-anatomical alignment is not capable of successfully bringing functional regions fully into spatial correspondence, fCBA integrates both curvature (anatomical information) and functionally defined regions to produce a group brain, which is well aligned in terms of anatomy and function.
In order to apply the general gradient-based alignment procedure to POIs (defined as a contigous region of vertices in a cortex mesh), the functionally labelled areas are smoothed with a Gaussian kernel (approximated by iterative neighborhood smoothing) providing a gradient into the area’s neighbourhood to be followed during the alignment procedure maximising overlap across subjects. The integrated alignment process follows the same stages of standard CBA including a coarse-to-fine alignment of curvature information to a dynamic group average. The (smoothed) POI information is simply added as another constraint operating simultaneously with the curvature maps. The relative strengths of macro-anatomical (curvature-based) vs functional (smoothed POI based) forces is controlled by relative weighting parameters (see below).
The execution of functionally informed CBA proceeds largely identical to a standard CBA process since fCBA merely adds additional constraints to the standard cortex-based alignment routine. Adding functional constraints to a CBA task is specified in the Cortex-Based Group Alignment Parameters dialog (see snapshot above) that can be invoked by clicking the Options button in the Align Group tab of the main Cortex-Based Alignment dialog. Below the CBA parameters table that contains (default) values for CBA parameters, a section is available for specifying fCBA parameters that is initially disabled (greyed out).
To turn on fCBA, check the Include functional constraints - fCBA option. This will enable the options below that are specific for fCBA. Furthermore a new column, "fCBA Weight" is added in the general CBA parameters table in the upper section of the dialog (see snapshot above). The weights in the added column specify how strong provided functional information (weight wF) should drive the alignment relative to macro-anatomical (curvature) information (weight wF) with the constraint: wA + wF = 1. A value of "1.0" for the functional weight indicates that only functional information drives the alignment whereas a value of "0.0" indicates that only macro-anatomical information is used. The default values start with modest functional influence at the beginning of the CBA routine (0.2 for smooth level 4) and increase during the coarse-to-fine alignment scheme (0.4 for smooth level 3, 0.6 for smooth level 2 and 0.8 for smooth level 1).
In addition to functional forces during the 4 major stages of CBA, an additional 5th alignment stage is added that continues from the moment when the stage defined in the first row of the CBA parameter table (smooth level 1) has been completed. The parameters of this additional stage can be specified in the single-row Extra-fCBA table that is located below the Include functional alignment option (see snapshot above). The reason for adding this step is to allow extra strong effects of functional constraints, i.e. if one wants to first run a rather standard macro-anatomical CBA and to add the fCBA step afterwards or in order to force alignment of homologue functionally defined regions tolerating some macro-anatomical misalignment. As has been turned out (Frost & Goebel, submitted), the balance of forces during CBA combined with the possibility of an extra step provides a rather fine-grained control to achieve desired results (see below).
Besides setting parameters that control functional alignment forces, the homologue regions that should be aligned need to be specified in subject-specific space. This is done by providing a POI file for each subject containing the definition of corresponding areas. If more than one area is defined, the different regions need to appear in the same order in the POI file. If, for example, a localizer run has been performed to map the fusiform face area (FFA) and the parahippocampal place area (PPA), a single .poi file with the FFA POI and the PPA POI needs to be created and stored in the same order. POIs can be created in the usual way, e.g. by running a MTC-GLM on individual SPH LH/RH cortex meshes followed by overlaying a "FFA vs PPA" contrast; after adjusting the threshold of the contrast t map appropriately, a target region (e.g. FFA) needs to be selected on the map as a contigous region by performing a CTRL-click (CMD-click on Mac); the marked region can then be converted into a Patch-Of-Interest (POI) by using the context menu entry Define Patch-Of-Interest (POI); if the order of the defined areas is not the same in each subject, the Up and Down arrow buttons in the Patch-Of-Interest Analysis dialog can be used to rearrange the POIs in a consistent order. The created subject-specific POI files need to be stored in the same folder that contains also the macro-anatomical alignment information, i.e. the curvature maps. Furthermore, the names of the POI files need to be identical except for the intial substring that must contain the subject ID - the same ID that is also used for the curvature maps. If, for example, a subject's curvature map has the name "AC_LH_SPH_CURVATURE.smp", the poi file need to also start with the "AC_" substring. In a simple example data set, three pois (FFA, hMT, FEF) were defined using appropriate localizer runs and the .poi files for three of the subjects included in the example were named "AC_FFA_MT_FEF.poi", "AZ_FFA_MT_FEF.poi" and "CG_FFA_MT_FEF.poi". In order to tell the program that these are the .poi names of the included subjects, only the file name sub-string following the subject-ID needs to be entered in the Name of pois text field without file extension (see lower red ellipse in the snapshot above); with this name "template", the program figures out the names for the individual subjects using the subject-ID from the curvature maps and it appends the extension ".poi". Furthermore, the number of POIs defined in each POI file needs to be set in the No. of functional maps/pois spin box (see upper red ellipse in the snapshot above). Since in our example there are 3 POIs, the value entered is "3". Note that the Put in one map option should always be turned on (default) at present; this option will create one surface map file containing all smoothed POIs at once per subject as opposed to creating a separate surface map per POI. Finally the No. of smooth steps value (for explanation of its function see below) can be set but it is recommended to keep the default value of "50". After all settings have been made, the Cortex-Based Group Alignment Parameters dialog can be closed by clicking the OK button.
When now a dynamic group average CBA process is started (with proper specification of the curvature maps of the corresponding hemispheres for the included subjects as usual), the functionally defined constraints will be effective during alignment. Since the alignment process needs gradient information to work successfully, the POIs are first spatially smoothed before CBA begins so that they are represented with a peak value in the middle that gradually falls off if one moves away from the peak in any direction. In order to perform this smoothing operation, POIs are first converted into surface maps internally and saved to disk. In the snapshot above, the smoothed FFA area is shown for the three included subjects in three different colors. The amount of iterative neighborhood smoothing is specified by the No. of smooth steps value (see above). The performed smoothing spreads the binary region occupied by the original POIs into the neighborhood increasing the chance that homologue regions across subjects (e.g. the FFA POIs of all subjects) will overlap; note that functional alignment forces only have an effect if homologue POIs overlap across subjects.
After fCBA has been completed, one can inspect the achieved POI (and macroa-anatomical) alignment in various ways. First, the program automatically puts the smoothed POIs in a final surface map that is displayed at the end of fCBA superimposed on a standard sphere or directly on the averaged folded group brain (if the individual "RECOSM_SPH" meshes can be located in the subject-specific sub-folders). The surface map is also stored to disk automatically and contains 4 macro-anatomical and 1 POI sub-map for each of the included subjects. In the overlaid surface map, only aligned curvature sub-maps are selected as default. In order to visualize the POI sub-maps, they need to be selected in the sub-maps list. Note that in order to visualize multiple maps simultaeneously, the Multiple selections option in the Advanced tab of the Surface Maps dialog needs to be turned on. In the snapshot above (right side), the POI sub-maps have been selected for each of the three subjects showing the nearly perfect alignment of the smoothed FFA POIs after fCBA. To visualize the different subjects better, the default orange sub-map color has been changed to a red, green and blue color value.
While the surface maps provide already insight in the quality of POI alignment, the alignment of the original POI shapes should be inspected by aligning the POIs using the standard CBA POI transformation tool. The snapshot above shows teh resulting "ALIGNED" POIs on the right side and the location of the POIs prior to fCBA alignment on the left side (both are shown on the average folded group brain for comparison). As is evident, the POIs are indeed located at the same region after fCBA.
The fCBA procedure offers a powerful means to align homologue functional regions across subjects. While the default functional weighting values provide good results, it may be necessary to adjust the values for specific purposes. While functional constraints should align homologue POIs as good as possible, it is also desirable to achieve a good as possible macro-anatomical (curvature-based) alignment. The provided default values are set to optimize POI overlap tolerating visible decrease of macro-anatomical alignment as compared to running pure macro-anatomical CBA. Recently, optimal functional constraint parameters have been reported that appear well suited to integrate both goals in a near-optimal way (Frost & Goebel, submitted). While not always leading to a perfect POI alignment (i.e. POIs do not exhibit complete overlap after fCBA), parameters are suggested that lead to near-optimal POI alignment with minimal macro-anatomical detoriations as compared to pure macro-anatomical CBA. The functional weights for this "balanced" weighting for the four coarse-to-fine stages were 0.2, 0.4, 0.6, 0.9 with 400, 300, 300 and 1000 iterations with 4, 4, 3 and 10 number of updates, respectively. The recommended parameters for the extra step were 0.2 functional weighting, 400 iterations and 4 updates. While the provided default values and the values suggested here are good candidates to either push POI alignment or to go for a more balanced functional-macro-anatomical alignment, it should be noted that some experimentation with the parameters might be needed, i.e. by running fCBA repeatedly with different parameter sets.
Copyright © 2014 Rainer Goebel. All rights reserved.