Skip to main content

Differential expression of the Nrf2-linked genes in pediatric septic shock

Abstract

Introduction

Experimental data from animal models of sepsis support a role for a transcription factor, nuclear erythroid-related factor 2 p45-related factor 2 (Nrf2), as a master regulator of antioxidant and detoxifying genes and intermediary metabolism during stress. Prior analysis of a pediatric septic shock transcriptomic database showed that the Nrf2 response is a top 5 upregulated signaling pathway in early pediatric septic shock.

Methods

We conducted a focused analysis of 267 Nrf2-linked genes using a multicenter, genome-wide expression database of 180 children with septic shock 10 years of age or younger and 53 healthy controls. The analysis involved RNA isolated from whole blood within 24 h of pediatric intensive care unit admission for septic shock and a false discovery rate of 5 %. We compared differentially expressed genes from (1) patients with septic shock and healthy controls and (2) across validated gene expression–based subclasses of pediatric septic shock (endotypes A and B) using several bioinformatic methods.

Results

We found upregulation of 123 Nrf2-linked genes in children with septic shock. The top gene network represented by these genes contained primarily enzymes with oxidoreductase activity involved in cellular lipid metabolism that were highly connected to the peroxisome proliferator activated receptor and the retinoic acid receptor families. Endotype A, which had higher organ failure burden and mortality, exhibited a greater downregulation of Nrf2-linked genes than endotype B, with 92 genes differentially regulated between endotypes.

Conclusions

Our findings indicate that Nrf2-linked genes may contribute to alterations in oxidative signaling and intermediary metabolism in pediatric septic shock.

Introduction

The prevalence of pediatric septic shock is rising [13], and multiorgan dysfunction syndrome (MODS) and secondary infections are leading causes of morbidity and mortality among children admitted to the pediatric intensive care unit (PICU) [4, 5]. Mitochondrial bioenergetic dysfunction and immunoparalysis are thought to play major roles in sepsis-associated deaths [69]. Metabolomic analyses of adults [10, 11] and children [12] suggest that cells fail to generate adequate energy to supply increased metabolic demands during sepsis, and this phenomenon is associated with mortality in adults [10, 11]. Oxidative stress is strongly linked to the mitochondria and is a hallmark of sepsis-associated MODS [1315].

Evidence for mitochondrial dysfunction in pediatric septic shock is supported by a recent report that nuclear-encoded mitochondrial genes are differentially expressed early in pediatric septic shock compared with healthy controls and across septic shock endotypes [16]. Moreover, direct measurements of mitochondrial respiration demonstrated decreased bioenergetic reserve and increased mitochondrial uncoupling in peripheral blood mononuclear cells of children with early sepsis [17]. Mitochondrial uncoupling is directly regulated by reactive oxygen species (ROS) production and glutathionylation status of the uncoupling proteins that regulate ROS production and influence cell signaling [18].

A mechanistic link between sensing of oxidative stress, immune dysregulation, and mitochondrial bioenergetic failure can be made with a nuclear transcription factor, nuclear erythroid-related factor 2 p45-related factor 2 (Nrf2) [19, 20]. Under conditions of increased oxidative stress, Nrf2 facilitates upregulation of genes involved in the antioxidant response [19]. Specifically, Nrf2 promotes a gene expression profile that directs glycolytic intermediates through the pentose phosphate pathway (PPP) to generate reducing equivalents, such as nicotinamide adenine dinucleotide phosphate hydrate (NADPH), used by glutathione (GSH) and thioredoxin to combat oxidative stress. In addition to restoring cellular redox balance, Nrf2 directly regulates cellular energy metabolism by modulating substrate availability for mitochondrial respiration [19, 20]. Prior studies have also shown that loss of Nrf2 leads to mitochondrial depolarization with a resultant decrease in cellular ATP levels, whereas genetic activation of Nrf2 increases the efficiency of mitochondrial oxidative phosphorylation [20]. Finally, in an Nrf2-knockout mouse model, exposure to lipopolysaccharide or cecal ligation and puncture induced higher levels of lung inflammation, proinflammatory gene expression, and mortality compared with wild-type mice [21]. Analysis of an established genome-wide expression database of children with septic shock [22] demonstrated that the Nrf2 oxidative stress response pathway is upregulated in children with septic shock [23].

We sought to further investigate the differential expression of Nrf2-linked genes by performing a focused analysis of a U.S. multisite genome-wide expression database of children with septic shock [22]. In our directed approach, we focused on 267 genes linked to the Nrf2 pathway. The genes were selected on the basis of existing literature [19, 2428]. We hypothesized that expression of whole blood–derived, Nrf2-linked genes would be differentially regulated between pediatric patients with septic shock within the first 24 h of presentation to the PICU and healthy controls. We further hypothesized that Nrf2 pathway genes would be differentially regulated between two endotypes of pediatric septic shock, defined by differential expression of genes involved in adaptive immunity and glucocorticoid receptor signaling and that have previously been shown to have distinct clinical phenotypes (endotypes A and B) [29].

Material and methods

Patients and data collection

The study protocol was approved by the institutional review boards of each participating institution. The name of each institutional review board corresponds to the name of each respective institution: Cincinnati Children’s Hospital Medical Center, The Children’s Hospital of Philadelphia, University of California Benioff Children’s Hospital Oakland, Penn State Hershey Children’s Hospital, Children’s Mercy Hospital, Children’s Hospital of Orange County, Akron Children’s Hospital, Children’s National Medical Center, Miami Children’s Hospital, Texas Children’s Hospital, and C.S. Mott Children’s Hospital at the University of Michigan. Children 10 years of age or younger admitted to the PICU who met pediatrics-specific criteria for septic shock were eligible for enrollment [30]. Age-matched controls were recruited from the ambulatory departments of participating institutions using published inclusion and exclusion criteria [31]. All subjects, including controls, and the data collection methods have been previously reported in microarray-based studies addressing hypotheses entirely different from that of the present study, and details of the study protocol were previously published [32]. All microarray data have been deposited in the National Center for Biotechnology Gene Expression Omnibus database (accession numbers [GEO:GSE26440] and [GEO:GSE26378]).

RNA extraction and microarray hybridization

Written informed consent to participate in this study was obtained from the parents or legal guardians of all children with septic shock and all control subjects. Blood samples were obtained within the first 24 h of meeting criteria for septic shock. Total RNA was isolated from whole blood using the PAXGene Blood RNA System (PreAnalytiX, Hombrechtikon, Switzerland). Microarray hybridization was performed as previously described using the GeneChip Human Genome U133 Plus 2.0 Array (Affymetrix, Santa Clara, CA, USA) [31].

Nrf2-linked gene selection

Nrf2-linked genes chosen for our focused analysis were derived from data compiled in independent studies that used microarray analyses and chromatin immunoprecipitation with massively parallel DNA sequencing (ChIP-Seq) experiments [19, 24–28; and references therein]. We identified 267 genes linked to Nrf2, as defined by Nrf2 target genes containing putative antioxidant response elements or xenobiotic response elements in their promoter regions and genes that directly or indirectly regulate expression of Nrf2. These 267 genes correspond to 566 gene probes on the GeneChip Human Genome U133 Plus 2.0 Array, as shown in Additional file 1: Table S1.

Data analysis

We analyzed existing normalized microarray data. The original analyses were performed using one patient sample per chip. Image files were captured using a GeneChip Scanner 3000 (Affymetrix). Raw data files (.CEL) were subsequently preprocessed using robust multiarray average (RMA) normalization and GeneSpring GX 7.3 software (Agilent Technologies, Palo Alto, CA, USA). All signal intensity–based data were used after RMA normalization, which specifically suppresses all but significant variation among lower-intensity probe sets [33]. All chips representing septic shock samples were then normalized to the respective median values of controls on a per-gene basis.

Differences in mRNA abundance between study groups were determined using Welch’s t test and corrections for multiple comparisons using a Benjamini-Hochberg false discovery rate (FDR) of 5 %. We did not include predetermined fold expression filters in the analysis, because the biological implications of a specific threshold change in our target gene set have not previously been established. Thus, to account for the possibility that even a modest change in the expression of genes from a common metabolic pathway could yield dramatic variability in flux through that pathway [34], we reported all significant differences in gene expression using the 5 % FDR in this analysis. For clarity, further details regarding microarray data analysis and gene list generation are provided in the Results section.

Gene lists of differentially regulated genes were analyzed using the Ingenuity Pathway Analysis (IPA) application (Ingenuity Systems, Redwood City, CA, USA) to explore potential associations with specific domains of Nrf2 function. IPA is a database generated from peer-reviewed scientific publications that provides a tool for discovery of signaling pathways and gene networks within the uploaded gene lists. Adjunct analyses of gene lists were conducted using the ToppGene application [35].

Gene expression mosaics representing the expression patterns of differentially regulated genes were generated using the Gene Expression Dynamics Inspector (GEDI) [36]. GEDI produces expression mosaics recognizable via human pattern recognition. The algorithm for creating the mosaics is a self-organizing map to depict complex genomic data.

Ordinal and continuous clinical variables not normally distributed were evaluated by analysis of variance on ranks. Dichotomous clinical variables were analyzed using a χ2 test (SigmaStat software; Systat Software, San Jose, CA, USA).

Results

Differential regulation of Nrf2-linked genes in patients with septic shock versus healthy controls

We first determined how many gene probes were differentially regulated between children with septic shock (n =180) and healthy control pediatric subjects (n =53). The demographic characteristics of the two study groups were reported previously by Weiss et al. [16]. We conducted a Welch’s t test starting with all 54,675 gene probes on the array and corrected for multiple comparisons using a Benjamini-Hochberg FDR of 5 %. There were 23,912 gene probes differentially regulated between children with septic shock and healthy control pediatric subjects.

We then conducted a Venn diagram analysis to determine how many of the 566 gene probes corresponding to the Nrf2-linked genes were found in this list of 23,912 differentially regulated gene probes. This analysis yielded 281 gene probes corresponding to 178 unique genes, of which 123 were upregulated and 55 were downregulated in children with septic shock relative to controls. The complete list of differentially expressed Nrf2-linked genes is provided in Additional file 2: Table S2.

Uploading the 281 differentially regulated Nrf2-linked gene probes to the IPA platform resulted in the top 5–scoring canonical pathways listed in Table 1. The gene network shown in Fig. 1 demonstrates that the differentially regulated Nrf2-linked genes are enriched for genes encoding enzymes involved in fatty acid metabolism and are connected to the peroxisome proliferator-activated receptor (PPAR) and retinoic acid receptor-α families. The network genes are listed in Additional file 3: Table S3. A majority of these genes are upregulated in septic shock relative to controls. Importing the network genes from IPA to the ToppGene platform returned “oxidoreductase activity” as the top molecular function and “cellular lipid metabolism” as the top biological process.

Table 1 Top 5 canonical pathways in pediatric septic shock vs. controls
Fig. 1
figure 1

Differentially regulated genes corresponding to a gene network having peroxisome proliferator-activated receptor-α (PPARA)-related and retinoic acid receptor-α (RXRA) genes as highly connected nodes. The degree of green intensity in a gene node corresponds to decreased expression, and the degree of red intensity in a given gene node corresponds to increased expression in the subjects with septic shock, relative to controls, respectively. The list of network genes is provided in Additional file 3: Table S3. HMG-CoA 3-hydroxy-3-methylglutaryl-coenzyme A synthase

Differential expression of Nrf2-linked genes across endotypes of pediatric septic shock

We next determined whether the Nrf2-linked genes were differentially regulated across previously validated gene expression-based subclasses of pediatric septic shock: endotype A and endotype B [29]. The clinical and demographic data for the patients in septic shock endotype A (n =60) and endotype B (n =160) have been previously published by Weiss et al. [16]. Patients with endotype A have a higher mortality rate, Pediatric Risk of Mortality score, and pediatric sepsis biomarker risk model–based mortality risk, as well as the maximum number of organ failures, compared with patients in endotype B.

We first determined how many gene probes were differentially regulated between endotype A and endotype B patients. We conducted a Welch’s t test starting with all 54,675 gene probes on the array and corrected for multiple comparisons using a Benjamini-Hochberg FDR of 5 %. There were 11,630 gene probes differentially regulated between endotype A and endotype B patients.

We then conducted a Venn diagram analysis to determine how many of the 566 gene probes corresponding to the Nrf2-linked genes were found in this list of 11,630 differentially regulated gene probes. This analysis yielded 138 gene probes corresponding to 92 unique genes, of which 40 were upregulated and 52 were downregulated in endotype A patients, relative to endotype B patients. The complete list of differentially expressed Nrf2-linked genes is provided in Additional file 4: Table S4.

We used the GEDI program to build gene expression mosaics to provide a global representation of the 92 differentially expressed genes in endotype A and endotype B patients. The GEDI mosaics in Fig. 2 illustrate that endotype A patients exhibited a greater repression of Nrf2-linked genes (higher proportion of blue color intensity) than endotype B patients.

Fig. 2
figure 2

Gene Expression Dynamics Inspector–generated mosaics of differentially expressed mitochondrial genes for the two previously defined septic shock endotypes. The 138 gene probes, corresponding to 92 unique genes, are depicted along the same coordinates across the two expression mosaics. Red intensity correlates with increased gene expression, and blue intensity correlates with decreased gene expression. Clear differences in color patterns illustrate differential expression of Nrf2-regulated genes across patient endotypes A and B, with general downregulation in endotype A. Endotype A subjects have higher illness severity, higher mortality, and higher organ failure burden than endotype B subjects

Among the 92 differentially expressed genes, the top 5 canonical pathways from IPA are shown in Table 2. The top scoring gene network among these differentially regulated genes is shown in Fig. 3. The network genes are listed in Additional file 5: Table S5. Uploading the network genes to the ToppGene platform analysis returned “glutathione binding” as the top molecular function and “response to oxidative stress” as the top biological process for network 1. The second top scoring gene network among these differentially regulated genes is shown in Fig. 4. The network genes are listed in Additional file 6: Table S6. Uploading these network genes to the ToppGene analysis platform returned “superoxide-generating NADPH oxidase activity” as the top molecular function and “superoxide metabolic process” as the top biological process for network 2.

Table 2 Top 5 canonical pathways in pediatric endotype A vs. endotype B
Fig. 3
figure 3

Differentially regulated genes corresponding to a gene network using glutathione to respond to oxidative stress as highly connected nodes. The degree of green intensity in a gene node corresponds to decreased expression, and the degree of red intensity in a given gene node corresponds to increased expression in the endotype A subjects with septic shock, relative to the endotype B subjects with septic shock, respectively. The list of network genes is provided in Additional file 5: Table S5

Fig. 4
figure 4

Differentially regulated genes corresponding to a gene network using superoxide-generating nicotinamide adenine dinucleotide phosphate hydrate (NADPH) oxidase activity to respond to reactive oxygen species as highly connected nodes. The degree of green intensity in a gene node corresponds to decreased expression, and the degree of red intensity in a given gene node corresponds to increased expression in the endotype A subjects with septic shock, relative to the endotype B subjects with septic shock, respectively. The list of network genes is provided in Additional file 6: Table S6. PI3K phosphatidylinositide 3-kinase

Discussion

The Nrf2 pathway has been shown in previous microarray analyses of pediatric sepsis to be in the top 5 signaling and metabolic pathways upregulated in whole blood–derived RNA, with 24 genes represented [23]. The present study delves more deeply into the discovery of Nrf2-linked genes as a critical pathway in pediatric sepsis. We found differential expression of Nrf2-linked genes in children with septic shock compared with controls and that the Nrf2-linked genes showed a greater degree of repression in endotype A subjects, who tend to have the most organ dysfunction and highest mortality.

Nrf2 plays a role in degradation of triglycerides and/or phospholipids and enzymes involved in fatty acid oxidation. In addition, Nrf2 negatively regulates many genes encoding enzymes involved in lipid biosynthesis, fatty acid desaturation, and fatty acid transport through unknown mechanisms [19]. Langley and colleagues performed a metabolomic study of adults with septic shock and demonstrated that fatty acid transport, β-oxidation, gluconeogenesis, and the citric acid cycle pathways are all differentially regulated in non-survivors compared with survivors of septic shock [11]. Specifically, lactate, pyruvate, α-ketoglutarate, oxaloacetate, and acyl-carnitine are all higher in the plasma of sepsis non-survivors than in sepsis survivors [11], purportedly owing to an inability of sepsis non-survivors to harness energetic substrates for aerobic catabolism [11]. Rogers et al also performed metabolomic analysis of critically ill adults using the same cohort as that analyzed by Langley et al and noted that, of the 31 metabolites that differentiated survivors from non-survivors, the six metabolites that were lower in the patients who died were a part of the lipid metabolism pathway [10]. Mickiewicz and colleagues performed nuclear magnetic resonance analysis of metabolites in children with sepsis compared with healthy control children and found an increase in three compounds (2-hydroxybutyrate, 2-hydroxyisovalerate, and lactate) that are associated with enhanced fat breakdown resulting in ketoacid and lactic acid production owing to increased energy demands during septic shock [12]. Our data are consistent with these previous studies and, overall, support the concept that pediatric septic shock is characterized by alterations in the expression of genes essential for oxidative stress responses and lipid metabolism.

Concurrent measures of oxidative stress and metabolite data are not available to determine how the observed Nrf2 pathway gene expression changes might alter the total antioxidant status, fatty acid metabolism, and ATP production in early pediatric sepsis. It is likely that changes in intermediary metabolism and oxidative state affect Nrf2-linked gene expression, and researchers in future studies should measure the effects of oxidative stress, metabolite concentration, and flux through various metabolic pathways influenced by Nrf2, such as glycolysis, the citric acid cycle, the PPP, fatty acid oxidation, and fatty acid synthesis.

The present study is constrained by the fact that we limited the analysis to 267 genes linked to Nrf2, as determined by previous human and murine studies. Linkage was defined as either genes regulated by Nrf2 or genes that regulate Nrf2 function. It is certainly possible that our selection method excluded some relevant genes. Limiting our working gene list in this manner is a potential source of bias. We attempted to mitigate this bias by conducting the initial statistical analyses, with corrections for multiple comparisons, using all gene probes available on the array. We subsequently searched the resulting gene lists for our selected Nrf2 genes and found that many of these genes were differentially regulated when we considered all available gene probes.

Conclusions

Nrf2-linked genes involved in cellular lipid metabolism were predominantly upregulated early in pediatric septic shock compared with healthy controls. Genes corresponding to the oxidative stress response mediated by GSH-binding enzymes and NADPH-oxidase generation of superoxide were predominantly downregulated in a distinct, previously defined endotype of pediatric septic shock having higher organ failure burden and mortality. Our findings indicate that Nrf2-linked genes may contribute to alterations in oxidative signaling and intermediary metabolism in pediatric septic shock. Studying the temporal changes in expression of Nrf2 pathway genes and correlation of gene expression with metabolomic analysis of intermediary metabolites should provide insight into the energy and oxidative states of children over the time course of sepsis. We propose that future work should be focused on the hypothesis that differential regulation of Nrf2-linked genes may be an important mechanism contributing to the host ability to combat oxidative stress and reprogram metabolic pathways to provide necessary energy to survive childhood sepsis and MODS.

Key messages

  • The oxidative stress transcription factor Nrf2-linked genes were differentially expressed in children with septic shock compared with healthy controls.

  • Nrf2-linked genes can be used to differentiate two previously validated subclasses of children with septic shock, with a greater degree of repression in children with a higher mortality and more organ failure.

  • Nrf2 may be an important transcription factor contributing to alterations in oxidative signaling and intermediary metabolism in pediatric sepsis.

Abbreviations

ChIP-Seq:

Chromatin immunoprecipitation with massively parallel DNA sequencing

FDR:

False discovery rate

GEDI:

Gene Expression Dynamics Inspector

GSH:

Glutathione

IL:

Interleukin

IPA:

Ingenuity Pathway Analysis

LPS:

Lipopolysaccharide

MODS:

Multiorgan dysfunction syndrome

NADPH:

Nicotinamide adenine dinucleotide phosphate hydrate

Nrf2:

Nuclear erythroid-related factor 2 p45-related factor 2

PI3K:

Phosphatidylinositide 3-kinase

PICU:

Pediatric intensive care unit

PPAR:

Peroxisome proliferator-activated receptor

PPP:

Pentose phosphate pathway

RMA:

Robust multiarray average

ROS:

Reactive oxygen species

RXR:

Retinoic acid receptor

References

  1. Ruth A, McCracken CE, Fortenberry JD, Hall M, Simon HK, Hebbar KB. Pediatric severe sepsis: current trends and outcomes from the Pediatric Health Information Systems database. Pediatr Crit Care Med. 2014;15:828–38.

    Article  Google Scholar 

  2. Hartman ME, Linde-Zwirble WT, Angus DC, Watson RS. Trends in the epidemiology of pediatric severe sepsis. Pediatr Crit Care Med. 2013;14:686–93.

    Article  Google Scholar 

  3. Balamuth F, Weiss SL, Neuman MI, Scott H, Brady PW, Paul R, et al. Pediatric severe sepsis in U.S. children’s hospitals. Pediatr Crit Care Med. 2014;15:798–805.

    Article  Google Scholar 

  4. Hall MW, Knatz NL, Vetterly C, Tomarello S, Wewers MD, Volk HD, et al. Immunoparalysis and nosocomial infection in children with multiple organ dysfunction syndrome. Intensive Care Med. 2011;37:525–32.

    Article  CAS  Google Scholar 

  5. Proulx F, Joyal JS, Mariscalco MM, Leteurtre S, Leclerc F, Lacroix J. The pediatric multiple organ dysfunction syndrome. Pediatr Crit Care Med. 2009;10:12–22.

    Article  Google Scholar 

  6. Crouser ED. Mitochondrial dysfunction in septic shock and multiple organ dysfunction syndrome. Mitochondrion. 2004;4:729–41.

    Article  CAS  Google Scholar 

  7. Fink MP. Bench-to-bedside review: cytopathic hypoxia. Crit Care. 2002;6:491–9.

    Article  Google Scholar 

  8. Singer M. Mitochondrial function in sepsis: acute phase versus multiple organ failure. Crit Care Med. 2007;35(9 Suppl):S441–8.

    Article  CAS  Google Scholar 

  9. Singer M, De Santis V, Vitale D, Jeffcoate W. Multiorgan failure is an adaptive, endocrine-mediated, metabolic response to overwhelming systemic inflammation. Lancet. 2004;364:545–8.

    Article  Google Scholar 

  10. Rogers AJ, McGeachie M, Baron RM, Gazourian L, Haspel JA, Nakahira K, et al. Metabolomic derangements are associated with mortality in critically ill adult patients. PLoS One. 2014;9:e87538.

    Article  Google Scholar 

  11. Langley RJ, Tsalik EL, van Velkinburgh JC, Glickman SW, Rice BJ, Wang C, et al. An integrated clinico-metabolomic model improves prediction of death in sepsis. Sci Transl Med. 2013;5:195ra195.

    Article  Google Scholar 

  12. Mickiewicz B, Vogel HJ, Wong HR, Winston BW. Metabolomics as a novel approach for early diagnosis of pediatric septic shock and its mortality. Am J Respir Crit Care Med. 2013;187:967–76.

    Article  CAS  Google Scholar 

  13. Trentadue R, Fiore F, Massaro F, Papa F, Iuso A, Scacco S, et al. Induction of mitochondrial dysfunction and oxidative stress in human fibroblast cultures exposed to serum from septic patients. Life Sci. 2012;91:237–43.

    Article  CAS  Google Scholar 

  14. Ware LB, Fessel JP, May AK, Roberts 2nd LJ. Plasma biomarkers of oxidant stress and development of organ failure in severe sepsis. Shock. 2011;36:12–7.

    Article  CAS  Google Scholar 

  15. Victor VM, Rocha M, Esplugues JV, De la Fuente M. Role of free radicals in sepsis: antioxidant therapy. Curr Pharm Des. 2005;11:3141–58.

    Article  CAS  Google Scholar 

  16. Weiss SL, Cvijanovich NZ, Allen GL, Thomas NJ, Freishtat RJ, Anas N, et al. Differential expression of the nuclear-encoded mitochondrial transcriptome in pediatric septic shock. Crit Care. 2014;18:623.

    Article  Google Scholar 

  17. Weiss SL, Selak MA, Tuluc F, Perales Villarroel J, Nadkarni VM, Deutschman CS, et al. Mitochondrial dysfunction in peripheral blood mononuclear cells in pediatric septic shock. Pediatr Crit Care Med. 2015;16:e4–e12.

    Article  Google Scholar 

  18. Mailloux RJ, Harper ME. Uncoupling proteins and the control of mitochondrial reactive oxygen species production. Free Radic Biol Med. 2011;51:1106–15.

    Article  CAS  Google Scholar 

  19. Hayes JD, Dinkova-Kostova AT. The Nrf2 regulatory network provides an interface between redox and intermediary metabolism. Trends Biochem Sci. 2014;39:199–218.

    Article  CAS  Google Scholar 

  20. Holmström KM, Baird L, Zhang Y, Hargreaves I, Chalasani A, Land JM, et al. Nrf2 impacts cellular bioenergetics by controlling substrate availability for mitochondrial respiration. Biol Open. 2013;2:761–70.

    Article  Google Scholar 

  21. Thimmulappa RK, Lee H, Rangasamy T, Reddy SP, Yamamoto M, Kensler TW, et al. Nrf2 is a critical regulator of the innate immune response and survival during experimental sepsis. J Clin Invest. 2006;116:984–95.

    Article  CAS  Google Scholar 

  22. Wong HR. Genome-wide expression profiling in pediatric septic shock. Pediatr Res. 2013;73:564–9.

    Article  CAS  Google Scholar 

  23. Wong HR, Freishtat RJ, Monaco M, Odoms K, Shanley TP. Leukocyte subset-derived genomewide expression profiles in pediatric septic shock. Pediatr Crit Care Med. 2010;11:349–55.

    PubMed  PubMed Central  Google Scholar 

  24. Malhotra D, Portales-Casamar E, Singh A, Srivastava S, Arenillas D, Happel C, et al. Global mapping of binding sites for Nrf2 identifies novel targets in cell survival response through ChIP-Seq profiling and network analysis. Nucleic Acids Res. 2010;38:5718–34.

    Article  CAS  Google Scholar 

  25. Taylor RC, Acquaah-Mensah G, Singhal M, Malhotra D, Biswal S. Network inference algorithms elucidate Nrf2 regulation of mouse lung oxidative stress. PLoS Comput Biol. 2008;4:e1000166.

    Article  Google Scholar 

  26. Agyeman AS, Chaerkady R, Shaw PG, Davidson NE, Visvanathan K, Pandey A, et al. Transcriptomic and proteomic profiling of KEAP1 disrupted and sulforaphane-treated human breast epithelial cells reveals common expression profiles. Breast Cancer Res Treat. 2012;132:175–87.

    Article  CAS  Google Scholar 

  27. Cho HY, Kleeberger SR. Nrf2 protects against airway disorders. Toxicol Appl Pharmacol. 2010;244:43–56.

    Article  CAS  Google Scholar 

  28. Chorley BN, Campbell MR, Wang X, Karaca M, Sambandan D, Bangura F, et al. Identification of novel NRF2-regulated genes by ChIP-Seq: influence on retinoid X receptor α. Nucleic Acids Res. 2012;40:7416–29.

    Article  CAS  Google Scholar 

  29. Wong HR, Cvijanovich NZ, Anas N, Allen GL, Thomas NJ, Bigham MT, et al. Developing a clinically feasible personalized medicine approach to pediatric septic shock. Am J Respir Crit Care Med. 2015;191:309–15.

    Article  Google Scholar 

  30. Goldstein B, Giroir B. Randolph A; International Consensus Conference on Pediatric Sepsis. International Pediatric Sepsis Consensus Conference: definitions for sepsis and organ dysfunction in pediatrics. Pediatr Crit Care Med. 2005;6:2–8.

    Article  Google Scholar 

  31. Wong HR, Shanley TP, Sakthivel B, Cvijanovich N, Lin R, Allen GL, et al. Genome-level expression profiles in pediatric septic shock indicate a role for altered zinc homeostasis in poor outcome. Physiol Genomics. 2007;30:146–55.

    Article  CAS  Google Scholar 

  32. Wong HR, Cvijanovich N, Lin R, Allen GL, Thomas NJ, Willson DF, et al. Identification of pediatric septic shock subclasses based on genome-wide expression profiling. BMC Med. 2009;7:34.

    Article  Google Scholar 

  33. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  Google Scholar 

  34. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.

    Article  CAS  Google Scholar 

  35. Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009;37(Web Server issue):W305–311.

    Article  CAS  Google Scholar 

  36. Eichler GS, Huang S, Ingber DE. Gene Expression Dynamics Inspector (GEDI): for integrative analysis of expression profiles. Bioinformatics. 2003;19:2321–2.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

Financial support was provided by National Institutes of Health (NIH) grants R01 GM108025 and R01 GM099773 (to HRW). JRG is supported by NIH grant T32GM095442. SLW is supported by National Institute of General Medical Sciences grant K23GM110496. The NIH had no role in the design, collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hector R. Wong.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JRG conceived and developed the study and wrote the manuscript. HRW obtained funding for the study, conducted analyses, and edited the manuscript. SLW, NZC, GLA, NJT, RJF, NA, KM, PAC, TPS, MTP, and JF enrolled subjects at the participating institutions and provided clinical data and biological samples. KHo and EF maintained the clinical database and coordinated all interinstitutional research activity. KHa maintained the biological repository and processed all biological samples. All authors read, edited, and approved the manuscript.

Additional files

Additional file 1: Table S1.

List of Nrf2-linked genes (n =267) used in the focused analysis. The 267 genes correspond to 566 gene probes on the Human Genome U133 Plus 2.0 GeneChip. The functional link is provided for each gene. (XLSX 42 kb)

Additional file 2: Table S2.

List of Nrf2-linked gene probes (n =281) corresponding to 178 unique genes differentially regulated between subjects with septic shock and healthy controls. (XLSX 35 kb)

Additional file 3: Table S3.

List of network genes shown in Fig. 1. (XLSX 13 kb)

Additional file 4: Table S4.

List of Nrf2-linked gene probes (n =138) corresponding to 92 unique genes differentially regulated between septic shock endotypes A and B subjects. (XLSX 20 kb)

Additional file 5: Table S5.

List of network genes shown in Fig. 3. (XLSX 12 kb)

Additional file 6: Table S6.

List of network genes shown in Fig. 4. (XLSX 10 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Grunwell, J.R., Weiss, S.L., Cvijanovich, N.Z. et al. Differential expression of the Nrf2-linked genes in pediatric septic shock. Crit Care 19, 327 (2015). https://doi.org/10.1186/s13054-015-1052-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-015-1052-0

Keywords