- Research
- Open access
- Published:
Mapping dynamic molecular changes in hippocampal subregions after traumatic brain injury through spatial proteomics
Clinical Proteomics volume 21, Article number: 32 (2024)
Abstract
Background
Traumatic brain injury (TBI) often results in diverse molecular responses, challenging traditional proteomic studies that measure average changes at tissue levels and fail to capture the complexity and heterogeneity of the affected tissues. Spatial proteomics offers a solution by providing insights into sub-region-specific alterations within tissues. This study focuses on the hippocampal sub-regions, analyzing proteomic expression profiles in mice at the acute (1Â day) and subacute (7Â days) phases of post-TBI to understand subregion-specific vulnerabilities and long-term consequences.
Methods
Three mice brains were collected from each group, including Sham, 1-day post-TBI and 7-day post-TBI. Hippocampal subregions were extracted using Laser Microdissection (LMD) and subsequently analyzed by label-free quantitative proteomics.
Results
The spatial analysis reveals region-specific protein abundance changes, highlighting the elevation of FN1, LGALS3BP, HP, and MUG-1 in the stratum moleculare (SM), suggesting potential immune cell enrichment post-TBI. Notably, established markers of chronic traumatic encephalopathy, IGHM and B2M, exhibit specific upregulation in the dentate gyrus bottom (DG2) independent of direct mechanical injury. Metabolic pathway analysis identifies disturbances in glucose and lipid metabolism, coupled with activated cholesterol synthesis pathways enriched in SM at 7-Day post-TBI and subsequently in deeper DG1 and DG2 suggesting a role in neurogenesis and the onset of recovery. Coordinated activation of neuroglia and microtubule dynamics in DG2 suggest recovery mechanisms in less affected regions. Cluster analysis revealed spatial variations post-TBI, indicative of dysregulated neuronal plasticity and neurogenesis and further predisposition to neurological disorders. TBI-induced protein upregulation (MUG-1, PZP, GFAP, TJP, STAT-1, and CD44) across hippocampal sub-regions indicates shared molecular responses and links to neurological disorders. Spatial variations were demonstrated by proteins dysregulated in both or either of the time-points exclusively in each subregion (ELAVL2, CLIC1 in PL, CD44 and MUG-1 in SM, and SHOC2, LGALS3 in DG).
Conclusions
Utilizing advanced spatial proteomics techniques, the study unveils the dynamic molecular responses in distinct hippocampal subregions post-TBI. It uncovers region-specific vulnerabilities and dysregulated neuronal processes, and potential recovery-related pathways that contribute to our understanding of TBI’s neurological consequences and provides valuable insights for biomarker discovery and therapeutic targets.
Introduction
Traumatic Brain Injury (TBI), a significant global health challenge, is a leading cause of death and disability worldwide, affecting an estimated 60 million individuals globally each year [1]. Despite extensive research spanning decades and numerous therapeutic strategies focusing on neuroprotection, neurovascular regeneration, and neuro-restoration, discovering a truly effective therapy for TBI remains a formidable challenge.
TBI, primarily caused by external forces impacting the head, results in extensive and lasting neuropathological damage through a complex cascade of secondary brain injury, leading to progressive neurological impairment [2]. The intricate nature of TBI is derived from an incomplete understanding of the complex molecular processes within the injured microenvironment. This knowledge gap includes changes at the molecular level in specific micro-regions of the brain after injury. Research has increasingly focused on the molecular changes specific to certain brain regions, particularly in relation to neurogenesis, synaptogenesis, angiogenesis, and axonal remodeling [3, 4]. TBI significantly influences synaptic plasticity and neurogenesis, two of the pivotal processes that govern brain function and recovery [5]. The hippocampus, a key area for memory, is notably impacted by TBI, particularly in maintaining synaptic plasticity and neurogenesis, often leading to impaired memory functions [6, 7]. The impact on the hippocampus, its subregions, and connected neural pathways can selectively influence specific sensory-modalities, consequently affecting certain memory-related behaviors [1]. These behaviors and altered memory plasticity are commonly reported as future clinical manifestations following TBI, such as parkinsonism [8, 9]. For example, immunohistochemistry results showed that VEGF and tau protein, which account for changes in synaptic plasticity were elevated in the dorsal hippocampus in TBI [10]. In addition, changes in expression levels of MAGI3 observed in hippocampal proteomics post TBI are linked to a wide series of molecular and biological system interactors. These changes may be related to short-term neuroplasticity phenomena or long-term neurodegenerative consequences such as parkinsonism [1]. As a site active in neurogenesis, immature new neurons in the hippocampal dentate gyrus (DG) subregion of mice have been shown to be vulnerable to injury, exhibiting significant degeneration as early as 2 days post-injury [11] and are likely to undergo apoptosis [12]. However, the pathophysiological mechanisms contributing to secondary brain injury in moderate/severe TBI remain largely elusive, primarily due to the tissue and cellular heterogeneity in affected brain regions associated with cognitive impairments and subsequent neurological disorders.
Traditional proteomics and immunohistochemistry have significantly advanced our understanding of TBI pathogenesis [13,14,15]. However, these methods often lack the spatial resolution necessary to dissect the complex changes that occur at the microenvironment level following TBI. To address this limitation, we have employed mass-spectrometry based spatial proteomics to scrutinize the molecular landscape, mapping the alterations in protein expression within distinct hippocampal subregions during the acute (1 day) and subacute (7 days) phases [16] which are particularly susceptible to long-term neurological deficits following TBI [17]. By unveiling spatially distinct protein expression patterns, our approach aims to illuminate TBI's impact on critical processes associated with the development of neurological disorders, particularly within the various hippocampal subregions. This investigation provides valuable insights that may help pave the way for the development of targeted and effective therapeutic interventions for TBI and the subsequent emergence of neurodegenerative disorders.
Methods
Animals
All animal experiments were carried out following National Institutes of Health (NIH) guidelines and all protocols were approved by the Institution Animal Care and Use Committee of the Tulane University School of Medicine (Protocol #2000). C57BL/6J male mice were purchased from Jackson Lab (Stock No: 000664). All experiments were conducted with biological replicates (n = 3). Our study included 9 animals including the three animals in Sham group and 6 experimental injury animals (three in 1-Day post-TBI and three in 7-Day post-TBI).
Controlled cortical impact model
Controlled Cortical Impact (CCI) models were induced in male C57BL/6J mice (10 weeks old, 25–27 g) using a pneumatically controlled cortical impactor device (Precision Systems and Instrumentation LLC, TBI-0310), following a previously described protocol with slight modification [18]. Any mice deviating from precise parameters such as having unusually small or large brain volumes were excluded from the study. Briefly, mice were anesthetized in a chamber containing 2.5% isoflurane (Anaquest, Memphis, TN, USA) in 70% N2O and 30% O2, fixed in a stereotaxic apparatus with a gas anesthesia mask using 2% isoflurane delivered through a Fluotec 3 vaporizer (Colonial Medical, Amherst, NH, USA). The scalp was opened to expose the skull and a portable trephine drill (Fine Science Tools, Foster City, CA, USA) was used to create a 5-mm craniotomy on the left cerebral hemisphere equidistant between bregma and lambda. Subsequently, a moderate-severe injury was induced in the ipsilateral cortex and hippocampus, accompanied by pericontusional axonal injury using a 3 mm flat-tip impounder (4.6 m/s impact velocity, 0.65 mm impact depth, and 500 ms duration). Following CCI, the incision was closed with interrupted 4–0 silk sutures, and mice were allowed to recover in their cages. The entire surgical procedure lasted approximately 8 min. Sham mice underwent the same surgical steps, including craniotomy, but without cortical impact, followed by suturing. Mice were anesthetized with 3% isoflurane for 1–2 min and transcardially perfused with 50 mL of ice-cold PBS followed by 30 mL of 4% paraformaldehyde (PFA) in PBS. The brains were rapidly removed and immersed in 4% PFA at 4 °C for 24 h.
Sectioning and slide preparation
Brains were sectioned using a vibrating microtome (Leica VT1000 E, USA). The dorsal brain region including the cerebellum and the brain stem was trimmed to approximately position upright in the specimen holder, leaving about 1 cm of the cut site above the holder and accessible for sectioning. The specimen holder with an immobilized brain was placed in a buffer tray that was filled with PBS to fully submerge the specimen. Sections of 30 µm thickness were made with the following settings: control knob of knife advancement speed set to 5 and control knob of knife vibration frequency set to 10. Cut sections were immediately transferred to PEN slides with PBS supplemented with 1 U/µl RiboLock RNase Inhibitor (Thermo Fisher Scientific, cat # EO0381).
Hematoxylin and eosin staining
Brain tissue was rehydrated through a gradient of ethanol dilutions (100%, 95%, 70%; 5 min each) and then washed in diH2O for 3 min. Hematoxylin and eosin (H&E) staining was performed using a commercial H&E staining kit (Abcam, #ab245880) following the manufacturer’s recommended protocol, with slight modification to account for slice thickness. Briefly, slides were stained with Mayer’s hematoxylin for 5 min, washed in two changes of diH2O, and differentiated in bluing solution for 15 s. Slides were again washed in two changes of diH2O, partially dehydrated through a gradient of ethanol dilutions (70%, 95%; 30 s each) and counterstained with eosin Y for 30 s. Excess stain was removed with 95% ethanol and slides were further dehydrated with 100% ethanol (3 × 2 min), cleared with xylene (3 × 5 min), and cover-slipped using Permount medium (Fisher Scientific, SP15-100). Whole slide scanning was performed using an Axio Scan.Z1 Slide Scanner (Zeiss, Oberkochen, Germany).
Nissl staining
Brain tissue was demyelinated through a series of ethanol dilutions (70%, 95%, 100%, 95%, 70%; 5 min each) and washed in diH2O for 3 min. Neuronal Nissl substance was stained using a commercial cresyl violet staining kit (Abcam, ab246816) per the manufacturer’s recommended protocol with minor modification. Briefly, slides were stained in cresyl violet staining solution for 7 min before differentiation in 95% ethanol until all visible cytoplasmic staining was removed. Slides were further dehydrated with 100% ethanol (30 dips) and cleared with xylene (3 × 5 min). Cover slipping and whole slide scanning was performed as above.
Validation of protein expression by immunofluorescence
Free floating sections were first washed in 1 × PBS for 10 min before being transferred to blocking solution (10% donkey serum, 0.3% Triton X-100 in 1 × TBS) for 1 h. The sections were incubated in primary antibody diluted in incubation solution (5% donkey serum, 0.3% Triton X-100 in 1 × TBS) overnight at 4 °C. Sections were then washed three times in 1 × TBS before being transferred to the proper secondary antibody diluted in incubation solution for one hour at room temperature. Then, sections were again washed three times in 1 × TBS before being mounted to slides. Slides were then cover-slipped using VECTASHIELD Antifade Mounting Medium with DAPI (Vector Labs, H-1200-10) and imaged with an Olympus IX-83 Confocal Laser Scanning Microscope (Olympus Life Sciences, Tokyo, Japan). Antibodies and their dilutions were as follows: mouse anti-GFAP (clone GA5, Millipore, MAB360, 1:100); rat anti-CD44 (clone IM7, Invitrogen, 14-0441-82, 1:200); donkey anti-mouse Alexafluor 488 (Invitrogen, A21202, 1:1000); and donkey anti-rat Alexafluor 488 (Invitrogen, A21208, 1:1000).
Laser microdissection
Sub-regions of hippocampus were harvested from dry sections of mice brains on PEN slides using the laser microdissection (Leica LMD7, USA) system with Leica Laser Microdissection (LMD)_V 8.3.1 software. The circumferences of regions to be harvested were drawn using the freehand tool. Cutting was done using the ‘Draw and Cut’ function with the following settings: 10× objective, 10% cut speed, 54% cut energy, 90% cut focus, aperture of 3 and specimen balance of 42. The selected tissue regions were catapulted into the cap of microfuge tubes (Non-stick RNase-Free 0.5 ml, Invitrogen, Texas, USA) placed in the tube collector at a working height of—15,356. RNAse-free 0.5 ml collection tubes. Images were taken with DMC45000 cameras, for brightfield images.
Sample preparation for proteomics
The preparation of microdissected tissue samples was optimized from a previously reported method [19]. Briefly, tissues were mixed with 10 µl lysis buffer (50% Trifluroethanol (TFE), 25 mM ammonium bicarbonate and 5 mM Dithiothreitol) and then reduced at 91 °C for 10 min in a Thermomixer (Eppendorf) and 5 min sonication. This step is followed by alkylation with 25 mM (final concentration) iodoacetamide by incubating for 20 min in the dark. The mixture was added with 50 µl of 25 mM ammonium bicarbonate and digested by trypsin at an enzyme to protein mass ratio of 1:50 overnight at 37 °C. The samples were acidified followed by drying by vacuum centrifugation. The dried pellet was resuspended in 2% acetonitrile and 0.1% formic acid.
Strong cation exchange (SCX) fractionation
For spectral library generation, the SCX stage tips were prepared, and samples were processed as previously mentioned [19]. Briefly, the tips were activated with 50 µl acetonitrile and the enriched peptides were loaded followed by their elution in 6 buffers: SCX buffer 1: 50 mM ammonium acetate, 20% acetonitrile, 0.5% formic acid; SCX buffer 2: 75 mM ammonium acetate, 20% acetonitrile, 0.5% formic acid; SCX buffer 3: 125 mM ammonium acetate, 20% acetonitrile, 0.5% formic acid; SCX buffer 4: 200 mM ammonium acetate, 20% acetonitrile, 0.5% formic acid; SCX buffer 5: 300 mM ammonium acetate, 20% acetonitrile, 0.5% formic acid; SCX buffer 6: 5% ammonium hydroxide, 80% acetonitrile. All the fractions were dried by vacuum centrifugation and suspended into 0.1% formic acid/2% acetonitrile before injection.
LC–MS/MS acquisition
Samples were measured using liquid chromatography tandem mass-spectrometry (LC–MS/MS) instrumentation consisting of an Vanquish Neo nanoLC system (Thermo Fisher Scientific) coupled via a nano-electrospray ion source (Thermo Fisher Scientific) to Orbitrap Eclipse™ Tribrid™ MS (Thermo Fisher Scientific). For generation of spectral library, the SCX fractions were loaded onto a 500 bar 300 um × 5 mm C18 trap column and PepMap™ Easy-spray™ C18 (75 um × 500 mm,) analytical column in buffer A (0.1% formic acid in water). The fractions were eluted with a linear gradient from 5 to 25% buffer B (80% acetonitrile, 0.1% formic acid) over 90 min. Following the linear separation, the system was set up to 40% buffer B over 30 min and finally set to 95% buffer B for 20 min, which was followed by re-equilibration to 5% buffer B prior to the subsequent injection. Data were acquired using data-dependent acquisition (DDA) set in 3 field asymmetry ion mobility spectrometry (FAIMS) experiments of − 40, − 55 and − 70. After adjusting each fraction to an estimated 0.5–1.0 μg on column, the fractions were measured at a cycle time of 2 s with 60 s dynamic exclusion. Precursor spectra were collected from 370–1500 m/z at 240,000 resolutions (customized automated gain control (AGC) target of 300, auto max injection time (IT)). The MS/MS were collected on + 2H to + 7H precursors achieving a standard AGC. MS/MS scans were collected using ion trap (with a dynamic IT) with an isolation width of 1.6 m/z and a NCE of 30.
For the DIA acquisitions of the individual sample with a FAIMS voltage at − 50, the peptides were injected using the same chromatography setup and column as mentioned above. The orbitrap was set to acquire 33 scan events of a biological sample pool using 120,000 precursor resolution and 30,000 fragment resolution. The AGC target was set to 3000, with auto IT, the stepped NCE was set to 25, 27, 30, and + 2H was assumed as the default charge state. The DIA was performed with one full MS event followed by 41 MS/MS windows in one cycle resulting in a cycle time of 3 s. Individual samples for proteome profiling acquisitions used single-injection DIA acquisitions was performed using 15 m/z precursor isolation windows in a staggered-window pattern with optimized window placements from 395 to 950 m/z.
Peptide identification and protein quantification
All DDA data were processed and combined into spectral library by Spectronaut (v17.3) against Mus musculus (Proteome ID UP000000589, 17,196 reviewed entries, accessed in February 2023). Protein, peptide and PSM FDR were set to 0.01. All tolerance parameter was set to dynamic. Trypsin was used as the cleaving enzyme and up to two missed cleavages were allowed. The search engine set cysteine carbamidomethylation as a fixed modification and N-acetylation and oxidation of methionine as variable modifications. The data were also searched against a decoy database so that protein identifications were accepted at a false discovery rate (FDR) of 1%. The DIA raw files were analyzed by Spectronaut using the spectral library generated. The XIC extraction window was set to dynamic. Calibration mode was automatic with MZ extraction strategy set to maximum intensity. Precursor cut off was set to 0.01 and single hit was defined by stripped sequence. Quantitation was performed at MS2 level with automatic normalization strategy used by the software.
Statistical analysis
All statistical analysis was performed in R Studio (version 2022.12.0) unless mentioned otherwise. One-way ANOVA and Fisher’s Least Significant Differences (LSD) test was used to compare groups in each hippocampal sub-regions and generate p values. Proteins were deemed significantly differentially expressed if fold change (FC) ≥ 2 or ≤ 0.5 and p < 0.05. GraphPad Prism (version 10.1.1) was used to generate graphs, heatmaps by using z scores of the protein intensities and dot plots by using number of DEPs and -log10 FDR values of each pathway. To investigate protein clusters representing different trends of expression, clustering analysis was performed using Fuzzy c means [20] within R studio, and the degree of clustering ambiguity was set to 1.5.
Results
Workflow and quantifying proteomic expression profiles on hippocampus with spatial context
The study workflow, outlined in Fig. 1, involved the microdissection of hippocampal subregions from three brains each from Sham, 1-Day post-TBI and 7-day post-TBI groups (Fig. 1A). Utilizing laser microdissection, specific hippocampal subregions were precisely collected, including the pyramidal layer (PL) extending to the stratum radiatum, stratum moleculare (SM) of cornu ammonis (CA1), and the DG divided into top (DG1) and bottom (DG2) subregions. Furthermore, CCI achieved reproducible, substantial damage in cerebral cortex post-TBI (Additional file 1: Fig. S1).
Optimization of sample processing was conducted using two methods: sodium dodecyl sulfate (SDS) and TFE initially applied to 3 pooled samples. These samples had a cut size of approximately 50,000 μm2 and a 30 μm thick brain section. Tandem mass spectrometry (MS/MS) acquisition methods were optimized using two compensation voltages (CVs) and two collision energies (stepped and fixed) (Additional file 1: Fig. S2A). Among the various methods, the TFE protocol with CV − 50 and stepped collision energy identified approximately 1500 proteins with an injection of approximately 20 ng from a cut of approximately 50,000 μm2 and a 30 μm thick brain section (S2A). This method, when applied to individual samples, yielded a comparable number of protein identifications with a ~ 20 ng injection (Additional file 1: Fig. S2A). In terms of peptide identification, the TFE lysate protocol with CV – 50 and stepped collision energy identified over 5000 peptides, which is at least five times higher than the other combinations of CVs and sample preparation strategies (Additional file 1: Fig. S2B). The difference in protein identifications between fixed and stepped collision energy for fragmentation was minimal. An investigation into the feasibility of injecting samples with different CVs was conducted because we employed FAIMS due to its increased selectivity by reducing background noise and allowing better detection of low abundant peptides [21]. Between CV − 35 and CV − 50, the latter identified a greater number of proteins (Additional file 1: Fig. S2C). In contrast, when comparing stepped and fixed collision energy, stepped CE identified more proteins (Additional file 1: Fig. S2D). Consequently, for all samples in the project, the TFE sample protocol with CV − 50 and stepped collision energy was chosen as the preferred method.
The sub-region cuts used for further DIA analysis, were between 200,000–500,000 um2 as shown in Additional file 1: Fig. S3A. For three sub-regions (SM, DG1 and DG2), the cut sizes are comparable across all 9 animals within each sub-region. In PL region, two times larger of PL cut were obtained from the Sham group compared to PL cuts from the two TBI groups (Additional file 1: Fig. S3A). Looking into the structure of PL layer, there was a notable shrinkage in TBI groups compared to other subregions (Additional file 1: Fig. S3B). After processing the samples using our optimized protocol, the analysis of all 36 samples was conducted using DIA acquisition. Additionally, six SCX fractions from each of the specified subregions were subjected to DDA to generate a spectral library which successfully identified a total of 4690 proteins (Fig. 1B). Subsequent DIA analysis of the samples identified approximately 3500 proteins from each of the targeted subregions (Fig. 1C, Additional file 2: Table S1). The Principal Component Analysis (PCA) performed on the dataset revealed discernible discrimination among brains from Sham, two TBI groups. The PCA analysis suggests that the proteomic profiles exhibited distinct variations in response to different experimental conditions, thus providing valuable insights into the molecular alterations associated with traumatic brain injury and the subsequent recovery process (Fig. 1D).
TBI induces differential protein expression pattern and signaling pathways in hippocampal sub-regions
A comprehensive analysis was conducted to map the alterations in protein expression across distinct subregions in the hippocampus-PL, SM, DG1 and DG2 at two distinct time points, 1-Day and 7-Day post-TBI, including Sham group as control.
Within the PL subregion, 24 proteins exhibited higher expression levels in both 1-Day and 7-Day post-TBI, compared to Sham, while 52 proteins were consistently downregulated at both time points following TBI compared to Sham. A subset of proteins, 10 and 69 respectively, were uniquely upregulated at 1-Day and 7-Day post-TBI. Conversely, exclusive downregulation was noted in 98 and 23 proteins at 1-Day and 7-Day post-TBI, respectively (Fig. 2A). In the SM subregion, 29 proteins increased at time points post-TBI relative to Sham, and 20 proteins decreased in comparison to Sham. In addition, 40 proteins were upregulated, and 24 proteins were demonstrated only at 7-Day post-TBI. There were also changes unique to the Day-1 post TBI, with 4 proteins upregulated and 32 downregulated (Fig. 2B). The analysis of DG1 revealed 40 proteins with higher expression and 41 proteins with lower expression at both 1-Day and 7-Day post-TBI. Furthermore, 64 and 96 proteins exhibited increased and decreased expression exclusively at 7-Day post-TBI, respectively, while 47 and 23 proteins displayed elevated and reduced expression solely at 1-Day post-TBI (Fig. 2C). In DG2, there were 42 proteins consistently increased, and 17 decreased at both time points post-TBI. Additionally, at 7-Day, 58 proteins were uniquely increased, and 38 were decreased, whereas at 1-Day, 41 proteins showed an increase, and 38 a decrease exclusively (Fig. 2D).
Using Ingenuity Pathway Analysis (IPA) to focus on metabolic and neurological pathways, our study identified distinct signaling pathway alterations across these hippocampal subregions post TBI. In PL, synaptogenesis, CREB signaling, glutaminergic and serotonin receptor signaling, along with NMDA receptor signaling, were consistently downregulated at both 1-Day and 7-Day post-TBI (Fig. 2E, Additional file 2: Table S2). This downregulation suggests impairments in synaptic formation and neuronal communication. Upregulated pathways associated with dementia and seizures indicate potential long-term neurological impacts. Disruptions in metabolic pathways were evident, with glucose metabolism disorder upregulated and fatty acid metabolism downregulated, hinting at energy production and utilization issues. In the SM subregions, acute phase response signaling, and amyloidosis were upregulated at both time points. Interestingly, metabolic pathways such as fatty acid metabolism, quantity of carbohydrates and transport of lipids were down-regulated only in 1-Day post-TBI. Notably, 7-days post-TBI saw the dysregulation of pathways such as DHCR24 and LXR/RXR signaling, along with the early onset of neurological disorder and prostaglandin metabolism, indicating secondary and potentially chronic developments (Fig. 2F, Additional file 2: Table S3).
In the two DG subregions, there are shared disruptions across DG1 and DG2, such as DHCR24 signaling, seizures, and reduced neural cells proliferation (Fig. 2G, H, Additional file 2: Tables S4, 5). At 7-Day post TBI, the decline in neural cell proliferation in both DG1 and DG2 suggests potential deficits in long-term potentiation. Additionally, each subregion exhibits unique metabolic pathway dysregulations, for example, DG1 with altered membrane lipid metabolism (Fig. 2G, Additional file 2: Table S4) and DG2 with impaired glucose metabolism, suggesting spatial metabolic variations in subregions (Fig. 2H, Additional file 2: Table S5). Pathways like amyloidosis and CREB signaling are predominant in DG1 (Fig. 2G, Additional file 2: Table S4), whereas DG2 is characterized by changes in microtubule dynamics and steroid transport (Fig. 2H, Additional file 2: Table S5). These variations emphasize spatial dynamics of protein expression and signaling pathway changes post-TBI.
Collectively, these findings contribute to an overall understanding of the temporal dynamics and regional specificity of protein expression alterations and subregion-specific nature of the signaling pathway alterations post-TBI. Other signaling and metabolic pathways marking sub-region differences at different time-points include estrogen receptor signaling, RHO GTPase cycle, MAPK signaling and serotonin receptor signaling (Additional file 2: Tables S2–5).
Trend analysis of spatial proteomics results based on DG1, DG2, CA1, and PL subregions of mouse brain slices from 1-day post-TBI and 7-day post-TBI models
We applied fuzzy c-means clustering algorithm [22] analysis on proteomic data at different time points post-brain injury to group proteins according to the changes in their expression levels that might be associated with TBI recovery or deterioration. The analysis resulted in 9 different clusters in each sub-region that allowed the construction of dynamic models capturing post-TBI processes, offering valuable insights into the development of various brain diseases (Additional file 1: Fig. S4).
To gain insights into the global change associated with TBI recovery or deterioration, acute and chronic response, we chose the following trends of expression in different sub-regions (clusters were shown in Additional file 1: Fig. S4): consistently up-regulated at both time points (cluster 1, 4 and 5 in PL, cluster 6 and 7 in SM, cluster 2 and 7 in DG1, and cluster 8 and 9 in DG2); up-regulated only at 1Â day post-TBI (cluster 8 in PL, cluster 8 in SM, cluster 5 and 6 in DG1, and cluster 3 in DG2) and up-regulated only in 7Â days post-TBI (cluster 7 in PL, cluster 9 in SM, cluster 3 in DG1, and cluster 4 and 6 in DG2). Similarly, the down-regulated clusters were grouped as consistently down-regulated in both time points (cluster 6 and 9 in PL, cluster 3, 4 and 5 in SM, cluster 4 and 9 in DG1, and cluster 2 and 7 in DG2); down-regulated only at 1Â day post-TBI (cluster 2 in PL, cluster 2 in SM, cluster 1 in DG1, and cluster 3 in DG2) and down-regulated only in 7Â days post-TBI (cluster 3 in PL, cluster 1 in SM, cluster 8 in DG1, and cluster 1 in DG2).
Subsequently, proteins were short-listed from these up and down-regulated clusters, based on their functional associations with synaptic plasticity, neurogenesis, and diseases such as Parkinson’s, dementia, and epilepsy from previous literature (Fig. 3A, Table 1) indicated by their trends of expression in the Sham vs. two TBI- time points (Fig. 3B) and spatially mapped based on their expression patterns (Fig. 3C). The roles of these proteins were also supported by studies conducted in mice models or samples collected from humans (e.g., cerebrospinal fluids) at genomics, transcriptomics and/or proteomics level (Table 1). This comprehensive analysis provides valuable insights into differentially expressed proteins in TBI and their potential links to susceptibility to subsequent neurological disorders. For example, proteins consistently up-regulated across all sub-regions, either at both time points (e.g., murinoglobulin-1 (MUG-1) and pregnancy zonal protein (PZP)) or progressively at 7 days post-TBI [e.g., glial fibrillary acidic protein (GFAP), tight junction protein ZO-1 (TJP-1), signal transducer and activator of transcription 1 (STAT-1), CD44 antigen (CD44), and ceruloplasmin (CP)] (Fig. 3C), hint at shared molecular mechanisms in response to TBI-induced stress. These proteins are also found associated with neurological disorders and their overexpression in a maximum time-point of seven days could potentially reflect the predisposition of TBI brain to the neurological disorders (Fig. 3A). To further confirm our results, GFAP and CD44 were selected for validation in tissue sections using immunofluorescence, which showed a consistent trend with proteomics results. These proteins were chosen because blood GFAP is an FDA authorized marker in individuals with suspected TBI [23] and both the proteins have distinguished roles in synaptic plasticity and neurogenesis, with reported associations in Parkinson’s and dementia [24,25,26,27,28,29]. CD44 was found to be upregulated in 7 days post-TBI compared to 1 day post-TBI and Sham in all sub-regions except SM where it showed upregulation in both the TBI time points as compared to Sham. GFAP consistently showed significant upregulation in 7 days post-TBI compared to Sham and 1 day post-TBI as shown by immunofluorescence and proteomics results (Additional file 1: Fig. S5).
Spatial variations were also evident among sub-regions with proteins being dysregulated in both or either of the time-points exclusively in each subregion (Additional file 2: Tables S6, 7). Among these proteins, several were found to influence synaptic plasticity or neurogenesis, two key factors of cognitive functioning. For instance, proteins with overexpression in both time points, such as ELAV like protein 2 (ELAVL2)[50] (Fig. 3B, Additional file 2: Table S6), chloride intracellular channel protein 1 (CLIC1) [78] (Additional file 2: Table S6) in PL, CD44 [79] and MUG-1 [43] in SM (Fig. 3B, Additional file 2: Table S6) and leucine-rich repeat protein SHOC-2 (SHOC2) [57], galectin 3 (LGALS3) [80] in DG1 (Fig. 3B, Additional file 2: Table S6), were associated with synaptic plasticity and subsequent neurodegeneration. Such examples were also evident among time-point-specific proteins such as gelsolin (GSN) [81] and protein—S100A9 [53] overexpressed in 1 day and nicastrin (NCSTN) [58] and teneurin-3 TENM3 [75] exclusively overexpressed at 7 days post-TBI in PL (Fig. 3B, Additional file 2: Table S6). A similar pattern was also observed in downregulated proteins such as voltage-dependent P/Q-type calcium channel subunit alpha-1A (CACNA1A) [82] and sodium- and chloride-dependent GABA transporter 1 (SLC6A1) [83], down-regulated in both time points from PL and SM (Fig. 3B, Additional file 2: Table S7), respectively, along with actinin-1 (ACTN1) [84] down-regulated only in 1 day post-TBI from PL, and diacylglycerol lipase-alpha (DAGLA) [85] downregulated in PL exclusively at 7 days post-TBI (Fig. 3B, Additional file 2: Table S7). These discerned spatial variations in protein expression patterns following TBI, notably those impacting synaptic plasticity and neurogenesis, contribute pivotal insights into the inclination of TBI towards subsequent neurological consequences.
Discussion
In clinical proteomics studies, tissue heterogeneity is often disregarded due to the requirement for large input material amounts [86]. Conversely, spatial proteomics proves advantageous in elucidating insights into tissue sub-region heterogeneity, thereby revealing molecular alterations specific to these sub-regions in various diseases. This study comprehensively analyzes proteomic expression profiles in spatially located subregions of the hippocampus, shedding light on molecular alterations associated with two time points (1 and 7 days) post-TBI. Sub-region protein expression patterns, unveil distinct variations at different post-TBI time points, potentially linked to long-term consequences. Following LCM, we present a sample preparation regime with MS acquisition workflow to detect protein expression differences with spatial sub-region context in mice brain tissue. Our method, which utilizes tissue sections of approximately 0.5 mm2 and a thickness of 30 μm, yields over 3500 protein IDs—a substantial improvement in resolution over a previous study on microproteomics landscape in brain injury, which identified around 2000 proteins from sections of 1 mm2 and a thickness of 200 μm [87].
The spatial analysis of protein expression patterns in distinct subregions following TBI provides valuable insights into the heterogeneous molecular response within the hippocampus. Noteworthy are the region-specific abundance changes observed in proteins like fibronectin 1 (FN1), LGALS3BP, haptoglobin (HP), and MUG-1, which were exclusively overexpressed in the SM subregion at both TBI day 1 and day 7 post-TBI. Of these proteins, LGALS3BP is a key player in astrocyte proliferation and neurosphere formation [88], and MUG-1 was previously implicated in neutrophil degranulation [89]. As also evidenced by our pathway analysis, acute phase response signaling was also enriched in SM as compared to other sub-regions. These spatially altered protein expression patterns could possibly indicate the immune cell enrichment in SM following TBI.
Conversely, the over-expression of immunoglobulin heavy constant mu (IGHM) and beta-2-microglobulin (B2M), recognized hallmarks of chronic traumatic encephalopathy [90, 91], specifically in DG2, the deepest region of the hippocampus, underscores the impact of TBI without direct mechanical injury. This observation in DG2 may signify a distinct molecular response to TBI, independent of the mechanical injury, further emphasizing the complexity and heterogeneity of the pathological processes associated with TBI. On a broader scale, the observation of several proteins, including GFAP, signal transducer and STAT1, PZP, CP, VTN, and KN, being overexpressed across all subregions following TBI is also potentially indicative of shared molecular mechanisms in response to TBI-induced stress.
The analysis of metabolic pathways in the aftermath of TBI has unveiled notable disturbances in glucose and lipid metabolism. However, a distinctive observation is the activation of cholesterol synthesis pathways, specifically involving DHCR24 and the LXR/RXR activation, in the SM on post-TBI day 7. Remarkably, these pathways persist in all deeper sub-regions at both TBI time points. This potentially suggests that regions experiencing less mechanical injury might initiate the recovery process, given the critical role of cholesterol synthesis in maintaining brain homeostasis [92]. This restorative trend is further substantiated by the observed activation of neuroglia and microtubule dynamics in DG2 on post-TBI day 7. The coordination of these processes in DG2 signifies a potential mechanism for neural repair and reconstruction, reinforcing the hypothesis that regions less affected by mechanical injury may actively engage in recovery-related pathways.
The cluster analysis has led to the identification of 20 proteins that serve as indicators of dysregulated neuronal plasticity and neurogenesis, crucial consequences of TBI, with potential implications for future neurological disorders like Parkinson’s, dementia, and epilepsy. One prominent protein among these is GFAP, consistently upregulated in all sub-regions at both time points post-TBI. The upregulation of GFAP aligns with its established role as a blood marker for TBI [23], providing robust support that our findings are reliable and confident. Moreover, GFAP has been implicated in various neurological disorders, including dementia, Parkinson’s, and epilepsy, further emphasizing its significance in the context of TBI [93, 94]. Spatial variation was observed by proteins like MUG-1, PZP, and ELAVL2 exclusively being expressed in the PL subregion while S100-A9 showed consistent upregulation in DG1 at both time points. Of these, MUG-1, known as an inhibitor of synaptic plasticity [43], potentially implies its role in disrupting crucial neuronal processes post TBI. On the other hand, PZP and ELAVL2 were known as pyramidal neuron-specific markers linked to Parkinson’s [95, 96]. This implies that the molecular changes identified post-TBI may not only indicate general neuronal dysfunction but also suggest involvement of neuro-pathways and markers in specific sub-regions.
Conversely, downregulated proteins such as SLC17A7 and synapsin-2 (SYN-2), crucial for maintaining synaptic plasticity and neurogenesis [97, 98], were identified. Previous associations of these proteins with diseases like Parkinson's, dementia, and epilepsy support their role in TBI-related processes. Notably, SLC17A7, previously found in hippocampal astrocytes [99] associated with TBI, reinforces the relevance of our discovery. In summary, the short-listed proteins provide valuable insights into the molecular signatures associated with dysregulated neuronal plasticity and neurogenesis following TBI. The links established between these proteins and neurological disorders along with their spatial distribution underscore their potential as biomarkers and therapeutic targets, advancing our understanding of the long-term consequences of traumatic brain injury.
The spatial proteome changes found by our study have the potential to advance the discovery of diagnostic and therapeutic approaches for an early-stage intervention in subsequent neurological disorders such as Parkinson’s disease or dementia, which are associated with TBI. Hippocampal neurons are known to be particularly vulnerable to TBI, leading to cognitive dysfunctions including epilepsy, due to neurodegeneration [100] and injuries to the hippocampus have been observed to impair neurogenesis [101]. Our study further provided localization of DEPs in the hippocampal microregions and their associations with injury pathways and subsequent neuroinflammatory responses. Understanding the spatial distribution of these DEPs can provide valuable insights for identifying therapeutic targets, including inflammation modulation and tissue repair promotion, aimed at ameliorating long-term neurological deficits in TBI patients [87].
Though our study is limited by sample size, it is in line with prior TBI research using mouse models [102, 103]. To properly balance the need for statistical significance with the limited throughput of the spatial proteomics method, we carefully decided to use the current sample size. However, this highlights the need for the development of a more streamlined method for spatial sample proteomics analysis. Our research sheds light on molecular alterations post-TBI during the well-recognized acute and sub-acute phases on day 1 and 7 [16]. However, it is important to acknowledge that the consequences of TBI are long-term, extending beyond these time points. Therefore, future investigations employing longitudinal samples over extended durations could provide deeper insights into the lasting effects of TBI. The scope of this study was centered on spatial proteome changes within the hippocampal sub-regions. Consequently, this focus reveals a limitation of our current study in assessing the effects on other brain regions impacted by TBI, such as the frontal and temporal lobes and/or the cerebral cortex. To address this gap, future research work will be necessary to explore the spatial resolution of protein changes throughout the entire brain post-TBI.
In conclusion, in our study, spatial proteomics unveiled nuanced molecular alterations in distinct hippocampal subregions following TBI at two time points of 1 and 7 days. Specific protein expression patterns in sub-regions highlight unique region-specific vulnerabilities, even in the deepest hippocampal sub-regions without undergoing any direct trauma. Utilizing precise LCM capture and MS, our methodology significantly advances protein identification over previous brain injury studies. Spatial analysis highlights region-specific immune cell enrichment in the SM post-TBI. Metabolic pathway analysis, especially the activation of cholesterol synthesis pathways in DG, suggests potential contribution to neurogenesis and recovery. Cluster analysis has identified 20 proteins indicative of dysregulated neuronal plasticity and neurogenesis post-TBI, including GFAP, an FDA-approved TBI blood marker, and region-specific proteins like MUG-1, PZP, and ELAVL2. This comprehensive analysis provides valuable insights into TBI's dynamic molecular landscape, offering potential biomarkers and therapeutic targets for mitigating long-term consequences.
Availability of data and materials
Raw mass spectrometry data and results are available in PRIDE repository with identifier PXD051093. All data generated and analyzed during this study are included in the additional files.
Abbreviations
- TBI:
-
Traumatic brain injury
- PL:
-
Pyrimidal layer
- SM:
-
Stratum moleculare
- CA1:
-
Cornu ammonis
- DG:
-
Dentate gyrus
- CCI:
-
Controlled cortical impact
- PFA:
-
Paraformaldehyde
- H&E:
-
Hematoxylin and eosin
- SCX:
-
Strong cation exchange chromatography
- TFE:
-
Trifluoroethanol
- LC–MS/MS:
-
Liquid chromatography–tandem mass-spectrometry
- DDA:
-
Data-dependent acquisition
- DIA:
-
Data-independent acquisition
- FAIMS:
-
Field asymmetry ion mobility spectrometry
- AGC:
-
Automatic gain control
- IT:
-
Ion injection time
- NCE:
-
Normalized collision energy
- PSM:
-
Peptide spectrum match
- FDR:
-
False discovery rate
- SDS:
-
Sodium dodecyl sulfate
- CV:
-
Compensation voltage
- PCA:
-
Principal component analysis
- IPA:
-
Ingenuity pathway analysis
References
Iacono D, Hatch K, Murphy EK, Cole RN, Post J, Leonessa F, Perl DP. Proteomic changes in the hippocampus after repeated explosive-driven blasts. J Proteome Res. 2023;23:397–408.
Jassam YN, Izzy S, Whalen M, McGavern DB, El Khoury J. Neuroimmunology of traumatic brain injury: time for a paradigm shift. Neuron. 2017;95(6):1246–65.
Mathew B, Bathla S, Williams KR, Nairn AC. Deciphering spatial protein-protein interactions in brain using proximity labeling. Mol Cell Proteomics. 2022;21(11): 100422.
Xiong Y, Zhang Y, Mahmood A, Chopp M. Investigational agents for treatment of traumatic brain injury. Expert Opin Investig Drugs. 2015;24(6):743–60.
Weston NM, Rolfe AT, Freelin AH, Reeves TM, Sun D. Traumatic brain injury modifies synaptic plasticity in newly-generated granule cells of the adult hippocampus. Exp Neurol. 2021;336: 113527.
Avchalumov Y, Mandyam CD. Plasticity in the hippocampus, neurogenesis and drugs of abuse. Brain Sci. 2021;11(3):404.
Massa F, Koehl M, Wiesner T, Grosjean N, Revest JM, Piazza PV, Abrous DN, Oliet SH. Conditional reduction of adult neurogenesis impairs bidirectional hippocampal synaptic plasticity. Proc Natl Acad Sci U S A. 2011;108(16):6644–9.
Chu HY. Synaptic and cellular plasticity in Parkinson’s disease. Acta Pharmacol Sin. 2020;41(4):447–52.
Rebelo D, Oliveira F, Abrunhosa A, Januario C, Lemos J, Castelo-Branco M. A link between synaptic plasticity and reorganization of brain activity in Parkinson’s disease. Proc Natl Acad Sci U S A. 2021;118(3):e2013962118.
Kovesdi E, Gyorgy AB, Kwon SK, Wingo DL, Kamnaksh A, Long JB, Kasper CE, Agoston DV. The effect of enriched environment on the outcome of traumatic brain injury; a behavioral, proteomics, and histological study. Front Neurosci. 2011;5:42.
Gao X, Deng-Bryant Y, Cho W, Carrico KM, Hall ED, Chen J. Selective death of newborn neurons in hippocampal dentate gyrus following moderate experimental traumatic brain injury. J Neurosci Res. 2008;86(10):2258–70.
Sun D, McGinn M, Hankins JE, Mays KM, Rolfe A, Colello RJ. Aging- and injury-related differential apoptotic response in the dentate gyrus of the hippocampus in rats following brain trauma. Front Aging Neurosci. 2013;5:95.
Abu Hamdeh S, Shevchenko G, Mi J, Musunuri S, Bergquist J, Marklund N. Proteomic differences between focal and diffuse traumatic brain injury in human brain tissue. Sci Rep. 2018;8(1):6807.
Shultz SR, Shah AD, Huang C, Dill LK, Schittenhelm RB, Morganti-Kossmann MC, Semple BD. Temporal proteomics of human cerebrospinal fluid after severe traumatic brain injury. J Neuroinflamm. 2022;19(1):291.
Sowers JL, Wu P, Zhang K, DeWitt DS, Prough DS. Proteomic changes in traumatic brain injury: experimental approaches. Curr Opin Neurol. 2018;31(6):709–17.
Zhao LR. Erratum. Brain repair by hematopoietic growth factors in the subacute phase of traumatic brain injury. J Neurosurg. 2019;130(4):1398.
Atkins CM. Decoding hippocampal signaling deficits after traumatic brain injury. Transl Stroke Res. 2011;2(4):546–55.
Liu N, Han J, Li Y, Jiang Y, Shi SX, Lok J, Whalen M, Dumont AS, Wang X. Recombinant annexin A2 inhibits peripheral leukocyte activation and brain infiltration after traumatic brain injury. J Neuroinflamm. 2021;18(1):173.
Maity S, Mayer MG, Shu Q, Linh H, Bao D, Blair RV, He Y, Lyon CJ, Hu TY, Fischer T, et al. Cerebrospinal fluid protein markers indicate neuro-damage in SARS-CoV-2-infected nonhuman primates. Mol Cell Proteomics. 2023;22(4): 100523.
Bhimavarapu U, Chintalapudi N, Battineni G. Brain tumor detection and categorization with segmentation of improved unsupervised clustering approach and machine learning classifier. Bioengineering (Basel). 2024;11(3):266.
Reilly L, Lara E, Ramos D, Li Z, Pantazis CB, Stadler J, Santiana M, Roberts J, Faghri F, Hao Y, et al. A fully automated FAIMS-DIA mass spectrometry-based proteomic pipeline. Cell Rep Methods. 2023;3(10): 100593.
Gong M, Liang Y, Shi J. Fuzzy C-means clustering with local information and kernel metric for image segmentation. IEEE Trans Image Process Publ IEEE Signal Process Soc. 2012;22:573–84.
Abdelhak A, Foschi M, Abu-Rumeileh S, Yue JK, D’Anna L, Huss A, Oeckl P, Ludolph AC, Kuhle J, Petzold A, et al. Blood GFAP as an emerging biomarker in brain and spinal cord disorders. Nat Rev Neurol. 2022;18(3):158–72.
Roszkowska M, Skupien A, Wojtowicz T, Konopka A, Gorlewicz A, Kisiel M, Bekisz M, Ruszczycki B, Dolezyczek H, Rejmak E, et al. CD44: a novel synaptic cell adhesion molecule regulating structural and functional plasticity of dendritic spines. Mol Biol Cell. 2016;27(25):4055–66.
Su W, Foster SC, Xing R, Feistel K, Olsen RH, Acevedo SF, Raber J, Sherman LS. CD44 transmembrane receptor and hyaluronan regulate adult hippocampal neural stem cell quiescence and differentiation. J Biol Chem. 2017;292(11):4434–45.
Shir D, Graff-Radford J, Hofrenning EI, Lesnick TG, Przybelski SA, Lowe VJ, Knopman DS, Petersen RC, Jack CR Jr, Vemuri P, et al. Association of plasma glial fibrillary acidic protein (GFAP) with neuroimaging of Alzheimer’s disease and vascular pathology. Alzheimers Dement (Amst). 2022;14(1): e12291.
Mochol M, Taubøll E, Aukrust P, Ueland T, Andreassen OA, Svalheim S, Ben-Menachem E. Serum markers of neuronal damage and astrocyte activity in patients with chronic epilepsy: elevated levels of glial fibrillary acidic protein. Acta Neurol Scand. 2023;2023:1–5.
Johnson ECB, Dammer EB, Duong DM, Ping L, Zhou M, Yin L, Higginbotham LA, Guajardo A, White B, Troncoso JC, et al. Large-scale proteomic analysis of Alzheimer’s disease brain and cerebrospinal fluid reveals early changes in energy metabolism associated with microglia and astrocyte activation. Nat Med. 2020;26(5):769–80.
Li S, Shao H, Chang L. The important role of perituberal tissue in epileptic patients with tuberous sclerosis complex by the transcriptome analysis. Biomed Res Int. 2020;2020:4980609.
Wruck W, Adjaye J. Meta-analysis of human prefrontal cortex reveals activation of GFAP and decline of synaptic transmission in the aging brain. Acta Neuropathol Commun. 2020;8(1):26.
Zwirner J, Lier J, Franke H, Hammer N, Matschke J, Trautz F, Tse R, Ondruschka B. GFAP positivity in neurons following traumatic brain injuries. Int J Legal Med. 2021;135(6):2323–33.
Liu Y, Namba T, Liu J, Suzuki R, Shioda S, Seki T. Glial fibrillary acidic protein-expressing neural progenitors give rise to immature neurons via early intermediate progenitors expressing both glial fibrillary acidic protein and neuronal markers in the adult hippocampus. Neuroscience. 2010;166(1):241–51.
Che N, Ou R, Li C, Zhang L, Wei Q, Wang S, Jiang Q, Yang T, Xiao Y, Lin J, et al. Plasma GFAP as a prognostic biomarker of motor subtype in early Parkinson’s disease. NPJ Parkinsons Dis. 2024;10(1):48.
Kronenberg G, Gertz K, Baldinger T, Kirste I, Eckart S, Yildirim F, Ji S, Heuser I, Schröck H, Hörtnagl H, et al. Impact of actin filament stabilization on adult hippocampal and olfactory bulb neurogenesis. J Neurosci. 2010;30(9):3419–31.
Jiang C, Hopfner F, Katsikoudi A, Hein R, Catli C, Evetts S, Huang Y, Wang H, Ryder JW, Kuhlenbaeumer G, et al. Serum neuronal exosomes predict and differentiate Parkinson’s disease from atypical parkinsonism. J Neurol Neurosurg Psychiatry. 2020;91(7):720–9.
Gamez-Valero A, Campdelacreu J, Rene R, Beyer K, Borras FE. Comprehensive proteomic profiling of plasma-derived extracellular vesicles from dementia with lewy bodies patients. Sci Rep. 2019;9(1):13282.
Bitsika V, Duveau V, Simon-Areces J, Mullen W, Roucard C, Makridakis M, Mermelekas G, Savvopoulos P, Depaulis A, Vlahou A. High-throughput LC-MS/MS proteomic analysis of a mouse model of mesiotemporal lobe epilepsy predicts microglial activation underlying disease development. J Proteome Res. 2016;15(5):1546–62.
Bendriem RM, Singh S, Aleem AA, Antonetti DA, Ross ME. Tight junction protein occludin regulates progenitor self-renewal and survival in developing cortex. Elife. 2019;8:e49376.
Forouzandeh M, Mostafavi H, Ghasemloo E, Mohammadi P, Hosseini M, Eskandari M. Increased expression of tight junction proteins and blood-brain barrier integrity in mcao rats following injection of miR-149-5p. Int J Mol Cell Med. 2022;11(3):223–35.
Yan BC, Xu P, Gao M, Wang J, Jiang D, Zhu X, Won MH, Su PQ. Changes in the blood-brain barrier function are associated with hippocampal neuron death in a kainic acid mouse model of epilepsy. Front Neurol. 2018;9:775.
Choudhury ME, Ozaki S, Miyaue N, Matsuura T, Mikami K, Islam A, Kubo M, Ando R, Yano H, Kunieda T, et al. Chloride intracellular channel protein 2 promotes microglial invasion: a link to microgliosis in the Parkinson’s disease brain. Brain Sci. 2022;13(1):55.
Shen KF, Yang XL, Liu GL, Zhu G, Wang ZK, Shi XJ, Wang TT, Wu ZF, Lv SQ, Liu SY, et al. The role of voltage-gated chloride channels in the epileptogenesis of temporal lobe epilepsy. EBioMedicine. 2021;70: 103537.
Kato K, Kishi T, Kamachi T, Akisada M, Oka T, Midorikawa R, Takio K, Dohmae N, Bird PI, Sun J, et al. Serine proteinase inhibitor 3 and murinoglobulin I are potent inhibitors of neuropsin in adult mouse brain. J Biol Chem. 2001;276(18):14562–71.
Sun J, Jiang T, Gu F, Ma D, Liang J. TMT-based proteomic analysis of plasma from children with rolandic epilepsy. Dis Markers. 2020;2020:8840482.
Garcia-Revilla J, Boza-Serrano A, Jin Y, Vadukul DM, Soldan-Hidalgo J, Camprubi-Ferrer L, Garcia-Cruzado M, Martinsson I, Klementieva O, Ruiz R, et al. Galectin-3 shapes toxic alpha-synuclein strains in Parkinson’s disease. Acta Neuropathol. 2023;146(1):51–75.
Nicolas CS, Peineau S, Amici M, Csaba Z, Fafouri A, Javalet C, Collett VJ, Hildebrandt L, Seaton G, Choi SL, et al. The Jak/STAT pathway is involved in synaptic plasticity. Neuron. 2012;73(2):374–90.
Kim JH, Choi DJ, Jeong HK, Kim J, Kim DW, Choi SY, Park SM, Suh YH, Jou I, Joe EH. DJ-1 facilitates the interaction between STAT1 and its phosphatase, SHP-1, in brain microglia and astrocytes: a novel anti-inflammatory function of DJ-1. Neurobiol Dis. 2013;60:1–10.
Hsu WL, Ma YL, Hsieh DY, Liu YC, Lee EH. STAT1 negatively regulates spatial memory formation and mediates the memory-impairing effect of Abeta. Neuropsychopharmacology. 2014;39(3):746–58.
Bolognani F, Qiu S, Tanner DC, Paik J, Perrone-Bizzozero NI, Weeber EJ. Associative and spatial learning and memory deficits in transgenic mice overexpressing the RNA-binding protein HuD. Neurobiol Learn Mem. 2007;87(4):635–43.
Berto S, Usui N, Konopka G, Fogel BL. ELAVL2-regulated transcriptional and splicing networks in human neurons link neurodevelopment and autism. Hum Mol Genet. 2016;25(12):2451–64.
Noureddine MA, Qin XJ, Oliveira SA, Skelly TJ, van der Walt J, Hauser MA, Pericak-Vance MA, Vance JM, Li YJ. Association between the neuron-specific RNA-binding protein ELAVL4 and Parkinson disease. Hum Genet. 2005;117(1):27–33.
Diouf I, Bush AI, Ayton S. Alzheimer’s disease Neuroimaging I: Cerebrospinal fluid ceruloplasmin levels predict cognitive decline and brain atrophy in people with underlying beta-amyloid pathology. Neurobiol Dis. 2020;139: 104810.
Qiu X, Ping S, Kyle M, Longo J, Chin L, Zhao L-R. S100 calcium-binding protein A9 knockout contributes to neuroprotection and functional improvement after traumatic brain injury. J Neurotrauma. 2020;37(7):950–65.
Tian Y, Cao R, Che B, Sun D, Tang Y, Jiang L, Bai Q, Liu Y, Morozova-Roche LA, Zhang C. Proinflammatory S100A9 regulates differentiation and aggregation of neural stem cells. ACS Chem Neurosci. 2020;11(21):3549–56.
Wang X, Liu H, Li J, Jiang J, Li T. Association study of S100A9 gene polymorphisms with Parkinson’s disease risk and age of disease onset. Acta Neurol Belg. 2024. https://doi.org/10.1007/s13760-024-02486-0.
Wang C, Klechikov AG, Gharibyan AL, Warmlander SK, Jarvet J, Zhao L, Jia X, Narayana VK, Shankar SK, Olofsson A, et al. The role of pro-inflammatory S100A9 in Alzheimer’s disease amyloid-neuroinflammatory cascade. Acta Neuropathol. 2014;127(4):507–22.
Norcross RG, Abdelmoti L, Rouchka EC, Andreeva K, Tussey O, Landestoy D, Galperin E. Shoc2 controls ERK1/2-driven neural crest development by balancing components of the extracellular matrix. Dev Biol. 2022;492:156–71.
Lee SH, Sharma M, Sudhof TC, Shen J. Synaptic function of nicastrin in hippocampal neurons. Proc Natl Acad Sci U S A. 2014;111(24):8973–8.
Tabuchi K, Chen G, Sudhof TC, Shen J. Conditional forebrain inactivation of nicastrin causes progressive memory impairment and age-related neurodegeneration. J Neurosci. 2009;29(22):7290–301.
Hernández-Villamizar LF, GarcÃa-Ayllón MS, Pérez-González R, Zetterberg H, Blennow K, Sáez-Valero J. Levels of nicastrin species are increased in Alzheimer’s disease cerebrospinal fluid. Alzheimers Dement. 2023;19(S24):e082732.
Kodam A, Ourdev D, Maulik M, Hariharakrishnan J, Banerjee M, Wang Y, Kar S. A role for astrocyte-derived amyloid beta peptides in the degeneration of neurons in an animal model of temporal lobe epilepsy. Brain Pathol. 2019;29(1):28–44.
Smajic S, Prada-Medina CA, Landoulsi Z, Ghelfi J, Delcambre S, Dietrich C, Jarazo J, Henck J, Balachandran S, Pachchek S, et al. Single-cell sequencing of human midbrain reveals glial activation and a Parkinson-specific neuronal state. Brain. 2022;145(3):964–78.
Schnizler MK, Schnizler K, Zha XM, Hall DD, Wemmie JA, Hell JW, Welsh MJ. The cytoskeletal protein alpha-actinin regulates acid-sensing ion channel 1a through a C-terminal interaction. J Biol Chem. 2009;284(5):2697–705.
Venugopal AK, Sameer Kumar GS, Mahadevan A, Selvan LD, Marimuthu A, Dikshit JB, Tata P, Ramachandra Y, Chaerkady R, Sinha S, et al. Transcriptomic profiling of medial temporal lobe epilepsy. J Proteomics Bioinform. 2012;5(2):1000210.
Terpolilli NA, Dolp R, Waehner K, Schwarzmaier SM, Rumbler E, Todorov B, Ferrari MD, van den Maagdenberg A, Plesnila N. Ca(V)2.1 channel mutations causing familial hemiplegic migraine type 1 increase the susceptibility for cortical spreading depolarizations and seizures and worsen outcome after experimental traumatic brain injury. Elife. 2022;11:e74923.
Jouvenceau A, Eunson LH, Spauschus A, Ramesh V, Zuberi SM, Kullmann DM, Hanna MG. Human epilepsy associated with dysfunction of the brain P/Q-type calcium channel. Lancet. 2001;358(9284):801–7.
Schurman LD, Carper MC, Moncayo LV, Ogasawara D, Richardson K, Yu L, Liu X, Poklis JL, Liu QS, Cravatt BF, et al. Diacylglycerol lipase-alpha regulates hippocampal-dependent learning and memory processes in mice. J Neurosci. 2019;39(30):5949–65.
Ludanyi A, Eross L, Czirjak S, Vajda J, Halasz P, Watanabe M, Palkovits M, Magloczky Z, Freund TF, Katona I. Downregulation of the CB1 cannabinoid receptor and related molecular elements of the endocannabinoid system in epileptic human hippocampus. J Neurosci. 2008;28(12):2976–90.
Zander JF, Munster-Wandowski A, Brunk I, Pahner I, Gomez-Lira G, Heinemann U, Gutierrez R, Laube G, Ahnert-Hilger G. Synaptic and vesicular coexistence of VGLUT and VGAT in selected excitatory and inhibitory synapses. J Neurosci. 2010;30(22):7634–45.
Johannesen KM, Gardella E, Linnankivi T, Courage C, de Saint MA, Lehesjoki AE, Mignot C, Afenjar A, Lesca G, Abi-Warde MT, et al. Defining the phenotypic spectrum of SLC6A1 mutations. Epilepsia. 2018;59(2):389–402.
Sun J, Bronk P, Liu X, Han W, Südhof TC. Synapsins regulate use-dependent synaptic plasticity in the calyx of Held by a Ca2+/calmodulin-dependent pathway. Proc Natl Acad Sci. 2006;103(8):2880–5.
Barbieri R, Contestabile A, Grazia Ciardo M, Forte N, Marte A, Baldelli P, Benfenati F, Onofri F. Synapsin I and Synapsin II regulate neurogenesis in the dentate gyrus of adult mice. Oncotarget. 2018;9(27):18760.
Stern S, Lau S, Manole A, Rosh I, Percia MM, Ben Ezer R, Shokhirev MN, Qiu F, Schafer S, Mansour AA, et al. Reduced synaptic activity and dysregulated extracellular matrix pathways in midbrain neurons from Parkinson’s disease patients. NPJ Parkinsons Dis. 2022;8(1):103.
Larson ME, Greimel SJ, Amar F, LaCroix M, Boyle G, Sherman MA, Schley H, Miel C, Schneider JA, Kayed R, et al. Selective lowering of synapsins induced by oligomeric alpha-synuclein exacerbates memory deficits. Proc Natl Acad Sci U S A. 2017;114(23):E4648–57.
Zhang X, Lin PY, Liakath-Ali K, Sudhof TC. Teneurins assemble into presynaptic nanoclusters that promote synapse formation via postsynaptic non-teneurin ligands. Nat Commun. 2022;13(1):2297.
Xu C, Li Z, Lyu C, Hu Y, McLaughlin CN, Wong KKL, Xie Q, Luginbuhl DJ, Li H, Udeshi ND et al: Molecular and Cellular Mechanisms of Teneurin Signaling in Synaptic Partner Matching. bioRxiv 2024:2024.2002.2023.581689.
Bastias-Candia S, Martinez M, Zolezzi JM, Inestrosa NC. Wnt signaling upregulates teneurin-3 expression via canonical and non-canonical Wnt pathway crosstalk. Front Neurosci. 2019;13:505.
Cuentas-Condori A, Chen S, Krout M, Gallik KL, Tipps J, Gailey C, Flautt L, Kim H, Mulcahy B, Zhen M, et al. The epithelial Na(+) channel UNC-8 promotes an endocytic mechanism that recycles presynaptic components to new boutons in remodeling neurons. Cell Rep. 2023;42(11): 113327.
Dzwonek J, Wilczynski GM. CD44: molecular interactions, signaling and functions in the nervous system. Front Cell Neurosci. 2015;9:175.
Soares LC, Al-Dalahmah O, Hillis J, Young CC, Asbed I, Sakaguchi M, O’Neill E, Szele FG. Novel galectin-3 roles in neurogenesis, inflammation and neurological diseases. Cells. 2021;10(11):3047.
Hlushchenko I, Hotulainen P. Chemical LTD, but not LTP, induces transient accumulation of gelsolin in dendritic spines. Biol Chem. 2019;400(9):1129–39.
Dedic N, Pohlmann ML, Richter JS, Mehta D, Czamara D, Metzger MW, Dine J, Bedenk BT, Hartmann J, Wagner KV, et al. Cross-disorder risk gene CACNA1C differentially modulates susceptibility to psychiatric disorders during development and adulthood. Mol Psychiatry. 2018;23(3):533–43.
Mermer F, Poliquin S, Rigsby K, Rastogi A, Shen W, Romero-Morales A, Nwosu G, McGrath P, Demerast S, Aoto J, et al. Common molecular mechanisms of SLC6A1 variant-mediated neurodevelopmental disorders in astrocytes and neurons. Brain. 2021;144(8):2499–512.
Pinho J, Marcut C, Fonseca R. Actin remodeling, the synaptic tag and the maintenance of synaptic plasticity. IUBMB Life. 2020;72(4):577–89.
Lee D, Kim E, Tanaka-Yamamoto K. Diacylglycerol kinases in the coordination of synaptic plasticity. Front Cell Dev Biol. 2016;4:92.
Davis S, Scott C, Ansorge O, Fischer R. Development of a sensitive, scalable method for spatial, cell-type-resolved proteomics of the human brain. J Proteome Res. 2019;18(4):1787–95.
Mallah K, Quanico J, Raffo-Romero A, Cardon T, Aboulouard S, Devos D, Kobeissy F, Zibara K, Salzet M, Fournier I. Mapping spatiotemporal microproteomics landscape in experimental model of traumatic brain injury unveils a link to Parkinson’s disease. Mol Cell Proteomics. 2019;18(8):1669–82.
Sirko S, Schichor C, Della Vecchia P, Metzger F, Sonsalla G, Simon T, Burkle M, Kalpazidou S, Ninkovic J, Masserdotti G, et al. Injury-specific factors in the cerebrospinal fluid regulate astrocyte plasticity in the human brain. Nat Med. 2023;29(12):3149–61.
Chalise U, Daseke MJ 2nd, Kalusche WJ, Konfrst SR, Rodriguez-Paar JR, Flynn ER, Cook LM, Becirovic-Agic M, Lindsey ML. Macrophages secrete murinoglobulin-1 and galectin-3 to regulate neutrophil degranulation after myocardial infarction. Mol Omics. 2022;18(3):186–95.
Gutierrez-Quiceno L, Dammer EB, Johnson AG, Webster JA, Shah R, Duong D, Yin L, Seyfried NT, Alvarez VE, Stein TD, et al. A proteomic network approach resolves stage-specific molecular phenotypes in chronic traumatic encephalopathy. Mol Neurodegener. 2021;16(1):40.
Huo Q, Dong W, Gao Y, Zhang Y, Liu X, Yang L, Nan D, Yang J. Effect of beta2-microglobulin in evaluating the severity and prognosis of brain injury: a clinical study. BMC Neurol. 2022;22(1):327.
Zhang J, Liu Q. Cholesterol metabolism and homeostasis in the brain. Protein Cell. 2015;6(4):254–64.
Cronje HT, Liu X, Odden MC, Moseholm KF, Seshadri S, Satizabal CL, Lopez OL, Bis JC, Djousse L, Fohner AE, et al. Serum NfL and GFAP are associated with incident dementia and dementia mortality in older adults: The cardiovascular health study. Alzheimers Dement. 2023;19(12):5672–80.
Tang Y, Han L, Li S, Hu T, Xu Z, Fan Y, Liang X, Yu H, Wu J, Wang J. Plasma GFAP in Parkinson’s disease with cognitive impairment and its potential to predict conversion to dementia. NPJ Parkinsons Dis. 2023;9(1):23.
Nijholt DAT, Ijsselstijn L, van der Weiden MM, Zheng P-P, Sillevis Smitt PAE, Koudstaal PJ, Luider TM, Kros JM. Pregnancy zone protein is increased in the Alzheimer’s disease brain and associates with senile plaques. J Alzheimers Dis. 2015;46:227–38.
Pastic A, Negeri O, Ravel-Chapuis A, Savard A, Trung MT, Palidwor G, Guo H, Marcogliese P, Taylor J, Okano H et al. LRRK2 phosphorylates neuronal elav RNA-binding proteins to regulate phenotypes relevant to Parkinson’s disease. In: bioRxiv; 2022.
Medrihan L, Cesca F, Raimondi A, Lignani G, Baldelli P, Benfenati F. Synapsin II desynchronizes neurotransmitter release at inhibitory synapses by interacting with presynaptic calcium channels. Nat Commun. 2013;4:1512.
Roberts BM, Doig NM, Brimblecombe KR, Lopes EF, Siddorn RE, Threlfell S, Connor-Robson N, Bengoa-Vergniory N, Pasternack N, Wade-Martins R, et al. GABA uptake transporters support dopamine release in dorsal striatum with maladaptive downregulation in a parkinsonism model. Nat Commun. 2020;11(1):4958.
Piao CS, Holloway AL, Hong-Routson S, Wainwright MS. Depression following traumatic brain injury in mice is associated with down-regulation of hippocampal astrocyte glutamate transporters by thrombin. J Cereb Blood Flow Metab. 2019;39(1):58–73.
Shetty AK. Hippocampal injury-induced cognitive and mood dysfunction, altered neurogenesis, and epilepsy: can early neural stem cell grafting intervention provide protection? Epilepsy Behav. 2014;38:117–24.
Redell JB, Maynard ME, Underwood EL, Vita SM, Dash PK, Kobori N. Traumatic brain injury and hippocampal neurogenesis: functional implications. Exp Neurol. 2020;331: 113372.
Nguyen TT, Dammer EB, Owino SA, Giddens MM, Madaras NS, Duong DM, Seyfried NT, Hall RA. Quantitative proteomics reveal an altered pattern of protein expression in brain tissue from mice lacking GPR37 and GPR37L1. J Proteome Res. 2020;19(2):744–55.
Park J, Lee SH, Shin D, Kim Y, Kim YS, Seong MY, Lee JJ, Seo HG, Cho WS, Ro YS, et al. Multiplexed quantitative proteomics reveals proteomic alterations in two rodent traumatic brain injury models. J Proteome Res. 2024;23(1):249–63.
Acknowledgements
We thank Dr. Deniz Ucar to assist us in operating the laser microdissection microscope.
Funding
This study was supported by NIH High-End Instrumentation Award (S10OD032453) to J.F. for Orbitrap Eclipse Tribrid Mass Spectrometer system, housed in the Tulane Proteomics Core. The authors would like to acknowledge the Tulane COBRE for Clinical and Translational Research in Cardiometabolic Diseases and NIH grants (P20GM109036) to J.F. and N.L.; American Heart Association Career Development Award (23CDA1055341) to N.L.; NIH grant (R21NS130542) to J. C.; NIH grant (R01NS126503-01A1) to X.W; NIH grant (R21AI169583) to B.N; The study was also supported by Tulane Brain Institute microscopy core funded by a Louisiana Board of Regents Enhancement Grant LEQSF(2018-23)-ENH-DE-15. Support was provided through Jim and Betty Karam postdoctoral scholarship funds, Tulane University to S.M.
Author information
Authors and Affiliations
Contributions
S.M., and J.F. contributed to experimental design; M.D.K., X.W., and N.L. performed mice perfusion and histological staining; S.M., Y.H., L.O.V., A.N.T. and M.J.G. performed tissue sectioning and laser capture microdissection. S.M., Y.H., N.L., and J.F. analyzed the data and prepare the figures; S.M., Y.H., and J.F. performed proteomics analysis; S.M., Y.H., X.X., J.C. performed the pathway analysis. S.M., J.F., and B.N., organized the data and drafted the manuscript; All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
All animal experiments were carried out following National Institutes of Health (NIH) guidelines and all protocols were approved by the Institution Animal Care and Use Committee of the Tulane University School of Medicine (Protocol #2000).
Consent for publication
This manuscript has been read and approved by all the authors to publish and is not submitted or under consideration for publication elsewhere.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Three brains, sham, 1 and 7 days post-TBI were stained with Nissl and hematoxylin and eosin (H & E) stain to show the morphological damage to the cortex of TBI as compared to sham. Figure S2. A. Comparison was showed between different sample processing strategy and MS acquisition methods. Of all the methods, sample processing with trifluroethanol followed by FAIMS CV voltage of -50 and stepped collision energy yielded the maximum number of proteins. B. The same protocol yielded a corresponding peptide number of 5000 with the individual samples peptide number ranging between 6000–8000. C. Venn diagram depicts the improvement in protein yield with CV 50 as compared to CV 35 and D. using stepped and fixed collision energies. Figure S3. A. Microdissection by LCM yielded approximately a cut of 40,000 um2 in PL, which was reduced to almost half in post TBI samples. On the contrary, area sections stayed consistent in rest of the sub- regions. B. The Pyramidal layer shows shrinkage on 1 and 7 days post-TBI. Figure S4. Clustering of stimulation time expression patterns in identified proteins from (A) DG1, (B) DG2, (C) SM1 and (D) PL subregions with 9 kinds of rising or falling patterns using the fuzzy c-means algorithm. Warm and cold colors indicate low and high deviation from the consensus profile, respectively. Figure S5. Immunofluorescence was performed by using CD44 and GFAP antibody on mice brain sections belonging to Sham, and 1-day post-TBI and 7-day post-TBI groups. CD44 showed significant upregulation at day 7 for all sub-regions except SM which also showed upregulation in day 1. This observation was also supported by our proteomics results from different sub-regions at different time-points. We observed similar correspondence between the proteomics and the IHC dataset of GFAP which showed prominent expression on day-7 as compared to day-1 and sham. ns – not significant; * p < 0.05.
Additional file 2: Table S1.
Proteins identified by Spectronaut. Table S2. IPA analysis of differentially expressed proteins in PL. Table S3. IPA analysis of differentially expressed proteins in SM. Table S4. IPA analysis of differentially expressed proteins in DG1. Table S5. IPA analysis of differentially expressed proteins in DG2. Table S6. Up-regulated proteins spatially identified in each sub-region and either at both or at each time-point. Table S7. Down-regulated proteins spatially identified in each sub-region and either at both or at each time-point.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Maity, S., Huang, Y., Kilgore, M.D. et al. Mapping dynamic molecular changes in hippocampal subregions after traumatic brain injury through spatial proteomics. Clin Proteom 21, 32 (2024). https://doi.org/10.1186/s12014-024-09485-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12014-024-09485-6