Skip to main content
Download PDF

Identification of potential biomarkers and therapeutic targets related to post-traumatic stress disorder due to traumatic brain injury

European Journal of Medical Research volume 29, Article number: 44 (2024) Cite this article

Abstract

Background

Post-traumatic stress disorder (PTSD), a disease state that has an unclear pathogenesis, imposes a substantial burden on individuals and society. Traumatic brain injury (TBI) is one of the most significant triggers of PTSD. Identifying biomarkers associated with TBI-related PTSD will help researchers to uncover the underlying mechanism that drives disease development. Furthermore, it remains to be confirmed whether different types of traumas share a common mechanism of action.

Methods

For this study, we screened the eligible data sets from the Gene Expression Omnibus (GEO) database, obtained differentially expressed genes (DEGs) through analysis, conducted functional enrichment analysis on the DEGs in order to understand their molecular mechanisms, constructed a PPI network, used various algorithms to obtain hub genes, and finally evaluated, validated, and analyzed the diagnostic performance of the hub genes.

Results

A total of 430 upregulated and 992 down-regulated differentially expressed genes were extracted from the TBI data set. A total of 1919 upregulated and 851 down-regulated differentially expressed genes were extracted from the PTSD data set. Functional enrichment analysis revealed that the differentially expressed genes had biological functions linked to molecular regulation, cell signaling transduction, cell metabolic regulation, and immune response. After constructing a PPI network and introducing algorithm analysis, the upregulated hub genes were identified as VNN1, SERPINB2, and ETFDH, and the down-regulated hub genes were identified as FLT3LG, DYRK1A, DCN, and FKBP8. In addition, by comparing the data with patients with other types of trauma, it was revealed that PTSD showed different molecular processes that are under the influence of different trauma characteristics and responses.

Conclusions

By exploring the role of different types of traumas during the pathogenesis of PTSD, its possible molecular mechanisms have been revealed, providing vital information for understanding the complex pathways associated with TBI-related PTSD. The data in this study has important implications for the design and development of new diagnostic and therapeutic methods needed to treat and manage PTSD.

Introduction

According to statistics, more than 50 million people worldwide suffer from traumatic brain injury (TBI) each year, and over the past 30 years, both the incidence and prevalence have increased [1, 2]. This type of injury imposes a significant burden on global healthcare institutions. In addition, recent studies have found that TBI will evolve over time [3]. Over the past decade, there has been a fundamental shift in understanding TBI, as experts now recognize it as a chronic condition affecting multiple bodily systems [4]. Patients often experience persistent cognitive, emotional, and functional impairments in their daily lives, greatly impacting their quality of life [5]. Among these impairments, post-traumatic stress disorder (PTSD) is one of the most typical presentations [6, 7]. According to the fifth edition of the diagnostic and statistical manual of mental disorders [8], PTSD encompasses four symptomatic clusters. These include the presence of intrusive memories related to the stressful event, persistent avoidance of stimuli associated with the event, negative alterations in mood or cognition, and hyperarousal or excessive vigilance. These symptoms must present within one month following the experience of a distressing or traumatic event. However, symptoms may have a delayed onset in some cases, occurring years later [9]. It is estimated that PTSD affects approximately 5–10% of the global population, with a gender ratio of approximately 1:2, favoring females over males [10]. In countries and regions with an internal conflict backdrop, the prevalence of PTSD is estimated to be even higher, potentially affecting more than 25% of the population [11]. PTSD is considered a heterogeneous disorder, with multiple pathological pathways believed to be involved in its onset and progression. In recent years, various studies have been exploring the possible biological mechanisms that underlie PTSD. These include the hypothalamic–pituitary–adrenal (HPA) axis, neurochemical factors, autonomic nervous system, inflammation, and immune dysregulation, which have all been proposed as potential neuroimmune biomarkers for PTSD [12,13,14,15,16]. However, recent literature review analyses have shown conflicting data on the role of some biomarkers in individuals with PTSD [17]. One of the reasons for such discrepancies stems from current molecular mechanism research concerning PTSD, which relies primarily on animal experiments that necessitate further investigation using human samples [18]. Conversely, the development of the disorder also depends on the characteristics of the trauma and individual risk factors [19]. To screen and explore potential biomarkers as well as the therapeutic targets used to treat PTSD that TBI causes, we conducted further analysis of the molecular mechanisms of TBI-associated PTSD. We utilized publicly available human specimen data sets related to TBI and PTSD and identified differentially expressed genes (DEGs) using bioinformatics methods. Subsequently, we performed GO and KEGG analyses and constructed protein–protein interaction (PPI) networks in order to elucidate the functional classification and the metabolic pathways of the DEGs. We also identified potential hub genes. To investigate whether different traumatic features can influence the molecular mechanisms of PTSD, we introduced a third human specimen data set related to fracture injury for validation. This allowed us to explore whether different traumatic features have shared mechanisms in their impact on PTSD, as well as to determine whether TBI-related PTSD has distinct characteristics that warrant further independent research.

Methods

Data set information

The keywords “trauma” and “PTSD” were used as screening terms to search the GEO database. The search criteria were further limited to “Homo sapiens” and “Expression profiling by array.” The search results were analyzed individually, and the original data set literature that met the research criteria was carefully read and reviewed. Finally, the GSE223245 and GSE81761 data sets were identified as subjects for further research. GSE223245 includes 12 TBI patients and 4 healthy controls, while GSE81761 includes 39 PTSD patients and 27 healthy controls. Both data sets collected venous blood for gene detection. The raw data of both data sets were downloaded using R software (version 4.2.1, http://r-project.org/). Then, the “affy” package [20] within the R software was used to read the raw data of the two data sets used for background correction and data normalization.

Differential gene expression analysis

The “limma” package [21] was used to analyze the two data sets separately. The disease group and control group within each data set were labeled. Differentially expressed genes between the two groups in each data set were obtained by performing comparative analysis. Since the molecular mechanisms of PTSD are not yet well understood, in order to screen for more valuable genes, the criteria were set as a |log2 (fold-change)|> 0.2 and a p value < 0.05 to include genes as differentially expressed genes. Those that met these two criteria were considered to have significant differences. Among the selected DEGs, genes with logFC > 0.2 were considered upregulated, while genes with logFC < − 0.2 were considered downregulated.

GO and KEGG pathway enrichment analysis

The GO database is a database established by the Gene Ontology Consortium and is primarily divided into three semantic categories: cellular component (CC), molecular function (MF), and biological process (BP). GO enrichment analysis can be used to determine the functions of genes and provide preliminary annotations by analyzing the sequence information of the genes [22]. The KEGG database is a database for systematic analysis of gene products and their metabolic pathways in cells and is commonly used for metabolic pathway analysis [23]. The “clusterProfiler” package in R software [24] is used for enrichment analysis of DEGs, followed by the use of the “org.Hs.eg.db” package for ID conversion [25]. Finally, the “ggplot2” package is used to visualize the enrichment analysis results.

PPI network construction

Analysis of the common differentially expressed genes obtained from the two data sets is uploaded to the STRING database [26] (http://string-db.org/). The STRING database is an online database for searching known protein–protein interaction relationships, which helps in the exploration of core regulatory genes. The output results are imported into the Cytoscape software [27]. Algorithms are used to explore the potential correlations among these common differentially expressed genes.

Hub gene selection

Hub genes play essential roles in biological processes and often act as key regulators of other genes in different molecular pathways. The Cyto-Hubba plugin [28] in the Cytoscape software can identify hub genes within the PPI network. To avoid selection bias caused by the different algorithms, the hub genes identified by the different algorithms intersect in order to define the final hub genes.

Preliminary analysis of hub genes

The PANTHER Classification System [29] (http://pantherdb.org/) can be used to classify proteins (and their genes). It is a comprehensive and annotated gene family system database. The Human Protein Atlas [30] (https://www.proteinatlas.org/) is a public database for querying gene expression profiles in different human organs. It can provide a list of the basic RNA and protein expression levels of specific genes. Selecting the “brain” section of the database, we compared the transcriptome data of 13 major brain regions to classify human brain regions and display the standardized RNA expression levels of the 13 brain regions. We can conduct a preliminary analysis of the hub genes using these two commonly utilized public databases.

Diagnostic performance of hub genes

ROC curve analysis was conducted on each hub gene data set to assess the diagnostic performance of hub genes. The area under the ROC curve (AUC) of the hub genes was calculated using the “pROC” package in R software. This analysis evaluated the accuracy of diagnosis and the ability to recognize diseases.

Validation of hub genes

PTSD is a heterogeneous disease with multiple pathophysiological pathways potentially involved in its onset and progression. It requires expanding the data set of the different types of trauma to assess whether different events trigger distinct pathophysiological responses [17]. A recent prospective study also indicated a significantly higher proportion of PTSD occurrence after a TBI compared to other injuries, such as fractures [31]. Therefore, in order to differentiate the molecular mechanisms that underlie the impact of different trauma characteristics on PTSD, we introduced a third human specimen data set focusing on fracture injuries for validation. This validation aims to explore whether different trauma characteristics have common mechanisms affecting PTSD and to confirm if the specificity of PTSD related to TBI warrants further independent research.

Results

Differential expression gene screening

After data processing of the two data sets, in the GSE223245 TBI data set, there was a total of 430 genes identified as upregulated DEGs that met the criteria of the log2(fold-change) > 0.2 and a p value < 0.05. In addition, 992 genes were identified as down-regulated DEGs that met the log2 (fold-change) criteria < 0.2 and a p value < 0.05. In the GSE81761 PTSD data set, there were a total of 1919 genes identified as upregulated DEGs that met the criteria of the log2 (fold-change) > 0.2 and a p value < 0.05 and 851 genes identified as down-regulated DEGs that met the criteria of the log2 (fold-change) < 0.2 and a p value < 0.05. The volcano plot was used to visualize the screening results of the differential expression genes (Fig. 1).

Fig. 1
figure 1

DEGs in the two data sets. A Volcano plot of GSE223245; red represents upregulated genes, blue represents downregulated genes, and gray represents nonsignificantly expressed genes. B Volcano plot of GSE81761. The criteria for statistically significant differences in DEGs were adjusted |log2(FC)|> 0.2 and p.adj < 0.05 in expression

Full size image

GO and KEGG pathway enrichment analysis

The “clusterProfiler” package was used to perform GO and KEGG pathway enrichment analysis on the TBI and PTSD data sets independently. The results showed that in the GSE223245 TBI data set, the DEGs were primarily enriched in biological processes (BP) such as viral life cycle regulation (GO:1903900), protein autophosphorylation (GO:0046777), the immune response-regulating signaling pathway (GO:0002764), viral process regulation (GO:0050792), and peptidyl-tyrosine phosphorylation (GO:0018108). In terms of molecular functions (MF), the DEGs were primarily associated with phosphatidylserine binding (GO:0001786), modified amino acid binding (GO:0072341), nuclease activity (GO:0004518), exodeoxyribonuclease activity (GO:0004529), and exodeoxyribonuclease activity that produces 5’-phosphomonoesters (GO:0016895). The enriched KEGG pathways were related to Th1 and Th2 cell differentiation (hsa04658) as well as pantothenate and CoA biosynthesis (hsa00770). In the GSE81761 PTSD data set, the DEGs were predominantly enriched in BP such as proteasomal protein catabolic processes (GO:0010498), proteasome-mediated ubiquitin-dependent protein catabolic processes (GO:0043161), mRNA processing (GO:0006397), positive regulation of cellular catabolic processes (GO:0031331), and the regulation of cellular protein catabolic processes (GO:1903362). In terms of cellular components (CC), the focus was on methyltransferase complexes (GO:0034708), ubiquitin ligase complexes (GO:0000151), nuclear speck (GO:0016607), spindle pole (GO:0000922), and nuclear pore nuclear baskets (GO:0044615). MF was associated with transcription coregulator activity (GO:0003712) and transcription coactivator activity (GO:0003713). The enrichment analysis results indicate that the BPs, CCs, and MFs are associated with different biological pathways, including molecular regulation, cell signaling transduction, cell metabolic regulation, and immune response, suggesting further refinement and classification to facilitate more in-depth research. Visualization was performed using circular plots (Fig. 2) to present the enrichment results clearly. The enrichment results of all the data sets were listed in a ternary table format (Tables 1, 2).

Fig. 2
figure 2

GO terms and KEGG pathway enrichment. A Details of GO terms and KEGG pathway enrichment in GSE223245. The inner ring is a bar plot where the height of the bar indicates the significance of the term (− log10 p value), and color corresponds to the z score. The outer ring displays scatterplots of the expression levels (logFC) for the genes in each term. B Details of GO terms and KEGG pathway enrichment in GSE81761

Full size image
Table 1 Details of GO terms and KEGG pathway enrichment in GSE223245
Full size table
Table 2 Details of GO terms and KEGG pathway enrichment in GSE81761
Full size table

PPI network construction

The overlapping upregulated and downregulated DEGs from the two data sets were analyzed separately using the “VennDiagram” package (Fig. 3A). This resulted in 17 commonly upregulated DEGs and 36 commonly downregulated DEGs. These genes were uploaded to the STRING database to construct the PPI network. The generated data was then exported and imported into the Cytoscape software to visualize and analyze the PPI network, resulting in a graphical representation of PPI (Fig. 3B, C).

Fig. 3
figure 3

Venn diagram and PPI network. A The Venn diagram shows the co-DEGs of the two data sets. B A PPI network was constructed using coupregulated genes. C A PPI network was constructed using codownregulated genes

Full size image

Hub gene selection

The Cyto-Hubba plugin in the Cytoscape software is a widely used network visualization tool that allows for the discovery of key target information for various types of biological data, such as gene regulation and signal transduction. The plugin assigns a value to each gene using topological network algorithms and ranks them to identify the hub genes. In order to avoid bias caused by different topological analysis methods, four methods were used for analysis: MCC, DMNC, MNC, and EPC. The results were then merged and analyzed using the “VennDiagram” package, resulting in the identification of the hub genes (Fig. 4). The upregulated hub genes identified were VNN1, SERPINB2, and ETFDH, while the downregulated hub genes identified were FLT3LG, DYRK1A, DCN, and FKBP8.

Fig. 4
figure 4

The hub genes were merged by the MCC, DMNC, MNC, and EPC algorithms. A VNN1, SERPINB2 and ETFDH were upregulated hub genes. B FLT3LG, DYRK1A, DCN and FKBP8 were downregulated hub genes

Full size image

Preliminary analysis of hub genes

First, the hub genes were subjected to a one-stop annotation analysis using the PANTHER Classification System (http://pantherdb.org/) to obtain an initial understanding of their functions (Table 3). Considering that the hub genes may exhibit differential expression in different brain regions during cranial brain injury processing, the Human Protein Atlas database (https://www.proteinatlas.org/) was used to compile the basic RNA and protein expression levels of the hub genes. These expression levels were then displayed across different brain regions based on standardized RNA expression levels. The color coding was based on the brain regions, and the bar chart depicts the highest expression levels within the included subregions (Fig. 5).

Table 3 The PANTHER classification system of hub genes
Full size table
Fig. 5
figure 5

The Human Protein Atlas was used to analyze the hub genes. Regional classification in the human brain by comparing transcriptome data aggregated to 13 major regions of the brain. Normalized RNA expression levels shown for the 13 brain regions. Color coding is based on brain region, and the bar shows the highest expression among the subregions included. A VNN1 B SERPINB2 C ETFDH D FLT3LG E DYRK1A F DCN G FKBP8

Full size image

Diagnostic performance of hub genes

To validate the diagnostic performance of the hub genes, we conducted a ROC curve analysis based on the two GEO data sets (Fig. 6). The results showed that, except for DCN, all of the hub genes exhibited good diagnostic performance within the TBI data set. Considering that the gene expression differences between the PTSD patients and the normal controls may be relatively small, thus, the diagnostic performance of the hub genes within the PTSD data set may not be as expected. However, the results indicated that the area under the ROC curve for all of the hub genes in the PTSD data set was > 60%. Therefore, we believe that the selected hub genes from this study are key genes in the pathogenesis of TBI-related PTSD and possess a significantly high diagnostic value.

Fig. 6
figure 6

GEO data set was used to verify the diagnostic performance of the hub genes. A ROC curve of hub genes in the GSE223245 data set. B ROC curve of hub genes in the GSE81761 data set

Full size image

Validation of hub genes

To speculate whether different trauma characteristics affect the molecular mechanisms of PTSD, we introduced a third data set for validation due to the involvement of different biological pathways containing different trauma characteristics[19]. Using the same selection criteria as before, we selected the GSE93138 fracture data set as the validation data set. After preprocessing the raw data, screening DEGs, constructing PPI networks, and identifying hub genes, we obtained the upregulated hub genes in the fracture and PTSD data sets: RBX1, COMMD8, COMMD6, DCUN1D1, DCUN1D5, TCEB1, COPS2, and CUL3. The downregulated hub genes identified were CD79B, CD79A, CD40, PAX5, FCER2, HNRNPH1, G3BP1, IL2RG, SRRM2, and APP. By comparing the hub genes between the different groups, it was observed that for PTSD, the target genes involved in the various trauma mechanisms are not identical. Thus, it can be reasonably inferred that there are differences in the molecular mechanisms during the occurrence and development of PTSD under different trauma characteristics. This also explains why some biomarkers show contradictory results in PTSD patients [17]. These validation results suggest that researchers should pay attention to the effects of the trauma characteristics when studying the mechanisms of PTSD, thus providing an increase in best practices used in guiding clinical practice and research.

Discussion

PTSD is a chronic mental disorder that leads to a decrease in patient quality of life and an increase in economic burden. Exposure to traumatic stressors is a triggering factor in the development of PTSD [32]. The first-line intervention for patients with PTSD is psychological therapy, either trauma-focused or non-trauma-focused [18]. However, these psychological therapies share common issues, such as prolonged treatment duration and significant individual variation. Many patients are unable to eliminate PTSD symptoms even after receiving treatment, leading to poor adherence to long-term therapy [33]. Currently, there are no specific drugs that have been proven to be effective in treating PTSD [34]. The two first-line medications approved by the US FDA for the clinical treatment of PTSD are sertraline and paroxetine; both are selective serotonin reuptake inhibitors (SSRIs), which are commonly used as antidepressants [35]. However, these medications have issues such as low efficacy rates, delayed onset of action, and significant adverse effects [36]. Therefore, conducting in-depth research on the pathogenesis of PTSD, identifying new potential treatment targets, and exploring novel treatment strategies are of significant importance in improving the treatment outcomes for PTSD patients. In recent years, there has been increased attention surrounding PTSD biomarker research [37]. These biomarkers may aid in screening and supporting early detection of PTSD, as well as identifying drug targets that will lead to timely intervention and better outcomes for individuals suffering from PTSD [38,39,40]. Based on our research findings, we have identified VNN1, SERPINB2, and ETFDH as upregulated hub genes and FLT3LG, DYRK1A, DCN, and FKBP8 as downregulated hub genes associated with PTSD following TBI. According to the review of previous literature, Vanin-1 (VNN1) is a glycosylphosphatidylinositol-anchored pantetheinase that catalyzes the hydrolysis of pantetheine into cysteamine and pantothenic acid [41]. Studies have shown that VNN1 plays a role in oxidative stress, inflammation, and cell migration [42,43,44,45]. In addition, VNN1’s cysteamine pantetheinase activity contributes to local tissue adaptation, injury resilience, and tissue stress tolerance [46]. A recent study demonstrated that gene knockout of VNN1 prevented organ damage from progressing from acute to chronic stages [47]. These mechanistic studies and animal experiments provide theoretical support for the role of VNN1 in PTSD. Serpin family B member 2 (SERPINB2), also known as plasminogen activator inhibitor-2 (PAI-2), is an inhibitor of the extracellular protease plasminogen activator and is expressed in various cell types [48]. SERPINB2 participates in multiple cellular functions, including cell survival, cell differentiation, inflammation, immunity, cell adhesion, migration, and extracellular matrix remodeling in multiple disease states via its interactions with intra and extracellular proteins [49,50,51]. Despite extensive research on SERPINB2, no consensus exists on its pathophysiological functions [52]. However, a recent study suggests that the downregulation of SERPINB2 may have potential therapeutic value in relieving chronic visceral pain syndromes, such as chronic pelvic pain syndrome within the urological system [53]. This improvement in chronic pain-related conditions also suggests that SERPINB2 may hold some value in the treatment of PTSD. Electron transfer flavoprotein dehydrogenase (ETFDH), also known as ETF–ubiquinone oxidoreductase (ETF–QO), is a protein located in the inner mitochondrial membrane and plays a central role in the electron transfer system [54]. It is closely associated with Multiple Acyl-CoA Dehydrogenase Deficiency (MADD). However, recent studies have also found that ETFDH plays an important role in the development of several diseases [55, 56], suggesting the need for future research in order to explore its association with biological pathways that drive PTSD.Fms-related receptor tyrosine kinase 3 ligand (FLT3LG) is a growth factor that binds to and forms non-covalent dimers with FLT3 (CD135), acting as an agonist [57]. FLT3LG stimulates the differentiation and proliferation of dendritic cells [58] and has a potential therapeutic value in treating autoimmune diseases and tumors [57, 59]. However, there is relatively limited research on FLT3LG, and further exploration is needed regarding its expression levels, functions, and clinical value in various diseases [60]. Dual-specificity tyrosine phosphorylation regulated kinase 1A (DYRK1A) is an evolutionarily conserved protein kinase and the most extensively studied member of the dual-specificity tyrosine-regulated kinase (DYRK) family. Studies have shown that it is involved in the development of many diseases, and both low and high expression levels of the DYRK1A protein contribute significantly to disease pathology [61]. A recent study suggests that elevated plasma levels of DYRK1A might help prevent neurodegenerative diseases [62], highlighting the potential value of this target in detecting, monitoring, and managing neurodegenerative disorders. As PTSD is a chronic mental disorder, there is a theoretical basis to further explore the role of DYRK1A in neuropsychiatric disorders. Decorin (DCN) is an extracellular matrix protein belonging to the family of small leucine-rich proteoglycans. As a multifunctional molecule, DCN assembles the extracellular matrix and regulates the biological activity of cellular growth factors. Studies [63,64,65,66,67] have shown that DCN acts as a ligand for various cytokines and growth factors by directly or indirectly interacting with signaling molecules involved in cell growth, differentiation, proliferation, adhesion, and migration. The multifunctional nature of DCN makes it a potential therapeutic target for various diseases and shows promising clinical and research applications [68]. FKBP prolyl isomerase 8 (FKBP8) is a member of the FK506-binding protein family, usually found in the mitochondria and the endoplasmic reticulum. FKBP8 plays a critical role in cellular functions, including protein transport and folding [69]. It has been found that FKBP8 plays a vital role in intracellular transport, protein folding, cell apoptosis, cell growth, and differentiation. It also participates in various cellular processes associated with immune regulatory effects [69, 70]. It is worth further investigation into whether FKBP8 plays a role in the occurrence and development of PTSD. At the same time, it should be emphasized that gene expression is generally regulated at transcription and translation levels. The correlation between the expression levels of mRNA and its corresponding protein depends on many regulatory factors and metabolic processes. The transcription process from DNA to mRNA and then to protein involves a set of expression regulation mechanisms that drive mRNA formation and include regulatory activities such as post-transcriptional regulation, epigenetic modifications (such as DNA methylation), and post-translational regulation. The abundance of mRNA and protein expression may not be consistent, and the differences between the two can also suggest more biological significance and regulatory mechanisms. Therefore, a detailed exploration of the specific biological gene expression processes is needed.

One of the limitations of our study is that only two data sets were included in this article. If there is bias in these two data sets, it may affect the results of the analysis. At the same time, the results of potential biomarkers and therapeutic targets screened by bioinformatics can only be used for theoretical guidance to facilitate researchers’ design of future clinical trials. Another limitation is that we did not analyze factors such as race, gender, education level, and dietary habits that may affect outcomes. However, the results of this study can be seen as a supplemental and validating contribution to existing theories. It informs researchers of the potential biological mechanisms and pathways that need to be considered in all areas of PTSD research. This can aid in developing intervention strategies that target these biomarkers as well as the related biological pathways. Our study also has the potential to help expand the pharmacological approaches used to prevent and treat individuals who experience traumatic stress.

Conclusions

This study utilized a bioinformatic analytic approach to screen potential biomarkers and therapeutic targets (VNN1, SERPINB2, ETFDH, FLT3LG, DYRK1A, DCN, FKBP8) related to PTSD due to TBI. The validation confirmed the potential value of these targets. In addition, by comparing data with other trauma-type patients, it was found that PTSD exhibits different mechanisms under the influence of different trauma characteristics. This provides new insights into the pathophysiology of PTSD as well as valuable references for future clinical practice.

Availability of data and materials

Publicly available data sets were analyzed in this study. The data sets generated and/or analyzed during the current study are available in the GEO repository [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE98793; https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164805; https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5972].

Abbreviations

TBI:

Traumatic brain injury

PTSD:

Post-traumatic stress disorder

GEO:

Gene expression omnibus

DEGs:

Differentially expressed genes

PPI:

Protein–protein interaction

CC:

Cellular component

BP:

Biological process

MF:

Molecular function

References

  1. Dams-O’Connor K, Juengst SB, Bogner J, et al. Traumatic brain injury as a chronic disease: insights from the United States Traumatic Brain Injury Model Systems Research Program. Lancet Neurol. 2023;22(6):517–28.

    Article  PubMed  Google Scholar 

  2. 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. Lancet Neurol. 2019;18(1):56–87.

  3. Wilson L, Stewart W, Dams-O’Connor K, et al. The chronic and evolving neurological consequences of traumatic brain injury. Lancet Neurol. 2017;16(10):813–25.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Livingston G, Huntley J, Sommerlad A, et al. Dementia prevention, intervention, and care: 2020 report of the Lancet Commission. Lancet (London, England). 2020;396(10248):413–46.

    Article  PubMed  Google Scholar 

  5. Lefkovits AM, Hicks AJ, Downing M, et al. Surviving the “silent epidemic”: a qualitative exploration of the long-term journey after traumatic brain injury. Neuropsychol Rehabil. 2021;31(10):1582–606.

    Article  PubMed  Google Scholar 

  6. Ashina H, Eigenbrodt AK, Seifert T, et al. Posttraumatic headache attributed to traumatic brain injury: classification, clinical characteristics, and treatment. Lancet Neurol. 2021;20(6):460–9.

    Article  PubMed  Google Scholar 

  7. Kong LZ, Zhang RL, Hu SH, et al. Military traumatic brain injury: a challenge straddling neurology and psychiatry. Mil Med Res. 2022;9(1):2.

    PubMed  PubMed Central  Google Scholar 

  8. Weathers FW, Bovin MJ, Lee DJ, et al. The Clinician-Administered PTSD Scale for DSM-5 (CAPS-5): development and initial psychometric evaluation in military veterans. Psychol Assess. 2018;30(3):383–95.

    Article  PubMed  Google Scholar 

  9. Bisson JI, Cosgrove S, Lewis C, et al. Posttraumatic stress disorder. BMJ (Clinical Research ed). 2015;351: h6161.

    PubMed  Google Scholar 

  10. Yehuda R, Hoge CW, McFarlane AC, et al. Posttraumatic stress disorder. Nat Rev Dis Primers. 2015;1:15057.

    Article  PubMed  Google Scholar 

  11. Hoppen TH, Priebe S, Vetter I, et al. Global burden of posttraumatic stress disorder and major depression in countries affected by war between 1989 and 2019: a systematic review and meta-analysis. BMJ Glob Health. 2021;6(7):e006303.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Miller AH. Norman Cousins Lecture. Mechanisms of cytokine-induced behavioral changes: psychoneuroimmunology at the translational interface. Brain Behav Immunity. 2009;23(2):149–58.

    Article  CAS  Google Scholar 

  13. Herman JP, McKlveen JM, Ghosal S, et al. Regulation of the hypothalamic-pituitary-adrenocortical stress response. Compr Physiol. 2016;6(2):603–21.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Geuze E, van Berckel BN, Lammertsma AA, et al. Reduced GABAA benzodiazepine receptor binding in veterans with posttraumatic stress disorder. Mol Psychiatry. 2008;13(1):74–83.

    Article  CAS  PubMed  Google Scholar 

  15. Hendrickson RC, Raskind MA, Millard SP, et al. Evidence for altered brain reactivity to norepinephrine in Veterans with a history of traumatic stress. Neurobiology of stress. 2018;8:103–11.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Naegeli C, Zeffiro T, Piccirelli M, et al. Locus coeruleus activity mediates hyperresponsiveness in posttraumatic stress disorder. Biol Psychiat. 2018;83(3):254–62.

    Article  PubMed  Google Scholar 

  17. Dell’Oste V, Fantasia S, Gravina D, et al. Metabolic and inflammatory response in posttraumatic stress disorder (PTSD): a systematic review on peripheral neuroimmune biomarkers. Int J Environ Res Public Health. 2023;20(4):2937.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Al Jowf GI, Ahmed ZT, Reijnders RA, et al. To predict, prevent, and manage posttraumatic stress disorder (PTSD): a review of pathophysiology, treatment, and biomarkers. Int J Mol Sci. 2023;24(6):5238.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Koenen KC, Ratanatharathorn A, Ng L, et al. Posttraumatic stress disorder in the World Mental Health Surveys. Psychol Med. 2017;47(13):2260–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Gautier L, Cope L, Bolstad BM, et al. affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20(3):307–15.

    Article  CAS  PubMed  Google Scholar 

  21. Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Altermann E, Klaenhammer TR. PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. BMC Genomics. 2005;6:60.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Yu G, Wang LG, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Walter W, Sánchez-Cabo F, Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics (Oxford, England). 2015;31(17):2912–4.

    CAS  PubMed  Google Scholar 

  26. Szklarczyk D, Gable AL, Nastou KC, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12.

    Article  CAS  PubMed  Google Scholar 

  27. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chin CH, Chen SH, Wu HH, et al. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Mi H, Muruganujan A, Huang X, et al. Protocol Update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat Protocols. 2019;14(3):703–21.

    Article  CAS  PubMed  Google Scholar 

  30. Uhlen M, Karlsson MJ, Zhong W, et al. A genome-wide transcriptomic analysis of protein-coding genes in human blood cells. Science (New York, NY). 2019;366(6472).

  31. Stein MB, Jain S, Giacino JT, et al. Risk of posttraumatic stress disorder and major depression in civilian patients after mild traumatic brain injury: a TRACK-TBI study. JAMA Psychiat. 2019;76(3):249–58.

    Article  Google Scholar 

  32. Roehr B. American Psychiatric Association explains DSM-5. BMJ (Clinical Research ed). 2013;346: f3591.

    PubMed  Google Scholar 

  33. Dewar M, Paradis A, Fortin CA. Identifying trajectories and predictors of response to psychotherapy for posttraumatic stress disorder in adults: a systematic review of literature. Can J Psychiatry Revue canadienne de psychiatrie. 2020;65(2):71–86.

    PubMed  Google Scholar 

  34. Williams T, Phillips NJ, Stein DJ, et al. Pharmacotherapy for post traumatic stress disorder (PTSD). Cochrane Database Syst Rev. 2022;3(3):CD002795.

    PubMed  Google Scholar 

  35. Stein DJ, Ipser JC, Seedat S. Pharmacotherapy for post traumatic stress disorder (PTSD). Cochrane Database Syst Rev. 2006;2006(1):CD002795.

    PubMed  PubMed Central  Google Scholar 

  36. Aspesi D, Pinna G. Animal models of posttraumatic stress disorder and novel treatment targets. Behav Pharmacol. 2019;30(2 and 3-Spec Issue):130–50.

    Article  CAS  PubMed  Google Scholar 

  37. Stein DJ, McLaughlin KA, Koenen KC, et al. DSM-5 and ICD-11 definitions of posttraumatic stress disorder: investigating “narrow” and “broad” approaches. Depress Anxiety. 2014;31(6):494–505.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Amstadter AB, Nugent NR, Yang BZ, et al. Corticotrophin-releasing hormone type 1 receptor gene (CRHR1) variants predict posttraumatic stress disorder onset and course in pediatric injury patients. Dis Markers. 2011;30(2–3):89–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Schmidt U, Kaltwasser SF, Wotjak CT. Biomarkers in posttraumatic stress disorder: overview and implications for future research. Dis Markers. 2013;35(1):43–54.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Michopoulos V, Norrholm SD, Jovanovic T. Diagnostic biomarkers for posttraumatic stress disorder: promising horizons from translational neuroscience research. Biol Psychiat. 2015;78(5):344–53.

    Article  CAS  PubMed  Google Scholar 

  41. Hu YW, Wu SG, Zhao JJ, et al. VNN1 promotes atherosclerosis progression in apoE-/- mice fed a high-fat/high-cholesterol diet. J Lipid Res. 2016;57(8):1398–411.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Jansen PA, Kamsteeg M, Rodijk-Olthuis D, et al. Expression of the vanin gene family in normal and inflamed human skin: induction by proinflammatory cytokines. J Invest Dermatol. 2009;129(9):2167–74.

    Article  CAS  PubMed  Google Scholar 

  43. Gensollen T, Bourges C, Rihet P, et al. Functional polymorphisms in the regulatory regions of the VNN1 gene are associated with susceptibility to inflammatory bowel diseases. Inflamm Bowel Dis. 2013;19(11):2315–25.

    Article  PubMed  Google Scholar 

  44. Berruyer C, Pouyet L, Millet V, et al. Vanin-1 licenses inflammatory mediator production by gut epithelial cells and controls colitis by antagonizing peroxisome proliferator-activated receptor gamma activity. J Exp Med. 2006;203(13):2817–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Millet V, Gensollen T, Maltese M, et al. Harnessing the Vnn1 pantetheinase pathway boosts short chain fatty acids production and mucosal protection in colitis. Gut. 2023;72(6):1115–28.

    Article  CAS  PubMed  Google Scholar 

  46. Naquet P, Pitari G, Duprè S, et al. Role of the Vnn1 pantetheinase in tissue tolerance to stress. Biochem Soc Trans. 2014;42(4):1094–100.

    Article  CAS  PubMed  Google Scholar 

  47. Chen J, Lu H, Wang X, et al. VNN1 contributes to the acute kidney injury-chronic kidney disease transition by promoting cellular senescence via affecting RB1 expression. FASEB J. 2022;36(9): e22472.

    Article  CAS  PubMed  Google Scholar 

  48. Piao YJ, Kim HS, Han W, et al. Transcriptome analysis of SerpinB2-deficient breast tumors provides insight into deciphering SerpinB2-mediated roles in breast cancer progression. BMC Genomics. 2022;23(1):479.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Medcalf RL, Stasinopoulos SJ. The undecided serpin. The ins and outs of plasminogen activator inhibitor type 2. FEBS J. 2005;272(19):4858–67.

    Article  CAS  PubMed  Google Scholar 

  50. Schroder WA, Le TT, Major L, et al. A physiological function of inflammation-associated SerpinB2 is regulation of adaptive immunity. J Immunol (Baltimore, Md: 1950). 2010;184(5):2663–70.

    Article  CAS  Google Scholar 

  51. Harris NLE, Vennin C, Conway JRW, et al. SerpinB2 regulates stromal remodelling and local invasion in pancreatic cancer. Oncogene. 2017;36(30):4288–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Teske NC, Engelen-Lee JY, Dyckhoff-Shen S, et al. The role of plasminogen activator inhibitor-2 in pneumococcal meningitis. Acta Neuropathol Commun. 2022;10(1):155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lien WC, Zhou XR, Liang YJ, et al. Therapeutic potential of nanoceria pretreatment in preventing the development of urological chronic pelvic pain syndrome: immunomodulation via reactive oxygen species scavenging and SerpinB2 downregulation. Bioeng Transl Med. 2023;8(1): e10346.

    Article  CAS  PubMed  Google Scholar 

  54. Missaglia S, Tavian D, Angelini C. ETF dehydrogenase advances in molecular genetics and impact on treatment. Crit Rev Biochem Mol Biol. 2021;56(4):360–72.

    Article  CAS  PubMed  Google Scholar 

  55. Meng L, Wang Z, Ji HF, et al. Causal association evaluation of diabetes with Alzheimer’s disease and genetic analysis of antidiabetic drugs against Alzheimer’s disease. Cell Biosci. 2022;12(1):28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Maurya NS, Kushwaha S, Chawade A, et al. Transcriptome profiling by combined machine learning and statistical R analysis identifies TMEM236 as a potential novel diagnostic biomarker for colorectal cancer. Sci Rep. 2021;11(1):14304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Maiorova V, Mollaev MD, Vikhreva P, et al. Mutated Flt3Lg provides reduced Flt3 recycling compared to wild-type Flt3Lg and retains the specificity of Flt3Lg-based CAR T-cell targeting in AML models. Int J Mol Sci. 2023;24(8):7626.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Barry KC, Hsu J, Broz ML, et al. A natural killer-dendritic cell axis defines checkpoint therapy-responsive tumor microenvironments. Nat Med. 2018;24(8):1178–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Ramos MI, Perez SG, Aarrass S, et al. FMS-related tyrosine kinase 3 ligand (Flt3L)/CD135 axis in rheumatoid arthritis. Arthritis Res Ther. 2013;15(6):R209.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Chen L, Huang Y, Dong B, et al. Low level FLT3LG is a novel poor prognostic biomarker for cervical cancer with immune infiltration. J Inflamm Res. 2022;15:5889–904.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Yang Y, Fan X, Liu Y, et al. Function and inhibition of DYRK1A: emerging roles of treating multiple human diseases. Biochem Pharmacol. 2023;212: 115521.

    Article  CAS  PubMed  Google Scholar 

  62. Delabar JM, Lagarde J, Fructuoso M, et al. Increased plasma DYRK1A with aging may protect against neurodegenerative diseases. Transl Psychiatry. 2023;13(1):111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Islam M, Gor J, Perkins SJ, et al. The concave face of decorin mediates reversible dimerization and collagen binding. J Biol Chem. 2013;288(49):35526–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Parkinson LG, Toro A, Zhao H, et al. Granzyme B mediates both direct and indirect cleavage of extracellular matrix in skin after chronic low-dose ultraviolet light irradiation. Aging Cell. 2015;14(1):67–77.

    Article  CAS  PubMed  Google Scholar 

  65. Zhu J, Li Y, Shen W, et al. Relationships between transforming growth factor-beta1, myostatin, and decorin: implications for skeletal muscle fibrosis. J Biol Chem. 2007;282(35):25852–63.

    Article  CAS  PubMed  Google Scholar 

  66. Seidler DG, Mohamed NA, Bocian C, et al. The role for decorin in delayed-type hypersensitivity. J Immunol (Baltimore, Md: 1950). 2011;187(11):6108–19.

    Article  CAS  Google Scholar 

  67. Stapor P, Wang X, Goveia J, et al. Angiogenesis revisited—role and therapeutic potential of targeting endothelial metabolism. J Cell Sci. 2014;127(Pt 20):4331–41.

    CAS  PubMed  Google Scholar 

  68. Zhang W, Ge Y, Cheng Q, et al. Decorin is a pivotal effector in the extracellular matrix and tumour microenvironment. Oncotarget. 2018;9(4):5480–91.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Aguilera MO, Robledo E, Melani M, et al. FKBP8 is a novel molecule that participates in the regulation of the autophagic pathway. Biochim Biophys Acta. 2022;1869(5): 119212.

    Article  CAS  Google Scholar 

  70. Edlich F, Lücke C. From cell death to viral replication: the diverse functions of the membrane-associated FKBP38. Curr Opin Pharmacol. 2011;11(4):348–53.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was supported by the National Key R&D Program of the Ministry of Science and Technology (2019YFC0119601), Special Program for military Nursing Innovation and cultivation Project (2021HL091) and Young Elite Scientist Sponsorship Program by CAST 2020QNRC001 (to MJ.H.).

Author information

Author notes

  1. Peng Qi and Mengjie Huang have contributed equally to this work and share first authorship.

Authors and Affiliations

  1. Department of Emergency, First Medical Center of Chinese, PLA General Hospital, 28 Fuxing Road, Beijing, 100853, China

    Peng Qi, Xuewen Ren, Yongzhi Zhai & Haiyan Zhu

  2. Department of Nephrology, First Medical Center of Chinese, PLA General Hospital, 28 Fuxing Road, Beijing, 100853, China

    Mengjie Huang

  3. Department of Orthopedics, Fourth Medical Center of Chinese, PLA General Hospital, Beijing, 100853, China

    Chen Qiu

Authors
  1. Peng Qi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mengjie Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xuewen Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yongzhi Zhai
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Chen Qiu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Haiyan Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

PQ and MH carried out bioinformatics and made significant contributions to the manuscript’s drafting. XR analyzed the data. The paper was reviewed by YZ. CQ and HZ provided writing guidance and gave the work a thorough assessment. The final manuscript was reviewed and approved by all writers.

Corresponding authors

Correspondence to Chen Qiu or Haiyan Zhu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qi, P., Huang, M., Ren, X. et al. Identification of potential biomarkers and therapeutic targets related to post-traumatic stress disorder due to traumatic brain injury. Eur J Med Res 29, 44 (2024). https://doi.org/10.1186/s40001-024-01640-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-024-01640-x

Keywords

European Journal of Medical Research

ISSN: 2047-783X

Contact us