Screening of Key Genes in Severe Burn Injury at Different Stages via Analyzing Gene Expression Data

Identification of Important Modules and Biomarkers That Are Related to Immune Infiltration Cells in Severe Burns Based on Weighted Gene Co-Expression Network Analysis

Background: Immunosuppression is an important trigger for infection and a significant cause of death in patients with severe burns. Nevertheless, the prognostic value of immune-related genes remains unclear. This study aimed to identify the biomarkers related to immunosuppression in severe burns. Methods: The gene expression profile and clinical data of 185 burn and 75 healthy samples were obtained from the GEO database. Immune infiltration analysis and gene set variation analysis were utilized to identify the disorder of circulating immune cells. A weighted gene co-expression network analysis (WGCNA) was carried out to select immune-related gene modules. Enrichment analysis and protein-protein interaction (PPI) network were performed to select hub genes. Next, LASSO and logistic regression were utilized to construct the hazard regression model with a survival state. Finally, we investigated the correlation between high- and low-risk patients in total burn surface area (TBSA), age, and inhalation injury. Results: Gene set variation analysis (GSVA) and immune infiltration analysis showed that neutrophils increased and T cells decreased in severe burns. In WGCNA, four modular differently expressed in burns and controls were related to immune cells. Based on PPI and enrichment analysis, 210 immune-related genes were identified, mainly involved in T-cell inhibition and neutrophil activation. In LASSO and logistic regression, we screened out key genes, including LCK, SKAP1 and GZMB, and LY9. In the ROC analysis, the area under the curve (AUC) of key genes was 0.945, indicating that the key genes had excellent diagnostic value. Finally, we discovered that the key genes were related to T cells, and the regression model performed well when accompanied by TBSA and age. Conclusion: We identified LCK, SKAP1, GZMB, and LY9 as good prognostic biomarkers that may play a role in post-burn immunosuppression against T-cell dysfunction and as potential immunotherapeutic targets for transformed T-cell dysfunction.

Identifying changes in immune cells and constructing prognostic models using immune-related genes in post-burn immunosuppression

BACKGROUND

Burn patients are prone to infection as well as immunosuppression, which is a significant cause of death. Currently, there is a lack of prognostic biomarkers for immunosuppression in burn patients. This study was conducted to identify immune-related genes that are prognosis biomarkers in post-burn immunosuppression and potential targets for immunotherapy.

METHODS

We downloaded the gene expression profiles and clinical data of 213 burn patients and 79 healthy samples from the Gene Expression Omnibus (GEO) database. Immune infiltration analysis was used to identify the proportion of circulating immune cells. Functional enrichment analyses were carried out to identify immune-related genes that were used to build miRNA-mRNA networks to screen key genes. Next, we carried out correlation analysis between immune cells and key genes that were then used to construct logistic regression models in GSE77791 and were validated in GSE19743. Finally, we determined the expression of key genes in burn patients using quantitative reverse transcription polymerase chain reaction (qRT-PCR).

RESULTS

A total of 745 differently expressed genes were screened out: 299 were up-regulated and 446 were down-regulated. The number of Th-cells (CD4+) decreased while neutrophils increased in burn patients. The enrichment analysis showed that down-regulated genes were enriched in the T-cell activation pathway, while up-regulated genes were enriched in neutrophil activation response in burn patients. We screened out key genes (NFATC2, RORA, and CAMK4) that could be regulated by miRNA. The expression of key genes was related to the proportion of Th-cells (CD4+) and survival, and was an excellent predictor of prognosis in burns with an area under the curve (AUC) value of 0.945. Finally, we determined that NFATC2, RORA, and CAMK4 were down-regulated in burn patients.

CONCLUSION

We found that NFATC2, RORA, and CAMK4 were likely prognostic biomarkers in post-burn immunosuppression and potential immunotherapeutic targets to convert Th-cell dysfunction.

The Genetic Evidence of Burn-Induced Cardiac Mitochondrial Metabolism Dysfunction

Burn-induced cardiac dysfunction is thought to involve mitochondrial dysfunction, although the mechanisms responsible are unclear. In this study, we used our established model of in vivo burn injury to understand the genetic evidence of burn-induced mitochondrial confusion dysfunction by describing cardiac mitochondrial metabolism-related gene expression after burn. Cardiac tissue was collected at 24 hours after burn injury. An O2K respirometer system was utilized to measure the cardiac mitochondrial function. Oxidative phosphorylation complex activities were determined using enzyme activity assays. RT Profiler PCR array was used to identify the differential regulation of genes involved in mitochondrial biogenesis and metabolism. The quantitative qPCR and Western blotting were applied to validate the differentially expressed genes. Burn-induced cardiac mitochondrial dysfunction was supported by the finding of decreased state 3 respiration, decreased mitochondrial electron transport chain activity in complex I, III, IV, and V, and decreased mitochondrial DNA-encoded gene expression as well as decreased levels of the corresponding proteins after burn injury. Eighty-four mitochondrial metabolism-related gene profiles were measured. The mitochondrial gene profile showed that 29 genes related to mitochondrial energy and metabolism was differentially expressed. Of these 29 genes, 16 were more than 2-fold upregulated and 13 were more than 2-fold downregulated. All genes were validated using qPCR and partial genes were correlated with their protein levels. This study provides preliminary evidence that a large percentage of mitochondrial metabolism-related genes in cardiomyocytes were significantly affected by burn injury.

Transcriptome modulation by hydrocortisone in severe burn shock: ancillary analysis of a prospective randomized trial

Background

Despite shortening vasopressor use in shock, hydrocortisone administration remains controversial, with potential harm to the immune system. Few studies have assessed the impact of hydrocortisone on the transcriptional response in shock, and we are lacking data on burn shock. Our objective was to assess the hydrocortisone-induced transcriptional modulation in severe burn shock, particularly modulation of the immune response.

Methods

We collected whole blood samples during a randomized controlled trial assessing the efficacy of hydrocortisone administration in burn shock. Using whole genome microarrays, we first compared burn patients (n = 32) from the placebo group to healthy volunteers to describe the transcriptional modulation induced by burn shock over the first week. Then we compared burn patients randomized for either hydrocortisone administration or placebo, to assess hydrocortisone-induced modulation.

Results

Study groups were similar in terms of severity and major outcomes, but shock duration was significantly reduced in the hydrocortisone group. Many genes (n = 1687) were differentially expressed between burn patients and healthy volunteers, with 85% of them exhibiting a profound and persistent modulation over seven days. Interestingly, we showed that hydrocortisone enhanced the shock-associated repression of adaptive, but also innate immunity.

Conclusions

We found that the initial host response to burn shock encompasses wide and persistent modulation of gene expression, with profound modulation of pathways associated with metabolism and immunity. Importantly, hydrocortisone administration may worsen the immunosuppression associated with severe injury. These data should be taken into account in the risk ratio of hydrocortisone administration in patients with inflammatory shock.

Trial registration

ClinicalTrials.gov, NCT00149123. Registered on 6 September 2005.

Electronic supplementary material

The online version of this article (doi:10.1186/s13054-017-1743-9) contains supplementary material, which is available to authorized users.