Transcriptomic profiles of multiple organ dysfunction syndrome phenotypes in pediatric critical influenza.

Background: Influenza virus is responsible for a large global burden of disease, especially in children. Multiple Organ Dysfunction Syndrome (MODS) is a life-threatening and fatal complication of severe influenza infection. Methods: We measured RNA expression of 469 biologically plausible candidate genes in children admitted to North American pediatric intensive care units with severe influenza virus infection with and without MODS. Whole blood samples from 191 influenza-infected children (median age 6.4 years, IQR: 2.2, 11) were collected a median of 27 hours following admission; for 45 children a second blood sample was collected approximately seven days later. Extracted RNA was hybridized to NanoString mRNA probes, counts normalized, and analyzed using linear models controlling for age and bacterial co-infections (FDR q<0.05). Results: Comparing pediatric samples collected near admission, children with Prolonged MODS for ≥7 days (n=38; 9 deaths) had significant upregulation of nine mRNA transcripts associated with neutrophil degranulation (RETN, TCN1, OLFM4, MMP8, LCN2, BPI, LTF, S100A12, GUSB) compared to those who recovered more rapidly from MODS (n=27). These neutrophil transcripts present in early samples predicted Prolonged MODS or death when compared to patients who recovered, however in paired longitudinal samples, they were not differentially expressed over time. Instead, five genes involved in protein metabolism and/or adaptive immunity signaling pathways (RPL3, MRPL3, HLA-DMB, EEF1G, CD8A) were associated with MODS recovery within a week. Conclusion: Thus, early increased expression of neutrophil degranulation genes indicated worse clinical outcomes in children with influenza infection, consistent with reports in adult cohorts with influenza, sepsis, and acute respiratory distress syndrome.

A robust host-response-based signature distinguishes bacterial and viral infections across diverse global populations.

Limited sensitivity and specificity of current diagnostics lead to the erroneous prescription of antibiotics. Host-response-based diagnostics could address these challenges. However, using 4,200 samples across 69 blood transcriptome datasets from 20 countries from patients with bacterial or viral infections representing a broad spectrum of biological, clinical, and technical heterogeneity, we show current host-response-based gene signatures have lower accuracy to distinguish intracellular bacterial infections from viral infections than extracellular bacterial infections. Using these 69 datasets, we identify an 8-gene signature to distinguish intracellular or extracellular bacterial infections from viral infections with an area under the receiver operating characteristic curve (AUROC) > 0.91 (85.9% specificity and 90.2% sensitivity). In prospective cohorts from Nepal and Laos, the 8-gene classifier distinguished bacterial infections from viral infections with an AUROC of 0.94 (87.9% specificity and 91% sensitivity). The 8-gene signature meets the target product profile proposed by the World Health Organization and others for distinguishing bacterial and viral infections.

Profiling chromatin accessibility responses in human neutrophils with sensitive pathogen detection.

Sepsis, sequela of bloodstream infections and dysregulated host responses, is a leading cause of death globally. Neutrophils tightly regulate responses to pathogens to prevent organ damage. Profiling early host epigenetic responses in neutrophils may aid in disease recognition. We performed assay for transposase-accessible chromatin (ATAC)-seq of human neutrophils challenged with six toll-like receptor ligands and two organisms; and RNA-seq after Escherichia coli exposure for 1 and 4 h along with ATAC-seq. ATAC-seq of neutrophils facilitates detection of pathogen DNA. In addition, despite similarities in genomic distribution of differential chromatin changes across challenges, only a fraction overlaps between the challenges. Ligands depict shared signatures, but majority are unique in position, function, and challenge. Epigenomic changes are plastic, only ∼120 are shared by E coli challenges over time, resulting in varied differential genes and associated processes. We identify three classes of gene regulation, chromatin access changes in the promoter; changes in the promoter and distal enhancers; and controlling expression through changes solely in distal enhancers. These and transcription factor footprinting reveal timely and challenge specific mechanisms of transcriptional regulation in neutrophils.

Transcriptomic similarities and differences in host response between SARS-CoV-2 and other viral infections.

The pandemic 2019 novel coronavirus disease (COVID-19) shares certain clinical characteristics with other acute viral infections. We studied the whole-blood transcriptomic host response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using RNAseq from 24 healthy controls and 62 prospectively enrolled patients with COVID-19. We then compared these data to non-COVID-19 viral infections, curated from 23 independent studies profiling 1,855 blood samples covering six viruses (influenza, respiratory syncytial virus (RSV), human rhinovirus (HRV), severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1), Ebola, dengue). We show gene expression changes in COVID-19 versus non-COVID-19 viral infections are highly correlated (r = 0.74, p < 0.001). However, we also found 416 genes specific to COVID-19. Inspection of top genes revealed dynamic immune evasion and counter host responses specific to COVID-19. Statistical deconvolution of cell proportions maps many cell type proportions concordantly shifting. Discordantly increased in COVID-19 were CD56bright natural killer cells and M2 macrophages. The concordant and discordant responses mapped out here provide a window to explore the pathophysiology of the host response to SARS-CoV-2.

The IL20 Genetic Polymorphism Is Associated with Altered Clinical Outcome in Septic Shock.

BACKGROUND: The IL10 family of genes includes crucial immune regulators. We tested the hypothesis that single nucleotide polymorphisms (SNPs) in IL10, IL19, IL20, and IL24 of the IL10 family gene cluster alter the clinical outcome of septic shock. METHODS: Patients with septic shock (n = 1,193) were genotyped for 13 tag SNPs of IL10, IL19, IL20, and IL24. IL20 gene expression was measured in genotyped lymphoblastoid cells in vitro. Cardiac surgical ICU patients (n = 981) were genotyped for IL20 rs2981573 A/G. The primary outcome variable was 28-day mortality. RESULTS: Patients with the G allele of IL20 rs2981573 had a significantly increased hazard of death over the 28-day period compared to patients with the A allele in the septic shock cohort (adjusted hazard ratio 1.27; 95% confidence interval 1.10-1.47; p = 8.0 × 10-4). Patients with the GG genotype had more organ dysfunction (p < 0.05). The GG genotype was associated with increased IL20 gene expression in stimulated lymphoblastoid cells in vitro (p < 0.05). The cardiac surgical ICU patients with the GG genotype had an increased length of ICU stay (p = 0.032). CONCLUSIONS: The GG genotype of IL20 rs2981573 SNP was associated with increased IL20 gene expression and increased adverse outcomes in patients with septic shock and following cardiac surgery.

TNFAIP2 Inhibits Early TNFα-Induced NF-x03BA;B Signaling and Decreases Survival in Septic Shock Patients.

During septic shock, tumor necrosis factor alpha (TNFα) is an early response gene and induces a plethora of genes and signaling pathways. To identify robust signals in genes reliably upregulated by TNFα, we first measured microarray gene expression in vitro and searched methodologically comparable, publicly available data sets to identify concordant signals. Using tag single-nucleotide polymorphisms in the genes common to all data sets, we identified a genetic variant of the TNFAIP2 gene, rs8126, associated with decreased 28-day survival and increased organ dysfunction in an adult cohort in the Vasopressin and Septic Shock Trial. Similar to this cohort, we found that an association with rs8126 and increased organ dysfunction is replicated in a second cohort of septic shock patients in the St. Paul's Hospital Intensive Care Unit. We found that TNFAIP2 inhibits NF-x03BA;B activity, impacting the downstream cytokine interleukin (IL)-8. The minor G allele of TNFAIP2 rs8126 resulted in greater TNFAIP2 expression, decreased IL-8 production and was associated with decreased survival in patients experiencing septic shock. These data suggest that TNFAIP2 is a novel inhibitor of NF-x03BA;B that acts as an autoinhibitor of the TNFα response during septic shock.

VPS13D Gene Variant Is Associated with Altered IL-6 Production and Mortality in Septic Shock.

BACKGROUND: Genetic variations contribute to septic shock mortality. To discover a novel locus, we performed in vitro genome-wide association studies (GWAS) and further tested the result in a cohort of septic shock patients. METHODS: Two in vitro GWAS using a quantitative trait locus analysis of stimulated IL-6 production in lymphoblastoid cells from 60 individuals of European ancestry were performed. VPS13D rs6685273 was genotyped in European ancestry patients (n = 498). The VPS13D gene was silenced in vitro. RESULTS: Two GWAS using lymphoblastoid cells identified the locus of VPS13D rs6685273 that was significant in the same direction in both GWAS. The VPS13D rs6685273 C allele was associated with increased IL-6 production. Patients with septic shock who had the VPS13D rs6685273 CC genotype had an increased 28-day mortality (p = 0.023) and more organ failure (p < 0.05) compared to the CT/TT genotypes. VPS13D in vitro gene silencing in the HeLa cell line increased IL-6 production. Furthermore, the rs6685273 genotype was associated with differential VPS13D splice variant expression. CONCLUSIONS: The VPS13D rs6685273 C allele was associated with increased IL-6 production in vitro. The patients with the VPS13D rs6685273 CC genotype had increased 28-day mortality and increased organ failure. VPS13D appears to regulate IL-6 production.

Identification of a nonsynonymous polymorphism in the SVEP1 gene associated with altered clinical outcomes in septic shock.

OBJECTIVES: Mortality from septic shock is highly heritable. The identification of causal genetic factors is insufficient. To discover key contributors, we first identified nonsynonymous single-nucleotide polymorphisms in conserved genomic regions that are predicted to have significant effects on protein function. We then test the hypothesis that these nonsynonymous single-nucleotide polymorphisms across the genome alter clinical outcome of septic shock. DESIGN: Genetic-association study plus in vitro experiment using primary cells plus in silico analysis using genomic DNA and protein database. SETTING: Twenty-seven ICUs at academic teaching centers in Canada, Australia, and the United States. PATIENTS: Patients with septic shock of European ancestry (n = 520). INTERVENTIONS: Patients with septic shock were genotyped for 843 nonsynonymous single-nucleotide polymorphisms in conserved regions of the genome and are predicted to have damaging effects from the protein sequence. MEASUREMENTS AND MAIN RESULTS: The primary outcome variable was 28-day mortality. Secondary outcome variables were organ dysfunction. Productions of adhesion molecules including interleukin-8, growth-regulated oncogene-α, monocyte chemoattractant protein-1, and monocyte chemoattractant protein-3 were measured in human umbilical vein endothelial cells after SVEP1 gene silencing by RNA interference. Patients with septic shock having the SVEP1 C allele of nonsynonymous single-nucleotide polymorphism, SVEP1 c.2080A>C (p. Gln581His, rs10817033), had a significant increase in the hazard of death over the 28 days (hazard ratio, 1.72; 95% CI, 1.31-2.26; p = 9.7 × 10-5) and increased organ dysfunction and needed more organ support (p < 0.05). Silencing SVEP1 significantly increased interleukin-8, growth-regulated oncogene-α, monocyte chemoattractant protein-1, monocyte chemoattractant protein-3 production in human umbilical vein endothelial cells under lipopolysaccharide stimulation (p < 0.01). CONCLUSIONS: C allele of SVEP1 c.2080A>C (p. Gln581His) single-nucleotide polymorphism, a non-synonymous single-nucleotide polymorphism in conserved regions and predicted to have damaging effects on protein structure, was associated with increased 28-day mortality and organ dysfunction of septic shock. SVEP1 appears to regulate molecules of the leukocyte adhesion pathway.

PCSK9 is a critical regulator of the innate immune response and septic shock outcome.

A decrease in the activity of proprotein convertase subtilisin/kexin type 9 (PCSK9) increases the amount of low-density lipoprotein (LDL) receptors on liver cells and, therefore, LDL clearance. The clearance of lipids from pathogens is related to endogenous lipid clearance; thus, PCSK9 may also regulate removal of pathogen lipids such as lipopolysaccharide (LPS). Compared to controls, Pcsk9 knockout mice displayed decreases in inflammatory cytokine production and in other physiological responses to LPS. In human liver cells, PCSK9 inhibited LPS uptake, a necessary step in systemic clearance and detoxification. Pharmacological inhibition of PCSK9 improved survival and inflammation in murine polymicrobial peritonitis. Human PCSK9 loss-of-function genetic variants were associated with improved survival in septic shock patients and a decrease in inflammatory cytokine response both in septic shock patients and in healthy volunteers after LPS administration. The PCSK9 effect was abrogated in LDL receptor (LDLR) knockout mice and in humans who are homozygous for an LDLR variant that is resistant to PCSK9. Together, our results show that reduced PCSK9 function is associated with increased pathogen lipid clearance via the LDLR, a decreased inflammatory response, and improved septic shock outcome.

A multibiomarker-based outcome risk stratification model for adult septic shock*.

OBJECTIVES: Clinical trials in septic shock continue to fail due, in part, to inequitable and sometimes unknown distribution of baseline mortality risk between study arms. Investigators advocate that interventional trials in septic shock require effective outcome risk stratification. We derived and tested a multibiomarker-based approach to estimate mortality risk in adults with septic shock. DESIGN: Previous genome-wide expression studies identified 12 plasma proteins as candidates for biomarker-based risk stratification. The current analysis used banked plasma samples and clinical data from existing studies. Biomarkers were assayed in plasma samples obtained from 341 subjects with septic shock within 24 hours of admission to the ICU. Classification and regression tree analysis was used to generate a decision tree predicting 28-day mortality based on a combination of both biomarkers and clinical variables. The derived tree was first tested in an independent cohort of 331 subjects, then calibrated using all subjects (n = 672), and subsequently validated in another independent cohort (n = 209). SETTING: Multiple ICUs in Canada, Finland, and the United States. SUBJECTS: Eight hundred eighty-one adults with septic shock or severe sepsis. INTERVENTION: None. MEASUREMENTS AND MAIN RESULTS: The derived decision tree included five candidate biomarkers, admission lactate concentration, age, and chronic disease burden. In the derivation cohort, sensitivity for mortality was 94% (95% CI, 87-97), specificity was 56% (50-63), positive predictive value was 50% (43-57), and negative predictive value was 95% (89-98). Performance was comparable in the test cohort. The calibrated decision tree had the following test characteristics in the validation cohort: sensitivity 85% (76-92), specificity 60% (51-69), positive predictive value 61% (52-70), and negative predictive value 85% (75-91). CONCLUSIONS: We have derived, tested, calibrated, and validated a risk stratification tool and found that it reliably estimates the probability of mortality in adults with septic shock.

Sepsis in transit: from clinical to molecular classification.

In the previous issue of Critical Care, Maslove and colleagues studied circulating neutrophil transcriptional expression to discover and validate a molecular subclassification of adult patients with sepsis. The authors divided patients into small derivation (n = 55) and validation (n = 71) cohorts. Their complex methodology included partitioning around medoid and hierarchical clustering methods to define two transcriptionally distinct subtypes of sepsis. Pathway analysis found that chemokine and cytokine pathways as well as Toll-like receptor signaling were enhanced. Investigation of specific drug target genes relevant to sepsis found significantly different expression levels between the two molecular subtypes. Interestingly, most patient characteristics did not differ between groups, except for an increase in the proportion of severe sepsis in molecular subtype 1. Possible confounders of this study were the small sample size, population stratification, and lack of information regarding drug interventions, all of which support the need for more studies with larger cohorts that include transcriptional profiles. This thought-provoking hypothesis-generating study could lead to a new neutrophil expression-based molecular classification of adult sepsis.

Association of angiotensin II type 1 receptor-associated protein gene polymorphism with increased mortality in septic shock.

OBJECTIVE: Angiotensin II and its postreceptor signaling are crucial in regulating vasomotor tone. The objective of this study was to test the hypothesis that single nucleotide polymorphisms in angiotensin II pathway genes alter outcome of septic shock. DESIGN: Genetic association study and in vitro experiment. SETTING: Intensive care units at academic teaching centers. PATIENTS: Derivation and validation septic shock cohorts (n = 589 and n = 616, respectively) and a coronary artery bypass surgery cohort (n = 551). INTERVENTIONS: Patients with septic shock in the derivation cohort were genotyped for tag single nucleotide polymorphisms: angiotensin-converting enzyme (six single nucleotide polymorphisms), angiotensin II receptor type 1 (five single nucleotide polymorphisms), and angiotensin II type 1 receptor-associated protein (three single nucleotide polymorphisms), which is a negative regulator of angiotensin II receptor type 1. Patients in the septic shock replication cohort and the coronary artery bypass graft cohort were genotyped for the angiotensin II type 1 receptor-associated protein rs11121816. MEASUREMENTS AND MAIN RESULTS: The primary outcome variable was 28-day mortality. Secondary outcome variables were blood pressure and heart rate. Angiotensin II type 1 receptor-associated protein messenger RNA expression was measured in genotyped lymphoblastoid cells in vitro. Patients with septic shock patients the GG genotype of angiotensin II type 1 receptor-associated protein rs11121816 had increased 28-day mortality in the derivation cohort (54.8% vs. 41.4%; adjusted hazard ratio, 1.46; 95% confidence interval, 1.09-1.93; p = .010 [all ethnicities]; p = .050 [white]) and in the replication cohort (43.8% vs. 32.3%; hazard ratio, 1.42; 95% confidence interval, 1.03-1.98; p = .035 [all ethnicities]; p = .037 [white]). Patients having the GG genotype had decreased mean arterial pressure (98.3% of other genotype, p = .058 [derivation cohort]; 97.7%, p = .00060 [replication cohort]) and increased heart rate (104.1%, p = .023 [derivation cohort], 102.9%, p = nonsignificant [replication cohort]). GG genotype patients undergoing coronary artery bypass grafting had decreased postoperative mean arterial pressure and increased postoperative heart rate (p < .05). GG genotype lymphoblastoid cells had 2.0-fold higher angiotensin II type 1 receptor-associated protein messenger RNA expression (p < .05). CONCLUSIONS: For angiotensin II type 1 receptor-associated protein, the negative regulator of angiotensin II receptor type 1, the GG genotype of rs11121816 was associated with increased angiotensin II type 1 receptor-associated protein expression, decreased blood pressure, and increased heart rate as well as increased 28-day mortality in septic shock.

Leucyl/cystinyl aminopeptidase gene variants in septic shock.

BACKGROUND: Vasopressin is an essential peptide hormone regulating cardiovascular homeostasis and an adjunctive vasopressor therapy for septic shock. METHODS: We tested for association between single nucleotide polymorphisms (SNPs) in vasopressin pathway genes and altered outcome in derivation (n = 589) and replication (n = 616) cohorts of patients with septic shock. The primary outcome was 28-day mortality and the secondary outcome was vasopressin clearance. In a third cardiac surgical cohort (n = 977), we tested for locus-specific heritability of serum sodium concentrations. RESULTS: Of 17 tested tag SNPs in five vasopressin pathway genes (arginine vasopressin [AVP], arginine vasopressin receptor 1A and 1B [AVPR1A, AVPR1B], leucyl/cystinyl aminopeptidase [LNPEP], and oxytocin receptor [OXTR]), rs18059 in LNPEP (also known as vasopressinase) was associated with 28-day mortality in the derivation cohort (P = .037). Therefore, we resequenced the 160-kb haplotype block encompassing the LNPEP gene, including rs18059, and genotyped the 230 identified SNPs in the derivation cohort. The strongest signal was found for LNPEP rs4869317 (adjusted P = .044). The rs4869317 TT genotype was associated with increased 28-day mortality in the derivation cohort (51.0% [TT] vs 34.5% [AA/AT]; adjusted hazard ratio [HR], 1.58; 95% CI, 1.21-2.06; P = .00073) and the replication cohort (38.6% vs 29.6%; HR, 1.36; 95% CI, 1.03-1.80; P = .030). We found that the TT genotype was associated with increased plasma vasopressin clearance (P = .028), and the rs4869317 genotype accounted for 80% of the variance of serum sodium concentrations (locus-specific heritability) in cardiac surgical patients. CONCLUSIONS: The genetic variation in LNPEP (vasopressinase) is associated with 28-day mortality in septic shock and is associated with biologic effects on vasopressin clearance and serum sodium regulation. Further confirmation in additional cohorts is required.

A single nucleotide polymorphism in NF-κB inducing kinase is associated with mortality in septic shock.

We tested the hypothesis that single nucleotide polymorphisms (SNPs) within genes of the NF-κB pathway are associated with altered clinical outcome of septic shock patients. We genotyped 59 SNPs in the NF-κB pathway in a discovery cohort of septic shock patients (St. Paul's Hospital [SPH], N = 589), which identified the C allele of rs7222094 T/C within MAP3K14 (NF-κB inducing kinase; NIK) associated with increased 28-d mortality (uncorrected p = 0.00024, Bonferroni corrected p = 0.014). This result was replicated in a second cohort of septic shock patients (Vasopressin and Septic Shock Trial [VASST; N = 616]) in which the CC genotype of rs7222094 was associated with increased 28-d mortality (Cox regression: SPH cohort hazard ratio [HR], 1.35; 95% confidence interval [CI], 1.12-1.64; p = 0.002 Caucasian only; and VASST cohort HR, 1.24; 95% CI, 1.00-1.52; p = 0.048 Caucasian only). Patients having the CC genotype of rs7222094 in SPH experienced more renal and hematological dysfunction (p = 0.003 and p = 0.011), while patients of the VASST cohort with the rs7222094 CC genotype showed the same trend toward more renal dysfunction. In lymphoblastoid cell lines, we found the rs7222094 genotype most strongly associated with mRNA expression of CXCL10, a chemokine regulated by NF-κB. Accordingly, we measured CXCL10 protein levels and found that the CC genotype of rs7222094 was associated with significantly lower levels than those of the TT genotype in lymphoblastoid cell lines (p < 0.05) and in septic shock patients (p = 0.017). This suggests that the CC genotype of NIK rs7222094 is associated with increased mortality and organ dysfunction in septic shock patients, perhaps due to altered regulation of NF-κB pathway genes, including CXCL10.

beta2-Adrenergic receptor gene polymorphism is associated with mortality in septic shock.

RATIONALE: The CysGlyGln haplotype of the beta(2)-adrenergic receptor gene (ADRB2) is functional and associated with altered responses to adrenergic agonists in patients with asthma. Whether this functional haplotype alters outcome in patients receiving adrenergic agonists in septic shock is unknown. OBJECTIVES: To determine whether genetic variation of ADRB2 influences outcome in septic shock. METHODS: Two cohorts of patients with septic shock were studied: a single center (St. Paul's Hospital [SPH]) cohort (n = 589) and the Vasopressin and Septic Shock Trial (VASST) cohort (n = 616). The A allele of the rs1042717 G/A polymorphism is in complete linkage disequilibrium with the CysGlyGln haplotype of ADRB2; therefore, rs1042717 was genotyped. Modulation by norepinephrine and salbutamol of IL-6 production by stimulated in vitro lymphoblastoid cells was measured by genotype. MEASUREMENTS AND MAIN RESULTS: Patients who had the AA genotype of rs1042717 displayed increased 28-day mortality in SPH (adjusted hazard ratio, 2.23; 95% confidence interval, 1.33-3.72; P = 0.0022), and this result was replicated in VASST (adjusted hazard ratio 2.82; 95% confidence interval, 1.56-5.09; P = 0.0006). This genotypic effect was eliminated in patients treated with acute low-dose corticosteroids. In all patients, the AA genotype was associated with more organ dysfunction. Patients with the AA genotype had a higher heart rate (SPH; P < 0.05; VASST; P < 0.05) and required a higher norepinephrine dose over Days 1 through 3 (VASST; P < 0.05). The AA genotype was associated with decreased norepinephrine and salbutamol inhibition of IL-6 production by stimulated lymphoblastoid cells in vitro (P < 0.05). CONCLUSIONS: The AA genotype of ADRB2 rs1042717, identifying homozygotes for the CysGlyGln haplotype, was associated with increased mortality and more organ dysfunction in septic shock.