PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Neurosci. Author manuscript; available in PMC 2010 August 4.
Published in final edited form as:
PMCID: PMC2915894
NIHMSID: NIHMS187483

Orbitofrontal Cortex Lesions Alter Anxiety-Related Activity in the Primate Bed Nucleus of Stria Terminalis

Abstract

In children, behavioral inhibition (BI) in response to potential threat predicts the development of anxiety and affective disorders and primate lesion studies suggest involvement of the orbtiofrontal cortex (OFC) in mediating BI. Lesion studies are essential for establishing causality in brain-behavior relationships, but should be interpreted cautiously because the impact of a discrete lesion on a complex neural circuit extends beyond the lesion location. Complementary functional imaging methods assessing how lesions influence other parts of the circuit can aid in precisely understanding how lesions affect behavior. Using this combination of approaches in monkeys, we found that OFC lesions concomitantly alter BI and metabolism in the Bed Nucleus of Stria Terminalis (BNST) region and that individual differences in BNST activity predict BI. Thus it appears that an important function of the OFC in response to threat is to modulate the BNST which may more directly influence the expression of BI.

Keywords: Pet (Positron Emission Tomography), Anxiety, Behavior, Prefrontal Cortex, Stria Terminalis, Primate

Introduction

In human children, laboratory-based measures of behavioral inhibition (BI) in response to unfamiliar individuals and/or novelty are markers of anxious temperament and are early predictors of the subsequent development of anxiety and affective disorders (N. A. Fox, Henderson, Marshall, Nichols, & Ghera, 2005). Because mechanistic studies cannot be performed in children, we developed and validated an animal model of trait-like BI in young rhesus monkeys (Kalin & Shelton, 2003). In a previous study, we demonstrated that OFC lesions decreased the expression of BI (Kalin, Shelton, & Davidson, 2007). To understand the specific role the OFC plays in the neural circuit that underlies BI, we assessed regional brain metabolic activity in freely behaving OFC-lesioned animals and cage-mate controls using high-resolution Positron Emission Tomography (PET).

The OFC is part of the neural circuit involved with emotion processing, and appears critical for guiding and maintaining emotional responses (Wallis, 2007; Milad & Rauch, 2007). Recent work has shown OFC neurons to encode and store the relative value of different choices, and to be required for reversing the value of a cue (Wallis, 2007). Additionally, the OFC is bi-directionally connected with the subcortical components of the circuit that mediates fear and anxiety, and is thought to play an important role in the pathophysiology of anxiety and affective disorders (Kalin & Shelton, 2003; Milad & Rauch, 2007). Studies in both OFC- as well as amygdala-lesioned monkeys exposed to the no-eye-contact condition (NEC) of the human intruder paradigm have demonstrated replicable reductions in threat-induced freezing, which is a measure of BI (Kalin et al., 2007; Kalin, Shelton, & Davidson, 2004; Elizabeth A Murray & Izquierdo, 2007; Machado & Bachevalier, 2008; though see: Izquierdo, Suda, & Elisabeth A Murray, 2005). The results of functional imaging studies in monkeys in relation to BI have not been entirely consistent with the results from the lesion studies. For example, [F-18]fluoro-deoxyglucose (FDG) PET imaging in intact animals failed to find a significant correlation between OFC metabolism and either NEC-associated freezing behavior or anxious temperament, a composite measure such that individuals with higher levels of anxious temperament have increased freezing behavior, decreased vocalizations and increased cortisol (Kalin, Shelton, A. S. Fox, Oakes, & Davidson, 2005; A. S. Fox, Shelton, Oakes, Davidson, & Kalin, 2008). Instead, these functional imaging studies highlighted subcortical stress-circuits, in which the amygdala and a region encompassing the Bed Nucleus of Stria Terminalis (BNST; a part of the extended amygdala) predicted individual differences in the duration of monkeys’ threat-associated freezing behavior and anxious temperament (Kalin et al., 2005; A. S. Fox et al., 2008). These data suggest the influences of the OFC on anxiety-related behavior may be indirect, as OFC lesions alter freezing but in intact animals brain activity in this region does not predict individual differences in freezing. The imaging data suggest that the BNST and amygdala are the proximate regions that mediate freezing. Based on these data, we hypothesized that lesions of the OFC alter freezing and BI by modulating the function of the amygdala and BNST, regions that more directly mediate anxiety and fear.

Methods

Subjects

Six male rhesus monkeys (Macaca mulatta) that received bilateral >60% OFC aspiration lesions (which included Areas 11, 12, 13 and 14) discussed in Kalin, Shelton & Davidson (2007) and 6 unoperated cage-mate controls were the focus of this study. Further characterizations of lesion extent based on Walker areas are not presented due to concerns about the reliability of using MRIs post-lesion to specifically quantitate damage to subregions of OFC. Each animal was injected with FDG, exposed to the NEC condition, anesthetized, and PET scanned at both time-1, where all animals were unoperated, and at time-2, which was post-lesion for the experimental group (Kalin et al., 2007). Of the 6 pairs of animals, one pair was excluded because we were unable to obtain quality FDG-PET data for the OFC-lesioned animal at time-1. The ten animals included in the analysis (age: mean=2.86 years) were matched for age, and were housed as pairs with each experimental animal living with a control animal. All testing was done in accordance with institutional guidelines.

Positron Emission Tomography and Behavioral Measurements

The six lesioned animals and their age-matched cage-mate controls were examined twice with FDG-PET (time-1, pre lesion for experimental subjects and time-2, post lesion for experimental subjects) during exposure to the NEC condition of the human intruder paradigm. Individual animals were injected with FDG intravenously and placed in a test cage, where they were exposed to a modified version of the NEC condition of the human intruder paradigm. In the modified human intruder paradigm a human experimenter entered the room and stood with their profile to the monkey, making no eye contact (NEC) with the animal 2.5 meters from the test cage. To prevent habituation during the 30-minutes of FDG uptake, the human remained in the test room for 10 minutes, left the test room for 5 minutes, returned for 5 minutes, left for five-minutes and returned for the remaining 5-minutes. Behavior was recorded on closed circuit video, scored by trained raters blind to condition and aggregated across the three periods of the NEC condition. Freezing was defined as a period of at least 3 seconds characterized by tense body posture without vocalizations and movement other than slow movements of the head. Cooing was defined as a vocalization made by rounding and pursing the lips with an increase then decrease in frequency and intensity. Locomotion was coded as of one or more full steps at any speed. The freezing results presented here are a subset of the results presented in Kalin et al., 2007(Kalin et al., 2007). After the 30-minutes of FDG uptake animals were anesthetized with ketamine (15 mg/kg) and were administered intramuscular atropine sulfate (.27 mg). They were then fitted with an endotracheal tube, to administer 1–2% isofluorane gas anesthesia. The animal was then placed in the microPET P4 scanner (Concorde Microsystems, Inc., Knoxville, Tennessee)(Tai et al., 2001). Sixty-minute emission PET scans were reconstructed using filtered back projection and reflect the integrated brain metabolism that occurred during the ~30-minutes of FDG uptake.

Pre-processing

Time-2 (post-surgical for the experimental group) PET scans were masked to exclude prefrontal regions where anatomy may have changed in the experimental group, and posterior brain regions which were misplaced outside of the PET scanner field of view in some of the pre-surgery images (see: Kalin et al., 2005 for additional information). Each animal’s time-2 masked images were transformed using a 6-parameter rigid-body transformation to match their pre-lesion PET scan. Transformations were applied to the unmasked time-2 PET scans and manually inspected by ASF to assure each subject’s time-2 PET scans was accurately aligned to the time-1 PET image. Non-linear transformations based on the time-1 images (computed using AIR5) were applied to the time-2 images, resulting in time-2 images in the standard space of Paxinos, Huang and Toga (Woods, Grafton, Holmes, Cherry, & Mazziotta, 1998; Paxinos, Huang, & Toga, 1999). Additional details on this registration can be found in Kalin et al., (2005). To facilitate inter-scan comparisons, images were globally scaled by adjusting the mean, based on a partial brain ROI created by intersecting the masked images used for alignment (Kalin et al., 2005; A. S. Fox et al., 2005). Each globally scaled image was smoothed using a 4mm3 Full-Width at Half Max (FWHM) gaussian kernel to account for slight across subject variation in anatomy and to enforce a normal distribution.

Statistical Analyses

Statistical analyses were performed using standard voxelwise regression techniques. To verify lesion efficacy in decreasing OFC metabolism, a region of interest based on the target lesion area was drawn on an MRI template, and relative glucose metabolism within this region was extracted for each subject. This effect was highly significant, suggesting our analysis was sensitive to the effects of the lesion (data not shown). In order to assess brain differences in non-OFC regions, analyses were performed using a voxelwise search within regions where data were accurately aligned to standard space. Posterior visual regions were excluded from the analyses, since this region was outside of the scanner FOV in some of the pre-lesion PET scans (see Kalin et al., 2005), and anterior regions (anterior to the genu of the corpus collosum) were excluded because we could not accurately account for lesion-induced deformations of brain structure within these regions. Within trusted regions, we identified regions greater than 8mm3 (the FWHM resolution of the scanner) that showed a significant lesion/no lesion by time-1/time-2 interaction (p<.05, two-tailed uncorrected) and a significant time-2 main effect of Group (lesion vs. no lesion) (p<.005, two-tailed uncorrected). Follow-up analyses were performed in SPSS, focusing on linear regressions on the mean values from the clusters identified in the voxelwise analysis. To test the effects of laterality we also extracted mean values from the opposing hemisphere of significant clusters, and performed a lesion/no lesion by time-1/time-2 by left-hemisphere/right-hemisphere interaction.

Results

We first examined the effects of OFC lesions on freezing behavior that occurred during the period of FDG uptake. Consistent with Kalin et al. (2007), we observed that OFC lesioned animals, compared to intact animals, exhibited significantly less freezing when tested during the time 2 FDG-uptake period (t(8)=2.371, p=.026; see supplementary note; Fig 1b)(Kalin et al., 2007). Further analyses revealed a significant increase in freezing duration within the control group (t(4)=4.848, p=.008), with no group differences at time-1 or any significant decrease in freezing behavior in the experimental group (p’s>.5; additional discussion can be found in supplementary information).

Figure 1
Orbital Frontal Cortex (OFC) Lesions Alter Freezing and Metabolism in the BNST Region

Based on the hypothesis that the OFC indirectly influences freezing by modulating amygdala and BNST activity, we predicted that OFC lesions would decrease activity in the amygdala and BNST. Results of the voxelwise analyses demonstrated that NEC-exposed OFC-lesioned animals had significantly (p<.005, two-tailed uncorrected time-2 effect & a p<.05, two-tailed uncorrected interaction) lower glucose metabolism in a region that overlapped with the BNST (Fig 1b) and additionally in the entorhinal cortex. The analysis also revealed increased activation in OFC- lesioned animals in regions of parietal cortex (LIP), mid-cingulate gyrus, and motor cortex (for a complete list of affected regions see table 1). Analyses of asymmetry revealed no group by time by hemisphere interaction in BNST metabolism (F(4,25)=1.44, p=.25), suggesting no significant differences in BNST laterality. Although we did not observe significant effects of the OFC lesions on amygdala metabolism, this null finding may result from a lack of statistical power as previously published associations between freezing behavior and amygdala metabolism suggest a relatively low zero-order correlation which gets stronger when combined with other measures not discussed in this manuscript (Kalin et al., 2005).

Table 1
Control minus Experimental Post-Lesion p<.005

To further examine the role of the BNST region in BI, we investigated the correlation between brain activity in the BNST region with freezing behavior at time 2. Because we hypothesized OFC lesions were mediated by the BNST, and would not alter the relationship between BNST and freezing, we predicted these data would replicate the previously observed positive correlation between BNST and freezing behavior (Kalin et al., 2005). Our findings demonstrated that mean metabolic activity in the BNST region that showed a group difference between lesioned and non lesioned animals strongly predicted freezing (spearman’s-rho = .782, p=0.008; controlling for age in both variables, Figure 2). To determine the specificity of this correlation to freezing behavior, we performed the same analyses correlating individual differences in BNST region activity with individual differences in affiliative vocalizations (coo-calls) and locomotion and found that neither were significantly correlated with BNST region metabolism (p’s>0.174). To ensure that the BNST region was accounting for the majority of the variance associated with freezing behavior, we entered each of the regions affected by the OFC lesions into a hierarchical linear regression with the BNST region. This analysis revealed that metabolic activity in no other region significantly accounted for additional variance in freezing behavior beyond that accounted for by the BNST region (F-Change’s<2.402, p’s>0.165). These analyses suggest that the BNST may be the critical brain region associated with the effect of the OFC lesions on decreasing freezing behavior. Additionally, this finding is consistent with previous research suggesting motor, parietal, and mid-cingulate regions may be related to specific operations associated with BI, such as alterations in locomotion and visual attention, rather than the altered motivation associated with BI itself (Gottlieb, 2007; Paus, 2001).

Figure 2
Bed Nucleus of Stria Terminalis (BNST) Region Metabolism Predicts Freezing

Discussion

Using a combined lesion and imaging approach in an established primate model, we demonstrate an interaction between the OFC and the BNST region that sheds light on how the OFC modulates trait-like anxiety-related responses in primates. OFC lesions altered freezing, a behavior that is critical and adaptive in helping an animal remain undiscovered in the face of a potential threat. The current data demonstrate that OFC lesion-induced changes in freezing could be accounted for by alterations in BNST activity.

Recent research has highlighted a number of important functions of the OFC that help elucidate its role in anxiety-related behavior. In particular, electrophysiological recordings in monkeys and functional imaging in humans have shown this region to track the relative value of choices and maintain this representation over time (Wallis, 2007; Kringelbach & Rolls, 2004). Additionally, lesion studies have demonstrated the OFC to be critical for updating the value of reward-related cues (Wallis, 2007). The OFC may respond similarly to threat-related cues such that in the NEC paradigm the OFC could be accumulating and maintaining the emotional representation of the human intruder to facilitate adaptive behavioral responses. Importantly, this is consistent with structural and functional differentiation of the OFC into “medial” and “orbital” OFC networks, with the orbital OFC (the region lesioned in this study) demonstrating more activation to cues that signal negative outcomes (Kringelbach & Rolls, 2004; Price, 2007). Based on the data from the current study, we further suggest that the behavioral actions induced by the perception of threat cues are implemented via the OFC’s influence on the BNST region.

The relation we observed between individual differences in BNST region metabolism and trait-like anxiety agrees with our previous work as well as rodent work examining the functional significance of BNST activity (Kalin et al., 2005; A. S. Fox et al., 2008; Walker et al., 2003). The BNST is part of the limbic system and basal forebrain, with the posterior edge near the central nucleus of the amygdala (CeA) and the anterior edge bordering on the nucleus accumbens (L Heimer et al., 1997). The similar development, cytoarchitecture and connectivity of the BNST and CeA have led researchers to group these regions into the broader construct of the “extended amygdala” (Lennart Heimer, 2003). Because the extended amygdala receives substantial input from other amygdala nuclei and projects to the brainstem and hypothalamic regions required to mount a stress-response, the extended amygdala is ideally situated to be the interface to the motor responses associated with emotional reactions (L Heimer et al., 1997). More recently, Walker and Davis have suggested a functional dissociation within the extended amygdala, such that the CeA is necessary for rapid-onset, short-duration, acute fear responding to specific threat, and the BNST is associated with slower-onset, longer-lasting responses associated with anxiety and BI (Walker et al., 2003). These data are consistent with the interpretation that OFC lesions alter metabolism in the BSNT-region, which in turn initiates BI measured in response to a relatively long (30-minute) exposure to the NEC condition. Importantly, our findings are also highly consistent with ongoing work in rodents demonstrating the BNST region to mediate individual differences in anxiety, and our work in monkeys suggesting that serotonin transporter availability in the BNST may represent a complementary mechanism for altering BNST activity (Duvarci, Bauer, & Paré, 2009; Oler et al., 2009).

Interestingly, anatomical tract tracing studies in primates have not observed the orbital OFC network to be directly connected to the BNST. Retrograde tracing studies measuring the direct projections to the orbital OFC have failed to label neurons in the BNST, and the anterograde studies investigating the projections from the OFC to the BNST have not been published (Price, 2007). Thus, it remains unclear if the OFC lesions are decreasing freezing behavior via direct or indirect connections to the BNST (e.g. via known thalamic, insular or amygdalar connections), and highlights the fact that these two regions represent only a portion of the functional network that underlies BI (Ghashghaei, Hilgetag, & Barbas, 2007).

Although this sample size is consistent with that of other studies investigating the effects of targeted lesions on behavior, it is substantially underpowered for analyses that attempt to understand the relation between individual differences in lesion extent with brain activity and behavior. Future work with larger sample sizes will facilitate an understanding of the relationships among lesion location, remaining functional activation within the lesioned areas, and lesion-induced behavioral alterations.

In conclusion, these data suggest that the OFC does not directly alter freezing per se, but rather the behavioral effects of OFC lesions result from the removal of descending modulatory inputs (direct or indirect) to the BNST region. It follows that alterations in OFC function may underlie difficulties in the adaptive regulation of anxiety, which is thought to be a cardinal feature of anxiety and affective disorders (Kalin & Shelton, 2003). In contrast to the strategy of using either selective lesion techniques or functional imaging to understand brain behavior relations, this study underscores the power of combining these methods to better understand the complex interactions among components of relevant neural circuits.

Supplementary Material

Supplementary Note

Acknowledgements

This work was supported by grants MH46729, MH69315, The HealthEmotions Research Institute, Meriter Hospital, and Training Grant T32-MH018931 to ASF. We are grateful to H. Van Valkenberg, T. Johnson, J. King, S. Mansavage, J. Droster, L. Greischar, A Jahn, T Johnstone, A Shackman, J Oler, and the staff at the Harlow Center for Biological Psychology and the National Primate Research Center at the University of Wisconsin for their technical support and Drs. Elisabeth Murray and Jocelyne Bachevalier for their neurosurgical advice.

References

  • Duvarci S, Bauer EP, Paré D. The bed nucleus of the stria terminalis mediates inter-individual variations in anxiety and fear. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience. 2009;29(33):10357–10361. doi: 10.1523/JNEUROSCI.2119-09.2009. [PMC free article] [PubMed]
  • Fox AS, Oakes TR, Shelton SE, Converse AK, Davidson RJ, Kalin NH. Calling for help is independently modulated by brain systems underlying goal-directed behavior and threat perception. Proceedings of the National Academy of Sciences of the United States of America. 2005;102(11):4176–4179. doi: 10.1073/pnas.0409470102. [PubMed]
  • Fox AS, Shelton SE, Oakes TR, Davidson RJ, Kalin NH. Trait-Like Brain Activity during Adolescence Predicts Anxious Temperament in Primates. PLoS ONE. 2008;3(7):e2570. doi: 10.1371/journal.pone.0002570. [PMC free article] [PubMed]
  • Fox NA, Henderson HA, Marshall PJ, Nichols KE, Ghera MM. Behavioral Inhibition: Linking Biology and Behavior within a Developmental Framework. Annual Review of Psychology. 2005;56(1):235–262. doi: 10.1146/annurev.psych.55.090902.141532. [PubMed]
  • Ghashghaei HT, Hilgetag CC, Barbas H. Sequence of information processing for emotions based on the anatomic dialogue between prefrontal cortex and amygdala. NeuroImage. 2007;34(3):905–923. doi: 10.1016/j.neuroimage.2006.09.046. [PMC free article] [PubMed]
  • Gottlieb J. From thought to action: the parietal cortex as a bridge between perception, action, and cognition. Neuron. 2007;53(1):9–16. doi: 10.1016/j.neuron.2006.12.009. [PubMed]
  • Heimer L, Harlan RE, Alheid GF, Garcia MM, de Olmos J. Substantia innominata: a notion which impedes clinical-anatomical correlations in neuropsychiatric disorders. Neuroscience. 1997;76(4):957–1006. [PubMed]
  • Heimer L. A new anatomical framework for neuropsychiatric disorders and drug abuse. The American Journal of Psychiatry. 2003;160(10):1726–1739. [PubMed]
  • Izquierdo A, Suda RK, Murray EA. Comparison of the effects of bilateral orbital prefrontal cortex lesions and amygdala lesions on emotional responses in rhesus monkeys. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience. 2005;25(37):8534–8542. doi: 10.1523/JNEUROSCI.1232-05.2005. [PubMed]
  • Kalin NH, Shelton SE. Nonhuman Primate Models to Study Anxiety, Emotion Regulation, and Psychopathology. Annals of the New York Academy of Sciences. 2003;1008:189–200. (Roots of Mental Illness in Children) doi: 10.1196/annals.1301.021. [PubMed]
  • Kalin NH, Shelton SE, Davidson RJ. Role of the Primate Orbitofrontal Cortex in Mediating Anxious Temperament. Biological Psychiatry. 2007;62(10):1134–1139. doi: 10.1016/j.biopsych.2007.04.004. [PubMed]
  • Kalin NH, Shelton SE, Davidson RJ. The Role of the Central Nucleus of the Amygdala in Mediating Fear and Anxiety in the Primate. J. Neurosci. 2004;24(24):5506–5515. doi: 10.1523/JNEUROSCI.0292-04.2004. [PubMed]
  • Kalin NH, Shelton SE, Fox AS, Oakes TR, Davidson RJ. Brain Regions Associated with the Expression and Contextual Regulation of Anxiety in Primates. Biological Psychiatry. 2005;58(10):796–804. doi: 10.1016/j.biopsych.2005.05.021. [PMC free article] [PubMed]
  • Kringelbach ML, Rolls ET. The functional neuroanatomy of the human orbitofrontal cortex: evidence from neuroimaging and neuropsychology. Progress in Neurobiology. 2004;72(5):341–372. doi: 10.1016/j.pneurobio.2004.03.006. [PubMed]
  • Machado CJ, Bachevalier J. Behavioral and hormonal reactivity to threat: effects of selective amygdala, hippocampal or orbital frontal lesions in monkeys. Psychoneuroendocrinology. 2008;33(7):926–941. doi: 10.1016/j.psyneuen.2008.04.012. [PMC free article] [PubMed]
  • Milad MR, Rauch SL. The Role of the Orbitofrontal Cortex in Anxiety Disorders. Annals of the New York Academy of Sciences. 2007;1121:546–561. (Linking affect to Action: Critical Contributions of the Orbitofrontal Cortex) doi: 10.1196/annals.1401.006. [PubMed]
  • Murray EA, Izquierdo A. Orbitofrontal Cortex and Amygdala Contributions to Affect and Action in Primates. Annals of the New York Academy of Sciences. 2007;1121:273–296. (Linking affect to Action: Critical Contributions of the Orbitofrontal Cortex) doi: 10.1196/annals.1401.021. [PubMed]
  • Oler JA, Fox AS, Shelton SE, Christian BT, Murali D, Oakes TR, et al. Serotonin transporter availability in the amygdala and bed nucleus of the stria terminalis predicts anxious temperament and brain glucose metabolic activity. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience. 2009;29(32):9961–9966. doi: 10.1523/JNEUROSCI.0795-09.2009. [PMC free article] [PubMed]
  • Paus T. Primate anterior cingulate cortex: where motor control, drive and cognition interface. Nature Reviews. Neuroscience. 2001;2(6):417–424. doi: 10.1038/35077500. [PubMed]
  • Paxinos G, Huang X, Toga A. The Rhesus Monkey Brain in Stereotaxic Coordinates. Academic Press; 1999.
  • Price JL. Definition of the Orbital Cortex in Relation to Specific Connections with Limbic and Visceral Structures and Other Cortical Regions. Annals of the New York Academy of Sciences. 2007;1121:54–71. (Linking affect to Action: Critical Contributions of the Orbitofrontal Cortex) doi: 10.1196/annals.1401.008. [PubMed]
  • Tai C, Chatziioannou A, Siegel S, Young J, Newport D, Goble RN, et al. Performance evaluation of the microPET P4: a PET system dedicated to animal imaging. Physics in Medicine and Biology. 2001;46(7):1845–1862. [PubMed]
  • Walker DL, Toufexis DJ, Davis M. Role of the bed nucleus of the stria terminalis versus the amygdala in fear, stress, and anxiety. European Journal of Pharmacology. 2003;463(1-3):199–216. [PubMed]
  • Wallis JD. Orbitofrontal Cortex and Its Contribution to Decision-Making. Annual Review of Neuroscience. 2007;30(1):31–56. doi: 10.1146/annurev.neuro.30.051606.094334. [PubMed]
  • Woods RP, Grafton ST, Holmes CJ, Cherry SR, Mazziotta JC. Automated image registration: I. General methods and intrasubject, intramodality validation. Journal of Computer Assisted Tomography. 1998;22(1):139–152. [PubMed]