limitations of phineas gage case study

Uncovering the Impact of Phineas Gage’s Accident on Psychology

Phineas Gage is a name that has become synonymous with studying psychology. His case has greatly interested researchers and students alike for years. Gage experienced a traumatic brain injury in 1848 when an iron rod was driven through his skull, destroying much of his frontal lobe. Despite the severity of his injury, Gage survived the accident, and his story has contributed significantly to our understanding of the brain and its functions.

Gage’s case is unique because it is one of the first recorded instances of a person surviving a severe brain injury. The accident profoundly impacted Gage’s personality, behavior, and cognitive abilities. Before the accident, Gage was described as a responsible and hard-working man. However, after the injury, he became impulsive, inconsistent, and unable to plan for the future. His case has been a subject of much debate and speculation, and it has contributed significantly to our understanding of the brain and its functions.

This article will explore Phineas Gage’s accident and its impact on psychology. We will examine his life before and after the injury, his personality and behavior changes, and the significance of his psychology case. We will also discuss how Gage’s case has contributed to our understanding of the brain and its functions.

The Life of Phineas Gage

Phineas Gage was born on July 9, 1823, in Grafton County, New Hampshire. He was the first of five children born to Jesse Eaton Gage and Hannah Trussell Gage. According to records, Phineas was a healthy and intelligent child who enjoyed playing outdoors and was well-liked by his peers.

Early Years

Phineas grew up in a modest household and attended school until he was 12. He then began working on his family’s farm, where he gained experience working with animals and machinery. In his late teens, he moved to Vermont to work as a farmhand and eventually found employment as a stagecoach driver.

Career as a Railroad Construction Foreman

In 1848, at 25, Phineas worked as a railroad construction foreman for the Rutland and Burlington Railroad Company. On September 13 that year, he suffered a traumatic brain injury when an iron rod measuring 43 inches long and 1.25 inches in diameter was accidentally driven through his skull. The rod entered below his left cheekbone and exited through the top of his head, destroying much of his frontal lobe.

Despite the severity of his injury, Phineas survived and could walk and talk within minutes of the accident. However, his personality and behavior were dramatically altered. He became impulsive, irritable, and unable to plan or make decisions. He also struggled with memory loss and had difficulty with social interactions.

After the accident, Phineas worked odd jobs and traveled around the country, becoming a celebrity due to his survival story. He eventually settled in San Francisco, where he worked as a longshoreman until his death on May 21, 1860.

Phineas Gage: The Accident

The incident.

We have all heard of Phineas Gage, who survived a devastating brain injury that changed his personality forever. But what exactly happened to him? In 1848, Gage was a 25-year-old railroad construction foreman working in Vermont. On September 13 that year, he was using a tamping iron to pack explosive powder into a hole when the powder ignited, sending the iron rocketing through his skull.

The iron, 43 inches long and weighing over 13 pounds, entered Gage’s head just below his left cheekbone and exited through the top of his skull, landing several yards away. Miraculously, Gage remained conscious throughout the ordeal and was able to speak within minutes of the accident.

Immediate Aftermath

The immediate aftermath of the accident was chaotic. Gage’s coworkers rushed to his aid and found him sitting up, blood pouring from his head. They took him to a nearby hotel, where a physician named Edward H. Williams examined him. Williams later described the scene: “I first noticed the wound upon the head before I alighted from my carriage, the brain pulsations being very distinct. Mr. Gage, during the time I was examining this wound, was relating how he was injured to the bystanders. I did not believe Mr. Gage’s statement then but thought he was deceived.”

Despite Williams’ skepticism, Gage’s account of the accident was accurate. He was eventually taken to his boarding house, where Dr. John Martyn Harlow attended. Harlow later wrote that “the iron entered the left side of the face, shattering the upper jaw, and passing back of the left eye, tearing the left lobe of the brain, and passing out at the top of the head, carrying with it a portion of the brain and other substances.”

Gage’s survival was miraculous, but it came at a significant cost. The iron had destroyed much of his left frontal lobe, which regulates emotions, personality, and decision-making. Gage’s behavior changed dramatically in the aftermath of the accident, and he became irritable, impulsive, and unpredictable. His story would become one of the most famous case studies in psychology and neuroscience, and it continues to fascinate researchers and laypeople today.

Medical and Psychological Impact

Phineas Gage’s accident had significant medical and psychological consequences. This section will discuss his injury’s physical consequences, behavioral changes, and long-term effects.

Physical Consequences

The iron rod that went through Phineas Gage’s skull damaged much of his frontal lobe, resulting in significant physical consequences. He lost his left eye and partially lost vision in his right eye. Additionally, he suffered from seizures and chronic headaches.

Behavioral Changes

Phineas Gage’s injury also resulted in significant behavioral changes. He became impulsive, irritable, and lacked empathy. He struggled with decision-making and planning, and his personality underwent a complete transformation. His friends and family reported that he was no longer the person he was before the accident.

Long-Term Effects

The long-term effects of Phineas Gage’s injury were significant. He could not return to his previous job as a railroad construction foreman and struggled to maintain employment. He became a circus attraction, traveling around the country as a living example of the effects of brain injury.

Phineas Gage’s case profoundly impacted the field of psychology, as it was one of the first documented cases of the link between brain damage and behavior. It helped researchers understand the role of the frontal lobe in decision-making, planning, and personality.

Significance in Psychology

Phineas Gage’s accident and its aftermath have significantly impacted the field of psychology. It has provided valuable insight into brain function, influenced neuropsychology development, and had implications for personality theory.

Insights into Brain Function

Gage’s injury provided a unique opportunity to study the relationship between the brain and behavior. It demonstrated that damage to specific brain areas can result in profound changes in personality and behavior. It also highlighted the importance of the prefrontal cortex in regulating social behavior, decision-making, and emotional regulation.

Influence on Neuropsychology

Gage’s case was one of the first documented cases of brain injury resulting in significant changes in behavior. It influenced the development of neuropsychology, a field that focuses on the relationship between brain function and behavior. Neuropsychologists use a variety of tests to assess cognitive function, including memory, attention, and language skills. They also use brain imaging techniques like magnetic resonance imaging (MRI) to study the brain’s structure and function.

Implications for Personality Theory

Gage’s case challenged the prevailing view of personality as fixed and unchanging. It demonstrated that brain injury can alter personality and that different parts of the brain are responsible for various aspects of personality. For example, damage to the prefrontal cortex can result in impulsivity, poor judgment, and emotional instability. This has led to theories that emphasize personality’s dynamic nature and the brain’s importance in shaping behavior.

Overall, Phineas Gage’s case has had a lasting impact on psychology. It has provided valuable insights into brain function, influenced neuropsychology development, and challenged prevailing views of personality.

Public Response and Legacy

Public perception.

After Phineas Gage’s accident, his story quickly spread throughout the public. People were fascinated by the idea that a simple iron rod could cause such a dramatic change in a person’s personality. However, the public’s understanding of Gage’s case was often oversimplified and exaggerated. Many people believed that Gage’s accident had turned him into a completely different person when, in reality, the changes were more subtle and complex.

Over time, Gage’s story became a cautionary tale about the dangers of head injuries. It was used to warn people about the potential consequences of traumatic brain injuries and to encourage them to take precautions to protect their heads.

Legacy in Science and Popular Culture

Phineas Gage’s case profoundly impacted the field of psychology and neuroscience. It was one of the first documented cases of a person with a brain injury that affected their personality and behavior. Gage’s case helped to establish the idea that different parts of the brain are responsible for various functions and that damage to these areas can cause specific changes in behavior.

Gage’s story has also had a lasting impact on popular culture. It has been referenced in countless books, movies, and TV shows and has become a symbol of the strange and mysterious workings of the human brain. In recent years, Gage’s case has even been used to promote the idea of neuroplasticity, which suggests that the brain can change and adapt throughout a person’s life.

Overall, Phineas Gage’s accident and its aftermath have left a lasting legacy in both the scientific and popular realms. While the public’s understanding of Gage’s case may be oversimplified, his story’s impact on psychology and neuroscience must be considered.

Frequently Asked Questions

Did phineas gage lose his eye.

No, Phineas Gage did not lose his eye in the accident. However, the iron rod that went through his skull damaged his left eye and caused him to lose vision in that eye.

When did Phineas Gage’s accident happen?

Phineas Gage’s accident happened on September 13, 1848. He was working on a railroad construction crew in Vermont when the iron rod he was using to tamp down blasting powder accidentally ignited the powder, causing the rod to shoot through his skull.

How long did Phineas Gage live after the accident?

Phineas Gage lived for another 12 years after the accident. He suffered from seizures and personality changes due to the damage to his brain, but he could continue living a relatively everyday life until he died in 1860.

How did Phineas Gage survive?

Phineas Gage’s survival is considered a medical miracle. The iron rod that went through his skull destroyed much of his frontal lobe, responsible for many essential functions such as decision-making, personality, and social behavior. However, the fact that the rod went through his brain in a relatively straight line and did not damage other vital areas likely contributed to his survival.

Why did Phineas Gage not feel pain?

It is unclear why Phineas Gage did not feel pain immediately after the accident. Some speculate that the damage to his brain may have affected his ability to perceive pain, while others believe that his body went into shock and he did not feel the pain at the time.

What happened to Phineas Gage, and how did it impact the field of psychology?

Phineas Gage’s accident and subsequent personality changes profoundly impacted the field of psychology. His case was one of the first to suggest a link between brain function and behavior, and it helped to establish the field of neuropsychology. Gage’s story also highlighted the importance of the frontal lobe in regulating personality and decision-making, and it continues to be studied by psychologists and neuroscientists today.

Leave a Comment Cancel reply

Save my name, email, and website in this browser for the next time I comment.

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here .

Loading metrics

Open Access

Peer-reviewed

Research Article

Mapping Connectivity Damage in the Case of Phineas Gage

* E-mail: [email protected]

Affiliation Laboratory of Neuro Imaging (LONI), Department of Neurology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America

Affiliation Surgical Planning Laboratory, Department of Radiology, Brigham & Women's Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

• John Darrell Van Horn, 
• Andrei Irimia, 
• Carinna M. Torgerson, 
• Micah C. Chambers, 
• Ron Kikinis, 
• Arthur W. Toga

• Published: May 16, 2012
• https://doi.org/10.1371/journal.pone.0037454
• Reader Comments

White matter (WM) mapping of the human brain using neuroimaging techniques has gained considerable interest in the neuroscience community. Using diffusion weighted (DWI) and magnetic resonance imaging (MRI), WM fiber pathways between brain regions may be systematically assessed to make inferences concerning their role in normal brain function, influence on behavior, as well as concerning the consequences of network-level brain damage. In this paper, we investigate the detailed connectomics in a noted example of severe traumatic brain injury (TBI) which has proved important to and controversial in the history of neuroscience. We model the WM damage in the notable case of Phineas P. Gage, in whom a “tamping iron” was accidentally shot through his skull and brain, resulting in profound behavioral changes. The specific effects of this injury on Mr. Gage's WM connectivity have not previously been considered in detail. Using computed tomography (CT) image data of the Gage skull in conjunction with modern anatomical MRI and diffusion imaging data obtained in contemporary right handed male subjects (aged 25–36), we computationally simulate the passage of the iron through the skull on the basis of reported and observed skull fiducial landmarks and assess the extent of cortical gray matter (GM) and WM damage. Specifically, we find that while considerable damage was, indeed, localized to the left frontal cortex, the impact on measures of network connectedness between directly affected and other brain areas was profound, widespread, and a probable contributor to both the reported acute as well as long-term behavioral changes. Yet, while significantly affecting several likely network hubs, damage to Mr. Gage's WM network may not have been more severe than expected from that of a similarly sized “average” brain lesion. These results provide new insight into the remarkable brain injury experienced by this noteworthy patient.

Citation: Van Horn JD, Irimia A, Torgerson CM, Chambers MC, Kikinis R, Toga AW (2012) Mapping Connectivity Damage in the Case of Phineas Gage. PLoS ONE 7(5): e37454. https://doi.org/10.1371/journal.pone.0037454

Editor: Olaf Sporns, Indiana University, United States of America

Received: August 3, 2011; Accepted: April 23, 2012; Published: May 16, 2012

Copyright: © 2012 Van Horn et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by 2U54EB005149-06 “NAMIC: Traumatic Brain Injury - Driving Biological Project” to JVH, 1RC1MH088194 “Informatics Meta-Spaces for the Exploration of Human Neuroanatomy” to JVH, and P41RR013642 “Computational Anatomy and Multidimensional Modeling” to AWT. This work was performed as part of the Human Connectome Project (HCP; www.humanconnectomeproject.org ). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The mapping of human brain connectivity through the use of modern neuroimaging methods has enjoyed considerable interest, examination, and application in recent years [1] , [2] . Through the use of diffusion weighted (DWI) and magnetic resonance imaging (MRI), it is possible to systematically assess white matter (WM) fiber pathways between brain regions to measure fiber bundle properties, their influence on behavior and cognition, as well as the results of severe brain damage. The potential for using combined DWI/MRI methods to understand network-level alterations resulting from neurological insult is among their major research and clinical advantages.

In this paper, we investigate the detailed connectomics of a noted example of severe traumatic brain injury (TBI) which has proved important to and controversial in the history of neuroscience. Few cases in the history of the medical sciences have been so important, interpreted, and misconstrued, as the case of Phineas P. Gage [3] , in whom a “tamping iron” was accidentally shot through his skull and brain, resulting in profound behavioral changes, and which contributed to his death 151 years ago. On September 13th, 1848, the 25-year old Phineas P. Gage was employed as a railroad construction supervisor near Cavendish, Vermont to blast and remove rock in preparation for the laying of the Rutland and Burlington Railroad. Having drilled a pilot hole into the rock and filling it partially with gunpowder, he instructed an assistant to pour sand into the hole atop the powder. Averting his attention for a moment to speak with his men, he apparently assumed the sand had been added. He then commenced dropping the end of a 110 cm long, 3.2 cm diameter iron rod into the hole in order to “tamp” down its contents. The 13 lb. iron struck the interior wall of the hole causing a spark to ignite the powder which, in turn, launched the pointed iron rod upwards, through the left cheek of Mr. Gage just under the zygomatic arch, passing behind his left eyeball, piercing his cranial vault under the left basal forebrain, passing through his brain, and then exiting the top and front of his skull near the sagittal suture. A large amount of brain tissue was expelled from the opening and the rod was found later “smeared with blood and brains”, washed in a stream, and, eventually, returned to him. After receiving treatment and care from Dr. John Martyn Harlow over subsequent weeks, Mr. Gage was able to recover sufficiently from his physical injuries and return to his family in nearby New Hampshire. However, reports of profound personality changes indicate that he was unable to return to his previous job and caused co-workers to comment that he was “no longer Gage.” Following several years of taking manual labor jobs and travelling throughout New England and eventually to Valparaiso, Chile, always in the company of “his iron”, he was reunited with his family in San Francisco whereupon Mr. Gage died on May 21, 1860, nearly 12 years after his injury – presumably due to the onset of seizures evidently originating from damage resulting from the tamping rod incident. Several years later, Dr. Harlow, upon learning of Gage's death, asked Gage's sister's family to exhume his body to retrieve his skull and rod for presentation to the Massachusetts Historical Society and deposition with Harvard Medical School where, to this day, it remains on display in the Warren Anatomical Museum in the Francis A. Countway Library of Medicine at Harvard Medical School ( Fig. 1a ).

• PPT PowerPoint slide
• PNG larger image
• TIFF original image

a) The skull of Phineas Gage on display at the Warren Anatomical Museum at Harvard Medical School. b) CT image volumes were reconstructed, spatially aligned, and manual segmentation of the individual pieces of bone dislodged by the tamping iron (rod), top of the cranium, and mandible was performed. Surface meshes for each individual element of the skull were created. Based upon observations from previous examinations of the skull as well as upon the dimensions of the iron itself, fiducial constraint landmarks were digitally imposed and a set of possible rod trajectories were cast through the skull. This figure shows the set of possible rod trajectory centroids which satisfied each of the anatomical constraints. The trajectory nearest the mean trajectory was considered the true path of the rod and was used in all subsequent calculations. Additionally, voxels comprising the interior boundary and volume of the cranial vault were manually extracted and saved as a digital endocast of Mr. Gage's brain cavity. c) A rendering of the Gage skull with the best fit rod trajectory and example fiber pathways in the left hemisphere intersected by the rod. Graph theoretical metrics for assessing brain global network integration, segregation, and efficiency [92] were computed across each subject and averaged to measure the changes to topological, geometrical, and wiring cost properties. d) A view of the interior of the Gage skull showing the extent of fiber pathways intersected by the tamping iron in a sample subject ( i.e. one having minimal spatial deformation to the Gage skull). The intersection and density of WM fibers between all possible pairs of GM parcellations was recorded, as was average fiber length and average fractional anisotropy (FA) integrated over each fiber.

https://doi.org/10.1371/journal.pone.0037454.g001

The amount of damage to Mr. Gage's left frontal cortical grey matter (GM) with secondary damage to surrounding GM has been considered by several authors with reference to Gage's reported change in temperament, character, etc [4] , [5] , [6] ( Table 1 ). With the aid of medical imaging technology, two previous published articles have sought to illustrate the impact of the rod on Mr. Gage's skull and brain. Most famously, Damasio et al. [7] illustrated that the putative extent of damage to the left frontal cortex would be commensurate with the disinhibition, failures to plan, memory deficiencies, and other symptoms noted in patients having frontal lobe injury. Ratiu et al. [8] sought to illustrate the trajectory of the tamping iron, characterize the pattern of skull damage, and explain potential brain damage using a single, example subject. However, while many authors have focused on the gross damage done by the iron to Gage's frontal cortical GM, little consideration has been given to the degree of damage to and destruction of major connections between discretely affected regions and the rest of his brain.

https://doi.org/10.1371/journal.pone.0037454.t001

WM fasciculi link activity between cortical areas of the brain [9] , [10] , become systematically myelinated through brain maturation [11] , govern fundamental cognitive systems [12] , and may be disrupted in neurological [13] and psychiatric disease [14] . Penetrative TBI in cases of wartime [15] , industrial [16] , gunshot [17] , or domestic [18] injury often result in significant damage to brain connectivity, loss of function, and often death. Yet, in some instances, recovery from objects penetrating WM [19] have been reported with minimal sequelae [20] . Neuroimaging studies of WM tracts in TBI have revealed not only significant acute damage to fiber pathways but also that measures of fiber integrity can show partial fiber recovery over time [21] , presumably due to cortical plasticity [22] in non-penetrative cases.

Given recent interest in the atlasing of the human WM connectome (e.g. http://www.humanconnectomeproject.org ), a detailed consideration of the putative damage to Mr. Gage's connectomics and implications for changes in behavior is provocative and compelling. Nerve damage is superficially evident through reports of eventual loss of sight in Gage's left eye, left eyelid ptosis [23] , and recognition of potential WM damage by other investigators [7] . Further examination of the extent of Gage's WM damage and of its effects on network topology and regional connectedness can offer additional context into putative behavioral changes. Due to the absence of original brain tissue and to the lack of a recorded autopsy from this case, one can only estimate the extent of damage from bony structures and can never be confident concerning which precise brain tissues were impacted. However, brain tissue in situ from a representative population can be considered and it can be assumed that Mr. Gage's anatomy would have been similar. In this examination, we obtained the original high-resolution CT data of the Gage skull used by Ratiu et al. , and computationally estimated the best-fit rod trajectory through the skull. Via multimodal analysis of T1-weighted anatomical MRI and DWI in N = 110 normal, right-handed males, aged 25–36, we quantify the extent of acute regional cortical loss and examine in detail the expected degree of damage to Mr. Gage's WM pathways.

Computationally projecting a model of the tamping iron through the T1 MRI anatomical volumes warped to the Gage skull geometry ( Table 2 ; Fig. 1b–c ; see also Methods ) in light of previously reported anatomical constraints ( Table 3 ) and healthy brain morphometry and connectivity ( Fig. 2 ), the average percentage of total cortical GM volume intersected was 3.97±0.29% (mean±SD), where the cortical regions most affected by the rod (>25% of their regional volumes) included (mean±SD): the left orbital sulcus (OrS; 90.86±6.97%), the left middle frontal sulcus (MFS; 80.33±10.01), the horizontal ramus of the anterior segment of the lateral sulcus (ALSHorp; 71.03±22.08%), the anterior segment of the circular sulcus of the insula (ACirInS; 61.81±18.14%), the orbital gyrus (OrG; 39.45±6.17%), the lateral orbital sulcus (LOrS; 37.96±20.24%), the superior frontal sulcus (SupFS; 36.29±12.16%), and the orbital part of the inferior frontal gyrus (InfFGOrp; 28.22±19.60%). While extensive damage occurred to left frontal, left temporal polar, and insular cortex, the best fit rod trajectory did not result in the iron crossing the midline as has been suggested by some authors (see Methods ). As a result, no direct damage appeared to occur in right frontal cortices as evident from our representative sample cohort. A complete list of all cortical areas experiencing damage is listed in Table 4 .

The outermost ring shows the various brain regions arranged by lobe (fr – frontal; ins – insula; lim – limbic; tem – temporal; par – parietal; occ- occipital; nc – non-cortical; bs – brain stem; CeB - cerebellum) and further ordered anterior-to-posterior based upon the centers-of-mass of these regions in the published Destrieux atlas [72] (see also Table 6 for complete region names, abbreviations, and FreeSurfer IDs, and Table 7 for the abbreviation construction scheme). The left half of the connectogram figure represents the left-hemisphere of the brain, whereas the right half represents the right hemisphere with the exception of the brain stem, which occurs at the bottom, 6 o'clock position of the graph. The lobar abbreviation scheme is given in the text. The color map of each region is lobe-specific and maps to the color of each regional parcellation as shown in Fig. S2 . The set of five rings (from the outside inward) reflect average i) regional volume, ii) cortical thickness, iii) surface area, and iv) cortical curvature of each parcellated cortical region. For non-cortical regions, only average regional volume is shown. Finally, the inner-most ring displays the relative degree of connectivity of that region with respect to WM fibers found to emanate from this region, providing a measure of how connected that region is with all other regions in the parcellation scheme. The links represent the computed degrees of connectivity between segmented brain regions. Links shaded in blue represent DTI tractography pathways in the lower third of the distribution of fractional anisotropy, green lines the middle third, and red lines the top third. Circular “color bars” at the bottom of the figure describe the numeric scale for each regional geometric measurement and its associated color on that anatomical metric ring of the connectogram.

https://doi.org/10.1371/journal.pone.0037454.g002

https://doi.org/10.1371/journal.pone.0037454.t002

https://doi.org/10.1371/journal.pone.0037454.t003

https://doi.org/10.1371/journal.pone.0037454.t004

The amount of total WM volume lost due to the tamping iron was 10.72±5.46% (mean±SD). Examination of lesioned connectivity matrices indicated that fiber bundles from nearly the entire extent of the left frontal cortex were impacted by the presence of the tamping iron (e.g. Fig. 1d ), which in turn affected most of that hemisphere as well as contralateral regions ( Fig. 3 ). The effect of this lesion on network properties was assessed 1) with respect to the healthy intact network, generally, as well as 2) in contrast to the average effects of similarly-sized lesions simulated elsewhere in the cortex, as related to local GM loss as well as distributed loss of connectivity ( Fig. 4a–c ). Metrics representative of three specific global network attributes were examined: characteristic path length (λ, measuring network integration), mean local efficiency (e, segregation), and small worldness (S) ( Table 5 ).

The lines in this connectogram graphic represent the connections between brain regions that were lost or damaged by the passage of the tamping iron. Fiber pathway damage extended beyond the left frontal cortex to regions of the left temporal, partial, and occipital cortices as well as to basal ganglia, brain stem, and cerebellum. Inter-hemispheric connections of the frontal and limbic lobes as well as basal ganglia were also affected. Connections in grayscale indicate those pathways that were completely lost in the presence of the tamping iron, while those in shades of tan indicate those partially severed. Pathway transparency indicates the relative density of the affected pathway. In contrast to the morphometric measurements depicted in Fig. 2 , the inner four rings of the connectogram here indicate (from the outside inward) the regional network metrics of betweenness centrality, regional eccentricity, local efficiency, clustering coefficient, and the percent of GM loss, respectively, in the presence of the tamping iron, in each instance averaged over the N = 110 subjects.

https://doi.org/10.1371/journal.pone.0037454.g003

WM fiber pathways intersected by the rod were pooled across all N = 110 subjects and examined for a) the relative lengths (w ij ) of affected pathways and b) the relative percentages of lost fiber density (g ij ); c) the bivariate distribution of g ij versus w ij indicating that local fiber pathways were affected, e.g. relatively short pathways proximal to the injury site, as well as damaging dense, longer-range fiber pathways, e.g. innervating regions some distance from the tamping iron injury (see “ Calculation of Pathology Effects upon GM/WM Volumetrics ” for further details).

https://doi.org/10.1371/journal.pone.0037454.g004

https://doi.org/10.1371/journal.pone.0037454.t005

Tables 6 and 7 provide details on the regional coding used for brain parcellation which were subjected to estimation of the effects of the tamping iron, lesion simulation modeling, and which encode the text on the outer-most rings of Figs. 2 and 3 . Differences in measures of network connectivity due to the rod's passage were apparent in terms of network integration, segregation, but not small worldness as compared to the unlesioned, healthy network. Specifically, when removing those cortical areas and fiber pathways intersected by the iron, characteristic path length was found to be significantly decreased in Gage compared to the intact network ( p ≤0.0001), mean local efficiency was decreased ( p ≤0.0001), while small worldness showed no statistical difference ( p ≤0.9467, ns). Regionally-specific network theoretical metrics in the affected regions and those to which they connect were also affected (see Fig. 5A ). This suggests that, not surprisingly, with significant loss of WM connectivity between left frontal regions and the rest of the brain, the surviving network of brain was likely to have been heavily impaired and its functions considerably compromised.

A) Cortical maps of regional graph theoretical properties. Regions affected by the passage of the tamping iron include those having relatively high betweenness centrality and clustering coefficients but relatively low mean local efficiency and eccentricity. B) A cortical surface schematic of the relative effects of systematic lesions of similar WM/GM attributes over the cortex for both network integration (i) and segregation (ii). For each mapping, colors represent the Z-score difference between systematic lesions of that area relative the average change in integration taken across all simulated lesions. C) Cortical maps of the differences/similarity between the effects on integration and segregation observed from the tamping iron lesion with that of each simulated lesion. Here black is most similar (e.g. the observed lesion is most similar to itself) whereas white is least similar to (e.g. most different from) the tamping iron's effects on these measures of network architecture.

https://doi.org/10.1371/journal.pone.0037454.g005

https://doi.org/10.1371/journal.pone.0037454.t006

https://doi.org/10.1371/journal.pone.0037454.t007

To further provide a baseline for comparison of the tamping iron lesion against similarly-sized lesions located elsewhere in the cortex, we conducted a systematic random simulation of 500 similarly-sized lesions across our N = 110 healthy subject cohort. The network containing the lesion due to the tamping iron was systematically compared against the distributions of the above mentioned metrics from the simulated lesion set. When paired t -statistics were computed to determine whether tamping iron lesion differed significantly from the standpoint of network metric values, standardized with respect to the intact network, as compared to other brain lesions of the same size, the characteristic path length (integration), mean local efficiency (segregation), and small worldness, while significantly different from that of the intact networks, were not found to be more severe than the average network properties of average similarly sized GM/WM lesion. These results are summarized in Fig. 5B and 5C . These indicate that alterations to network integration resulting from the tamping iron lesion resulted in greater average path length than that of the intact network but which was less than the average effects of other equally sized lesions. Likewise, segregation, as measured using mean local efficiency, was reduced compared to the intact network, but greater than the average effects of the simulated lesions. These results suggest that Mr. Gage's lesion, while severe and certain to have affected WM connectivity in his left cerebral hemisphere and throughout his brain, could have been considerably more severe had the tamping iron pierced other areas of his brain.

The case of Phineas Gage is among the most famous and infamous in the history of brain science. The interpretations of his incredible injury and attempts to characterize it have been ongoing since soon after it occurred. Our consideration sought to provide a modern connectomic understanding of Mr. Gage's injury and put it into context as involving brain WM in addition to the GM damage discussed by other authors. While we, too, are constrained by the relics left from Mr. Gage's life and what evidence can be gleaned from them, work detailed in this article differs considerably from previous examinations of this case and topic in several key areas: we 1) precisely model the trajectory of the tamping iron through high resolution computed tomographic data of Mr. Gage's skull - a rare imaging data set that, until now, had been lost to science for over a decade; 2) geometrically fit N = 110 age, gender, and handedness matched modern subject MRI brain volumes into the Gage cranial vault to assess average cortical metrics and their degree of variability; 3) in so doing, illustrate that while ∼4% of the cortex was intersected by the rod's passage, ∼11% of total white matter was also damaged, and provide estimates of the degree of damage experienced under a well-established brain parcellation scheme; 4) map high angular resolution diffusion neuroimaging tractography into the same space to measure damage to the pair-wise connections between atlas-defined cortical regions; and 5) compare the graph theoretical properties of the observed lesion against those expected from theoretically similar lesions systematically located throughout the brain. In what follows, we comment on our approach and findings.

Trajectory of the Tamping Iron

Various descriptions of the trajectory of the tamping iron through the Mr. Gage's skull have been given, which has understandably led to differing opinions about which parts of his brain were subjected to damage. Harlow, the physician responsible for Gage's initial treatment, documented that only the left hemisphere had been affected while the right remained unaffected [24] . In contrast, Bigelow maintained that some right-sided damage must have occurred. Dupuy [25] agreed with the left sidedness of the trajectory but placed it more posterior, claiming that motor and language areas had been destroyed – supporting the anti-localizationist arguments popular of the era. Ferrier [26] , illustrating that the motor and language areas had been spared, concluded that damage was limited only to the left hemisphere – a conclusion later echoed by Cobb [27] . In their measurements, Damasio et al. estimated the damage to be more frontal and right sided, whereas Ratiu and colleagues concluded that damage was limited to the left frontal lobe and that it did not cross the midline. Central to these differences in interpretation is likely to be how mandible position has been considered. To satisfy the observed anatomical constraints with the mouth closed would result in a greater right-sided inclination of the rod. Yet, as Harlow originally noted, Gage was in the act of speaking to his men at the moment of the injury and, thus, his mouth was likely open. We observe that with the jaw opened, the best-fit rod trajectory satisfying all constraints does not intersect or cross the superior sagittal sulcus and the injury is specific to the left frontal lobe. Thus, our conclusions are congruent with those of Harlow, Ferrier, as well as with those of Ratiu and Talos and, given the detailed computational approach taken, seem to provide the most likely reconstruction of the acute damage caused by the tamping iron.

Alterations of Network Connectivity Due to the Tamping Iron

The loss of ∼11% total WM volume in the left frontal lobe suggests that the iron's effects on Mr. Gage's brain extended well beyond the loss of left frontal GM alone. Overall differences in metrics of network integration as well as segregation were observed relative to intact connectivity, suggesting widespread disruption of networks involving damage to the left frontal and temporal pathways. Alterations indicate major changes to global network topology which affected network-wide efficiency. In the healthy cohort examined here, the region-to-region WM connectedness when in the presence of the rod was found to be associated with several important fiber bundles. Specifically, connectivity was affected between the frontal lobes and the basal ganglia, the insula, limbic, and other major lobes of the left hemisphere, in addition to right frontal, insular, and limbic areas. This severed portions of the uncinate fasciculus (UF) - connecting parts of the limbic system such as the hippocampus and amygdala in the temporal lobe with frontal regions such as the orbito-frontal cortex. The cingulum bundle - the collection of WM fibers projecting from the cingulate gyrus to the entorhinal cortex, allowing for communication between components of the limbic system – was also damaged. Additionally, the superior longitudinal fasciculus (SLF) was impacted – the long bi-directional bundles of neurons connecting the rostral and caudal aspects of the cerebrum in which each association fiber bundle is lateral to the centrum ovale linking the frontal, occipital, parietal, and temporal lobes. Fibers here pass from the frontal lobe through the operculum to the posterior end of the lateral sulcus, where numerous processes radiate into the occipital lobe while others turn downward and forward around the putamen and project to anterior portions of the temporal lobe. The occipito-temporal projection [28] in humans connects the temporal lobe and occipital lobe, running along the lateral walls of the inferior and posterior cornua of the lateral ventricle. The connectivity of the orbital cortex with temporal lobe regions via the UF which is among the last to complete myelination in development [29] , has been shown to be particularly affected in patients with mental illness [30] , and to be related to cognitive deficits in TBI [31] , [32] . WM fascicular damage in these instances was likely an important factor in Gage's reported post-injury symptomatology as well as in his reported and putative behavioral issues.

The obtained results suggest that GM damage had wider reaching influence than previously described and compromised several aspects of Gage's network of WM connectivity. Regions whose connectivity within and between cerebral hemispheres were affected included: the left frontal lobe (the transverse fronto-polar gyrus, fronto-marginal gyrus, middle frontal gyrus, lateral orbital sulcus, orbital sulcus, oribital part of the inferior frontal gyrus, triangular part of the inferior frontal gyrus, inferior frontal sulcus, medial orbital sulcus, orbital gyri, superior frontal gyrus, and opercular part of the inferior frontal gyrus); left insular cortex (horizontal and vertical ramus of the anterior segment of the lateral sulcus/fissure, the anterior/inferior/superior segments of the circular sulcus of the insula, short insular gyri, and long insular gyrus and central insular sulcus); and the left temporal lobe (the temporal pole and polar plane of the superior temporal gyrus) ( Fig. 3 ). The marginal and bivariate probability distributions of average brain-normalized WM fiber bundle length (w ij ) and proportion of GM density lost (g ij ; Fig. 4a–c ) unsurprisingly indicated a considerable number of relatively short connections being affected locally by the presence of the rod while, additionally, a considerable number of longer fiber bundles connecting relatively large regions of cortex were also impacted.

Alterations to these connections contribute to the significant reductions in characteristic path length and mean local efficiency of the remaining network after removal of the affected fibers. That no significant difference was observed concerning the small worldness of the tamping iron network as compared to the intact network suggest that a lesion of this size and scope, while severe, may not had appreciable effects on the degree of clustering of unaffected nodes Mr. Gage's brain relative to randomly degree-equivalent versions of that network. On the other hand, the average simulated lesion did show a significant reduction in small worldness, indicating that regions other than those affected may have more influence over the degree of measured network clustering. Thus, Mr. Gage's unaffected network may have still maintained its small world architecture of nodal clustering and presumed functional integrity, despite loss of major frontal and temporal lobe participation in the system resulting in deficits to measures network integration and segregation.

Several previous articles have precisely investigated the direct effects of node deletions of various size on network connectivity and architecture [33] , [34] , [35] . In particular, the paper by Alstott et al. provides a detailed examination of simulated lesion effects on brain networks both in terms of lesion location and extent. In their study of structural and functional connectivity data from N = 5 healthy subjects, they found that lesions to midline areas resulted in more profound effects on various network metrics than do more lateral brain regions. As might be expected, the magnitude of change was dependent upon the number of nodes removed from the network and the manner in which they were removed. Their observations indicate that networks may be insensitive to lesions involving random node removal or where node removal was based only upon a node's degree of connectedness. However, network lesioning based upon the targeted removal of nodes having high betweenness centrality - a measurement of the number of shortest paths from all vertices to all others which pass through a given node - resulted in greater network vulnerability as evident from significant reductions in global efficiency in contrast to random lesioning. This result is particularly compelling in regards to assessing the robustness of cortical architecture in the face of brain damage to major network hubs localized proximal to the cortical midline.

There is little doubt that a tamping iron injury to central nodes of the frontal lobe would have severely impacted Gage's brain connectivity. Fig. 3 shows the extent of white matter damage and the effects on several measures of network connectivity, including regional betweenness centrality, local efficiency, clustering coefficient, and eccentricity. Note that this illustration differs from Fig. 2 in that the inner-most rings are now colored according to the respective average nodal connectivity metrics in the presence of cortical loss incurred from the tamping iron. Additionally, Fig. 5A illustrates the spatial distribution of these regional connectivity metrics over the cortex when pooled across the intact healthy networks from our sample (sub-cortex not shown). We note that, as observed by Alstott et al., areas of relatively high betweenness centrality tended to be located along the frontal midline. Other metrics show similar regional concentrations ( Fig. 5Aii–iv ). However, while intact frontal areas of both hemispheres show high betweenness centrality ( Fig. 5Ai ), the regions of tamping iron damage encompassed many other regions as well having relatively less betweenness centrality, e.g. TrFPoG/S, RG, SbCaG, TPo. Removal of these areas, as illustrated by the various metric rings in the left frontal segment of the connectogram in Fig. 3 , has wide ranging effects on the regionally-specific network metrics in unaffected brain regions.

It is evident that removal of these areas produce significant effects on global metrics of network segregation and integration. However, from systematic lesion simulation using a similar extent of GM/WM involvement, the effects on Mr. Gage's network integration and segregation were not found to be more severe that that observed from the “average” lesion. Clearly, a larger lesion would have affected a greater number of network nodes including various hubs resulting in further deleterious effects on network integration and segregation. Moreover, a different lesion altogether would have possibly resulted in more outwardly obvious sensorimotor deficits. Located in occipital cortex, for instance, the lesion might have resulted in sensory-specific changes in connectivity (e.g. blindness), or one involving more of the sub-cortex and brain stem could have been more clinically serious and resulted in death. Nevertheless, the observed damage illustrates that severe network insult affecting the majority of left hemisphere connectivity as well as right hemispheric inter-connections, was experienced. Such damage can be expected to have had its influence over the normal functioning of many regions non-local to the injury and their subsequent connectivity as well.

Therefore, in light of these observations, it would be safe to conclude that 1) Mr. Gage's injury very likely destroyed portions of the central hub structure in left frontal midline structures as well as temporal pole and limbic structures which have extensive connectivity throughout the left hemisphere as well as inter-hemispherically, 2) that the tamping iron's passage did not specifically remove only the most central network hubs but a host of regions having a range of network properties, and 3) that such damage to important network hubs connection to other brain regions having secondary levels of centrality, clustering, etc. are likely to have combined to give rise to the behavioral and cognitive symptomatology originally reported by Harlow. Knowledge of Gage's affected connectivity help provide clarity and context for symptomatologies subsequently only inferred by others.

Implications for Gage's Reported Behavioral Changes

Traumatic brain injury of the frontal cortices is often associated with profound behavioral alterations, changes mood [36] , working memory [37] and planning deficits [38] , [39] , social functioning [40] , among other cognitive symptoms [41] , [42] , [43] , [44] . Alterations to functional connectivity have also been reported [45] , [46] which, in addition to cortical damage, likely related to accompanying diffuse axonal injury [47] , [48] . It is also worth noting neurodegenerative diseases, such as the leukodystrophies [49] , Alzheimer's Disease (AD) [50] , [51] , and early-stage frontotemporal dementia (FTD) [52] , also have effects on brain networks involving connectivity of the frontal lobe. Altered structural connectivity in these disorders illustrates changes in large-scale brain network organization deviating from healthy network organization [53] , with possible effects on resting state connectivity [54] . Disruptions of WM connectivity are also known to underlie elements of psychiatric illness [55] , [56] , [57] which are associated with behavioral alterations not dissimilar to those reported in Mr. Gage.

In particular, network damage, predominantly of the left basal forebrain and of its connections throughout the left as well as into right frontal cortices, was particularly extensive. Processing of emotion stimuli have been associated with connectivity of the frontal cortex and amygdala, in particular involving the connectivity of the uncinate fasciculi [58] . Thus, in addition to disinhibition symptoms considered by Damasio et al., with evidence of potentially greater degree of WM rather than cortical injury, there is also similarity between Mr. Gage's behavioral changes and network alterations observed in FTD and related WM degenerative syndromes. This suggests that network topological changes may have been the source of Mr. Gage having not only executive function deficits but also problems resulting from damage to connections associated with the encoding of episodic memory as well as the processing of emotion – consistent with reports on changes in his personality.

Historical Implications of Gage's WM Damage

While observations of severe network damage and their resulting affects may not be surprising given that which has been documented of Mr. Gage's accident and behavioral changes, one can only speculate upon the possible contribution to Gage's survival, recovery, and the uniqueness of changes to his WM networks. Macmillan [3] has noted that many reports on Gage's behavioral changes are anecdotal, largely in error, and that what we formally know of Mr. Gage's post-accident life comes largely from the follow-up report of Harlow [23] according to which Gage, despite the description of him having some early difficulties, appeared to adjust moderately well for someone experiencing such a profound injury. Indeed, the recent discovery of daguerreotype portraits of Mr. Gage show a “handsome…well dressed and confident, even proud” man [59] in the context of 19 th century portraiture. That he was any form of vagrant following his injury is belied by these remarkable images. While certainly neuroanatomically profound, the changes to his cognitive capacities were much more subtle upon his full recovery than may have been otherwise described. In spite of recovering from severe brain trauma, his mental state appears to have eventually stabilized sufficiently for him to travel throughout New England, take on several (some might say menial) forms of employment, travel through South America for several years, and to return to his family in the Western US, before succumbing to epilepsy which was presumably related to the injuries directly affecting his WM connectivity. That his network damage, though extensive, was not apparently more severe than an “average” brain lesion would incur may help to explain his ability to have sufficiently recovered in spite of the residual behavioral changes reported by Harlow.

Limitations of our Study

We have worked to provide a detailed, accurate, and comprehensive picture of the extent of damage from this famous brain injury patient and its effect on network connectivity. While the approach used here to model the tamping iron's trajectory is precise and the computation of average volume lost across our population of subjects is reflective of the acute level of damage, we acknowledge that there was likely more damage than that caused by its presence alone. The iron likely propelled unrecovered bone fragments through the brain. The resulting hemorrhage from the wound was also considerable. Subsequent infection and a large abscess took further toll. Consequently, more GM and WM tissue may have been lost than estimated here. Like Damasio et al. and Ratiu et al. , we make the assumption that Gage's brain and its position within the skull can be estimated from the structure of the skull itself, and that its sub-regions, WM, and connective anatomy can be localized through population averaging. Such a supposition may have its limitations and could be open to debate. Nevertheless, ours represents the best current estimation as to the extent of brain damage likely to have occurred at the level of both cortex and WM fiber pathways. We also have no way of assessing the biochemical cascade of changes to biomarker proteins measureable post-injury in modern TBI patients which may also have influenced the trajectory of Mr. Gage's recovery.

Another potential criticism is that we compare the loss of GM, WM, and connectivity in Mr. Gage by computationally casting the tamping iron through the WM fibers of healthy age- and gender-matched subjects and measuring the resulting changes in network topology. We also systematically lesion the brains of our healthy cohort to derive “average” network metrics and compare the observed values with respect to them – an approach that has been recommended elsewhere [35] . This technique is helpful for creating a representative expectation of inter-regional connectivity against which to compare observed or hypothetical lesions. However, some might consider this approach to be misguided in this instance due to the fact that Mr. Gage's brain was damaged in such a way that he survived the injury whereas a host of other lesions resulting from penetrative missile wounds would likely have resulted in death. Indeed, as noted originally by Harlow, the trajectory of the 110 cm long, 3.2 cm thick, 13 lb. tamping iron was likely along the only path that it could have taken without killing Mr. Gage. Thus, any distribution of lesioned topological values might not provide a useful foundation for comparison because the majority of these penetrative lesions would, in reality, be fatal. We recognize these concerns and the practical implications for subject death which would also be a caveat of other network theoretical applications of targeted or random network lesioning. Indeed, such considerations are something to be taken into account generally in such investigations. Nevertheless, our simulations provide supporting evidence for the approximate neurological impact of the tamping iron on network architecture and form a useful basis for comparison beyond utilizing the intact connectivity of our normal sample in assessing WM connectivity damage. So, while this might be viewed as a limitation of our study, especially given the absence of the actual brain for direct inspection, the approach taken provides an appropriate and detailed assessment of the probable extent of network topological change. All the same, we look forward to further work by graph theoreticians to develop novel approaches for assessing the effects of lesioned brain networks.

Conclusions

In as much as earlier examinations have focused exclusively on GM damage, the study of Phineas Gage's accident is also a study of the recovery from severe WM insult. Extensive loss of WM connectivity occurred intra- as well as inter-hemispherically, involving direct damage limited to the left cerebral hemisphere. Such damage is consistent with modern frontal lobe TBI patients involving diffuse axonal injury while also being analogous to some forms of degenerative WM disease known to result in profound behavioral change. Not surprisingly, structural alterations to network connectivity suggest major effects on Mr. Gage's overall network efficiency. Connections lost between left-frontal, left-temporal, right-frontal cortices as well as left limbic structures likely had considerable impact on executive as well as emotional functions. Consideration of WM damage and connectivity loss is, therefore, an essential consideration when interpreting and discussing this famous case study and its role in the history of neuroscience. While, finally, the quantification of connectomic change might well provide insights regarding the extent of damage and potential for clinical outcome in modern day brain trauma patients.

Ethics Statement

No new neuroimaging data was obtained in carrying out this study. All MRI data were drawn from the LONI Integrated Data Archive (IDA; http://ida.loni.ucla.edu ) from large-scale projects in which subjects provided their informed written consent to project investigators in line with the Declaration of Helsinki, U.S. 45 CFR 46, and approval by local ethics committees at their respective universities and research centers. Research neuroimaging data sets deposited with the LONI IDA and made available to the public are fully anonymized with respect to all identifying labels and linked meta-data for the purposes of data sharing, re-use, and re-purposing. IDA curators do not maintain linked coding or keys to subject identity. Therefore, in accordance with the U.S. Health Insurance Portability and Accountability Act (HIPAA; http://www.hhs.gov/ocr/privacy ), our study does not involve human subjects' materials.

Medical Imaging of the Gage Skull

Medical imaging technology has been applied to the Gage skull on three known occasions to model the trajectory of the tamping iron, infer extent of GM damage, and theorize about the changes in personality which a patient with such an injury might have incurred. In an influential study, Damasio and coworkers [7] used 2D X-rays to obtain the dimensions of the skull itself and to compute the trajectory of the iron bar through the regions of frontal cortex based on independently obtained CT data from a normal subject. Prior to this, CT scanning of the skull had been obtained by Tyler and Tyler in 1982 for presentation and discussion at a neurological scientific meeting. The location of the raw CT data files from this imaging session is unknown but the data were last reproduced in An Odd Kind of Fame (Appendix E), though they were not part of any other scientific publication of which we are aware. The most recent occurrence of scanning on record was performed on June 12 th , 2001 through the Surgical Planning Laboratory (SPL) at Brigham and Women's Hospital, Harvard Medical School. A series of two high-resolution CT image series were obtained of the skull: one covering the portion of the jaw up to approximately the bridge of the nose, and another covering the cranial vault (see details below). These data were used by Ratiu et al. [8] , [60] to digitally reconstruct and animate the passage of the tamping iron through the skull. An additional CT image of the Gage life-mask, a plaster likeness presumed to have been commissioned by Dr. Bigelow during one of Gage's visits to Harvard Medical School, was also obtained and used to create a surface model of Mr. Gage's face, scalp, and neck. New CT or other medical imaging of the skull specimen is unlikely to be performed in the future due to the age and fragile state of the specimen.

Documented Extent of Neurological Damage

In the book An Odd Kind of Fame (2000, pg 85), Macmillan conveniently summarizes the reports from various anatomists on the damage to Gage's brain. We reproduce these summaries here and also add the findings of Ratiu et al. [8] which appeared after the publication of An Odd Kind of Fame .

Skull CT Data Processing

Due to a variety of circumstances, the raw and processed digital imaging data from the 2001 CT imaging session at Brigham and Women's Hospital were improperly archived and effectively lost to science. However, these image volumes were subsequently recovered by the authors and represent the highest quality data/resolution available (0.5 mm slice thickness) for modeling the skull of this noted patient and for use in the modeling of affected anatomy and connectivity. The scan data were originally obtained with the superior, cut portion of the calvarium and the mandible in the correct anatomical position on a Siemens Somatom CAT scanner (Siemens AG, Erlangen, Germany), in the Department of Radiology, Brigham and Women's Hospital (Boston, MA) [8] . These data were converted from ECAT format to the NIFTI file format ( http://nifti.nimh.nih.gov ) using the program “mri_convert” – part of the FreeSurfer neuroimaging data analysis software package (surfer.nmr.mgh.harvard.edu/fswiki/mri_convert). The CT images were systematically segmented and masked by hand using MRICron ( http://www.cabiatl.com/mricro/mricron/index.html ) and seg3D ( http://www.sci.utah.edu/cibc/software/42-seg3d.html ) to isolate the skull cap (the portion of the skull created by its being cut with a saw upon deposition at the Warren Museum by Dr. Harlow), each piece of remaining/healed bone fragments, the left frontal/temporal portion of the skull along the readily evident fracture lines, and the lower jaw, and separate 3-D surface mesh models were generated for each segment using 3D Slicer ( http://www.slicer.org ). An additional binary image volume was created by hand-filling the space of the cranium that contained Gage's brain. This volume represents a digital version of the standard endocast often used in the analysis of paleontological specimens [61] , [62] , [63] . Use of the Gage skull and life mask CT data is courtesy of the SPL and the Warren Anatomical Museum at Harvard Medical School.

The LONI Pipeline Workflow Environment

For all major image processing operations (e.g. bias field correction, skull stripping, image alignment, etc.) we employed the LONI Pipeline Workflow Environment ( http://pipeline.loni.ucla.edu ; Fig. S1 ). This program is a graphical environment for construction, validation, and execution of advanced neuroimaging data analysis protocols. It enables automated data format conversion, leverages Grid computer systems, facilitates data provenance, and provides a significant library of computational tools [64] , [65] , [66] .

For instance, employing LONI Pipeline, we used the Brainsfit software package ( http://www.nitrc.org/projects/multimodereg/ ) to register the T1 anatomical MRI volumes to the endocast template. Diffusion gradient image data were processed in native subject space using Diffusion Toolkit ( http://trackvis.org ) to reconstruct the fiber tracts. Data processing workflows to compute inter-regional connectivity matrices were constructed using purpose-built software. Fig. S2 illustrates an example connectivity matrix displayed using Matlab (Mathworks, Natick, MA, USA).

Measurements of the Skull

Consistent with Damasio et al. , the physical dimensions of the Gage skull were measured as follows in Table 2 using the Slicer software program. Additionally, the following landmarks were identified on the Gage skull: Entrance of the Left Auditory Canal: (49.56, 219.46, −807.75 mm); Entrance of the Right Auditory Canal: (175.04, 212.26, −802.85 mm); and the Middle of Crease Between Frontal Bone Plate and Nasal Bone: (117.04, 301.73, −800.72 mm). Given these landmarks, all the other points can be accurately positioned.

Measurements of the Tamping Iron

One of our team (MCC) visited the Warren Anatomical Museum and, working with lead curator Dominic Hall, obtained the following measurements of the iron using a SPI Digimax caliper (Model: 30440-2): 110 cm in length, 9.5 cm circumference, and 2.88 cm diameter at tail. The rear taper is approximately 19 cm long, the maximum diameter (between the rear and tip taper) is 10.5 cm circumference (3.2 cm diameter), the taper beginning at the tip is 27 cm long, and the diameter at the rod's tip is 72 mm.

The Trajectory of the Tamping Iron

The trajectory of the tamping iron through Mr. Gage's skull and brain has been the subject of much debate and several attempts have been made to infer the relationship between putative damage on the one hand and the lore surrounding Gage's personality and behavioral changes resulting from his accident on the other. Bigelow [67] first attempted to formally model the trajectory of the rod by drilling a hole through another “common” skull (pg. 21), and noted that “a considerable portion of the brain must have been carried away; that while a portion of its lateral substances may have remained intact, the whole central part of the anterior lobe, and the front of the sphenoidal or middle lobe must have been lacerated and destroyed”. Importantly, Damasio [7] and coworkers provided a detailed analysis of the rod trajectory through the skull attempting to identify which brain regions were impacted by the flight of the iron and what effect this impact had on the patient's post-injury behavior. While this study has been well cited, their methodology for determining the rod trajectory has been subsequently questioned [3] .

Ratiu et al. [60] constrained their modeling of the rod trajectory by noting bony injuries to the skull, and by more closely aligning the rod with the clinical information provided by both Harlow and Bigelow. Ratiu et al. inserted the brain of a single normal subject into Mr. Gage's cranial cavity to examine which structures might have been affected. Their reconstruction shows that the path of the iron passed left of the superior sagittal sinus (their Fig. 4b,d ). This is corroborated by the fact that damage to the superior sagittal sinus would have almost certainly caused air embolism and/or significant blood loss, resulting in Mr. Gage's death. In addition, their reconstruction shows, in their normal subject's brain, that the iron's trajectory was also anterior to the cingulate gyrus and to the left lateral ventricle (their Fig. 4 e,f ). No rhinoliquorhea or other indication of post-traumatic CSF fistula was reported, nor that Gage developed ventriculitis, a condition which very likely would have been lethal - especially in the 1840's before the use of antibiotics in common medical practice. However, there is little way of being empirically precise with respect to location of major structures when employing only a single, example subject to represent Mr. Gage's unknown neuroanatomy.

To address this issue, we fit the T1 anatomical and diffusion images from the N = 110 normal, right handed subjects, aged 25–36 into the space of Phineas Gage's cranial vault to map the probability to regional injury and the effects of the tamping rod on WM fiber connectivity. The process of morphing data into the Gage skull is described in the following sections.

Determining the Trajectory of the Tamping Iron

Using the measurements of the original tamping iron [3] , [8] , [24] , [67] , on display at the Warren Museum, a 3-D model of the tamping iron was generated using Matlab and stored as an VTK surface ( http://www.vtk.org ) for visualization using 3D Slicer and for processing using the segmented brain regions and fiber tracts morphed into the space of the Gage 3D cranial endocast volume model.

To constrain the trajectory of the rod through the Gage skull, we examined the work of previous authors to identify noteworthy statements on the condition of the skull, particular patterns of breakage, chips in the bone, and other prominent features that could be used as landmarks to restrict the possible paths which the rod might have taken ( Fig. S3 ). For instance, the left maxillary molar is missing and osteological analysis by the Warren Museum states that it was lost ante-mortem (Object File WAM 00949, Warren Anatomical Museum, Francis A. Countway Library of Medicine). While Harlow and/or Bigelow do not specifically mention the loss of this tooth, it is likely that the rod made contact with it after passing through Gage's cheek, and was either dislodged completely or knocked loose and lost sometime during his recovery. Additionally, for the zygomatic arch the Warren Museum records (also WAM 00949) indicate “Maxilla: ante-mortem sharp force trauma remodeling” but are not more specific about the potential for complete breakage of the zygomatic process which was suspected by Ratiu et al. Still, it can be assumed that some contact was made between the iron and the interior portions of the arch. A collection of previously reported observations contributing to the set of applied constraints are noted in Table 3 .

In particular, we concur with Ratiu et al. that Mr. Gage had his jaw open at the moment of the accident. Harlow reports Gage looking over his right shoulder and saying something to his crew at critical moment of the blast. In the casting of possible rod trajectories, the most likely position of the jaw was determined to be −15° in pitch (downward) and 5° in yaw (to the right) relative to the closed position of the jaw. This position allowed the unhindered passage of 1.303×10 3 out of 1×10 9 viable rod trajectories inclusive through the skull. With this jaw position, in contrast to the suspicion of Bigelow, we noted no contact between the rod and that of Mr. Gage's coronoid process. Jaw rotations at greater pitch angles were inconsequential to our results. Therefore, these values represent the minimal angular jaw deflections needed to allow the maximal number of rod passage scenarios without jaw intersection. Additionally, these values are typical for the acts of speaking and mastication in which the maximum typical jaw pitch extension in males is ∼30° [68] . Assuming the jaw to be in a completely closed position forces rod trajectories to incline more toward the right hemisphere in order to avoid contact with the jaw and breaking it - as may result from the trajectories identified by Damasio et al. Having the jaw open provides a greater number of possible paths which are closer to the vertical axis, which thus does not enforce an intersection of the rod with the right hemisphere ( Fig. S4A , B, D; Fig. S5 A–D). The rod's intersection with white matter fiber tractography was thereby determined ( Fig. S6 ). Movie S1 illustrates the path of the tamping iron through Mr. Gage's skull and the white matter fiber pathways of his left hemisphere.

Normal Subjects

T1 anatomical MRI and 64-direction diffusion tensor images (DTI) from N = 110 right-handed male subjects between the ages of 25 and 36 were selected from the LONI Integrated Data Archive (IDA; http://ida.loni.ucla.edu ). The age range was specifically selected to match the age at which Mr. Gage received his injury (25 years old) as well as the age at which he succumbed as a presumed result of the brain damage he experienced (36 years old). Subjects were all healthy “normals” with no neurological or history of psychiatric illnesses.

Segmentation and Parcellation

Segmentation and regional parcellation were performed using FreeSurfer [69] , [70] , [71] following the nomenclature described in [72] . For each hemisphere, a total of 74 cortical structures were identified in addition to 7 subcortical structures and to the cerebellum. The 82 cortical and sub-cortical label names were assigned per hemisphere to each brain based upon the nomenclature described in Destrieux et al. [72] . Regional parcellation was performed using FreeSurfer [73] , [74] , [75] , [76] (see also above). The numbers of hemispheric partitions in the segmentation was as follows – frontal (21), insula (8), limbic (8), temporal (12), parietal (11), occipital (14), basal ganglia (8), and brain stem (1). The complete coding scheme is as presented describing the parcellation scheme naming convention ( Table 6 ) and their abbreviations ( Table 7 ), which can be used to identify the regional labels in Figs. 2a and 3 .

Connectogram Design

Neuroanatomical structure and connectivity information were graphically depicted in a circular diagram format using freely available Circos software ( [77] , www.cpan.org/ports ). Briefly, Circos is a cross-platform Perl-based application which employs a circular layout to facilitate the representation of relationships between pairs of positions by the use of various graphical elements, including links and heat maps. While traditionally used to render genomic information, Circos can be effectively adapted to the exploration of data sets involving complex relationships between large numbers of factors. In our case, cortical parcellations were represented as a circular array of 165 radially aligned elements representing the left and right cerebral hemispheres, each positioned symmetrically with respect to the vertical axis. We term this representation a “connectogram”. The brain stem was positioned at the most inferior extremity of the Circos ring as a consequence of its inclusion as the only midline structure. In this manner, Circos' ability to illustrate chromosomes was modified for lobar depiction, while its functionality for illustrating cytogenetic bands was modified to represent cortical parcellations. As previously described, each parcellation was assigned an arbitrary but unique RGB color (see below). Parcellations were arranged within each lobe in the order of their location along the antero-posterior axis of the cortical surface associated with the published FreeSurfer normal population atlas [72] . To determine this ordering, the center of mass was computed for the GM surface portion associated with each parcellation, and the order of all parcellations was determined based on the locations of these centers of mass as their distance from the frontal pole increased along the antero-posterior coordinate axis. A LONI Pipeline workflow for the creation of the connectogram images using parcellation and connectivity matrix information is available upon request from the authors. A complete description of the methods for connectogram construction can be found in [78] with applied examples in [79] .

Color Coding Schemes

Each cortical lobe was assigned a unique color scheme: black to red to yellow (Fro), charlotte to turquoise to forest green (Ins), primrose to lavender rose (Lim), pink to lavender to rosebud cherry (Tem), lime to forest green (Par), and lilac to indigo (Occ). Each structure was assigned its unique RGB color based on esthetic considerations; e.g. subcortical structures were colored light gray to black. Color scheme choice and assignment to each lobe were made by taking into account the arrangement and adjacency of lobes on the cortical surface, with the goal of avoiding any two adjacent lobes from having overlapping or similar color schemes which were too similar. The individual colors of the scheme associated with any particular lobe were assigned to every parcellation within that lobe in such a way as to create a distinct contrast when displayed on cortical surfaces ( Fig. S2 ) or on the connectogram graphics ( Figs. 2 and 3 ). The particular regional color mappings employed in this article can be considered arbitrary and are not intended to convey any universal or standard regional color scheme, per se .

Representation of Cortical Metrics

Within the circular framework representing the cortical parcellations, five circular heat maps were generated, each encoding one of five structural measures associated with the corresponding parcellation. Proceeding inward towards the center of the circle in Fig. 2 , these measures were: total GM volume, total area of the surface associated with the GM-WM interface (forming the base of the cortical ribbon), mean cortical thickness, mean curvature and connectivity per unit volume. For subject-level analysis, these measures were computed over the entire volumetric (or areal, as appropriate) extent of each parcellation; for the population-level analysis, they were averaged over all subjects.

Values for each measure were mapped to colors, using a scheme that ranged from the minimum to the maximum of the data set. For example, the cortical thickness t with values ranging from t min to t max was normalized as t 1  = ( t − t min )/( t max − t min ). The latter value was mapped onto a unique color from the color map of choice. Thus, for example, hues at color map extremities correspond to t min and t max , as required. For subcortical structures, brain stem and cerebellum, three measures (area, thickness and curvature) were unavailable on a parcellation-by-parcellation basis; their corresponding heat map entries were consequently left blank.

The connectogram in Fig. 3 , illustrating the effects of the tamping iron lesion, represents the individual regionally-specific network metrics (i.e. betweenness centrality, eccentricity, mean local efficiency, and clustering coefficient) and are colored distinctly to be consistent with the cortical maps of the same but unaffected network metrics presented in Fig. 5A . The inner-most ring of the connectogram in Fig. 3 represents the average proportion of regional GM loss taken across subjects.

Connectivity Calculation

To compute connectivity between regions for each subject, the location of each fiber tract extremity within the brain was identified, while the GM volume associated with each parcellation was also delineated. For those fibers which both originated as well as terminated within any two distinct parcellations of the 165 available, each fiber extremity was associated with the appropriate parcellation. For each such fiber, the corresponding entry in the connectivity matrix (e.g. Fig. S2 ) of the subject's brain was appropriately updated to reflect an increment in fiber count [80] , [81] . Each subject's connectivity matrix was normalized over the total number of fibers within that subject; for population-level analysis, all connectivity matrices were pooled across subjects and averaged to compute probabilistic connection probabilities.

Connectivity Representation

For subject-level connectograms, links were generated between any two parcellations whenever a WM tract existed between them. In population-level analyses, the former was done whenever there was a non-vanishing probability for a WM tract to exist between the two regions ( Fig. 2 ). Links were color-coded by the average fractional anisotropy (FA) value associated with the fibers between the two regions connected by the link, as follows. The lowest and highest FA values over all links ( FA min and FA max , respectively) were first computed. For any given connection i where i  = 1, …, N ( N being the total number of connections), the FA value FA i associated with that connection was normalized as FA′ i  = ( FA i − FA min )/( FA max − FA min ), where the prime indicates the FA i value after normalization. After this normalization, FA′ i values were distributed in the interval 0 to 1, where 0 corresponds to FA min and 1 corresponds to FA max . The interval 0 to 1 was then divided into three subintervals (bins) of equal size, namely 0 to 1/3, 1/3 to 2/3, and 2/3 to 1. For every i  = 1, …, N , link i was color-coded in either blue, green or red, depending on whether its associated FA′ i value belonged to the first, second, or third bin above, respectively. Thus, these bins represent low, medium, and high FA. In addition to encoding FA in the link's color as described, relative fiber density (the proportion of fibers for each connection out of the total number of fibers) was also encoded as link transparency. Thus, within each of the three FA bins described, the link associated with the highest fiber density within that bin was rendered as completely opaque, whereas the link with the lowest fiber density was colored as transparent as possible without rendering it invisible. For example, the link with FA′ i  = 1/3 was colored as opaque blue, whereas the link with the lowest FA′ i value was colored as most transparent blue. Similarly, the link with FA′ i  = 2/3 was colored as opaque green, and the link with the lowest value of FA′ i greater than 1/3 was colored as faintest green. The links associated with the lowest fiber densities were drawn first, and links with progressively larger relative fiber densities were drawn on top of the former. The process was successively repeated by drawing links with higher fiber densities on top of links with lower fiber densities. Thus, links associated with the largest fiber densities were drawn “on top” of all other links.

Representation of Connectivity Affected by Pathology

Links associated with fibers affected by pathology were designed to encode fiber density using the same transparency coding scheme as described in the previous subsection. In contrast with the case of healthy fibers, however, two different color schemes were used to encode pathology. Whenever fibers existed between one cortical region that was affected by pathology and another that was not, the color used to draw the corresponding link was brown. By contrast, links between parcellations that were both affected by pathology were drawn using the color gray. This allows one to visually distinguish between connections that involve only one affected region (brown links) and connections that involve two regions that were both affected (grayscale links) ( Fig. 3 ).

Calculation of Pathology Effects upon GM/WM Volumetrics

The calculation described above estimated the amount of GM that was directly affected by the passage of the rod. To compute the total amount of GM that was affected by pathology, however, it is not sufficient to compute the sum of directly lesioned GM parcellation volumes because pathology-affected GM includes cells with intact somas whose axons were nevertheless injured in at least one location along their paths. In other words, a population of neurons whose GM axons were destroyed or affected in spite of their somas being outside the volume of direct injury should also be taken into account when computing the amount of affected GM. Furthermore, the destruction of fibers originating in some parcellated region r 1 that had been directly affected by pathology could also have affected the GM in parcellations to which r 1 is connected by WM fibers originating in r 1 . Consequently, an appropriate calculation of the total GM volume affected by pathology must take into account available quantitative information concerning the extent to which WM fibers affected by pathology could indirectly affect GM as well. To obtain and interpret such information meaningfully, one can use the measures of GM and WM atrophy described below:

Average Percentages of Brain Regions Intersected by the Rod

The average percentage regional volumes (and their standard deviations) intersected by the rod pooled over N = 110 subjects are listed in Table 7 and illustrated graphically in the connectogram of Fig. 3 .

Network Analysis

Because network theory can provide essential insight into the structural properties of cortical connectivity networks in both health and disease [83] , several network metrics of particular significance were computed for each subject, starting with the degree of each node. In our case, nodes were denoted by parcellated regions and edges were represented by fiber tracts. Nodal degree is the number of edges connected to a node and its calculation has fundamental impact upon many network measures; moreover, node degree distributions are highly informative of network architecture. The entry indexed by i and j in the distance matrix of the graph contains the minimum weighted physical length of the path connecting vertices i and j and was computed using the algebraic shortest paths algorithm [84] . Degree of connectivity is represented as the inner-most ring in Fig. 2 , though was not analyzed further beyond its being utilized in the computations of some of the overall network metrics detailed below.

The measurement of network attributes can be generally broken down into the examination of overall network integration – the measurement of path lengths between nodes in a network and the extent of network-wide interaction and ease of communication between distinct regions; segregation – the extent to which nodes of the network group themselves into separate communities; and small worldness – the quantification of the generally shorter path lengths and higher clustering observed in many biological and technological networks with respect to randomly connected systems [85] . To specifically measure these overall network properties, we chose to focus on three particular metrics. To assess network integration from each subject's connectivity matrix we measured the characteristic path length, a measurement of the global average of a graph's distance matrix [86] . Appropriate to our application, the weighted characteristic path length of a network may be altered as a result of brain trauma [87] . To measure the degree of segregation, we computed the mean local efficiency of each network. Investigating network segregation can be important because it can reveal how much information brain regions are able to exchange as well as the extent to which such regions remain structurally segregated from each other. In this instance, reduced efficiency might be expected as a result of a severe penetrating head wound. Finally, we measured network small worldness , i.e. the ratio comprised of the observed characteristic path length relative to that observed in a random network having the same degree distribution and the observed clustering coefficient relative to that observed in a random network.

Additionally, to characterize the regionally-specific effects of the tamping iron lesion, we also computed several additional graph theoretical measurements for each parcellated brain region. These included 1) betweenness centrality, measuring the number of shortest paths from all vertices to all others that pass through that node, 2) local efficiency, the mean shortest absolute path length of at that node, 3) clustering coefficient, measuring the degree to which a node is nodes in a graph is a member of a cluster or clique, and 4) eccentricity, representing the greatest geodesic distance between that node and any other vertex in the graph. Metrics were computed for each subject and averaged with respect to weighting by subject-wise regional parcellation volume. To be consistent with other studies reporting these regionally-specific values, we chose not to normalize them with respect to those obtained in equivalent random networks. Averages of these metrics are illustrated in Fig. 5a(i–iv) along with linear colorbars indicating the ranges of observed mean values. Effects on these metrics in the presence of the tamping iron can be seen as the first four of the inner-most rings of the connectogram presented in Fig. 3 .

Several additional global as well as local graph metrics were computed but not reported here due to potentially excessive colinearlity, imprecision, or due to recognized difficulty with interpretation. For instance, network modularity [88] was not considered due to the heuristic nature of its computation and tendency to provide unreliable values upon repeated estimation. While many of these other network metrics are well known and have their unique advantages [83] , the ones chosen parsimoniously capture the overall changes in network architecture for this patient and the extent to which his injury would compare to similarly-sized lesions in other areas of the cortex. The Brain Connectivity Toolbox (BCT; https://sites.google.com/a/brain-connectivity-toolbox.net/bct/Home ) was used for all weighted and unweighted connection density- and path-length related graph theoretical computations [84] .

For each of the global graph theory measures described above, the mean and standard deviation was computed for each subject in both intact (healthy) and pathology-affected scenarios (the tamping iron lesion as well as simulated lesions over the brain). As an additional basis, we also performed a degree-preserving randomization process using the BCT for each subject's intact network, computed the aforementioned network measurements, and report these averaged across subjects. Such normalization has been recently advised by Rubinov and Sporns [84] . In our case, this involved 10,000 “rewiring” iterations of the BCT null_model_und_sign (compiled C-code version of the Matlab code from the “the bct-cpp project”; http://code.google.com/p/bct-cpp ) algorithm per region by subject. To accommodate the computational cost of performing such a randomization process, we utilized fully the 1200 node Linux cluster based at the Laboratory of Neuro Imaging (LONI) at UCLA to randomize subjects and regions in parallel. Incidentally, normalization of each network type by its own randomized version has the effect of scaling out differences between networks – lesioned or otherwise – and thus makes the metrics largely insensitive to the effects of network damage. So, to provide a common frame of reference across each network type, the observed metrics for the intact, tamping iron, and simulated lesions were normalized with respect to the degree-preserving randomization of the intact network. Finally, to specifically test the differences between the intact and the tamping iron-lesioned networks between subjects, paired Student's t-tests were applied for each normalized measure to identify significant differences between means at p≤0.01. Results are summarized in Table 5 . Further details on the lesion simulation are provided in the section below.

Equivalent Lesion Simulation and Comparison

To examine the tamping iron lesion's specificity to changes in network structure, we investigated whether changes Gage's brain network properties were significantly different from those that would be expected by chance for the same amount of GM loss located in other regions of the brain. To address this, network properties were computed for a set of simulated lesions systematically positioned over the cortex (excluding the tamping iron lesion itself) and Mr. Gage's network measurements were compared to the distribution of the average metric values taken over subjects and lesions. Specifically, we adopted an approach similar to that of Alstott et al. [89] , who simulated the effects on functional connectivity of targeted lesions distributed in various regions of the cerebral cortex. In our extension of this method, localized area removal was performed by deleting all nodes and their connections within regions consisting of contiguous anatomic parcellations as defined using the methods of Destrieux et al. [72] . In contrast to Alstott et al., however, our structural connectivity simulations also sought to account for additional lesion effects upon WM by modeling the removal of so-called “fibers of passage”. To do so, connectivity network edges between anatomic parcellations neighboring the GM lesion were removed without deleting the corresponding nodes connected by these edges, unless these nodes also belonged to the GM portion of the lesion itself.

The details of our simulation are as follows: 500 distinct lesions were simulated by first populating the cortical surface with 500 distinct sets of contiguous parcellations. Each of these sets was subsequently used as a synthetic “lesion”, subject to the constraints that the percentages of WM and GM lost due to the lesion were the same as had been estimated for Gage's tamping iron injury. This process was repeated until 500 distinct lesions were created uniformly across the brain, and the procedure was repeated for all 110 subjects included in the study. To ensure that each of the lesions had approximately the same position in each subject, lesion configurations were defined using the cortical atlas of Fischl, Dale et al. [71] , and the corresponding location of every lesion in each subjects was identified by mapping the lesion configuration from the atlas to each subject's cortical surface using existing/published FreeSurfer methodology [70] , [90] , [91] . Thus, by the process described above, 500 distinct lesions that were identical in size to Gage's from the standpoint of percentage WM and GM loss were created uniformly over the brain in each of the 110 subjects. Subsequently, each lesion's effect on overall network properties was computed. Global network metrics were then pooled over all subjects and simulations so as to obtain the average (i.e. most probable) value of every metric for each of the 500 simulated lesioned networks.

Finally, we compared the observed effects of the tamping iron lesion on the random network normalized graph theory measures of integration and segregation against that observed for all remaining lesions. Computed as Z-statistics, the results of these comparisons are illustrated graphically for network integration and segregation in Fig. 5c (i and ii) , respectively, and are colored to show those effects most similar to the tamping iron lesion (black), moderately similar (orange), and most dissimilar (white). Generally, as one moves posteriorly away from the Gage lesion site, similarity on network effects tends to be reduced. However, exceptions exist in bilateral post-central gyrus and the left superior and posterior portion of the parahippocampal gyrus.

Supporting Information

The LONI Pipeline Workflow Environment. We applied the LONI Pipeline [93] , [94] for segmentation and registration of the input MRI image volume data, the processing of all DTI tractography, and computation of tract statistics. This grid-based solution provides validation and distribution of new computational tools, and an intuitive graphical interface for developing and executing parallel volumetric processing software. See http://pipeline.loni.ucla.edu for additional details.

https://doi.org/10.1371/journal.pone.0037454.s001

Views of the cortical parcellation of a sample subject. Top rows show the lateral, anterior, and dorsal surfaces; second row shows medial, posterior, and ventral pial surfaces, while the bottom two rows show the same orientations but as inflated pial surfaces to more adequately present the extent of regional parcellations and their color coding. The arbitrarily chosen regional colors are the same as those of the outer-most ring in Figure 2 and whose RGB values are referenced Table 5 are shared by the outer most ring of brain regions on the connectogram images permitting rapid cross-reference.

https://doi.org/10.1371/journal.pone.0037454.s002

Connectivity Matrix. Each row and each column represent distinct parcellated regions where in each cell i,j was computed the number of fibers that were found to begin or end in each region pair, the average FA, and the average fiber length over subjects.

https://doi.org/10.1371/journal.pone.0037454.s003

Modeling of the Skull Fragmentation and the Rod. a) Models of the eyeballs were placed in to the ocular cavities in order to use them as constraints for the trajectory of the tamping iron. According to Harlow's account, the left orbit was extended outward “by half its diameter”. b) The bones of the skull representing the major breakages were systematically labeled and can be independently manipulated using Slicer. The mandible was also rotated downward and laterally in order to allow the tamping iron not to impinge on it and also to comply with Harlow's account that Gage was in the act of speaking to his men at the moment of the blast. c) The surface model of the Gage skull, with closed mandible, along with the surface of the life mask commissioned by Bigelow. d) A view looking superiorly along the tamping iron's computed trajectory noting how the iron displaced the left anterior frontal bone as it passed.

https://doi.org/10.1371/journal.pone.0037454.s004

Illustrating the Intersection of the Rod and the Brain. a) A figure showing the passage of the rod through the skull with the bones above the cranial “cap” cut at Harlow's direction, and its intersection with the left anterior white matter fiber pathways of an example subject. The complementary hemisphere is displayed to illustrate that the rod did not intersect that hemisphere. b) A view of the rod displacing the bones of the skull. c) A close up, coxial view of the inferior portion of the iron along its trajectory. d) The intersection of the tamping iron with the left frontal cortex with each major bone fragment removed.

https://doi.org/10.1371/journal.pone.0037454.s005

The Effects of the Tamping Iron on White Matter Fiber Tractography. a) A view of the Gage skull with the white matter fiber tracts of an example subject warped to the space. In this view, fibers which intersect the rod's pathway have been removed. b) A transaxial view of the DTI fiber pathways remaining after those which were intersected by the rod had been removed. c) The fibers intersected by the rod connect areas of cortex throughout the left cerebral hemisphere as well as between hemispheres. d) A sagittal view of the fibers experiencing damage by the tamping iron. All bone fragments and the cranial “cap” have been removed.

https://doi.org/10.1371/journal.pone.0037454.s006

Movie of The Effects of the Tamping Iron on White Matter Fiber Tractography. This movie rendering illustrates the passage of the tamping iron through the Gage skull and its intersection with left hemispheric white matter fiber pathways. The right hemispheric cortical surface model is displayed to illustrate that the rod did not cross the midline to damage right frontal cortex. The rendering was created using 3D Slicer ( http://slicer.org ).

https://doi.org/10.1371/journal.pone.0037454.s007

Acknowledgments

The authors wish to acknowledge the assistance of Dominic Hall, Curator, Warren Anatomical Museum, Center for the History of Medicine, Francis A. Countway Library of Medicine 10 Shattuck Street, Boston, MA 02115 for access to Mr. Gage's skull, life mask, and tamping iron. We also express our gratitude to Marianna Jakab of the Surgical Planning Laboratory at Harvard Medical School for assistance with the CT image volumes, and to Drs. Danielle Bassett (Department of Physics, University of California Santa Barbara), Randal McIntosh (Rotman Institute, Toronto, Canada), and Paul M. Thompson (Department of Neurology, University of California Los Angeles) for their input and guidance on our network theoretical analyses. We are also extremely grateful for the rigorous and thorough comments of two anonymous reviewers on earlier versions of this article. Finally, we are indebted to the dedicated staff of the Laboratory of Neuro Imaging (LONI) at UCLA.

Author Contributions

Conceived and designed the experiments: JVH AWT RK. Performed the experiments: MCC CMT AI. Analyzed the data: MCC AI CMT. Contributed reagents/materials/analysis tools: AWT RK. Wrote the paper: JVH AI MCC. Provided computational resources and database access needed for neuroimaging data analysis: AWT. Provided access to data essential for the study: RK.

• 1. Rosen B, Wedeen V, Van Horn JD, Fischl B, Buckner R, et al. (2010) The Human Connectome Project. Organization for Human Brain Mapping Annual Meeting. Barcelona, Spain.
• View Article
• Google Scholar
• 3. Macmillan M (2000) An Odd Kind of Fame: Stories of Phineas Gage. Boston: MIT Press. 576 p.
• 4. Damasio AR (1995) 336 p. Descartes' Error: Harper Perennial.
• 83. Sporns O (2011) Networks of the Brain. Cambridge, MA: MIT Press.