Elsevier

NeuroImage

Volume 28, Issue 3, 15 November 2005, Pages 635-662
NeuroImage

A Population-Average, Landmark- and Surface-based (PALS) atlas of human cerebral cortex

https://doi.org/10.1016/j.neuroimage.2005.06.058Get rights and content

Abstract

This report describes a new electronic atlas of human cerebral cortex that provides a substrate for a wide variety of brain-mapping analyses. The Population-Average, Landmark- and Surface-based (PALS) atlas approach involves surface-based and volume-based representations of cortical shape, each available as population averages and as individual subject data. The specific PALS-B12 atlas introduced here is derived from structural MRI volumes of 12 normal young adults. Accurate cortical surface reconstructions were generated for each hemisphere, and the surfaces were inflated, flattened, and mapped to standard spherical configurations using SureFit and Caret software. A target atlas sphere was generated by averaging selected landmark contours from each of the 24 contributing hemispheres. Each individual hemisphere was deformed to this target using landmark-constrained surface registration. The utility of the resultant PALS-B12 atlas was demonstrated using a variety of analyses. (i) Probabilistic maps of sulcal identity were generated using both surface-based registration (SBR) and conventional volume-based registration (VBR). The SBR approach achieved markedly better consistency of sulcal alignment than did VBR. (ii) A method is introduced for ‘multi-fiducial mapping’ of volume-averaged group data (e.g., fMRI data, probabilistic architectonic maps) onto each individual hemisphere in the atlas, followed by spatial averaging across the individual maps. This yielded a population-average surface representation that circumvents the biases inherent in choosing any single hemisphere as a target. (iii) Surface-based and volume-based morphometry applied to maps of sulcal depth and sulcal identity demonstrated prominent left–right asymmetries in and near the superior temporal sulcus and Sylvian fissure. Moreover, shape variability in the temporal lobe is significantly greater in the left than the right hemisphere. The PALS-B12 atlas has been registered to other surface-based atlases to facilitate interchange of data and comparison across atlases. All data sets in the PALS-B12 atlas are accessible via the SumsDB database for online and offline visualization and analysis.

Introduction

Atlases play an increasingly important role in providing a spatial and structural framework for visualizing and analyzing many aspects of brain structure, function, and development (Van Essen, 2002, Toga and Thompson, 2002, Mazziotta et al., 2001). In recent decades, numerous human brain atlases have been developed, each with its particular advantages and attractive features but also with various limitations and drawbacks.

For human cerebral cortex, the major challenges in generating a suitable atlas arise from the fact that cortical convolutions are very complex and are highly variable from one individual to the next. To address these challenges, strategic choices must be made on four major fronts: (i) representing cortical morphology and its variability; (ii) visualizing experimental data and underlying cortical structure; (iii) selecting coordinate systems; and (iv) choosing algorithms for registering data from individuals to the atlas. On each of these fronts, a key issue is whether the approach is based on surfaces, volumes, or a combination of both. (i) Representing morphology and individual variability. The core of a brain atlas is a representation of brain morphology derived from one or more individuals, captured by in vivo imaging or postmortem histological methods. Structural magnetic resonance imaging (sMRI) is particularly attractive as a substrate for a human cortical atlas because it can be readily obtained in a population of normal subjects, and it allows for isotropic volume representations that reveal the full complexity of cortical folds. Some human cortical atlases are based on an individual brain or hemisphere, which allows accurate representation of cortical shape in the selected individual but fails to represent variability across the population (Talairach and Tournoux, 1988, Roland et al., 1994, Van Essen and Drury, 1997, Geyer et al., 2001). Other atlases are based on averaging structural MRI volumes across a population of subjects, which capture some shape consistencies across the population, but invariably lead to extensive blurring of cortical convolutions (Evans et al., 1994, Friston et al., 1995, Mazziotta et al., 2001). (ii) Visualization and analysis options. Cortical structure is sufficiently complex that any single visualization mode is inadequate when used in isolation. It is particularly valuable to be able to view atlas data using a combination of volume and surface visualization. (iii) Coordinate systems. Spatial coordinates provide a concise and objective way to represent locations in human cortex. The most basic distinction is between 3D stereotaxic (x,y,z) coordinates and surface-based coordinates (e.g., latitude and longitude on a spherical map). Within each type, numerous specific coordinate systems are in common use, and these differ from one another in important respects. (iv) Registration algorithms. The objective of registering individual brains to an atlas is to compensate for multiple types of variability, including the initial position and orientation of the individual brain, the overall dimensions of each hemisphere, the pattern of cortical convolutions, and the relationship of functional subdivisions (cortical areas) to these convolutions (Toga, 1999). Given the nature of cortical structure and function, it is not possible to remove all variability, i.e., to place all geographic and functional subdivisions into precise correspondence throughout an atlas. Volume-based registration (VBR) algorithms are fundamentally limited because they fail to respect the topology of the cortical sheet and hence can lead to topologically incorrect mappings. Surface-based registration (SBR) avoids this limitation and allows more flexible application of relevant geographic and/or functional constraints, using explicit landmarks (Van Essen et al., 2005) or using metrics that reflect shape characteristics throughout the cortical sheet (Fischl et al., 1999b). However, there are as yet few evaluations of different SBR methods using objective performance measures (Desai et al., 2005).

The present report introduces the Population-Average, Landmark- and Surface-based (PALS) approach and illustrates attractive solutions it offers for each of the aforementioned challenges in generating human cerebral cortical atlases. The specific atlas illustrated here is based on structural MRI volumes from 12 normal subjects and is accordingly identified as the PALS-B12 atlas to distinguish it from other PALS atlases derived from a different set of brains. The PALS-B12 atlas data set includes (i) structural MRI volumes from individual brains and the population average; (ii) cortical segmentations and fiducial surface reconstructions for all contributing hemispheres; (iii) inflated, flattened, and spherical configurations for each hemisphere; (iv) all hemispheres registered to a target atlas surface whose landmarks reflect shape characteristics of the population rather than any individual; (v) individual surfaces resampled to a standard mesh to facilitate simultaneous visualization; (vi) multiple coordinate systems, both stereotaxic and surface-based, for expressing cortical location; (vii) maps of sulcal depth and other shape characteristics for individuals and group averages; (viii) maps of sulcal identity for individuals and the population average; and (ix) many additional types of experimental data mapped from other hemispheres to the PALS-B12 atlas.

Central to the PALS approach is the choice of geographic landmarks used to constrain the registration from individual hemispheres to the atlas. The PALS atlas uses only six well-defined and consistent geographic landmark contours, because other gyral and sulcal features in human cortex are so variable as to be unsuitable (or at least problematic) for use as landmarks. The effectiveness of this approach is demonstrated here using a variety of quantitative analyses related to sulcal shape and sulcal identity. One set of analyses relates to the issue of asymmetries between the left and right hemispheres. Although there is an extensive literature on left–right morphological asymmetries (Toga and Thompson, 2003, Ochiai et al., 2004; see Discussion), much remains unresolved about the nature and magnitude of these differences. The surface-based morphometry (SBM) methods developed in the present study provide valuable additional information regarding left–right asymmetries, especially in the temporal lobe. They are complementary to other SBM analyses (MacDonald et al., 2000, Chung et al., 2003, Salat et al., 2004) that have focused on comparing maps of cortical thickness in different groups rather than the pattern of cortical convolutions.

Surface-based atlases are especially useful for analyzing and comparing complex functional activation patterns obtained in many thousands of fMRI studies involving a wide variety of behavioral tasks and analysis paradigms. Currently, the great majority of fMRI studies obtain statistical power by averaging across subjects using volume-based registration. Mapping such volume-based population-average data onto an individual atlas surface facilitates visualization and comparisons across studies (Van Essen, 2002), but it inevitably introduces significant biases related to the particular convolutions of the target hemisphere. Here, a method for circumventing this bias is introduced, in which fMRI data are mapped separately onto each individual fiducial surface in the PALS-B12 atlas population, then averaged across all maps. This yields an objective representation of complex activation patterns and avoids the bias of a particular individual atlas target.

All of the data sets used in this study are accessible from the SumsDB database (http://sumsdb.wustl.edu:8081/sums). These data can be readily visualized online (using WebCaret) and can be downloaded for offline visualization and analysis (using Caret software). A ‘scene display’ option allows immediate viewing (in either Caret or WebCaret) of the precise screen displays used in generating each of the Caret-derived panels of each figure in this paper.

Section snippets

Methods

Structural MRI volumes were obtained from 24 normal human subjects, all right-handed young adults (ages 18–24; data from Buckner et al., 2004, Head et al., in press). Data from six males and six females were used for the PALS-B12 target atlas. Multiple (three or four) T1-weighted MP-RAGE scans (1.5 T Siemens Vision scanner, 1 × 1 × 1.25 mm voxels) were acquired from each individual. Volumes were registered to the Washington University 711-2C version of stereotaxic atlas space. The 711-2C target

Results

The results are presented in eight sections that address three broad topics. One topic concerns the diversity of cortical shape characteristics; the results can be viewed in a population of individuals (Diversity of sulcal patterns section) or using population-average representations of sulcal depth (Population-based shape representations section) and sulcal identity (Probabilistic surface and volume maps of sulcal identity section). A second topic focuses on methodological considerations,

Discussion

This study was carried out with six major objectives in mind: (i) to introduce the PALS atlas approach as a general way to cope with the complexity of human cortical convolutions and their variability across individuals; (ii) to make the PALS-B12 atlas data set available, including a wide variety of surface and volume data, along with a flexible set of online and offline visualization tools; (iii) to evaluate and compare the quality of surface-based and volume-based registration; (iv) to

Acknowledgments

I thank J. Harwell, D. Hanlon, and J. Dickson for superlative software development; E. Reid for excellent data analysis; S. Danker for manuscript preparation; R. Buckner for MRI data sets; J. Hegde for advice on statistical methods; K. Zilles, K. Amunts, and H. Mohlberg for providing probabilistic architectonic data sets; and A. Snyder, M. Raichle, R. Buckner, and D. Fair for valuable discussions. This project was supported by NIH Research Grant R01-MH-60974, funded by the National Institute of

References (65)

  • D. MacDonald et al.

    Automated 3-D extraction of inner and outer surfaces of cerebral cortex from MRI

    NeuroImage

    (2000)
  • P. Morosan et al.

    Human primary auditory cortex: cytoarchitectonic subdivisions and mapping into a spatial reference system

    Neuroimage

    (2001)
  • T. Ochiai et al.

    Sulcal pattern and morphology of the superior temporal sulcus

    Neuroimage

    (2004)
  • J.G. Ojemann et al.

    Anatomic localization and quantitative analysis of gradient refocused echo-lanaar fMRI suceptibility artifacts

    Neuroimage

    (1997)
  • J. Rademacher et al.

    Probabilistic mapping and volume measurement of human primary auditory cortex

    Neuroimage

    (2001)
  • D. Riviere et al.

    Automatic recognition of cortical sulci of the human brain using a congregation of neural networks

    Med. Image Anal.

    (2002)
  • A.W. Toga et al.

    New approaches in brain morphometry

    Am. J. Geriatr. Psychiatry

    (2002)
  • D.C. Van Essen

    Windows on the brain. The emerging role of atlases and databases in neuroscience

    Curr. Opin. Neurobiol.

    (2002)
  • K. Zilles et al.

    Quantitative analysis of cyto- and receptor architecture of the human brain

  • K. Amunts et al.

    Advances in cytoarchitectonic mapping of the human cerebral cortex

  • K. Amunts et al.

    Broca's region revisited: cytoarchitecture and intersubject variability

    J. Comp. Neurol.

    (1999)
  • K. Brodmann

    Vergleichende Lokalisationslehre der Grosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues

    (1909)
  • R.L. Buckner et al.

    Dissociation of human prefrontal cortical areas across different speech production tasks and gender groups

    J. Neurophysiol.

    (1995)
  • M. Corbetta et al.

    Control of goal-directed and stimulus-driven attention in the brain

    Nat. Rev. Neurosci.

    (2002)
  • M. Corbetta et al.

    Voluntary orienting is dissociated from target detection in human posterior parietal cortex

    Nat. Neurosci.

    (2000)
  • R.F. Dougherty et al.

    Visual field representations and locations of visual areas V1/2/3 in human visual cortex

    J. Vis.

    (2003)
  • H.A. Drury et al.

    Computerized mappings of the cerebral cortex. A multiresolution flattening method and a surface-based coordinate system

    J. Cogn. Neurosci.

    (1996)
  • H.A. Drury et al.

    Surface-based analyses of the human cerebral cortex

  • A.C. Evans et al.

    Three-dimensional correlative imaging: applications in human brain mapping

  • D.J. Felleman et al.

    Distributed hierarchical processing in primate cerebral cortex

    Cereb. Cortex

    (1991)
  • B. Fischl et al.

    High-resolution intersubject averaging and a coordinate system for the cortical surface

    Hum. Brain Mapp.

    (1999)
  • K.J. Friston et al.

    Statistical parametric maps in functional imaging: a general linear approach

    Hum. Brain Mapp.

    (1995)
  • Cited by (985)

    View all citing articles on Scopus
    View full text