Original Research ARTICLE
A Moderate Blast Exposure Results in Dysregulated Gene Network Activity Related to Cell Death, Survival, Structure, and Metabolism
- 1National Institute of Nursing Research, National Institutes of Health, Bethesda, MD, United States
- 2The Henry M. Jackson Foundation for the Advancement of Military Medicine, Bethesda, MD, United States
- 3Wake Forest School of Medicine, Wake Forest University, Winston-Salem, NC, United States
- 4School of Nursing, University of Texas at Austin, Austin, TX, United States
- 5Department of Neurology, University of Texas, Austin, TX, United States
- 6College of Nursing, University of Iowa, Iowa City, IA, United States
- 7Walter Reed Army Institute of Research, Silver Spring, MD, United States
- 8Naval Medical Research Center, Silver Spring, MD, United States
- 9CNRM Co-Director Biomarkers Core, Uniformed Services University of the Health Sciences, Bethesda, MD, United States
Blast exposure is common in military personnel during training and combat operations, yet biological mechanisms related to cell survival and function that coordinate recovery remain poorly understood. This study explored how moderate blast exposure influences gene expression; specifically, gene-network changes following moderate blast exposure. On day 1 (baseline) of a 10-day military training program, blood samples were drawn, and health and demographic information collected. Helmets equipped with bilateral sensors worn throughout training measured overpressure in pounds per square inch (psi). On day 7, some participants experienced moderate blast exposure (peak pressure ≥5 psi). On day 10, 3 days post-exposure, blood was collected and compared to baseline with RNA-sequencing to establish gene expression changes. Based on dysregulation data from RNA-sequencing, followed by top gene networks identified with Ingenuity Pathway Analysis, a subset of genes was validated (NanoString). Five gene networks were dysregulated; specifically, two highly significant networks: (1) Cell Death and Survival (score: 42), including 70 genes, with 50 downregulated and (2) Cell Structure, Function, and Metabolism (score: 41), including 69 genes, with 41 downregulated. Genes related to ubiquitination, including neuronal development and repair: UPF1, RNA Helicase and ATPase (UPF1) was upregulated while UPF3 Regulator of Nonsense Transcripts Homolog B (UPF3B) was downregulated. Genes related to inflammation were upregulated, including AKT serine/threonine kinase 1 (AKT1), a gene coordinating cellular recovery following TBIs. Moderate blast exposure induced significant gene expression changes including gene networks involved in (1) cell death and survival and (2) cellular development and function. The present findings may have implications for understanding blast exposure pathology and subsequent recovery efforts.
When an individual is in close proximity to a blast, the resulting overpressure (i.e., shock wave) can cause injury to the brain and/or body (1). The increased use of improvised explosives, sophisticated weaponry, and explosive entry techniques has led to increased risk of blast exposure. Specifically, in military personnel who deployed to recent conflicts of Operation Iraqi Freedom (3) and Operation Enduring Freedom (OEF), an estimated 300,000 service members were exposed to at least one blast from adversary attack (2), and blast overpressure from firing weapons is increasing commensurate with increases in weaponry power. High intensity blast exposure events can damage connective tissues, including the central nervous system, resulting in cerebrovascular damage and blood-brain barrier disruption. Significant overpressure can result in tearing of the long axons of neurons (diffuse axonal injury) leading to the associated deficits and comorbidities of a traumatic brain injury (TBI) (4, 5). Although there is evidence suggesting blast-induced TBI (biTBI) has distinct features from blunt-force or penetrating TBI (6), it is difficult to evaluate the consequences of blast in isolation using human subjects because there is often concomitant blunt force or penetrating TBI when objects are propelled and contact the skull (e.g., shrapnel) or the individual is thrown. These challenges contribute to the relatively poor understanding of the pathophysiologic responses to blast and lack of therapies to treat blast-exposed individuals. Moreover, the response and subsequent recovery from blast exposure represent an important line of research, which remains to be further explored and may elucidate the biological mechanisms associated with blast.
Differential gene expression is reported in a small number of clinical TBI studies (7–10), with few studies relevant to blast TBI (11–13). Gene expression regulation is imperative to appropriate cellular response to external mechanical, environmental, or biological stimuli, and the nuclear factor kappa-light-chain enhancer of activated B cells (NF-κB) complex is a main transcription factor of these adaptive gene expression changes (14). More specifically, the NF-κB complex is a transcription factor central to numerous cellular pathways influencing cell survival and proliferation, including inflammatory and immune responses, gene activation, and ubiquitination (14). Animal models demonstrate that the NF-κB complex regulates the innate immune response through upregulation of proinflammatory cytokines including tumor necrosis factor (TNF) (15), interleukin 1 (IL-1) (16), and interleukin 6 (IL-6) (17). In addition, mutations and epigenetic changes within the NF-κB pathway have been linked to immune and inflammatory diseases (18). Cytokines are among a number of factors that may activate NF-κB. NF-κB becomes activated when ubiquitin degrades its inhibitory protein, IκK, freeing NF-κB to enter the nucleus and activate gene transcription (17, 19). Study of gene expression changes following blast exposure may elucidate some of these complexities surrounding the roles and relationships of ubiquitin and inflammatory cytokines following blast exposure.
Within clinical studies of TBI, changes in the NF-κB network are reported in a limited number of studies (7–10); however, they have not yet been examined in biTBI. Preclinical studies of blast exposures have demonstrated altered gene expression, including cognitive impairment (20, 21) and immune function (22). Recent work in military training that involves personnel exposure to blast has demonstrated that ubiquitin carboxy-terminal hydrolase-1 (UCH-L1) is weakly correlated with repeated exposure to low-level blast (12), which is consistent with previous work in TBI (23) and blast exposure (11, 24). In particular, Heinzelmann et al. (11) found protein ubiquitination genes (associated with neuronal recovery, central regulator in IPA) to be downregulated in military personnel with chronic symptoms following blast head injury. UCH-L1 is predominately expressed in the neurons and neuroendocrine cells within the brain (25, 26) and is an enzyme responsible for protein degradation, thus providing a role in ubiquitin stability within neurons and maintaining neuronal health (27). In animal models, a mutation in the UCH-L1 gene causing a truncated protein is associated with neurodegeneration, likely due to the buildup of ubiquitin and subsequent lack of protein clearance (28). Given this limited number of clinical studies, this study sought to further examine differential gene expression pathways in a blast-exposed population. The purpose of this study was to examine gene networks involving cell death and survival as well as cell structure, function, and metabolism to investigate the role of these networks specific to biTBI.
Materials and Methods
To address the gaps in knowledge surrounding the consequences of exposure to isolated blast, a unique cohort of military personnel engaged in training on advanced techniques for breaching buildings with controlled explosives was utilized. The breaching activities were conducted under close supervision and with personal protective equipment and established safety procedures, eliminating the chance of concomitant blunt-force or penetrating TBI. Moreover, recruiting from a training environment (as opposed to real-world combat) facilitated accurate measurement of isolated blast exposures using helmets equipped with pressure sensors (see “Blast Measurement”). This novel sampling also facilitated a collection of baseline data, including pre-exposure blood draws to support assessment of gene expression changes after blast. During the 2 week training program, some participants (n = 29) experienced a moderate blast exposure with peak pressure exceeding 5 pounds per square inch (psi), which exceeded the training range limit of 4 psi and was more than 200% greater than typical exposures measured in such training [e.g., Carr et al. (12)]. The moderate blast exposure was an isolated event, and blast exposures remained ≤2 psi on all other training days. These 29 cases were studied for gene expression changes related to cell death and survival as well as cell structure, function, and metabolism from training day 1 to 10. Unbiased RNA-sequencing (RNA-seq) was used to detect dysregulated genes (13). Ingenuity pathway analysis (29) of dysregulated genes was used to identify gene networks, two of which were validated in the present study using NanoString's nCounter system.
All study protocols were reviewed and approved by the Institutional Review Boards (IRBs) at the Naval Medical Research Center and Walter Reed Army Institute of Research (NMRC#2011.0002; WRAIR#1796) as described in a past publication (12). Prior to study participation, each participant provided informed consent. The parent study from which the present study is drawn was comprised of (N = 108) male active-duty military service members who were engaged in 2 week blast training programs (as either a student or instructor). The goal of the course was to teach advanced techniques for explosive breaching, a tactic used to gain access into secured structures. All participants provided demographic and health history data at baseline, as well as blood samples. For the present study, participants (n = 29) examined were those who experienced a moderate blast exposure (≥5 psi). These 29 individuals, who provided blood samples at the end of training (day 10), were used in the present study to examine gene expression changes from baseline to 3 days post-moderate blast exposure.
Self-reported data were provided by participants at baseline included demographic, health, and blast-history information. Demographic data included age, military rank, and educational status; health information collected included smoking status and history of TBI (Table 1). Previous blast exposure data were also obtained through self-reports on how many blast exposures had been experienced during breaching and artillery fires using the following ordinal scale: 0, 1–9, 10–39, 40–99, 100–199, 200–399, and 400+ blast exposures. Details regarding the surveys used to collect data have been previously described (12).
Objective blast data were collected using standard Army combat helmets equipped with bilateral sensors capable of measuring blast parameters greater than a threshold of 0.4 psi on either sensor. Helmets were worn throughout training and the average of the right and left sensors was used as data to approximate levels of explosive blast each participant experienced. The sensitivity of the sensors is based on the technological specifications of the device itself (micro Data Acquisition System, μDAS; Applied Research Associates, Inc., Albuquerque, NM) as well as considerations for signal-to-noise ratios and effects on data interpretation.
Blood samples were collected at baseline and at the end of 2 week training, which was 3 days after moderate blast. Blood was collected in PAXgene tubes and stored in a −80°C freezer until the time of batch processing.
Random fragmentation of complementary deoxyribonucleic acid (cDNA) followed by 5′ and 3′ adapter ligation was used to create a cDNA library. Average fragment length was 150–170 bp. RNA integrity was assessed using Aligent Technologies 2100 Bioanalyzer and the mean value was 8.9 with standard error of 0.05. Samples from 29 participants on day 1 and 10 were sequenced for mRNA using the Illumina HiSeq-2500 Next Generation Sequencing system (Illumina Inc., San Diego, CA). Using this system, we performed RNA-seq to read paired-ends and read 101 bases per each end. Sequencing data used in the study were deposited in the Gene Expression Omnibus (GEO) with GEO ID GSE89866.
Ingenuity Pathway Analysis
Dysregulated genes were further explored using IPA software (build version 389077M, content version 27821452, released 2016-06-14, Qiagen, Redwood City, CA). Two pathways of interest were identified (see section Results for details and Figures 1, 2).
Figure 1. Network shows dysregulated cell death and survival pathway following moderate blast. Genes described in the text included: UPF1, UPF3B, ARRB1, ZBTB7B, flt3, HIBCH, RPL6, RPL35, MRPL1, MRPL3, MRPL36, and MRPL50. Red indicates increased measurement; green represents decreased measurement; with increased color saturation representing more extreme measurement in dataset. Solid lines represent direct interactions, non-targeting interactions, or correlations between chemicals, proteins, or RNA. Dotted lines represent indirect interactions. Purple lines denote the commonly found relationship among the two merged networks. Arrowed lines represent activation, causation, expression, localization, membership, modification, molecular cleavage, phosphorylation, protein-DNA interactions, protein-RNA interactions, regulation of binding, transcription. Shapes represent molecule type (double circle, complex/group; square, cytokine; diamond, enzyme; inverted triangle, kinase; circle, other; triangle, phosphatase; oval, transcription regulator; trapezoid, transporter). The network was generated through the use of IPA (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis).
Figure 2. Network shows dysregulated structure, function, and development pathway following moderate blast. Genes described in the text include: TRIP12, NAE1, AKT1, MBIP, COA5, and EPB41L3. Red indicates increased measurement; green represents decreased measurement; with increased color saturation representing more extreme measurement in dataset. Solid lines represent direct interactions, non-targeting interactions, or correlations between chemicals, proteins, or RNA. Dotted lines represent indirect interactions. Purple lines denote the commonly found relationship among the two merged networks. Arrowed lines represent activation, causation, expression, localization, membership, modification, molecular cleavage, phosphorylation, protein-DNA interactions, protein-RNA interactions, regulation of binding, transcription. Shapes represent molecule type (double circle, complex/group; square, cytokine; diamond, enzyme; inverted triangle, kinase; circle, other; triangle, phosphatase; oval, transcription regulator; trapezoid, transporter). The network was generated through the use of IPA (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis).
A subset of genes examined in RNA-seq data were selected to validate gene expression changes using a direct digital detection system (NanoString Technologies, Seattle, WA). In selecting genes to validate, the extent of dysregulation, biological plausibility, and the position of the protein within the IPA pathway diagrams were considered. Two pathways were identified: one focused on cell death and survival and another focused on basic structure, function, and development. A panel was designed for each pathway to include 50 markers of interest, plus a total of nine reference genes for data normalization (Tables 2, 3). Reference genes met the following criteria: (1) not dysregulated in the RNA-seq data for the same samples; (2) not clearly implicated in traumatic brain injury, blast exposure, or a similar condition; and (4) no published evidence that this was an unstable reference gene in human blood. Probes for the 50 genes of interest and the reference genes were designed and manufactured by NanoString Technologies. NanoString was used to determine the mean copy number of each mRNA probe of interest based on manufacturer's protocol. The standard manufacturer protocol was followed for sample preparation, hybridization, and detection.
The Statistical Package for the Social Sciences (SPSS; version 22; IBM Corporation, Armonk, NY) and NanoString's nSolver Analysis Software (version 3; NanoString Technologies, Seattle, WA) were usedfor analyses.
The moderate blast exposed cases (n = 29) met QC criteria based on the RNA Integrity Number (RIN) and were subsequently sequenced. In total, between 52.5 and 75.5 million read counts were completed for each sample; in 94.95% of base calls, an accuracy of at least Q30 was achieved. To establish bioinformatics quality control (QC), FastQC (version 0.11.5, Babraham Bioinformatics, Cambridgeshire, UK) was used. Data were aligned to a reference genome (hg19) using an open-source aligner, STAR (version 2.5) (30). To count the number of reads mapped to genes, HTSeq software was used (version 0.6.1p1) (31). DESeq2 (version 1.12.3) (32) was used to identify differentially expressed genes, with cutoff value of false discovery rate (FDR) of 0.05. This cut off value of the false discovery rate is considered to be a conservative balance in minimizing both false positives (type I errors) and false negatives (type II errors) (33, 34).
NanoString Validation Analysis
RNA sequencing data were validated with 50 genes from each network (100 genes total), plus reference genes, using a custom NanoString codeset. Raw data were analyzed using nSolver 3.0 digital analyzer software using standard settings and quality control parameters. Raw data were normalized against reference genes. Gene expression differences before and after blast exposure were measured using t-tests, with the Benjamini-Yekutieli FDR method for multiple comparison corrections to conservatively balance the minimization of type I and type II errors (35). Log2 fold changes and adjusted p-values were calculated for samples before and after blast exposure, with statistical significance defined at the level of p < 0.05, demonstrating congruency with RNA sequencing findings.
Participants in the study were male military service members with a mean age of 30.42 and a mean length of service of 9.94 years (Table 1). Almost half of participants (46.3%) had a history of >40 prior blast exposures. No significant differences based on demographic information were noted among the cohort [described in previous report by Gill et al. (13)].
Results of the RNA-seq analysis demonstrated significantly dysregulated gene activity changes following moderate blast exposure. These genes were entered into IPA software (29). IPA analysis identified five significantly dysregulated networks before and after moderate blast exposure (Table 4). Two sets of two networks shared overlapping functions and were subsequently merged to form two networks of interest in the present study: (1) Cell Death and Survival and (2) Cell Structure, Function, and Metabolism (Table 4). NanoString analysis confirmed 32 significantly differentially expressed genes in the Cell Death and Survival network (p < 0.05) and 35 significantly differentially expressed genes in the Cell Structure, Function, and Metabolism network (p < 0.05), validating differential expression of these two gene networks following blast exposure. Here, we report on the most significant gene network activity changes in the Cell Death and Survival and the Cell Structure, Function, and Metabolism networks (Figures 1, 2).
Cell Death and Survival Network
One merged pathway centered on cell death and survival. Genes identified in this pathway are defined in Table 2. This pathway was comprised of genes implicated in apoptosis, necrosis, autophagy, mitophagy, ferroptosis, survival, regeneration, and recovery, with an IPA score of 42 (Figure 1). Significantly dysregulated genes include UPF1, RNA helicase and ATPase (UPF1), UPF3 regulator of nonsense transcripts homolog B (UPF3B), arrestin β1 (ARRB1), zinc finger and BTB domain containing 7B (ZBTB7B), fms related tyrosine kinase 3 (flt3), 3-hydroxyisobutyryl-CoA hydrolase (HIBCH), ribosomal proteins (RPL6, -L35), and mitochondrial ribosomal proteins (MRPL1, -L3, -L36, and -L50). Major hubs within this network include growth factor receptor bound protein 2 (GRB2) and staufen double-stranded RNA binding protein 1 (STAU1).
Cell Structure, Function, and Metabolism Network
The second merged pathway focused on development, metabolism, and cell structure/function. Genes identified in this pathway are defined in Table 3. This pathway consisted of genes involved in cytoskeleton, organelles, cellular metabolism, lipid metabolism, heat shock, cell motion, cell growth, and differentiation, with an IPA score of 41 (Figure 2). Significant genes within this network include tripartite motif containing 12 (TRIP12), NEDD8 activating enzyme E1 subunit 1 (NAE1), cytochrome C oxidase assembly factor 5 (COA5), and erythrocyte membrane protein band 4.1-like 3 (EPB41L3). AKT serine/threonine kinase 1 (AKT1), amyloid precursor protein (APP), and MAP3K12 binding inhibitory protein 1 (MBIP) are the major hubs in this network.
Activity changes in two gene networks were found after moderate blast exposure in military personnel engaged in training. Differentially regulated networks after blast included (1) Cell Death and Survival (Figure 1) and (2) Cell Structure, Function, and Development (Figure 2). Genes within these two networks relate to ubiquitination, nonsense mediated decay (NMD), apoptosis, as well as activity related to ribosomes, mitochondria, and inflammation. Findings provide novel insights for understanding the biological changes that occur following blast, which for some individuals, may result in biological changes that increase their risk for neurological or behavioral symptoms and deficits. These findings may ultimately contribute to characterizing the cellular mechanisms of blast exposure to improve diagnosis, monitoring, and prognosis of military personnel exposed to blast.
In this study, genes related to ubiquitination were increased in activity following blast exposure, including tripartite motif containing 12 (TRIP12), which is an E3 ubiquitin-protein ligase involved in ubiquitin fusion degradation. Protein ubiquitination initiates the removal of oxidized and misfolded proteins following injury, and its processes can protect neurons from reactive oxidative species (ROS) that accumulate following blast exposure in pre-clinical models (36). Our findings provide further evidence of increased UCH-L1, the primary protein for ubiquitination, following repeated low-level blast (12). This finding suggests that there may also be overlap with the biological mechanisms related to recovery from TBIs in civilians, as UCHL1 increases are one of the most often reported changes following a TBI (37, 38). However, these findings are in contrast to another previous report in which the activity of genes related to ubiquitin were lower in activity in military personnel with TBIs, with many related to blast exposures, and chronic symptoms (11). Therefore, it may be that ubiquitin activity is critical to acute recovery from biTBIs, and that in some individuals, there is a reduction in activity that may place them at higher risk for chronic symptoms. In support of this, pre-clinical studies show that reductions or inactivation of ubiquitin activity results in poor outcomes, including behavioral deficits, possibly indicating long-term neurodegenerative processes (39).
Genes that may relate to neuronal recovery were altered in activity following a moderate blast. Specifically, we report gene activity changes within the NMD pathway, including UPF1 and UPF3B, which are responsible for neuronal specific cell development and repair through a reciprocal pattern of activity (40). Previous studies show an interaction in the activity of these two genes, such that when one gene is less active, the other gene will compensate, preserving the activity of this network; our findings mirror this. Here we report that UPF1 was increased in activity, whereas UPF3B was downregulated. These findings suggest that in response to the blast, injury mechanisms may have been initiated (inflammation, aberrant cellular formation, and cell death), and this initiation may result in an upregulation of UPF1, in an effort to preserve the activity of the NMD pathway. Subsequently, here we report that the expression of UPF3B is suppressed, hindering possible detrimental neurological effects. These findings suggest complex gene-activity changes following blast exposure that may be occurring to promote recovery; additional studies are needed to increase understanding of the temporal relationship of these changes and their relation to neuronal recovery.
A downregulated gene within the structure, function, and development pathway was NAE1 (NEDD8 Activating Enzyme E1 Subunit 1), a protein associated with the neddylation pathway. Vogl et al. showed that neddylation was a critical regulator of dendritic spine development, reporting that in NAE1 knockout mice, there were cognitive deficits as well as synaptic and neurotransmitter impairments (41). The down regulation observed in our military population could suggest similarly that exposure to blast hinders the neddylation pathway, which might suggest a marker of injury resulting directly from blast exposure. Additionally, recent in vitro work suggests Il-1β may inhibit NEDD8 and neddylation in conjunction with increased ubiquitination; while activation of NEDD8 downregulates the nuclear factor kappa light-chain enhancer of activated B cells (NF-κB) pathway (42). This finding is of interest because we also found that genes within the NF-κB network show activity changes, with most genes becoming more active. The NF-κB network is a dominant activator of the immune system following TBI, and this activity is essential because it initiates secondary injury mechanisms required for neuronal recovery. However, if activity of this pathway is too high, or too long-lasting, it can be detrimental to neuronal recovery (43). One such gene is ARRB1(arrestin β1), which is increased following blast exposure. This gene has been reported to play a role in the beta-adrenergic receptor kinase (BARK) mediated desensitization of beta-adrenergic receptors. In TBI patients, catecholamine surge after injury has been linked to immunosuppression and greater mortality risk that is reversed through β-blocker treatment (44). ZBTB7B (zinc finger and BTB domain containing 7B) is also upregulated after blast and linked to reductions in CD8-cytotoxic activity (45), which could be a mechanism to prevent further cellular damage after blast injury.
In addition to the activation of genes within the NF-κB network, AKT1 is another immune-related gene with increased activity. AKT1 is a hub with 14 connections in the structure, function, and development network. AKT1 encodes for a serine-threonine protein kinase (AKT1), which is known to regulate a vast number of cellular processes including neuronal survival, glucose uptake, protein and fatty acid synthesis, cell proliferation, and the previously mentioned role in apoptosis (46, 47). Additionally, AKT1 may function in the inflammatory response as an upstream activator of NF-κB (48). Interestingly, in this population, significantly elevated levels of the cytokines tumor necrosis factor alpha (TNFα) and interleukin 6 (IL-6) have been reported during the acute period following moderate blast (49). This finding is relevant because NF-κB is recognized as a master regulator of cytokines including TNFα and IL-6 (47, 50). The NF-κB pathway has been implicated in the regulation of proinflammatory cytokines during meningitis (51) and in blood-brain barrier permeability (52). Moreover, the NF-κB pathway has been found to be dysregulated in clinical studies of acute and subacute TBI (7–10). Upregulation of AKT1 in this sample suggests activation of the NF-κB pathway; a finding that supports these prior studies, although the specific role of AKT1 in blast effects on the central nervous system remains to be examined.
Also showing increased activity in this study was the Flt3 (dimer), which encodes for a receptor tyrosine kinase and is also related to NF-κB pathway. Flt3 is implicated in multiple signaling pathways including regulation of the proliferation and survival of hematopoietic cells, which ultimately relates to the number of intermediate monocytes (53). Because intermediate monocytes promote production of inflammatory cytokines within the NF-κb network, including TNF-α and Il-1β (54), these findings suggest the possibility of a pro-inflammatory response through increased production of intermediate monocytes.
Our findings of changes in activity of apoptosis-related genes following blast are of interest because preclinical models show that blast exposure results in astrocytic and microglial activation, oxidative stress, axonal and vascular damage, and inflammation, which ultimately contribute to programmed cell-death (55–58). Specifically, we found activation of caspase complexes, a family of cysteine-dependent proteases, which have been previously associated with neuronal and oligo-dendroglial cell death in both pre-clinical and human brain injuries (59). Otherwise referred to as apoptosis executioners, caspase-3 and -7 are both indirectly activated by MBIP, a major hub of the cell structure, function, and metabolism network. Increased expression in caspase-3 and -7 complexes have also been previously linked to TBIs in pre-clinical models (60, 61) and to mortality in patients with severe TBIs (62). We also report increased activity of other apoptosis genes following blast including EPB41L3, or erythrocyte membrane protein band 4.1-like 3, and EPB41L3, a tumor suppressor gene strongly expressed in the brain that promotes apoptotic pathways and inhibits cellular proliferation (63). These findings suggest that moderate blast results in expression of apoptosis inducing genes, and that mitigating these activities may be protective.
Lastly, several mitochondrial genes and genes connected with the mitochondrial gene network are dysregulated, including COA5, HIBCH, RPL6, RPL35, as well as mitochondrial ribosomal genes MRPL50, MRPL1, MRPL3, and MRPL46. Although the function of mitochondria is not yet well-understood in blast exposures, it is worth noting that mitochondrial dysfunction has been implicated in preclinical TBI pathology. Previous studies have indicated that following TBI an influx of intracellular calcium leads to disruption of the mitochondrial membrane potential, impairing ATP production and creating ROS, activating cell death pathways and leading to neuronal damage associated with cognitive impairments (64, 65). The biological mechanisms specific to blast effects on the central nervous system in the context of mitochondrial genes are not yet known.
Alterations in these gene expression pathways may have translational implications for blast-related neuropathology if replicated. In this military training population, dysregulation in gene expression pathways related to ubiquitination, NMD, apoptosis, inflammation, and ribosomal and mitochondrial activity were observed, suggesting a role for these pathways following acute blast exposures. Examination of related, downstream proteins may indicate potential diagnostic or prognostic biomarkers. Mapping temporal changes in gene expression, downstream proteins, and symptomology may shed light on the role of these pathways in underlying neuropathological processes and clinical outcomes.
Although these initial findings provide novel insights into gene-activity changes following blast, this study has a number of limitations. First, the secondary data analysis in this study precluded the comparison of control personnel who were engaged in blast training, but did not sustain a moderate blast; therefore, gene expression changes related to normal daily training activities or other activities of daily living such as circadian rhythm or diet cannot be determined. In addition, a sub-portion of the moderate blast cohort were experienced trainers with previous blast exposures, and thus were not naïve to blast. However, this sample represents typical training cohorts. Second, blast exposure may affect cell types throughout the body in addition to brain-related pathways. Preclinical models demonstrate that the CNS is affected by blast exposure (66, 67); however, at this time translating this information to clinical populations presents difficulty as accessing the CNS requires intensive, invasive procedures. Third, due to the exploratory nature of this work, clinical symptomology measures were not collected.
These limitations represent important considerations for future work. Additional studies in clinical blast exposure will need to collect samples from unaffected control cohorts, not only a priori matched samples, in order to differentiate possible causative and confounding agents such as blast exposure vs. training effects. We are addressing this design limitation in future studies. Next, considering the limitation in directly studying the CNS in clinical populations, an exciting future direction for this work would be the ability to measure brain-specific peripheral biomarkers, in order to differentiate brain-related pathways from other cellular processes influenced by blast exposures. Finally, the significant gene expression changes found in this study warrant further research into possible symptomology that may be experienced following blast exposure utilizing validated clinical measures and mapping these symptoms to biomarker changes. Notably, this study identified altered gene expression pathways, which may be the focus of more specific biological mechanisms significant to blast exposure in future studies.
The initial findings reported here show that there are robust gene activity changes following a moderate blast exposure in a sample of military personnel. Notably, this study's findings are important to distinguish the effects of blast exposure without the co-occurring impact of blunt-force injuries that take place in combat stations. Findings from this study suggest that additional studies are needed to examine gene-activity related to blast exposure, in order to examine the impact of previous blast exposures on biomarker changes, relationships to neurological impacts, and subsequent clinical symptoms. Thus, these findings provide novel insights into gene network dysregulation observed following objectively measured blast exposure that warrant future clinical studies to advance the understanding of neuropathology related to blast exposure.
Data Availability Statement
Sequencing data used in the study were deposited in the Gene Expression Omnibus (GEO) with GEO ID GSE89866.
This study protocol was reviewed and approved by the Institutional Review Boards at the Naval Medical Research Center and Walter Reed Army Institute of Research (NMRC#2011.0002; WRAIR#1796). All subjects gave written informed consent prior to participation in the study.
KE, JG, NO, and H-SK contributed to the conception or design of the study. KE, VM, SY, Y-EC, and CL contributed to lab analysis, interpretation of the data, and drafting of the work. KE, SY, Y-EC, CL, KD, WC, PW, SA, ML, AY, AT, and JG contributed to the collection, interpretation of the data, and drafting and editing of the work. All authors contributed to critical revision of the manuscript, read and approved the submitted version, and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
This work was supported by the National Institute of Nursing Research Intramural Program, the US Army Medical Research and Materiel Command, and the US Navy Bureau of Medicine (NMRC#2011.0002; WRAIR#1796). While producing this article, KE was supported as a Jonas Nurse Leader Scholar.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The opinions and assertions in this manuscript are those of the authors and are not to be construed as official as reflecting true views of the Department of Navy, Department of the Army, Department of Defense, the Uniformed Services University of the Health Sciences, U.S. Government, the Center for Neuroscience and Regenerative Medicine or any other agency of the U.S. government. This study was funded by National Institutes of Health, Intramural Department of Research, and the US Army Medical Research and Material Command and the US Navy Bureau of Medicine (NMRC#2011.0002; WRAIR#1796). The study protocol was approved by the Walter Reed Army Institute of Research/Naval Medical Research Center Institutional Review Boards in compliance with AR 70-25 and all applicable Federal regulations governing the protection of human subjects. Some of the authors are military service members or employees of the U.S. Government. This work was prepared as part of their official duties. Title 17 U.S.C. §105 provides that “Copyright protection under this title is not available for any work of the United States Government.” Title 17 U.S.C. §101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person's official duties. Material has been reviewed by the Walter Reed Army Institute of Research. There is no objection to its presentation and/or publication.
We thank the leadership and Soldiers of the military units studied for their service to our nation and their participation in the study. The authors acknowledge the helpful comments provided by Dr. Kathleen Valentine, Dr. Julia Eggert, and Dr. Mary Beth Steck of the Clemson University Healthcare Genetics Program, and by Dr. Sheila Alexander of the University of Pittsburgh School of Nursing. The authors thank Dr. Jim McDonell of Clemson University for his statistical advice.
2. Tenielien T, Jaycox L. Invisible Wounds of War: Psychological and Cognitive Injuries, Their Consequences, and Services to Assist Recovery. Santa Monica, CA: RAND Center for Military Health Policy Research (2008). doi: 10.1037/e527612010-001
3. Meerlo P, Sgoifo A, Suchecki D. Restricted and disrupted sleep: effects on autonomic function, neuroendocrine stress systems and stress responsivity. Sleep Med Rev. (2008) 12:197–210. doi: 10.1016/j.smrv.2007.07.007
4. Yurgil KA, Barkauskas DA, Vasterling JJ, Nievergelt CM, Larson GE, Schork NJ, et al. Association between traumatic brain injury and risk of posttraumatic stress disorder in active-duty Marines. JAMA Psychiatry. (2014) 71:149–57. doi: 10.1001/jamapsychiatry.2013.3080
5. Mac Donald CL, Johnson AM, Wierzechowski L, Kassner E, Stewart T, Nelson EC, et al. Outcome trends after US military concussive traumatic brain injury. J Neurotrauma. (2016) 34:2206–9. doi: 10.1089/neu.2016.4434
6. Courtney A, Courtney M. The complexity of biomechanics causing primary blast-induced traumatic brain injury: a review of potential mechanisms. Front Neurol. (2015) 6:221. doi: 10.3389/fneur.2015.00221
7. Livingston WS, Gill J, Cota M, Olivera A, O'Keefe JL, Martin C, et al. Differential Gene Expression Associated with Meningeal Injury in Acute Mild TBI. J Neurotrauma (2016) 34:853–86. doi: 10.1089/neu.2016.4479
8. Cho YE, Latour LL, Kim H, Turtzo LC, Olivera A, Livingston WS, et al. Older age results in differential gene expression after mild traumatic brain injury and is linked to imaging differences at acute follow-up. Front Aging Neurosci. (2016) 8:168. doi: 10.3389/fnagi.2016.00168
9. Gill J, Merchant-Borna K, Lee H, Livingston WS, Olivera A, Cashion A, et al. Sports-related concussion results in differential expression of nuclear factor-kappaB pathway genes in peripheral blood during the acute and subacute periods. J Head Trauma Rehabil. (2016) 31:269–76. doi: 10.1097/HTR.0000000000000191
10. Merchant-Borna K, Lee H, Wang D, Bogner V, van Griensven M, Gill J, et al. Genome-wide changes in peripheral gene expression following sports-related concussion. J Neurotrauma (2016) 33:1576–85. doi: 10.1089/neu.2015.4191
11. Heinzelmann M, Reddy SY, French LM, Wang D, Lee H, Barr T, et al. Military personnel with chronic symptoms following blast traumatic brain injury have differential expression of neuronal recovery and epidermal growth factor receptor genes. Front Neurol. (2014) 5:198. doi: 10.3389/fneur.2014.00198
12. Carr W, Yarnell AM, Ong R, Walilko T, Kamimori GH, da Silva U, et al. Ubiquitin carboxy-terminal hydrolase-l1 as a serum neurotrauma biomarker for exposure to occupational low-level blast. Front Neurol. (2015) 6:49. doi: 10.3389/fneur.2015.00049
13. Gill J, Cashion A, Osier N, Arcurio L, Motamedi V, Dell KC, et al. Moderate blast exposure alters gene expression and levels of amyloid precursor protein. Neurol Genet. (2017) 3:e186. doi: 10.1212/NXG.0000000000000186
15. Bohuslav J, Kravchenko VV, Parry GC, Erlich JH, Gerondakis S, Mackman N, Ulevitch RJ. Regulation of an essential innate immune response by the p50 subunit of NF-kappaB. J Clin Invest. (1998) 102:1645–52. doi: 10.1172/JCI3877
21. Tweedie D, Rachmany L, Rubovitch V, Li Y, Holloway HW, Lehrmann E, et al. Blast traumatic brain injury-induced cognitive deficits are attenuated by preinjury or postinjury treatment with the glucagon-like peptide-1 receptor agonist, exendin-4. Alzheimers Dement. (2016) 12:34–48. doi: 10.1016/j.jalz.2015.07.489
22. Struebing FL, King R, Li Y, Chrenek MA, Lyuboslavsky PN, Sidhu CS, et al. Transcriptional changes in the mouse retina following ocular blast injury: a role for the immune system. J Neurotrauma. (2017) 35:118–29. doi: 10.1089/neu.2017.5104
23. Papa L, Lewis LM, Silvestri S, Falk JL, Giordano P, Brophy GM, et al. Serum levels of ubiquitin C-terminal hydrolase distinguish mild traumatic brain injury from trauma controls and are elevated in mild and moderate traumatic brain injury patients with intracranial lesions and neurosurgical intervention. J Trauma Acute Care Surg. (2012) 72:1335–44. doi: 10.1097/TA.0b013e3182491e3d
24. Tate CM, Wang KK, Eonta S, Zhang Y, Carr W, Tortella FC, et al. Serum brain biomarker level, neurocognitive performance, and self-reported symptom changes in soldiers repeatedly exposed to low-level blast: a breacher pilot study. J Neurotrauma. (2013) 30:1620–30. doi: 10.1089/neu.2012.2683
25. Doran JF, Jackson P, Kynoch PA, Thompson RJ. Isolation of PGP 9.5, a new human neurone-specific protein detected by high-resolution two-dimensional electrophoresis. J Neurochem. (1983) 40:1542–7. doi: 10.1111/j.1471-4159.1983.tb08124.x
27. Osaka H, Wang YL, Takada K, Takizawa S, Setsuie R, Li H, et al. Ubiquitin carboxy-terminal hydrolase L1 binds to and stabilizes monoubiquitin in neuron. Hum Mol Genet. (2003) 12:1945–58. doi: 10.1093/hmg/ddg211
28. Saigoh K, Wang YL, Suh JG, Yamanishi T, Sakai Y, Kiyosawa H, et al. Intragenic deletion in the gene encoding ubiquitin carboxy-terminal hydrolase in gad mice. Nat Genet. (1999) 23:47–51. doi: 10.1038/12647
33. Glickman ME, Rao SR, Schultz MR. False discovery rate control is a recommended alternative to Bonferroni-type adjustments in health studies. J Clin Epidemiol. (2014) 67:850–7. doi: 10.1016/j.jclinepi.2014.03.012
34. Li J, Witten DM, Johnstone IM, Tibshirani R. Normalization, testing, and false discovery rate estimation for RNA-sequencing data. Biostatistics. (2012) 13:523–38. doi: 10.1093/biostatistics/kxr031
36. Kochanek PM, Dixon CE, Shellington DK, Shin SS, Bayir H, Jackson EK, et al. Screening of biochemical and molecular mechanisms of secondary injury and repair in the brain after experimental blast-induced traumatic brain injury in rats. J Neurotrauma. (2013) 30:920–37. doi: 10.1089/neu.2013.2862
38. Diaz-Arrastia R, Wang KK, Papa L, Sorani MD, Yue JK, Puccio AM, et al. Acute biomarkers of traumatic brain injury: relationship between plasma levels of ubiquitin C-terminal hydrolase-L1 and glial fibrillary acidic protein. J Neurotrauma. (2014) 31:19–25. doi: 10.1089/neu.2013.3040
39. Svetlov SI, Prima V, Kirk DR, Gutierrez H, Curley KC, Hayes RL, et al. Morphologic and biochemical characterization of brain injury in a model of controlled blast overpressure exposure. J Trauma. (2010) 69:795–804. doi: 10.1097/TA.0b013e3181bbd885
41. Vogl AM, Brockmann MM, Giusti SA, Maccarrone G, Vercelli CA, Bauder CA, et al. Neddylation inhibition impairs spine development, destabilizes synapses and deteriorates cognition. Nat Neurosci. (2015) 18:239–51. doi: 10.1038/nn.3912
42. Yan F, Guan J, Peng Y, Zheng X. MyD88 NEDDylation negatively regulates MyD88-dependent NF-kappaB signaling through antagonizing its ubiquitination. Biochem Biophys Res Commun. (2017) 482:632–7. doi: 10.1016/j.bbrc.2016.11.084
44. Schroeppel TJ, Fischer PE, Zarzaur BL, Magnotti LJ, Clement LP, Fabian TC, et al. Beta-adrenergic blockade and traumatic brain injury: protective? J Trauma. (2010) 69:776–82. doi: 10.1097/TA.0b013e3181e981b8
45. Wang L, Wildt KF, Castro E, Xiong Y, Feigenbaum L, Tessarollo L, et al. The zinc finger transcription factor Zbtb7b represses CD8-lineage gene expression in peripheral CD4+ T cells. Immunity. (2008) 29:876–87. doi: 10.1016/j.immuni.2008.09.019
48. Lian S, Xia Y, Khoi PN, Ung TT, Yoon HJ, Kim NH, et al. Cadmium induces matrix metalloproteinase-9 expression via ROS-dependent EGFR, NF-small ka, CyrillicB, and AP-1 pathways in human endothelial cells. Toxicology. (2015) 338:104–16. doi: 10.1016/j.tox.2015.10.008
49. Gill J, Motamedi V, Osier N, Dell K, Arcurio L, Carr W, et al. Moderate blast exposure results in increased IL-6 and TNFalpha in peripheral blood. Brain Behav Immun. (2017) 65:90–4. doi: 10.1016/j.bbi.2017.02.015
51. Barichello T, Generoso JS, Simoes LR, Elias SG, Quevedo J. Role of oxidative stress in the pathophysiology of pneumococcal meningitis. Oxid Med Cell Longev. (2013) 2013:371465. doi: 10.1155/2013/371465
52. Merrill JE, Murphy SP. Inflammatory events at the blood brain barrier: regulation of adhesion molecules, cytokines, and chemokines by reactive nitrogen and oxygen species. Brain Behav Immun. (1997) 11:245–63. doi: 10.1006/brbi.1997.0496
53. Zawada AM, Schneider JS, Michel AI, Rogacev KS, Hummel B, Krezdorn N, et al. DNA methylation profiling reveals differences in the 3 human monocyte subsets and identifies uremia to induce DNA methylation changes during differentiation. Epigenetics. (2016) 11:259–72. doi: 10.1080/15592294.2016.1158363
54. Wong KL, Chen W, Balakrishnan T, Toh YX, Fink K, Wong SC. Susceptibility and response of human blood monocyte subsets to primary dengue virus infection. PLoS ONE. (2012) 7:e36435. doi: 10.1371/journal.pone.0036435
55. Saljo A, Mayorga M, Bolouri H, Svensson B, Hamberger A. Mechanisms and pathophysiology of the low-level blast brain injury in animal models. Neuroimage. (2011) 54(Suppl 1):S83–8. doi: 10.1016/j.neuroimage.2010.05.050
57. Agoston DV, Gyorgy A, Eidelman O, Pollard HB. Proteomic biomarkers for blast neurotrauma: targeting cerebral edema, inflammation, and neuronal death cascades. J Neurotrauma. (2009) 26:901–11. doi: 10.1089/neu.2008.0724
58. Goodrich JA, Kim JH, Situ R, Taylor W, Westmoreland T, Du F, et al. Neuronal and glial changes in the brain resulting from explosive blast in an experimental model. Acta Neuropathol Commun. (2016) 4:124. doi: 10.1186/s40478-016-0395-3
60. Larner SF, McKinsey DM, Hayes RL, Wang WKK. Caspase 7: increased expression and activation after traumatic brain injury in rats. J Neurochem. (2005) 94:97–108. doi: 10.1111/j.1471-4159.2005.03172.x
61. Clark RS, Kochanek PM, Watkins SC, Chen M, Dixon CE, Seidberg NA, et al. Caspase-3 mediated neuronal death after traumatic brain injury in rats. J Neurochem. (2000) 74:740–53. doi: 10.1046/j.1471-4159.2000.740740.x
62. Zhang X, Alber S, Watkins SC, Kochanek PM, Marion DW, Graham SH, Clark RS. Proteolysis consistent with activation of caspase-7 after severe traumatic brain injury in humans. J Neurotrauma. (2006) 23:1583–90. doi: 10.1089/neu.2006.23.1583
63. Li X, Zhang Y, Zhang H, Liu X, Gong T, Li M, et al. miRNA-223 promotes gastric cancer invasion and metastasis by targeting tumor suppressor EPB41L3. Mol Cancer Res. (2011) 9:824–33. doi: 10.1158/1541-7786.MCR-10-0529
65. Ohta M, Higashi Y, Yawata T, Kitahara M, Nobumoto A, Ishida E, et al. Attenuation of axonal injury and oxidative stress by edaravone protects against cognitive impairments after traumatic brain injury. Brain Res. (2013) 1490:184–92. doi: 10.1016/j.brainres.2012.09.011
67. Bugay V, Bozdemir E, Vigil FA, Holstein DM, Chun SH, Elliot W, et al. A mouse model of repetitive blast traumatic brain injury reveals post-trauma seizures and increased neuronal excitability. J Neurotrauma. (2019) 37:248–61. doi: 10.1089/neu.2018.6333
Keywords: blast, overpressure, gene expression, RNA-sequencing, NanoString
Citation: Edwards KA, Motamedi V, Osier ND, Kim H-S, Yun S, Cho Y-E, Lai C, Dell KC, Carr W, Walker P, Ahlers S, LoPresti M, Yarnell A, Tschiffley A and Gill JM (2020) A Moderate Blast Exposure Results in Dysregulated Gene Network Activity Related to Cell Death, Survival, Structure, and Metabolism. Front. Neurol. 11:91. doi: 10.3389/fneur.2020.00091
Received: 24 January 2019; Accepted: 27 January 2020;
Published: 26 February 2020.
Edited by:Henrik Zetterberg, University of Gothenburg, Sweden
Reviewed by:Francisco Capani, University of Buenos Aires, Argentina
Nicholas Ashton, University of Gothenburg, Sweden
At least a portion of this work is authored by Katie A. Edwards, Vida Motamedi, Nicole D. Osier, Hyung-Suk Kim, Young-Eun Cho, Chen Lai, Kristine C. Dell, Walter Carr, Peter Walker, Stephen Ahlers, Matthew LoPresti, Angela Yarnell, Anna Tschiffley and Jessica M. Gill on behalf of the U.S. Government and, as regards Katie A. Edwards, Vida Motamedi, Nicole D. Osier, Hyung-Suk Kim, Young-Eun Cho, Chen Lai, Kristine C. Dell, Walter Carr, Peter Walker, Stephen Ahlers, Matthew LoPresti, Angela Yarnell, Anna Tschiffley and Jessica M. Gill and the U.S. Government, is not subject to copyright protection in the United States. Foreign and other copyrights may apply. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katie A. Edwards, email@example.com