Patient characteristics
The patient characteristics on study inclusion, hospital length of stay (LOS), and hospital mortality as well as results of statistical patient subgroup comparisons are tabulated for the total cohort and each analytical subcohort in Additional File 1: Tables S2–S8. The patient-level information on subcohort membership, SIRS etiology (type of surgery), infectious focus and agent as well as the SOFA score is available from heiDATA (https://heidata.uni-heidelberg.de/privateurl.xhtml?token=1fe2d50e-57f1-4e52-977a-93797d68b9c9).
Sixty patients in total were included: 24 SIRS (of which 21 were post-surgical), 25 with sepsis, and 11 with septic shock. A detailed description of the patient characteristics for the total cohort is given in Additional File 1: Text S1. Briefly, the average (standard deviation) age was 65.4 (12.9) years. Forty percent were female. General surgery was the most frequently admitting department (55.0%), and the mean hospital LOS was 24.6 (22.5) days. There were no statistically significant subgroup differences across the total cohort in demographics, admitting department, and hospital LOS. As expected, patients with septic shock had higher disease severity as measured by the SOFA score (median SOFA = 10) and hospital mortality (81.8%) than those with SIRS (median SOFA = 4, hospital mortality = 12.5%) and sepsis (median SOFA = 4, hospital mortality = 16.0%). Tumor disease was more frequent in SIRS (62.5%), where it was a common reason for surgery, than in sepsis (28.0%) and septic shock (18.2%). Higher levels of the inflammatory blood markers WBC (white blood cell) count and CRP in the latter two subgroups compared to SIRS agreed with the presence of infection. Noteworthily, 1.43-fold higher mean WBC counts [18.0 (9.86) versus 12.6 (3.36) × 109/L, p = 0.0154] and 2.18-fold higher mean levels of CRP [211 (99.8) versus 96.8 (62.5) mg/L, p < 0.0001] were the only statistically significant differences between sepsis and SIRS besides more frequent tumor disease in SIRS.
The clinical characteristics of our three patient subgroups in the total cohort and each of the analytical subcohorts are juxtaposed as radar charts in Additional File 1: Fig. S2. In all subcohorts, septic shock remained associated with the worst values (Additional File 1: Tables S2–S8). The following summary of clinical characteristics in the subcohorts thus focuses on the comparison of sepsis and SIRS. In the 51- and 42-patient flow cytometry subcohorts used to determine precursor blood counts and the relative subpopulation abundance, respectively, in high- and low-density granulocyte fractions, subgroup differences (Additional File 1: Tables S3 and S4) were highly similar to those in the total cohort. For both subcohorts, a higher mean value for CRP in sepsis than SIRS was the only statistically significant difference between these two subgroups. This difference was also seen in the 53-patient multiplex immunoassay subcohort used to determine granule proteins in plasma. In this subcohort, the frequency of tumor disease and the mean WBC count were higher and lower, respectively, in SIRS than sepsis, as in the total cohort, and the mean blood pH was lower in sepsis than SIRS (Additional File 1: Table S5). Despite high relative variations in the small 15-patient QuantiGene™ Plex subcohort with only five patients per subgroup, only a lower mean blood pH in sepsis reached statistical significance (Additional File 1: Table S6) as in the multiplex subcohort. In both the 30-patient qRT-PCR subcohort used to analyze the total granulocytes as well as in the 25-patient qRT-PCR subcohort used for the HD and LD granulocytes, fully contained in the first qRT-PCR subcohort, only higher mean CRP values in sepsis than SIRS were statistically significant. Notably, the p-values for the sepsis versus SIRS comparison of the SOFA scores ranged between 0.29 for the qRT-PCR (HD and LD granulocytes) subcohort and 1.00 for the QuantiGene™ Plex subcohort. This argues against important differences in disease severities between patients with sepsis and SIRS across our analytical subcohorts.
Granulocyte precursor blood counts
Across all precursor populations and the polymorphonuclear neutrophils (PMNs), blood counts were highest in the septic shock and lowest in the SIRS patients (Fig. 2). The count difference between the two patient subgroups reached statistical significance for late promyelocytes, metamyelocytes, and band cells with median differences of 6.5-fold, 37.8-fold and 11.7-fold, respectively. For the intermediate sepsis subgroup, only a 2.9-fold higher median band cell count compared to the SIRS patients reached statistical significance. These results are consistent with previous reports of higher granulocyte precursor blood counts in sepsis than SIRS on ICU admission [19, 31].
Next, we assessed whether our granulocyte precursor counts correlated with disease severity measured by the concurrent SOFA score (Additional File 1: Fig. S3). Correlations were only observed for band cells (p = 0,0024, ρ = 0,513) and metamyelocytes (p = 0,0267, ρ = 0,386) when the septic shock and sepsis subgroups were combined. They persisted when the SIRS subgroup was added. However, counts of neither the most mature (PMNs) nor least mature population (late promyelocytes myelocytes) showed any correlation with SOFA.
Relative abundance of precursor and mature granulocytes in the HD and LD fractions
The median proportions of the granulocyte precursor subpopulations were generally higher in the LD than HD fractions across all patient subgroups (Fig. 3A–C) consistent with their expected enrichment in the LD fraction. For all precursors, the difference was overall most pronounced in SIRS, reaching 111-fold for late promyelocytes (Fig. 3A), followed by sepsis (Fig. 3B) and septic shock (Fig. 3C). Due to the low patient number in the latter subgroup, statistical significance was reached only for late promyelocytes. In all subgroups, mature PMNs were still by far superior in counts in either fraction, ranging from a median of 95.4% in the LD fraction from septic shock patients to 99.8% in the HD fraction from SIRS patients. PMNs did not appear depleted in the LD fraction. Nevertheless, 2%, 3%, and 4% lower proportions in the LD than the HD fractions in SIRS, sepsis, and septic shock, respectively, all reached statistical significance.
There were no differences between sepsis and SIRS in relative abundance for any granulocyte maturation stage in the LD or HD fraction (Fig. 3D).
Neutrophil granule proteins in blood plasma
Blood plasma levels of the azurophilic granule enzyme myeloperoxidase were previously shown to be elevated in sepsis and especially septic shock compared to SIRS [32] suggesting a higher degree of granulocyte degranulation in sepsis than SIRS. Degranulation is expected to lower the buoyant density of granulocytes, potentially leading to a shift from the HD to the LD fraction [33]. Therefore, we determined the plasma levels of myeloperoxidase and four additional prototypical neutrophil granule proteins: elastase 2 (azurophilic granules), NGAL and lactotransferrin (both specific granules), and cathepsin D (ficolin-containing granules) (Fig. 4).
Throughout, the median plasma levels of all five proteins were significantly higher in septic shock than in both sepsis and SIRS. They were, however, only marginally higher in sepsis than SIRS, ranging from only 1.1-fold for cathepsin D to 1.5-fold for lactotransferrin, not reaching statistical significance in any case. This supports more blood granulocyte degranulation in septic shock than in sepsis and SIRS without an important difference, however, between the latter two.
Transcriptional profiling of granulocytic differentiation in total, LD, and HD granulocytes by QuantiGene™ Plex analysis
According to our previous results, gene signatures of early terminal granulocytic differentiation, determined in total whole blood granulocytes, distinguish sepsis from SIRS on ICU admission [19]. Because all granulocyte precursor populations were enriched in the LD over the HD fraction (Fig. 3), we predicted that the sepsis-SIRS contrast in these gene signatures was more pronounced in the LD fraction and less in the HD fraction each compared to total granulocytes. To test this, we assembled a panel of signature genes that characterize the granulocytic series from promyelocytes to PMNs. We applied the following rules to select these genes based on the results of our previous transcriptome analysis [19]. First, for each of the seven canonical promyelocyte- and myelocyte-restricted pathways showing the strongest enrichment in sepsis compared to SIRS (i.e., OxPhos, Proteasome, Ribosome, Lysosome, Cell Cycle, TCA Cycle, and Fatty Acid Metabolism), the ten genes with the highest running enrichment scores were included. We additionally included here the sepsis-associated Carbon Metabolism pathway [19]. Second, signature genes of granule biogenesis pathways, differentially expressed in sepsis and SIRS as reported [19], were included if the log2 mean difference for sepsis versus SIRS was > 0.4. The pathways comprised Azurophilic Granules (promyelocytes), Specific Granules (myelocytes), Gelatinase-Containing Granules (metamyelocytes), Ficolin-Containing Granules (band cells), Secretory Vesicles (segmented neutrophils), and Cell Membrane (PMNs). Third, eight differentially expressed lysosomal/endolysosomal genes were included, for which we had independently validated higher expression in sepsis than in SIRS CD15+ cells [19]. An overview of the resultant 135-gene QuantiGene™ Plex test panel with pathway memberships is given in Additional file 1: Table S1. Total granulocytes and the derived HD and LD fractions from five patients each with SIRS, sepsis, and septic shock were subjected to QuantiGene™ Plex analysis.
Comparisons of total, HD, and LD granulocytes
We conducted pairwise comparisons of the obtained results for the total, HD, and LD granulocyte preparations. Due to the low number of five patients per subgroup, the lowest p-value obtained in any comparison was 0.0625, which applied to altogether 112 genes (83% of the test panel) summarized in Fig. 5.
Overall, the expression levels were most similar across all cell preparations and patients within the SIRS subgroup and showed the highest inter-patient variations in septic shock. In SIRS compared to sepsis HD granulocytes, the relatively high expression levels of the genes encoding the ribosomal proteins L18, L24, L27, S15, and S15A stood out. However, none of the expected enrichments of early granule biogenesis and promyelocyte- and myelocyte-restricted canonical pathways in the LD fraction was apparent in SIRS other than partially for Fatty Acid Metabolism and OxPhos. In sepsis contrarily, the expression profile overall agreed with our previous results [19]. Namely, the levels for the selected endolysosomal, granule biogenesis except for the Cell Membrane pathway, and canonical pathway genes were highest in the LD and lowest in the HD fractions. The intermediate position of the total granulocyte preparation in sepsis appeared clearest for the Ribosome pathway. This data indicates that the higher expression of signature genes of early terminal granulocytic differentiation in sepsis than SIRS was more pronounced in the LD and less in the HD fraction compared to the parent total granulocyte preparation.
Comparisons of patient subgroups
Next, we considered the pairwise comparisons of gene expression in total, HD and LD granulocytes between the three patient subgroups (Fig. 6 and Additional File 1: Fig. S4). With two exceptions, the expression levels for altogether 42 differentially expressed genes (DEGs) (31.1% of the test panel) across all subgroup comparisons were lower in SIRS than in sepsis or septic shock, and lower in sepsis than in septic shock. For sepsis versus SIRS and septic shock versus SIRS, we counted 24 and 28 DEGs, respectively, distributed over all three granulocyte preparations in each comparison. For sepsis versus SIRS, there was no overlap in the DEGs between total and HD granulocytes, and little overlap between LD and total as well as LD and HD granulocytes (Additional File 1: Fig. S4A). Contrarily, 46.4% of the DEGs from the septic shock versus SIRS comparison were shared between total and HD granulocytes, and 32.1% between all three preparations (Additional File 1: Fig. S4B). In the comparison of septic shock versus sepsis, there were merely seven DEGs each in the total and the HD granulocytes with four shared between these two preparations (Additional File 1: Fig. S4C).
Promyelocyte- and myelocyte-restricted canonical pathways were represented among the DEGs across all comparisons and granulocyte preparations (Fig. 6). Notably, the only granule biogenesis pathway genes, known to be restricted to these two earliest precursor stages [25] and indeed represented in the LD granulocytes, were CTSA and CD63 (both Azurophilic Granules) in the sepsis versus SIRS comparison, and CTSA and LAIR1 (the latter Specific Granules) in the septic shock versus SIRS comparison. As the above-mentioned exception, the two PMN-expressed Cell Membrane pathway genes, CXCR2 and CEACAM4, showed higher expression levels in HD granulocytes from patients with SIRS than sepsis. Given the absence of any differences between sepsis and SIRS in PMN blood counts (Fig. 2) and their proportions in both HD and LD granulocytes (Fig. 3D), PMNs appear to express lower levels of CXCR2 and CEACAM4 in sepsis than SIRS.
Differential expression profiles of SeptiCyte™ LAB genes and CD63 prompt their selection for RT-PCR validation
We noted that all four genes of the SeptiCyte™ LAB test [11, 12] had been selected as bona fide granule biogenesis pathway genes [25] in our 135-gene QuantiGene™ Plex test panel (Additional File 1: Table S2). These were the promyelocyte-restricted Azurophilic Granule pathway genes LAMP1 and PLAC8 as well as the PMN-restricted Cell Membrane pathway genes CEACAM4 and PLA2G7. LAMP1 and PLAC8 both appeared more highly expressed in LD than total granulocytes of SIRS patients and in LD than HD granulocytes in sepsis patients (Fig. 5), consistent with granulocyte precursors as their source, but no patient subgroup differences were seen for either gene (Fig. 6). CEACAM4, however, appeared not only more highly expressed in total and HD granulocytes than in the LD fraction of septic shock patients (Fig. 5), consistent with PMNs as their source, but also showed higher levels in SIRS than sepsis HD granulocytes (Fig. 6). The latter led us to select CEACAM4 for validation by qRT-PCR in HD and total granulocytes. Because the SeptiCyte™ LAB score is based on the expression level ratio of the Azurophilic Granule to the Cell Membrane pathway genes [11], we further included one of the former, PLAC8. Additionally, we selected the Azurophilic Granule gene CD63 [34, 35] because its levels were exclusively higher in sepsis than SIRS LD granulocytes, which would also have been our expectation for PLAC8 and LAMP1. We hypothesized that forming the expression level ratio of each PLAC8 and CD63 from LD granulocytes to CEACAM4 from the HD fraction improves the discrimination between sepsis and SIRS compared to the LD levels of these genes considered separately.
Elevated expression of CD63 but not PLAC8 in sepsis is restricted to LD granulocytes
The distributions of the qRT-PCR data are displayed in Fig. 7A–C. Higher expression levels of CD63 and PLAC8 in sepsis than SIRS reached statistical significance in both total and LD granulocytes. This difference was most pronounced for PLAC8 and only marginal for CD63 in the total granulocytes. After referencing to CEACAM4, the difference in total granulocytes persisted for PLAC8 but not CD63 (Fig. 7B). In the LD granulocytes by contrast, the sepsis-SIRS differences in CD63 and PLAC8 expression were comparable (Fig. 7C). As hypothesized above, this difference was enhanced for both genes after referencing the measurements from the LD fractions to CEACAM4 determined in the corresponding HD fractions. Taken together, this indicates that elevated expression of CD63 but not PLAC8 in sepsis is restricted to LD granulocytes and that a concomitant reduction of CEACAM4 expression in HD granulocytes improves the ratiometric discrimination between sepsis and SIRS.
Among the clinical characteristics, a higher blood value of the infection marker CRP in patients with sepsis than SIRS was the only significant patient subgroup difference in the qRT-PCR subcohort (Additional File 1: Table S7). The median CRP concentration in sepsis was 2.94-fold higher than in SIRS (Fig. 7D). This was only exceeded by a 6.20-fold higher median level of PLAC8 expression in total granulocytes (Fig. 7A).