Hostname: page-component-78c5997874-dh8gc Total loading time: 0 Render date: 2024-11-05T03:15:44.890Z Has data issue: false hasContentIssue false

ErbB Signaling Pathway Genes Are Differentially Expressed in Monozygotic Twins Discordant for Sports-Related Concussion

Published online by Cambridge University Press:  26 May 2022

Tyler J. Stahl
Affiliation:
University of Rochester School of Arts and Sciences, Rochester, NY, USA
Rohith Palli
Affiliation:
Medical Scientist Training Program, Biophysics, Structural, and Computational Biology Program, Rochester, NY, USA
Jessica M. Gill
Affiliation:
Johns Hopkins University, Baltimore, MD, USA
Elizabeth L. Saionz
Affiliation:
Medical Scientist Training Program, Biophysics, Structural, and Computational Biology Program, Rochester, NY, USA
Viktoria Bogner-Flatz
Affiliation:
Department of Trauma Surgery, Ludwig-Maximilians-University, Munich, Germany
Arokoruba Cheetham-West
Affiliation:
University of Rochester School of Arts and Sciences, Rochester, NY, USA
Kian Merchant-Borna
Affiliation:
Department of Emergency Medicine, University of Rochester School of Medicine and Dentistry, Rochester, NY, USA
Dzung L. Pham
Affiliation:
Center for Neuroscience and Regenerative Medicine, National Institutes of Health, Bethesda, MD, USA
Lachlan Strike
Affiliation:
Queensland Brain Institute, University of Queensland, Brisbane, QLD, Australia
Margaret J. Wright
Affiliation:
Queensland Brain Institute, University of Queensland, Brisbane, QLD, Australia
Paul Thompson
Affiliation:
Institute for Neuroimaging and Informatics, University of Southern California, Los Angeles, CA, USA
Jeffrey J. Bazarian*
Affiliation:
Department of Emergency Medicine, University of Rochester School of Medicine and Dentistry, Rochester, NY, USA
*
Author for correspondence: Jeffrey Bazarian, Email: [email protected]

Abstract

Transcriptional changes involved in neuronal recovery after sports-related concussion (SRC) may be obscured by inter-individual variation in mRNA expression and nonspecific changes related to physical exertion. Using a co-twin study, the objective of this study was to identify important differences in mRNA expression among a single pair of monozygotic (MZ) twins discordant for concussion. A pair of MZ twins were enrolled as part of a larger study of concussion biomarkers among collegiate athletes. During the study, Twin A sustained SRC, allowing comparison of mRNA expression to the nonconcussed Twin B. Twin A clinically recovered by Day 7. mRNA expression was measured pre-injury and at 6 h and 7 days postinjury using Affymetrix HG-U133 Plus 2.0 microarray. Changes in mRNA expression from pre-injury to each postinjury time point were compared between the twins; differences >1.5-fold were considered important. Kyoto Encyclopedia of Genes and Genomes identified biologic networks associated with important transcripts. Among 38,000 analyzed genes, important changes were identified in 153 genes. The ErbB (epidermal growth factor receptor) signaling pathway was identified as the top transcriptional network from pre-injury to 7 days postinjury. Genes in this pathway with important transcriptional changes included epidermal growth factor (2.41), epiregulin (1.73), neuregulin 1 (1.54) and mechanistic target of rapamycin (1.51). In conclusion, the ErbB signaling pathway was identified as a potential regulator of clinical recovery in a MZ twin pair discordant for SRC. A co-twin study design may be a useful method for identifying important gene pathways associated with concussion recovery.

Type
Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of International Society for Twin Studies

Researchers have long relied on gene expression changes identified in animal models to develop targeted therapy for traumatic brain injury (TBI), but so far these therapies have not proved effective in human trials (Xiong et al., Reference Xiong, Mahmood and Chopp2013). While these preclinical models have provided valuable mechanistic insights, they might not reflect the unique pathophysiology of concussion in humans (Boomsma et al., Reference Boomsma, Busjahn and Peltonen2002). Identifying important and potentially modifiable gene expression changes in humans diagnosed with TBI may be essential for developing effective therapies.

The inability to study expressional changes in injured brain tissue has complicated efforts to identify important gene changes in humans with TBI and in particular with mild TBI (concussion). This obstacle has been largely overcome by the recognition that gene changes in peripheral blood monocyte cells (PBMCs) are an acceptable albeit imperfect reflection of brain transcriptional changes (Sullivan et al., Reference Sullivan, Fan and Perou2006). However, an equally important obstacle, one that has not yet been overcome, is the large degree of genetic ‘noise’ that can obscure identification of a TBI-related transcriptional ‘signal’, when genetically unrelated individuals with and without injury are compared. This noise comes primarily from natural inter-individual variation in gene transcription as well as variation in environmental exposures affecting transcription.

Comparing mRNA expression before and after injury on a subject-specific level is one way to minimize interference from this genetic noise. Athletes are a population particularly well-suited to such studies, as they are readily identifiable at baseline and have a known, measurable risk for concussion. Another way to minimize genetic noise is to compare gene changes in genetically related individuals, such as twins. Monozygotic (MZ) twins have less variation in gene expression (0−1.76%) compared to genetically unrelated individuals (up to 14.13%; Sharma et al., Reference Sharma, Sharma, Horn-Saban, Lancet, Ramachandran and Brahmachari2005). Twin studies examining outcome differences within a twin pair discordant for exposure (co-twin control study) have been used to investigate a variety of potentially harmful brain exposures (e.g. combat stress, trauma, childhood sexual abuse; Brown et al., Reference Brown, Berenz, Aggen, Gardner, Knudsen, Reichborn-Kjennerud, Kendler and Amstadter2014; Kendler et al., Reference Kendler, Bulik, Silberg, Hettema, Myers and Prescott2000; Koenen et al., Reference Koenen, Lyons, Goldberg, Simpson, Williams, Toomey, Elsen, True and Tsuang2003; Mitchell et al., Reference Mitchell, Wolf, Lyons, Goldberg and Magruder2018), but not head injury. The combination of a co-twin control study design with a before-and-after study design would theoretically maximize the identification of gene changes unique to concussion. The serendipitous occurrence of a concussion occurring in one member of a pair of MZ twins participating in larger study of biomarkers in concussed collegiate athletes afforded us the unique opportunity to identify novel transcriptional changes with reduced genetic noise. Because the twins reported here were roommates and teammates, genetic noise may have been further reduced by their exposure to a very similar daily environment.

The objective of the current study was to identify important and potentially novel differences in messenger RNA (mRNA) expression among this single pair of MZ twins discordant for sport-related concussion (SRC).

Materials and Methods

A single pair of right-handed MZ twins was recruited as part of a National Institutes of Health funded study of National Collegiate Athletic Association (NCAA) athletes at the University of Rochester and Rochester Institute of Technology designed to validate blood-based markers of concussion. During the study, one twin (Twin A) sustained a concussion while playing soccer, affording us the rare opportunity to compare mRNA expression to the uninjured twin (Twin B).

The study methods have been described in detail elsewhere (Merchant-Borna et al., Reference Merchant-Borna, Lee, Wang, Bogner, van Griensven, Gill and Bazarian2016). In brief, athletes were eligible for inclusion if they were ≥18 years of age and played an NCAA varsity contact sport. A baseline, preseason (T0) assessment included a venous blood sample, cognitive performance (ImPACT test) and postural stability (Balance Error Scoring System) on all participating athletes. Athletes were followed prospectively for development of SRC as defined by the Third International Conference on Concussion in Sport (McCrory et al., Reference McCrory, Meeuwisse, Johnston, Dvorak, Aubry, Molloy and Cantu2009). After Twin A sustained a SRC, a second blood sample was obtained from both twins within 6 h of the injury (T1 or ‘acutely’ postinjury). White matter integrity was measured with diffusion tensor imaging (DTI) at 2 days and 1-year postinjury. Cognitive performance, balance and symptoms were remeasured 3 days postinjury. A third blood sample, repeat cognitive performance and balance measures were obtained at 7 days postinjury (T2 or ‘subacutely’ postinjury).

The institutional review boards at University of Rochester and Rochester Institute of Technology approved this protocol; written informed consent was obtained from study participants prior to subject participation. Supplementary Material is available on the Cambridge Core website

PBMC and mRNA Isolation

The methods to isolate mRNA from PBMCs have been previously described (Merchant-Borna et al., Reference Merchant-Borna, Lee, Wang, Bogner, van Griensven, Gill and Bazarian2016). In summary, PBMCs were isolated within 1 h of blood collection (Fuss et al., Reference Fuss, Kanof, Smith and Zola2009), and the resulting pellets were suspended in RPMI-10 medium and stored at −190°C until analysis. RNA was isolated from PBMC pellets using TRIzol® Plus RNA Purification Kits (Life Technologies, Grand Island, NY) followed by DNase I-Amplification Grade Kits (Life Technologies). A NanoDrop DN-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE) was used to verify the purity and concentration of the RNA samples. The integrity of the RNA samples was measured on an Agilent bioanalyzer 2100 using the RNA 6000 Nano Kit (Agilent, Sana Clara, CA). Only samples with an RNA integrity number (RIN) greater than 7.0 were used in analysis; samples below this threshold were excluded. Each RNA underwent reverse transcription using a GeneChip® 3’ IVT Expression Kit, converting the RNA strands to biotinylated complementary RNA, which was hybridized to Affymetrix HG-U133 Plus 2.0 microarray (Affymetrix, Sana Clara, CA), stained with streptavidin-phycoerythrin and washed. Raw data were generated by laser scanning imaging.

DTI Methods

Both twins were scanned on a Siemens 3T Tim Trio scanner with a 32-channel head coil using a single-shot spin-echo echoplanar imaging (60 diffusion directions with b = 1200 s/mm2, 10 images with b = 0 s/mm2, TR = 8900 ms, TE = 86 ms, FOV = 256 Å∼ 256 mm2, matrix = 128 Å∼ 128, voxel size = 2 mm3 by 2 mm3 by 2 mm3, 70 axial slices). A double-echo gradient echo field map sequence (echo time difference = 2.46 ms, EPI dwell time = 0.75 ms) was acquired with the same resolution as the DTI sequence and was used to correct for distortion caused by B0 inhomogeneity. Diffusion-weighted images in conjunction with an atlas-based Markov random field representation were used to define 31 WM structures (Bazin et al., Reference Bazin, Ye, Bogovic, Shiee, Reich, Prince and Pham2011). FA values were measured in these 31 regions. For each ROI, the FA difference between Twin A and Twin B was calculated and compared to the 95% CI for the mean FA difference among 65 uninjured female MZ twins (mean age 22.9 + 2.1 yrs) scanned between 2007 and 2012 as part of the Queensland Twin Imaging Study (QTIM; De Zubicaray et al., Reference De Zubicaray, Chiang, McMahon, Shattuck, Toga, Martin, Wright and Thompson2008; Kochunov et al., Reference Kochunov, Jahanshad, Sprooten, Nichols, Mandl, Almasy, Booth, Brouwer, Curran, de Zubicaray, Dimitrova, Duggirala, Fox, Hong, Landman, Lemaitre, Lopez, Martin, McMahon and de Zubicaray2014). ROIs in which FA differences between the twins were either above or below this 95% CI were considered indicative of axonal injury. The QTIM used the Johns Hopkins University white matter parcellation atlas (Mori et al., Reference Mori, Oishi, Jiang, Jiang, Li, Akhter, Hua, Faria, Mahmood, Woods, Toga, Pike, Neto, Evans, Zhang, Huang, Miller, van Zijl and Mazziotta2008) to define 41 ROIs, 15 of which were identical to the ROIs defined in the study twins using the atlas-based Markov method.

Analysis

Partek Genomics Suite Software, version 6.6 (Partek Inc., St. Louis, MO), was used to analyze the microarray data. The robust multiarray average method was used to correct for background signaling once interrogating probes were imported. Additional corrections were applied to the GC content probes. Quantile normalization was used to standardize the probesets, and the data distribution was normalized by transforming gene expression on a log-2 scale.

In order to isolate the effect of concussion on the gene expression in Twin A, we first calculated the change in mRNA expression for each twin across each time point; pre-injury to 6 h postinjury (T1-T0), pre-injury to 7 days postinjury (T2-T0) and 6 h postinjury to 7 days postinjury (T2-T1) (Supplementary Table S1). The difference in the longitudinal change of mRNA expression was then determined between the twins. Taking the difference in the longitudinal change maximizes detection of gene changes unique to brain injury. Important genes were identified based on a fold-change threshold value of greater than 1.5 or less than −1.5 for all three time point comparisons.

Pathway Enrichment Analysis

Gene set enrichment analysis by hypergeometric test was performed in R (R Core Team, 2013), version 1.2.1335, using Biocarta (Nishimura, Reference Nishimura2001), KEGG(Ogata et al., Reference Ogata, Goto, Sato, Fujibuchi, Bono and Kanehisa1999), and Reactome (Vastrik et al., Reference Vastrik, D’Eustachio, Schmidt, Gopinath, Croft, de Bono, Gillespie, Jassal, Lewis, Matthews, Wu, Birney and Stein2007) gene sets from msigdb (Liberzon et al., Reference Liberzon, Birger, Thorvaldsdottir, Ghandi, Mesirov and Tamayo2015) and gmt utilities from qusage (Liberzon et al., Reference Liberzon, Birger, Thorvaldsdottir, Ghandi, Mesirov and Tamayo2015) with a conservative estimate of 20,000 for size of universe. P values were adjusted using the Benjamini–Hochberg method (Benjamini & Hochberg, Reference Benjamini and Hochberg1995) and adjusted p values < .2 output. Cytoscape (Shannon et al., Reference Shannon, Markiel, Ozier, Baliga, Wang, Ramage, Amin, Schwikoswki and Ideker2003) was used to visualize ErbB KEGG pathway. The gene set enrichment analysis was confined to the baseline-to-subacute time interval (T2-T0) in order to understand the genetic changes that may be associated with clinical recovery.

Results

Clinical and White Matter Changes

Twin A was diagnosed with an SRC while participating in collegiate soccer. She was struck in the face by a ball kicked at close range and was witnessed to have a loss of consciousness of several seconds. She was amnestic for the event and reported headache, dizziness, lightheadedness. She was immediately removed from play. Both twins underwent DTI scanning 2 days postinjury. Of 31 ROIs, we were able to compare 15 ROIs to 65 female MZ twin pairs from the QTIM study. Of these 15 ROIs, FA differences between the study twins were outside the regional 95% CI in 5 ROIs (33%). Clinical testing 3 days postinjury revealed a marked worsening in verbal memory, visual motor speed and reaction time — but not postural stability — and increase in concussive symptoms, compared to her personal baseline and compared to testing done simultaneously on Twin B (Figure 1).

Fig. 1. Longitudinal changes in cognitive performance (A-D), postconcussive symptom score (E) and balance scores (F) between concussed (Twin A) and nonconcussed (Twin B) MZ twins. Absolute FA differences between twins (Twin A-Twin B) at 2 day postinjury, compared to 95% CI in FA differences in 15 ROIs among 65 uninjured females MZ twin pairs from the Queensland Twin Imaging Study (G). BCC-body of corpus callosum, CGC-cingulate gyrus, CST-corticospinal tract, GCC-genu of corpus callosum, IFO-inferior fronto-occipital fasciculus, SCC-splenium of corpus callosum, SFO-superior fronto-occipital fasciculus, SLF-superior longitudinal fasciculus, UNC-uncinate fasciculus

By 7 days postinjury, Twin A’s symptoms, verbal memory, visual motor speed and reaction time had returned to baseline levels suggesting clinical recovery. The improvement in Twin A’s visual memory score between 3 and 7 days postinjury is also consistent with concussion recovery. However, why Twin A’s baseline visual memory score was lower than that of Twin B and why Twin B’s score declined between 3 and 7 days in the absence of concussive injury are unclear. Possibilities include sleep deprivation, pain and emotional stress. By Day 7, Twin A was cleared to begin the gradual return to play progression, which she successfully completed on day 13 postinjury. Twin B continued sport participation throughout the duration of this study. By 1-year postinjury, FA differences between the study twins returned to within the 95% CI in 4 of the 5 ROIs that had been outside the 95% CI at 2 days (not shown).

Messenger RNA Expression Changes

54,675 probesets representing 38,000 genes were compared for important differences in longitudinal mRNA expression between the twins. Changes in 158 probes representing 153 unique genes had either a greater than 1.5 or less than −1.5 difference in the longitudinal change between the twins and were thus considered important (Figure 2). The majority (144) of these changed genes were protein-coding; however, 2 antisense RNA transcripts and 7 noncoding RNA transcripts also were changed. Considering all time points, more probesets were upregulated (111 probesets) than downregulated (68 probesets).

Fig. 2. Heat map of the difference in longitudinal change in mRNA expression between single MZ twin pair discordant for concussion across the three time intervals (X-axis): baseline to acute (T0 to T1), baseline to subacute (T0 to T2) and acute to subacute (T1 to T2). The Y-axis represents individual genes, color-coded by upregulation (red) or downregulation (blue)

Baseline to acute

From pre-injury baseline (T0) to the acute postinjury time point (T1), 79 probesets had important changes between the twins, representing 77 unique genes; 38 probesets were upregulated, while 38 probesets were downregulated (Supplementary Tables S2a and S2b). Of these 79 probesets, 64 were unique to this time interval, while 15 probesets had important changes in both baseline-to-acute and baseline-to-subacute time intervals. Fibrillin 2 (FBN2; 2.06-fold), Ras and Rab interactor 2 (RIN2; 2.04-fold), and Ikaros family zinc protein finger 1 (IKZF1; 2.01-fold) were found to have the largest differential change of the upregulated genes. Pregnancy-specific beta 1-glycoprotein 6 (PSG6; −2.0-fold), immunoglobulin lambda constant 1 (IGLC1; −2.0-fold) and proyl 3-hydroxylase (P3H2; −2.25-fold) had the largest downregulation. Another notable differentially expressed gene during this time interval was calcium channel, voltage-dependent, L type and alpha 1F subunit (CACNA1F, −1.90-fold).

Baseline to subacute

From pre-injury baseline (T0) to the subacute time point (T2), 71 probesets had important changes in expression, representing 68 unique genes. Of these probesets, 50 were upregulated and 21 were downregulated (Supplementary Tables S3a and S3b). Fifty probesets were unique to this time interval, 15 overlapped with the baseline-to-acute time interval, and 6 overlapped with the acute to subacute time interval. Mitochondrial ribosomal protein L19 (MRPL1; 2.13-fold) and G protein-coupled receptor 155 (GPR155; 2.12-fold) had the largest differential change among the upregulated genes. Important changes in the downregulated genes included transmembrane protein with EGF-like and two follistatin-like domains 2 (TMEFF2; −2.80-fold), epidermal growth factor receptor (EGFR; −2.41-fold; Probe ID:1565483_at), immunoglobulin kappa constant (IGKC; −2.19-fold; Probe ID: 214768_x_at) and prolyl 3-hydroxylase 2 (P3H2; −2.01-fold). Other notable differentially expressed genes during this time interval included neuregulin 1 (NRG1, 1.54-fold) and glutamate-cysteine ligase, modifier subunit (GCLM, 1.75-fold).

Acute to subacute

From the acute (T1) to subacute time points (T2), 29 probesets had important changes in expression, representing 29 genes. Of these, 23 were upregulated and 6 were downregulated (Supplementary Tables S4a and S4b). 23 probesets were unique to the acute to subacute time interval, while 6 probesets were also importantly changed from the baseline-to-acute time interval. Glutamine and serine-rich 1 (QSER1; 2.19-fold) and poly(rC) binding protein 2 (PCBP2; 1.92-fold) had the greatest differential upregulated change. Mitotic spindle organizing protein 1 (MZT1; −1.76-fold) and pleckstrin homology-like domain, family B, member 3 (PHLDB3; −1.55-fold) had the largest differential downregulated change during this time period. Another notable differentially expressed gene during this time interval was N-methyl D-aspartate 2D (GRIN2D, 1.52-fold).

Pathway enrichment analysis: Baseline to subacute

The top pathway of gene changes from baseline pre-injury (T0) to subacute postinjury (T2) was related to the ErbB signaling pathway (p-adj = .010826; Figure 3). Individual genes in this pathway that were importantly changed (<−1.5 or >1.5) in the MZ twins included EGFR, epiregulin (EREG), neuregulin 1 (NRG1), and mechanistic target of rapamycin (MTOR). EREG had a total increased expression from baseline pre-injury to subacute postinjury of 1.73-fold. NRG1 was upregulated 1.54-fold, while MTOR was upregulated 1.51-fold. Conversely, EGFR(Probe ID: 1565483_at) was downregulated, with a total -2.41-fold change from baseline pre-injury to subacute postinjury.

Fig. 3. ErbB signaling pathway from KEGG.(Ogata et al., Reference Ogata, Goto, Sato, Fujibuchi, Bono and Kanehisa1999) Genes with important differential change between the twins discordant for concussion from baseline to subacute (T0 to T2) are shown in either blue (downregulation) or red (upregulation). Genes of the ErbB pathway that did not have important upregulation or downregulation are indicated in green. Receptors of the ErbB pathway are either homodimerized or heterodimerized, demonstrated in the figure by the blocks containing two gene symbols (e.g. ErbB4/ErbB4). Important upregulated genes (red) included epiregulin (EREG), neuregulin 1 (NRG1) and mechanistic target of rapamycin (MTOR). Epidermal growth factor (EGFR) was downregulated (blue). The +P symbol indicates a gene in which the encode protein is responsible for a phosphorylation event in the ErbB pathway

Discussion

In the current study, we identified the ErbB signaling pathway as a potential regulator of recovery in single pair of MZ twins discordant for concussion. The chance occurrence of concussion in one member of a pair of MZ twins provided the rare and unique opportunity to identify novel concussion-related transcriptional changes in the setting of reduced genetic noise. An important study strength was the acquisition of PBMCs from each twin before injury, due to their participation in a larger NIH-funded study. This allowed us to further isolate the effects of concussion by comparing pre- and postinjury transcriptional changes in each twin. Another key strength was the use of structural neuroimaging to verify the diagnosis of concussion. Although this diagnosis was based on Twin A’s symptoms at the time of injury, the finding of abnormal FA difference between the twins in five white matter brain regions acutely postinjury increased confidence in the validity of a concussion diagnosis and suggests her injury resulted in reduced white matter integrity. The resolution of the FA differences in four of these five regions 1 year after clinical recovery further supports that these ROIs were injured acutely.

Although we observed important changes in single transcripts at both the acute and subacute time points, the pathway enrichment analysis of changes at the subacute time point, when Twin B was back to baseline clinically, is likely a more robust indicator of biologic pathways involved in concussion recovery. Pathway enrichment analysis has several advantages over analyzing single genes. Statistical power is improved by reducing the dimensionality of large gene lists into a smaller number of systems or pathways (Reimand et al., Reference Reimand, Isserlin, Voisin, Kucera, Tannus-Lopes, Rostamianfar, Wadi, Meyer, Wong, Xu, Merico and Bader2019). Moreover, pathway analysis simplifies data interpretation by constructing familiar concepts (e.g. biological processes) from the gene set, facilitating identification of putative causal mechanisms and drug targets (Reimand et al., Reference Reimand, Isserlin, Voisin, Kucera, Tannus-Lopes, Rostamianfar, Wadi, Meyer, Wong, Xu, Merico and Bader2019). Previous concussion research has taken advantage of pathway analysis to better interpret and understand genetic data (Johnson et al., Reference Johnson, Loeffert, Stokes, Olympia, Bramley and Hicks2018; LaRocca et al., Reference LaRocca, Barns, Hicks, Brindle, Williams, Uhlig, Johnson, Neville and Middleton2019; Merchant-Borna et al., Reference Merchant-Borna, Lee, Wang, Bogner, van Griensven, Gill and Bazarian2016).

Our pathway analysis of mRNA expression between pre-injury baseline and 7 days postinjury suggested the ErbB signaling pathway to be a potential regulator of the biologic processes involved in recovery following SRC. Specific genes in this pathway with important longitudinal changes in the twins included EREG, NRG1, MTOR and EGFR. Although insufficient ErbB signaling is known to be associated with neurodegenerative diseases, such as Alzheimer’s disease and multiple sclerosis, its role in TBI recovery has been less frequently reported (Bublil & Yarden, Reference Bublil and Yarden2007). In an animal model of TBI, Erlich et al. (Reference Erlich, Shohami and Pinkas-Kramarski2000) demonstrated significant upregulation of ErbB4 receptor (which binds NRG1) in neuronal and glial cells at the site of injury. LaRocca et al. (Reference LaRocca, Barns, Hicks, Brindle, Williams, Uhlig, Johnson, Neville and Middleton2019) reported increased expression of salivary microRNAs targeting transcripts of the ErbB pathway among 50 amateur mixed martial artists exposed to subconcussive head impacts during a fight (Atif & Hicks, Reference Atif and Hicks2019).

ErbB signaling consists of four homodimerized, or heterodimerized, receptors that regulate a variety of pathways including PI3K/AKT/mTOR, p38/MAPK, calcium signaling and others (Yarden & Shilo, Reference Yarden and Shilo2007; Kataria et al., Reference Kataria, Alizadeh and Karimi-Abdolrezaee2019). Previous TBI studies have targeted the PI3K/AKT/mTOR and p38/MAPK pathways post-TBI in order to improve functional outcome (Bachstetter et al., Reference Bachstetter, Rowe, Kaneko, Goulding, Lifshitz and Van Eldik2013; Erlich et al., Reference Erlich, Alexandrovich, Shohami and Pinkas-Kramarski2007). Postmortem brain tissues of CTE patients have demonstrated a number of genes related to the MAPK and calcium signaling pathway to be significantly downregulated (Seo et al., Reference Seo, Lee, Shin, Hwang, Cho, Yoo, Kim, Lim, Kim, Hwang, Kim, Kim, Hyeon, Yun, Kim, Kim, Alvarez, Stein, Lee and McKee2017). Yet, studies have not extensively investigated the role of the ErbB signaling in regulating these downstream pathways postconcussion.

We identified NRG1, a growth factor that binds ErbB receptors, to be upregulated from pre-injury baseline to 7 days postinjury. The exact role this gene plays in concussion recovery cannot be determined from our results, but evidence from other disease processes (spinal cord injury, brain ischemia, Alzheimer’s disease, Parkinson’s disease) suggests that NRG1 controls several critical functions that could potentially promote concussion recovery including myelination, blood–brain barrier (BBB) function, synaptic function and immune response.(Kataria et al., Reference Kataria, Alizadeh and Karimi-Abdolrezaee2019)

NRG1 has been shown to protect microvasculature and reduce permeability of the BBB via signaling in brain endothelial cells (Lok et al., Reference Lok, Sardi, Guo, Besancon, Ha, Rosell, Kim, Corfas and Lo2009; Lok et al., Reference Lok, Zhao, Leung, Seo, Navaratna, Wang, Whalen and Lo2012), which could mitigate disruption of the BBB following TBI (Blyth et al., Reference Blyth, Farhavar, Gee, Hawthorn, He, Nayak, Stöcklein and Bazarian2009). In the postconcussion period, white matter tracts are susceptible to injury and demyelination, which is characterized by loss of myelin sheath and oligodendrocyte cell death (Mierzwa et al., Reference Mierzwa, Marion, Sullivan, McDaniel and Armstrong2015; Shi et al., Reference Shi, Hu, Leak, Shi, An, Suenaga, Chen and Gao2015). NRG1 plays an important role in the development, maturation and survival of oligodendrocytes and seems to affect the quality of myelination in the CNS (Shahriary et al., Reference Shahriary, Kataria and Karimi-Abdolrezaee2019). Additionally, NRG1 serves a neuroprotective role through anti-inflammatory actions, including the regulation of microglia, astrocyte and immune cell response after in vitro CNS injury (Alizadeh et al., Reference Alizadeh, Dyck, Kataria, Shahriary, Nguyen, Santhosh and Karimi-Abdolrezaee2017; Shahriary et al., Reference Shahriary, Kataria and Karimi-Abdolrezaee2019). These studies provide evidence for a possible role of NRG1 in the brain’s recovery process following concussion through the ErbB signaling pathway.

Limitations

An obvious limitation of the study is sample size, with only one twin pair evaluated, which increases the uncertainty that the important genetic changes we identified are unique to concussion. Additionally, as the MZ twins were athletes, the generalizability of results must be carefully considered, as physical exertion itself has been found to alter gene expression (Connolly et al., Reference Connolly, Caiozzo, Zaldivar, Nemet, Larson, Hung, Heck, Hatfield and Cooper2004). While the concussed twin was prohibited from returning to play until recovered, the control twin continued to engage in physical activity. Therefore, some of the important gene changes we identified may reflect not only processes related to brain injury but also to physical exertion. Furthermore, epigenetic divergence of MZ twins could contribute to the differences in gene expression observed (Fraga et al., Reference Fraga, Ballestar, Paz, Ropero, Setien, Ballestar, Heine-Suñer, Cigudosa, Urioste, Benitez, Boix-Chornet, Sanchez-Aguilera, Ling, Carlsson, Poulsen, Vaag, Stephan, Spector and Esteller2005), with additional genetic variability in female twin pairs (such as the pair in the present study) caused by random x-inactivation and/or reproductive hormone cycling.

Conclusions

By evaluating a set of MZ twins discordant for concussion and reduced white matter integrity, we identified the ErbB signaling pathway as the top differentially expressed gene network from pre-injury to 7 days postinjury. Since the concussed twin appeared clinically recovered by 7 days, we believe that ErbB signaling may be important for concussion recovery. Understanding genetic regulation of concussion recovery can help identify biologic processes that may be modifiable and suitable targets for future therapies. Because MZ twins have less inter-individual difference in gene expression than the general population, focusing on them as study subjects may facilitate detection of gene networks specific to concussion. A co-twin study design may thus be a new avenue for identifying modifiable targets with an improved chance of leading to effective concussion treatments.

Supplementary Material

To view supplementary material for this article, please visit https://doi.org/10.1017/thg.2022.15.

Data

The full gene data set used in this study has been publicly available on Gene Expression Omnibus (GEO) repository since December 31, 2021, under the following code/name: GSE176388 - Expression data from MZ Twins discordant for concussion

Acknowledgments

The authors wish to acknowledge the head athletic trainer at the Rochester Institute of Technology, Ben Emke, for his assistance and support. The authors gratefully acknowledge the assistance of RIT’s athletic training staff, as well as Kirsten Ross, Stephanie Amalfe and Nikita Bourque for their coordination efforts.

Author contributions

JB, JG, KMB, ES were responsible for the design and conceptualization of the study. JB, JG, KMB were responsible for the acquisition of data. TS, RP, KMB, ES, JG, JB contributed to the analysis of the data. TS, JB, ACW contributed to the interpretation of the data. JB drafted the manuscript. All authors contributed to the revision of the manuscript for intellectual content.

Financial support

This work was supported by NICHD K24HD064754 (JJB), University of Rochester CTSA award number TL1 TR000096 from the National Center for Advancing Translational Sciences (ELS) and from the NIH Intramural Research Program (JG).

Conflict of interest

J. J. Bazarian: reports research support from BrainScope and scientific advisory board for Abbott Point of Care Diagnostics. There is no conflict of interest for all other authors.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional guides on the care and use of laboratory animals.

References

Alizadeh, A., Dyck, S. M., Kataria, H., Shahriary, G. M., Nguyen, D. H., Santhosh, K. T., & Karimi-Abdolrezaee, S. (2017). Neuregulin-1 positively modulates glial response and improves neurological recovery following traumatic spinal cord injury. Glia, 65, 11521175.CrossRefGoogle ScholarPubMed
Atif, H., & Hicks, S. D. (2019). A review of microRNA biomarkers in traumatic brain injury. Journal of Experimental Neuroscience, 13, 1179069519832286.CrossRefGoogle ScholarPubMed
Bachstetter, A. D., Rowe, R. K., Kaneko, M., Goulding, D., Lifshitz, J., & Van Eldik, L. J. (2013). The p38alpha MAPK regulates microglial responsiveness to diffuse traumatic brain injury. Journal of Neuroscience, 33, 61436153.CrossRefGoogle ScholarPubMed
Bazin, P. L., Ye, C., Bogovic, J. A., Shiee, N., Reich, D. S., Prince, J. L., & Pham, D. L. (2011). Direct segmentation of the major white matter tracts in diffusion tensor images. NeuroImage, 58, 458468.CrossRefGoogle ScholarPubMed
Benjamini, Y., & Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological), 57, 289300.CrossRefGoogle Scholar
Blyth, B. J., Farhavar, A., Gee, C., Hawthorn, B., He, H., Nayak, A., Stöcklein, V., & Bazarian, J. J. (2009). Validation of serum markers for blood-brain barrier disruption in traumatic brain injury. Journal of Neurotrauma, 26, 14971507.CrossRefGoogle ScholarPubMed
Boomsma, D., Busjahn, A., & Peltonen, L. (2002). Classical twin studies and beyond. Nature Reviews Genetics, 3, 872882.CrossRefGoogle ScholarPubMed
Brown, R. C., Berenz, E. C., Aggen, S. H., Gardner, C. O., Knudsen, G. P., Reichborn-Kjennerud, T., Kendler, K. S., & Amstadter, A. B. (2014). Trauma exposure and axis i psychopathology: A co-twin control analysis in Norwegian young adults. Psychological Trauma: Theory, Research, Practice, and Policy, 6, 652660.CrossRefGoogle ScholarPubMed
Bublil, E. M., & Yarden, Y. (2007). The EGF receptor family: Spearheading a merger of signaling and therapeutics. Current Opinion in Cell Biology, 19, 124134.CrossRefGoogle ScholarPubMed
Connolly, P. H., Caiozzo, V. J., Zaldivar, F., Nemet, D., Larson, J., Hung, S. P., Heck, J. D., Hatfield, G. W., & Cooper, D. M. (2004). Effects of exercise on gene expression in human peripheral blood mononuclear cells. Journal of Applied Physiology (1985), 97, 14611469.CrossRefGoogle ScholarPubMed
De Zubicaray, G. I., Chiang, M.-C., McMahon, K. L., Shattuck, D. W., Toga, A. W., Martin, N. G., Wright, M. J., & Thompson, P. M. (2008). Meeting the challenges of neuroimaging genetics. Brain imaging and Behavior, 2, 258263.CrossRefGoogle ScholarPubMed
Erlich, S., Alexandrovich, A., Shohami, E., & Pinkas-Kramarski, R. (2007). Rapamycin is a neuroprotective treatment for traumatic brain injury. Neurobiology of Disease, 26, 8693.CrossRefGoogle ScholarPubMed
Erlich, S., Shohami, E., & Pinkas-Kramarski, R. (2000). Closed head injury induces up-regulation of ErbB-4 receptor at the site of injury. Molecular and Cellular Neuroscience, 16, 597608.CrossRefGoogle ScholarPubMed
Fraga, M. F., Ballestar, E., Paz, M. F., Ropero, S., Setien, F., Ballestar, M. L., Heine-Suñer, D., Cigudosa, J. C., Urioste, M., Benitez, J., Boix-Chornet, M., Sanchez-Aguilera, A., Ling, C., Carlsson, E., Poulsen, P., Vaag, A., Stephan, Z., Spector, T. D., … Esteller, M. (2005). Epigenetic differences arise during the lifetime of monozygotic twins. Proceedings of the National Academy of Sciences of the United States of America, 102, 1060410609.CrossRefGoogle ScholarPubMed
Fuss, I. J., Kanof, M. E., Smith, P. D., & Zola, H. (2009). Isolation of whole mononuclear cells from peripheral blood and cord blood. Current Protocols in Immunology, 85, 7-1.CrossRefGoogle Scholar
Johnson, J. J., Loeffert, A. C., Stokes, J., Olympia, R. P., Bramley, H., & Hicks, S. D. (2018). Association of salivary microRNA changes with prolonged concussion symptoms. JAMA Pediatrics, 172, 6573.CrossRefGoogle ScholarPubMed
Kataria, H., Alizadeh, A., & Karimi-Abdolrezaee, S. (2019). Neuregulin-1/ErbB network: An emerging modulator of nervous system injury and repair. Progress in Neurobiology, 180, 101643.CrossRefGoogle ScholarPubMed
Kendler, K. S., Bulik, C. M., Silberg, J., Hettema, J. M., Myers, J., & Prescott, C. A. (2000). Childhood sexual abuse and adult psychiatric and substance use disorders in women: An epidemiological and cotwin control analysis. Archives of General Psychiatry, 57, 953959.CrossRefGoogle ScholarPubMed
Kochunov, P., Jahanshad, N., Sprooten, E., Nichols, T. E., Mandl, R. C., Almasy, L., Booth, T., Brouwer, R. M., Curran, J. E., de Zubicaray, G. I., Dimitrova, R., Duggirala, R., Fox, P. T., Hong, L. E., Landman, B. A., Lemaitre, H., Lopez, L. M., Martin, N. G., McMahon, K. L., … de Zubicaray, G. I. (2014). Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: Comparing meta and megaanalytical approaches for data pooling. NeuroImage, 95, 136150.CrossRefGoogle ScholarPubMed
Koenen, K. C., Lyons, M. J., Goldberg, J., Simpson, J., Williams, W. M., Toomey, R., Elsen, S. A., True, W., & Tsuang, M. T. (2003). Co-twin control study of relationships among combat exposure, combat-related PTSD, and other mental disorders. Journal of Traumatic Stress, 16, 433438.CrossRefGoogle ScholarPubMed
LaRocca, D., Barns, S., Hicks, S. D., Brindle, A., Williams, J., Uhlig, R., Johnson, P., Neville, C., & Middleton, F. A. (2019). Comparison of serum and saliva miRNAs for identification and characterization of mTBI in adult mixed martial arts fighters. PLoS One, 14, e0207785e0207785.CrossRefGoogle ScholarPubMed
Liberzon, A., Birger, C., Thorvaldsdottir, H., Ghandi, M., Mesirov, J. P., & Tamayo, P. (2015). The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Systems, 1, 417425.CrossRefGoogle ScholarPubMed
Lok, J., Sardi, S. P., Guo, S., Besancon, E., Ha, D. M., Rosell, A., Kim, W. J., Corfas, G., & Lo, E. H. (2009). Neuregulin-1 signaling in brain endothelial cells. Journal of Cerebral Blood Flow & Metabolism, 29, 3943.CrossRefGoogle ScholarPubMed
Lok, J., Zhao, S., Leung, W., Seo, J. H., Navaratna, D., Wang, X., Whalen, M. J., & Lo, E. H. (2012). Neuregulin-1 effects on endothelial and blood-brain-barrier permeability after experimental injury. Translational Stroke Research, 3, S119124.CrossRefGoogle ScholarPubMed
McCrory, P., Meeuwisse, W., Johnston, K., Dvorak, J., Aubry, M., Molloy, M., & Cantu, R. (2009). Consensus statement on concussion in sport: The 3rd International Conference on Concussion in Sport held in Zurich, November 2008. Journal of Athletic Training, 44, 434448.CrossRefGoogle Scholar
Merchant-Borna, K., Lee, H., Wang, D., Bogner, V., van Griensven, M., Gill, J., & Bazarian, J. J. (2016). Genome-wide changes in peripheral gene expression following sports-related concussion. Journal of Neurotrauma, 33, 15761585.CrossRefGoogle ScholarPubMed
Mierzwa, A. J., Marion, C. M., Sullivan, G. M., McDaniel, D. P., & Armstrong, R. C. (2015). Components of myelin damage and repair in the progression of white matter pathology after mild traumatic brain injury. Journal of Neuropathology & Experimental Neurology, 74, 218232.CrossRefGoogle ScholarPubMed
Mitchell, K. S., Wolf, E. J., Lyons, M. J., Goldberg, J., & Magruder, K. M. (2018). A co-twin control study of the association between combat exposure, PTSD and obesity in male veterans. Psychological Medicine, 48, 24492452.CrossRefGoogle ScholarPubMed
Mori, S., Oishi, K., Jiang, H., Jiang, L., Li, X., Akhter, K., Hua, K., Faria, A. V., Mahmood, A., Woods, R., Toga, A. W., Pike, G. B., Neto, P. R., Evans, A., Zhang, J., Huang, H., Miller, M. I., van Zijl, P., & Mazziotta, J. (2008). Stereotaxic white matter atlas based on diffusion tensor imaging in an ICBM template. NeuroImage, 40, 570582.CrossRefGoogle Scholar
Nishimura, D. (2001). BioCarta. Biotech Software & Internet Report, 2, 117120.CrossRefGoogle Scholar
Ogata, H., Goto, S., Sato, K., Fujibuchi, W., Bono, H., & Kanehisa, M. (1999). KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Research, 27, 2934.CrossRefGoogle ScholarPubMed
R Core Team. (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing.Google Scholar
Reimand, J., Isserlin, R., Voisin, V., Kucera, M., Tannus-Lopes, C., Rostamianfar, A., Wadi, L., Meyer, M., Wong, J., Xu, C., Merico, D., & Bader, G. D. (2019). Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and EnrichmentMap. Nature Protocols, 14, 482517.CrossRefGoogle ScholarPubMed
Seo, J.-S., Lee, S., Shin, J.-Y., Hwang, Y. J., Cho, H., Yoo, S.-K., Kim, Y., Lim, S., Kim, Y. K., Hwang, E. M., Kim, S. H., Kim, C. H., Hyeon, S. J., Yun, J. Y., Kim, J., Kim, Y., Alvarez, V. E., Stein, T. D., Lee, J., … McKee, A. C. (2017). Transcriptome analyses of chronic traumatic encephalopathy show alterations in protein phosphatase expression associated with tauopathy. Experimental & Molecular Medicine, 49, e333e333.CrossRefGoogle ScholarPubMed
Shahriary, G. M., Kataria, H., & Karimi-Abdolrezaee, S. (2019). Neuregulin-1 fosters supportive interactions between microglia and neural stem/progenitor cells. Stem Cells International, 2019, 8397158.CrossRefGoogle ScholarPubMed
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., Amin, N., Schwikoswki, B., & Ideker, T. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Research, 13, 24982504.CrossRefGoogle ScholarPubMed
Sharma, A., Sharma, V. K., Horn-Saban, S., Lancet, D., Ramachandran, S., & Brahmachari, S. K. (2005). Assessing natural variations in gene expression in humans by comparing with monozygotic twins using microarrays. Physiological Genomics, 21, 117123.CrossRefGoogle ScholarPubMed
Shi, H., Hu, X., Leak, R. K., Shi, Y., An, C., Suenaga, J., Chen, J., & Gao, Y. (2015). Demyelination as a rational therapeutic target for ischemic or traumatic brain injury. Experimental Neurology, 272, 1725.CrossRefGoogle ScholarPubMed
Sullivan, P. F., Fan, C., & Perou, C. M. (2006). Evaluating the comparability of gene expression in blood and brain. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 141B, 261268.CrossRefGoogle ScholarPubMed
Vastrik, I., D’Eustachio, P., Schmidt, E., Gopinath, G., Croft, D., de Bono, B., Gillespie, M., Jassal, B., Lewis, S., Matthews, L., Wu, G., Birney, E., & Stein, L. (2007). Reactome: A knowledge base of biologic pathways and processes. Genome Biology, 8, R39.CrossRefGoogle ScholarPubMed
Xiong, Y., Mahmood, A., & Chopp, M. (2013). Animal models of traumatic brain injury. Nature Reviews Neuroscience, 14, 128142.CrossRefGoogle ScholarPubMed
Yarden, Y., & Shilo, B. Z. (2007). SnapShot: EGFR signaling pathway. Cell, 131, 1018.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Longitudinal changes in cognitive performance (A-D), postconcussive symptom score (E) and balance scores (F) between concussed (Twin A) and nonconcussed (Twin B) MZ twins. Absolute FA differences between twins (Twin A-Twin B) at 2 day postinjury, compared to 95% CI in FA differences in 15 ROIs among 65 uninjured females MZ twin pairs from the Queensland Twin Imaging Study (G). BCC-body of corpus callosum, CGC-cingulate gyrus, CST-corticospinal tract, GCC-genu of corpus callosum, IFO-inferior fronto-occipital fasciculus, SCC-splenium of corpus callosum, SFO-superior fronto-occipital fasciculus, SLF-superior longitudinal fasciculus, UNC-uncinate fasciculus

Figure 1

Fig. 2. Heat map of the difference in longitudinal change in mRNA expression between single MZ twin pair discordant for concussion across the three time intervals (X-axis): baseline to acute (T0 to T1), baseline to subacute (T0 to T2) and acute to subacute (T1 to T2). The Y-axis represents individual genes, color-coded by upregulation (red) or downregulation (blue)

Figure 2

Fig. 3. ErbB signaling pathway from KEGG.(Ogata et al., 1999) Genes with important differential change between the twins discordant for concussion from baseline to subacute (T0 to T2) are shown in either blue (downregulation) or red (upregulation). Genes of the ErbB pathway that did not have important upregulation or downregulation are indicated in green. Receptors of the ErbB pathway are either homodimerized or heterodimerized, demonstrated in the figure by the blocks containing two gene symbols (e.g. ErbB4/ErbB4). Important upregulated genes (red) included epiregulin (EREG), neuregulin 1 (NRG1) and mechanistic target of rapamycin (MTOR). Epidermal growth factor (EGFR) was downregulated (blue). The +P symbol indicates a gene in which the encode protein is responsible for a phosphorylation event in the ErbB pathway

Supplementary material: PDF

Stahl et al. supplementary material

Tables S1-S4

Download Stahl et al. supplementary material(PDF)
PDF 573.5 KB