ALE meta-analysis workflows via the BrainMap database: progress towards a probabilistic functional brain atlas
1
Research Imaging Center, University of Texas Health Science Center, San Antonio, TX, USA
2
Institute for Neuroscience and Medicine (INM – 2), Research Center Jülich, Jülich, Germany
3
Department of Psychiatry and Psychotherapy, RWTH Aachen University, Aachen, Germany
4
Jülich Aachen Research Alliance–Translational Brain Medicine, Germany
5
Department of Psychiatry and Human Behavior, University of California, Irvine, CA, USA
6
Neuroscience Institute, Scott and White Memorial Hospital, Temple, TX, USA
With the ever-increasing number of studies in human functional brain mapping, an abundance of data has been generated that is ready to be synthesized and modeled on a large scale. The BrainMap database archives peak coordinates from published neuroimaging studies, along with the corresponding metadata that summarize the experimental design. BrainMap was designed to facilitate quantitative meta-analysis of neuroimaging results reported in the literature and supports the use of the activation likelihood estimation (ALE) method. In this paper, we present a discussion of the potential analyses that are possible using the BrainMap database and coordinate-based ALE meta-analyses, along with some examples of how these tools can be applied to create a probabilistic atlas and ontological system of describing function–structure correspondences.
Over the last three decades, neuroimaging research has produced an enormous amount of data localizing the neural effects of specific mental operations in both healthy and diseased populations. Community-wide standards of spatial normalization and the reporting of peak activation locations in stereotactic coordinates allow researchers to compare results across studies when the primary data are unavailable or difficult to obtain. Due to the nearly universal adherence to these standards, the BrainMap project was designed to create tools for large-scale data mining and meta-analysis of the brain mapping literature (Fox and Lancaster, 2002
; Laird et al., 2005a
).
BrainMap is a community accessible database
1
that allows a user to relate behavioral functions to specific brain locations through retrieval and visualization of peak coordinates and their associated metadata. These metadata allow each coordinate to be linked with how the observed activation was experimentally derived, a formulation that lends itself to very rich data mining. BrainMap was originally conceived by Peter Fox in 1987 and received its original funding from the James S. McDonnell Foundation (1988–1990). Continued BrainMap development was funded by the Office of Naval Research (1991–1992), the EJLB Foundation (1992–1996), and the National Library of Medicine (2000–2003). BrainMap is currently funded by the Human Brain Project of the National Institute of Mental Health.
There are three desktop applications (Scribe, Sleuth, and GingerALE) and one web application (BrainMapWeb) that allow interaction with the BrainMap database, all of which are coded in Java. The desktop applications run in the Java Runtime Environment on Macintosh, Windows, Linux, and Unix operating systems, while the web application uses Java server-side technologies. Scribe
2
is used to code papers for entry into BrainMap. Peer-reviewed publications can be submitted to the database by the original authors (uncommon) or by investigators performing a meta-analysis (very common). Most data fields have candidate responses presented in scrollable lists. Coordinate tables of peak locations can be imported from a tab-delimited file or entered by hand. Upon insertion into the database, each x,y,z, coordinate is assigned an anatomical location using the Talairach Daemon
3
(Lancaster et al., 2000
). All entries are reviewed for quality control by BrainMap staff and faculty before being entered into the database to ensure the accuracy and consistency of coding. The Sleuth application
4
allows a user to search the BrainMap database and retrieve data, which can then be filtered and visualized on a standard Talairach atlas image. Specific locations may be searched for according to user-defined, coordinate-based regions of interest (defined by Talairach or MNI coordinates) or anatomical labels from the Talairach Daemon nomenclature. BrainMap queries may also be implemented via an internet browser using BrainMapWeb
5
, which includes query functions that are similar to those of Sleuth, but lack 3D brain visualizations. Data view and manipulation capabilities are much more restricted than in Sleuth, which can output data in publication-ready graphics, text files of Talairach or MNI coordinates, or workspaces that specify search rules and filters for meta-analyses. By archiving coordinates of activation locations rather than raw image data, BrainMap focuses on discoveries derived from coordinate-based meta-analyses of functional neuroimaging data. The last BrainMap application, GingerALE
6
, is used for performing activation likelihood estimation (ALE) meta-analyses on sets of coordinates extracted from the database in Talairach or MNI space. Scribe, Sleuth, and GingerALE are closely integrated to transition seamlessly from database submission to search refinement to meta-analysis results (Figure 1
).
Figure 1. BrainMap software applications. There are three desktop applications that allow interaction with the BrainMap database: Scribe (data entry), Sleuth (search and retrieval), and GingerALE (meta-analysis). These Java applications are freely downloadable from the BrainMap website (http://brainmap.org
). Database statistics are current as of 04/01/09.
To summarize the experimental design and results of a published study for inclusion in the database, BrainMap utilizes a rigorous taxonomy that is composed chiefly of structured keywords. BrainMap entries include descriptions on the scanned subjects and experimental conditions, including the presented stimuli, instructions, and responses. A complete metadata listing can be seen in Figure 2
. To facilitate meta-analysis, several hierarchically structured keywords have been developed that categorize the nature of each experimental contrast to allow rapid, comprehensive retrieval of results. “Context” broadly categorizes the purpose of the work; for example, normal mapping (the comparison of different experimental conditions in a group of healthy subjects), age effects, disease effects, or drug effects. “Behavioral Domain” classifies the research in terms of the neural systems studied according to six main categories and their related subcategories: cognition, action, perception, emotion, interoception, or pharmacology (Figure 3
). “Paradigm Class” categorizes the challenge presented, preferably in the jargon of the field, such as anti-saccades, Stroop, delayed match to sample, or mental rotation tasks (Figure 4
).
Figure 2. The BrainMap coding scheme. The BrainMap metadata coding scheme follows a hierarchy naturally occurring within the brain mapping literature. Every Paper reports experimental results drawn from one or more subject groups, the members of which have been functionally imaged during one or more behavioral Conditions. In BrainMap, an Experiment is defined by the production of a statistical parametric image. From each Experiment one or more functional activations (Locations) are extracted in the form of peak coordinates.
Figure 3. BrainMap behavioral domains. In the BrainMap coding scheme, the behavioral domain is recorded for each statistical contrast, which describes the research in terms of the neural systems studied according to 6 main categories: cognition, action, perception, emotion, interoception, or pharmacology. The yellow background represents that this field lies at the Experiment level.
Figure 4. BrainMap paradigm classes. The paradigm class is recorded for each statistical contrast in BrainMap. This field categorizes the task or challenge presented, preferably in the jargon of the field. BrainMap currently contains 80 paradigm classes. The yellow background represents that this field lies at the Experiment level.
While a given a given paradigm class (e.g., n-back) is often immediately associated with a given behavioral domain (e.g., working memory), this is not always the case. Each experiment must be evaluated based on the conditions contrasted since many comparisons are designed to elicit processing in domains not directly linked to the paradigm class employed. For example, n-back and working memory are appropriate choices when comparing an n-back condition to a control condition, but if two n-back conditions are contrasted that differ only in the modality of stimulus presentation (e.g., visual or auditory), then additional perceptual domains should be coded for that experiment. We have found that the context, behavioral domain, and paradigm class represent the three most critical components of a functional neuroimaging study; their orthogonality fully defines and gives contextual meaning to the coordinates archived in the database.
Typically, investigators use citation indexing services such as PubMed to search the literature according to user-defined keywords. This results in the identification of a subset of desired studies that explicitly match the search criteria. BrainMap’s strategy has been to design paradigm class entries in order to pool similar studies, rather than segregate them according to domain-specific keywords. Categories and sub-categories are created and refined only as needed, based on the demands of the literature and the continued development of functional brain imaging. For example, the first study entered into BrainMap that utilized the Wisconsin Card Sorting Test was initially coded under the paradigm class of “Deductive Reasoning”. This was done to classify the study in the same set as other similar papers in the database (i.e., “best fit”), a practice that the classes are rich enough to be useful. Once four papers of the same sub-category are entered into the database, a new class is created and defined (e.g., “Wisconsin Card Sorting Test”). Then all entries matching this new class are manually searched for and updated to reflect the new designation. This procedure requires continuous and labor-intensive maintenance of the database, yet yields a high-quality database full of rich metadata categories and provides an evolving and flexible tool.
Given the sheer volume of neuroimaging data that is currently being produced, it is rapidly becoming overwhelming for an investigator to reconcile new results to those previously published, particularly when studies pertain to different areas of research. Derrfuss and Mar (2009)
estimated that BrainMap contains approximately one-fifth of the relevant published studies, making it the largest coordinate-based database in functional neuroimaging to date. In BrainMap, an ROI of 1 cm3 currently contains an average of 23 experiments, and includes results from 15 paradigm classes. Databases designed to simply retrieve studies reporting activations in proximate locations result in subsequent manual filtering and interpretation. While BrainMap’s data entry procedure is labor-intensive (Laird et al., 2005a
), the depth of the current coding strategy is what provides diverse data mining opportunities and establishes the overall value of the database. BrainMap was structured with the goal of not only retrieving studies returned by regional searches without domain-specific biases, but also allowing the results to be synthesized. Specifically, BrainMap development of neuroinformatics tools focuses on knowledge discovery that is made possible by coordinate-based function–location meta-analysis (Fox et al., 1998
).
Activation likelihood estimation (ALE) is a method of coordinate-based meta-analysis that is supported within the BrainMap software environment. It is a useful tool for integrating the neuroimaging literature wherein consistent regions of activation are identified across a collection of studies. In particular, peak coordinates are collected from studies that share a similar feature of interest, which can be a specific task (e.g., go/no-go) or a more generalized cognitive process (e.g., inhibition). In ALE, coordinates are then modeled with a Gaussian function to accommodate the spatial uncertainty associated with a reported coordinate and are analyzed for where they converge.
Since its introduction (Chein et al., 2002
; Turkeltaub et al., 2002
), ALE has been applied in many aspects of normal brain function (Costafreda et al., 2008
; Decety and Lamm, 2007
; Eickhoff et al., 2006a
; Grosbras et al., 2005
; Soros et al., 2009
; Spreng et al., 2009
), as well as in studies of neuropsychiatric and neurological disorders, such as schizophrenia (Glahn et al., 2005
; Minzenberg et al., 2009
; Ragland et al., 2009
), obsessive-compulsive disorder (Menzies et al., 2008
), depression (Fitzgerald et al., 2008
), and developmental stuttering (Brown et al., 2005
). Recently, ALE has been extended to voxel-based morphometry (Ellison-Wright et al., 2008
; Glahn et al., 2008
; Schroeter et al., 2007
) and diffusion tensor imaging studies (Ellison-Wright and Bullmore, 2009
).
The most interesting ALE applications do not merely merge previous results in a retrospective fashion, but instead generate or test a new hypothesis (Eickhoff et al., 2009a
; Price et al., 2005
), identify a previously unspecified region (Derrfuss et al., 2005
), resolve conflicting views (Laird et al., 2005b
; Petacchi et al., 2005
), or validate a new paradigm (McMillan et al., 2007
). Several studies have used ALE as a preliminary step, followed by an analysis of network co-occurrences (Lancaster et al., 2005
; Neumann et al., 2005
, 2008
; Toro et al., 2008
) or structural equation modeling (Laird et al., 2008
). Each of these novel meta-analytic applications was carried out using typical ALE analysis parameters, and many of them involved the comparison of multiple meta-analyses. For example, Price et al. (2005)
examined the results of a picture naming meta-analysis and found that use of a high-level baseline condition to control for speech production and perceptual processing resulted in increased sensitivity to activation in areas associated with semantic processing, visual-speech integration, and response selection. A prospective fMRI study was performed to test this hypothesis, which subsequently allowed the picture naming system to be decomposed into its perceptual, semantic, and phonological components. In a different application of the ALE method, a meta-analysis was performed on studies in which transcranial magnetic stimulation was applied to left motor cortex (Laird et al., 2008
). The results of this meta-analysis were used to determine the location of regions of interest in a prospective study of TMS/PET data that examined the effective connectivity of the motor system using structural equation modeling. These examples of how meta-analysis results have been applied to guide analyses in newly acquired experimental data demonstrate the power of the ALE method and provide evidence of its efficacy beyond that of a purely retrospective tool.
The ALE method was originally developed and validated by Turkeltaub et al. (2002)
in a meta-analysis of single word reading. BrainMap developers obtained the algorithm from the Georgetown University CSL group, ported the code into Java, and created a graphical user interface (GingerALE). A cluster analysis script was added that identifies ALE clusters (areas of high activation likelihood) and returns the cluster extent above a user-specified threshold, x-y-z coordinates of the weighted center-of-mass and peak locations, and an anatomical label assigned by the Talairach Daemon (Lancaster et al., 2000
). A coordinate conversion utility was also included to convert MNI coordinates to Talairach space (Lancaster et al., 2007
). Two extensions of the original ALE method, a correction for multiple comparisons and a method for computing statistical contrasts of pairs of ALE images, were also added (Laird et al., 2005c
).
In the original implementation of the ALE method, several limitations were known to exist: (1) the size of the modeled Gaussian distribution was user-specified and therefore subjective, (2) the permutation test for significance was not anatomically constrained, leading to some modeled activation in white matter, and (3) the analysis tested for above-chance clustering of individual coordinates (a fixed-effects analysis), preventing the generalization of results that is possible in a random-effects analysis (Wager et al., 2007
, 2009
). Recent advances in the ALE technique have overcome these limitations to provide a more valid and statistically reliable meta-analysis framework (Eickhoff et al., 2009b
). Rather than relying on user-dependent Gaussian distributions, quantitative estimates of the between-subject and between-template variability were empirically determined in order to more explicitly model the spatial uncertainty associated with each coordinate (a correction that also includes a weighting of each study by the number of included subjects). In addition, the permutation test was limited to regions of gray matter and modified to test for the above-chance clustering between experiments, resulting in a transition from a fixed-effects to a random-effects method of statistical inference. By progressing from an analysis based on the clustering across coordinates to the clustering across experiments, ALE results no longer may potentially be driven by a single study. The new ALE formulation was validated against the classical algorithm and experimental data (Grefkes et al., 2008
) and found to increase the specificity of results without losing the sensitivity of the original approach. These improvements have been implemented in the most recent version of GingerALE, which is currently available for beta testing on the BrainMap website.
The ALE meta-analysis method can be applied in a variety of ways to answer specific research questions. Most frequently, ALE is applied to sets of neuroimaging studies that share some similar aspect of evoked brain function. These function-based meta-analyses usually involve pooling studies with similar experimental designs, and these studies may also be segregated into different collections to evaluate the functional specificity of the task or process being investigated. In the BrainMap coding scheme, experimental conditions are described according to the presented stimuli, instructions given, and response requested. Each of these fields has candidate entries to choose from, which collectively capture the essence of the scanned conditions (Figure 5
). These three condition axes (stimulus, response, and instructions) provide a structure in which differential activation patterns can be systematically probed for variations of a given task.
Figure 5. BrainMap conditions. In each paper archived in BrainMap, the scanned experimental conditions can be segregated into components based on the stimulus presented (modality and type), instructions given, and required response (modality and type). The values in the BrainMap metadata scheme are standardized to as few options as possible, so as to group similar studies. The blue background represents that this field lies at the Paper level.
For example, we performed an ALE meta-analysis of all studies in BrainMap that were coded with a paradigm class of “word generation” (66 papers, 197 experiments, 1552 coordinates), a widely used test of neuropsychological function. The meta-analytic results revealed extensive convergence in large portions of the left inferior frontal gyrus, centering on Brodmann area 44/45 (Broca’s area), and the left dorsolateral prefrontal cortex (DLPFC, BA 46/9), regions commonly known to be associated with word retrieval and executive function. ALE clusters were also observed in the bilateral insula, anterior cingulate cortex (ACC, BA 32), supplementary motor area (SMA, BA6), precuneus (BA 7), posterior cingulate cortex (PCC, BA 31), left posterior temporal cortex (Wernicke’s area, BA 22), left inferior parietal cortex (BA 40), right posterior cerebellum, and left thalamus. These areas are generally understood to be involved with the production of language and executive processing that is characteristic of verbal fluency tasks (Heim et al., 2008
; Petersen et al., 1988
; Warburton et al., 1996
; Wise et al., 1991
).
We then examined BrainMap metadata for these studies and found that the imaged tasks varied according to the modality of responses (covert or overt), the modality of stimulus presentation (visual or auditory), or the stimulus type (words or letters). Supplemental ALE meta-analyses were then performed according to each of these task variations, yielding differential patterns of activation likelihood (Figure 6
). Covert word generation yielded more extensive engagement of lateral and medial prefrontal areas (similar to Basho et al., 2007
), precuneus, and posterior cingulate cortex, while analysis of overt tasks revealed distinct concordance in the right cerebellum. ALE results of studies using visual stimuli were observed in visual cortex (BA 17/18) and the fusiform gyrus (BA 37), while auditory stimuli were localized to left auditory cortex (BA 41), Wernicke’s area (BA 22), and precuneus. Generation of words in response to a presented word (visual or auditory) was associated with the middle frontal gyrus (BA 9). Semantic verbal fluency was also associated with the ventral portion of the left inferior frontal gyrus, in agreement with previous meta-analysis results (Costafreda et al., 2006
), as well as BA 45 (Amunts et al., 2004
). The observation that Wernicke’s area was preferentially involved auditory stimulation, particularly with words, likely reflects this region’s role in auditory processing and the comprehension of spoken words. This word generation meta-analysis highlights the rich data mining that is possible using the BrainMap coding scheme.
Figure 6. Function-based meta-analysis of word generation. Studies utilizing a word generation paradigm were downloaded from the BrainMap database. Papers were segregated according to stimulus and response features. Separate ALE meta-analyses were performed for each set of coordinates: overt speech responses (30 papers, 78 experiments, 604 coordinates), covert speech responses (35 papers, 92 experiments, 729 coordinates), visual stimuli (39 papers, 115 experiments, 906 coordinates), auditory stimuli (32 papers, 86 experiments, 621 coordinates), visual words (25 papers, 88 experiments, 638 coordinates), visual letters (14 papers, 38 experiments, 373 coordinates), auditory words (39 papers, 124 experiments, 1001 coordinates), and auditory letters (14 papers, 35 experiments, 279 coordinates). Results are displayed at P < 0.05, corrected. Left lateral (x = −44) and medial (x = 0) sagittal slices display the comparative meta-analytic results. In all images, the overlap between meta-analysis results is shown in yellow.
While function-based meta-analyses tend to dominate the literature, structure-based meta-analyses offer an alternative meta-analytic strategy. Instead of pooling studies that share a common experimental design, structure-based meta-analyses focus on a specific anatomical region and look for global coactivation patterns across a diverse range of tasks. The theory behind this type of meta-analysis is that groups of coordinates that coactivate across experiments can be pooled to identify functionally connected networks in the brain. Like other methods of analyzing functional connectivity (Cordes et al., 2000
; Rogers et al., 2007
; Xiong et al., 1999
), structure-based meta-analyses are based on the co-occurrence of spatially separate neurophysiological events. Koski and Paus (2000)
used this technique to study the meta-analytic connectivity of the anterior cingulate cortex, although their analysis was limited to the frontal lobe. In their study, the authors manually collected and filtered data from 107 studies tasks to examine regional co-occurrences and found evidence for functional heterogeneity within the ACC. A similar meta-analysis was performed on basal ganglia activation using 126 published studies to determine the functional connectivity between cortex and striatum (Postuma and Dagher, 2006
).
Increasing both the size and diversity of structure-based meta-analyses adds to the generalizability of the results. When used in conjunction with the BrainMap database, the procedure is more automated, resulting in larger meta-analyses that include decades of neuroimaging data from a diverse range of paradigms and behavioral domains. Recently, a large-scale meta-analysis of the functional connectivity of the amygdala was carried out in which the ROIs for left and right amygdala were defined according the Harvard-Oxford structural probability atlas distributed with FSL (Smith et al., 2004
) and seeded in BrainMap (Robinson et al., 2009
). This anatomically defined meta-analysis of 240 papers (326 experiments with 3842 coordinates) revealed that the amygdala plays a integrative role in both emotion and cognition, and was validated according to the non-human primate database, CoCoMac (Stephan et al., 2001
).
To illustrate the use of structure-based meta-analyses, we note that the function-based meta-analysis of verbal fluency studies (Figure 6
) was characterized by extensive activation of the left inferior frontal gyrus, centered on Broca’s area in BA 44 and 45. To determine what connectivity differences exist between these two cytoarchitectonic regions, a location query was performed within BrainMap for studies activating these regions, using ROIs of the cytoarchitectonically defined areas (Amunts et al., 1999
) as distributed with the SPM Anatomy Toolbox (Eickhoff et al., 2005
, 2006b
). The returned studies were then analyzed in a structure-based meta-analysis of both left BA 44 and BA 45 (Figure 7
). For both regions, extensive bilateral connectivity was observed across the inferior frontal gyrus, precentral gyrus, inferior parietal lobule, fusiform gyrus, and insula, as well as medial ACC and SMA. Joint connectivity was also observed in the left thalamus and right cerebellum. In contrast to the word generation meta-analysis, no connectivity was observed in the precuneus or posterior cingulate cortex, likely reflecting a memory retrieval component (Cabeza and Nyberg, 2000
) of verbal fluency processing. Much overlap was observed between the two images; however, comparison of the maps for BA 44 and BA 45 revealed strong dissociation in subcortical regions. BA 44 exhibited extensive connectivity with bilateral thalamus, caudate, and putamen in agreement with Eickhoff et al. (2009a)
, while only the left thalamus was observed for BA 45. In addition, only the BA 44 map returned connectivity in Wernicke’s area, reflecting preferential engagement of this region in comparison to BA 45. Structure-based meta-analyses can therefore find distinct differences in functionally connected networks even for regions that lie very close to each other, such as BA 44 and 45.
Figure 7. Structure-based meta-analysis of Brodmann areas 44 and 45. ALE meta-analyses were performed for all experiments in BrainMap that reported activation in either BA 44 or BA 45 to determine the meta-analytic functional connectivity of these cytoarchitectonically-defined regions. Results are displayed at P < 0.05, corrected. The top slices (left and right) are centered at x = −54, y = 18, z = 24. The upper left slices display the ROIs used to search BrainMap for BA 44 (red) and BA 45 (blue), which were obtained from the SPM Anatomy Toolbox (Eickhoff et al., 2005
; Eickhoff et al., 2006b
). ALE results for these regions are shown on the top right in axial, sagittal, and coronal slices, as well as on the bottom panel in axial slices from z = 58 to z = −26. The overlap between meta-analyses for BA 44 and 45 is shown in purple.
Toro et al. (2008)
expanded upon the idea of structure-based meta-analyses and developed an algorithm to test the likelihood of a functional connection between regions, yielding a 3D meta-coactivation map for every voxel in the brain. In their analysis of 825 papers (3402 experiments with 27,909 coordinates), Toro et al. observed distinct and recognizable functional networks that are commonly associated with processes such as attention, motor function, and the resting state. Given that meaningful networks were extracted from the coordinates contained in BrainMap via a coactivation analysis, a recent study proposed that known functional networks of the brain during explicit activation could be derived using independent component analysis (ICA) of BrainMap data (Smith et al., 2009
). These networks, when compared to resting state networks (RSNs) obtained by ICA of resting state fMRI data (Damoiseaux et al., 2006
) were virtually identical. That is, the set of major covarying activation networks identified from a massive-scale meta-analysis (1687 papers, 7342 experiments, 58,620 coordinates) matched the set of networks that are present in the resting brain. These results provide strong evidence that RSNs reflect functional neural networks, and that these dynamic networks are engaged even at rest (Fox and Raichle, 2007
). Given the independent nature of these two analyses on fundamentally different types of data, as well as the heterogeneity of data contained in BrainMap due to differences in subjects, scanners, analyses, and paradigms, it is remarkable that such strong correspondence was observed between the resting state and meta-analytic results. In sum, this study supports the validity of using BrainMap and coordinate-based meta-analyses to identify functional neural networks on a large scale.
One of the broad goals of functional neuroimaging research is to determine function–structure relationships in the brain. A concrete deliverable of this aim is a probabilistic functional atlas, in which specific mental operations are mapped to discrete networks of brain regions. Price and Friston (2005)
point out that the relationship between a brain region and a mental function is not a one-to-one mapping. Instead this relationship is a many-to-many mapping, as a single region can be involved in many cognitive processes, and a single process usually activates multiple regions. Evaluating these mappings will require collating the immense amount of neuroimaging data that has been acquired thus far and continuing the development of advanced meta-analytic techniques in order to efficiently and effectively synthesize all of these data.
As the development of comprehensive neuroinformatics tools progresses, the need for comprehensive data ontologies increases. An ontology is a machine-interpretable description of concepts and their relationships with the purpose of sharing of ideas and information in a manner facilitated by semantic interoperability (Stevens et al., 2000
). Until a foundational ontology for neuroimaging is established and adopted, the communication within and between databases will be limited, hindering the creation of a functional brain atlas. In the field of neuroimaging, ontology development is proceeding rapidly in the domains of representing neuroanatomical findings [e.g., NeuroNames, Bowden and Dubach, 2003
; Bowden and Martin, 1995 and the Foundational Model of Anatomy (FMA), Rosse and Mejino, 2003
, describing imaging acquisition strategies (e.g., RadLex, Langlotz, 2006
; Rubin, 2008
), and identifying clinical assessments [e.g., the Systematized Nomenclature of Medicine (SNOMED), Coté and Robboy, 1980
. In addition, the Neuroscience Information Framework Standardized (NIFSTD) Ontology, developed by BIRN and the NIF, is a collection of these and other neuroscience ontologies (Bug et al., 2008
).
There is currently no accepted ontology for describing the range of mental operations performed by the human brain, although need for such an ontology is increasingly being discussed (Binder et al., 2009
; Poldrack, 2006
; Price and Friston, 2005
; Toga, 2002
). A research question such as “find the data examining the relationship between hippocampal volume in Alzheimer’s disease” requires knowledge about clinical diagnosis, neuroanatomy, and MR scan acquisition. As mentioned previously, there are ongoing ontological efforts designed to curate knowledge in these domains. However, a query such as “find the data examining neural activations observed during sustained attention” involves knowledge about cognitive processing, for which no ontology exists to date. Creating a realistic functional brain atlas will require a systematic description of mental operations that are reported in the literature so that informative and standardized labels can be applied to different brain networks.
Investigators frequently utilize alternate and sometimes competitive terminologies when referring to the cognitive processes elicited by specific tasks. When different words are employed to represent the same concepts, the grouping of related ideas across different resources is impeded. For example, a semantically interoperable ontology will allow the linking of data designated with terms such as “declarative” and “explicit” memory types and be capable of relating “working memory” to either “memory” or “executive processing”. Cognition represents the most difficult domain to explicate as it contains the most intricate of all concepts, such as language, attention, and memory. But while much imaging research focuses on cognitive processing, important results are also being published in the areas of perception, action, emotion, and autonomic functions. To comprehensively describe the many-to-many mappings of structure and function being investigated in neuroimaging research, a complete mental ontology must be developed.
Poldrack (2008)
suggests that the many-to-many mapping dilemma may be complicated by our current understanding of what cognitive processes exist and how they are defined across functional neuroimaging experiments. Optimally, an appropriate and useful cognitive ontology will not merely be a catalogue of various mental operations parsed down to very fine detail in accord with current theories of cognitive psychology. While the consideration of competing theories often results in new knowledge discovery, the development of such a top-down ontology would be so continuously and vehemently debated that it could never reach a sufficient degree of consensus in order to be considered adoptable by the neuroimaging community. In contrast, a biologically based ontology that is driven by the way we observe the brain to operate in imaging experiments may reveal a cognitive architecture that has not previously been considered.
The synergy that exists between the BrainMap database and the ALE meta-analysis method was designed to facilitate the creation of a functional brain atlas. BrainMap’s search capabilities can support different types of queries, such as “for a given function, what regions are typically engaged?”, “for a given region, what tasks elicit activation?”, or “for a given region, what other regions are coactivated?” These questions highlight the value of meta-analytic results in comparison to results from individual studies. We believe that these correspondences (function-to-regions, region-to-tasks, or region-to-network) must be constructed according to a bottom-up strategy, using knowledge gleaned from data-driven analyses. Current probabilistic structural atlases (e.g., the Harvard-Oxford structural probability atlas, Smith et al., 2004
, or the Jülich cytoarchitectonic atlas, Eickhoff et al., 2005
, 2006a
) have proven to be useful in testing region-to-tasks or region-to-network associations, as shown in Figure 7
. Theoretically, any method of determining regions of interest can be used to query the BrainMap database, whether structurally or functionally defined.
BrainMap is also capable of generating function-to-regions associations, albeit at a coarse resolution. Whole brain meta-analytic maps can be created for each behavioral domain category in BrainMap, which can then be decomposed into sub-networks based on different levels of the domain hierarchy. To illustrate, ALE meta-analyses were performed on nine different behavioral domain categories and sub-categories: action, cognition, emotion, perception, interoception, language, memory, vision, and audition (Figure 8
). Each ALE image provides a unique mapping of the neural network associated with the relevant domain. Many regions are observed in multiple domain maps, and some maps are very similar, but none are identical. However, a more detailed domain structure is needed to fully characterize the range of human cognition. We propose that applying high-level filters from the entire BrainMap coding scheme to these meta-maps, in the way condition-based filters were applied in the word generation meta-analysis (Figure 6
), can be an effective strategy for refining the spatial specificity of these images. Thus, while paradigm class and behavioral domain are important metadata fields in the BrainMap coding scheme, all fields have the potential to assist in unraveling the brain’s systems and their interactions.
Figure 8. BrainMap’s functional atlas development strategy. ALE meta-analyses were performed on nine arbitrarily selected behavioral domain categories and sub-categories: action, cognition, emotion, perception, interoception, language, memory, vision, and audition. These results illustrate the mapping of function-to-regions correspondences that are currently capable using the BrainMap database. Refinement of these maps can be accomplished through high-level filtering and mining of the database. Results are displayed at P < 0.05, corrected. The total number of coordinates in each meta-analysis is listed after the domain heading for each map in parentheses.
As the BrainMap database increases in size, these results will evolve and grow more powerful, perhaps leading to a multi-layered, multi-modal probabilistic functional brain atlas derived from many different large-scale coordinate-based meta-analyses. This approach would likely be enhanced by the development of a standardized mental ontology. Differences in competitive terminology must be resolved to allow for the union of experimentally similar data sets. Perhaps the best strategy would be to combine all of BrainMap’s data-driven methods in establishing function–structure relationships with other ontology initiatives, such at the Cognitive Atlas
7
(Bilder et al., 2009
), the NIFSTD ontology (Bug et al., 2008
), or the Neural ElectroMagnetic Ontologies (NEMO) (Frishkoff et al., 2009
). Given the complex nature of human brain function, it is reasonable to suggest that no single approach will be powerful enough to solve the fundamental challenges associated with mapping the mind, but rather a joint effort will be required.
The authors report no competing interests.
This work was supported by the NIMH (R01-MH074457-03; PI = Peter Fox and R01-MH084812-01; PI = Angela Laird and Jessica Turner), NINDS (T35-NS051166-04; PI = Peter Fox), NIDCD (F32-DC009116-02; PI = Matthew Cykowski), NCRR (U24-RR021992; PI = Steven Potkin), and the Helmholz Initiative on Systems-Biology (SBE).
Amunts, K., Weiss, P. H., Mohlberg, H., Pieperhoff, P., Eickhoff, S., Gurd, J. M., Marshall, J. C., Shah, N. J., Fink, G. R., and Zilles, K. (2004). Analysis of neural mechanisms underlying verbal fluency in cytoarchitectonically defined stereotactic space – the roles of Brodmann areas 44 and 45. Neuroimage 22, 42–56.
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., Bannister, P. R., De Luca, M., Drobnjak, I., Flitney, D. E., Niazy, R. K., Saunders, J., Vickers, J., Zhang, Y., De Stefano, N., Brady, J. M., and Matthews, P. M. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23(Suppl. 1), S208–S219.