Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 19 February 2020
Sec. Microbes and Innate Immunity

Extensive Changes in Transcriptomic “Fingerprints” and Immunological Cells in the Large Organs of Patients Dying of Acute Septic Shock and Multiple Organ Failure Caused by Neisseria meningitidis

\nBerit Sletbakk Brusletto,
Berit Sletbakk Brusletto1,2*Else Marit Lberg,Else Marit Løberg2,3Bernt Christian HellerudBernt Christian Hellerud4Ingeborg Lstegaard Goverud,Ingeborg Løstegaard Goverud2,3Jens Petter Berg,Jens Petter Berg1,2Ole Kristoffer OlstadOle Kristoffer Olstad1Unni Gopinathan,Unni Gopinathan1,2Petter Brandtzaeg,,Petter Brandtzaeg1,2,5Reidun 
vstebReidun Øvstebø1
  • 1Department of Medical Biochemistry, Oslo University Hospital, Oslo, Norway
  • 2Institute of Clinical Medicine, University of Oslo, Oslo, Norway
  • 3Department of Pathology, Oslo University Hospital, Oslo, Norway
  • 4Institute of Immunology, Oslo University Hospital, Oslo, Norway
  • 5Department of Pediatrics, Oslo University Hospital, Oslo, Norway

Background: Patients developing meningococcal septic shock reveal levels of Neisseria meningitidis (106-108/mL) and endotoxin (101-103 EU/mL) in the circulation and organs, leading to acute cardiovascular, pulmonary and renal failure, coagulopathy and a high case fatality rate within 24 h.

Objective: To investigate transcriptional profiles in heart, lungs, kidneys, liver, and spleen and immunostain key inflammatory cells and proteins in post mortem formalin-fixed, paraffin-embedded (FFPE) tissue samples from meningococcal septic shock patients.

Patients and Methods: Total RNA was isolated from FFPE and fresh frozen (FF) tissue samples from five patients and two controls (acute non-infectious death). Differential expression of genes was detected using Affymetrix microarray analysis. Lung and heart tissue samples were immunostained for T-and B cells, macrophages, neutrophils and the inflammatory markers PAI-1 and MCP-1. Inflammatory mediators were quantified in lysates from FF tissues.

Results: The transcriptional profiles showed a complex pattern of protein-coding and non-coding RNAs with significant regulation of pathways associated with organismal death, cell death and survival, leukocyte migration, cellular movement, proliferation of cells, cell-to-cell signaling, immune cell trafficking, and inflammatory responses in an organ-specific clustering manner. The canonical pathways including acute phase response-, EIF2-, TREM1-, IL-6-, HMBG1-, PPAR signaling, and LXR/RXR activation were associated with acute heart, pulmonary, and renal failure. Fewer genes were regulated in the liver and particularly in the spleen. The main upstream regulators were TNF, IL-1β, IL-6, RICTOR, miR-6739-3p, and CD3. Increased numbers of inflammatory cells (CD68+, MPO+, CD3+, and CD20+) were found in lungs and heart. PAI-1 inhibiting fibrinolysis and MCP-1 attracting leukocyte were found significantly present in the septic tissue samples compared to the controls.

Conclusions: FFPE tissue samples can be suitable for gene expression studies as well as immunostaining of specific cells or molecules. The most pronounced gene expression patterns were found in the organs with highest levels of Neisseria meningitidis DNA. Thousands of protein-coding and non-coding RNA transcripts were altered in lungs, heart and kidneys. We identified specific biomarker panels both protein-coding and non-coding RNA transcripts, which differed from organ to organ. Involvement of many genes and pathways add up and the combined effect induce organ failure.

Introduction

Neisseria meningitidis is feared among health care personnel and lay people owing to its propensity to cause acute meningitis and fulminant septicemia, often in clusters of the population (van Deuren et al., 2000; Rosenstein et al., 2001; Stephens et al., 2007). The case fatality rate of systemic meningococcal disease (SMD) has remained fairly stable around 10% since the introduction of sulfa therapy in the 1930s and penicillin in the 1940s (Barquet et al., 1999; Brooks et al., 2006; Stanton et al., 2011; Stoof et al., 2015).

Septic shock is the principal cause of death in patients with SMD in industrialized countries (Gedde-Dahl et al., 1983; Halstensen et al., 1987; de Greeff et al., 2008; Brandtzaeg and Van Deuren, 2012). Persistent septic shock in SMD is a consequence of the rapid proliferation of N. meningitidis in the circulation and in the vasculature of various organs. This is reflected by very high levels of N. meningitidis (106-108) in blood and tissues (104-109) as measured by the number of meningococcal DNA molecules (NmDNA), levels of meningococcal lipopolysaccharides (LPS, endotoxin) in the range of 101-103 endotoxin units (EU)/ mL plasma and an accompanying immune response which is dose-related to the LPS levels in plasma (Brandtzaeg et al., 1989a; van Deuren and Brandtzaeg, 2000; Hackett et al., 2002; Brandtzaeg and Van Deuren, 2012).

The septic shock is characterized by a profound vasodilation, reduced ventricular contractility, vascular leakage reducing the circulation blood volume and a gradual development of inflammation-induced cardiac failure, manifested by increased left ventricle volume and decreasing ejection fraction (Boucek et al., 1984; Mercier et al., 1988; Hagmolen of ten Have et al., 2000; Thiru et al., 2000; Brandtzaeg, 2006; Martin et al., 2019). The septic cardiac failure and vasodilation have been associated with increased production of nitric oxide (NO) in endothelial cells, high levels of tumor necrosis factor (TNF), interleukin 1β (IL-1β), IL-6, vasoactive intestinal polypeptide and persistent high levels of anaphylatoxins C3a and C5a in plasma (Waage et al., 1989a; Baines et al., 1999; Pathan et al., 2004; Brandtzaeg, 2006; Kobsar et al., 2011).

Concomitantly, the coagulation system is triggered by multiple factors and leads to disseminated intravascular coagulation (DIC) and formation of thrombi and hemorrhages in different tissues and organs including the adrenals, lungs, kidneys, skin, and extremities (Ferguson and Chapman, 1948). Formation of thrombi in capillaries and larger vessels is the hallmark of meningococcal septic shock. The balance between coagulation and fibrinolysis is tilted toward coagulation. Monocytes and monocyte derived microvesicles carrying functional active tissue factor appear to initiate the process in plasma (Osterud and Flaegstad, 1983; Ovstebo et al., 2012, 2014; Hellum et al., 2014). Complement factor 5, which is massively activated to C5a and C5b in these patients, augments the coagulation cascade and thrombus formation (Brandtzaeg et al., 1989b; Ovstebo et al., 2014). The fibrinolytic system is severely impaired by plasminogen activator inhibitor 1 (PAI-1). The low levels of protein C and to a lesser extent antithrombin, are associated with purpura fulminance i.e., diffuse thrombotic lesions of the skin. The endothelial cells reveal concurrently reduced levels of thrombomoduline and protein C receptors. Adherence of leukocytes, primarily neutrophils, and meningococci to the endothelial cells may alter the surface and facilitate thrombus formation in the skin, various organs and peripheral extremities (Brandtzaeg et al., 1989c, 1990; Hazelzet et al., 1996; de Kleijn et al., 1998; Faust et al., 2001; Melican et al., 2013). Meningococcal septic shock is often described as overwhelming owing to rapid development of the symptoms, high levels of the bacterium in the circulation and large organs as well as endotoxin levels in plasma higher than measured in any other human infections. This clinical entity is clearly different from most cases of septic shock and multiple organ failure as observed in intensive care units in present day hospitals. However, studies of the gene regulation in lungs, heart, kidneys, liver and spleen of meningococcal infection may provide valuable information in our understanding of multiple organ failure induced by other organisms and in other clinical settings.

Previous studies have applied transcriptomic methods to investigate changes in the gene expression underlying the inflammatory response in meningococcal sepsis and other forms of sepsis (Maslove and Wong, 2014). These studies have primarily been conducted in whole blood or leukocyte fractions (Maslove and Wong, 2014). Studies of meningococcal sepsis in an in vivo porcine model and in patients suggest that the concentration of N. meningitidis may be significantly higher in the major organs as compared to blood (Hellerud et al., 2015; Brusletto et al., 2017). These findings suggest that the local inflammatory response in different organs contributes quantitatively to the organ dysfunction. Yet no studies have so far investigated changes in transcriptional profiles and the corresponding activation of signaling pathways in major organs affected by massive proliferation of N. meningitidis.

The primary aim of this study was to investigate gene expression changes in details induced in patients who died of meningococcal septic shock and to identify activation of signaling pathways that might explain the development of organ dysfunction in these patients. We used both paraffinized and fresh frozen tissues collected during routine autopsies in which the bacterial load in each organ previously had been determined (Hellerud et al., 2015; Brusletto et al., 2017). To our knowledge this has not been done before. A secondary aim was to visualize the immune reaction by immunohistochemistry and quantify various immune cells and specific molecules that might play crucial roles in the pathophysiology which develops in the heart and lungs in lethal meningococcal septic shock. We also compared results between paraffinized and fresh frozen tissues to generate knowledge about the influence of storage methods on human material.

Materials and Methods

Ethics Approval and Consent to Participate

The study was approved by the Regional Medical Ethical Committee of South East Norway (2011/1413C “Translational research, meningococcal disease” and 2011/753 “Studies of invasive meningococcal and pneumococcal disease”). The patients' samples were collected after informed consent from patient parents or relatives and according to the Helsinki declaration. The Director of Public Prosecutions approved the use of forensic material for this research.

Clinical Definitions

Systemic meningococcal disease (SMD) was present if N. meningitidis was cultivated and or confirmed by polymerase chain reaction (PCR) in blood, cerebrospinal fluid (CSF) or by organ tissue examination (Brandtzaeg et al., 1989a; Ovstebo et al., 2004; Brusletto et al., 2017). Severe septic shock was defined as persistent hypotension because of bacterial infection, with a systolic blood pressure <90 mm Hg in adults (≥12 year) and <70 mmHg in children (<12 year), that required fluid therapy and treatment with vasoactive drugs (dopamine, epinephrine, norepinephrine) for at least 24 h or until death (Brandtzaeg et al., 1989a).

Multiple organ failure was defined as: (1) reduced pulmonary function requiring artificial ventilation to maintain an adequate arterial oxygenation and (2) renal failure as reduced creatinine clearance (<60 mL/minute per 1.73 m2 body surface) or pathologically elevated serum creatinine (related to age and collected within 12 hours after admission).

Time Between Hospital Admission and Death

Patients 1–4 (Table 1) died within 270 min of hospital admission. Patient 5 was found dead in bed.

TABLE 1
www.frontiersin.org

Table 1. Overview of patients, clinical data, storage age of tissue, and type of tissue.

Autopsy Procedure

The autopsy procedures were carried out at the Department of Pathology, Oslo University Hospital, Department of Pathology, Stavanger University Hospital, or at the section for Forensic Pediatric Pathology, Oslo University Hospital. Tissue samples from different organs were fixed in 4% buffered–neutral formalin at room temperature for 6–48 h, dehydrated, cleared, embedded in paraffin and cut in 4 μm thick sections. All tissue sections were routinely stained with hematoxylin and eosin (H&E). The FFPE sections were stored at room temperature.

Subjects

Post Mortem Findings

Macro- and microscopic findings have previously been reported in patients No 1–5 (Brusletto et al., 2017). All five patients had hemorrhagic skin lesions and hemorrhagic adrenals.

Formalin-Fixed, Paraffin-Embedded (FFPE) Tissue From Patients With Meningococcal Shock and Multiple Organ Failure (Patients No 1–5)

The formalin-fixed, paraffin-embedded tissues were selected according to histopathological findings; presence of neutrophilic inflammatory infiltrates or thrombi. Small tissue specimens from five lungs, five hearts, four livers, four kidneys, and three spleens were available.

The samples were collected during the routine post mortem examination within 24 h after the patient died. The storage times of the FFPE tissue samples were 11, 10, 6, 5, and 2 years (Table 1).

Fresh Frozen (FF) Tissue From Patients With Meningococcal Shock and Multiple Organ Failure (Patients No 3–5)

Three lungs, three hearts, two livers, three kidneys and three spleens were collected in parallel with the routine post mortem examination and frozen at −80°C for later analysis. The storage times of the FF tissue were 6, 5, and 2 years. The samples had been partially thawed once and examined before this analysis (Hellerud et al., 2015).

Formalin-Fixed, Paraffin-Embedded (FFPE) Tissue From Patients With Acute Non-infectious Death (Controls) (Patients No 6 and 7)

FFPE tissues from two patients with a non-inflammatory disease were used as negative controls; patient No 6 died of a head injury and patient No 7 had a cerebral hemorrhage. The storage time of the specimens was 15 and 3 years, respectively. Tissue from lungs, heart, liver, spleen and kidneys were analyzed. No pathological findings were found in the tissues. The organ samples were collected at routine post mortem examination 24–48 h after death.

Immunohistochemical Staining of FFPE Lungs and Heart Tissue Samples

IHC staining was performed on samples from patients 2, 3, and 7 (Table 1). Sections from lungs and heart were analyzed for presence of macrophages (CD68), neutrophils (MPO), T cells (CD3), and B cells (CD20). The inflammatory markers, plasminogen activator inhibitor-1 (PAI-1), presently also denoted serpine-1, and monocyte chemoattractant protein (MCP-1/CCL2), were also examined. The tissue samples were deparaffinized, rehydrated and demasked in a microwave oven for 24 min in Target Retrieval Solution or Tris/Edta buffer pH 9.1 for ab CD 68 and CD3. The primary antibodies used were anti CD3 1:50 (rabbit polyclonal ab 5690, Abcam), anti CD20 1:2500 (confirm anti-CD20 (L26) monoclonal mouse anti CD ready-to-use (RTU) cat. no 760-2531 Roche Diagnostics), anti CD68 1:3000 (monoclonal mouse anti human CD68 cat. no M0814 Agilent Technologies), anti MPO 1:1000 (rabbit polyclonal cat. no A0398, Agilent Technologies), anti PAI-1 1:400 (rabbit polyclonal ab 66705, Abcam) and anti MCP-1 1:400 (rabbit polyclonal ab9669, Abcam). Antigen-antibody reactions were visualized with DAKO EnVision horse radish peroxidase system (Agilent Technologies or DAKO Cytomation for MCP-1) using 3, 3′-diaminobenzidin as the chromogen. All tissue sections were counterstained with hematoxylin.

Quantification of Immune Cells and Inflammatory Markers

Immunostained sections from heart and lungs from the patients and controls were examined using an Olympus BX51 microscope (obj. x20). Positively stained inflammatory cells were counted in 10 pictures (Diagnostic instruments. Inc, model 11.2 camera) from representative areas, a total area of 2 mm2 per section. The inflammatory markers (PAI-1 and MCP-1) were registered as present or not present.

RNA Extraction

FFPE tissue samples for RNA extraction were cut using a microtome. A new sterile blade was used for each patient and washed with 70% ethanol between each organ block. The first sections of the tissue samples were discarded before cutting sections for RNA extraction. Freshly cut slices of two 10 μm-thick sections were isolated in parallel and immediately placed in 160 μL of deparaffinization Solution (cat. no: 19093 Qiagen) in a microcentrifuge tube. The miRNeasy FFPE kit (Qiagen, Hilden, Germany) was further used for extraction of total RNA in the QIAcube robot www.qiagen.com/MyQIAcube according to manufacturer's instructions. RNase–Free DNase I digestion step was added to remove DNA contamination and highly fragmented molecules in the RNA samples. The RNA was eluted in 20 μL RNase-free water. A negative control (sample without tissue sample) was subjected in parallel for isolation to check for contaminations. The total RNA samples were stored at −80°C before further analysis. All reactions were performed in an RNase-free environment; benches, instruments, and pipettes were cleaned and treated with RNaseZap solution (Ambion Inc., Austin TX) and RNase-free tips and microtubes were used. Total RNA isolated in parallel sections from FFPE tissue samples gave almost similar yields and purities (260/280 ratio) (data not shown). Therefore, a single tissue sample from each tissue was used for microarray analysis.

Approximately 30 mg of FF tissue samples were cut by a sterile scalpel and subjected for homogenization for 2 min in 700 μL QIAzol using a TissueLyzer (Qiagen, Hilden, Germany). Total RNA was furthermore extracted using the miRNeasy Mini Kit (Qiagen) and Phase Lock Gel™ Heavy (5 PRIME GmbH, Hamburg, Germany) and the QIAcube robot (www.qiagen.com/MyQIAcube) according to manufacturer's instructions. Total RNA was eluted in 50 μL RNase-free water and stored at −80°C before further analysis.

Total RNA concentration and purity (260/280 ratio) in each sample was determined with the NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA) and the quality assessed on Agilent 2100 Bioanalyzer RNA 6000 Nano Kits (Agilent Technologies, Palo Alto, CA).

Microarray Analyses

Microarray analyses were performed using the Affymetrix GeneChip Human Transcriptome 2.0 Arrays (Affymetrix, Santa Clara, CA, USA). Total RNA (100 ng) was subjected to the Sensation Plus™ FFPE Amplification and WT Labeling Kit, following the manufacturer's protocol for whole–genome gene expression analysis. Biotinylated and fragmented single-stranded cDNAs were hybridized to the arrays. The arrays were washed and stained using an FS-450 fluidics station (Affymetrix, fluidics protocol FS450_0001). Signal intensities were detected by a Hewlett Packard (Palo Alto, CA, USA) 30007G gene array scanner. The scanned images were processed using the AGCC (Affymetrix GeneChip Command Console) software, and the CEL files were imported into Partek® Genomics Suite™ software (Partek, St.Louis, MO, USA) for statistical analysis. The Robust Multichip Analysis (RMA) algorithm was applied for generation of signal values and normalization. On each array 44,699 protein coding genes (transcript clusters) and 22,829 non–protein coding genes (including immature microRNAs) could be detected. Gene transcripts with maximal signal values of < 5 (log2) across all arrays were removed to filter for low and non-expressed genes. For expression comparisons of different groups, profiles were compared using a one-way ANOVA model. The results were expressed as fold changes (FC). Genes with FC ≥ |±2| and a P-value < 0.05 were regarded as significantly regulated. Partek Genomics Suite software was used to generate principal component analysis (PCA), Venn diagram and a table for top up- and down-regulated transcripts of the gene expression data. To minimize experimental artifacts RNA extraction, sample amplification and labeling, hybridization and washing, and scanning procedures, were carried out by the same operator.

The Affymetrix Synthesis

Staggered concentrations of internal positive controls, poly A RNA controls, which monitor the entire labeling process, and hybridization controls, which monitor the hybridization process, were spiked into the Sensation Plus™ FFPE Amplification and WT Labeling Kit, and found as Present calls in all samples.

Gene Expression Data Analysis

To summarize the information obtained in the microarray analysis and to check data quality, we used Principal Component Analysis (PCA) which identifies the directions (principal components) in which the variation in the data is maximal, and enables us to visualize this variation in a plot.

A Venn diagram illustrating the number of differentially expressed genes (one-way ANOVA, p < 0.05, FDR 5%) in FFPE tissue samples (lungs, heart, kidneys, liver and spleen) from meningococcal septic shock patients compared with controls, was performed using: http://bioinformatics.psb.ugent.be/webtools/Venn/.

The top 10 up- and down-regulated gene transcripts, differentially expressed in large organs from patients with meningococcal septic shock vs. controls were calculated and displayed in a table.

Functional Analysis by Ingenuity Pathway Analysis (IPA)

Bioinformatics analysis was conducted on significantly regulated genes to identify biological functions/pathways that were most significantly associated with the data set by Ingenuity Pathways Analysis (IPA) (Ingenuity Systems, Redwood City, CA, USA). Briefly, the data set containing gene identifiers and corresponding fold changes and p-values was uploaded into the web-delivered application where each gene identifier was mapped to its corresponding gene object in the IPA Knowledge Base (IPKB). Fisher's exact test was performed to calculate p-values determining the probability that each biological function and/or disease assigned to the data set was due to chance alone. The program then computed a score for each network according to the fit of the network to the set of focus genes. The score was derived from a p-value which indicates the likelihood of the focus genes in a network being found together because of random chance. The Z-score indicates predicted activation state of the biofunctions, canonical pathways, and upstream regulators. The data sets were mined for significant pathways with the IPA library of canonical pathways and networks were generated by using IPA as graphical representations of the molecular relationships between genes and gene products.

Enrichment analysis was performed using IPA's ≪core analysis≫ for each tissue sample and a ≪comparison analysis≫ between the tissue samples. These functions have the ability to identify significantly activated biological functions and pathways, molecular functions, and relationships in our dataset of genes. A right-tailed Fisher's exact test calculated p-values corrected for multiple testing by the Benjamini-Hochberg method. In addition an “Upstream regulator” analysis was used to identify the cascade of upstream transcriptional regulators that were involved in MSS patients and whether they were likely activated or inhibited to obtain the observed gene expression profile changes in our datasets.

Canonical Pathway Analysis in IPA

There are two main groups of canonical pathways in IPA: metabolic and signaling. These pathways are hierarchically grouped according to a number of sub-categories. In order to identify those pathways most relevant to the cell-types and disease of our study, disease- and cell-specific pathways not deemed relevant to our study were excluded from the analysis. Canonical pathways significantly enriched by the differentially expressed genes in the datasets were identified with the right-tailed Fisher's exact test, which calculates a P-value determining the probability that the canonical pathway is associated with the data set due to random chance alone. The P-values were corrected for multiple testing using the Benjamini-Hochberg method for correcting the false discovery rate (Benjamini and Hochberg, 1995).

Validation of Genes by Real Time Quantitative Reverse Transcription PCR (qRT-PCR)

To validate the gene expression data from FFPE and FF tissue samples analyzed by microarray analysis, selected differentially expressed genes with the significantly highest fold change value were quantified by qRT-PCR (TaqMan gene expression assays and the Applied Biosystems ViiA7 sequence detection system). Total RNA (100 ng) was reverse transcribed using SuperScript™ VILO™ cDNA Syntheisis Kit (cat. no: 11754250) (ThermoFisher Scientific). In a total volume of 20 μL; 20 ng cDNA (2 μL), and 1 μL of either RT-PCR primer cat 4331182 (CCL2 Hs 00234140-m1, IL1R Hs 01073300-m1, CXCL8 Hs 00174103-m1, IL1B Hs 01555410-ml, IL6 Hs 00174131-mL, SERPINE 1 Hs 00167155-m1, MT1A Hs00831826-s1, RPL9 Hs 01552541-g1, PLA2G2A Hs 00179898-m1, HAMP Hs 00221783 –m1, DYRK2 Hs 00705109-s1) were added to 10 μL TaqMan Fast Advanced mastermix (Applied Biosystems, www.appliedbiosystems.com) and 7 μL of H2O. The relative changes of each transcript using the mean of DYRK2 (Hs 00705109-s1) as endogenous control, were calculated using the ViiA™7 Software v1.1 and the ΔΔCT method (Livak and Schmittgen, 2001).

Validation of Proteins by Multiplex Assays

Homogenized FF tissue samples were used to examine whether mRNAs are translated into proteins. Fifty mg of tissue sample, 495 μL CytoBuster Protein Extraction Reagent (Novagen, San Diego, CA) and 5 μL Protease Inhibitor cocktail set I (Calbiochem, Darmstadt, Germany) was homogenized with Xiril Dispomix. After completion, the samples were incubated for 5 min on ice and thereafter centrifuged at 2500 x g for 20 min at 4°C. The supernatants were transferred to Nunc tubes (Thermo Fisher Scientific Inc., Waltham, MA) and stored at −80°C before analysis. The supernatants were analyzed using multiplex cytokine assay (Bio-Plex Pro™ Human Cytokine assay), (Bio-Rad Laboratories, Hercules, CA) by Luminex® Technology.

Quantification of N. meningitidis DNA and LPS in Plasmas/Serum/CSF From Patients With Meningococcal Disease in Samples Collected on Hospital Admission (Table 1)

Heparin-blood was collected in LPS-free vacuum tubes, centrifuged, plasma pipetted off and aliquoted as described in detail earlier (Brandtzaeg et al., 1989a; van Deuren and Brandtzaeg, 2000). Quantification of N. meningitidis DNA was performed as previously described in detail (Ovstebo et al., 2004; Gopinathan et al., 2012). The detection limit was 103 N. meningitidis DNA copies/mL.

Quantification of LPS in plasma/serum/ CSF was initially performed with an in house developed limulus amebocyte lysate (LAL) assay and later with Chromo-LAL (Associates of Cape Cod, USA) with a detection limit of 0.2 EU/mL. The serum level is on average 63% of the plasma level (Brandtzaeg et al., 1992; van Deuren and Brandtzaeg, 2000).

Quantification of N. meningitidis DNA in FFPE and FF Tissue From Patients With Meningococcal Disease in Samples Collected Post Mortem (Table 1)

All FFPE tissue samples were prepared according to routine procedures at the Department of Pathology. FF tissue were collected in parallel with the routine post mortem examination and frozen at −80°C for later analysis. Quantification of N. meningitidis DNA was performed as previously described in detail (Brusletto et al., 2017).

Statistical Analysis

GraphPad Prism, version 7 was used for Standard error of the mean (SEM), medians and correlation plots (GraphPad Software, Inc.). Fischer's exact test was used when appropriate.

Results

Evaluation of RNA and Affymetrix Synthesis

RNA Extraction: Yield and Purity

We compared the concentration, purity and RIN of RNA from the different tissues samples. The concentration of RNA isolated from FFPE tissue samples (2 × 10 μm) from five patients with meningococcal septic shock ranged from 59 to 624 ng/μL, (n = 21 median 119 ng/μL). The purity ranged from 1.62 to 2.00 (260/280 ratio) (n = 21 and median 260/280 ratio = 1.89). RIN ranged from 2.3 to 2.6 (median 2.5 in RIN).

The concentration of RNA isolated from FFPE tissue samples from controls (non-infectious disease) ranged from 41 to 439 ng/μL (n = 8 and median 62 ng/μL). The purity ranged from 1.66 to 1.91 (260/280 ratio) (n = 8 and median 260/280 ratio = 1.8. RIN ranged from 1.9 to 2.5 (median 2.4 in RIN). The concentration of RNA isolated from FF tissue samples (30 mg) from three patients with meningococcal disease ranged from 24 to 2762 ng/μL (n = 15 and median 716 ng/μL). The purity ranged from 1.96 to 2.08 (260/280 ratio) (n = 15 and median 260/280 ratio = 2.05). RIN ranged from 2 to 8.5 (median 3.5 in RIN).

Overall, these samples showed adequate RNA quantity and quality for further amplification and labeling using Affymetrix Sensation Plus FFPE kit allowing RIN values > 1.4 to be acceptable.

Evaluation of the Gene Expression Profiles

Principal Component Analysis: Sources of Variability and Effect of Storage Time

Principal component analysis was used to determine the significant sources of variability, identify patterns in the data sets, and if storage time had effect on gene expression patterns. The PCA analysis (Additional File 1: Figure S1) indicated that the two groups, meningococcal septic shock patients and controls, cluster with minor overlap, indicating different gene expression patterns. The storage time of the tissue samples at the time of analysis varied greatly within the study population. The effect of degradation of the nucleic acids was evaluated using PCA. Despite low RIN values, the PCA plot (Additional File 1: Figure S1) indicates that the age of the tissue samples had no impact on the gene expression analysis, showing no clustering of young or old tissue samples.

Gene Expression Data in FFPE Tissue Samples vs. FF Tissue Samples

The gene expression profiles from fresh frozen (FF n = 3) and formalin-fixed, paraffin embedded tissues (FFPE n = 5) tissue samples (lungs, heart, kidneys, liver, and spleen) from meningococcal shock patients were compared to identify if gene expression profiles were different between the storage methods. Correlation plots (Figure 1) of the gene expression profiles showed correlations ranging from r = 0.88 to 0.97, indicating good consistency between the methods.

FIGURE 1
www.frontiersin.org

Figure 1. Correlation plots of the gene expression profiles between FF and FFPE tissue samples from meningococcal septic shock patients. The mean signal intensities for unfiltered gene expression (25000 transcripts) with NM (m-RNA) and NR (ncRNA) annotations on FF (x-axis) and FFPE (Y-axis) tissue samples for respectively lungs, heart, kidneys, liver and spleen. Pearson correlation coefficients (r) for each scatter plot are shown on the top.

The selected transcripts (CCL2, IL1R, CXCL8, SERPINE 1, MT1A, IL-6, RPL9, PLA2G2A, HAMP, IL-1B) quantified by qRT-PCR in FFPE and FF, showed correlation coefficients for lungs r = 0.896, heart r = 0.816, kidneys r = 0.936, liver r = 0.413 and spleen r = 0.521.

Immunohistochemical Staining of Lungs and Heart

Overall, the lungs had a higher number of infiltrative inflammatory immune cells than the heart. CD68-positive macrophages, MPO-positive neutrophils, CD3-positive T-lymphocytes and CD20–positive B-lymphocytes were present in higher numbers in the lungs and the heart of meningococcal septic shock patients as compared to the control (Figure 2 and Table 2). In the lungs, large CD68-positive macrophages were observed in the alveolar spaces whereas smaller CD68-positive cells were present within the alveolar walls and small vessels (Figure 2 and Table 2). More endothelial and inflammatory cells were significantly positive for PAI-1/SERPINEI and MCP-1 in the lungs and heart from the meningococcal septic shock patients compared with the controls (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Immunohistochemical examination of FFPE tissue samples (lungs and heart) from one meningococcal septic shock patient and one control (acute non-infectious death). Immunohistochemical detection of CD68, MPO, MCP-1, and PAI-1 positive cells (labeled brown) in sections from lungs and heart from meningococcal septic shock patient and control patient with acute non-infectious death (Obj. x20) (#The control patient, elderly smoker).

TABLE 2
www.frontiersin.org

Table 2. Immunohistochemically examination of FFPE tissue samples (lungs and heart) from two patients with meningococcal septic shock and one control patient (acute non-infectious death).

Gene Expression Profiles in FFPE Tissue Samples From Meningococcal Septic Shock Patients and Controls

PCA was used to visualize variations in gene expression profiles from the different organs from meningococcal septic shock patients (n = 5) and controls (n = 2). Overall, the PCA plot demonstrated an organ specific gene clustering of the data sets (Additional File 2: Figure S2).

A comparison of the lists from the transcripts from lungs, heart, kidneys, liver and spleen using a Venn diagram (Figure 3) showed 827 specific transcripts for lungs, 982 for heart, 837 for kidneys, 559 for liver and 182 for the spleen. Five transcripts including 1 mapped transcript; DMTF1, a transcription regulator in the senescent pathway and 4 transcripts with unmapped ID in IPA, were found to be common for these five organs.

FIGURE 3
www.frontiersin.org

Figure 3. Comparison of differentially expressed gene transcripts in FFPE tissue samples (lungs, heart, kidneys, liver, and spleen) from meningococcal septic shock patients vs. control patients (acute non-infectious death). Gene transcripts with signal values of <5 (log2) across all arrays were removed to filter for low and non-expressed genes. For expression comparisons of different groups, the profiles were compared using a one way ANOVA model. P-value <0.05, correction FDR 5%. The Venn diagram illustrate that 827 gene transcripts were specific for the lungs of the 2039 differentially expressed gene transcripts, in heart 982 of 2029, in kidneys 837 of 2231, in liver 559 of 1531 and in spleen 182 of 435. Five gene transcripts were common for lungs, heart, kidneys, liver, and spleen.

The top 10 most up-and down-regulated differentially expressed transcripts in the organs (lungs, heart, kidneys, liver, and spleen) from patients with meningococcal septic shock vs. controls also showed an organ specific gene expression where transcripts for different classes of non-coding RNAs such as transfer RNA and piRNA were most pronounced in most organs (Table 3). In lungs and heart, also transcripts resulting in proteins involved in inflammation, were on the top 10 up-regulated list.

TABLE 3
www.frontiersin.org

Table 3. Top 10 differentially expressed up-regulated (A) and down-regulated (B) gene transcripts in large organs from patients with meningococcal septic shock vs. controls.

Descriptive Analysis and Predicted Biological Functions of Differentially Expressed Genes in FFPE Tissue Samples From Meningococcal Septic Shock Patients

Table 4A shows the numbers of differentially expressed genes in the meningococcal septic shock tissue samples compared to control tissue samples and their association with molecular, cellular, and physiological functions that were identified in lungs, heart, kidneys, liver, and spleen. Based on the number of differentially expressed genes, all tissue samples were affected, with least affection in the spleen. The highest number of affected molecules (Table 4B) was associated to cellular movement, cell to cell signaling and interaction and cell death and survival with lungs and heart as most influenced.

TABLE 4
www.frontiersin.org

Table 4. Numbers of genes differentially expressed (A) and functional analysis (B) in organs from patients with meningococcal septic shock (n = 5) performed by IPA.

When gene expression changes in all the tissue samples from meningococcal septic shock patients were compared, the IPA “comparison analysis” predicted the down-regulated biofunctions to be organismal death, migration of cells, cell movement, necrosis and cell death (Figure 4A). The top regulated canonical pathways were predicted to be the acute phase response signaling, EIF2 signaling, TREM1 signaling, IL-6 signaling and HMGB1 signaling whereas PPAR signaling, LXR/RXR activation, PPARα/RXRα activation and NF-κB to be the most down-regulated (Figure 4B). Among predicted genes in the top five up-regulated canonical pathways, SERPINE 1 was on top in both acute phase signaling and HMGB1 signaling while MCP-1 was on top in TREM 1 and HMBG1 signaling (Figure 5A). Top predicted genes in the top down-regulated canonical pathways were IL1RL1, MCP-1, TNFAIP3, NFKB1A, and genes for the heat shock protein 90 family (Figure 5B). To identify the cascade of upstream transcriptional regulators that can explain the gene expression changes in our data sets, an upstream regulator ≪core analysis≫ was performed in IPA. The ≪core analysis≫, showed highly significant p-values (ranging from E35 to E3) (Table 5) of the top upstream regulators in the tissue samples. A ≪comparison analysis≫ of upstream regulators in the different tissue samples predicted TNF, IL1B, IL6, IFNG, and IL1A to be on top concerning up-regulation (Figure 4C) while RICTOR, CD3, and many miRNAs were down-regulated (Figure 4C).

FIGURE 4
www.frontiersin.org

Figure 4. Predicted biofunctions (A), canonical pathways (B), and upstream regulators (C) in FFPE tissue samples from patients with meningococcal septic shock. ≪Comparison analysis≫ of biofunctions, canonical pathways, and upstream regulators significantly enriched in FFPE tissue samples from lungs, heart, kidneys, liver and spleen from meningococcal septic shock patients vs. controls patients (acute non-infectious death). The figure shows the most up-regulated biofunctions (A), canonical pathways (B), and upstream regulators (C) ranked according to expression levels in FFPE tissues from lungs. The Z-score indicates predicted activation state of the biofunctions, canonical pathways, and upstream regulators. Z-score value >|±1| are displayed. Blue color or lighter shades of blue indicate a negative Z-score and down-regulation of a biofunction, canonical pathways, and upstream regulators. Orange or lighter shades of orange indicates a positive Z-score and up-regulation of a biofunction, canonical pathways, and upstream regulators. Note that only the top pathways are shown.

FIGURE 5
www.frontiersin.org

Figure 5. Predicted gene signaling pathways from the top up-regulated (A) and down-regulated (B) canonical pathways in FFPE tissue samples from patients with meningococcal septic shock vs. controls. The Z-score indicates predicted activation state of canonical pathways. Orange or lighter shades of orange indicate a positive Z-score and up-regulation of the pathway. Blue color or lighter shades of blue indicate a negative Z-score and down-regulation of the pathway. The transcripts in the gene signaling pathway are expressed as Fold Change (FC) values. Red or lighter shades of red indicates positive FC-values and up-regulation of transcripts, green color or lighter shades of green indicates negative FC- values and down-regulation of transcripts. Color gray indicates that a predicted activation state of a gene/transcript in the canonical pathway signaling network is not affected.

TABLE 5
www.frontiersin.org

Table 5. Top upstream regulators differentially expressed in organs from five patients with meningococcal septic shocka vs. controls.

An IPA ≪core analysis≫ was performed separately for each organ to fine tune the top enriched pathways for each organs (Figures 6, 7A–E).

FIGURE 6
www.frontiersin.org

Figure 6. Transcriptional profiles of canonical pathways in FFPE tissue samples from meningococcal septic shock patients. The figure displays the top canonical pathways enriched in each organ. The mean values of N. meningitidis DNA/μg human DNA was in lungs: 2.5 × 108, heart 5.4 × 107, kidneys 4.6 × 106 liver 8.3 × 107, and spleen 9.4 × 104. A ≪core analysis≫ was performed separately for each organ. Significantly enriched canonical pathways were identified with a right-tailed Fisher's Exact Test (p < 0.05, after correction for multiple testing using the Benjamini-Hochberg method). Ratio denotes the number of significantly expressed genes compared with the total number of genes associated with the canonical pathway. The Z-score=|±1| indicates predicted activation state of canonical pathway. Blue color or lighter shades of blue indicate a negative Z-score and down-regulation of the pathway, and orange or lighter shades of orange indicate a positive Z-score and up-regulation of the pathway. Gray color indicates no activity pattern available. Z-score value >|±1| are displayed.

FIGURE 7
www.frontiersin.org

Figure 7. Interaction networks of top score canonical pathways and transcripts affected in FFPE tissue samples from meningococcal septic shock patients. The interaction networks in lungs (A), heart (B), kidneys (C), liver (D), and spleen (E) are generated through the use of IPA ≪core analysis≫. In addition, ≪Canonical Pathways Analysis≫ was performed to determine genes that were involved in well-documented canonical signal transduction or metabolic pathways, from the library of canonical pathways in IPA. Benjamini–Hochberg procedure for multiple testing corrections was performed. Significantly enriched genes in the dataset were overlayed into the specific canonical pathway. The analysis was performed in silico using Molecular Activity Predictor analysis of IPA. The shapes represent the molecular classes of the ≪gene≫. In the figures, red represents upregulation and green down-regulation, and color intensity represents the relative magnitude of change in gene expression. Gray color indicates no activity pattern available. Direct and indirect interactions are indicated by solid and dashed lines, respectively. The network diagram shows the biological relationship between the indicated genes lines: — represents direct physical interactions; —– represents indirect functional interactions; → represents activation; ⊣ represents inhibition. The blue lines indicate that the direction of regulation is consistent with IPA prediction. In contrast, yellow lines indicate that the regulation observed is inconsistent with expectations, while gray lines indicate lack of pre-existing data to formulate expectations. Nodes are displayed using various shapes that represent the functional class of the genes.

Our mRNA data in the lung tissue samples (Figures 6, 7A and Additional File 3: Table S1) showed a significant up-regulation of the Toll-like receptor-, acute phase response-, TREM1-, IL-6-, HMBG1-, p38MAPK-signaling, LXR/RXR activation, IL-1 signaling, PPARα/RXRα activation and NF-κB signaling pathways, resulting in striking up-regulations of several cytokines; CCL2, IL-6, IL1B, CCL3 (MIP-1α), transmembrane receptors; IL1RL, ITGA5, transcription regulators; NFKBIA, kinases; MAP2K3, IRAK3 peptidases; PLAT, enzymes; SOD2, FTL, TNFAIP3 and other molecules; such as IL1, SERPINE1, SERPINA1, HAMP, MKK3/6, TNFAIP6, and TNIP1.

The mRNA data in the heart tissue samples (Figures 6, 7B and Additional File 3: Table S2) showed a significant up-regulation of the EIF-2, TREM1-, P13/AKT-and HMGB1-signaling pathways and significant down-regulation of the PPAR signaling pathway (Figure 6) which led to activation of several cytokines; CCL2, CCL3, CXCL8, transmembrane proteins SELE (E-selectin), IL1RL1 (Interleukin 1 receptor –like 1) (ST2), ITGA5 (Integrin alpha-5), ICAM (Intracellular Adhesion Molecule), transcription regulator; STAT3, translation regulator; PABPC1, enzymes; HNRNPA1, HSPA5, and other molecules; Hsp90, IL1, SERPINE1, and YWHAE.

mRNA data in the kidneys tissue samples revealed a significant up-regulation of the EIF2 signaling pathways and the acute phase response (Figures 6, 7C and Additional File 3: Table S3) resulting in strikingly up-regulations of several peptidases such as; C3, HP, SOD2, transcription regulators; EIF1, EIF2, FOS, EIF2S3, and other molecules; such as SERPINE1, SERPINA1, SERPINA3 and 43S Translation Preinitiation.

Validation of Gene Expression Profiles With qRT–PCR

We validated the microarray results with real time qRT-PCR, for genes that were significantly up-regulated in each tissue. Some of the genes were only up-regulated in one tissue, where as other genes were up-regulated in many tissues. The correlations were calculated from the fold changes in all the meningococcal shock patients compared with the controls, determined by qRT-PCR and microarray analysis giving the following results: for lungs r = 0.884 (MT1A, CCL2, SERPINE1), heart r = −0.997 (CCL2, HAMP, IL1RL1), kidneys r = −0.767 (CCL2, RPL9, CXCL8), liver r = 0.977 (PLA2G2A, SERPINE1, CXCL8) and spleen r = 0.376 (SERPINE1, MT1A, CCL2).

Quantification of Selected Proteins in FF Tissue Samples From Meningococcal Septic Shock Patients

Proteins were selected for quantitation by Luminex multiplex analysis based high gene expression FC values when comparing tissues from MSS patients and controls (Figure 8). Tissue sample concentrations for 10 of the cytokines (TNF-α, IL-1β, Il-6, IL-8, IL-10, IP-10, G-CSF, IL-17, MIP-1β, and MCP-1) have previously been published in Hellerud et al. (2015) whereas IL-1ra, RANTES, G-CSF, IL-17, M-CSF, ICAM, MIF, and PAI-1 were analyzed for this publication. High levels of most cytokines, with individual variations, were found in the tissue samples (Hellerud et al., 2015).

FIGURE 8
www.frontiersin.org

Figure 8. Cytokine concentration in fresh frozen tissue samples from patients with meningococcal septic shock (n = 3). The concentration unit is ng/g for the organ samples.

Quantification of N. meningitidis DNA and LPS in Plasma/Serum or CSF From Patients With Systemic Meningococcal Disease and Shock

The number of N. meningitidis/mL in the circulation of the patients with meningococcal shock ranged from 3.0 × 107/mL to 2.8 × 108/mL in plasma/serum (Table 1). LPS in plasma or serum ranged from 271 EU/mL to 3800 EU/mL (Table 1).

Quantification of N. meningitidis DNA in FFPE and FF Tissue Samples From Patients With Systemic Meningococcal Disease and Shock

Copies of N. meningitidis DNA/μg human DNA in FFPE tissue samples in patients with systemic meningococcal disease and shock ranged from 8.1 × 104 to 1.2 × 109 (Table 1) and from 4.2 × 106 to 2.4 × 108 in FF tissue samples (Table 1) (Brusletto et al., 2017). The mean values of N. meningitidis DNA /μg human DNA was in lungs: 2.5 × 108, heart 5.4 × 107, kidneys 4.6 × 106 liver 8.3 × 107 and spleen 9.4 × 104.

Discussion

Post mortem tissue samples (FFPE) from patients with meningococcal septic shock provide an invaluable resource for profiling of gene expression and identifying affected molecular mechanisms. Our material is unique since both tissue samples frozen immediately after the autopsies as well paraffinized tissues stored at room temperature from the same organs in three of the five patients are included. This made it possible to study the influence of storage on the transcriptional profiles in samples stored from 2 to 6 years (Figure 1). The results suggest that FFPE tissues reflect fairly accurately the transcriptional profiles at the time of autopsy. Such studies have been lacking in gene expression studies of human bacterial infectious diseases. Specific proteins detected in the FFPE and FF tissue samples with immunohistochemistry techniques (Figure 2) or multiprotein assays (Figure 8), respectively manifested the regulation of the specific molecular pathways at protein levels.

Nucleic acids have historically been difficult to extract from FFPE tissue samples due to the need to remove paraffin and to counteract covalent protein-DNA interactions that results from the fixation process (Krafft et al., 1997; Farragher et al., 2008). The quality of nucleic acids from FFPE specimens may furthermore be affected by several other factors, such as variability in the time to tissue handling after onset of death (24–48 h), degree of autolysis (Gupta et al., 2012), archival storage time and various tissue sources (Chung et al., 2008; van Maldegem et al., 2008; Abdueva et al., 2010; Ferruelo et al., 2011). Isolating RNA has also been challenging since the FFPE archival methods may lead to chemical modifications and a partial degradation of the RNA (up to 50% of the RNA may not contain an intact poly-A-tail) (Farragher et al., 2008; Mittempergher et al., 2011) complicating established downstream analysis (Farragher et al., 2008). However, commercially available RNA extraction kits now enables to pick up fragmented RNA/cDNA isolated from FFPE tissue samples in a way that satisfy the quality and quantity needed for microarray analysis (Xiang et al., 2003; Zhu et al., 2016).

The quality control of isolated RNA (Boeckx et al., 2011; Ludyga et al., 2012; Seiler et al., 2016) from all the tissue samples in this study showed RNA integrity numbers (RIN) ≥1.9. The FFPE tissue samples ranged from 1.9 to 2–6 while the RIN of FF tissue samples ranged from 2 to 8.5. A RIN value of 1 indicates completely degraded RNA, whereas an RIN value of 10 indicates intact RNA. Our data indicate that the RNA isolated from the FFPE tissue samples from both meningococcal septic shock patients and controls, is heavily degraded. However, a RIN value of 1.4 and RNA purity measured as the 260/280 ratio of ≥1.8, has been considered as a minimum for successful completion of further microarray analyzes using the Affymetrix Sensation Plus FFPE amplification assay (Schroeder et al., 2006; Ribeiro-Silva et al., 2007; von Ahlfen et al., 2007).

Our study included FFPE tissue samples or fresh frozen tissue samples from five different organs from the same three meningococcal septic shock patients. To investigate the effect of the storage methods, the mean signaling intensity of the gene expression profiles from FFPE and FF tissue samples were compared and a correlation analysis resulted in r values ranging from 0.88 to 0.97 (Figure 1). These data indicate that RNA isolated from FFPE–and FF tissue samples give similar results and that microarray analysis enables to pick up and analyze the RNA molecules for gene expression profiles despite differences in storage time and low RIN values in the tissue samples. This is in line with other gene expression studies on FFPE stored cancer tissue samples (Linton et al., 2009; Zhu et al., 2016).

FFPE tissue samples from two persons with acute non–infectious death stored for 3 and 15 years were available as controls. A principal component analysis including information on storage times and storage methods (Additional File 1: Figure S1) indicated a clustering between the patient groups and that storage times and methods had minor effect on gene expression methodology, consistent with other studies (Thomas et al., 2013; Tyekucheva et al., 2015; Webster et al., 2015; Zhu et al., 2016). A PCA plot, a Venn diagram and a top 10 list of up-and down-regulated transcripts also suggests an organ specific transcriptional pattern clustering (Additional File 2: Figure S2, Figure 3 and Table 3). A ≪comparison analysis≫ in IPA was furthermore used to identify affected biofunctions and pathways in the different FFPE tissue samples and a disease associated specific tissue clustering was found in the septic-lungs, -hearts, -kidneys, -livers, and -spleens (Figures 4A–C) compared to controls. This in line with what is presented in Genotype-Tissue Expression (GTEx) - a database collecting gene expression data from numerous healthy and diseased tissue samples which observe tissue specific gene expression profiles (Carithers et al., 2015).

When comparing the multitude of gene expression results we see the contours of a very complex process resulting in a profound reprogramming with changes of gene expression in thousands of genes, both protein-coding and non-coding RNA transcripts. Many of these genes are related to the host's inflammatory response. However, also genes regulating metabolic processes, protein synthesis and folding, transmembrane calcium transport and mitochondrial functions, not directly related to inflammation, are altered probably as a consequence of inflammatory molecules with their multitude of effects. Notably, several non-coding RNAs, including long ncRNA, transfer RNA, piRNA and miRNA, were identified significantly regulated. These genetics elements are suggested to play essential roles in transcriptional and post-transcriptional regulation in innate immunity, mitochondrial- and organ dysfunction in trauma and septic patients (Ho et al., 2016a; Zhang et al., 2017). Given this complexity, it is impossible in this study, to pinpoint one single mechanism or genetic pathway which is the dominant cause of the organ dysfunction in the tissue samples. Our results suggest that involvement of many different genes and pathways may add up and the combined effect induce organ failure. The hypothesis that dysfunction of one “major” pathway leads to organ failure, may be too simplistic. We also find that the differential gene expression of these pathways differ from organ to organ (Figure 4). Some of our findings have previously been documented in a porcine meningococcal septic shock model (Hellerud et al., 2015), but not to our knowledge in human material (Ricote et al., 1998; Castrillo and Tontonoz, 2004; Pascual et al., 2005; Kidani and Bensinger, 2012; Soares et al., 2014; Khan et al., 2015; Hotchkiss et al., 2016; Balmer and Hess, 2017).

The main up-regulated bio-functions were involved in migration of cells, cell movement, cellular homeostasis, leukocyte migration, cell survival, chemotaxis, cell viability, and inflammatory responses while the main repressed cellular functions were associated with organismal death, necrosis, cell death and apoptosis (Table 4B, Figure 4A). Our data show that the pathophysiology observed in meningococcal septic shock tissue samples with different levels of N. meningitidis (Brusletto et al., 2017) may involve a complex interaction between the N. meningitidis and the host's immune system, most likely through recognition of LPS molecules embedded in the outer membrane of the bacteria. However, non-LPS bacterial molecules may also contribute to the inflammatory response (Hellerud et al., 2010).

The main up-regulated canonical pathways in our data set were acute phase response-, EIF-2-, TREM1-, IL-6- and HMBG1 signaling pathways (Figure 4B). Our study identified immunomodulating molecules such as SERPINE1 (PAI-1), CCL2 (MCP-1), CXCL8 (IL-8), ICAM, IL-6, SOD-2, NFKB1A, STAT3, IL1RL1, and RPL (Figure 5A) responding to the microbial invasion consistent with previous studies investigating dysregulated systemic inflammatory responses in sepsis/septic shock (Khan et al., 2015; Hotchkiss et al., 2016).

Among the most interesting findings was the down-regulation of pathways involving several lipid-activated transcription factors, namely nuclear receptors involved in activation of PPAR, LXR/RXR and PPARα/RXRα canonical pathways (Figures 4B, 5B). These are important for physiological lipid and cholesterol metabolism as well as for inflammation (Castrillo and Tontonoz, 2004). Short-term and mild changes in metabolism can positively modulate immune responses to eliminate pathogens and protect the host via disease tolerance (Soares et al., 2014; Van Wyngene et al., 2018). However, uncontrolled and severe disturbances of metabolic homeostasis are unfavorable (Balmer and Hess, 2017). The pathways are expressed in numerous tissues (Ricote et al., 1998) and have been shown to be implicated in the negative regulation of inflammatory responses to bacterial infections via anti-inflammatory or phagocytic roles by binding of transcription factors to co-repressor complexes and thereby maintaining inflammatory genes in a repressed state (Pascual et al., 2005; Kidani and Bensinger, 2012). Our results show that these pathways are variably down-regulated by N. meningitidis in the different organs (Figure 4B), indicating that the down-regulation might be a mechanism to dampen local inflammatory responses.

Quantification of cytokines in lysates, obtained from FF tissue samples from lungs, heart, kidneys, liver and spleen, showed that the organs under these circumstances might synthesize both pro- and anti -inflammatory cytokines in large quantities. Our data (Figure 8) demonstrated increased levels of selected cytokines in all tissue samples, with surprisingly high levels in the spleen. In the spleen, the transcriptional profiles changed minimally as compared with the transcriptional profiles in the other organs, given at the strictly chosen criteria in this study. This may possibly indicate a more profound posttranslatory regulation of the cytokine production in the spleen than in the other organs examined in this study. Most organs probably release these inflammatory mediators into the circulation. High levels of TNF, IL-1, IL-6, IL-10, and chemokines have been measured in plasma from patients with meningococcal sepsis (Waage et al., 1989b; Brandtzaeg et al., 1996b; Hazelzet et al., 1996; Moller et al., 2005; Brandtzaeg, 2006). A notable observation in this study is that comparatively much lower quantities of IL-10, the major inhibitor of proinflammatory cytokines are detected in different tissues as compared with blood levels (Lehmann et al., 1995; Brandtzaeg et al., 1996b; Hellerud et al., 2015). This may imply that interactions between pro- and anti-inflammatory molecules operating at tissue levels are quantitatively different in organ tissues as compared with the circulation. Importantly, many cytokine receptors are released in plasma during meningococcal septic shock which presumably dampen the proinflammatory immune reaction (van Deuren et al., 1995, 1997). However, we have no knowledge about the release of these receptors in the different tissues and the functional inhibitory effects of such release.

Transcriptional Profiles in Lung Tissue Samples

Pulmonary function is heavily affected during meningococcal disease with capillary leakage leading to increased intra-alveolar fluid, pulmonary edema, tachypnea, and respiratory failure (Ferguson and Chapman, 1948; Pathan et al., 2003). The lung capillary cells are “bombarded” by a range of proinflammatory mediators including cytokines and chemokines possibly partly produced locally or conveyed via plasma. Our results suggest that the lungs appear to be key organs for intravascular leukocyte adherence during meningococcal sepsis. This is in line with the results from a porcine model (Hellerud et al., 2015). A striking observation was that the neutrophils in the meningococcal septic shock patients were detected within the alveolar walls and small vessels, as compared to patients with bacterial pneumonias where the neutrophils are located in the alveolar spaces. Leukopenia, particularly a low neutrophils count, is a striking laboratory observation in the meningococcal septic shock patients (Hazelzet et al., 1996; Brandtzaeg et al., 2001). Immunohistochemical analysis of lung tissues in this study documents the accumulation on neutrophils, T-lymphocytes and macrophages in the lung capillaries (Table 2 and Figure 2). A previous study of meningococcal shock patients suggests that the neutrophils are activated simultaneously with an intravascular adherence reflected by high plasma levels of neutrophil-specific elastase (Brandtzaeg and Kierulf, 1992). SERPINA1 (alpha-1-antitrypsin), a serine protease inhibitor which is up-regulated in this study, may dampen the proteolytic effect of the neutrophil-specific elastase (Janciauskiene et al., 2018). However, in the complex interplay between proinflammatory and prothrombotic factors, the capillary integrity appears to be negatively influenced resulting in an increased transcapillary flux of filtrated plasma, accumulating in the alveoli and causing pulmonary edema. PAI-1 (SERPINE 1) transcripts, involved in fibrinolysis, were markedly up-regulated as a part of the inflammatory response (Figure 7A). This inhibition of the fibrinolysis was counter balanced by other molecules such as PLAT (tPA) (Figure 7A and Additional File 3: Table S1). However, the balance between formation of thrombi vs. fibrinolysis appears to be tilted toward thrombosis in capillaries in different organs which was observed in the lungs of patient 2 (Brusletto et al., 2017). Our results are in line with an endotoxemic animal study (Semeraro and Colucci, 2000-2013) showing increased levels of PAI-1 mRNA in multiple organs.

CCL2 (MCP-1) detected in the septic lung tissue samples (Figure 7A and Additional File 3: Table S1), is a key chemoattractant protein that regulates migration and infiltration of monocytes/macrophages, and is synthesized after activation of the TREM1, HMBG1 signaling, and LXR/RXR activation pathways. The transmembrane receptors IL1RL1 and ITGA5, previously found to be regulating chemokines as CCL2 and CCL3, are also implicated in the pathogenesis of lung injury (Oshikawa et al., 2002; Akhabir and Sandford, 2010; Sarangi et al., 2012; Hellerud et al., 2015; Hotchkiss et al., 2016). Several inhibitors/modulators of the proinflammatory response such as NFKB1A-, TNFAIP3-, TNFAIP6-, TNIP1-, and IRAK3 mRNA (Figure 7A and Additional File 3: Table S1) were found to be regulated in lungs of the studied patients, suggesting modulation of the macrophages in the lungs from pro- to an anti-inflammatory phenotype. These molecules are associated with poor outcome (Pino-Yanes et al., 2011; Mittal et al., 2016; Jimenez-Sousa et al., 2019).

The MAP2K3 gene in the p38 MAPK signaling pathway, was significantly up-regulated in the lung tissue samples (Figures 6, 7A). This pathway has widespread effects in the pathophysiology of multiorgan dysfunction in septic shock in lungs (Asaduzzaman et al., 2008) as well as myocardium in meningococcal septic shock (Pathan et al., 2011). Autophagy by binding of LPS to TLR4 and inducing MAPK/p38 signaling, is found to be a pathway which is protective against multiple organ injuries in a murine sepsis models by preventing apoptosis, maintaining a balance between the productions of pro- and anti-inflammatory cytokines, and preserving mitochondrial functions (Ho et al., 2016b).

Iron is an essential factor required for the N. meningitidis bacteria to colonize and cause disease in humans (Ali et al., 2017). Iron hemostasis in the lungs was affected via up-regulation of FTL in our study (Figure 7A and Additional File 3: Table S1), an iron storage protein as well as an important regulator which diminishes inflammation and increases the anti-inflammatory response (Fan et al., 2014). Furthermore, the up-regulated genes SOD2 and HAMP, also have roles against oxidative stress and inflammatory cytokines from incoming neutrophils, as well as to preventing iron from being requisitioned by invading bacteria (Schmidt, 2015).

The extent of infiltration of inflammatory cells varies among tissue samples, but overall we found the greatest influx of cells in the lungs (Figure 2 and Table 2). This is consistent with our gene expression results that showed extensively pro-inflammatory responses in the lung tissue samples. The results also support our previous data which demonstrated a higher number of N. meningitidis in lung tissue samples compared with the other tissues (Brusletto et al., 2017). Other studies based on FFPE tissue samples have confirmed attachment of meningococci to the endothelial cells in different FFPE tissue organs using immunohistochemistry and PCR assays (Guarner et al., 2004; Mairey et al., 2006). The adherence of monocytes, which have turned into inflammatory macrophages, and presence of neutrophils, are facilitated by up-regulation of various adhesion molecules (ICAM, VCAM, E-selectin) on the endothelial cells and circulating leukocytes (E-selectin) (Hellerud et al., 2015). We found that PAI-1 and MCP-1, were located to the epithelium (Figure 2). Our results are consistent with an E.coli induced sepsis model in baboons which detected neutrophil granulocytes and macrophages on the endothelial cells in the lungs (Tang et al., 2007). Concomitantly PAI-1 was up-regulated in lungs (Silasi-Mansat et al., 2010). CCR2, the receptor for MCP-1 (CCL2) has also been found up-regulated in macrophages in lungs from patients with sepsis-induced lung injury (An et al., 2009).

Transcriptional Profiles in Heart Tissue Samples

In industrialized countries heart failure combined with persistent hypotension are the primary causes of death in meningococcal disease. Acute heart failure has previously been attributed to inflammatory foci in the myocardium (Ferguson and Chapman, 1948; Dacosta et al., 1991; Neveling and Kaschula, 1993; Garcia et al., 1999). Troponin I, a specific marker of myocardial impairment and cell death, increased within 48 h of hospital admission of meningococcal shock patients and was related to a clinical severity score (Thiru et al., 2000). During the last 40 years, the circulatory collapse in fulminant meningococcal septicemia has been regarded as a consequence of declining vascular resistance and hypovolemia combined with an acute inflammation-induced cardiac failure (van Deuren et al., 2000; Stephens et al., 2007; Brandtzaeg and Van Deuren, 2012).

Myocardial depression is a well-recognized manifestation of organ dysfunction in sepsis. Excessive formation of nitric oxide (NO), reactive oxygen species (ROS) or nitrogen radicals, and transcriptional and metabolic changes have been proposed to explain the dysfunction (Martin et al., 2019) by either affecting a transient rise in cytosolic calcium (Ca2+) or a decrease in the cardiac contractile forces that defines septic cardiomyopathy (Martin et al., 2019). Endoplasmic reticulum (ER) stress due to activation of the unfolded protein response (UPR) via EIF-2 signaling pathway may also lead to myocardial apoptosis contributing to cardiac contractile dysfunction (Ceylan-Isik et al., 2010). The EIF signaling pathway, extensively affected in the septic heart tissue in our study (Figures 6, 7B) is known to moderate the synthesis of multiple molecules that are produced during a bacterial infection and to dampen ER stress (Nakayama et al., 2010; Shrestha et al., 2012). Several long non-coding RNAs (lncRNAs) (data not shown), were also significantly regulated in the septic heart samples. Such lncRNAs are in a rat model, suggested play an interaction and regulation role in the pathogenesis of sepsis-induced myocardial depression (Zhang T. N. et al., 2019).

The enzymes HSPA5 (BiP), HNRNPA1 and the transcription regulator PABPC1 previously found to be significantly expressed in response to ER stress (Dudek et al., 2009) and during gram-negative sepsis, may all induce autophagy (Xu et al., 2007; Ho et al., 2016b) and were substantially up-regulated in our meningococcal shock tissue samples (Figure 7B and Additional File 3: Table S2).

Both the TREM-1 and the HMBG1signaling pathways were up-regulated in the heart (Figures 6, 7B). They converge to amplify an inflammatory response by increasing the synthesis of pro-inflammatory mediators such as MCP-1, IL-8, IL1R1, and ICAM1. These molecules may also lead to binding of leukocytes to endothelial surfaces and disruption of the structures that maintain the integrity of the endothelium in a tissue. The transmembrane receptors SELE (E-selectin), ITGA5 and IL1R1 (ST2), all known (Hotchkiss et al., 2016; Hakanpaa et al., 2018; Szekely and Arbel, 2018) to recruit and enable leukocytes to roll along surfaces, were significantly up-regulated in our data set (Figure 7B and Additional File 3: Table S2). Studies of patients with meningococcal septic shock, has previously shown elevated levels of ICAM1 in the blood (Baines et al., 1999) while increased levels of IL1R1 correlates with poor prognosis in cardiovascular failure (Sabatine et al., 2008). Increased presence of macrophages and neutrophils in the heart tissue samples expressing IL-8, MCP-1 and PAI-1 (SERPINE1 gene) was confirmed by quantification of IL8 by Luminex and MCP-1 and PAI-1 by immunohistochemically methods (Figures 2, 8 and Table 2). In vitro studies of LPS-treated myocardium suggest that ICAM1 and VCAM1 contribute to myocardial dysfunction independent of neutrophil accumulation (Raeburn et al., 2002).

A dysregulation of the fibrinolytic system in the heart tissue samples due to significant up-regulation of SERPINE1 (PAI-1) was observed (Figures 2, 7B, 8 and Additional File 3: Table S2). No fibrin clots were observed histochemically. However, the presence of PAI-1 was found in cells in the heart tissue samples, suggesting that a reduced removal of fibrin may occur, leading to a deposition of fibrin clots in small blood vessels and inadequate tissue perfusion and organ failure (Brandtzaeg et al., 1990; Madoiwa et al., 2006). Previous studies indicate that elevated PAI-1 in plasma or serum is a significant predictor of disease severity and mortality, first discovered in plasma samples from meningococcal sepsis patients (Brandtzaeg et al., 1990; Kornelisse et al., 1996; Tipoe et al., 2018).

Damage-associated molecular pattern (DAMP) proteins such as High Mobility Group Box 1 (HMGB1), S100 proteins (S100A8/A9) and IL-1α, shown in our data to be significantly regulated, may all contribute to a host's defense by interacting through pattern recognition receptors (PRRs) such as RAGE and TLR4. HMGB1 has a pleiotropic role in inflammation and may be either secreted by activated macrophages and monocytes at a late pro-inflammatory stage, or released through cell necrosis and apoptosis (Wang et al., 1999; Bertheloot and Latz, 2017). HMBG1 has also an important role in the pathogenesis of cardiac dysfunction through increased ROS levels (Zhang et al., 2014).

The PI3K/Akt pathway (Figures 6, 7B), significantly up-regulated in the heart tissue samples, may negatively regulate NF-kB and probably limit pro-inflammatory and apoptotic events in monocytes/macrophages to protect the myocardium (Pourrajab et al., 2015). In a study with endotoxemic mice the P13K/AKT signaling pathway suppressed LPS–induced inflammation and coagulation (Schabbauer et al., 2004). In addition, YWHAE (14-3-3ε) was found up-regulated (Figure 7B and Additional File 3: Table S2). This is a protein that tightly regulates cellular and tissue homeostasis (Kosaka et al., 2012), and has been shown to have a role in inhibiting apoptosis in cardiomyocytes (Xing et al., 2000).

The LXR/RXR signaling pathway was significantly down-regulated in the heart tissue samples (Figures 4B, 5B). It involves the gene S100A8 that acts together with S100A9 to form the calprotectin protein, with both pro–and anti-inflammatory properties (Wang et al., 2018), as well as it exerts antimicrobial function by the ability to bind and control Zn2++ and Mn2++ required for bacterial growth (Damo et al., 2013). Calprotectin has a role in cardiovascular disease, and acts as an amplifier of the inflammatory response (Schiopu and Cotoi, 2013). The median calprotectin level in plasma from 13 patients in a meningococcal septic shock study was 20-fold higher than controls (Johne et al., 1997).

The PPAR signaling pathway was also significantly down-regulated in heart tissue samples. The nuclear receptors, particularly PPARs regulate cardiac fatty acid oxidation (FAO). In sepsis-mediated cardiac dysfunction a reduced energy production due in part to compromised FAO (Drosatos et al., 2011, 2013) and glucose catabolism (Tessier et al., 2003) will occur. In response to LPS and bacterial infections, down-regulation of PPARs and reduced FAO levels have been found in heart, liver and kidneys (Beigneux et al., 2000; Feingold et al., 2008; Drosatos et al., 2013). In children with septic shock decreased PPARα expression in whole blood correlates with severity (Standage et al., 2012).

Transcriptional Profiles in Kidney Tissue Samples

Severe renal failure sets in gradually during bacterial sepsis and significantly contributes to the mortality of late-phase sepsis (White et al., 2013; Hotchkiss et al., 2016). Immunohistochemical examination of postmortem kidney tissues from patients with severe sepsis, reveals heterogeneous and rather nonspecific findings (Takasu et al., 2013). The cellular and molecular pathways are not well-understood (Hato et al., 2019).

In meningococcal septic shock, activation of the coagulation system and inhibition of the fibrinolytic system are so pronounced that thrombosis of the glomerular capillaries appears to be the principal cause of the acute renal failure. Patient 1, 2 and 4 (kidney samples from patient 5 was not microscopically examined) had thrombotic lesions of the glomeruli leading to acute reduced glomerular filtration and necrosis of the proximal tubuli (Brusletto et al., 2017). Several serine proteases were significantly up-regulated in the kidneys (Figures 6, 7C and Additional File 3: Table S3). Most strikingly were SERPINE1 (PAI-1) inhibiting fibrionolysis (Malgorzewicz et al., 2013), and SERPINA3 (α1-antichymotrypsin) and SERPINA1 (α1-antitrypsin), both known to affect the renin-angiotensin system in various ways (Schmaier, 2002; Zhu et al., 2007).

Our results may differ from septic shock studies caused by other pathogens with less pronounced coagulopathy (Brandtzaeg, 2006). The microarray results are in line with other studies in which the kidneys are exposed to endoplasmic reticulum (ER) stress and unfolded protein response (UPR) during severe pathological and inflammatory conditions (Khan et al., 2015). In a gram-negative animal sepsis model, the EIF1/EIF2a transcription regulators were mediators of translation initiation block in late-phase sepsis when transcriptomic changes were examined at multiple time points (Hato et al., 2019). Recent research indicates that ER stress is a major factor in renal tubular cell apoptosis resulting from ischemic acute kidney injury (Xu et al., 2016). The hypotension combined with reduced cardiac output will have a dramatic negative effect on renal function because of inadequate perfusion (Dickhout et al., 2011).

Our data also show activation of the complement system via C3, a component at the convergence of all three activation pathways in the complement system, leading to lysis and opsonophagocytosis when C3 is cleaved to C3b. The complement system is massively up-regulated in blood in meningococcal septic shock patients and directly associated with a lethal outcome (Brandtzaeg et al., 1989b, 1996a).

The SOD2 transcript was up-regulated in the kidneys. SOD2 is one of the oxidative stress genes that may counteract the effects of incoming neutrophils in the tissue (Macdonald et al., 2003). Up-regulation of this gene was observed in the lung tissue in a E.coli sepsis baboon model (Zhu et al., 2007). Haptoglobin (HP) is another oxidative stress gene found up-regulated. This high-affinity hemoglobin-binding protein and antioxidant may prevent degradative enzymes from gaining access to hemoglobin and thereby prevent loss of iron through the kidneys (Zager et al., 2012).

Transcriptional Profiles in Liver Tissue Samples

The liver has critical roles during infections associated with removal of intravascular bacteria and LPS, i.e., bacteria scavenging, detoxifications by acyloxyacyl hydrolase (AOAH) and alkaline phosphatase and synthesizing acute-phase proteins (Poelstra et al., 1997; Shao et al., 2007; Deng et al., 2013). The response to the PAMPS is a highly regulated process involving several cells present in the liver such as hepatocytes, Kupffer cells, sinusoidal endothelial cells, primarily under the influence of cytokines such as TNF, IL-6 and IL-1β, (Nesseler et al., 2012). A combination of the cytokines are widely involved in acute phase protein production in the liver (Dhainaut et al., 2001).

The mRNA data in the liver tissue samples (Figures 6, 7D and Additional File 3: Table S4), showed a significant down-regulation of the LXR/RXR activation signaling pathway involving HDL, ApoA5 and LDL. These genes, involved in lipid uptake and cholesterol efflux, may be controlled by peroxisomes proliferation-activated receptors (PPARs) and liver X receptors (LXRs) (Zelcer and Tontonoz, 2006). Top upstream regulators in the liver tissue samples (Figure 4C) were TNF, IL1B and IL6. mRNA of AOAH and alkaline phosphatase were not significantly altered.

Transcriptional Profiles in Spleen Tissue Samples

The genes of the spleen involved in septic shock (Figures 6, 7E and Additional File 3: Table S5) were to our surprise altered less than the genes in the other organs examined. The numbers of N.meningitidis in spleen tissue were also lower than in the other septic tissue samples (Table 1). Oxidative phosphorylation and mitochondrial dysfunction were affected and this may be explained by multiple organ failure following sepsis, possibly representing an adaptive state during which the organs “shut down” their normal metabolic functions to protect themselves from the overwhelming and prolonged insult. The decrease in energy supply due to mitochondrial inhibition may also trigger a “hibernation like-state” (Levy, 2007; Azevedo, 2010). Our top-ten list from the spleen tissue samples also showed an up-regulation of several piRNAs, molecules found to be involved in gene silencing and regulatory mechanisms (Parhad and Theurkauf, 2019; Zhang X. et al., 2019). Is the low number of genes regulated an indication of a dysfunctional spleen? Clinical experience suggests that these patients are not more inclined to contract other types of life threatening invasive bacterial infections before or after contracting N.meningitidis leading to meningococcal septic shock. Could the low number of regulated genes in the spleen represent an evolutionary adaptation mechanism since a key function of the spleen is removal of live and disintegrated bacteria from the circulation?

Upstream Regulators

Prediction of top upstream transcription regulators in our tissue samples were found to be TNF, IL1B, IL6, IFNG, NFκB (complex), RICTOR, several miRNAs and CD3 (Figure 4C and Table 5), regulating signaling pathways in a sequential way. TNF, IL1B, IL6, IFNG, and NFκB (complex) are predicted to regulate cell survival signaling pathways, proliferation and metabolic processes. In addition, TNF can trigger programmed cell death (Varfolomeev and Vucic, 2018) (Figures 4C, 9A). The production of these upstream regulators is consistent with upstream regulators found in lung tissue samples in an E.coli sepsis baboon model (Zhu et al., 2007).

FIGURE 9
www.frontiersin.org

Figure 9. Predicted gene signaling network from the top up-regulated (A) and down-regulated (B) upstream regulators in FFPE tissue samples from patients with meningococcal septic shock vs. controls. The Z-score indicates predicted activation state of genes in the upstream regulator signaling network. Orange or lighter shades of orange indicate a positive Z-score and up-regulation of the upstream regulator. Blue color or lighter shades of blue indicate a negative Z-score and down-regulation of the upstream regulator. The transcripts in the signaling network are expressed as Fold Change (FC) values. Red or lighter shades of red indicates positive FC-values and up-regulation of transcripts, green color or lighter shades of green indicates negative FC-values and down-regulation of transcripts. Color gray indicates that a predicted activation state of a gene/transcript in the upstream regulator signaling network is not affected.

RICTOR, involved in regulation of cell growth and survival, may be inhibited due to overproduction of proinflammatory cytokines and cell survival (Cheng et al., 2016) (Figures 4C, 9B). Regulation of microRNA has been shown in several sepsis studies, reducing inflammatory cytokines and inhibition of cellular apoptosis (Ho et al., 2016a). Our data show several immature microRNAs to be down-regulated in our tissue samples. miR-6739-3p the most regulated in our data set, involved e.g., up-regulation of the chemokine CXCL1 (GRO-α) (Figures 4C, 9B), by recruiting and activating neutrophils for microbial killing at the tissue site (Sawant et al., 2016). miR-182-5p was found to be associated with down-regulation of multiple proteins in several tissue samples, however, most extensively in heart (Figures 4C, 9B). Regulation of miR-182-5p has been observed both in metabolic pathways, in cardiomyocytes and in sepsis (Vasilescu et al., 2009; Zhang et al., 2016, 2018).

Down-regulation of the upstream regulator CD3 was most strikingly observed in the lung tissue samples. CD3, associated with T-cell receptors that activate T-lymphocytes, has also been reported to be down-regulated in blood samples from other septic shock studies (Venet et al., 2005; Cazalis et al., 2014). However slightly higher numbers of CD3 positive immunostained cells were detected in both lungs and heart from meningococcal septic shock patients as compared to controls (Table 2).

Conclusions

Our study demonstrates that organismal death and multiple defense mechanisms, both pro-and anti-inflammatory, were heavily activated in tissue samples from our patients. Genes associated with several non-inflammatory pathways related to basic metabolism and energy production showed reduced expression. These observations possibly represent a complex preprogrammed response to the massive amounts of meningococci in order to balance the homeostasis in the inflamed organs.

Septic lungs showed the most extensive gene expression changes whereas tissues from the spleen revealed, unexpectedly, fewer regulated genes. The different organs showed a quite specific transcriptional signature. These signatures reflect the activation of thousands of specific genes in resident cells and of immune cells attracted to the different organs from the circulation. Collectively they represent the cellular mRNA “fingerprint” associated with a lethal outcome in meningococcal septic shock and multiple organ failure induced by the massive proliferation of N. meningitidis. The results suggest that the acute organ failure is a combined effect of multiple genes being up- or down-regulated in concert and not changes in one single pathway. Future studies will add details that may explain more clearly the contribution of protein-coding and non-coding RNA transcripts in the pathophysiology of meningococcal septic shock and possibly identify new intervention methods.

Data Availability Statement

The datasets supporting the conclusions of this article are available in the Gene Expression Omnibus (GEO) repository https://www.ncbi.nlm.nih.gov/geo/ under the identifier GSE141864 in accordance with minimum information about a microarray experiment (MIAME) standards.

Ethics Statement

The study was approved by the Regional Medical Ethical Committee of South East Norway (2011/1413C Translational research, meningococcal disease and 2011/753 Studies of invasive meningococcal and pneumococcal disease). The patients' samples were collected after informed consent from patient parents or relatives and according to the Helsinki declaration. The Director of Public Prosecutions approved the use of forensic material for this research.

Author Contributions

BB, RØ, PB, BH, JB, UG, OO, EL, and IG conceptualized and designed the study. EL and BH contributed with patient data and paraffin blocks from their hospital. BB, IG, EL, and RØ performed the laboratory experiments. BB, RØ, and PB performed the statistical analysis and drafted the manuscript. BB and OO performed the data analysis. BB, RØ, PB, JB, BH, UG, EL, and IG critically revised the manuscript. RØ, PB, and JB supervised the study. All authors read and approved the final manuscript.

Funding

This research was funded by South–Eastern Norway Regional Health Authority program, pro 2009062, meningococcal infection.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to the patients and their families that allowed us to use materials for research studies. Jack Levin at Department of Laboratory Medicine, School of Medicine, University of California San Francisco, USA. Åshild Vege at Section for Forensic Pediatric Pathology, Department of Forensic Sciences, Oslo University Hospital, Norway, Hege Ulland Dirdal at Division of Medical Service, Department of Pathology, Stavanger University Hospital, Norway. Anne-Marie Siebke Trøseid, Blood Cell Research Group, Section for Research, Development and Innovation, Department of Medical Biochemistry, Oslo University Hospital, Norway. Lilly Alice Steffensen, Blood Cell Research Group Department of Medical Biochemistry, Oslo University Hospital, Norway.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00042/full#supplementary-material

Abbreviations

ApoA5, Apolipoprotein A5; C3, Complement component 3; CD, cluster of differentiation; CREB, cAMP response element-binding protein; CSF, Cerebrospinal fluid; CXCL2, Chemokine (C-X-C motif) ligand 2 (alias MIP-2α); CXCL8, chemokine (C-X-C motif) ligand 8 (alias IL-8); CXCL10, (C-X-C) motif chemokine 10 (alias IP-10); DAMPs, Damage-associated molecular patterns (alias alarmin); DIC, disseminated intravascular coagulation; DYRK2, Dual specificity tyrosine-phosphorylation-regulated kinase 2; EIF2, Eukaryotic Initiation Factor 2; ER, endoplasmic reticulum; EU, Endotoxin unit; FF, fresh frozen; FFPE, formalin-fixed, paraffin-embedded; FOS, Fos proto-oncogene, AP-1 transcription factor subunit; FTL, Ferritin light chain; G-CSF, Granulocyte-colony stimulating factor (alias CFR 3); HAMP, Hepcidin; HDL, High-density lipoprotein; HMBG1, High mobility group box 1 protein; HNRNPA1, Heterogeneous nuclear ribonucleoprotein A1; HP, Haptoglobin; Hsp90, heat shock protein 90; HSP90AA1, Heat Shock Protein 90 Alpha Family Class A Member 1; HSP90AB1, Heat Shock Protein 90 Alpha Family Class B Member 1; HSP90B1, Heat shock protein 90kDa beta member 1; HSPA5, heat shock 70 kDa protein 5 (alias BiP); ICAM, intercellular adhesion molecules (alias CD 54); IFNG, Interferon gamma; IHC, Immunohistochemistry; IL, interleukin; IL1RL1, Interleukin 1 receptor-like 1 (ST2); IL1RN, interleukin-1 receptor antagonist (IL-1RA); IP-10, Interferon gamma-induced protein 10; IPA, Ingenuity Pathway Analysis; IRAK-3, Interleukin-1 receptor-associated kinase 3; ITGA5, Integrin alpha-5; ITGB1, Integrin beta-1; LAL, limulus amebocyte lysate; LDL, Low-density lipoprotein; lncRNA, long non-coding RNA; LPS, lipopolysaccharides; LXR, liver X receptor; MAP2K3, Dual specificity mitogen-activated protein kinase kinase 3; MCP-1, Monocyte chemoattractant protein-1 (alias CCL2); M-CSF, macrophage colony-stimulating factor (alias CSF1); MIF, Macrophage migration inhibitory factor; MIP, Macrophage Inflammatory Proteins; miR, microRNA; miRNA, microRNA; MOF, multiple organ failure; MPO, Myeloperoxidase; MSS, meningococcal septic shock; mTOR, mammalian target of rapamycin; MYCN, N-myc proto-oncogene protein; MyD88, Myeloid differentiation primary response 88; NFKB1A, NF-Kappa-B Inhibitor Alpha; NF-κB, nuclear factor kappa-light-chain-enhancer of activated B cells; NmDNA, N. meningitidis DNA; NO, nitric oxide; P13/AKT, phosphoinositide 3-kinase (PI3K)/AKT; PABPC1, Polyadenylate-binding protein 1; PAI-1, Plasminogen activator inhibitor-1 (alias SEPINE 1); PAMP, Pathogen-associated molecular pattern; PCA, Principal Component Analysis; piRNA, piwi-interacting RNA; PLAT, Tissue plasminogen activator (alias tPA); PPAR, peroxisome proliferator-activated receptors; PPARA, Peroxisome proliferator-activated receptor alpha; PPARGC1A, Peroxisome proliferator-activated receptor gamma coactivator 1-alpha; RAGE, receptor for advanced glycation endproducts; RANTES, regulated on activation, normal T cell expressed and secreted (alias CCL5); RICTOR, Rapamycin-insensitive companion of mammalian target of rapamycin; RIN, RNA integrity number; ROS, reactive oxygen species; RPL, L ribosomal proteins; RXR, retinoid X receptor; RXRA, Retinoid X receptor alpha; S100A8, calcium-binding protein A8 (alias Calgranulin A); S100A9, calcium-binding protein A9; SELE, E-selectin, SERPINA1, Serpin Family A Member 1 (alias Alpha-1-antitrypsin); SERPINA3, Serpin Family A Member 3 (alias Alpha 1-antichymotrypsin); SMD, systemic meningococcal disease; SOD-2, Superoxide dismutase 2, mitochondrial; STAT3, Signal transducer and activator of transcription 3; TLR4, Toll-like receptor 4; TLRs, Toll-like receptor family; TNF, tumor necrosis factor; TNFAIP3, Tumor necrosis factor, alpha-induced protein 3; TNFAIP6, Tumor Necrosis Factor-Inducible Gene 6 Protein (alias TSG-6); TNIP1, TNFAIP3-interacting protein 1; TREM1, Triggering receptor expressed on myeloid cells 1; UPR, unfolded protein response; VCAM-1, vascular cell adhesion molecule 1; YWHAE, Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Epsilon (alias 14-3-3 protein epsilon).

References

Abdueva, D., Wing, M., Schaub, B., Triche, T., and Davicioni, E. (2010). Quantitative expression profiling in formalin-fixed paraffin-embedded samples by affymetrix microarrays. J. Mol. Diagn. 12, 409–417. doi: 10.2353/jmoldx.2010.090155

PubMed Abstract | CrossRef Full Text | Google Scholar

Akhabir, L., and Sandford, A. (2010). Genetics of interleukin 1 receptor-like 1 in immune and inflammatory diseases. Curr. Genomics 11, 591–606. doi: 10.2174/138920210793360907

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, M. K., Kim, R. Y., Karim, R., Mayall, J. R., Martin, K. L., Shahandeh, A., et al. (2017). Role of iron in the pathogenesis of respiratory disease. Int. J. Biochem. Cell Biol. 88, 181–195. doi: 10.1016/j.biocel.2017.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

An, J. L., Ishida, Y., Kimura, A., Tsokos, M., and Kondo, T. (2009). Immunohistochemical detection of CCR2 and CX3CR1 in sepsis-induced lung injury. Forensic Sci. Int. 192, e21–25. doi: 10.1016/j.forsciint.2009.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Asaduzzaman, M., Wang, Y., and Thorlacius, H. (2008). Critical role of p38 mitogen-activated protein kinase signaling in septic lung injury. Crit. Care Med. 36, 482–488. doi: 10.1097/01.CCM.0B013E31816204FA

PubMed Abstract | CrossRef Full Text | Google Scholar

Azevedo, L. C. (2010). Mitochondrial dysfunction during sepsis. Endocr. Metab. Immune Disord. Drug Targets 10, 214–223. doi: 10.2174/187153010791936946

PubMed Abstract | CrossRef Full Text | Google Scholar

Baines, P. B., Marzouk, O., Thomson, A. P., Sills, J. A., Riordan, F. A., and Hart, C. A. (1999). Endothelial cell adhesion molecules in meningococcal disease. Arch. Dis. Child. 80, 74–76 doi: 10.1136/adc.80.1.74

PubMed Abstract | CrossRef Full Text | Google Scholar

Balmer, M. L., and Hess, C. (2017). Starving for survival-how catabolic metabolism fuels immune function. Curr. Opin. Immunol. 46, 8–13. doi: 10.1016/j.coi.2017.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Barquet, N., Domingo, P., Cayla, J. A., Gonzalez, J., Rodrigo, C., Fernandez-Viladrich, P., et al. (1999). Meningococcal disease in a large urban population (Barcelona, 1987-1992): predictors of dismal prognosis. Barcelona Meningococcal Disease Surveillance Group. Arch. Intern. Med. 159, 2329–2340. doi: 10.1001/archinte.159.19.2329

PubMed Abstract | CrossRef Full Text | Google Scholar

Beigneux, A. P., Moser, A. H., Shigenaga, J. K., Grunfeld, C., and Feingold, K. R. (2000). The acute phase response is associated with retinoid X receptor repression in rodent liver. J. Biol. Chem. 275, 16390–16399. doi: 10.1074/jbc.M000953200

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Bertheloot, D., and Latz, E. (2017). HMGB1, IL-1alpha, IL-33 and S100 proteins: dual-function alarmins. Cell. Mol. Immunol. 14, 43–64. doi: 10.1038/cmi.2016.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Boeckx, C., Wouters, A., Pauwels, B., Deschoolmeester, V., Specenier, P., Lukaszuk, K., et al. (2011). Expression analysis on archival material: comparison of 5 commercially available RNA isolation kits for FFPE material. Diagn. Mol. Pathol. 20, 203–211. doi: 10.1097/PDM.0b013e3182230937

PubMed Abstract | CrossRef Full Text | Google Scholar

Boucek, M. M., Boerth, R. C., Artman, M., Graham, T. P. Jr., and Boucek, R. J. Jr. (1984). Myocardial dysfunction in children with acute meningococcemia. J. Pediatr. 105, 538–542. doi: 10.1016/S0022-3476(84)80416-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P. (2006). “Pathogenesis and pathophysiology of invasive meningococcal disease,” in Handbook of Meningococcal Disease: Infection Biology, Vaccination, Clinical Management, eds M. Frosch and M. C. J. Maiden (Weinheim: Wiley-VCH Verlag GmbH & Co.), 427–480. doi: 10.1002/3527608508.ch21

CrossRef Full Text | Google Scholar

Brandtzaeg, P., Bjerre, A., Ovstebo, R., Brusletto, B., Joo, G. B., and Kierulf, P. (2001). Neisseria meningitidis lipopolysaccharides in human pathology. J. Endotoxin Res. 7, 401–420. doi: 10.1179/096805101101533016

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Hogasen, K., Kierulf, P., and Mollnes, T. E. (1996a). The excessive complement activation in fulminant meningococcal septicemia is predominantly caused by alternative pathway activation. J. Infect. Dis. 173, 647–655. doi: 10.1093/infdis/173.3.647

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Joo, G. B., Brusletto, B., and Kierulf, P. (1990). Plasminogen activator inhibitor 1 and 2, alpha-2-antiplasmin, plasminogen, and endotoxin levels in systemic meningococcal disease. Thromb. Res. 57, 271–278. doi: 10.1016/0049-3848(90)90326-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., and Kierulf, P. (1992). “Endotoxin and meningococcemia. Intravascular inflammation induced by native endotoxin in man,” in Bacterial Endotoxic Lipopolysaccharides: 2: Immunopharmacology and Pathophysiology, eds. D. C. Morrison and J. L. Ryan (Boca Raton: CRC Press, 327–346.

Google Scholar

Brandtzaeg, P., Kierulf, P., Gaustad, P., Skulberg, A., Bruun, J. N., Halvorsen, S., et al. (1989a). Plasma endotoxin as a predictor of multiple organ failure and death in systemic meningococcal disease. J. Infect. Dis. 159, 195–204. doi: 10.1093/infdis/159.2.195

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Mollnes, T. E., and Kierulf, P. (1989b). Complement activation and endotoxin levels in systemic meningococcal disease. J. Infect. Dis. 160, 58–65. doi: 10.1093/infdis/160.1.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Osnes, L., Ovstebo, R., Joo, G. B., Westvik, A. B., and Kierulf, P. (1996b). Net inflammatory capacity of human septic shock plasma evaluated by a monocyte-based target cell assay: identification of interleukin-10 as a major functional deactivator of human monocytes. J. Exp. Med. 184, 51–60. doi: 10.1084/jem.184.1.51

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Ovsteboo, R., and Kierulf, P. (1992). Compartmentalization of lipopolysaccharide production correlates with clinical presentation in meningococcal disease. J. Infect. Dis. 166, 650–652. doi: 10.1093/infdis/166.3.650

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., Sandset, P. M., Joo, G. B., Ovstebo, R., Abildgaard, U., and Kierulf, P. (1989c). The quantitative association of plasma endotoxin, antithrombin, protein C, extrinsic pathway inhibitor and fibrinopeptide A in systemic meningococcal disease. Thromb. Res. 55, 459–470. doi: 10.1016/0049-3848(89)90054-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandtzaeg, P., and van Deuren, M. (2012). Classification and pathogenesis of meningococcal infections. Methods Mol. Biol. 799, 21–35. doi: 10.1007/978-1-61779-346-2_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, R., Woods, C. W., Benjamin, D. K. Jr., and Rosenstein, N. E. (2006). Increased case-fatality rate associated with outbreaks of Neisseria meningitidis infection, compared with sporadic meningococcal disease, in the United States, 1994-2002. Clin. Infect. Dis. 43, 49–54. doi: 10.1086/504804

PubMed Abstract | CrossRef Full Text | Google Scholar

Brusletto, B. S., Hellerud, B. C., Loberg, E. M., Goverud, I. L., Vege, A., Berg, J. P., et al. (2017). Traceability and distribution of Neisseria meningitidis DNA in archived post mortem tissue samples from patients with systemic meningococcal disease. BMC Clin. Pathol. 17:10. doi: 10.1186/s12907-017-0049-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Carithers, L. J., Ardlie, K., Barcus, M., Branton, P. A., Britton, A., Buia, S. A., et al. (2015). A novel approach to high-quality postmortem tissue procurement: the GTEx project. Biopreserv. Biobank. 13, 311–319. doi: 10.1089/bio.2015.0032

PubMed Abstract | CrossRef Full Text | Google Scholar

Castrillo, A., and Tontonoz, P. (2004). Nuclear receptors in macrophage biology: at the crossroads of lipid metabolism and inflammation. Annu. Rev. Cell Dev. Biol. 20, 455–480. doi: 10.1146/annurev.cellbio.20.012103.134432

PubMed Abstract | CrossRef Full Text | Google Scholar

Cazalis, M. A., Lepape, A., Venet, F., Frager, F., Mougin, B., Vallin, H., et al. (2014). Early and dynamic changes in gene expression in septic shock patients: a genome-wide approach. Intensive Care Med. Exp. 2:20. doi: 10.1186/s40635-014-0020-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceylan-Isik, A. F., Zhao, P., Zhang, B., Xiao, X., Su, G., and Ren, J. (2010). Cardiac overexpression of metallothionein rescues cardiac contractile dysfunction and endoplasmic reticulum stress but not autophagy in sepsis. J. Mol. Cell. Cardiol. 48, 367–378. doi: 10.1016/j.yjmcc.2009.11.003

CrossRef Full Text | Google Scholar

Cheng, S. C., Scicluna, B. P., Arts, R. J., Gresnigt, M. S., Lachmandas, E., Giamarellos-Bourboulis, E. J., et al. (2016). Broad defects in the energy metabolism of leukocytes underlie immunoparalysis in sepsis. Nat. Immunol. 17, 406–413. doi: 10.1038/ni.3398

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, J. Y., Braunschweig, T., Williams, R., Guerrero, N., Hoffmann, K. M., Kwon, M., et al. (2008). Factors in tissue handling and processing that impact RNA obtained from formalin-fixed, paraffin-embedded tissue. J. Histochem. Cytochem. 56, 1033–1042. doi: 10.1369/jhc.2008.951863

PubMed Abstract | CrossRef Full Text | Google Scholar

Dacosta, D. F., Neville, E., and Raman, G. V. (1991). Association of myocarditis and acute renal failure complicating meningococcal septicaemia. J. Infect. 22, 259–261. doi: 10.1016/S0163-4453(05)80008-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Damo, S. M., Kehl-Fie, T. E., Sugitani, N., Holt, M. E., Rathi, S., Murphy, W. J., et al. (2013). Molecular basis for manganese sequestration by calprotectin and roles in the innate immune response to invading bacterial pathogens. Proc. Natl. Acad. Sci. U.S.A. 110, 3841–3846. doi: 10.1073/pnas.1220341110

PubMed Abstract | CrossRef Full Text | Google Scholar

de Greeff, S. C., De Melker, H. E., Schouls, L. M., Spanjaard, L., and Van Deuren, M. (2008). Pre-admission clinical course of meningococcal disease and opportunities for the earlier start of appropriate intervention: a prospective epidemiological study on 752 patients in the Netherlands, 2003-2005. Eur. J. Clin. Microbiol. Infect. Dis. 27, 985–992. doi: 10.1007/s10096-008-0535-1

PubMed Abstract | CrossRef Full Text | Google Scholar

de Kleijn, E. D., Hazelzet, J. A., Kornelisse, R. F., and De Groot, R. (1998). Pathophysiology of meningococcal sepsis in children. Eur. J. Pediatr. 157, 869–880. doi: 10.1007/s004310050958

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, M., Scott, M. J., Loughran, P., Gibson, G., Sodhi, C., Watkins, S., et al. (2013). Lipopolysaccharide clearance, bacterial clearance, and systemic inflammatory responses are regulated by cell type-specific functions of TLR4 during sepsis. J. Immunol. 190, 5152–5160. doi: 10.4049/jimmunol.1300496

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhainaut, J. F., Marin, N., Mignon, A., and Vinsonneau, C. (2001). Hepatic response to sepsis: interaction between coagulation and inflammatory processes. Crit. Care Med. 29, S42–47. doi: 10.1097/00003246-200107001-00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickhout, J. G., Carlisle, R. E., and Austin, R. C. (2011). Interrelationship between cardiac hypertrophy, heart failure, and chronic kidney disease: endoplasmic reticulum stress as a mediator of pathogenesis. Circ. Res. 108, 629–642. doi: 10.1161/CIRCRESAHA.110.226803

PubMed Abstract | CrossRef Full Text | Google Scholar

Drosatos, K., Drosatos-Tampakaki, Z., Khan, R., Homma, S., Schulze, P. C., Zannis, V. I., et al. (2011). Inhibition of c-Jun-N-terminal kinase increases cardiac peroxisome proliferator-activated receptor alpha expression and fatty acid oxidation and prevents lipopolysaccharide-induced heart dysfunction. J. Biol. Chem. 286, 36331–36339. doi: 10.1074/jbc.M111.272146

PubMed Abstract | CrossRef Full Text | Google Scholar

Drosatos, K., Khan, R. S., Trent, C. M., Jiang, H., Son, N. H., Blaner, W. S., et al. (2013). Peroxisome proliferator-activated receptor-gamma activation prevents sepsis-related cardiac dysfunction and mortality in mice. Circ. Heart Fail. 6, 550–562. doi: 10.1161/CIRCHEARTFAILURE.112.000177

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudek, J., Benedix, J., Cappel, S., Greiner, M., Jalal, C., Muller, L., et al. (2009). Functions and pathologies of BiP and its interaction partners. Cell. Mol. Life Sci. 66, 1556–1569. doi: 10.1007/s00018-009-8745-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Zhang, J., Cai, L., Wang, S., Liu, C., Zhang, Y., et al. (2014). The effect of anti-inflammatory properties of ferritin light chain on lipopolysaccharide-induced inflammatory response in murine macrophages. Biochim. Biophys. Acta 1843, 2775–2783. doi: 10.1016/j.bbamcr.2014.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Farragher, S. M., Tanney, A., Kennedy, R. D., and Paul Harkin, D. (2008). RNA expression analysis from formalin fixed paraffin embedded tissues. Histochem. Cell Biol. 130, 435–445. doi: 10.1007/s00418-008-0479-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Faust, S. N., Levin, M., Harrison, O. B., Goldin, R. D., Lockhart, M. S., Kondaveeti, S., et al. (2001). Dysfunction of endothelial protein C activation in severe meningococcal sepsis. N. Engl. J. Med. 345, 408–416. doi: 10.1056/NEJM200108093450603

PubMed Abstract | CrossRef Full Text | Google Scholar

Feingold, K. R., Wang, Y., Moser, A., Shigenaga, J. K., and Grunfeld, C. (2008). LPS decreases fatty acid oxidation and nuclear hormone receptors in the kidney. J. Lipid Res. 49, 2179–2187. doi: 10.1194/jlr.M800233-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferguson, J. H., and Chapman, O. D. (1948). Fulminating meningococcic infections and the so-called Waterhouse-Friderichsen syndrome. Am. J. Pathol. 24, 763–795.

PubMed Abstract | Google Scholar

Ferruelo, A., El-Assar, M., Lorente, J. A., Nin, N., Penuelas, O., Fernandez-Segoviano, P., et al. (2011). Transcriptional profiling and genotyping of degraded nucleic acids from autopsy tissue samples after prolonged formalin fixation times. Int. J. Clin. Exp. Pathol. 4, 156–161.

PubMed Abstract | Google Scholar

Garcia, N. S., Castelo, J. S., Ramos, V., Rezende, G. S., and Pereira, F. E. (1999). Frequency of myocarditis in cases of fatal meningococcal infection in children: observations on 31 cases studied at autopsy. Rev. Soc. Bras. Med. Trop. 32, 517–522. doi: 10.1590/S0037-86821999000500008

PubMed Abstract | CrossRef Full Text | Google Scholar

Gedde-Dahl, T. W., Hoiby, E. A., Schillinger, A., Lystad, A., and Bovre, K. (1983). An epidemiological, clinical and microbiological follow-up study of incident meningococcal disease cases in Norway, winter 1981-1982. Material and epidemiology in the MenOPP project. NIPH Ann. 6, 155–168.

PubMed Abstract | Google Scholar

Gopinathan, U., Ovstebo, R., Olstad, O. K., Brusletto, B., Dalsbotten Aass, H. C., Kierulf, P., et al. (2012). Global effect of interleukin-10 on the transcriptional profile induced by Neisseria meningitidis in human monocytes. Infect. Immun. 80, 4046–4054. doi: 10.1128/IAI.00386-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Guarner, J., Greer, P. W., Whitney, A., Shieh, W. J., Fischer, M., White, E. H., et al. (2004). Pathogenesis and diagnosis of human meningococcal disease using immunohistochemical and PCR assays. Am. J. Clin. Pathol. 122, 754–764. doi: 10.1309/3489075U03LMK9AE

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, S., Halushka, M. K., Hilton, G. M., and Arking, D. E. (2012). Postmortem cardiac tissue maintains gene expression profile even after late harvesting. BMC Genomics 13:26. doi: 10.1186/1471-2164-13-26

PubMed Abstract | CrossRef Full Text | Google Scholar

Hackett, S. J., Guiver, M., Marsh, J., Sills, J. A., Thomson, A. P., Kaczmarski, E. B., et al. (2002). Meningococcal bacterial DNA load at presentation correlates with disease severity. Arch. Dis. Child. 86, 44–46. doi: 10.1136/adc.86.1.44

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagmolen of ten Have, W., Wiegman, A., Van Den Hoek, G. J., Vreede, W. B., and Derkx, H. H. (2000). Life-threatening heart failure in meningococcal septic shock in children: non-invasive measurement of cardiac parameters is of important prognostic value. Eur. J. Pediatr. 159, 277–282. doi: 10.1007/s004310050070

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakanpaa, L., Kiss, E. A., Jacquemet, G., Miinalainen, I., Lerche, M., Guzman, C., et al. (2018). Targeting beta1-integrin inhibits vascular leakage in endotoxemia. Proc. Natl. Acad. Sci. U.S.A. 115, E6467–e6476. doi: 10.1073/pnas.1722317115

PubMed Abstract | CrossRef Full Text | Google Scholar

Halstensen, A., Pedersen, S. H., Haneberg, B., Bjorvatn, B., and Solberg, C. O. (1987). Case fatality of meningococcal disease in western Norway. Scand. J. Infect. Dis. 19, 35–42. doi: 10.3109/00365548709032375

PubMed Abstract | CrossRef Full Text | Google Scholar

Hato, T., Maier, B., Syed, F., Myslinski, J., Zollman, A., Plotkin, Z., et al. (2019). Bacterial sepsis triggers an antiviral response that causes translation shutdown. J. Clin. Invest. 129, 296–309. doi: 10.1172/JCI123284

PubMed Abstract | CrossRef Full Text | Google Scholar

Hazelzet, J. A., Risseeuw-Appel, I. M., Kornelisse, R. F., Hop, W. C., Dekker, I., Joosten, K. F., et al. (1996). Age-related differences in outcome and severity of DIC in children with septic shock and purpura. Thromb. Haemost. 76, 932–938. doi: 10.1055/s-0038-1650688

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellerud, B. C., Nielsen, E. W., Thorgersen, E. B., Lindstad, J. K., Pharo, A., Tonnessen, T. I., et al. (2010). Dissecting the effects of lipopolysaccharides from nonlipopolysaccharide molecules in experimental porcine meningococcal sepsis. Crit. Care Med. 38, 1467–1474. doi: 10.1097/CCM.0b013e3181de8c94

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellerud, B. C., Olstad, O. K., Nielsen, E. W., Troseid, A. M., Skadberg, O., Thorgersen, E. B., et al. (2015). Massive organ inflammation in experimental and in clinical meningococcal septic shock. Shock 44, 458–469. doi: 10.1097/SHK.0000000000000441

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellum, M., Ovstebo, R., Brusletto, B. S., Berg, J. P., Brandtzaeg, P., and Henriksson, C. E. (2014). Microparticle-associated tissue factor activity correlates with plasma levels of bacterial lipopolysaccharides in meningococcal septic shock. Thromb. Res. 133, 507–514. doi: 10.1016/j.thromres.2013.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, J., Chan, H., Wong, S. H., Wang, M. H., Yu, J., Xiao, Z., et al. (2016a). The involvement of regulatory non-coding RNAs in sepsis: a systematic review. Crit. Care 20:383. doi: 10.1186/s13054-016-1555-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, J., Yu, J., Wong, S. H., Zhang, L., Liu, X., Wong, W. T., et al. (2016b). Autophagy in sepsis: degradation into exhaustion? Autophagy 12, 1073–1082. doi: 10.1080/15548627.2016.1179410

PubMed Abstract | CrossRef Full Text | Google Scholar

Hotchkiss, R. S., Moldawer, L. L., Opal, S. M., Reinhart, K., Turnbull, I. R., and Vincent, J. L. (2016). Sepsis and septic shock. Nat. Rev. Dis. Primers 2:16045. doi: 10.1038/nrdp.2016.45

PubMed Abstract | CrossRef Full Text | Google Scholar

Janciauskiene, S., Wrenger, S., Immenschuh, S., Olejnicka, B., Greulich, T., Welte, T., et al. (2018). The multifaceted effects of Alpha1-antitrypsin on neutrophil functions. Front. Pharmacol. 9:341. doi: 10.3389/fphar.2018.00341

PubMed Abstract | CrossRef Full Text | Google Scholar

Jimenez-Sousa, M. A., Fadrique, A., Liu, P., Fernandez-Rodriguez, A., Lorenzo-Lopez, M., Gomez-Sanchez, E., et al. (2019). TNFAIP3, TNIP1, and MyD88 polymorphisms predict septic-shock-related death in patients who underwent major surgery. J. Clin. Med. 8:E283. doi: 10.3390/jcm8030283

PubMed Abstract | CrossRef Full Text | Google Scholar

Johne, B., Fagerhol, M. K., Lyberg, T., Prydz, H., Brandtzaeg, P., Naess-Andresen, C. F., et al. (1997). Functional and clinical aspects of the myelomonocyte protein calprotectin. MP. Mol. Pathol. 50, 113–123. doi: 10.1136/mp.50.3.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, M. M., Yang, W. L., and Wang, P. (2015). Endoplasmic reticulum stress in sepsis. Shock 44, 294–304. doi: 10.1097/SHK.0000000000000425

PubMed Abstract | CrossRef Full Text | Google Scholar

Kidani, Y., and Bensinger, S. J. (2012). Liver X receptor and peroxisome proliferator-activated receptor as integrators of lipid homeostasis and immunity. Immunol. Rev. 249, 72–83. doi: 10.1111/j.1600-065X.2012.01153.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobsar, A., Siauw, C., Gambaryan, S., Hebling, S., Speer, C., Schubert-Unkmeir, A., et al. (2011). Neisseria meningitidis induces platelet inhibition and increases vascular endothelial permeability via nitric oxide regulated pathways. Thromb. Haemost. 106, 1127–1138. doi: 10.1160/TH11-07-0491

PubMed Abstract | CrossRef Full Text | Google Scholar

Kornelisse, R. F., Hazelzet, J. A., Savelkoul, H. F., Hop, W. C., Suur, M. H., Borsboom, A. N., et al. (1996). The relationship between plasminogen activator inhibitor-1 and proinflammatory and counterinflammatory mediators in children with meningococcal septic shock. J. Infect. Dis. 173, 1148–1156. doi: 10.1093/infdis/173.5.1148

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosaka, Y., Cieslik, K. A., Li, L., Lezin, G., Maguire, C. T., Saijoh, Y., et al. (2012). 14-3-3ε plays a role in cardiac ventricular compaction by regulating the cardiomyocyte cell cycle. Mol. Cell. Biol. 32, 5089–5102. doi: 10.1128/MCB.00829-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Krafft, A. E., Duncan, B. W., Bijwaard, K. E., Taubenberger, J. K., and Lichy, J. H. (1997). Optimization of the isolation and amplification of RNA from formalin-fixed, paraffin-embedded tissue: the armed forces institute of pathology experience and literature review. Mol. Diagn. 2, 217–230. doi: 10.1016/S1084-8592(97)80032-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehmann, A. K., Halstensen, A., Sornes, S., Rokke, O., and Waage, A. (1995). High levels of interleukin 10 in serum are associated with fatality in meningococcal disease. Infect. Immun. 63, 2109–2112. doi: 10.1128/IAI.63.6.2109-2112.1995

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, R. J. (2007). Mitochondrial dysfunction, bioenergetic impairment, and metabolic down-regulation in sepsis. Shock 28, 24–28. doi: 10.1097/01.shk.0000235089.30550.2d

PubMed Abstract | CrossRef Full Text | Google Scholar

Linton, K., Hey, Y., Dibben, S., Miller, C., Freemont, A., Radford, J., et al. (2009). Methods comparison for high-resolution transcriptional analysis of archival material on Affymetrix Plus 2.0 and Exon 1.0 microarrays. BioTechniques 47, 587–596. doi: 10.2144/000113169

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

CrossRef Full Text | Google Scholar

Ludyga, N., Grunwald, B., Azimzadeh, O., Englert, S., Hofler, H., Tapio, S., et al. (2012). Nucleic acids from long-term preserved FFPE tissues are suitable for downstream analyses. Virchows Arch. 460, 131–140. doi: 10.1007/s00428-011-1184-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Macdonald, J., Galley, H. F., and Webster, N. R. (2003). Oxidative stress and gene expression in sepsis. Br. J. Anaesth. 90, 221–232. doi: 10.1093/bja/aeg034

PubMed Abstract | CrossRef Full Text | Google Scholar

Madoiwa, S., Nunomiya, S., Ono, T., Shintani, Y., Ohmori, T., Mimuro, J., et al. (2006). Plasminogen activator inhibitor 1 promotes a poor prognosis in sepsis-induced disseminated intravascular coagulation. Int. J. Hematol. 84, 398–405. doi: 10.1532/IJH97.05190

PubMed Abstract | CrossRef Full Text | Google Scholar

Mairey, E., Genovesio, A., Donnadieu, E., Bernard, C., Jaubert, F., Pinard, E., et al. (2006). Cerebral microcirculation shear stress levels determine Neisseria meningitidis attachment sites along the blood-brain barrier. J. Exp. Med. 203, 1939–1950. doi: 10.1084/jem.20060482

PubMed Abstract | CrossRef Full Text | Google Scholar

Malgorzewicz, S., Skrzypczak-Jankun, E., and Jankun, J. (2013). Plasminogen activator inhibitor-1 in kidney pathology (Review). Int. J. Mol. Med. 31, 503–510. doi: 10.3892/ijmm.2013.1234

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, L., Derwall, M., Al Zoubi, S., Zechendorf, E., Reuter, D. A., Thiemermann, C., et al. (2019). The septic heart: current understanding of molecular mechanisms and clinical implications. Chest 155, 427–437. doi: 10.1016/j.chest.2018.08.1037

PubMed Abstract | CrossRef Full Text | Google Scholar

Maslove, D. M., and Wong, H. R. (2014). Gene expression profiling in sepsis: timing, tissue, and translational considerations. Trends Mol. Med. 20, 204–213. doi: 10.1016/j.molmed.2014.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Melican, K., Michea Veloso, P., Martin, T., Bruneval, P., and Dumenil, G. (2013). Adhesion of Neisseria meningitidis to dermal vessels leads to local vascular damage and purpura in a humanized mouse model. PLoS Pathog. 9:e1003139. doi: 10.1371/journal.ppat.1003139

PubMed Abstract | CrossRef Full Text | Google Scholar

Mercier, J. C., Beaufils, F., Hartmann, J. F., and Azema, D. (1988). Hemodynamic patterns of meningococcal shock in children. Crit. Care Med. 16, 27–33. doi: 10.1097/00003246-198801000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, M., Tiruppathi, C., Nepal, S., Zhao, Y. Y., Grzych, D., Soni, D., et al. (2016). TNFalpha-stimulated gene-6 (TSG6) activates macrophage phenotype transition to prevent inflammatory lung injury. Proc. Natl. Acad. Sci. U.S.A. 113, E8151–E8158. doi: 10.1073/pnas.1614935113

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittempergher, L., De Ronde, J. J., Nieuwland, M., Kerkhoven, R. M., Simon, I., Rutgers, E. J., et al. (2011). Gene expression profiles from formalin fixed paraffin embedded breast cancer tissue are largely comparable to fresh frozen matched tissue. PLoS ONE 6:e17163. doi: 10.1371/journal.pone.0017163

PubMed Abstract | CrossRef Full Text | Google Scholar

Moller, A. S., Bjerre, A., Brusletto, B., Joo, G. B., Brandtzaeg, P., and Kierulf, P. (2005). Chemokine patterns in meningococcal disease. J. Infect. Dis. 191, 768–775. doi: 10.1086/427514

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakayama, Y., Endo, M., Tsukano, H., Mori, M., Oike, Y., and Gotoh, T. (2010). Molecular mechanisms of the LPS-induced non-apoptotic ER stress-CHOP pathway. J. Biochem. 147, 471–483. doi: 10.1093/jb/mvp189

PubMed Abstract | CrossRef Full Text | Google Scholar

Nesseler, N., Launey, Y., Aninat, C., Morel, F., Malledant, Y., and Seguin, P. (2012). Clinical review: the liver in sepsis. Crit Care 16:235. doi: 10.1186/cc11381

PubMed Abstract | CrossRef Full Text | Google Scholar

Neveling, U., and Kaschula, R. O. (1993). Fatal meningococcal disease in childhood: an autopsy study of 86 cases. Ann. Trop. Paediatr. 13, 147–152. doi: 10.1080/02724936.1993.11747638

PubMed Abstract | CrossRef Full Text | Google Scholar

Oshikawa, K., Yanagisawa, K., Tominaga, S., and Sugiyama, Y. (2002). ST2 protein induced by inflammatory stimuli can modulate acute lung inflammation. Biochem. Biophys. Res. Commun. 299, 18–24. doi: 10.1016/S0006-291X(02)02578-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Osterud, B., and Flaegstad, T. (1983). Increased tissue thromboplastin activity in monocytes of patients with meningococcal infection: related to an unfavourable prognosis. Thromb. Haemost. 49, 5–7. doi: 10.1055/s-0038-1657303

PubMed Abstract | CrossRef Full Text | Google Scholar

Ovstebo, R., Aass, H. C., Haug, K. B., Troseid, A. M., Gopinathan, U., Kierulf, P., et al. (2012). LPS from Neisseria meningitidis is crucial for inducing monocyte- and microparticle-associated tissue factor activity but not for tissue factor expression. Innate Immun. 18, 580–591. doi: 10.1177/1753425911428230

CrossRef Full Text | Google Scholar

Ovstebo, R., Brandtzaeg, P., Brusletto, B., Haug, K. B., Lande, K., Hoiby, E. A., et al. (2004). Use of robotized DNA isolation and real-time PCR to quantify and identify close correlation between levels of Neisseria meningitidis DNA and lipopolysaccharides in plasma and cerebrospinal fluid from patients with systemic meningococcal disease. J. Clin. Microbiol. 42, 2980–2987. doi: 10.1128/JCM.42.7.2980-2987.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ovstebo, R., Hellum, M., Aass, H. C., Troseid, A. M., Brandtzaeg, P., Mollnes, T. E., et al. (2014). Microparticle-associated tissue factor activity is reduced by inhibition of the complement protein 5 in Neisseria meningitidis-exposed whole blood. Innate Immun. 20, 552–560. doi: 10.1177/1753425913502099

PubMed Abstract | CrossRef Full Text | Google Scholar

Parhad, S. S., and Theurkauf, W. E. (2019). Rapid evolution and conserved function of the piRNA pathway. Open Biol. 9:180181. doi: 10.1098/rsob.180181

PubMed Abstract | CrossRef Full Text | Google Scholar

Pascual, G., Fong, A. L., Ogawa, S., Gamliel, A., Li, A. C., Perissi, V., et al. (2005). A SUMOylation-dependent pathway mediates transrepression of inflammatory response genes by PPAR-gamma. Nature 437, 759–763. doi: 10.1038/nature03988

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathan, N., Faust, S. N., and Levin, M. (2003). Pathophysiology of meningococcal meningitis and septicaemia. Arch. Dis. Child. 88, 601–607. doi: 10.1136/adc.88.7.601

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathan, N., Franklin, J. L., Eleftherohorinou, H., Wright, V. J., Hemingway, C. A., Waddell, S. J., et al. (2011). Myocardial depressant effects of interleukin 6 in meningococcal sepsis are regulated by p38 mitogen-activated protein kinase. Crit. Care Med. 39, 1692–1711. doi: 10.1097/CCM.0b013e3182186d27

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathan, N., Hemingway, C. A., Alizadeh, A. A., Stephens, A. C., Boldrick, J. C., Oragui, E. E., et al. (2004). Role of interleukin 6 in myocardial dysfunction of meningococcal septic shock. Lancet 363, 203–209. doi: 10.1016/S0140-6736(03)15326-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Pino-Yanes, M., Ma, S. F., Sun, X., Tejera, P., Corrales, A., Blanco, J., et al. (2011). Interleukin-1 receptor-associated kinase 3 gene associates with susceptibility to acute lung injury. Am. J. Respir. Cell Mol. Biol. 45, 740–745. doi: 10.1165/rcmb.2010-0292OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Poelstra, K., Bakker, W. W., Klok, P. A., Hardonk, M. J., and Meijer, D. K. (1997). A physiologic function for alkaline phosphatase: endotoxin detoxification. Lab. Invest. 76, 319–327.

PubMed Abstract | Google Scholar

Pourrajab, F., Yazdi, M. B., Zarch, M. B., Zarch, M. B., and Hekmatimoghaddam, S. (2015). Cross talk of the first-line defense TLRs with PI3K/Akt pathway, in preconditioning therapeutic approach. Mol.Cell Ther. 3:4. doi: 10.1186/s40591-015-0041-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Raeburn, C. D., Calkins, C. M., Zimmerman, M. A., Song, Y., Ao, L., Banerjee, A., et al. (2002). ICAM-1 and VCAM-1 mediate endotoxemic myocardial dysfunction independent of neutrophil accumulation. Am. J. Physiol. Regul. Integr. Comp. Physiol. 283, R477–R486. doi: 10.1152/ajpregu.00034.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro-Silva, A., Zhang, H., and Jeffrey, S. S. (2007). RNA extraction from ten year old formalin-fixed paraffin-embedded breast cancer samples: a comparison of column purification and magnetic bead-based technologies. BMC Mol. Biol. 8:118. doi: 10.1186/1471-2199-8-118

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricote, M., Li, A. C., Willson, T. M., Kelly, C. J., and Glass, C. K. (1998). The peroxisome proliferator-activated receptor-gamma is a negative regulator of macrophage activation. Nature 391, 79–82. doi: 10.1038/34178

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenstein, N. E., Perkins, B. A., Stephens, D. S., Popovic, T., and Hughes, J. M. (2001). Meningococcal disease. N. Engl. J. Med. 344, 1378–1388. doi: 10.1056/NEJM200105033441807

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabatine, M. S., Morrow, D. A., Higgins, L. J., Macgillivray, C., Guo, W., Bode, C., et al. (2008). Complementary roles for biomarkers of biomechanical strain ST2 and N-terminal prohormone B-type natriuretic peptide in patients with ST-elevation myocardial infarction. Circulation 117, 1936–1944. doi: 10.1161/CIRCULATIONAHA.107.728022

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarangi, P. P., Hyun, Y. M., Lerman, Y. V., Pietropaoli, A. P., and Kim, M. (2012). Role of beta1 integrin in tissue homing of neutrophils during sepsis. Shock 38, 281–287. doi: 10.1097/SHK.0b013e31826136f8

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawant, K. V., Poluri, K. M., Dutta, A. K., Sepuru, K. M., Troshkina, A., Garofalo, R. P., et al. (2016). Chemokine CXCL1 mediated neutrophil recruitment: role of glycosaminoglycan interactions. Sci. Rep. 6:33123. doi: 10.1038/srep33123

PubMed Abstract | CrossRef Full Text | Google Scholar

Schabbauer, G., Tencati, M., Pedersen, B., Pawlinski, R., and Mackman, N. (2004). PI3K-Akt pathway suppresses coagulation and inflammation in endotoxemic mice. Arterioscler. Thromb. Vasc. Biol. 24, 1963–1969. doi: 10.1161/01.ATV.0000143096.15099.ce

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiopu, A., and Cotoi, O. S. (2013). S100A8 and S100A9: DAMPs at the crossroads between innate immunity, traditional risk factors, and cardiovascular disease. Mediators Inflamm. 2013:828354. doi: 10.1155/2013/828354

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmaier, A. H. (2002). The plasma kallikrein-kinin system counterbalances the renin-angiotensin system. J. Clin. Invest. 109, 1007–1009. doi: 10.1172/JCI0215490

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, P. J. (2015). Regulation of Iron Metabolism by Hepcidin under Conditions of Inflammation. J. Biol. Chem. 290, 18975–18983. doi: 10.1074/jbc.R115.650150

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, A., Mueller, O., Stocker, S., Salowsky, R., Leiber, M., Gassmann, M., et al. (2006). The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol. 7:3. doi: 10.1186/1471-2199-7-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Seiler, C., Sharpe, A., Barrett, J. C., Harrington, E. A., Jones, E. V., and Marshall, G. B. (2016). Nucleic acid extraction from formalin-fixed paraffin-embedded cancer cell line samples: a trade off between quantity and quality? BMC Clin. Pathol. 16:17. doi: 10.1186/s12907-016-0039-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Semeraro, N., and Colucci, M. (2000-2013). “Endothelial cell perturbation disseminated intravascular coagulation,” in: Madame Curie Bioscience Database (Austin, TX: Landes Bioscience).

Google Scholar

Shao, B., Lu, M., Katz, S. C., Varley, A. W., Hardwick, J., Rogers, T. E., et al. (2007). A host lipase detoxifies bacterial lipopolysaccharides in the liver and spleen. J. Biol. Chem. 282, 13726–13735. doi: 10.1074/jbc.M609462200

PubMed Abstract | CrossRef Full Text | Google Scholar

Shrestha, N., Bahnan, W., Wiley, D. J., Barber, G., Fields, K. A., and Schesser, K. (2012). Eukaryotic initiation factor 2 (eIF2) signaling regulates proinflammatory cytokine expression and bacterial invasion. J. Biol. Chem. 287, 28738–28744. doi: 10.1074/jbc.M112.375915

PubMed Abstract | CrossRef Full Text | Google Scholar

Silasi-Mansat, R., Zhu, H., Popescu, N. I., Peer, G., Sfyroera, G., Magotti, P., et al. (2010). Complement inhibition decreases the procoagulant response and confers organ protection in a baboon model of Escherichia coli sepsis. Blood 116, 1002–1010. doi: 10.1182/blood-2010-02-269746

PubMed Abstract | CrossRef Full Text | Google Scholar

Soares, M. P., Gozzelino, R., and Weis, S. (2014). Tissue damage control in disease tolerance. Trends Immunol. 35, 483–494. doi: 10.1016/j.it.2014.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Standage, S. W., Caldwell, C. C., Zingarelli, B., and Wong, H. R. (2012). Reduced peroxisome proliferator-activated receptor alpha expression is associated with decreased survival and increased tissue bacterial load in sepsis. Shock 37, 164–169. doi: 10.1097/SHK.0b013e31823f1a00

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanton, M. C., Taylor-Robinson, D., Harris, D., Paize, F., Makwana, N., Hackett, S. J., et al. (2011). Meningococcal disease in children in Merseyside, England: a 31 year descriptive study. PLoS ONE 6:e25957. doi: 10.1371/journal.pone.0025957

PubMed Abstract | CrossRef Full Text | Google Scholar

Stephens, D. S., Greenwood, B., and Brandtzaeg, P. (2007). Epidemic meningitis, meningococcaemia, and Neisseria meningitidis. Lancet 369, 2196–2210. doi: 10.1016/S0140-6736(07)61016-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoof, S. P., Rodenburg, G. D., Knol, M. J., Rumke, L. W., Bovenkerk, S., Berbers, G. A., et al. (2015). Disease burden of invasive meningococcal disease in the Netherlands between June 1999 and June 2011: a subjective role for serogroup and clonal complex. Clin. Infect. Dis. 61, 1281–1292. doi: 10.1093/cid/civ506

PubMed Abstract | CrossRef Full Text | Google Scholar

Szekely, Y., and Arbel, Y. (2018). A Review of Interleukin-1 in Heart Disease: Where Do We Stand Today? Cardiol. Ther. 7, 25–44. doi: 10.1007/s40119-018-0104-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Takasu, O., Gaut, J. P., Watanabe, E., To, K., Fagley, R. E., Sato, B., et al. (2013). Mechanisms of cardiac and renal dysfunction in patients dying of sepsis. Am. J. Respir. Crit. Care Med. 187, 509–517. doi: 10.1164/rccm.201211-1983OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Ivanciu, L., Popescu, N., Peer, G., Hack, E., Lupu, C., et al. (2007). Sepsis-induced coagulation in the baboon lung is associated with decreased tissue factor pathway inhibitor. Am. J. Pathol. 171, 1066–1077. doi: 10.2353/ajpath.2007.070104

PubMed Abstract | CrossRef Full Text | Google Scholar

Tessier, J. P., Thurner, B., Jungling, E., Luckhoff, A., and Fischer, Y. (2003). Impairment of glucose metabolism in hearts from rats treated with endotoxin. Cardiovasc. Res. 60, 119–130. doi: 10.1016/S0008-6363(03)00320-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiru, Y., Pathan, N., Bignall, S., Habibi, P., and Levin, M. (2000). A myocardial cytotoxic process is involved in the cardiac dysfunction of meningococcal septic shock. Crit. Care Med. 28, 2979–2983. doi: 10.1097/00003246-200008000-00049

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, M., Poignee-Heger, M., Weisser, M., Wessner, S., and Belousov, A. (2013). An optimized workflow for improved gene expression profiling for formalin-fixed, paraffin-embedded tumor samples. J. Clin. Bioinforma. 3:10. doi: 10.1186/2043-9113-3-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Tipoe, T. L., Wu, W. K. K., Chung, L., Gong, M., Dong, M., Liu, T., et al. (2018). Plasminogen activator inhibitor 1 for predicting sepsis severity and mortality outcomes: a systematic review and meta-analysis. Front. Immunol. 9:1218. doi: 10.3389/fimmu.2018.01218

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyekucheva, S., Martin, N. E., Stack, E. C., Wei, W., Vathipadiekal, V., Waldron, L., et al. (2015). Comparing platforms for messenger rna expression profiling of archival formalin-fixed, paraffin-embedded tissues. J. Mol. Diagn. 17, 374–381. doi: 10.1016/j.jmoldx.2015.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

van Deuren, M., and Brandtzaeg, P. (2000). Parents' and GPs' key role in diagnosis of meningococcal septicaemia. Lancet 356, 954–955. doi: 10.1016/S0140-6736(00)02705-7

PubMed Abstract | CrossRef Full Text | Google Scholar

van Deuren, M., Brandtzaeg, P., and Van Der Meer, J. W. (2000). Update on meningococcal disease with emphasis on pathogenesis and clinical management. Clin. Microbiol. Rev. 13, 144–166. doi: 10.1128/CMR.13.1.144

PubMed Abstract | CrossRef Full Text | Google Scholar

van Deuren, M., Van Der Ven-Jongekrijg, J., Bartelink, A. K., Van Dalen, R., Sauerwein, R. W., and Van Der Meer, J. W. (1995). Correlation between proinflammatory cytokines and antiinflammatory mediators and the severity of disease in meningococcal infections. J. Infect. Dis. 172, 433–439. doi: 10.1093/infdis/172.2.433

PubMed Abstract | CrossRef Full Text | Google Scholar

van Deuren, M., Van Der Ven-Jongekrijg, J., Vannier, E., Van Dalen, R., Pesman, G., Bartelink, A. K., et al. (1997). The pattern of interleukin-1beta (IL-1beta) and its modulating agents IL-1 receptor antagonist and IL-1 soluble receptor type II in acute meningococcal infections. Blood 90, 1101–1108. doi: 10.1182/blood.V90.3.1101

PubMed Abstract | CrossRef Full Text | Google Scholar

van Maldegem, F., De Wit, M., Morsink, F., Musler, A., Weegenaar, J., and Van Noesel, C. J. (2008). Effects of processing delay, formalin fixation, and immunohistochemistry on RNA Recovery From Formalin-fixed Paraffin-embedded Tissue Sections. Diagn. Mol. Pathol. 17, 51–58. doi: 10.1097/PDM.0b013e31814b8866

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Wyngene, L., Vandewalle, J., and Libert, C. (2018). Reprogramming of basic metabolic pathways in microbial sepsis: therapeutic targets at last? EMBO Mol. Med. 10:e8712. doi: 10.15252/emmm.201708712

PubMed Abstract | CrossRef Full Text | Google Scholar

Varfolomeev, E., and Vucic, D. (2018). Intracellular regulation of TNF activity in health and disease. Cytokine 101, 26–32. doi: 10.1016/j.cyto.2016.08.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasilescu, C., Rossi, S., Shimizu, M., Tudor, S., Veronese, A., Ferracin, M., et al. (2009). MicroRNA fingerprints identify miR-150 as a plasma prognostic marker in patients with sepsis. PLoS ONE 4:e7405. doi: 10.1371/journal.pone.0007405

PubMed Abstract | CrossRef Full Text | Google Scholar

Venet, F., Bohe, J., Debard, A. L., Bienvenu, J., Lepape, A., and Monneret, G. (2005). Both percentage of gammadelta T lymphocytes and CD3 expression are reduced during septic shock. Crit. Care Med. 33, 2836–2840. doi: 10.1097/01.CCM.0000189745.66585.AE

PubMed Abstract | CrossRef Full Text | Google Scholar

von Ahlfen, S., Missel, A., Bendrat, K., and Schlumpberger, M. (2007). Determinants of RNA quality from FFPE samples. PLoS ONE 2:e1261. doi: 10.1371/journal.pone.0001261

PubMed Abstract | CrossRef Full Text | Google Scholar

Waage, A., Brandtzaeg, P., Halstensen, A., Kierulf, P., and Espevik, T. (1989a). The complex pattern of cytokines in serum from patients with meningococcal septic shock. Association between interleukin 6, interleukin 1, and fatal outcome. J. Exp. Med. 169, 333–338. doi: 10.1084/jem.169.1.333

PubMed Abstract | CrossRef Full Text | Google Scholar

Waage, A., Halstensen, A., Shalaby, R., Brandtzaeg, P., Kierulf, P., and Espevik, T. (1989b). Local production of tumor necrosis factor alpha, interleukin 1, and interleukin 6 in meningococcal meningitis. Relation to the inflammatory response. J. Exp. Med. 170, 1859–1867. doi: 10.1084/jem.170.6.1859

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Bloom, O., Zhang, M., Vishnubhakat, J. M., Ombrellino, M., Che, J., et al. (1999). HMG-1 as a late mediator of endotoxin lethality in mice. Science 285, 248–251. doi: 10.1126/science.285.5425.248

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Song, R., Wang, Z., Jing, Z., Wang, S., and Ma, J. (2018). S100A8/A9 in Inflammation. Front. Immunol. 9:1298. doi: 10.3389/fimmu.2018.01298

PubMed Abstract | CrossRef Full Text | Google Scholar

Webster, A. F., Zumbo, P., Fostel, J., Gandara, J., Hester, S. D., Recio, L., et al. (2015). Mining the archives: a cross-platform analysis of gene expression profiles in archival formalin-fixed paraffin-embedded tissues. Toxicol. Sci. 148, 460–472. doi: 10.1093/toxsci/kfv195

PubMed Abstract | CrossRef Full Text | Google Scholar

White, L. E., Hassoun, H. T., Bihorac, A., Moore, L. J., Sailors, R. M., Mckinley, B. A., et al. (2013). Acute kidney injury is surprisingly common and a powerful predictor of mortality in surgical sepsis. J. Trauma Acute Care Surg. 75, 432–438. doi: 10.1097/TA.0b013e31829de6cd

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, C. C., Chen, M., Ma, L., Phan, Q. N., Inman, J. M., Kozhich, O. A., et al. (2003). A new strategy to amplify degraded RNA from small tissue samples for microarray studies. Nucleic Acids Res. 31:e53. doi: 10.1093/nar/gng053

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, H., Zhang, S., Weinheimer, C., Kovacs, A., and Muslin, A. J. (2000). 14-3-3 proteins block apoptosis and differentially regulate MAPK cascades. EMBO J. 19, 349–358. doi: 10.1093/emboj/19.3.349

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Guo, M., Jiang, W., Dong, H., Han, Y., An, X. F., et al. (2016). Endoplasmic reticulum stress and its effects on renal tubular cells apoptosis in ischemic acute kidney injury. Ren. Fail. 38, 831–837. doi: 10.3109/0886022X.2016.1160724

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Jagannath, C., Liu, X. D., Sharafkhaneh, A., Kolodziejska, K. E., and Eissa, N. T. (2007). Toll-like receptor 4 is a sensor for autophagy associated with innate immunity. Immunity 27, 135–144. doi: 10.1016/j.immuni.2007.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Zager, R. A., Vijayan, A., and Johnson, A. C. (2012). Proximal tubule haptoglobin gene activation is an integral component of the acute kidney injury “stress response”. Am. J. Physiol. Renal Physiol. 303, F139–148. doi: 10.1152/ajprenal.00168.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zelcer, N., and Tontonoz, P. (2006). Liver X receptors as integrators of metabolic and inflammatory signaling. J. Clin. Invest. 116, 607–614. doi: 10.1172/JCI27883

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Mo, M., Ding, W., Liu, W., Yan, D., Deng, J., et al. (2014). High-mobility group box 1 (HMGB1) impaired cardiac excitation-contraction coupling by enhancing the sarcoplasmic reticulum (SR) Ca(2+) leak through TLR4-ROS signaling in cardiomyocytes. J. Mol. Cell. Cardiol. 74, 260–273. doi: 10.1016/j.yjmcc.2014.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Li, Y., Yao, X., Wang, H., Zhao, L., Jiang, H., et al. (2016). miR-182 regulates metabolic homeostasis by modulating glucose utilization in muscle. Cell Rep. 16, 757–768. doi: 10.1016/j.celrep.2016.06.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T. N., Goodwin, J. E., Liu, B., Li, D., Wen, R., Yang, N., et al. (2019). Characterization of long noncoding RNA and mRNA profiles in sepsis-induced myocardial depression. Mol. Ther. Nucleic Acids 17, 852–866. doi: 10.1016/j.omtn.2019.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T. N., Li, D., Xia, J., Wu, Q. J., Wen, R., Yang, N., et al. (2017). Non-coding RNA: a potential biomarker and therapeutic target for sepsis. Oncotarget 8, 91765–91778. doi: 10.18632/oncotarget.21766

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Liang, Z., Zhang, Y., Zhu, M., Zhu, Y., Li, S., et al. (2019). Specific PIWI-interacting small noncoding RNA expression patterns in pulmonary tuberculosis patients. Epigenomics 11, 1779–1794. doi: 10.2217/epi-2018-0142

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Fang, J., and Ma, H. (2018). Inhibition of miR-182-5p protects cardiomyocytes from hypoxia-induced apoptosis by targeting CIAPIN1. Biochem. Cell Biol. 96, 646–654. doi: 10.1139/bcb-2017-0224

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H., Tang, Y., Ivanciu, L., Centola, M., Lupu, C., Taylor, F. B. Jr., et al. (2007). Temporal dynamics of gene expression in the lung in a baboon model of E. coli sepsis. BMC Genomics 8:58. doi: 10.1186/1471-2164-8-58

CrossRef Full Text | Google Scholar

Zhu, J., Deane, N. G., Lewis, K. B., Padmanabhan, C., Washington, M. K., Ciombor, K. K., et al. (2016). Evaluation of frozen tissue-derived prognostic gene expression signatures in FFPE colorectal cancer samples. Sci. Rep. 6:33273. doi: 10.1038/srep33273

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: FFPE, IHC, MOF, Neisseria meningitidis, septic shock, transcriptional profiles

Citation: Brusletto BS, Løberg EM, Hellerud BC, Goverud IL, Berg JP, Olstad OK, Gopinathan U, Brandtzaeg P and Øvstebø R (2020) Extensive Changes in Transcriptomic “Fingerprints” and Immunological Cells in the Large Organs of Patients Dying of Acute Septic Shock and Multiple Organ Failure Caused by Neisseria meningitidis. Front. Cell. Infect. Microbiol. 10:42. doi: 10.3389/fcimb.2020.00042

Received: 07 November 2019; Accepted: 22 January 2020;
Published: 19 February 2020.

Edited by:

Prajwal Gurung, The University of Iowa, United States

Reviewed by:

Roger Derek Pechous, University of Arkansas for Medical Sciences, United States
Carlos J. Orihuela, The University of Alabama at Birmingham, United States

Copyright © 2020 Brusletto, Løberg, Hellerud, Goverud, Berg, Olstad, Gopinathan, Brandtzaeg and Øvstebø. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Berit Sletbakk Brusletto, YmVyaXQuYnJ1c2xldHRvJiN4MDAwNDA7bWVkaXNpbi51aW8ubm8=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.