Gene expression patterns following unilateral traumatic brain injury reveals a local pro-inflammatory and remote anti-inflammatory response
BMC Genomics volume 14, Article number: 282 (2013)
Traumatic brain injury (TBI) results in irreversible damage at the site of impact and initiates cellular and molecular processes that lead to secondary neural injury in the surrounding tissue. We used microarray analysis to determine which genes, pathways and networks were significantly altered using a rat model of TBI. Adult rats received a unilateral controlled cortical impact (CCI) and were sacrificed 24 h post-injury. The ipsilateral hemi-brain tissue at the site of the injury, the corresponding contralateral hemi-brain tissue, and naïve (control) brain tissue were used for microarray analysis. Ingenuity Pathway Analysis (IPA) software was used to identify molecular pathways and networks that were associated with the altered gene expression in brain tissues following TBI.
Inspection of the top fifteen biological functions in IPA associated with TBI in the ipsilateral tissues revealed that all had an inflammatory component. IPA analysis also indicated that inflammatory genes were altered on the contralateral side, but many of the genes were inversely expressed compared to the ipsilateral side. The contralateral gene expression pattern suggests a remote anti-inflammatory molecular response. We created a network of the inversely expressed common (i.e., same gene changed on both sides of the brain) inflammatory response (IR) genes and those IR genes included in pathways and networks identified by IPA that changed on only one side. We ranked the genes by the number of direct connections each had in the network, creating a gene interaction hierarchy (GIH). Two well characterized signaling pathways, toll-like receptor/NF-kappaB signaling and JAK/STAT signaling, were prominent in our GIH.
Bioinformatic analysis of microarray data following TBI identified key molecular pathways and networks associated with neural injury following TBI. The GIH created here provides a starting point for investigating therapeutic targets in a ranked order that is somewhat different than what has been presented previously. In addition to being a vehicle for identifying potential targets for post-TBI therapeutic strategies, our findings can also provide a context for evaluating the potential of therapeutic agents currently in development.
It is estimated that 3.17-3.32 million persons in the United States are living with long-term or lifelong effects of traumatic brain injury (TBI) . Approximately 1.7 million new TBIs are sustained each year  resulting in 53,000 deaths  and as many as 125,000 additional people that have long-term behavioral deficits . Efforts to implement preventative measures for the leading causes of TBI (motor vehicle-traffic accidents, falls, and assaults) are critical [2, 3] but will not eliminate TBI as a major public health problem. Development of effective clinical treatment protocols post-TBI and potential prophylactic agents for military usage are necessary to address those TBIs that cannot be prevented.
Injury to the brain results in irreversible damage at the site of impact and initiates cellular and molecular processes that lead to delayed or secondary neural injury in the surrounding tissue [5, 6]. Neuroprotective strategies target these processes in an attempt to halt the progression of or prevent the delayed injury . These processes include inflammation, damage to the blood brain barrier, release of excitatory amino acids, oxidative stress, cerebral edema, reduced cerebral blood flow, hypoxia, and ischemia [5, 6]. In the current study, we used gene microarray analysis to determine which genes, pathways and networks were significantly altered after unilateral controlled cortical impact (CCI), an experimental model of TBI.
Microarray technology is a very powerful tool that allows the user to examine expression profiles for thousands of genes at one time. The true power of microarray technology can be maximized when the expression profiles can be attributed to significant alterations in biological functions and molecular pathways. Fortunately, the recent development of advanced bioinformatic analysis tools has made the utilization of microarray data more practical and allows for easier replication. One such tool is the Ingenuity Pathway Analysis (IPA) software program which uses a database built from published scientific literature to draw direct and indirect interactions between genes and to assign genes to specific biological functions and canonical pathways . IPA was used here for functional, canonical pathway, and network analysis of the genes that were altered by TBI. We examined gene expression on both sides of the brain in order to understand the alterations both locally (ipsilateral) and remotely on the opposite side of the brain (contralateral). We observed that the contralateral side of the brain, which has often been used as a control in similar experiments, exhibited significant alterations in gene expression following TBI.
The bioinformatic analysis tools we used also allowed for both identification of key molecules and elucidation of their interactions with each other. We used these interactions to identify molecules and molecular pathways central in the response to TBI that could provide novel targets for therapeutic strategies. Understanding how local and remote gene expression profiles change following TBI, when compared to non-injured brain, will provide valuable insight into delayed neuronal injury mechanisms as well as intrinsic neuroprotective processes.
Principal component analysis
The ipsilateral, contralateral, and naïve gene datasets generated 24 hours post-TBI were analyzed by principal component analysis (PCA) to assess the variability of our microarray data. Nine principal components were generated and the first 3 principal components explain 91.267% of the variance in our microarray data. The 3D score plot generated from these identified injury status as the major source of variability (Figure 1A). Similarity between data points in the PCA was determined by the distance between the points, with shorter distance indicating increased similarity. Each data point represents one animal's gene expression profile. Ipsilateral, contralateral and naïve all clustered together by injury status and each group was well isolated from the other two groups.
A total of 69 biological functions met IPA's threshold and cutoff conditions for the TBI-I (ipsilateral vs naïve) dataset. Inspection of the top 15 biological functions associated with the dataset in IPA showed that 7 of these functions are explicitly cellular in nature; cellular movement, cellular growth and proliferation, cellular development, cell death, cell cycle, cell morphology, and immune cell trafficking (Figure 1B). Of the other top functions, 5 are disease and disorder related; cancer, neurological disease, genetic disorder, skeletal and muscular disorders, and dermatological diseases and conditions. We posited that the majority of these functional categories had an inflammatory component and inflammatory response (IR) itself was ranked as the 15th most significant biological function for our dataset. To further investigate this, we calculated the percentage of gene overlap each function had with inflammatory response. This overlap ranged from 26% to 93% with 9 biological functions falling in the 35% to 55% range (Table 1). The cumulative overlap of IR with the other 14 top biological functions was 39.3%. Therefore, we chose to focus our subsequent analysis techniques on the inflammatory response genes in this study because of the significant involvement of inflammatory genes in all of the top ranked biological functions and the role of the IR in secondary neural injury [5, 6, 9, 10].
Histology and immunohistochemistry
Histology and immunostaining were used to assess damage to and the inflammatory state of the brain. Figure 2 demonstrates cortical damage to the brain tissue ipsilateral to the injury by Fluoro-Jade B (FJB) and DAPI staining (Figure 2A,C). Near the site of impact, the integrity of the cortex is disrupted and there is obvious tissue loss. ED-1 immunostaining showed that many activated microglia/macrophages were present in the area of tissue damage (Figure 2D). The contralateral side of the brain showed no overt structural damage or ED-1 immunoreactivity (Figure 2B,E,F).
The hippocampal region did not display overt structural damage on either side of the brain (Figure 3C,F) but did show cellular damage ipsilaterally when stained with FJB (Figure 3A). ED-1 staining similarly showed only macrophage/microglial activation in the hippocampus on the ipsilateral side of the brain (Figure 3D). No cortical or hippocampal labeling was seen in negative control sections (data not shown). Any apparent staining in the contralateral hemisphere represented autofluorescence and was similar to levels detected in negative controls.
In contrast to ED-1, CD11b stains all microglia regardless of activation state. When the microglia transform into activated macrophages they also undergo a morphological transformation from a resting ramified state to an amoeboid state. Figure 4 demonstrates that this morphological shift only occurred in areas of cellular damage. There were no amoeboid shaped microglia detected on the contralateral side (Figure 4A,B). Ipsilaterally, both microglial morphologies were observed. Ramified microglia were seen in subcortical regions that did not exhibit cellular damage (Figure 4C) and amoeboid microglia were present at the site of impact (Figure 4D) and in subcortical regions where cellular damage had occurred. This duality can be seen definitively in the ipsilateral hippocampal region where there is cellular damage interspersed between structurally intact tissue (Figure 4E,F). No immunostaining was seen in the cortex or hippocampus in negative control sections (data not shown).
Inflammatory gene expression patterns
Focusing on the inflammatory response genes in our datasets, we determined that 372 IR genes had a greater than 2-fold change in expression. Of these genes, 146 genes changed on both the ipsilateral and contralateral sides of the brain. In order to determine whether these common genes changed differently on one side of the brain compared to the other, we calculated the ratio of the TBI-I fold change to the TBI-C (contralateral vs naïve) fold change. Those genes that had a TBI-I/TBI-C ratio greater than 1.75 were determined to have change differently. We observed that 109 of the common IR genes (75%) changed similarly (TBI-I/TBI-C ratio < 1.75; Figure 5A). Of the genes that changed similarly, 79 genes (54%) increased in expression and 30 genes (21%) decreased in expression. The remaining 37 common IR genes (25%) were changed differently (TBI-I/TBI-C ratio > 1.75) (Figure 5A). Table 2 shows the 37 common IR genes that change differently. These genes span all cellular compartments (extracellular space, plasma membrane, cytoplasm, and nucleus) with diverse molecule types. The expression of all these genes was lower on the contralateral side of the brain. Because of their different expression patterns, these genes became our first group of genes of interest (GOI). Notable genes identified included CCND1, SPP1, ERAP1, LYN, THRA, TIMP1, the transcription regulators STAT3, CEBPD, and CBL, and the plasma membrane receptors and signaling molecules IL6ST, CD44, EGFR, ITGA5, and SDC1.
There were 188 IR genes that changed uniquely on the ipsilateral side of the brain. 179 of those genes (95%) increased while 9 genes (5%) decreased in expression (Figure 5B). Only 38 IR genes change uniquely on the contralateral side of the brain and, in contrast to what we observed on the ipsilateral side, only 10 genes (26%) increased while 28 genes (74%) decreased in expression (Figure 5C).
Canonical pathway analysis
We used canonical pathway and network analysis in IPA to identify genes in our datasets that were most relevant to the observed IR. Because IPA displays only the relative expression values, we defined GOI, in this context, as those genes that either increased or decreased on one side of the brain and showed no change in expression on the other side or genes that had opposite changes in expression. These genes were easily identified by side-by-side comparison of the canonical pathway and gene networks overlaid with the expression values of the TBI-I and TBI-C IR datasets. Canonical pathways in IPA are well-characterized metabolic and cell signaling pathways derived from information found in specific journal articles, review articles, text books, and KEGG Ligand . Figure 6 shows the IL-6 signaling canonical pathway. IPA determined that this pathway was highly associated with the TBI-I dataset. Additionally, this pathway also includes elements of IL-1, TNF-α, and lipopolysaccharide (LPS) signaling. IL-1 and TNF-α were previously associated with TBI and inflammation [9, 12–14] and TBI induced inflammation has been shown to be exacerbated by LPS challenge . By overlaying the relative expression values for all TBI-I IR genes (Figure 6A) and all TBI-C IR genes (Figure 6B), we were able to identify a number of GOI that were increased in TBI-I and were either unchanged or decreased in TBI-C, including IL1B, several transmembrane cytokine receptors and the transcriptional regulators NFkB, STAT3, CEBPB (NF-IL6), and FOS.
Gene network analysis
In contrast to canonical pathways, which are relatively immutable in IPA, gene networks are generated de novo in IPA based on the list of genes that are imported. IPA takes “seed” molecules from the gene list, searches the Ingenuity Knowledge Base, and uses a network algorithm to draw connections between molecules based on biological function . In order to generate networks that included IR genes that changed on both sides of the brain, we combined the TBI-I and TBI-C datasets and performed an IPA core analysis on that union dataset. IPA scores the networks in order to rank them according to their degree of relevance to the network eligible molecules in your dataset . Figure 7 shows the highest scored network associated with our union dataset. For this network, we overlaid the relative expression values for the unique TBI-I IR genes (Figure 7A) and unique TBI-C IR genes (Figure 7B) and were able to identify 29 GOI, such as the chemokine CXCL10, a number of cytokine and toll-like receptors, heat shock proteins, and transcriptional regulators that were increased in TBI-I and were either unchanged or decreased in TBI-C. Interestingly, IRF2 was upregulated on the contralateral side of the brain.
Figure 8 shows another network that was scored in the top three networks associated with our union dataset. For this network, we overlaid the relative expression values for all TBI-I IR genes (Figure 8A) and all TBI-C IR genes (Figure 8B) and were able to identify 15 GOI including the cytokine CXCL1, a number of genes in the JAK/STAT pathway, NFkB and several genes associated with apoptosis.
Because the IR is, in part, humoral in nature, we created a network within IPA by seeding it with IR cytokines and growth factors expressed uniquely in the TBI-I dataset and “growing them” (making direct functional connections) to the TBI-I/TBI-C union dataset. Figure 9 shows the resultant network. By overlaying the relative expression values for all TBI-I IR genes (Figure 9A) and all TBI-C IR genes (Figure 9B), we were able to identify 23 more GOI. Notable genes were IL1B, chemokines CCL7, CCL13, CXCL13 and CCL4, peptidases MMP13, MMP3, and MMP9 and many inflammation associated transcription factors.
Compiling the gene interaction hierarchy (GIH)
By combining the GOI identified through canonical pathway and network analysis with those in Table 2, we identified a total of 114 GOI. In order to determine which genes might be most relevant in the IR, we ranked these genes relative to each other by the number of direct interactions each had with the other GOI. Our analysis showed that 95 of the GOI formed an interconnected network, leaving 19 “orphan” genes (see Additional file 1). Genes having more than 10 connections (1st order) within the main GOI network were considered “primary” in this analysis (see Additional file 2 for an example). Genes having 5–10 connections were considered “secondary” (see Additional file 3 for an example) and those with less than 5 connections were considered “peripheral”. The resultant GIH is displayed in Table 3.
In this study, we used microarray technology to examine cellular and molecular mechanisms associated with secondary brain injury following TBI. Our findings indicated that the inflammatory response and its associated genes and pathways are significant in the post-TBI molecular response. This is consistent with our published studies indicating that inflammation is involved with delayed, secondary neuronal injury following other acute brain injuries (ABI’s) including stroke and neurotoxin exposure [17–20]. Activated microglia, astrocytes and macrophages have been shown to be the source of several of the inflammatory molecules identified in this study [10, 21–23]. The presence of activated inflammatory cells ipsilaterally and their absence contralaterally was confirmed by immunohistochemistry. This cellular inflammation offers further evidence for pro-inflammatory molecules being produced in the region of tissue damage following TBI.
Despite the absence of activated inflammatory cells on the contralateral side of the brain, our results also showed that TBI resulted in a significant alteration of the inflammatory gene response on the both sides of the brain. The observed IR gene expression pattern suggests that there is a baseline IR throughout the whole brain due to unilateral CCI. This is demonstrated by the nearly 30% of all IR genes that change similarly on both sides of the brain (109 of 372 IR genes). Above this baseline of inflammation, a distinct IR gene expression pattern emerges. The expression level is higher on the ipsilateral side for all the remaining genes expressed on both sides of the brain. 95% of the IR genes expressed only in ipsilateral tissues increase in expression while 74% of the IR genes expressed only in contralateral tissues decrease in expression. Examination of the contralateral gene expression in our GIH showed that 11 genes (3 primary tier, 3 secondary tier, 3 peripheral tier and 2 orphan) decreased in expression contralaterally while either remaining unchanged or increasing in expression ipsilaterally. Only 3 genes (1 secondary and 2 peripheral tier) increased contralaterally and remained unchanged ipsilaterally. This demonstrated that the most significant contralateral IR genes in this analysis show decreased expression.
Taken together with the cellular inflammation data, we can surmise that increased expression of the majority of IR genes ipsilaterally results in the development of functional inflammation that can contribute to secondary neural injury. Similarly, it is likely that the suppression of the majority of IR genes contralaterally prevents this development in brain regions remote to the injury. This summary of the overall IR gene expression does not take into account that some IR genes may be classified as anti-inflammatory in certain contexts. However, any counter effects of anti-inflammatory molecules seems to be negligible as the cellular inflammation pattern follows the IR gene expression pattern assuming the majority of IR genes are pro-inflammatory. The contralateral inflammatory response without detectable cellular inflammation also fits well with the idea that some inflammation may actually be beneficial following TBI .
Through canonical pathway and network analysis combined with identification of the common genes that change differently on each side of the brain, we identified 114 GOI. Many of these genes have been previously associated with the IR following ABI’s (i.e., TBI, stroke, nerve agent exposure) including CCL13/CCL2, CCL4, CCL6, CCL7, CCR5, CD14, CD44, CDKN1A, CEBPB, CEBPD, CXCL13, CXCL9, CSF1, FOS, HBEGF, HSPA1A/HSPA1B, ICAM1, IER3, IL-1β, IRF1, IRF2, JAK2, LCN2, MMP3, MMP9, MYD88, NF-κB, PSMB8, S100A4, SPP1, STAT3, TLR2, TLR4, and TNFRSF1A [10, 13, 21, 24–34]. These results support the utility of our methods for identifying the significant genes related to this biological function. Additionally, several IR genes that appear to be novel in the context of TBI were also identified, providing new targets for future study. These genes included HSP90AA1, ERAP1, PSMB9, CBL, BTK, RORA, THRA, and ITGA5. We wanted to take our analysis one step further and determine which genes were likely the most critical in the observed IR. We accomplished this by creating a network of these genes and determining how many 1st order connections each gene had with the other genes in the network. A GIH was created based on these numbers and there were some intriguing findings in terms of which tier certain molecule types predominated.
Not surprisingly, a large number of transcription regulators were included in the primary and secondary tiers as these molecules are the point of convergence for many of the inflammatory pathways and the regulatory step in the production of new proteins. Somewhat surprising was the large number of cytokines that fell into the peripheral tier since considerable focus has been placed on cytokines as mediators of inflammation and targets for therapeutic intervention [10, 21, 24]. This result may be due to a near 1-to-1 relationship that these cytokines have with their receptors, limiting the 1st order connections in the GOI network. While cytokines clearly play an important role in initiating the IR, they may not be the most critical molecules in modulating the IR. Only one cytokine, IL1B, is in the primary tier of our hierarchy while CCL7, CXCL10, and SPP1 are in the secondary tier. This suggests that neuroprotective strategies directed at one of the lower tier cytokines may not be sufficient for limiting the IR . Further, while therapeutic hypothermia after TBI has been shown to improve outcome, inflammatory cytokine levels were unaffected, implying minimal cytokine involvement in the observed neuroprotection . However, targeting the cytokines in the higher tiers may produce effective modification of the IR . Our GIH shows that regulation of transcription, phosphorylation (kinases), extracellular matrix/cell adhesion (FN1, MMPs, ICAM1), and receptors (transmembrane and G-coupled) figure prominently (primary and secondary tiers) in the post-TBI IR. These molecules may, therefore, be more efficient targets for therapeutic strategies to combat post-TBI inflammation because they are activated regardless of the initiating factor.
Two well characterized signaling pathways stood out in our GIH: toll-like receptor/NF-κB signaling and JAK/STAT signaling. There was some concern that our choice of canonical pathway may have skewed our analysis towards these signaling pathways (toll-like receptor/NF-κB and JAK/STAT) because large portions of both of these pathways are included in IL-6 signaling (Figure 4). While 12 GOI are included in that pathway, 8 of these genes were also identified by other analyses used to identify GOI, further supporting their importance to the post-TBI IR. Additionally, interconnection of the genes in IL-6 signaling could not account for the total number of 1st order connections for these genes in the GOI network used to build the GIH. This was especially true for genes in the primary tier. 19 GOI were either a part of (CEBPB, FOS, NFKB2, EGFR, IKBKB, IL1B, TLR4, MYD88, TLR2, TNFRSF1A, CD14, IL1R2), a product of (IRF1, FN1, CCND1, CASP3), or associated with (IRF7, LYN, HSP90AA1) toll-like receptor/NF-κB signaling . 15 of these genes were in the primary tier of the GIH. Secondarily, 5 GOI (JAK2, STAT3, CDKN1A, IL6R, IL6ST) were a part of JAK/STAT signaling , and three of these genes were in the primary tier. Based on our analysis, these signaling pathways are likely candidates for the induction of post-TBI IR. Therefore, therapeutic strategies aimed at the molecules in these pathways may reduce TBI-induced inflammation and, by extension, post-TBI neuronal death.
It should be noted that this study only takes a snapshot at 24 hours post-injury of a dynamic and evolving molecular process initiated by TBI . Therefore, significant molecular events that precede and follow this time point are not reflected in these data. For a more complete examination of the molecular response to TBI, a study including multiple time points before and after 24 hours will be necessary and are underway in our laboratory. Further, confirmation of the biological relevance of any observed gene expression profile following TBI is a critical next step to exploring potential therapies for brain trauma [26, 38].
Microarray analysis is a powerful tool that allows for the analysis of thousands of genes simultaneously. We demonstrated that TBI was associated with a powerful pro-inflammatory response in ipsilateral brain tissues. We also noted the distinct IR gene expression pattern that suggests a remote anti-inflammatory response. The use of multiple network and pathway analyses to identify GOI aided in making our datasets manageable and revealed 2 distinct pathways, toll-like receptor/NF-κB signaling and JAK/STAT signaling, associated with post-TBI secondary neural injury. Our GIH provides a starting point for investigating therapeutic targets in a ranked order that is somewhat different than what has been presented previously in microarray studies. In addition to being a vehicle for identifying potential targets for post-TBI therapeutic strategies, our GIH can also provide a context for evaluating the potential of therapeutic agents currently in development.
All animals used in these studies were treated humanely and with regard for alleviation of suffering and pain and all protocols involving animals were approved by the IACUCs of Morehouse School of Medicine and/or The Georgia Institute of Technology prior to the initiation of experimentation. Adult male Sprague–Dawley rats (290-300 g; Charles River Laboratories International, Inc., USA) were housed individually in standard plastic cages in a temperature-controlled room (22 ± 2°C) on a 12 h reverse light–dark cycle. Food and water were provided ad libitum.
Controlled cortical impact
Under isoflurane anesthesia, rats received a unilateral controlled cortical impact (CCI/TBI) using the Pittsburgh Precision Instruments, Inc. device. A craniotomy was made with the center 4 mm posterior and 3–4 mm lateral to bregma using a 6 mm diameter trephan drill bit. The impact was done at an angle of 15° from vertical with a velocity of 3 m/s to a depth of 2 mm using a 5 mm diameter impact tip. The rats were sacrificed 24 h post-injury and the brains were removed for RNA isolation or histology.
RNA preparation and GeneChip analysis
The ipsilateral hemi-brain tissue at the site of the injury, the corresponding contralateral hemi-brain tissue, and naïve (control) brain tissue (n = 3 for each) were used for RNA isolation. Total RNA was extracted with TRIzol Reagent (Life Technologies, Rockville, MD, USA) and cleaned (RNAqueous Kit, Ambion, Austin, TX, USA). The RNA was prepared for microarray hybridization with the GeneChip® 3′ IVT Express Kit (Affymetrix Inc., Santa Clara, CA, USA) aRNA amplification procedure. Briefly, total RNA was reverse transcribed to synthesize first-strand cDNA containing a T7 promoter sequence. The single-stranded cDNA was converted into a double-stranded DNA template for transcription. The reaction employed DNA polymerase and RNase H to simultaneously degrade the RNA and synthesize second-strand cDNA. In vitro transcription generated multiple copies of biotin-modified aRNA from the double-stranded cDNA templates (this was the amplification step). aRNA Purification removed unincorporated NTPs, salts, enzymes, and inorganic phosphate to improve the stability of the biotin-modified aRNA. Finally, the labeled aRNA was fragmented to prepare the target for hybridization to GeneChip® 3′ expression arrays . Following fragmentation, 15 μg of the biotinylated cRNA was hybridized to an Affymetrix Rat Genome 230 2.0 GeneChip. The chips were hybridized at 45°C for 16 h, and then washed, stained with streptavidin–phycoerythrin and scanned according to manufacturing guidelines.
Microarray data analysis
Data analysis was performed using Affymetrix Expression Console™ software that supports probe set summarization and CHP file generation of 3′ expression using the MAS5 Statistical algorithm. Affymetrix microarrays contain the hybridization, labeling and housekeeping controls that help determine the success of the hybridizations. The Affymetrix Expression Analysis algorithm uses the Tukey’s biweight estimator to provide a robust mean Signal value and the Wilcoxon’s rank test to calculate a significance or p-value and Detection call for each probe set. The Detection p-value is calculated using a Discrimination Score [R] for all probes. The Discrimination Score is a basic property of a probe pair that describes its ability to detect its intended target. It measures the target-specific intensity differences of the probe pair (perfect match (PM) – mismatch (MM)) relative to its overall hybridization intensity (PM + MM). Background estimation is provided by a weighted average of the lowest 2% of the feature intensities. Mismatch probes are utilized to adjust the perfect match (PM) intensity. Linear scaling of the feature level intensity values, using the trimmed mean, is the default to make the means equal for all arrays being analyzed. False-negative and false-positive rates are minimized by subtracting nonspecific signal from the PM probe intensities and performing an intensity-dependent normalization at the probe set level. Three chips were used for each experimental group: ipsilateral, contralateral and naïve control. The datasets produced by the Affymetrix software contain gene identifiers and corresponding expression values. The datasets used for this study can be accessed in the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) of the National Center for Biotechnology Information (NCBI) with accession number GSE45997. These data were analyzed in Microsoft Excel for calculation of fold change and whether the genes were confirmed as present in the tissue sample (as determined by the Affymetrix software). Genes in the injured brain that increased or decreased in expression by 2-fold or more compared to controls and were present in either all 3 ipsilateral samples or all 3 contralateral samples were identified. The gene datasets that were generated were ipsilateral vs naïve (TBI-I) and contralateral vs naïve (TBI-C) fold changes.
Principal component analysis
PCA was carried out using the Gene Expression Similarity Investigation Suite software (Genesis; Graz University of Technology). The microarray datasets were analyzed and a matrix was constructed to determine distribution of variants. To compute the principal components (PCs), Genesis calculates the n eigenvalues and their corresponding eigenvectors are calculated from the (n x n) distance matrix using Singular Value Decomposition (SVD), where n = number of genes . The score for each experimental dataset per PC was determined by multiplying the initial expression value by the Genesis calculated eigenvector for each gene and adding up the resulting values. The scores were plotted on the X, Y, and Z axes resulting in a 3D scores plot where the coordinates are (PC1 score, PC2 score, PC3 score).
Ingenuity pathway analysis
The gene datasets were analyzed using Ingenuity Pathway Analysis (Ingenuity® Systems, http://www.ingenuity.com) and overlaid onto a global molecular network developed from information contained in the Ingenuity Knowledge Base. Fischer’s exact test was used to calculate a p-value determining the probability that each biological function and/or disease assigned to that network is due to chance alone. The functions, canonical pathways, and gene networks that were most significant to the dataset were identified. Gene expression profiles were overlaid on the canonical pathway and gene network figures to reveal similarities and dissimilarities in their gene expression patterns. Gene networks were also created using Ingenuity Knowledge Base to further understand specific interactions between our genes of interest.
We used the following formulas to calculate the ratio of TBI-I to TBI-C fold changes: (1) Gene increased on both sides: ratio = (TBI-I)/(TBI-C); (2) Gene decreased on both sides: ratio = 1/[(TBI-I)/(TBI-C)]; and (3) Gene increased ipsilaterally and decreased contralaterally: ratio = (TBI-I)/-[1/(TBI-C)].
Histology and immunohistochemistry
At 24 h post injury, rats were anesthetized with an intraperitoneal injection of a ketamine:xylazine:acetylpromazine cocktail (50:10:1.67 mg/kg respectively) and perfused transcardially with saline followed by cold 4% paraformaldehyde solution in PBS for 30 min. Brains were quickly removed and cryoprotected in 30% sucrose. The brains were then frozen in OCT mounting medium and stored until sectioning. Coronal sections of 20 μm thickness were cryosectioned from the perilesional brain area of each animal. Sections were mounted on slides which were stored at -80°C until further processed. Fluoro-Jade® B (FJB; AG310, Millipore, Billerica, MA) labeling was performed as previously described . Immunohistochemical localization of macrophages and activated microglia was performed using antibodies against ED-1 (1:500, MAB1435, Millipore) and CD11b (1:500; CBL1512, Millipore). After rinsing in 0.01 M PBS, sections were blocked with buffer containing 5% normal goat serum and 0.3% triton-x 100 for 1 h at 4°C and then incubated for 1 h at 37°C with the primary antibodies. Sections incubated with antibodies to ED-1 and CD11b were washed with PBS and incubated with DyLight 594 and DyLight 488 conjugated goat anti-mouse secondary antibodies, respectively (1:400; Jackson ImmunoResearch Laboratory, West Grove, PA) for 1 h at room temperature. Negative control sections for immunohistochemistry were incubated with the secondary antibody only (no primary antibodies). A Zeiss microscope equipped with a CCD camera (Carl Zeiss Microimaging Inc, Thornwood, NY) was used to capture digital images of the sections.
Controlled cortical impact
Gene interaction hierarchy
Genes of interest
Ingenuity Pathway Analysis
Principal component analysis
Traumatic brain injury
Contralateral vs naïve gene dataset
Ipsilateral vs naïve gene dataset.
Zaloshnja E, Miller T, Langlois JA, Selassie AW: Prevalence of long-term disability from traumatic brain injury in the civilian population of the United States, 2005. J Head Trauma Rehabil. 2008, 23 (6): 394-400. 10.1097/01.HTR.0000341435.52004.ac.
Faul M, Xu L, Wald MM, Coronado VG: Traumatic brain injury in the United States: emergency department visits, hospitalizations and deaths 2002–2006. 2010, Atlanta (GA): Centers for Disease Control and Prevention, National Center for Injury Prevention and Control
Coronado VG, Xu L, Basavaraju SV, McGuire LC, Wald MM, Faul MD, Guzman BR, Hemphill JD: Surveillance for traumatic brain injury-related deaths–United States, 1997–2007. MMWR Surveill Summ. 2011, 60 (5): 1-32.
Selassie AW, Zaloshnja E, Langlois JA, Miller T, Jones P, Steiner C: Incidence of long-term disability following traumatic brain injury hospitalization, United States, 2003. J Head Trauma Rehabil. 2008, 23 (2): 123-131. 10.1097/01.HTR.0000314531.30401.39.
McIntosh TK, Saatman KE, Raghupathi R, Graham DI, Smith DH, Lee VM, Trojanowski JQ: The Dorothy Russell Memorial Lecture. The molecular and cellular sequelae of experimental traumatic brain injury: pathogenetic mechanisms. Neuropathol Appl Neurobiol. 1998, 24 (4): 251-267. 10.1046/j.1365-2990.1998.00121.x.
Bramlett HM, Dietrich WD: Pathophysiology of cerebral ischemia and brain trauma: similarities and differences. J Cereb Blood Flow Metab. 2004, 24 (2): 133-150.
Loane DJ, Faden AI: Neuroprotection for traumatic brain injury: translational challenges and emerging therapeutic strategies. Trends Pharmacol Sci. 2010, 31 (12): 596-604. 10.1016/j.tips.2010.09.005.
Ingenuity® Systems Web Site. http://www.ingenuity.com,
Morganti-Kossmann MC, Rancan M, Stahel PF, Kossmann T: Inflammatory response in acute traumatic brain injury: a double-edged sword. Curr Opin Crit Care. 2002, 8 (2): 101-105. 10.1097/00075198-200204000-00002.
Ghirnikar RS, Lee YL, Eng LF: Inflammation in traumatic brain injury: role of cytokines and chemokines. Neurochem Res. 1998, 23 (3): 329-340. 10.1023/A:1022453332560.
IPA® Canonical Pathways Overview. http://ingenuity.force.com/ipa/IPATutorials?id=kA250000000TN3aCAG,
Fan L, Young PR, Barone FC, Feuerstein GZ, Smith DH, McIntosh TK: Experimental brain injury induces expression of interleukin-1[beta] mRNA in the rat brain. Brain Res Mol Brain Res. 1995, 30 (1): 125-130. 10.1016/0169-328X(94)00287-O.
Rostworowski M, Balasingam V, Chabot S, Owens T, Yong VW: Astrogliosis in the neonatal and adult murine brain post-trauma: elevation of inflammatory cytokines and the lack of requirement for endogenous interferon-gamma. J Neurosci. 1997, 17 (10): 3664-3674.
Shojo H, Kaneko Y, Mabuchi T, Kibayashi K, Adachi N, Borlongan CV: Genetic and histologic evidence implicates role of inflammation in traumatic brain injury-induced apoptosis in the rat cerebral cortex following moderate fluid percussion injury. Neuroscience. 2010, 171 (4): 1273-1282. 10.1016/j.neuroscience.2010.10.018.
Hang CH, Shi JX, Tian J, Li JS, Wu W, Yin HX: Effect of systemic LPS injection on cortical NF-kappaB activity and inflammatory response following traumatic brain injury in rats. Brain Res. 2004, 1026 (1): 23-32. 10.1016/j.brainres.2004.07.090.
IPA® Network Generation. http://ingenuity.force.com/ipa/IPATutorials?id=kA250000000TNBZCA4,
Ford G, Xu Z, Gates A, Jiang J, Ford BD: Expression Analysis Systematic Explorer (EASE) analysis reveals differential gene expression in permanent and transient focal stroke rat models. Brain Res. 2006, 1071 (1): 226-236. 10.1016/j.brainres.2005.11.090.
Xu Z, Ford GD, Croslan DR, Jiang J, Gates A, Allen R, Ford BD: Neuroprotection by neuregulin-1 following focal stroke is associated with the attenuation of ischemia-induced pro-inflammatory and stress gene expression. Neurobiol Dis. 2005, 19 (3): 461-470. 10.1016/j.nbd.2005.01.027.
Xu Z, Jiang J, Ford G, Ford BD: Neuregulin-1 is neuroprotective and attenuates inflammatory responses induced by ischemic stroke. Biochem Biophys Res Commun. 2004, 322 (2): 440-446. 10.1016/j.bbrc.2004.07.149.
Li Y, Lein PJ, Liu C, Bruun DA, Tewolde T, Ford G, Ford BD: Spatiotemporal pattern of neuronal injury induced by DFP in rats: a model for delayed neuronal cell death following acute OP intoxication. Toxicol Appl Pharmacol. 2011, 253 (3): 261-269. 10.1016/j.taap.2011.03.026.
Israelsson C, Bengtsson H, Kylberg A, Kullander K, Lewen A, Hillered L, Ebendal T: Distinct cellular patterns of upregulated chemokine expression supporting a prominent inflammatory role in traumatic brain injury. J Neurotrauma. 2008, 25 (8): 959-974. 10.1089/neu.2008.0562.
Chavarria A, Alcocer-Varela J: Is damage in central nervous system due to inflammation?. Autoimmun Rev. 2004, 3 (4): 251-260. 10.1016/j.autrev.2003.09.006.
Laird MD, Vender JR, Dhandapani KM: Opposing roles for reactive astrocytes following traumatic brain injury. Neurosignals. 2008, 16 (2–3): 154-164.
Di Pietro V, Amin D, Pernagallo S, Lazzarino G, Tavazzi B, Vagnozzi R, Pringle A, Belli A: Transcriptomics of traumatic brain injury: gene expression and molecular pathways of different grades of insult in a rat organotypic hippocampal culture model. J Neurotrauma. 2010, 27 (2): 349-359. 10.1089/neu.2009.1095.
Herx LM, Rivest S, Yong VW: Central nervous system-initiated inflammation and neurotrophism in trauma: IL-1 beta is required for the production of ciliary neurotrophic factor. J Immunol. 2000, 165 (4): 2232-2239.
Kochanek PM, Berger RP, Bayir H, Wagner AK, Jenkins LW, Clark RS: Biomarkers of primary and evolving damage in traumatic and ischemic brain injury: diagnosis, prognosis, probing mechanisms, and therapeutic decision making. Curr Opin Crit Care. 2008, 14 (2): 135-141. 10.1097/MCC.0b013e3282f57564.
Khan M, Im YB, Shunmugavel A, Gilg AG, Dhindsa RK, Singh AK, Singh I: Administration of S-nitrosoglutathione after traumatic brain injury protects the neurovascular unit and reduces secondary injury in a rat model of controlled cortical impact. J Neuroinflammation. 2009, 6: 32-10.1186/1742-2094-6-32.
Wang Q, Tang XN, Yenari MA: The inflammatory response in stroke. J Neuroimmunol. 2007, 184 (1–2): 53-68.
Kong Y, Le Y: Toll-like receptors in inflammation of the central nervous system. Int Immunopharmacol. 2011, 11 (10): 1407-1414. 10.1016/j.intimp.2011.04.025.
Rall JM, Matzilevich DA, Dash PK: Comparative analysis of mRNA levels in the frontal cortex and the hippocampus in the basal state and in response to experimental brain injury. Neuropathol Appl Neurobiol. 2003, 29 (2): 118-131. 10.1046/j.1365-2990.2003.00439.x.
Fahlenkamp AV, Coburn M, Czaplik M, Ryang YM, Kipp M, Rossaint R, Beyer C: Expression analysis of the early chemokine response 4 h after in vitro traumatic brain injury. Inflamm Res. 2011, 60 (4): 379-387. 10.1007/s00011-010-0281-6.
von Gertten C, Flores MA, Holmin S, Mathiesen T, Nordqvist AC: Genomic responses in rat cerebral cortex after traumatic brain injury. BMC Neurosci. 2005, 6: 69-10.1186/1471-2202-6-69.
Kobori N, Clifton GL, Dash P: Altered expression of novel genes in the cerebral cortex following experimental brain injury. Brain Res Mol Brain Res. 2002, 104 (2): 148-158. 10.1016/S0169-328X(02)00331-5.
Poulsen CB, Penkowa M, Borup R, Nielsen FC, Caceres M, Quintana A, Molinero A, Carrasco J, Giralt M, Hidalgo J: Brain response to traumatic brain injury in wild-type and interleukin-6 knockout mice: a microarray analysis. J Neurochem. 2005, 92 (2): 417-432. 10.1111/j.1471-4159.2004.02877.x.
Ley EJ, Clond MA, Singer MB, Shouhed D, Salim A: IL6 deficiency affects function after traumatic brain injury. J Surg Res. 2011, 170 (2): 253-256. 10.1016/j.jss.2011.03.006.
Clausen F, Hanell A, Bjork M, Hillered L, Mir AK, Gram H, Marklund N: Neutralization of interleukin-1beta modifies the inflammatory response and improves histological and cognitive outcome following traumatic brain injury in mice. Eur J Neurosci. 2009, 30 (3): 385-396. 10.1111/j.1460-9568.2009.06820.x.
Feuerstein GZ, Wang X, Barone FC: Inflammatory gene expression in cerebral ischemia and trauma. Potential new therapeutic targets. Ann N Y Acad Sci. 1997, 825: 179-193. 10.1111/j.1749-6632.1997.tb48428.x.
Li HH, Lee SM, Cai Y, Sutton RL, Hovda DA: Differential gene expression in hippocampus following experimental brain trauma reveals distinct features of moderate and severe injuries. J Neurotrauma. 2004, 21 (9): 1141-1153.
Affymetrix: GeneChip® 3′IVT Express Kit User Manual.vol. P/N 702646 Rev. 8. 2010, http://www.affymetrix.com/support/technical/manuals.affx,
Sturn A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of microarray data. Bioinformatics. 2002, 18 (1): 207-208. 10.1093/bioinformatics/18.1.207.
The authors would like to thank Brock Wester, Kaelin Brewster, and Samantha Simon for their technical assistance. This work was supported by National Institutes of Health (NIH) grants U01 NS 057993 (BDF), U54 NS060659 (BDF), Department of Defense Contract #W81XWH-10-2-0055 (BDF), the W.M. Keck Foundation (BDF); BRIC NIH grant #5P20M006131-02 (GDF), Howard Hughes Medical Institute grant #52006306 (GDF) and PHS Grant (UL1 RR025008, KL2 RR025009 or TL1 RR025010) from the Clinical and Translational Science Award program, NIH, National Center for Research Resources (NCRR) (MCL). The project described was supported by Morehouse School of Medicine Grants Number U54 RR026137, G12RR003034 and S21MD000101 from the NCRR, a component of NIH, and its contents are solely the responsibility of the authors and do not necessarily represent the official views of NCRR or NIH.
The authors declare that they have no competing interests.
TEW, GDF and BDF were responsible for overall study design and execution. TEW, MCS and MCL were responsible for the animal models, sample preparation and histological analysis. TEW, GDF, ASG and BDF were responsible for carrying out microarray studies and bioinformatic data analysis. The manuscript was written by TEW and BDF. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: GOI (genes of interest) network created in IPA. This is the resulting network when IPA connected our 114 GOI using only direct (1st order) connections between the genes. 95 of the GOI formed an interconnected network, leaving 19 “orphan” genes. (PDF 7 MB)
Additional file 2: Primary tier example for how direct connections were counted. This figure shows an example of how we calculated the number of direct connections for a gene in our GOI network. In IPA, the gene in question was selected (JAK2 in this example). Then, its direct connections were selected by right clicking on JAK2 and using the "select nearest neighbors" option (highlighted in blue). A list of the selected genes was exported and JAK2 was removed from the list (upper right corner). The remaining genes were counted (15 in this example) and JAK2 was ranked in the gene interaction hierarchy (primary tier) by this number. (PDF 5 MB)
Additional file 3: Secondary tier example for how direct connections were counted. This figure shows an example of how we calculated the number of direct connections for a gene in our GOI network. In IPA, the gene in question was selected (IRF2 in this example). Then, its direct connections were selected by right clicking on IRF2 and using the "select nearest neighbors" option (highlighted in blue). A list of the selected genes was exported and IRF2 was removed from the list (upper right corner). The remaining genes were counted (7 in this example) and IRF2 was ranked in the gene interaction hierarchy (secondary tier) by this number. (PDF 5 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
White, T.E., Ford, G.D., Surles-Zeigler, M.C. et al. Gene expression patterns following unilateral traumatic brain injury reveals a local pro-inflammatory and remote anti-inflammatory response. BMC Genomics 14, 282 (2013). https://doi.org/10.1186/1471-2164-14-282
- Gene interaction hierarchy
- Traumatic brain injury