The Gene Coexpression Analysis Identifies Functional Modules Dynamically Changed After Traumatic Brain Injury

Selective neuronal expression of progranulin is sufficient to provide neuroprotective and anti-inflammatory effects after traumatic brain injury

Progranulin (PGRN), which is produced in neurons and microglia, is a neurotrophic and anti-inflammatory glycoprotein. Human loss-of-function mutations cause frontotemporal dementia, and PGRN knockout (KO) mice are a model for dementia. In addition, PGRN KO mice exhibit severe phenotypes in models of traumatic or ischemic central nervous system (CNS) disorders, including traumatic brain injury (TBI). It is unknown whether restoration of progranulin expression in neurons (and not in microglia) might be sufficient to prevent excessive TBI-evoked brain damage. To address this question, we generated mice with Nestin-Cre-driven murine PGRN expression in a PGRN KO line (PGRN-KONestinGrn) to rescue PGRN in neurons. PGRN expression analysis in primary CNS cell cultures from naïve mice and in (non-) injured brain tissue from PGRN-KONestinGrn revealed expression of PGRN in neurons but not in microglia. After experimental TBI, examination of the structural brain damage at 5 days post-injury (dpi) showed that the TBI-induced loss of brain tissue and hippocampal neurons was exacerbated in PGRN-KOGrnflfl mice (PGRN knockout with the mGrn fl-STOP-fl allele, Cre-negative), as expected, whereas the tissue damage in PGRN-KONestinGrn mice was similar to that in PGRN-WT mice. Analysis of CD68+ immunofluorescent microglia and Cd68 mRNA expression showed that excessive microglial activation was rescued in PGRN-KONestinGrn mice, and the correlation of brain injury with Cd68 expression suggested that Cd68 was a surrogate marker for excessive brain injury caused by PGRN deficiency. The results show that restoring neuronal PGRN expression was sufficient to rescue the exacerbated neuropathology of TBI caused by PGRN deficiency, even in the absence of microglial PGRN. Hence, endogenous microglial PGRN expression was not essential for the neuroprotective or anti-inflammatory effects of PGRN after TBI in this study.

Regulatory T cells subgroups in the tumor microenvironment cannot be overlooked: Their involvement in prognosis and treatment strategy in melanoma

BACKGROUND

Melanoma, the most lethal form of skin cancer, presents substantial challenges despite effective surgical interventions for in situ lesions. Regulatory T cells (Tregs) wield a pivotal immunomodulatory influence within the tumor microenvironment, yet their impact on melanoma prognosis and direct molecular interactions with melanoma cells remain elusive. This investigation employs single-cell analysis to unveil the intricate nature of Tregs in human melanoma.

METHODS

Single-cell RNA and bulk sequencing data, alongside clinical information, were obtained from public repositories. Initially, GO and GSEA analyses were employed to delineate functional disparities among distinct cell subsets. Pseudotime and cell-cell interconnection analyses were conducted, followed by an endeavor to construct a prognostic model grounded in Treg-associated risk scores. This model's efficacy was demonstrated via PCA and K-M analyses, with multivariate Cox regression affirming its independent prognostic value in melanoma patients. Furthermore, immune infiltration analysis, immune checkpoint gene expression scrutiny, and drug sensitivity assessments were performed to ascertain the clinical relevance of this prognostic model.

RESULTS

Following batch effect correction, 80 025 cells partitioned into 31 clusters, encompassing B cells, plasma cells, endothelial cells, fibroblasts, melanoma cells, monocytes, macrophages, and T_NK cells. Within these, 4240 CD4+ T cells were subclassified into seven distinct types. Functional analysis underscored the immunomodulatory function of Tregs within the melanoma tumor microenvironment, elucidating disparities among Treg subpopulations. Notably, the ITGB2 signaling pathway emerged as a plausible molecular nexus linking Tregs to melanoma cells. Our prognostic signature exhibited robust predictive capacities for melanoma prognosis and potential implications in evaluating immunotherapy response.

CONCLUSION

Tregs exert a critical role in immune suppression within the melanoma tumor microenvironment, revealing a potential molecular-level association with melanoma cells. Our innovative Treg-centered signature introduces a promising prognostic marker for melanoma, holding potential for future clinical prognostic assessments.

Single-cell RNA sequencing explored potential therapeutic targets by revealing the tumor microenvironment of neuroblastoma and its expression in cell death

BACKGROUND

Neuroblastoma (NB) is the most common extracranial solid tumor in childhood and is closely related to the early development and differentiation of neuroendocrine (NE) cells. The disease is mainly represented by high-risk NB, which has the characteristics of high mortality and difficult treatment. The survival rate of high-risk NB patients is not ideal. In this article, we not only conducted a comprehensive study of NB through single-cell RNA sequencing (scRNA-seq) but also further analyzed cuproptosis, a new cell death pathway, in order to find clinical treatment targets from a new perspective.

MATERIALS AND METHODS

The Seurat software was employed to process the scRNA-seq data. This was followed by the utilization of GO enrichment analysis and GSEA to unveil pertinent enriched pathways. The inferCNV software package was harnessed to investigate chromosomal copy number variations. pseudotime analyses involved the use of Monocle 2, CytoTRACE, and Slingshot software. CellChat was employed to analyze the intercellular communication network for NB. Furthermore, PySCENIC was deployed to review the profile of transcription factors.

RESULT

Using scRNA-seq, we studied cells from patients with NB. NE cells exhibited superior specificity in contrast to other cell types. Among NE cells, C1 PCLAF + NE cells showed a close correlation with the genesis and advancement of NB. The key marker genes, cognate receptor pairing, developmental trajectories, metabolic pathways, transcription factors, and enrichment pathways in C1 PCLAF + NE cells, as well as the expression of cuproptosis in C1 PCLAF + NE cells, provided new ideas for exploring new therapeutic targets for NB.

CONCLUSION

The results revealed the specificity of malignant NE cells in NB, especially the key subset of C1 PCLAF + NE cells, which enhanced our understanding of the key role of the tumor microenvironment in the complexity of cancer progression. Of course, cell death played an important role in the progression of NB, which also promoted our research on new targets. The scrutiny of these findings proved advantageous in uncovering innovative therapeutic targets, thereby bolstering clinical interventions.

Unveiling the cellular landscape: insights from single-cell RNA sequencing in multiple myeloma

Objective

The aim of this research was to gain a thorough understanding of the processes involved in cell communication and discover potential indicators for treating multiple myeloma (MM) through the use of single-cell RNA sequencing (scRNA-seq). And explored the expression of multiple myeloma-related subgroups on metal ion-related pathways to explore the relationship between MM and metal ions.

Methods

We performed a fair examination using single-cell RNA sequencing on 32 bone marrow specimens collected from 22 individuals at different points of MM advancement and 9 individuals without any health issues. To analyze the scRNA-seq data, we employed advanced computational algorithms, including Slingshot, Monocle2, and other methodologies. Specifically, Slingshot and Monocle2 enabled us to simulate the biological functionalities of different cell populations and map trajectories of cell developmental pathways. Additionally, we utilized the UMAP algorithm, a powerful dimension reduction technique, to cluster cells and identify genes that were differentially expressed across clusters.

Results

Our study revealed distinct gene expression patterns and molecular pathways within each patient, which exhibited associations with disease progression. The analysis provided insights into the tumor microenvironment (TME), intra- and inter-patient heterogeneity, and cell-cell interactions mediated by ligand-receptor signaling. And found that multiple myeloma-related subgroups were expressed higher levels in MMP and TIMP pathways, there were some associations.

Conclusion

Our study presents a fresh perspective for future research endeavors and clinical interventions in the field of MM. The identified gene expression patterns and molecular pathways hold immense potential as therapeutic targets for the treatment of multiple myeloma. The utilization of scRNA-seq technology has significantly contributed to a more precise understanding of the complex cellular processes and interactions within MM. Through these advancements, we are now better equipped to unravel the underlying mechanisms driving the development and progression of this complex disease.

Elucidating the role of tumor-associated ALOX5+ mast cells with transformative function in cervical cancer progression via single-cell RNA sequencing

Background

Cervical cancer (CC) is the fourth most common malignancy among women globally and serves as the main cause of cancer-related deaths among women in developing countries. The early symptoms of CC are often not apparent, with diagnoses typically made at advanced stages, which lead to poor clinical prognoses. In recent years, numerous studies have shown that there is a close relationship between mast cells (MCs) and tumor development. However, research on the role MCs played in CC is still very limited at that time. Thus, the study conducted a single-cell multi-omics analysis on human CC cells, aiming to explore the mechanisms by which MCs interact with the tumor microenvironment in CC. The goal was to provide a scientific basis for the prevention, diagnosis, and treatment of CC, with the hope of improving patients' prognoses and quality of life.

Method

The present study acquired single-cell RNA sequencing data from ten CC tumor samples in the ArrayExpress database. Slingshot and AUCcell were utilized to infer and assess the differentiation trajectory and cell plasticity of MCs subpopulations. Differential expression analysis of MCs subpopulations in CC was performed, employing Gene Ontology, gene set enrichment analysis, and gene set variation analysis. CellChat software package was applied to predict cell communication between MCs subpopulations and CC cells. Cellular functional experiments validated the functionality of TNFRSF12A in HeLa and Caski cell lines. Additionally, a risk scoring model was constructed to evaluate the differences in clinical features, prognosis, immune infiltration, immune checkpoint, and functional enrichment across various risk scores. Copy number variation levels were computed using inference of copy number variations.

Result

The obtained 93,524 high-quality cells were classified into ten cell types, including T_NK cells, endothelial cells, fibroblasts, smooth muscle cells, epithelial cells, B cells, plasma cells, MCs, neutrophils, and myeloid cells. Furthermore, a total of 1,392 MCs were subdivided into seven subpopulations: C0 CTSG+ MCs, C1 CALR+ MCs, C2 ALOX5+ MCs, C3 ANXA2+ MCs, C4 MGP+ MCs, C5 IL32+ MCs, and C6 ADGRL4+ MCs. Notably, the C2 subpopulation showed close associations with tumor-related MCs, with Slingshot results indicating that C2 subpopulation resided at the intermediate-to-late stage of differentiation, potentially representing a crucial transition point in the benign-to-malignant transformation of CC. CNVscore and bulk analysis results further confirmed the transforming state of the C2 subpopulation. CellChat analysis revealed TNFRSF12A as a key receptor involved in the actions of C2 ALOX5+ MCs. Moreover, in vitro experiments indicated that downregulating the TNFRSF12A gene may partially inhibit the development of CC. Additionally, a prognosis model and immune infiltration analysis based on the marker genes of the C2 subpopulation provided valuable guidance for patient prognosis and clinical intervention strategies.

Conclusions

We first identified the transformative tumor-associated MCs subpopulation C2 ALOX5+ MCs within CC, which was at a critical stage of tumor differentiation and impacted the progression of CC. In vitro experiments confirmed the inhibitory effect of knocking down the TNFRSF12A gene on the development of CC. The prognostic model constructed based on the C2 ALOX5+MCs subset demonstrated excellent predictive value. These findings offer a fresh perspective for clinical decision-making in CC.

Single-cell RNA sequencing reveals that MYBL2 in malignant epithelial cells is involved in the development and progression of ovarian cancer

Background

Ovarian carcinoma (OC) is a prevalent gynecological malignancy associated with high recurrence rates and mortality, often diagnosed at advanced stages. Despite advances in immunotherapy, immune exhaustion remains a significant challenge in achieving optimal tumor control. However, the exploration of intratumoral heterogeneity of malignant epithelial cells and the ovarian cancer tumor microenvironment is still limited, hindering our comprehensive understanding of the disease.

Materials and methods

Utilizing single-cell RNA sequencing (scRNA-seq), we comprehensively investigated the cellular composition across six ovarian cancer patients with omental metastasis. Our focus centered on analysis of the malignant epithelial cells. Employing CytoTRACE and slingshot pseudotime analyses, we identified critical subpopulations and explored associated transcription factors (TFs) influencing ovarian cancer progression. Furthermore, by integrating clinical factors from a large cohort of bulk RNA sequencing data, we have established a novel prognostic model to investigate the impact of the tumor immune microenvironment on ovarian cancer patients. Furthermore, we have investigated the condition of immunological exhaustion.

Results

Our study identified a distinct and highly proliferative subgroup of malignant epithelial cells, known as C2 TOP2A+ TCs. This subgroup primarily consisted of patients who hadn't received neoadjuvant chemotherapy. Ovarian cancer patients with elevated TOP2A expression exhibited heightened sensitivity to neoadjuvant chemotherapy (NACT). Moreover, the transcription factor MYBL2 in this subgroup played a critical role in ovarian cancer development. Additionally, we developed an independent prognostic indicator, the TOP2A TCs Risk Score (TTRS), which revealed a correlation between the High TTRS Group and unfavorable outcomes. Furthermore, immune infiltration and drug sensitivity analyses demonstrated increased responsiveness to Paclitaxel, Cisplatin, and Gemcitabine in the Low TTRS Group.

Conclusion

This research deepens our understanding of malignant epithelial cells in ovarian cancer and enhances our knowledge of the ovarian cancer immune microenvironment and immune exhaustion. We have revealed the heightened susceptibility of the C2 TOP2A+ TCs subgroup to neoadjuvant chemotherapy and emphasized the role of MYBL2 within the C2 subgroup in promoting the occurrence and progression of ovarian cancer. These insights provide valuable guidance for the management of ovarian cancer treatment.

Predicting prognosis in colorectal cancer patients with curative resection using albumin, lymphocyte count and RAS mutations

Colorectal cancer (CRC) poses a significant global health challenge, demanding reliable prognostic tools to guide treatment decisions. This study introduces a novel prognostic scoring system, the albumin-total lymphocyte count-RAS index (ALRI), integrating serum albumin, lymphocyte count, and RAS gene mutations. A cohort of 445 stage I-III CRC patients undergoing curative resection was analyzed, revealing ALRI's association with clinicopathological factors, including age, tumor location, and invasion depth. The ALRI demonstrated superior prognostic value, with a cutoff value of 2 distinguishing high and low-risk groups. The high-ALRI group exhibited elevated rates of recurrence. Univariate and multivariate analyses identified ALRI as an independent predictor for both 5 year recurrence-free survival (RFS) and overall survival (OS). Kaplan-Meier curves illustrated significant differences in RFS and OS between high and low-ALRI groups, emphasizing ALRI's potential as a prognostic marker. Importantly, ALRI outperformed existing nutritional indices, such as controlling nutritional status and neutrophil-to-lymphocyte ratio, in predicting overall survival. The study underscores the comprehensive insight provided by ALRI, combining inflammatory, nutritional, and genetic information for robust prognostication in CRC patients. This user-friendly tool demonstrates promise for preoperative prognosis and personalized treatment strategies, emphasizing the crucial role of inflammation and nutrition in CRC outcomes.

Macrophage migration inhibitor factor (MIF): Potential role in cognitive impairment disorders

Macrophage migration inhibitory factor (MIF) is a cytokine in the immune system, participated in both innate and adaptive immune responses. Except from immune cells, MIF is also secreted by a variety of non-immune cells, including hematopoietic cells, endothelial cells (ECs), and neurons. MIF plays a crucial role in various diseases, such as sepsis, rheumatoid arthritis, acute kidney injury, and neurodegenerative diseases. The role of MIF in the neuropathogenesis of cognitive impairment disorders is emphasized, as it recruits multiple inflammatory mediators, leading to activating microglia or astrocyte-derived neuroinflammation. Furthermore, it contributes to the cell death of neurons and ECs with the binding of apoptosis-inducing factor (AIF) through parthanatos-associated apoptosis-inducing factor nuclease (PAAN) / MIF pathway. This review comprehensively delves into the relationship between MIF and the neuropathogenesis of cognitive impairment disorders, providing a series of emerging MIF-targeted pharmaceuticals as potential treatments for cognitive impairment disorders.

Integrative Single-Cell Analysis of Cardiomyopathy Identifies Differences in Cell Stemness and Transcriptional Regulatory Networks among Fibroblast Subpopulations

Background

Cardiomyopathy encompasses a broad spectrum of diseases affecting myocardial tissue, characterized clinically by abnormalities in cardiac structure, heart failure, and/or arrhythmias. Clinically heterogeneous, major types include dilated cardiomyopathy (DCM), hypertrophic cardiomyopathy (HCM), restrictive cardiomyopathy (RM), ischemic cardiomyopathy (ICM), among which DCM is more prevalent, while ICM exhibits higher incidence and mortality rates. Myocardial injury during cardiomyopathy progression may lead to myocardial fibrosis. Failure to intervene early and inhibit the process of myocardial fibrosis may culminate in heart failure. Cardiac fibroblasts constitute crucial cellular components determining the extent and quality of myocardial fibrosis, with various subpopulations exerting diverse roles in cardiomyopathy progression. Despite this, understanding of the cellular plasticity and transcriptional regulatory networks of cardiac fibroblasts in cardiomyopathy remains limited. Therefore, in this study, we conducted comprehensive single-cell analysis of cardiac fibroblasts in cardiomyopathy to explore differences in cellular plasticity and transcriptional regulatory networks among fibroblast subpopulations, with the aim of providing as many useful references as possible for the diagnosis, prognosis, and treatment of cardiomyopathy.

Materials and Methods

Cells with mitochondrial gene expression comprising >20% of total expressed genes were excluded. Differential expression genes (DEGs) and stemness genes within cardiac fibroblast subpopulations were subjected to Gene Ontology (GO) analysis of biological processes (BP) and AUCell analysis. Monocle software was employed to analyze the pseudo-temporal trajectory of cardiac fibroblasts in cardiomyopathy. Additionally, the Python package SCENIC was utilized to assess enrichment of transcription factors and activity of regulators within cardiac fibroblast subpopulations in cardiomyopathy.

Results

Following batch effect correction, 179,927 cells were clustered into 32 clusters, designated as T_NK cells, endothelial cells, myeloid cells, fibroblasts, pericytes, SMCs, CMs, proliferating cells, EndoCs, and EPCs. Among them, 8148 fibroblasts were further subdivided into 4 subpopulations, namely C0 THBS4+ Fibroblasts, C1 LINC01133+ Fibroblasts, C2 FGF7+ Fibroblasts, and C3 AGT + Fibroblasts. Results from GO_BP and AUCell analyses suggest that C3 AGT + Fibroblasts may be associated with immune response activation, protein transport, and myocardial contractile function, correlating with disease progression in cardiomyopathy. Transcription factor enrichment analysis indicates that FOS is the most significant TF in C3 AGT + Fibroblasts, also associated with the M1 module, possibly implicated in protein hydrolysis, intracellular DNA replication, and cell proliferation. Moreover, correlation analysis of transcriptional regulatory activity between fibroblast subpopulations reveals a more pronounced heterogeneity within C3 AGT + Fibroblasts in cardiomyopathy.

Conclusion

C3 AGT + Fibroblasts exhibit increased sensitivity towards adverse outcomes in cardiomyopathy, such as myocardial fibrosis and impaired cardiac contractile function, compared to other cardiac fibroblast subpopulations. The differential cellular plasticity and transcriptional regulatory activity between C3 AGT + Fibroblasts and other subgroups offer new perspectives for targeting fibroblast subpopulation activity to treat cardiomyopathy. Additionally, stemness genes EPAS1 and MYC, along with the regulator FOS, may play roles in modulating the biological processes of cardiac fibroblasts in cardiomyopathy.

Identification of necroptosis-related gene signatures for predicting the prognosis of ovarian cancer

Ovarian cancer (OC) is one of the most prevalent and fatal malignant tumors of the female reproductive system. Our research aimed to develop a prognostic model to assist inclinical treatment decision-making.Utilizing data from The Cancer Genome Atlas (TCGA) and copy number variation (CNV) data from the University of California Santa Cruz (UCSC) database, we conducted analyses of differentially expressed genes (DEGs), gene function, and tumor microenvironment (TME) scores in various clusters of OC samples.Next, we classified participants into low-risk and high-risk groups based on the median risk score, thereby dividing both the training group and the entire group accordingly. Overall survival (OS) was significantly reduced in the high-risk group, and two independent prognostic factors were identified: age and risk score. Additionally, three genes-C-X-C Motif Chemokine Ligand 10 (CXCL10), RELB, and Caspase-3 (CASP3)-emerged as potential candidates for an independent prognostic signature with acceptable prognostic value. In Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, pathways related to immune responses and inflammatory cell chemotaxis were identified. Cellular experiments further validated the reliability and precision of our findings. In conclusion, necroptosis-related genes play critical roles in tumor immunity, and our model introduces a novel strategy for predicting the prognosis of OC patients.

XRCC1: a potential prognostic and immunological biomarker in LGG based on systematic pan-cancer analysis

X-ray repair cross-complementation group 1 (XRCC1) is a pivotal contributor to base excision repair, and its dysregulation has been implicated in the oncogenicity of various human malignancies. However, a comprehensive pan-cancer analysis investigating the prognostic value, immunological functions, and epigenetic associations of XRCC1 remains lacking. To address this knowledge gap, we conducted a systematic investigation employing bioinformatics techniques across 33 cancer types. Our analysis encompassed XRCC1 expression levels, prognostic and diagnostic implications, epigenetic profiles, immune and molecular subtypes, Tumor Mutation Burden (TMB), Microsatellite Instability (MSI), immune checkpoints, and immune infiltration, leveraging data from TCGA, GTEx, CELL, Human Protein Atlas, Ualcan, and cBioPortal databases. Notably, XRCC1 displayed both positive and negative correlations with prognosis across different tumors. Epigenetic analysis revealed associations between XRCC1 expression and DNA methylation patterns in 10 cancer types, as well as enhanced phosphorylation. Furthermore, XRCC1 expression demonstrated associations with TMB and MSI in the majority of tumors. Interestingly, XRCC1 gene expression exhibited a negative correlation with immune cell infiltration levels, except for a positive correlation with M1 and M2 macrophages and monocytes in most cancers. Additionally, we observed significant correlations between XRCC1 and immune checkpoint gene expression levels. Lastly, our findings implicated XRCC1 in DNA replication and repair processes, shedding light on the precise mechanisms underlying its oncogenic effects. Overall, our study highlights the potential of XRCC1 as a prognostic and immunological pan-cancer biomarker, thereby offering a novel target for tumor immunotherapy.

Potential Progression Mechanism and Key Genes in Early Stage of mTBI

Chronic traumatic encephalopathy (CTE) is a neurodegenerative disease caused by repetitive mild traumatic brain injury (rmTBI), and the lack of sensitive diagnostic and prognostic biomarkers for rmTBI leads to long-term sequelae after injury. The purpose of this study is to identify key genes of rmTBI and find the potential progression mechanism in early stage of mTBI. We downloaded the gene expression profiles of GSE2871 from Gene Expression Omnibus (GEO) datasets. Differentially expressed genes (DEGs) were screened from the cerebral cortex of rats 24 hours after smTBI, and these DEGs were then subjected to GO enrichment analysis, KEGG pathway analysis, PPI analysis, and hub analysis. Key genes were identified as the most significantly expressed genes and had a higher degree of connectivity from hub genes. By using homemade metal pendulum impact equipment and a multiple regression discriminant equation to assess the severity of rats after injury, smTBI and rmTBI rat models were established in batches, and q-PCR analyses were performed to verify the key genes. The main KEGG pathways were cytokine-cytokine receptor interaction and neuroactive ligand-receptor interaction. SPP1 and C3 were the most significant DEGs, and their connectivity degree was the highest 24 hours after smTBI (logFC > 4; connectivity degree ＞15). The q-PCR analyses were performed 24 hours and 14 days after mTBI. The results showed that SPP1 and C3 were significantly upregulated in smTBI and in rmTBI at 24 hours after injury compared with their levels in sham-injured rats, and the phenomenon persisted 14 days after injury. Notably, 14 days after injury, both of these genes were significantly upregulated in the rmTBI group compared with the smTBI. These pathways and genes identified could help understanding the development in mTBI.

Transcriptional Pathology Evolves over Time in Rat Hippocampus after Lateral Fluid Percussion Traumatic Brain Injury

Traumatic brain injury (TBI) causes acute and lasting impacts on the brain, driving pathology along anatomical, cellular, and behavioral dimensions. Rodent models offer an opportunity to study the temporal progression of disease from injury to recovery. Transcriptomic and epigenomic analysis were applied to evaluate gene expression in ipsilateral hippocampus at 1 and 14 days after sham (n = 2 and 4, respectively per time point) and moderate lateral fluid percussion injury (n = 4 per time point). This enabled the identification of dynamic changes and differential gene expression (differentially expressed genes; DEGs) modules linked to underlying epigenetic response. We observed acute signatures associated with cell death, astrocytosis, and neurotransmission that largely recovered by 2 weeks. Inflammation and immune signatures segregated into upregulated modules with distinct expression trajectories and functions. Whereas most down-regulated genes recovered by 14 days, two modules with delayed and persistent changes were associated with cholesterol metabolism, amyloid beta clearance, and neurodegeneration. Differential expression was paralleled by changes in histone H3 lysine residue 4 trimethylation at the promoters of DEGs at 1 day post-TBI, with the strongest changes observed for inflammation and immune response genes. These results demonstrate how integrated genomics analysis in the pre-clinical setting has the potential to identify stage-specific biomarkers for injury and/or recovery. Though limited in scope here, our general strategy has the potential to capture pathological signatures over time and evaluate treatment efficacy at the systems level.