IMR Press / JIN / Volume 22 / Issue 2 / DOI: 10.31083/j.jin2202044
Open Access Original Research
Identification of Key Genes and Pathways in the Hippocampus after Traumatic Brain Injury: Bioinformatics Analysis and Experimental Validation
Hai Zeng1,†Sisi Zhao1,†Zhao Pang2,†Shiqi Wang1Luxi Cao1,*Yimin Zhang1,*
Show Less
1 Department of Acupuncture and Moxibustion, College of Traditional Chinese Medicine, Jinan University, 510632 Guangzhou, Guangdong, China
2 Medical Administration Division, The First Affiliated Hospital of Jinan University, 510630 Guangzhou, Guangdong, China
*Correspondence: (Luxi Cao); (Yimin Zhang)
These authors contributed equally.
J. Integr. Neurosci. 2023, 22(2), 44;
Submitted: 25 June 2022 | Revised: 27 September 2022 | Accepted: 30 September 2022 | Published: 20 February 2023
Copyright: © 2023 The Author(s). Published by IMR Press.
This is an open access article under the CC BY 4.0 license.

Background: Traumatic brain injury (TBI) is a common brain injury with a high morbidity and mortality. The complex injury cascade triggered by TBI can result in permanent neurological dysfunction such as cognitive impairment. In order to provide new insights for elucidating the underlying molecular mechanisms of TBI, this study systematically analyzed the transcriptome data of the rat hippocampus in the subacute phase of TBI. Methods: Two datasets (GSE111452 and GSE173975) were downloaded from the Gene Expression Omnibus (GEO) database. Systematic bioinformatics analyses were performed, including differentially expressed genes (DEGs) analysis, gene set enrichment analysis (GSEA), Gene Ontology (GO) enrichment analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, protein-protein interaction (PPI) network construction, and hub gene identification. In addition, hematoxylin and eosin (HE), Nissl, and immunohistochemical staining were performed to assess the injured hippocampus in a TBI rat model. The hub genes identified by bioinformatics analyses were verified at the mRNA expression level. Results: A total of 56 DEGs were shared in the two datasets. GSEA results suggested significant enrichment in the MAPK and PI3K/Akt pathways, focal adhesion, and cellular senescence. GO and KEGG analyses showed that the common DEGs were predominantly related to immune and inflammatory processes, including antigen processing and presentation, leukocyte-mediated immunity, adaptive immune response, lymphocyte-mediated immunity, phagosome, lysosome, and complement and coagulation cascades. A PPI network of the common DEGs was constructed, and 15 hub genes were identified. In the shared DEGs, we identified two transcription co-factors and 15 immune-related genes. The results of GO analysis indicated that these immune-related DEGs were mainly enriched in biological processes associated with the activation of multiple cells such as microglia, astrocytes, and macrophages. HE and Nissl staining results demonstrated overt hippocampal neuronal damage. Immunohistochemical staining revealed a marked increase in the number of Iba1-positive cells in the injured hippocampus. The mRNA expression levels of the hub genes were consistent with the transcriptome data. Conclusions: This study highlighted the potential pathological processes in TBI-related hippocampal impairment. The crucial genes identified in this study may serve as novel biomarkers and therapeutic targets, accelerating the pace of developing effective treatments for TBI-related hippocampal impairment.

traumatic brain injury
bioinformatics analysis
hub gene
1. Introduction

Traumatic brain injury (TBI) is a lethal and disabling injury and is a major challenge to public health and social development worldwide [1]. The global incidence of TBI has been continuously increasing over the past several decades [1, 2]. In 2016, there were more than 27 million new cases of TBI, and the total number of prevalent cases exceeded 55 million [2]. TBI is caused by the application of an external mechanical force to the head, followed by a series of primary and secondary pathological mechanisms, such as cell death in the brain parenchyma, disruption of the neurovascular unit, axonal injury, neuroinflammation, and neurodegeneration [3]. Due to the complex pathological mechanisms of TBI, developing effective neuroprotective therapeutics for afflicted patients has been challenging [3]. In recent years, increasing emphasis has been placed on the progressive and chronic consequences of TBI. Many preclinical and clinical studies have suggested that TBI is a significant risk factor for neuropsychiatric and cognitive disorders such as posttraumatic stress disorder (PTSD), chronic traumatic encephalopathy (CTE), and Alzheimer’s disease (AD) [4, 5, 6, 7]. The hippocampus is the main locus responsible for cognitive processing involving learning and memory [6, 8]. It has been shown that hippocampus dysfunction is associated with cognitive deficits following a TBI [6, 8, 9]. Existing treatments for TBI primarily focus on immediate neurosurgical procedures and long-term behavioral rehabilitation, but few treatments are sufficient to improve TBI-induced cognitive deficits [10, 11]. TBI differs from other diseases in that it involves different periods after the injury. However, subacute phase studies of TBI functional recovery have received substantially less attention than acute phase studies [12, 13]. Since effective therapeutic strategies for the subacute phase of TBI are limited, it is essential to explore optimal therapeutic targets at this stage [12, 13].

Understanding the molecular mechanisms underlying pathological processes is a critical step toward developing effective treatments for TBI. Based on the systematic bioinformatics analysis, microarray and high throughput sequencing techniques have emerged as essential tools for exploring the pathogenesis of complex human diseases [14, 15]. Many bioinformatics studies have been conducted in the field of oncology and some non-tumor diseases, but relatively few have focused explicitly on traumatic nerve injuries such as TBI [16]. The Gene Expression Omnibus (GEO) is a publicly available database consisting of an invaluable resource of mass gene expression data that can be systematically analyzed for biomarker or therapeutic target discovery [17]. In this study, we used two datasets from the GEO database to screen the differentially expressed genes (DEGs) of the hippocampus in the subacute phase of TBI. Gene functional enrichment analyses, protein-protein interaction (PPI) network construction, and hub gene identification were also performed. Furthermore, we observed the main histopathological features of the injured hippocampus and selected several DEGs for verification in vivo. The schematic workflow of our study is presented in Fig. 1. This study may provide new insight into potential biomarkers and therapeutic targets for cognitive deficits induced by TBI.

Fig. 1.

Schematic illustration of bioinformatics analyses and experimental validation in the study.

2. Materials and Methods
2.1 Dataset Source

We downloaded two mRNA expression datasets (GSE111452 and GSE173975) from the GEO database ( [18, 19]. In two experiments, adult male Sprague-Dawley (SD) rats were used as the experimental models. The rats were anesthetized by inhaling 1.5–3.0% isoflurane; then, a circular craniotomy was made on the right side of the parietal bone. The rats in the model group were subjected to a fluid percussion injury (FPI). The fluid pressure pulses were 2.3 and 2.12 atm in the GSE111452 and GSE173975 datasets. The rats in the sham group underwent a similar surgical preparation except for the FPI. The rat hippocampus tissues on the side of the lesion were harvested on the 14th day (subacute phase) after TBI modeling. Each dataset contained four separate hippocampus tissue samples in the sham and TBI groups. Detailed information of the two datasets can be found in Supplementary Table 1.

2.2 DEGs Identification

We utilized the “limma” and “edgeR” packages to identify DEGs between the sham and TBI groups in the GSE111452 and GSE173975 datasets, respectively. The criterion for defining differential gene expression was a p-value < 0.05 and a |log2 (fold change)|>1.00. The common DEGs in the two datasets were identified using the “VennDiagram” package.

2.3 Gene Enrichment Analysis

Gene set enrichment analysis (GSEA), Gene Ontology (GO) enrichment analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted by the “clusterProfiler” package [20]. For GSEA, a |standardized enrichment score (NES)|>1.00 and an adjusted p-value < 0.05 were considered as the criteria for significant enrichment. For the GO and KEGG analyses, the significant threshold was set as an adjusted p-value < 0.05. The GO enrichment analysis contained three categories: biological process (BP), cellular components (CC), and molecular function (MF).

2.4 Construction of a PPI Network and Identification of Hub Genes

A PPI network was drawn by the Search Tool for the Retrieval of Interacting Genes (STRING) database ( The threshold for the minimum interaction score was set to 0.4. The PPI network was plotted in Cytoscape. The hub genes were identified using the MCC algorithm of the CytoHubba plug-in.

2.5 Transcription Factor (TF) and Co-Factor Analysis

The list of TF and co-factor of rats were downloaded from the AnimalTFDB database [21]. The differentially expressed TF and co-factors shared in common between the two datasets were identified using the “VennDiagram” package.

2.6 Immune-Related Genes (IRGs) Analysis

The list of IRGs was obtained from the InnateDB database [22, 23]. The overlapped differentially expressed IRGs in the two datasets were identified using the “VennDiagram” package.

2.7 Animal Experiment

Specific pathogen-free (SPF) Male SD rats (7 weeks old, 220–240 g body weight) were purchased from Beijing Huafukang Bioscience (Beijing, China). All animals were housed in a standard environment (12-h light/dark cycle, room temperature: 23–25 °C, relative humidity: 50%–60%) with water and food ad libitum. After adaptive feeding for ten days, fourteen rats were randomly allocated into the sham and TBI groups (n = seven per group). The experimental procedures were approved by the Laboratory Animal Ethics Committee of Jinan University (ethical number: 14075) and followed the Principles of Laboratory Animal Care.

2.7.1 TBI Model

TBI model was produced with a controlled cortical impact (CCI) device (YHC199, Wuhan Yihong Science & Technology Co., Ltd., Wuhan, Hubei, China). Briefly, rats were anesthetized with 2% pentobarbital sodium (40 mg/kg) with an intraperitoneal injection. Then, they were fixed in a stereotaxic frame, and a skin incision was made on the scalp to expose the skull. A craniotomy (diameter of 5.0 mm) was performed over the right parietal bone (2.0 mm lateral to the sagittal suture and 1.0 mm posterior to the bregma) to expose the dura and cerebral cortex. A flat-tipped impactor (diameter of 4 mm) was utilized to impact the exposed dura (impact parameters: depth of 2.5 mm, dwell time of 250 ms, and impact velocity of 3.5 m/s) [24]. After surgery, scalp incisions were sutured, and the rats were moved into clean cages with heated pads for recovery. The rats in the sham group were subjected to the same procedures without the CCI.

2.7.2 Tissue Collection and Preparation

To collect the hippocampal tissue, on the 14th day after TBI, four rats in each group were anesthetized with 2% pentobarbital sodium (40 mg/kg) with an intraperitoneal injection, followed by transcardial perfusion with ice-cold 0.9% normal saline. Then, the brain was removed, and the hippocampal tissue in the ipsilateral (lesion) side was rapidly removed. The collected samples were flash-frozen in liquid nitrogen and stored at –80 °C for subsequent analysis. For histopathological detection, three rats in each group were anesthetized in the same manner; then, the brain was harvested after transcardial perfusion with ice-cold 0.9% normal saline followed by 4% paraformaldehyde. The tissue blocks were embedded in paraffin after fixation in 4% paraformaldehyde for 24 hours.

2.8 Hematoxylin and Eosin (HE) Staining

HE staining was conducted using a HE staining kit (Servicebio, Wuhan, Hubei, China) according to the instruction manual. Brain sections (5 um) were subjected to deparaffinization with dimethylbenzene, ethanol hydration, hematoxylin staining, eosin staining, and ethanol dehydration. Finally, the sections were sealed with neutral resin. Tissue sections were observed by a light microscope (Olympus, Tokyo, Japan).

2.9 Nissl Staining

Nissl staining was performed using Nissl staining solution (Servicebio, Wuhan, Hubei, China) according to the manufacturer’s instructions. After deparaffinization and hydration, brain slices were stained in Nissl staining solution and then differentiated in 0.1% glacial acetic acid. After sealing the slices, the staining was visualized under a light microscope (Olympus, Tokyo, Japan).

2.10 Immunohistochemical Staining

After deparaffinization and hydration, brain sections underwent high-temperature antigen repair using citrate buffer (pH 6.0). The sections were cooled naturally after antigen retrieval. To block endogenous peroxidases, 3% hydrogen peroxide solution was applied to the sections for 25 minutes, and nonspecific bindings were blocked with 3% bovine serum albumin (BSA) solution for 30 minutes. Then, the sections were incubated with primary antibody mouse anti-Iba1 (1:500, Servicebio, Wuhan, Hubei, China) overnight at 4 °C. After removing the primary antibody, the sections were incubated in anti-mouse HRP-labeled secondary antibody (1:500, Servicebio, Wuhan, Hubei, China) at room temperature for one hour. Diaminobenzidine (DAB) staining was performed using the DAB chromogenic agent (Servicebio, Wuhan, Hubei, China). The sections were sealed with coverslips after nuclear counterstaining and dehydration. The positive staining cells were observed using a light microscope (Olympus, Tokyo, Japan).

2.11 Quantitative Real-Time PCR (RT-qPCR)

Total RNA was extracted from the hippocampal tissues using RNAiso Plus (Takara, Ohtsu, Japan). A NanoDrop spectrophotometer (Thermo Scientific, Waltham, MA, USA) was used to measure RNA concentration and purity. RT-qPCR was performed using the SYBR Green RT-qPCR kit (Takara, Ohtsu, Japan) and the CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). We used β-actin as a reference gene and calculated the relative mRNA expression level by the 2-ΔΔCt method. Reactions were performed in triplicate. The sequences of primers used are provided in Supplementary Table 2.

2.12 Statistical Analysis

All data were presented as mean ± standard error of the mean (SEM). Two groups of data were compared using Student’s t-test. Statistical analyses were performed with SPSS (version 23.0, IBM Corp., Chicago, IL, USA) software. Plots were generated using GraphPad Prism (version 6.0, GraphPad Software Inc., San Diego, CA, USA) and R (version 4.1.1, the R Foundation, Vienna, Austria) softwares. A p-value < 0.05 indicated statistical significance.

3. Results
3.1 Bioinformatics Analyses of the GEO Datasets
3.1.1 Identification of DEGs

A total of 445 genes were differentially expressed between the sham and TBI groups in the GSE111452 dataset, including 276 upregulated genes and 169 downregulated genes (Fig. 2A,C). In the GSE173975 dataset, 802 DEGs were screened between the sham and TBI groups, including 718 upregulated genes and 84 downregulated genes (Fig. 2B,D). The Venn diagram was utilized to identify the intersections of DEGs from the GSE111452 and GSE173975 datasets. Overall, 56 upregulated DEGs were shared in common between the two datasets.

Fig. 2.

The DEGs between the sham and TBI groups in the GSE111452 and GSE173975 datasets. (A) Volcano plot of the DEGs in the GSE111452 dataset. (B) Volcano plot of the DEGs in the GSE173975 dataset. Red dots represent upregulated genes, and blue dots represent downregulated genes. (C) Heatmap of the top 50 DEGs in the GSE111452 dataset. (D) Heatmap of the top 50 DEGs in the GSE173975 dataset.

3.1.2 GSEA Analysis

To avoid missing some genes with crucial biological relevance during the screening of DEGs, GSEA was conducted based on all gene expression data of the two datasets. We extracted four main pathways from the overlapped pathways in the two datasets, including the mitogen-activated protein kinase (MAPK) signaling pathway, phosphatidylinositol 3-kinases/protein kinase B (PI3K/Akt) signaling pathway, focal adhesion, and cellular senescence. The enrichment plots of the GSEA results are shown in Fig. 3.

Fig. 3.

The main results of GSEA in the GSE111452 (A) and GSE173975 (B) datasets.

3.1.3 GO and KEGG Analyses of DEGs

For the BP terms, the DEGs were mainly involved in antigen processing and presentation, leukocyte-mediated immunity, adaptive immune response, and lymphocyte-mediated immunity.

For the CC terms, the DEGs were primarily enriched in lytic vacuole, lysosome, late endosome, external side of plasma membrane, phagocytic vesicle, lysosomal membrane, lytic vacuole membrane, MHC class II protein complex, collagen trimer, and vacuolar membrane.

For the MF terms, the DEGs were significantly enriched in immune receptor activity, immunoglobulin receptor activity, peptide binding, cysteine-type endopeptidase activity, amide binding, immunoglobulin binding, amyloid-beta binding, cysteine-type peptidase activity, peptide antigen binding, and endopeptidase activity.

KEGG pathway enrichment analysis found that the DEGs were mainly involved in antigen processing and presentation, phagosome, lysosome, and complement and coagulation cascades, etc. The results of the GO and KEGG analyses are shown in Figs. 4,5.

Fig. 4.

GO analyses of the common DEGs. (A) The bar plot of the top 10 BP, CC, and MF terms. (B) The bubble plot of the top 10 BP, CC, and MF terms.

Fig. 5.

KEGG analyses of the common DEGs. (A) The bar plot of the top 10 KEGG pathways. (B) The bubble plot of the top 10 KEGG pathways.

3.1.4 PPI Analysis and Hub Gene Identification

The PPI network was constructed to reveal the potential connection among the DEGs. As shown in Fig. 6A, the PPI network consisted of 52 nodes and 220 linkages. We further identified the top 15 hub genes according to the MCC method. Next, we performed K-means clustering and grouped these hub genes into four distinct clusters as follows: cluster 1: allograft inflammatory factor 1 (Aif1), colony-stimulating factor 1 receptor (Csf1r), triggering receptor expressed on myeloid cells 2 (Trem2), and TYRO protein tyrosine kinase-binding protein (Tyrobp); cluster 2: C1qa, C1qb, C1qc, Cd74, and Cathepsin S (Ctss); cluster 3: Rho GDP dissociation inhibitor-β (Arhgdib), Cd53, lysosomal-associated protein transmembrane 5 (Laptm5); and cluster 4: Fc epsilon receptor Ig (Fcer1g), Fc gamma receptor IIb (Fcgr2b), and protein tyrosine phosphatase non-receptor type 6 (Ptpn6) (Fig. 6B).

Fig. 6.

Construction of the PPI networks of the DEGs (A) and the top 15 hub genes (B).

3.1.5 TF, TF Co-Factors, and IRGs Analysis

We identified two overlapped TF co-factors among the DEGs shared in the two datasets (Fig. 7A). They were adipocyte enhancer-binding protein 1 (Aebp1) and PYD and CARD domain containing (Pycard) genes. Furthermore, we also found 15 immune-related DEGs, including C1qa, interferon γ receptor 1 (Ifngr1), Csf1r, Unc-93 homolog B1 (Unc93b1), Ctss, C1qc, Ptpn6, C1qb, legumain (Lgmn), Tyrobp, Trem2, granulin precursor (Grn), annexin A2 (Anxa2), interferon-stimulated gene 15 (Isg15), and Pycard. The GO enrichment analysis of these immune-related DEGs indicated that the significant BP terms included microglial cell activation, cell junction disassembly, astrocyte activation, antigen processing and presentation, macrophage activation, positive regulation of response to external stimulus, myeloid cell differentiation, myeloid leukocyte activation, and astrocyte development (Fig. 7B).

Fig. 7.

Identification of differentially expressed TF, TF co-factors, and IRGs. (A) The Venn diagram presenting overlapped genes among the common DEGs, TF, TF co-factors, and IRGs. (B) The circle plot of the top 10 BP terms of the immune-related DEGs.

3.2 Experimental Validation in a Rat Model of TBI
3.2.1 Histopathological Features

HE staining was used to observe the cell morphology of neurons. In the sham group, neurons in the hippocampus were arranged regularly and showed clear and intact cellular structure (Fig. 8A). However, the TBI group showed distinct morphological alterations, nuclear pyknosis, and neuronal necrosis (Fig. 8A). Nissl staining was used to detect neuronal loss. The results of Nissl staining also demonstrated severe structural damage in the neurons and apparent reductions in the number of intact neurons in the TBI group (Fig. 8B). No such pathological abnormalities were observed in the sham group (Fig. 8B). Moreover, microglia was visualized by staining with anti-Iba1 antibody. The results of Iba1 staining suggested that the number of Iba1-positive cells was markedly increased in the TBI group (Fig. 8C).

Fig. 8.

Representative histopathological images of the hippocampus from the sham and TBI groups (scale bar = 50 um). (A) HE staining. (B) Nissl staining. (C) Immunohistochemical staining for Iba1. Arrows point to Iba1-positive cells.

3.2.2 The mRNA Expression Levels of the Hub Genes

To validate the transcriptome data, we selected the 15 hub genes (Aif1, Csf1r, Trem2, Tyrobp, C1qa, C1qb, C1qc, Cd74, Ctss, Arhgdib, Cd53, Laptm5, Fcer1g, Fcgr2b, and Ptpn6) for RT-qPCR (Fig. 9). Compared with the sham group, the mRNA expression levels of these genes were significantly higher in the TBI group (p < 0.05). Overall, these RT-qPCR results were in accordance with the transcriptome data.

Fig. 9.

Validation of the 15 hub genes at the mRNA level between the sham and TBI groups. n = four per group, data were presented as mean ± SEM. *p < 0.05, **p < 0.01, and ***p < 0.001.

4. Discussion

TBI is one of the most common traumatic injuries which can cause long-lasting cerebral damage and cognitive deficits [25]. Elucidating the potential molecular mechanisms and therapeutic targets for improving cognitive function helps to improve the clinical prognosis and alleviate the burden on patients suffering from TBI [26]. Transcriptomic profiling is a powerful tool that can identify the genes differentially expressed in specific physiological and pathological states and discover new diagnostic or therapeutic targets [27]. In this study, we conducted a systematic bioinformatics analysis of the two GEO datasets (GSE111452 and GSE173975) to reveal the potential molecular mechanisms responsible for hippocampal damage in rats in the subacute phase of TBI.

Initially, we performed GSEA on the basis of all transcriptome data of the two GEO datasets, separately. The GSEA results revealed notable enrichment in the MAPK signaling pathway, PI3K/Akt signaling pathway, focal adhesion, and cellular senescence. The MAPK family in mammals consists of c-Jun N-terminal kinase, p38 MAPK, and extracellular signal-regulated kinase [28]. Previous findings suggested that inhibiting the activation of the MAPK signaling pathway might produce neuroprotective effects in TBI [29, 30]. As a critical regulator of organismal growth and multiple cellular processes, the PI3K/Akt signaling pathway plays an important role under various pathophysiological conditions [31]. Preclinical studies of TBI revealed that the PI3K/Akt signaling pathway exerted positive effects on neural repair by regulating neuronal apoptosis and neurogenesis [32, 33]. Focal adhesions are protein complexes that mediate cell adhesion by linking the extracellular matrix to the cytoskeleton [34]. A recent in vitro study demonstrated that inhibiting focal adhesion kinase phosphorylation could mitigate astrocyte activation induced by mechano-stimulation [35]. Cellular senescence refers to a cell state triggered by various stress factors and some specific physiological conditions [36]. It is characterized by proliferative arrest, secretory phenotype, macromolecular damage, and alternation in metabolism [36, 37]. Several studies suggested that DNA damage-induced cellular senescence was a possible driver of TBI-associated sequelae such as cognitive impairment [38, 39, 40]. Among the common DEGs, several genes might play a role in the response to DNA damage. Du et al. [41] found that the activation of complement C3 contributed to retinal DNA damage and C3 deficiency alleviated alkylation-induced retinal degeneration in mice. Raso et al. [42] reported that Isg15 upregulation promoted DNA replication fork progression, leading to extensive DNA damage and chromosomal aberrations in a human osteosarcoma cell line. Xie et al. [43] found that downregulating the expression of CX3C chemokine receptor 1 (Cx3cr1) accelerated double-strand DNA damage in irradiated ovarian cancer cells. These experiments were primarily in the field of oncology, and the studies that investigated how these genes regulated DNA damage response to impact the progression of TBI are still limited.

To further explore potential biological functions and signaling pathways involved in hippocampal damage, we performed GO and KEGG analyses based on the common DEGs of the two datasets. We found that the DEGs were mainly involved in an array of critical immune and inflammatory processes, including antigen processing and presentation, leukocyte and lymphocyte mediated immunity, adaptive immune response, phagosome, lysosome, and complement and coagulation cascades. Following the initial insult triggered by TBI, peripheral leukocytes and lymphocytes were successively recruited to the injury site [44]. Antigen processing and presentation by antigen-presenting cells promote the transition between innate and adaptive immune responses [45]. Moderate immune responses promote the removal of cellular debris, repair, and regeneration in the central nervous system (CNS) post-TBI, but excessive immune responses result in reactions causing destructive neuroinflammation and maladaptive secondary injuries [3, 44]. In view of the enrichment of phagosome- and lysosome-related genes after TBI, it is reasonable to speculate that active phagocytosis plays a critical role in neural injury and repair [46]. Microglia act as important phagocytes in the brain and can rapidly remove cellular and myelin debris after brain injury [47]. Moderate microglial phagocytosis is beneficial for remodeling the brain microenvironment and alleviating secondary injury [47, 48], whereas excessive phagocytosis by microglia of neurons, synapses, or myelin may impede tissue and functional recovery [49, 50]. Current evidence indicates extensive crosstalk between the complement and coagulation systems, which have fundamental clinical implications for inflammation, immunity, and tissue damage [51]. A previous proteomics study suggested that complement and coagulation cascades were activated during the acute and subacute phases of TBI, which was consistent with our findings [52].

In addition to the functional enrichment and pathway analyses described above, we proceeded to construct the PPI network of the DEGs and then identified the top 15 hub genes which were upregulated in both datasets. Among the hub genes, we found several microglia signature genes such as Csf1r, Aif1, and Trem2. Csf1r is a crucial regulator of myeloid lineage cells and is also necessary for maintaining the cell viability of microglia [53]. Several experiments adopted the method of Csf1r inhibition to attenuate microglia-mediated neuroinflammation in rodent models of TBI [54, 55]. However, given that microglia plays a key role in maintaining homeostasis and normal function of the brain [56], inhibiting Csf1r may increase the risk of failure in clinical translation. In contrast, the infiltrating peripheral macrophages in the brain after TBI may be promising targets for developing potential treatment strategies [57]. Aif1 is a calcium-binding protein that regulates immune and inflammatory responses [58, 59]. Increased expression of Aif1 is an indication of microglia and macrophage activation in response to brain trauma [60, 61]. A previous experimental study identified Trem2 and Tyrobp as significant hub genes in TBI mice expressing human APOE [62]. Trem2 is an innate immune receptor mainly expressed on various tissue macrophages, such as microglia in the brain [63]. It transmits intracellular signals through binding with the adaptor protein Tyrobp, which in turn regulates many crucial biological processes, including phagocytosis, chemotaxis, inflammatory responses, and lipid metabolism [56, 63, 64]. It was reported that Trem2 deficiency alleviated the acute peripheral macrophage infiltration and attenuated chronic hippocampal atrophy and cognitive deficits in a murine TBI model [65]. Given the vital role of Trem2 in the progression of neurodegeneration, it was worthwhile to explore it as a potential therapeutic target in managing TBI-mediated neurodegenerative pathologies [56, 57, 65].

Ctss is a cysteine protease expressed by diverse immune cells that is responsible for cleaving certain extracellular matrix proteins and cell adhesion molecules to promote immune cell motility [66, 67]. In the acute phase after TBI, the upregulation of Ctss expression was mainly observed in microglia; inhibiting Ctss could reduce the level of inflammatory factors, alleviate brain edema, and improve neurobehavioral function [68]. The Cd74 molecule is the cognate receptor of macrophage migration inhibitory factor that is recognized to be a pleiotropic inflammatory cytokine [69, 70]. Tobin et al. [45] found that TBI-induced neurodegeneration depended on antigen processing and presentation that required Cd74. Cd74 deficiency could decrease peripheral lymphocyte activation and neurodegeneration following TBI, suggesting that targeting the switch between innate and adaptive immunity might be a potential therapeutic strategy [45]. C1qa, C1qb, and C1qc are three distinct subunits of complement component C1q, which is a crucial protein in the complement cascade within the innate immune system [71, 72]. One recent clinical study showed that elevated serum C1q levels were significantly correlated with traumatic severity and could serve as an independent prognostic factor for long-term outcomes after TBI [73]. After TBI, microglia and neurons were the primary source of C1q [74, 75]. The increased expression of C1q might lead to neuron loss and chronic neuroinflammation and correlate with sleep spindle loss and epileptic spikes. Blocking C1q was beneficial to counteract these effects, indicating that C1q might be a potential target for treating the devastating long-term outcomes of TBI [75].

We also discovered two upregulated Fc receptor genes. It was reported that Fcer1g might be a potential microglial biomarker related to aging and neurodegeneration [72, 76]. Among the Fc receptors, Fcgr2b could serve as one marker of the pro-inflammatory phenotype of microglia after TBI [77, 78, 79]. There was some evidence that Laptm5 had pro-inflammatory properties in macrophages [80]. In a chronic constriction injury in rats, Laptm5 was identified as a significant upregulated neuroinflammation-related gene that influenced neuropathic pain behavior [81]. However, due to the lack of relevant literature, it is unknown whether Laptm5 plays a role in the pathological processes of TBI.

TF co-factors play an essential role in diverse biological processes by interacting with TFs to suppress or activate gene transcription [21, 82]. Among the DEGs shared in common between the GSE111452 and GSE173975 datasets, we found two inflammation-related TF co-factors (Aebp1 and Pycard). Aebp1 has been shown to promote inflammatory processes in macrophages by activating the nuclear factor kappa B [83, 84]. Shijo et al. [85] found that Aebp1 was highly expressed in hippocampal neurons and glial cells in AD patients, implying its potential role in the progression of AD pathology. The Pycard gene encodes a critical adaptor protein in activating inflammasomes [86]. Recent research has shown that inflammasome proteins such as Pycard are elevated in the blood of TBI patients and are closely associated with injury severity and pathological outcomes [87]. Experimental inhibition of inflammasomes exerted anti-inflammatory and neuroprotective effects in rodent models of TBI, implying that developing inflammasome-targeting therapeutics might be a potential direction for future research [88, 89, 90].

CTE refers to a progressive Tau-dependent neurodegenerative disease that remains relatively understudied [91, 92]. In several animal models of TBI, the immunocytochemical signature for CTE has been presented, which is in agreement with clinical studies [4, 7]. Among the common DEGs identified in this study, several genes (e.g., Trem2, Tyrobp, and Cx3cr1) might be involved in the progression of tauopathy based on the current literature. Jiang et al. [93] found that Trem2 attenuated tau kinase activity by inhibiting neuroinflammation and Trem2 deficiency exacerbated tau pathology and neurodegeneration in P301S tau transgenic mice. Another study reported that Trem2 deletion enhanced tau dispersion through microglia exosomes [94]. In a tauopathy mouse model, silencing of Tyrobp reduced C1q levels and improved learning behavior and electrophysiological properties despite increasing tau phosphorylation [95]. Cx3cr1 was deemed as a potential target for treating tauopathy, since its deficiency impaired the phagocytosis and internalization of extracellular tau by microglia [96, 97].

Recent evidence suggested that endoplasmic reticulum (ER) stress resulted in extensive behavioral changes related to a CTE-like phenotype after TBI [91, 98]. Inhibiting ER stress ameliorated cognitive deficits by reducing tau hyperphosphorylation in rodent models of CTE [91, 98]. Several genes within the common DEGs might be involved in ER stress response. It was reported that Glycoprotein nonmetastatic melanoma protein B (Gpnmb) increased the expression of chaperone protein BiP to increase the protein-folding capacity of the ER and attenuate cell death caused by ER stress [99]. In a cerebral ischemia mouse model, the neuronal damage of Gpnmb-transgenic mice was significantly alleviated compared with wild-type mice [99]. Kam et al. [100] found that Fcgr2b mediated amyloid-β neurotoxicity by activating ER stress and caspase-12. It has been reported that the translocator protein (Tspo) is a potential biomarker in multiple inflammatory and neurodegenerative diseases [101, 102]. Loss of Tspo in hepatocytes caused free cholesterol accumulation, which subsequently induced ER stress [103]. More importantly, recent clinical studies suggested that PET imaging with tracers binding to Tspo could be applied to detect the neuroinflammatory changes in subjects with CTE [102, 104]. Interestingly, a multi-tracer PET Study showed that tau deposition and increased Tspo expression were co-localized in the cortex and hippocampus in a tau transgenic mouse model; however, the complex interaction between tauopathy and Tspo expression remains to be further studied [105]. Neuroimaging modalities could be conducted longitudinally at various time points on the same subject, showing an advantage in observing dynamic change during the progression of CTE [92]. In future studies, molecular imaging through PET technology such as Tspo-PET may be used as a novel tool to study CTE on preclinical models and help to elucidate the complex pathophysiological mechanisms leading from TBI to symptomatic CTE [102, 104]. In addition, Gouna and colleagues [106] demonstrated that Trem2 deficiency in microglia/macrophages resulted in ER stress due to impaired lipid droplet formation. The potential critical role of these genes in both tauopathy progression and ER stress implied that they might be of considerable value in the field of CTE modeling, although the specific molecular mechanisms remain to be determined in future studies. Tau transgenic rodent models combined with targeted gene deletion might be instrumental for the study of CTE modeling and help to reveal the potential pathophysiological processes of CTE.

In summary, our research was a comprehensive bioinformatics study analyzing the transcriptome data of rat hippocampus in the subacute phase of TBI. We used multiple tools for enrichment analysis to reveal the underlying biological mechanisms involved in TBI-related hippocampal impairment. The enrichment analysis showed that TBI might result in a significant dysregulation of multiple inflammation and immune processes. We further identified the hub genes, TF co-factors, and IRGs, which likely play a critical role in the progression of hippocampal impairment.

In the validation experiment, we found some apparent histopathologic changes in the hippocampus during the subacute phase of TBI, including neuronal necrosis and reductions in the number of intact neurons. A previous study supports the pathological changes described in our model [107]. Furthermore, we used immunohistochemical staining to observe the number of microglia. The staining results suggested that the level of microglia activation was markedly elevated in the hippocampus of rats in the TBI group when compared with the sham group. Through immunofluorescence staining analysis, Luo et al. [52] also found that the number of microglia was significantly increased in the hippocampus during the subacute phase of TBI, which was consistent with our results. Similarly, the immune-related DEGs identified by bioinformatic analyses were significantly enriched in BP terms associated with the activation of microglia and macrophages. Microglia is the major type of glial cells involved in neuroinflammation in the mammalian CNS [108]. The peripheral immune cells (especially macrophages) are also essential responders to tissue injuries in the CNS [3, 109]. These findings suggest that these inflammatory cells might be valuable targets in the therapy of hippocampal impairment in the subacute phase after TBI. Moreover, we also detected the mRNA expression levels of the 15 hub genes identified by bioinformatics analyses. Overall, the RT-qPCR results were consistent with the transcriptome data, indicating that the transcriptome data were reliable. Overall, traumatic homogeneity between the datasets and our validation experiment was relatively high. Both the datasets and our experiment used adult male SD rats as the experimental objects, and the location and extent of the brain injury were similar. In addition, the tissue type and disease stage were consistent.

However, there were several limitations of this study. First, it specifically focused on the hippocampus rather than other brain regions. The principal reason was that hippocampal impairment acts as a key player in the pathogenesis of TBI-induced cognitive decline [6, 8, 9], which was one of the most dominant clinical manifestations in the subacute phase [110, 111, 112]. Developing effective treatments to promote cognitive rehabilitation in TBI patients during the subacute phase is of high clinical importance [110]. The bioinformatics analyses and experimental validation performed in this study mainly focused on the subacute phase (14 days), which was important for assessing cognitive function in the rodent model of TBI [13, 113]. Notably, there was some evidence that several potential biomarkers were significantly upregulated in the acute phase. For instance, Saber et al. [65] found that both the protein and mRNA expressions of Trem2 were significantly increased in the ipsilateral cortex in mice on the third day after experimental TBI. The study by Luo and colleagues [52] revealed that the multiple proteins (e.g., C1qa, Pycard, Aif1, Ptpn6, and Fcer1g) in the ipsilateral hippocampus were markedly upregulated in a rat model of TBI on the third day. In addition, Wang et al. [54] reported that the gene expression levels of C1qa, Aif1, and Csf1r were prominently increased in the perilesional area on the fifth day post injury in TBI mice. Therefore, these biomarker genes might become important therapeutic targets during the early phase of TBI. Further study on these crucial genes might provide useful insight for treating TBI-induced cognitive deficits in fundamental, translational, and clinical research.

Second, in this study, only male rats were used in the experiments. Previous observations have suggested that there may be a link between sex and clinical outcomes after TBI [114, 115]. Some studies demonstrated that females had a lower risk of mortality and cognitive dysfunction than males, while others failed to find such advantages [114, 115]. Given the predominant role of microglia in posttraumatic neuroinflammation, it is worthwhile to determine whether sex differences exist in microglial activation after TBI [114]. Several studies have revealed that male rodents may exhibit more significant microglial activation than female rodents [114, 116, 117]. However, there has also been conflicting evidence about the influence of sex on microglial activation in rodent models of TBI [79, 114]. Regrettably, there is currently an apparent male bias in the field of TBI research, both in preclinical rodent studies and clinical trials [116, 118]. Thus, future studies are needed before drawing definite conclusions about the role of sex in TBI. Finally, traditional microarray and mRNA sequencing techniques mainly investigate the global gene expression patterns of the same tissues. Single-cell mRNA sequencing techniques may help to reveal the exact role of specific cell types in the complex pathophysiological mechanisms underlying TBI.

5. Conclusions

In conclusion, this study identified several critical genes and key signaling pathways involved in hippocampal impairment induced by TBI. Several genes involved in neuroinflammation (e.g., Trem2, C1q, and Pycard) might represent potential targets for treating TBI-induced cognitive deficits.


TBI, Traumatic brain injury; PTSD, posttraumatic stress disorder; CTE, chronic traumatic encephalopathy; AD, Alzheimer’s disease; GEO, Gene Expression Omnibus; DEGs, differentially expressed genes; PPI, protein-protein interaction; FPI, fluid percussion injury; SD, Sprague-Dawley; GSEA, gene set enrichment analysis; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; NES, standardized enrichment score; BP, biological process; CC, cellular components; MF, molecular function; STRING, Search Tool for the Retrieval of Interacting Genes; TF, Transcription factor; IRGs, immune-related genes; SPF, specific pathogen-free; CCI, controlled cortical impact; HE, hematoxylin and eosin; BSA, bovine serum albumin; DAB, Diaminobenzidine; RT-qPCR, quantitative Real-Time PCR; SEM, standard error of the mean; MAPK, mitogen-activated protein kinase; PI3K/Akt, phosphatidylinositol 3-kinases/protein kinase B; Ctss, Cathepsin S; Aif1, allograft inflammatory factor 1; Csf1r, colony-stimulating factor 1 receptor; Trem2, triggering receptor expressed on myeloid cells 2; Tyrobp, TYRO protein tyrosine kinase-binding protein; Arhgdib, Rho GDP dissociation inhibitor-β; Laptm5, lysosomal-associated protein transmembrane 5; Fcer1g, Fc epsilon receptor Ig; Fcgr2b, Fc gamma receptor IIb; Ptpn6, protein tyrosine phosphatase non-receptor type 6; Aebp1, adipocyte enhancer-binding protein 1; Pycard, PYD and CARD domain containing; Ifngr1, interferon γ receptor 1; Unc93b1, Unc-93 homolog B1; Lgmn, legumain; Grn, granulin precursor; Anxa2, annexin A2; Isg15, interferon-stimulated gene 15; Cx3cr1, CX3C chemokine receptor 1; CNS, central nervous system; ER, endoplasmic reticulum; Gpnmb, Glycoprotein nonmetastatic melanoma protein B; Tspo, translocator protein.

Availability of Data and Materials

The transcriptomic datasets are available in the GEO database

Author Contributions

HZ and YZ contributed to the study design. SZ, SW, and ZP participated in data collation and statistical analysis. HZ, SZ, SW, and ZP finished the experimental validation. HZ and LC wrote the manuscript. YZ and SW revised the manuscript. All authors read and approved the final manuscript.

Ethics Approval and Consent to Participate

The animal experiment was approved by the Laboratory Animal Ethics Committee of Jinan University (ethical number: 14075).


The authors appreciate the GEO database and contributors for providing their meaningful datasets.


This study was supported by the National Natural Science Foundation of China (nos. 81873362, 82174483, and 82205249), the Natural Science Foundation of Guangdong Province (nos. 2021A1515011219 and 2021A1515110146), China Postdoctoral Science Foundation (2022M710912), and the Administration of Traditional Chinese Medicine of Guangdong Province (no. 20231086).

Conflict of Interest

The authors declare no conflict of interest.

Maas AIR, Menon DK, Adelson PD, Andelic N, Bell MJ, Belli A, et al. Traumatic brain injury: integrated approaches to improve prevention, clinical care, and research. The Lancet. Neurology. 2017; 16: 987–1048.
GBD 2016 Traumatic Brain Injury and Spinal Cord Injury Collaborators. Global, regional, and national burden of traumatic brain injury and spinal cord injury, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. The Lancet. Neurology. 2019; 18: 56–87.
Corps KN, Roth TL, McGavern DB. Inflammation and neuroprotection in traumatic brain injury. JAMA Neurology. 2015; 72: 355–362.
Bramlett HM, Dietrich WD. Long-Term Consequences of Traumatic Brain Injury: Current Status of Potential Mechanisms of Injury and Neurological Outcomes. Journal of Neurotrauma. 2015; 32: 1834–1848.
Wu Z, Wang ZH, Liu X, Zhang Z, Gu X, Yu SP, et al. Traumatic brain injury triggers APP and Tau cleavage by delta-secretase, mediating Alzheimer’s disease pathology. Progress in Neurobiology. 2020; 185: 101730.
Arneson D, Zhang G, Ahn IS, Ying Z, Diamante G, Cely I, et al. Systems spatiotemporal dynamics of traumatic brain injury at single-cell resolution reveals humanin as a therapeutic target. Cellular and Molecular Life Sciences: CMLS. 2022; 79: 480.
Simon DW, McGeachy MJ, Bayır H, Clark RSB, Loane DJ, Kochanek PM. The far-reaching scope of neuroinflammation after traumatic brain injury. Nature Reviews. Neurology. 2017; 13: 171–191.
Keightley ML, Sinopoli KJ, Davis KD, Mikulis DJ, Wennberg R, Tartaglia MC, et al. Is there evidence for neurodegenerative change following traumatic brain injury in children and youth? A scoping review. Frontiers in Human Neuroscience. 2014; 8: 139.
Gao F, Hu M, Zhang J, Hashem J, Chen C. TDP-43 drives synaptic and cognitive deterioration following traumatic brain injury. Acta Neuropathologica. 2022; 144: 187–210.
Khodaei S, Avramescu S, Wang DS, Sheng H, Chan NK, Lecker I, et al. Inhibiting α5 Subunit-Containing γ-Aminobutyric Acid Type A Receptors Attenuates Cognitive Deficits After Traumatic Brain Injury. Critical Care Medicine. 2020; 48: 533–544.
Chou A, Krukowski K, Jopson T, Zhu PJ, Costa-Mattioli M, Walter P, et al. Inhibition of the integrated stress response reverses cognitive deficits after traumatic brain injury. Proceedings of the National Academy of Sciences of the United States of America. 2017; 114: E6420–E6426.
Marklund N, Bellander BM, Godbolt AK, Levin H, McCrory P, Thelin EP. Treatments and rehabilitation in the acute and chronic state of traumatic brain injury. Journal of Internal Medicine. 2019; 285: 608–623.
Yang ZY, Tang T, Li PF, Li XX, Wu Y, Feng DD, et al. Systematic analysis of tRNA-derived small RNAs reveals therapeutic targets of Xuefu Zhuyu decoction in the cortexes of experimental traumatic brain injury. Phytomedicine: International Journal of Phytotherapy and Phytopharmacology. 2022; 102: 154168.
Tang YL, Fang LJ, Zhong LY, Jiang J, Dong XY, Feng Z. Hub genes and key pathways of traumatic brain injury: bioinformatics analysis and in vivo validation. Neural Regeneration Research. 2020; 15: 2262–2269.
Hui J, Liu R, Zhang H, He S, Wei A. Screening and identification of critical biomarkers in erectile dysfunction: evidence from bioinformatic analysis. PeerJ. 2020; 8: e8653.
Fang S, Zhong L, Wang AQ, Zhang H, Yin ZS. Identification of Regeneration and Hub Genes and Pathways at Different Time Points after Spinal Cord Injury. Molecular Neurobiology. 2021; 58: 2643–2662.
Toro-Domínguez D, Martorell-Marugán J, López-Domínguez R, García-Moreno A, González-Rumayor V, Alarcón-Riquelme ME, et al. ImaGEO: integrative gene expression meta-analysis from GEO database. Bioinformatics (Oxford, England). 2019; 35: 880–882.
Catta-Preta R, Zdilar I, Jenner B, Doisy ET, Tercovich K, Nord AS, et al. Transcriptional Pathology Evolves over Time in Rat Hippocampus after Lateral Fluid Percussion Traumatic Brain Injury. Neurotrauma Reports. 2021; 2: 512–525.
Boone DR, Weisz HA, Willey HE, Torres KEO, Falduto MT, Sinha M, et al. Traumatic brain injury induces long-lasting changes in immune and regenerative signaling. PLoS ONE. 2019; 14: e0214741.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Cambridge (Mass.)). 2021; 2: 100141.
Hu H, Miao YR, Jia LH, Yu QY, Zhang Q, Guo AY. AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors. Nucleic Acids Research. 2019; 47: D33–D38.
Li KP, Yuan M, Wu YL, Pineda M, Zhang CM, Chen YF, et al. A High-Fat High-Fructose Diet Dysregulates the Homeostatic Crosstalk Between Gut Microbiome, Metabolome, and Immunity in an Experimental Model of Obesity. Molecular Nutrition & Food Research. 2022; 66: e2100950.
Breuer K, Foroushani AK, Laird MR, Chen C, Sribnaia A, Lo R, et al. InnateDB: systems biology of innate immunity and beyond–recent updates and continuing curation. Nucleic Acids Research. 2013; 41: D1228–33.
Wu J, Li H, He J, Tian X, Luo S, Li J, et al. Downregulation of microRNA-9-5p promotes synaptic remodeling in the chronic phase after traumatic brain injury. Cell Death & Disease. 2021; 12: 9.
Sanchez E, Blais H, Duclos C, Arbour C, Van Der Maren S, El-Khatib H, et al. Sleep from acute to chronic traumatic brain injury and cognitive outcomes. Sleep. 2022; 45: zsac123.
Mahncke HW, DeGutis J, Levin H, Newsome MR, Bell MD, Grills C, et al. A randomized clinical trial of plasticity-based cognitive training in mild traumatic brain injury. Brain: a Journal of Neurology. 2021; 144: 1994–2008.
Chambers DC, Carew AM, Lukowski SW, Powell JE. Transcriptomics and single-cell RNA-sequencing. Respirology (Carlton, Vic.). 2019; 24: 29–36.
Kim EK, Choi EJ. Compromised MAPK signaling in human diseases: an update. Archives of Toxicology. 2015; 89: 867–882.
Sharma HS, Sahib S, Tian ZR, Muresanu DF, Nozari A, Castellani RJ, et al. Protein kinase inhibitors in traumatic brain injury and repair: New roles of nanomedicine. Progress in Brain Research. 2020; 258: 233–283.
Li G, Duan L, Yang F, Yang L, Deng Y, Yu Y, et al. Curcumin suppress inflammatory response in traumatic brain injury via p38/MAPK signaling pathway. Phytotherapy Research: PTR. 2022; 36: 1326–1337.
Huang X, Liu G, Guo J, Su Z. The PI3K/AKT pathway in obesity and type 2 diabetes. International Journal of Biological Sciences. 2018; 14: 1483–1496.
Xu L, Xing Q, Huang T, Zhou J, Liu T, Cui Y, et al. HDAC1 Silence Promotes Neuroprotective Effects of Human Umbilical Cord-Derived Mesenchymal Stem Cells in a Mouse Model of Traumatic Brain Injury via PI3K/AKT Pathway. Frontiers in Cellular Neuroscience. 2019; 12: 498.
Zhou D, Liu J, Hang Y, Li T, Li P, Guo S, et al. TMT-based proteomics analysis reveals the protective effects of Xuefu Zhuyu decoction in a rat model of traumatic brain injury. Journal of Ethnopharmacology. 2020; 258: 112826.
Legerstee K, Houtsmuller AB. A Layered View on Focal Adhesions. Biology. 2021; 10: 1189.
Hlavac N, Guilhaume-Corrêa F, VandeVord PJ. Mechano-stimulation initiated by extracellular adhesion and cationic conductance pathways influence astrocyte activation. Neuroscience Letters. 2020; 739: 135405.
Gorgoulis V, Adams PD, Alimonti A, Bennett DC, Bischof O, Bishop C, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019; 179: 813–827.
Varela-Eirín M, Demaria M. Cellular senescence. Current Biology: CB. 2022; 32: R448–R452.
Schwab N, Leung E, Hazrati LN. Cellular Senescence in Traumatic Brain Injury: Evidence and Perspectives. Frontiers in Aging Neuroscience. 2021; 13: 742632.
Schwab N, Ju Y, Hazrati LN. Early onset senescence and cognitive impairment in a murine model of repeated mTBI. Acta Neuropathologica Communications. 2021; 9: 82.
Arun P, Rossetti F, Wilder DM, Sajja S, Van Albert SA, Wang Y, et al. Blast Exposure Leads to Accelerated Cellular Senescence in the Rat Brain. Frontiers in Neurology. 2020; 11: 438.
Du L, Peng GH. Complement C3 deficiency alleviates alkylation-induced retinal degeneration in mice. Eye and Vision (London, England). 2022; 9: 22.
Raso MC, Djoric N, Walser F, Hess S, Schmid FM, Burger S, et al. Interferon-stimulated gene 15 accelerates replication fork progression inducing chromosomal breakage. The Journal of Cell Biology. 2020; 219: e202002175.
Xie J, Gurler Main H, Sacks JD, Muralidhar GG, Barbolina MV. Regulation of DNA damage repair and lipid uptake by CX_3CR1 in epithelial ovarian carcinoma. Oncogenesis. 2018; 7: 37.
Alam A, Thelin EP, Tajsic T, Khan DZ, Khellaf A, Patani R, et al. Cellular infiltration in traumatic brain injury. Journal of Neuroinflammation. 2020; 17: 328.
Tobin RP, Mukherjee S, Kain JM, Rogers SK, Henderson SK, Motal HL, et al. Traumatic brain injury causes selective, CD74-dependent peripheral lymphocyte activation that exacerbates neurodegeneration. Acta Neuropathologica Communications. 2014; 2: 143.
Chen H, Zhang Y, Li X, Zhang W, He H, Du B, et al. Transcriptome Changes and Potential Immunotoxicity Analysis in RAW264.7 Macrophages Caused by Bisphenol F. Frontiers in Pharmacology. 2022; 13: 846562.
Herzog C, Pons Garcia L, Keatinge M, Greenald D, Moritz C, Peri F, et al. Rapid clearance of cellular debris by microglia limits secondary neuronal cell death after brain injury in vivo. Development (Cambridge, England). 2019; 146: dev174698.
Zhou H, Hu L, Li J, Ruan W, Cao Y, Zhuang J, et al. AXL kinase-mediated astrocytic phagocytosis modulates outcomes of traumatic brain injury. Journal of Neuroinflammation. 2021; 18: 154.
Wang K, Li J, Zhang Y, Huang Y, Chen D, Shi Z, et al. Central nervous system diseases related to pathological microglial phagocytosis. CNS Neuroscience & Therapeutics. 2021; 27: 528–539.
Yanuck SF. Microglial Phagocytosis of Neurons: Diminishing Neuronal Loss in Traumatic, Infectious, Inflammatory, and Autoimmune CNS Disorders. Frontiers in Psychiatry. 2019; 10: 712.
Amara U, Flierl MA, Rittirsch D, Klos A, Chen H, Acker B, et al. Molecular intercommunication between the complement and coagulation systems. Journal of Immunology (Baltimore, Md.: 1950). 2010; 185: 5628–5636.
Luo W, Yang Z, Zhang W, Zhou D, Guo X, Wang S, et al. Quantitative Proteomics Reveals the Dynamic Pathophysiology Across Different Stages in a Rat Model of Severe Traumatic Brain Injury. Frontiers in Molecular Neuroscience. 2022; 14: 785938.
Elmore MRP, Najafi AR, Koike MA, Dagher NN, Spangenberg EE, Rice RA, et al. Colony-stimulating factor 1 receptor signaling is necessary for microglia viability, unmasking a microglia progenitor cell in the adult brain. Neuron. 2014; 82: 380–397.
Wang Y, Wernersbach I, Strehle J, Li S, Appel D, Klein M, et al. Early posttraumatic CSF1R inhibition via PLX3397 leads to time- and sex-dependent effects on inflammation and neuronal maintenance after traumatic brain injury in mice. Brain, Behavior, and Immunity. 2022; 106: 49–66.
Henry RJ, Ritzel RM, Barrett JP, Doran SJ, Jiao Y, Leach JB, et al. Microglial Depletion with CSF1R Inhibitor During Chronic Phase of Experimental Traumatic Brain Injury Reduces Neurodegeneration and Neurological Deficits. The Journal of Neuroscience: the Official Journal of the Society for Neuroscience. 2020; 40: 2960–2974.
Hickman S, Izzy S, Sen P, Morsett L, El Khoury J. Microglia in neurodegeneration. Nature Neuroscience. 2018; 21: 1359–1369.
Puntambekar SS, Saber M, Lamb BT, Kokiko-Cochran ON. Cellular players that shape evolving pathology and neurodegeneration following traumatic brain injury. Brain, Behavior, and Immunity. 2018; 71: 9–17.
Fu Y, Wang X, Zhang L, Ren Y, Hao L. Allograft inflammatory factor-1 enhances inflammation and oxidative stress via the NF-κB pathway in diabetic kidney disease. Biochemical and Biophysical Research Communications. 2022; 614: 63–69.
Sikora M, Kopeć B, Piotrowska K, Pawlik A. Role of allograft inflammatory factor-1 in pathogenesis of diseases. Immunology Letters. 2020; 218: 1–4.
Coelho FS, Rodpai R, Miller A, Karinshak SE, Mann VH, Dos Santos Carvalho O, et al. Diminished adherence of Biomphalaria glabrata embryonic cell line to sporocysts of Schistosoma mansoni following programmed knockout of the allograft inflammatory factor. Parasites & Vectors. 2020; 13: 511.
Hiskens MI, Schneiders AG, Vella RK, Fenning AS. Repetitive mild traumatic brain injury affects inflammation and excitotoxic mRNA expression at acute and chronic time-points. PLoS ONE. 2021; 16: e0251315.
Castranio EL, Mounier A, Wolfe CM, Nam KN, Fitz NF, Letronne F, et al. Gene co-expression networks identify Trem2 and Tyrobp as major hubs in human APOE expressing mice following traumatic brain injury. Neurobiology of Disease. 2017; 105: 1–14.
Li RY, Qin Q, Yang HC, Wang YY, Mi YX, Yin YS, et al. TREM2 in the pathogenesis of AD: a lipid metabolism regulator and potential metabolic therapeutic target. Molecular Neurodegeneration. 2022; 17: 40.
Colonna M, Wang Y. TREM2 variants: new keys to decipher Alzheimer disease pathogenesis. Nature Reviews. Neuroscience. 2016; 17: 201–207.
Saber M, Kokiko-Cochran O, Puntambekar SS, Lathia JD, Lamb BT. Triggering Receptor Expressed on Myeloid Cells 2 Deficiency Alters Acute Macrophage Distribution and Improves Recovery after Traumatic Brain Injury. Journal of Neurotrauma. 2017; 34: 423–435.
van Dalen FJ, Bakkum T, van Leeuwen T, Groenewold M, Deu E, Koster AJ, et al. Application of a Highly Selective Cathepsin S Two-step Activity-Based Probe in Multicolor Bio-Orthogonal Correlative Light-Electron Microscopy. Frontiers in Chemistry. 2021; 8: 628433.
Fuchs N, Meta M, Schuppan D, Nuhn L, Schirmeister T. Novel Opportunities for Cathepsin S Inhibitors in Cancer Immunotherapy by Nanocarrier-Mediated Delivery. Cells. 2020; 9: 2021.
Xu J, Wang H, Ding K, Lu X, Li T, Wang J, et al. Inhibition of cathepsin S produces neuroprotective effects after traumatic brain injury in mice. Mediators of Inflammation. 2013; 2013: 187873.
Jin C, Shao Y, Zhang X, Xiang J, Zhang R, Sun Z, et al. A Unique Type of Highly-Activated Microglia Evoking Brain Inflammation via Mif/Cd74 Signaling Axis in Aged Mice. Aging and Disease. 2021; 12: 2125–2139.
Wirtz TH, Saal A, Bergmann I, Fischer P, Heinrichs D, Brandt EF, et al. Macrophage migration inhibitory factor exerts pro-proliferative and anti-apoptotic effects via CD74 in murine hepatocellular carcinoma. British Journal of Pharmacology. 2021; 178: 4452–4467.
Li S, Zhu Y, Wei C, Li C, Chen W, Jiang S, et al. Identification of Molecular Correlations Between DHRS4 and Progressive Neurodegeneration in Amyotrophic Lateral Sclerosis By Gene Co-Expression Network Analysis. Frontiers in Immunology. 2022; 13: 874978.
Xu J, Zhou H, Xiang G. Identification of Key Biomarkers and Pathways for Maintaining Cognitively Normal Brain Aging Based on Integrated Bioinformatics Analysis. Frontiers in Aging Neuroscience. 2022; 14: 833402.
Yan XJ, Li YB, Liu W, Wu HY, Yu GF. Elevated Serum Complement C1q Levels After Traumatic Brain Injury and Its Association with Poor Prognosis. Neuropsychiatric Disease and Treatment. 2022; 18: 47–55.
Ciechanowska A, Ciapała K, Pawlik K, Oggioni M, Mercurio D, De Simoni MG, et al. Initiators of Classical and Lectin Complement Pathways Are Differently Engaged after Traumatic Brain Injury-Time-Dependent Changes in the Cortex, Striatum, Thalamus and Hippocampus in a Mouse Model. International Journal of Molecular Sciences. 2020; 22: 45.
Holden SS, Grandi FC, Aboubakr O, Higashikubo B, Cho FS, Chang AH, et al. Complement factor C1q mediates sleep spindle loss and epileptic spikes after mild brain injury. Science (New York, N.Y.). 2021; 373: eabj2685.
Mukherjee S, Klaus C, Pricop-Jeckstadt M, Miller JA, Struebing FL. A Microglial Signature Directing Human Aging and Neurodegeneration-Related Gene Networks. Frontiers in Neuroscience. 2019; 13: 2.
Toledano Furman NE, Prabhakara KS, Bedi S, Cox CS, Jr, Olson SD. OMIP-041: Optimized multicolor immunofluorescence panel rat microglial staining protocol. Cytometry. Part a: the Journal of the International Society for Analytical Cytology. 2018; 93: 182–185.
Toledano Furman N, Gottlieb A, Prabhakara KS, Bedi S, Caplan HW, Ruppert KA, et al. High-resolution and differential analysis of rat microglial markers in traumatic brain injury: conventional flow cytometric and bioinformatics analysis. Scientific Reports. 2020; 10: 11991.
Scott MC, Prabhakara KS, Walters AJ, Olson SD, Cox CS, Jr. Determining Sex-Based Differences in Inflammatory Response in an Experimental Traumatic Brain Injury Model. Frontiers in Immunology. 2022; 13: 753570.
Song Z, Wang X, He L, Chen L, Ren Z, Song S. Suppression of lysosomal-associated protein transmembrane 5 ameliorates cardiac function and inflammatory response by inhibiting the nuclear factor-kappa B (NF-κB) pathway after myocardial infarction in mice. Experimental Animals. 2022; 71: 415–425.
Du H, Shi J, Wang M, An S, Guo X, Wang Z. Analyses of gene expression profiles in the rat dorsal horn of the spinal cord using RNA sequencing in chronic constriction injury rats. Journal of Neuroinflammation. 2018; 15: 280.
Schmeier S, Alam T, Essack M, Bajic VB. TcoF-DB v2: update of the database of human and mouse transcription co-factors and transcription factor interactions. Nucleic Acids Research. 2017; 45: D145–D150.
Hou Y, Sun B, Liu W, Yu B, Shi Q, Luo F, et al. Targeting of glioma stem-like cells with a parthenolide derivative ACT001 through inhibition of AEBP1/PI3K/AKT signaling. Theranostics. 2021; 11: 555–566.
Bogachev O, Majdalawieh A, Pan X, Zhang L, Ro HS. Adipocyte enhancer-binding protein 1 (AEBP1) (a novel macrophage proinflammatory mediator) overexpression promotes and ablation attenuates atherosclerosis in ApoE (-/-) and LDLR (-/-) mice. Molecular Medicine (Cambridge, Mass.). 2011; 17: 1056–1064.
Shijo M, Honda H, Suzuki SO, Hamasaki H, Hokama M, Abolhassani N, et al. Association of adipocyte enhancer-binding protein 1 with Alzheimer’s disease pathology in human hippocampi. Brain Pathology (Zurich, Switzerland). 2018; 28: 58–71.
Jansen K, Wirz OF, van de Veen W, Tan G, Mirer D, Sokolowska M, et al. Loss of regulatory capacity in Treg cells following rhinovirus infection. The Journal of Allergy and Clinical Immunology. 2021; 148: 1016–1029.e16.
Johnson NH, Hadad R, Taylor RR, Rodríguez Pilar J, Salazar O, Llompart-Pou JA, et al. Inflammatory Biomarkers of Traumatic Brain Injury. Pharmaceuticals (Basel, Switzerland). 2022; 15: 660.
Bortolotti P, Faure E, Kipnis E. Inflammasomes in Tissue Damages and Immune Disorders After Trauma. Frontiers in Immunology. 2018; 9: 1900.
Wang J, Xu C, Zhang K, Shi J, Liu Z, Wang X, et al. Inhibition of ASC enhances the protective role of salvianolic acid A in traumatic brain injury via inhibition of inflammation and recovery of mitochondrial function. Folia Neuropathologica. 2021; 59: 50–66.
Kuwar R, Rolfe A, Di L, Xu H, He L, Jiang Y, et al. A novel small molecular NLRP3 inflammasome inhibitor alleviates neuroinflammatory response following traumatic brain injury. Journal of Neuroinflammation. 2019; 16: 81.
Lucke-Wold BP, Turner RC, Logsdon AF, Nguyen L, Bailes JE, Lee JM, et al. Endoplasmic reticulum stress implicated in chronic traumatic encephalopathy. Journal of Neurosurgery. 2016; 124: 687–702.
Turner RC, Lucke-Wold BP, Logsdon AF, Robson MJ, Lee JM, Bailes JE, et al. Modeling Chronic Traumatic Encephalopathy: The Way Forward for Future Discovery. Frontiers in Neurology. 2015; 6: 223.
Jiang T, Tan L, Zhu XC, Zhou JS, Cao L, Tan MS, et al. Silencing of TREM2 exacerbates tau pathology, neurodegenerative changes, and spatial learning deficits in P301S tau transgenic mice. Neurobiology of Aging. 2015; 36: 3176–3186.
Zhu B, Liu Y, Hwang S, Archuleta K, Huang H, Campos A, et al. Trem2 deletion enhances tau dispersion and pathology through microglia exosomes. Molecular Neurodegeneration. 2022; 17: 58.
Audrain M, Haure-Mirande JV, Wang M, Kim SH, Fanutza T, Chakrabarty P, et al. Integrative approach to sporadic Alzheimer’s disease: deficiency of TYROBP in a tauopathy mouse model reduces C1q and normalizes clinical phenotype while increasing spread and state of phosphorylation of tau. Molecular Psychiatry. 2019; 24: 1383–1397.
Chidambaram H, Das R, Chinnathambi S. Interaction of Tau with the chemokine receptor, CX3CR1 and its effect on microglial activation, migration and proliferation. Cell & Bioscience. 2020; 10: 109.
Bolós M, Llorens-Martín M, Perea JR, Jurado-Arjona J, Rábano A, Hernández F, et al. Absence of CX3CR1 impairs the internalization of Tau by microglia. Molecular Neurodegeneration. 2017; 12: 59.
Lucke-Wold BP, Logsdon AF, Turner RC, Huber JD, Rosen CL. Endoplasmic Reticulum Stress Modulation as a Target for Ameliorating Effects of Blast Induced Traumatic Brain Injury. Journal of Neurotrauma. 2017; 34: S62–S70.
Noda Y, Tsuruma K, Takata M, Ishisaka M, Tanaka H, Nakano Y, et al. GPNMB Induces BiP Expression by Enhancing Splicing of BiP Pre-mRNA during the Endoplasmic Reticulum Stress Response. Scientific Reports. 2017; 7: 12160.
Kam TI, Song S, Gwon Y, Park H, Yan JJ, Im I, et al. FcγRIIb mediates amyloid-β neurotoxicity and memory impairment in Alzheimer’s disease. The Journal of Clinical Investigation. 2013; 123: 2791–2802.
Barichello T, Simões LR, Collodel A, Giridharan VV, Dal-Pizzol F, Macedo D, et al. The translocator protein (18kDa) and its role in neuropsychiatric disorders. Neuroscience and Biobehavioral Reviews. 2017; 83: 183–199.
Filippi L, Schillaci O, Palumbo B. Neuroimaging with PET/CT in chronic traumatic encephalopathy: what nuclear medicine can do to move the field forward. Expert Review of Molecular Diagnostics. 2022; 22: 149–156.
Li Y, Chen L, Li L, Sottas C, Petrillo SK, Lazaris A, et al. Cholesterol-binding translocator protein TSPO regulates steatosis and bile acid synthesis in nonalcoholic fatty liver disease. IScience. 2021; 24: 102457.
Stormezand GN, Doorduin J, Rakers SE, Spikman JM, van der Naalt J, García DV, et al. Imaging of neuroinflammation due to repetitive head injury in currently active kickboxers. European Journal of Nuclear Medicine and Molecular Imaging. 2022; 49: 3162–3172.
Endepols H, Anglada-Huguet M, Mandelkow E, Schmidt Y, Krapf P, Zlatopolskiy BD, et al. Assessment of the In Vivo Relationship Between Cerebral Hypometabolism, Tau Deposition, TSPO Expression, and Synaptic Density in a Tauopathy Mouse Model: a Multi-tracer PET Study. Molecular Neurobiology. 2022; 59: 3402–3413.
Gouna G, Klose C, Bosch-Queralt M, Liu L, Gokce O, Schifferer M, et al. TREM2-dependent lipid droplet biogenesis in phagocytes is required for remyelination. The Journal of Experimental Medicine. 2021; 218: e20210227.
Yang ZY, Wu Y, Li X, Tang T, Wang Y, Huang ZB, et al. Bioinformatics Analysis of miRNAs and mRNAs Network-Xuefu Zhuyu Decoction Exerts Neuroprotection of Traumatic Brain Injury Mice in the Subacute Phase. Frontiers in Pharmacology. 2022; 13: 772680.
Lai JQ, Shi YC, Lin S, Chen XR. Metabolic disorders on cognitive dysfunction after traumatic brain injury. Trends in Endocrinology and Metabolism: TEM. 2022; 33: 451–462.
Kolaczkowska E, Kubes P. Neutrophil recruitment and function in health and inflammation. Nature Reviews. Immunology. 2013; 13: 159–175.
Rodríguez-Rajo P, García-Rudolph A, Sánchez-Carrión R, Aparicio-López C, Enseñat-Cantallops A, García-Molina A. Computerized social cognitive training in the subacute phase after traumatic brain injury: A quasi-randomized controlled trial. Applied Neuropsychology. Adult. 2022; 1–14. (Online ahead of print)
Wang ML, Wei XE, Yu MM, Li WB. Cognitive impairment in mild traumatic brain injury: a diffusion kurtosis imaging and volumetric study. Acta Radiologica (Stockholm, Sweden: 1987). 2022; 63: 504–512.
Karlsen RH, Saksvik SB, Stenberg J, Lundervold AJ, Olsen A, Rautio I, et al. Examining the Subacute Effects of Mild Traumatic Brain Injury Using a Traditional and Computerized Neuropsychological Test Battery. Journal of Neurotrauma. 2021; 38: 74–85.
Fronczak KM, Li Y, Henchir J, Dixon CE, Carlson SW. Reductions in Synaptic Vesicle Glycoprotein 2 Isoforms in the Cortex and Hippocampus in a Rat Model of Traumatic Brain Injury. Molecular Neurobiology. 2021; 58: 6006–6019.
Caplan HW, Cox CS, Bedi SS. Do microglia play a role in sex differences in TBI? Journal of Neuroscience Research. 2017; 95: 509–517.
Mollayeva T, Mollayeva S, Pacheco N, Colantonio A. Systematic Review of Sex and Gender Effects in Traumatic Brain Injury: Equity in Clinical and Functional Outcomes. Frontiers in Neurology. 2021; 12: 678971.
Blaya MO, Raval AP, Bramlett HM. Traumatic brain injury in women across lifespan. Neurobiology of Disease. 2022; 164: 105613.
Villapol S, Loane DJ, Burns MP. Sexual dimorphism in the inflammatory response to traumatic brain injury. Glia. 2017; 65: 1423–1438.
Späni CB, Braun DJ, Van Eldik LJ. Sex-related responses after traumatic brain injury: Considerations for preclinical modeling. Frontiers in Neuroendocrinology. 2018; 50: 52–66.

Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Back to top