Neuroimage

In this study, we have assessed the validity and reliability of an automated labeling system that we have developed for subdividing the human cerebral cortex on magnetic resonance images into gyral based regions of interest (ROIs). Using a dataset of 40 MRI scans we manually identified 34 cortical ROIs in each of the individual hemispheres. This information was then encoded in the form of an atlas that was utilized to automatically label ROIs. To examine the validity, as well as the intra-and inter-rater reliability of the automated system, we used both intraclass correlation coefficients (ICC), and a new method known as mean distance maps, to assess the degree of mismatch between the manual and the automated sets of ROIs. When compared with the manual ROIs, the automated ROIs were highly accurate, with an average ICC of 0.835 across all of the ROIs, and a mean distance error of less than 1 mm. Intra-and inter-rater comparisons yielded little to no difference between the sets of ROIs. These findings suggest that the automated method we have developed for subdividing the human cerebral cortex into standard gyral-based neuroanatomical regions is both anatomically valid and reliable. This method may be useful for both morphometric and functional studies of the cerebral cortex as well as for clinical investigations aimed at tracking the evolution of disease-induced changes over time, including clinical trials in which MRI-based measures are used to examine response to treatment. D

Active contour segmentation and its robust implementation using level set methods are well-established theoretical approaches that have been studied thoroughly in the image analysis literature. Despite the existence of these powerful segmentation methods, the needs of clinical research continue to be fulfilled, to a large extent, using slice-by-slice manual tracing. To bridge the gap between methodological advances and clinical routine, we developed an open source application called ITK-SNAP, which is intended to make level set segmentation easily accessible to a wide range of users, including those with little or no mathematical expertise. This paper describes the methods and software engineering philosophy behind this new tool and provides the results of validation experiments performed in the context of an ongoing child autism neuroimaging study. The validation establishes SNAP intrarater and interrater reliability and overlap error statistics for the caudate nucleus and finds that SNAP is a highly reliable and efficient alternative to manual tracing. Analogous results for lateral ventricle segmentation are provided. D

Voxel-based-morphometry (VBM) is a whole-brain, unbiased technique for characterizing regional cerebral volume and tissue concentration differences in structural magnetic resonance images. We describe an optimized method of VBM to examine the effects of age on grey and white matter and CSF in 465 normal adults. Global grey matter volume decreased linearly with age, with a significantly steeper decline in males. Local areas of accelerated loss were observed bilaterally in the insula, superior parietal gyri, central sulci, and cingulate sulci. Areas exhibiting little or no age effect (relative preservation) were noted in the amygdala, hippocampi, and entorhinal cortex. Global white matter did not decline with age, but local areas of relative accelerated loss and preservation were seen. There was no interaction of age with sex for regionally specific effects. These results corroborate previous reports and indicate that VBM is a useful technique for studying structural brain correlates of ageing through life in humans.

Neuroimaging studies with positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) have begun to describe the functional neuroanatomy of emotion. Taken separately, specific studies vary in task dimensions and in type(s) of emotion studied and are limited by statistical power and sensitivity. By examining findings across studies, we sought to determine if common or segregated patterns of activations exist across various emotional tasks. We reviewed 55 PET and fMRI activation studies (yielding 761 individual peaks) which investigated emotion in healthy subjects. Peak activation coordinates were transformed into a standard space and plotted onto canonical 3-D brain renderings. We divided the brain into 20 nonoverlapping regions, and characterized each region by its responsiveness across individual emotions (positive, negative, happiness, fear, anger, sadness, disgust), to different induction methods (visual, auditory, recall/imagery), and in emotional tasks with and without cognitive demand. Our review yielded the following summary observations: (1) The medial prefrontal cortex had a general role in emotional processing; (2) fear specifically engaged the amygdala; (3) sadness was associated with activity in the subcallosal cingulate; (4) emotional induction by visual stimuli activated the occipital cortex and the amygdala; (5) induction by emotional recall/imagery recruited the anterior cingulate and insula; (6) emotional tasks with cognitive demand also involved the anterior cingulate and insula. This review provides a critical comparison of findings across individual studies and suggests that separate brain regions are involved in different aspects of emotion. © 2002 Elsevier Science (USA)

A component based method (CompCor) for the reduction of noise in both blood oxygenation level-dependent (BOLD) and perfusionbased functional magnetic resonance imaging (fMRI) data is presented. In the proposed method, significant principal components are derived from noise regions-of-interest (ROI) in which the time series data are unlikely to be modulated by neural activity. These components are then included as nuisance parameters within general linear models for BOLD and perfusion-based fMRI time series data. Two approaches for the determination of the noise ROI are considered. The first method uses high-resolution anatomical data to define a region of interest composed primarily of white matter and cerebrospinal fluid, while the second method defines a region based upon the temporal standard deviation of the time series data. With the application of CompCor, the temporal standard deviation of resting-state perfusion and BOLD data in gray matter regions was significantly reduced as compared to either no correction or the application of a previously described retrospective image based correction scheme (RETROICOR). For both functional perfusion and BOLD data, the application of CompCor significantly increased the number of activated voxels as compared to no correction. In addition, for functional BOLD data, there were significantly more activated voxels detected with CompCor as compared to RETRO-ICOR. In comparison to RETROICOR, CompCor has the advantage of not requiring external monitoring of physiological fluctuations.

The question of the self has intrigued philosophers and psychologists for a long time. More recently, distinct concepts of self have also been suggested in neuroscience. However, the exact relationship between these concepts and neural processing across different brain regions remains unclear. This article reviews neuroimaging studies comparing neural correlates during processing of stimuli related to the self with those of non-self-referential stimuli. All studies revealed activation in the medial regions of our brains' cortex during self-related stimuli. The activation in these so-called cortical midline structures (CMS) occurred across all functional domains (e.g., verbal, spatial, emotional, and facial). Cluster and factor analyses indicate functional specialization into ventral, dorsal, and posterior CMS remaining independent of domains. Taken together, our results suggest that self-referential processing is mediated by cortical midline structures. Since the CMS are densely and reciprocally connected to subcortical midline regions, we advocate an integrated cortical -subcortical midline system underlying human self. We conclude that self-referential processing in CMS constitutes the core of our self and is critical for elaborating experiential feelings of self, uniting several distinct concepts evident in current neuroscience. D

We present an approach to characterizing the differences among event-related hemodynamic responses in functional magnetic resonance imaging that are evoked by different sorts of stimuli. This approach is predicated on a linear convolution model and standard inferential statistics as employed by statistical parametric mapping. In particular we model evoked responses, and their differences, in terms of basis functions of the peri-stimulus time. This facilitates a characterization of the temporal response profiles that has a high effective temporal resolution relative to the repetition time. To demonstrate the technique we examined differential responses to visually presented words that had been seen prior to scanning or that were novel. The form of these differences involved both the magnitude and the latency of the response components. In this paper we focus on bilateral ventrolateral prefrontal responses that show deactivations for previously seen words and activations for novel words. 1998 Academic Press

Typically in neuroimaging we are looking to extract some pertinent information from imperfect, noisy images of the brain. This might be the inference of percent changes in blood flow in perfusion FMRI data, segmentation of subcortical structures from structural MRI, or inference of the probability of an anatomical connection between an area of cortex and a subthalamic nucleus using diffusion MRI. In this article we will describe how Bayesian techniques have made a significant impact in tackling problems such as these, particularly in regards to the analysis tools in the FMRIB Software Library (FSL). We shall see how Bayes provides a framework within which we can attempt to infer on models of neuroimaging data, while allowing us to incorporate our prior belief about the brain and the neuroimaging equipment in the form of biophysically informed or regularising priors. It allows us to extract probabilistic information from the data, and to probabilistically combine information from multiple modalities. Bayes can also be used to not only compare and select between models of different complexity, but also to infer on data using committees of models. Finally, we mention some analysis scenarios where Bayesian methods are impractical, and briefly discuss some practical approaches that we have taken in these cases.

Paradigms drawn from cognitive psychology have provided new insight into covert stages of action. These states include not only intending actions that will eventually be executed, but also imagining actions, recognizing tools, learning by observation, or even understanding the behavior of other people. Studies using techniques for mapping brain activity, probing cortical excitability, or measuring the activity of peripheral effectors in normal human subjects and in patients all provide evidence of a subliminal activation of the motor system during these cognitive states. The hypothesis that the motor system is part of a simulation network that is activated under a variety of conditions in relation to action, either self-intended or observed from other individuals, will be developed. The function of this process of simulation would be not only to shape the motor system in anticipation to execution, but also to provide the self with information on the feasibility and the meaning of potential actions.

We present a set of techniques for embedding the physics of the imaging process that generates a class of magnetic resonance images (MRIs) into a segmentation or registration algorithm. This results in substantial invariance to acquisition parameters, as the effect of these parameters on the contrast properties of various brain structures is explicitly modeled in the segmentation. In addition, the integration of image acquisition with tissue classification allows the derivation of sequences that are optimal for segmentation purposes. Another benefit of these procedures is the generation of probabilistic models of the intrinsic tissue parameters that cause MR contrast (e.g., T1, proton density, T2*), allowing access to these physiologically relevant parameters that may change with disease or demographic, resulting in nonmorphometric alterations in MR images that are otherwise difficult to detect. Finally, we also present a high band width multiecho FLASH pulse sequence that results in high signal-to-noise ratio with minimal image distortion due to B0 effects. This sequence has the added benefit of allowing the explicit estimation of T2* and of reducing test-retest intensity variability. D

Whereas some older adults show significant cognitive deficits, others perform as well as young adults. We investigated the neural basis of these different aging patterns using positron emission tomography (PET). In PET and functional MRI (fMRI) studies, prefrontal cortex (PFC) activity tends to be less asymmetric in older than in younger adults (Hemispheric Asymmetry Reduction in Old Adults or HAROLD). This change may help counteract age-related neurocognitive decline (compensation hypothesis) or it may reflect an age-related difficulty in recruiting specialized neural mechanisms (dedifferentiation hypothesis). To compare these two hypotheses, we measured PFC activity in younger adults, low-performing older adults, and high-performing older adults during recall and source memory of recently studied words. Compared to recall, source memory was associated with right PFC activations in younger adults. Low-performing older adults recruited similar right PFC regions as young adults, but high-performing older adults engaged PFC regions bilaterally. Thus, consistent with the compensation hypothesis and inconsistent with the dedifferentiation hypothesis, a hemispheric asymmetry reduction was found in high-performing but not in low-performing older adults. The results suggest that low-performing older adults recruited a similar network as young adults but used it inefficiently, whereas high-performing older adults counteracted age-related neural decline through a plastic reorganization of neurocognitive networks. © 2002 Elsevier

Precise localization of sulco-gyral structures of the human cerebral cortex is important for the interpretation of morpho-functional data, but requires anatomical expertise and is time consuming because of the brain's geometric complexity. Software developed to automatically identify sulco-gyral structures has improved substantially as a result of techniques providing topologically correct reconstructions permitting inflated views of the human brain. Here we describe a complete parcellation of the cortical surface using standard internationally accepted nomenclature and criteria. This parcellation is available in the FreeSurfer package. First, a computerassisted hand parcellation classified each vertex as sulcal or gyral, and these were then subparcellated into 74 labels per hemisphere. Twelve datasets were used to develop rules and algorithms (reported here) that produced labels consistent with anatomical rules as well as automated computational parcellation. The final parcellation was used to build an atlas for automatically labeling the whole cerebral cortex. This atlas was used to label an additional 12 datasets, which were found to have good concordance with manual labels. This paper presents a precisely defined method for automatically labeling the cortical surface in standard terminology.

Diffusion-weighted (DW) MR images contain information about the orientation of brain white matter fibres that potentially can be used to study human brain connectivity in vivo using tractography techniques. Currently, the diffusion tensor model is widely used to extract fibre directions from DW-MRI data, but fails in regions containing multiple fibre orientations. The spherical deconvolution technique has recently been proposed to address this limitation. It provides an estimate of the fibre orientation distribution (FOD) by assuming the DW signal measured from any fibre bundle is adequately described by a single response function. However, the deconvolution is ill-conditioned and susceptible to noise contamination. This tends to introduce artefactual negative regions in the FOD, which are clearly physically impossible. In this study, the introduction of a constraint on such negative regions is proposed to improve the conditioning of the spherical deconvolution. This approach is shown to provide FOD estimates that are robust to noise whilst preserving angular resolution. The approach also permits the use of super-resolution, whereby more FOD parameters are estimated than were actually measured, improving the angular resolution of the results. The method provides much better defined fibre orientation estimates, and allows orientations to be resolved that are separated by smaller angles than previously possible. This should allow tractography algorithms to be designed that are able to track reliably through crossing fibre regions.

There is great interest in estimating brain "networks" from FMRI data. This is often attempted by identifying a set of functional "nodes" (e.g., spatial ROIs or ICA maps) and then conducting a connectivity analysis between the nodes, based on the FMRI timeseries associated with the nodes. Analysis methods range from very simple measures that consider just two nodes at a time (e.g., correlation between two nodes' timeseries) to sophisticated approaches that consider all nodes simultaneously and estimate one global network model (e.g., Bayes net models). Many different methods are being used in the literature, but almost none has been carefully validated or compared for use on FMRI timeseries data. In this work we generate rich, realistic simulated FMRI data for a wide range of underlying networks, experimental protocols and problematic confounds in the data, in order to compare different connectivity estimation approaches. Our results show that in general correlation-based approaches can be quite successful, methods based on higher-order statistics are less sensitive, and lag-based approaches perform very poorly. More specifically: there are several methods that can give high sensitivity to network connection detection on good quality FMRI data, in particular, partial correlation, regularised inverse covariance estimation and several Bayes net methods; however, accurate estimation of connection directionality is more difficult to achieve, though Patel's τ can be reasonably successful. With respect to the various confounds added to the data, the most striking result was that the use of functionally inaccurate ROIs (when defining the network nodes and extracting their associated timeseries) is extremely damaging to network estimation; hence, results derived from inappropriate ROI definition (such as via structural atlases) should be regarded with great caution. © 2010 Elsevier Inc. All rights reserved. ICA dimensionality). If a low number of components is estimated , then it makes more sense to think of each component itself as a network. This will often include several noncontiguous regions, all having the same timecourse according to the ICA model, and hence within-component network analysis is not possible without further processing, such as splitting the components and reestimating each resulting node's timeseries. Furthermore, betweencomponent network analysis is quite possibly not reasonable, as each component will in itself constitute a gross, complex functional system. However, if a higher number of components is estimated , these are more likely to be smaller, isolated regions (functional parcels), which can more sensibly be then considered as nodes for use in network analysis.

225 SA-AM What you did to others you did to me? Brain activation during altruistic punishment in a first person versus a third party dictator game, P. Kirsch, 239 SA-AM Self-control in decision-making involves modulation of the vmPFC valuation system, TA Hare, CF Camerer, A Rangel, Caltech, Pasadena, CA, USA 241 SA-AM Beliefs, guilt, and the brain: Neural mechanisms underlying social cooperative behavior, L.LANGUAGE Language Acquisition 245 SA-AM Global effects of three remediation programs on children with developmental dyslexia : an fMRI study, R. Nenert, C. Leveque, M.T. LeNormand, P. Evrard, S. De Schonen, Inserm U825,

Longitudinal image analysis has become increasingly important in clinical studies of normal aging and neurodegenerative disorders. Furthermore, there is a growing appreciation of the potential utility of longitudinally acquired structural images and reliable image processing to evaluate disease modifying therapies. Challenges have been related to the variability that is inherent in the available cross-sectional processing tools, to the introduction of bias in longitudinal processing and to potential over-regularization. In this paper we introduce a novel longitudinal image processing framework, based on unbiased, robust, within-subject template creation, for automatic surface reconstruction and segmentation of brain MRI of arbitrarily many time points. We demonstrate that it is essential to treat all input images exactly the same as removing only interpolation asymmetries is not sufficient to remove processing bias. We successfully reduce variability and avoid over-regularization by initializing the processing in each time point with common information from the subject template. The presented results show a significant increase in precision and discrimination power while preserving the ability to detect large anatomical deviations; as such they hold great potential in clinical applications, e.g. allowing for smaller sample sizes or shorter trials to establish disease specific biomarkers or to quantify drug effects.

This work reports the use of diffusion tensor magnetic resonance tractography to visualize the three-dimensional (3D) structure of the major white matter fasciculi within living human brain. Specifically, we applied this technique to visualize in vivo (i) the superior longitudinal (arcuate) fasciculus, (ii) the inferior longitudinal fasciculus, (iii) the superior fronto-occipital (subcallosal) fasciculus, (iv) the inferior frontooccipital fasciculus, (v) the uncinate fasciculus, (vi) the cingulum, (vii) the anterior commissure, (viii) the corpus callosum, (ix) the internal capsule, and (x) the fornix. These fasciculi were first isolated and were then interactively displayed as a 3D-rendered object. The virtual tract maps obtained in vivo using this approach were faithful to the classical descriptions of white matter anatomy that have previously been documented in postmortem studies. Since we have been able to interactively delineate and visualize white matter fasciculi over their en...

Myelin damage, as seen in multiple sclerosis (MS) and other demyelinating diseases, impairs axonal conduction and can also be associated with axonal degeneration. Accurate assessments of these conditions may be highly beneficial in evaluating and selecting therapeutic strategies for patient management. Recently, an analytical approach examining diffusion tensor imaging (DTI) derived parameters has been proposed to assess the extent of axonal damage, demyelination, or both. The current study uses the well-characterized cuprizone model of experimental demyelination and remyelination of corpus callosum in mouse brain to evaluate the ability of DTI parameters to detect the progression of myelin degeneration and regeneration. Our results demonstrate that the extent of increased radial diffusivity reflects the severity of demyelination in corpus callosum of mouse brain affected by cuprizone treatment. Subsequently, radial diffusivity decreases with the progression of remyelination. Furthermore, radial diffusivity changes were specific to the time course of changes in myelin integrity as distinct from axonal injury, which was detected by BAPP immunostaining and shown to be most extensive prior to demyelination. Radial diffusivity offers a specific assessment of demyelination and remyelination, as distinct from acute axonal damage. D

The advent of functional neuroimaging has allowed tremendous advances in our understanding of brain-language relationships, in addition to generating substantial empirical data on this subject in the form of thousands of activation peak coordinates reported in a decade of language studies. We performed a large-scale meta-analysis of this literature, aimed at defining the composition of the phonological, semantic, and sentence processing networks in the frontal, temporal, and inferior parietal regions of the left cerebral hemisphere. For each of these language components, activation peaks issued from relevant component-specific contrasts were submitted to a spatial clustering algorithm, which gathered activation peaks on the basis of their relative distance in the MNI space. From a sample of 730 activation peaks extracted from 129 scientific reports selected among 260, we isolated 30 activation clusters, defining the functional fields constituting three distributed networks of frontal and temporal areas and revealing the functional organization of the left hemisphere for language. The functional role of each activation cluster is discussed based on the nature of the tasks in which it was involved. This meta-analysis sheds light on several contemporary issues, notably on the fine-scale functional architecture of the inferior frontal gyrus for phonological and semantic processing, the evidence for an elementary audio-motor loop involved in both comprehension and production of syllables including the primary auditory areas and the motor mouth area, evidence of areas of overlap between phonological and semantic processing, in particular at the location of the selective human voice area that was the seat of partial overlap of the three language components, the evidence of a cortical area in the pars opercularis of the inferior frontal gyrus dedicated to syntactic processing and in the posterior part of the superior temporal gyrus a region selectively activated by sentence and text processing, and the hypothesis that different working memory perception-actions loops are identifiable for the different language components. These results argue for large-scale architecture networks rather than modular organization of language in the left hemisphere.

Neuroimaging and neuropsychological studies have implicated left inferior prefrontal cortex (LIPC) in both semantic and phonological processing. In this study, functional magnetic resonance imaging was used to examine whether separate LIPC regions participate in each of these types of processing. Performance of a semantic decision task resulted in extensive LIPC activation compared to a perceptual control task. Phonological processing of words and pseudowords in a syllable-counting task resulted in activation of the dorsal aspect of the left inferior frontal gyrus near the inferior frontal sulcus (BA 44/45) compared to a perceptual control task, with greater activation for nonwords compared to words. In a direct comparison of semantic and phonological tasks, semantic processing preferentially activated the ventral aspect of the left inferior frontal gyrus (BA 47/45). A review of the literature demonstrated a similar distinction between left prefrontal regions involved in semantic processing and phonological/lexical processing. The results suggest that a distinct region in the left inferior frontal cortex is involved in semantic processing, whereas other regions may subserve phonological processes engaged during both semantic and phonological tasks. 1999

Diffusion-weighted magnetic resonance imaging can provide information related to the arrangement of white matter fibers. The diffusion tensor is the model most commonly used to derive the orientation of the fibers within a voxel. However, this model has been shown to fail in regions containing several fiber populations with distinct orientations. A number of alternative models have been suggested, such as multiple tensor fitting, q-space, and Q-ball imaging. However, each of these has inherent limitations. In this study, we propose a novel method for estimating the fiber orientation distribution directly from high angular resolution diffusion-weighted MR data without the need for prior assumptions regarding the number of fiber populations present. We assume that all white matter fiber bundles in the brain share identical diffusion characteristics, thus implicitly assigning any differences in diffusion anisotropy to partial volume effects. The diffusion-weighted signal attenuation measured over the surface of a sphere can then be expressed as the convolution over the sphere of a response function (the diffusion-weighted attenuation profile for a typical fiber bundle) with the fiber orientation density function (ODF). The fiber ODF (the distribution of fiber orientations within the voxel) can therefore be obtained using spherical deconvolution. The properties of the technique are demonstrated using simulations and on data acquired from a volunteer using a standard 1.5-T clinical scanner. The technique can recover the fiber ODF in regions of multiple fiber crossing and holds promise for applications such as tractography. D

This article discusses general modeling of multisubject and/or multisession FMRI data. In particular, we show that a two-level mixed-effects model (where parameters of interest at the group level are estimated from parameter and variance estimates from the single-session level) can be made equivalent to a single complete mixed-effects model (where parameters of interest at the group level are estimated directly from all of the original single sessions' time series data) if the (co-)variance at the second level is set equal to the sum of the (co-)variances in the single-level form, using the BLUE with known covariances. This result has significant implications for group studies in FMRI, since it shows that the group analysis requires only values of the parameter estimates and their (co-)variance from the first level, generalizing the well-established "summary statistics" approach in FMRI. The simple and generalized framework allows different prewhitening and different first-level regressors to be used for each subject. The framework incorporates multiple levels and cases such as repeated measures, paired or unpaired t tests and F tests at the group level; explicit examples of such models are given in the article. Using numerical simulations based on typical first-level covariance structures from real FMRI data we demonstrate that by taking into account lower-level covariances and heterogeneity a substantial increase in higher-level Z score is possible.

Tractography based on diffusion tensor imaging (DTI) allows visualization of white matter tracts. In this study, protocols to reconstruct eleven major white matter tracts are described. The protocols were refined by several iterations of intra-and inter-rater measurements and identification of sources of variability. Reproducibility of the established protocols was then tested by raters who did not have previous experience in tractography. The protocols were applied to a DTI database of adult normal subjects to study size, fractional anisotropy (FA), and T 2 of individual white matter tracts. Distinctive features in FA and T 2 were found for the corticospinal tract and callosal fibers. Hemispheric asymmetry was observed for the size of white matter tracts projecting to the temporal lobe. This protocol provides guidelines for reproducible DTI-based tract-specific quantification.

Alerting, orienting, and executive control are widely thought to be relatively independent aspects of attention that are linked to separable brain regions. However, neuroimaging studies have yet to examine evidence for the anatomical separability of these three aspects of attention in the same subjects performing the same task. The attention network test (ANT) examines the effects of cues and targets within a single reaction time task to provide a means of exploring the efficiency of the alerting, orienting, and executive control networks involved in attention. It also provides an opportunity to examine the brain activity of these three networks as they operate in a single integrated task. We used event-related functional magnetic resonance imaging (fMRI) to explore the brain areas involved in the three attention systems targeted by the ANT. The alerting contrast showed strong thalamic involvement and activation of anterior and posterior cortical sites. As expected, the orienting contrast activated parietal sites and frontal eye fields. The executive control network contrast showed activation of the anterior cingulate along with several other brain areas. With some exceptions, activation patterns of these three networks within this single task are consistent with previous fMRI studies that have been studied in separate tasks. Overall, the fMRI results suggest that the functional contrasts within this single task differentially activate three separable anatomical networks related to the components of attention. D

Intersubject variability and subtle differences in experimental design can lead to variable results in studies of cognitive processes such as reading. To accurately identify the neural processes associated with cognition and sensorimotor processing, meta-analytic methods capable of identifying areas of consistent activation among studies are useful. This paper describes a novel approach for combining published neuroimaging results from multiple studies, designed to maximize the quantification of interstudy concordance while minimizing the subjective aspects of meta-analysis. In this method, a localization probability distribution was modeled for each activation focus obtained from 11 PET studies of reading single words aloud, and the union of these distributions was taken to yield an activation likelihood estimate map for the brain. Significance was assessed via permutation analysis of randomly generated sets of foci. Regions of significant concordance were identified in bilateral motor and superior temporal cortices, pre-SMA, left fusiform gyrus, and the cerebellum. These meta-analytic results were validated by comparison with new fMRI data on aloud word reading in normal adult subjects. Excellent correspondence between the two statistical maps was observed, with fMRI maxima lying close to all meta-analysis peaks and statistical values at the peaks identified by the two techniques correlating strongly. This close correspondence between PET meta-analysis and fMRI results also demonstrates the validity of using fMRI for the study of language tasks involving overt speech responses. Advantages of this automated meta-analysis technique include quantification of the level of concordance at all brain locations and the provision for use of a threshold for statistical significance of concordance.

Brain registration to a stereotaxic atlas is an effective way to report anatomic locations of interest and to perform anatomic quantification. However, existing stereotaxic atlases lack comprehensive coordinate information about white matter structures. In this paper, white matter specific atlases in stereotaxic coordinates are introduced. As a reference template, the widely-used ICBM-152 was used. The atlas contains fiber orientation maps and hand-segmented white matter parcellation maps based on diffusion tensor imaging (DTI). Registration accuracy by linear and nonlinear transformation was measured, and automated template-based white matter parcellation was tested. The results showed high correlation between the manual ROI-based and the automated approaches for normal adult populations. The atlases are freely available and believed to be a useful resource as a target template and for automated parcellation methods.

Bayesian model selection (BMS) is a powerful method for determining the most likely among a set of competing hypotheses about the mechanisms that generated observed data. BMS has recently found widespread application in neuroimaging, particularly in the context of dynamic causal modelling (DCM). However, so far, combining BMS results from several subjects has relied on simple (fixed effects) metrics, e.g. the group Bayes factor (GBF), that do not account for group heterogeneity or outliers. In this paper, we compare the GBF with two random effects methods for BMS at the between-subject or group level. These methods provide inference on model-space using a classical and Bayesian perspective respectively. First, a classical (frequentist) approach uses the log model evidence as a subject-specific summary statistic. This enables one to use analysis of variance to test for differences in log-evidences over models, relative to intersubject differences. We then consider the same problem in Bayesian terms and describe a novel hierarchical model, which is optimised to furnish a probability density on the models themselves. This new variational Bayes method rests on treating the model as a random variable and estimating the parameters of a Dirichlet distribution which describes the probabilities for all models considered. These probabilities then define a multinomial distribution over model space, allowing one to compute how likely it is that a specific model generated the data of a randomly chosen subject as well as the exceedance probability of one model being more likely than any other model. Using empirical and synthetic data, we show that optimising a conditional density of the model probabilities, given the log-evidences for each model over subjects, is more informative and appropriate than both the GBF and frequentist tests of the log-evidences. In particular, we found that the hierarchical Bayesian approach is considerably more robust than either of the other approaches in the presence of outliers. We expect that this new random effects method will prove useful for a wide range of group studies, not only in the context of DCM, but also for other modelling endeavours, e.g. comparing different source reconstruction methods for EEG/MEG or selecting among competing computational models of learning and decision-making.

In vivo MRI-derived measurements of human cerebral cortex thickness are providing novel insights into normal and abnormal neuroanatomy, but little is known about their reliability. We investigated how the reliability of cortical thickness measurements is affected by MRI instrument-related factors, including scanner field strength, manufacturer, upgrade and pulse sequence. Several data processing factors were also studied. Two test -retest data sets were analyzed: 1) 15 healthy older subjects scanned four times at 2-week intervals on three scanners; 2) 5 subjects scanned before and after a major scanner upgrade. Within-scanner variability of global cortical thickness measurements was <0.03 mm, and the point-wise standard deviation of measurement error was approximately 0.12 mm. Variability was 0.15 mm and 0.17 mm in average, respectively, for cross-scanner (Siemens/GE) and cross-field strength (1.5 T/3 T) comparisons. Scanner upgrade did not increase variability nor introduce bias. Measurements across field strength, however, were slightly biased (thicker at 3 T). The number of (single vs. multiple averaged) acquisitions had a negligible effect on reliability, but the use of a different pulse sequence had a larger impact, as did different parameters employed in data processing. Sample size estimates indicate that regional cortical thickness difference of 0.2 mm between two different groups could be identified with as few as 7 subjects per group, and a difference of 0.1 mm could be detected with 26 subjects per group. These results demonstrate that MRI-derived cortical thickness measures are highly reliable when MRI instrument and data processing factors are controlled but that it is important to consider these factors in the design of multi-site or longitudinal studies, such as clinical drug trials. D

Detecting artifacts produced in EEG data by muscle activity, eye blinks and electrical noise is a common and important problem in EEG research. It is now widely accepted that independent component analysis (ICA) may be a useful tool for isolating artifacts and/or cortical processes from electroencephalographic (EEG) data. We present results of simulations demonstrating that ICA decomposition, here tested using three popular ICA algorithms, Infomax, SOBI, and FastICA, can allow more sensitive automated detection of small nonbrain artifacts than applying the same detection methods directly to the scalp channel data. We tested the upper bound performance of five methods for detecting various types of artifacts by separately optimizing and then applying them to artifact-free EEG data into which we had added simulated artifacts of several types, ranging in size from thirty times smaller (-50 dB) to the size of the EEG data themselves (0 dB). Of the methods tested, those involving spectral thresholding were most sensitive. Except for muscle artifact detection where we found no gain of using ICA, all methods proved more sensitive when applied to the ICA-decomposed data than applied to the raw scalp data: the mean performance for ICA was higher and situated at about two standard deviations away from the performance distribution obtained on raw data. We note that ICA decomposition also allows simple subtraction of artifacts accounted for by single independent components, and/or separate and direct examination of the decomposed non-artifact processes themselves.

Brain maturation is a complex process that continues well beyond infancy, and adolescence is thought to be a key period of brain rewiring. To assess structural brain maturation from childhood to adulthood, we charted brain development in subjects aged 5 to 30 years using diffusion tensor magnetic resonance imaging, a novel brain imaging technique that is sensitive to axonal packing and myelination and is particularly adept at virtually extracting white matter connections. Age-related changes were seen in major white matter tracts, deep gray matter, and subcortical white matter, in our large (n = 202), age-distributed sample. These diffusion changes followed an exponential pattern of maturation with considerable regional variation. Differences observed in developmental timing suggest a pattern of maturation in which areas with fronto-temporal connections develop more slowly than other regions. These in vivo results expand upon previous postmortem and imaging studies and provide quantitative measures indicative of the progression and magnitude of regional human brain maturation.

We report a functional neuroimaging study with positron emission tomography (PET) in which six healthy adult volunteers were scanned while watching silent computer-presented animations. The characters in the animations were simple geometrical shapes whose movement patterns selectively evoked mental state attribution or simple action description.

We present a novel skull-stripping algorithm based on a hybrid approach that combines watershed algorithms and deformable surface models. Our method takes advantage of the robustness of the former as well as the surface information available to the latter. The algorithm first localizes a single white matter voxel in a T1-weighted MRI image, and uses it to create a global

We performed quantitative meta-analyses on 65 neuroimaging studies of emotion. In an earlier report (NeuroImage 16 (2002), 331). we examined the effects of induction method, specific emotions, and cognitive demand in emotional tasks. This paper focuses on the effects of emotional valence (positive vs negative and approach vs withdrawal) and gender on regional brain activations, with particular emphasis on hypotheses concerning lateralization of brain function in emotion. Overall, we found no support for the hypothesis of overall right-lateralization of emotional function, and limited support for valence-specific lateralization of emotional activity in frontal cortex. In addition, we found that males showed more lateralization of emotional activity, and females showed more brainstem activation in affective paradigms. The study provides evidence that lateralization of emotional activity is more complex and region-specific than predicted by previous theories of emotion and the brain.

Human total brain size is consistently reported to be ~8-10% larger in males, although consensus on regionally-specific differences is weak. Here, in the largest longitudinal pediatric neuroimaging study reported to date (829 scans from 387 subjects, ages 3 to 27 years), we demonstrate the importance of examining size-by-age trajectories of brain development rather than group averages across broad age ranges when assessing sexual dimorphism. Using magnetic resonance imaging (MRI) we found robust male/female differences in the shapes of trajectories with total cerebral volume peaking at age 10.5 in females and 14.5 in males. White matter increases throughout this 24 year period with males having a steeper rate of increase during adolescence. Both cortical and subcortical gray matter trajectories follow an inverted U shaped path with peak sizes 1 to 2 years earlier in females. These sexually dimorphic trajectories confirm the importance of longitudinal data in studies of brain development and underline the need to consider sex matching in studies of brain development.

Increased intraindividual variability (IIV) is a hallmark of disorders of attention. Recent work has linked these disorders to abnormalities in a "default mode" network, comprising brain regions routinely deactivated during goal-directed cognitive tasks. Findings from a study of the neural basis of attentional lapses suggest that a competitive relationship between the "task-negative" default mode network and regions of a "task-positive" attentional network is a potential locus of dysfunction in individuals with increased IIV. Resting state studies have shown that this competitive relationship is intrinsically represented in the brain, in the form of a negative correlation or antiphase relationship between spontaneous activity occurring in the two networks. We quantified the negative correlation between these two networks in 26 subjects, during active (Eriksen flanker task) and resting state scans. We hypothesized that the strength of the negative correlation is an index of the degree of regulation of activity in the default mode and task-positive networks and would be positively related to consistent behavioral performance. We found that the strength of the correlation between the two networks varies across individuals. These individual differences appear to be behaviorally relevant, as interindividual variation in the strength of the correlation was significantly related to individual differences in response time variability: the stronger the negative correlation (i.e., the closer to 180°antiphase), the less variable the behavioral performance. This relationship was moderately consistent across resting and task conditions, suggesting that the measure indexes moderately stable individual differences in the integrity of functional brain networks. We discuss the implications of these findings for our understanding of the behavioral significance of spontaneous brain activity, in both healthy and clinical populations.

Longitudinal and multi-site clinical studies create the imperative to characterize and correct technological sources of variance that limit image reproducibility in high-resolution structural MRI studies, thus facilitating precise, quantitative, platform-independent, multi-site evaluation. In this work, we investigated the effects that imaging gradient non-linearity have on reproducibility of multi-site human MRI. We applied an image distortion correction method based on spherical harmonics description of the gradients and verified the accuracy of the method using phantom data. The correction method was then applied to the brain image data from a group of subjects scanned twice at multiple sites having different 1.5 T platforms. Within-site and across-site variability of the image data was assessed by evaluating voxel-based image intensity reproducibility. The image intensity reproducibility of the human brain data was significantly improved with distortion correction, suggesting that this method may offer improved reproducibility in morphometry studies. We provide the source code for the gradient distortion algorithm together with the phantom data. D

This meta-analysis explores the role of the mirror and mentalizing systems in the understanding of other people's action goals. Based on over 200 fMRI studies, this analysis demonstrates that the mirror systemconsisting of the anterior intraparietal sulcus and the premotor cortexis engaged when one perceives articulated motions of body parts irrespective of their sensory (visual or auditory) or verbal format as well as when the perceiver executes them. This confirms the matching role of the mirror system in understanding biological action. Observation of whole-body motions and gaze engage the posterior superior temporal sulcus and most likely reflects an orientation response in line with the action or attention of the observed actor. In contrast, the mentalizing systemconsisting of the temporo-parietal junction, the medial prefrontal cortex and the precuneusis activated when behavior that enables inferences to be made about goals, beliefs or moral issues is presented in abstract terms (e.g., verbal stories or geometric shapes) and there is no perceivable biological motion of body parts. A striking overlap of brain activity at the temporo-parietal junction between social inferences and other, non-social observations (e.g., Posner's cuing task) suggests that this area computes the orientation or direction of the behavior in order to predict its likely end-state (or goal). No conclusions are drawn about the specific functionality of the precuneus and the medial prefrontal cortex. Because the mirror and mentalizing systems are rarely concurrently active, it appears that neither system subserves the other. Rather, they are complementary. There seems, however, to be a transition from the mirror to the mentalizing system even when body-part motions are observed by perceivers who are consciously deliberating about the goals of others and their behavioral executions, such as when perceived body motions are contextually inconsistent, implausible or pretended.

Atlas normalization, as commonly used by functional data analysis, provides an automated solution to the widely encountered problem of correcting for head size variation in regional and whole-brain morphometric analyses, so long as an age-and population-appropriate target atlas is used. In the present article, we develop and validate an atlas normalization procedure for head size correction using manual total intracranial volume (TIV) measurement as a reference. The target image used for atlas transformation consisted of a merged young and old-adult template specifically created for cross age-span normalization. Automated atlas transformation generated the Atlas Scaling Factor (ASF) defined as the volume-scaling factor required to match each individual to the atlas target. Because atlas normalization equates head size, the ASF should be proportional to TIV. A validation analysis was performed on 147 subjects to evaluate ASF as a proxy for manual TIV measurement. In addition, 19 subjects were imaged on multiple days to assess test -retest reliability. Results indicated that the ASF was (1) equivalent to manual TIV normalization (r = 0.93), (2) reliable across multiple imaging sessions (r = 1.00; mean absolute percentage of difference = 0.51%), (3) able to correct between-gender head size differences, and (4) minimally biased in demented older adults with marked atrophy. Hippocampal volume differences between nondemented (n = 49) and demented (n = 50) older adults (measured manually) were equivalent whether corrected using manual TIV or automated ASF (effect sizes of 1.29 and 1.46, respectively). To provide normative values, ASF was used to automatically derive estimated TIV (eTIV) in 335 subjects aged 15 -96 including both clinically characterized nondemented (n = 77) and demented (n = 90) older adults. Differences in eTIV between nondemented and demented groups were negligible, thus failing to support the hypothesis that large premorbid brain size moderates Alzheimer's disease. Gender was the only robust factor that influenced eTIV. Men showed an approximately f12% larger eTIV than women. These results demonstrate that atlas normalization using appropriate template images provides a robust, automated method for head size correction that is equivalent to manual TIV correction in studies of aging and dementia. Thus, atlas normalization provides a common framework for both morphometric and functional data analysis.

There is a growing appreciation of the importance of nonlinearities in evoked responses in fMRI, particularly with the advent of event-related fMRI. These nonlinearities are commonly expressed as interactions among stimuli that can lead to the suppression and increased latency of responses to a stimulus that are incurred by a preceding stimulus. We have presented previously a model-free characterization of these effects using generic techniques from nonlinear system identification, namely a Volterra series formulation. At the same time Buxton et al. (1998) described a plausible and compelling dynamical model of hemodynamic signal transduction in fMRI. Subsequent work by Mandeville et al. (1999) provided important theoretical and empirical constraints on the form of the dynamic relationship between blood flow and volume that underpins the evolution of the fMRI signal. In this paper we combine these system identification and model-based approaches and ask whether the Balloon model is sufficient to account for the nonlinear behaviors observed in real time series. We conclude that it can, and furthermore the model parameters that ensue are biologically plausible. This conclusion is based on the observation that the Balloon model can produce Volterra kernels that emulate empirical kernels. To enable this evaluation we had to embed the Balloon model in a hemodynamic input-state-output model that included the dynamics of perfusion changes that are contingent on underlying synaptic activation. This paper presents (i) the full hemodynamic model (ii), how its associated Volterra kernels can be derived, and (iii) addresses the model's validity in relation to empirical nonlinear characterisations of evoked responses in fMRI and other neurophysiological constraints.

Several tracing studies have established a topographical distribution of fiber connections to the cortex in midsagittal cross-sections of the corpus callosum (CC). The most prominent example is Witelson's scheme, which defines five vertical partitions mainly based on primate data. Conventional MRI of the human CC does not reveal morphologically discernable structures, although microscopy techniques identified myelinated axons with a relatively small diameter in the anterior and posterior third of the CC as opposed to thick fibers in the midbody and posterior splenium. Here, we applied diffusion tensor imaging (DTI) in conjunction with a tract-tracing algorithm to gain cortical connectivity information of the CC in individual subjects. With DTI-based tractography, we distinguished five vertical segments of the CC, containing fibers projecting into prefrontal, premotor (and supplementary motor), primary motor, and primary sensory areas as well as into parietal, temporal, and occipital cortical areas. Striking differences to Witelson's classification were recognized in the midbody and anterior third of the CC. In particular, callosal motor fiber bundles were found to cross the CC in a much more posterior location than previously indicated. Differences in water mobility were found to be in qualitative agreement with differences in the microstructure of transcallosal fibers yielding the highest anisotropy in posterior regions of the CC. The lowest anisotropy was observed in compartments assigned to motor and sensory cortical areas. In conclusion, DTI-based fiber tractography of healthy human subjects suggests a modification of the widely accepted Witelson scheme and a new classification of vertical CC partitions.

There is growing interest regarding the role of the right inferior frontal gyrus (RIFG) during a particular form of executive control referred to as response inhibition. However, tasks used to examine neural activity at the point of response inhibition have rarely controlled for the potentially confounding effects of attentional demand. In particular, it is unclear whether the RIFG is specifically involved in inhibitory control, or is involved more generally in the detection of salient or task relevant cues. The current fMRI study sought to clarify the role of the RIFG in executive control by holding the stimulus conditions of one of the most popular response inhibition tasks–the Stop Signal Task–constant, whilst varying the response that was required on reception of the stop signal cue. Our results reveal that the RIFG is recruited when important cues are detected, regardless of whether that detection is followed by the inhibition of a motor response, the generation of a motor response, or no external response at all.