key: cord-0874758-vxcxdm8v
authors: Reyes, Miguel; Filbin, Michael R.; Bhattacharyya, Roby P.; Sonny, Abraham; Mehta, Arnav; Billman, Kianna; Kays, Kyle R.; Pinilla-Vera, Mayra; Benson, Maura E.; Cosimi, Lisa A.; Hung, Deborah T.; Levy, Bruce D.; Villani, Alexandra-Chloe; Sade-Feldman, Moshe; Baron, Rebecca M.; Goldberg, Marcia B.; Blainey, Paul C.; Hacohen, Nir
title: Plasma from patients with bacterial sepsis or severe COVID-19 induces suppressive myeloid cell production from hematopoietic progenitors in vitro
date: 2021-06-16
journal: Sci Transl Med
DOI: 10.1126/scitranslmed.abe9599
sha: 69f410413819a7aeba49683a6a8f106480d0b05a
doc_id: 874758
cord_uid: vxcxdm8v

Bacterial sepsis and severe COVID-19 share similar clinical manifestations and are both associated with dysregulation of the myeloid cell compartment. We previously reported an expanded CD14(+) monocyte state, MS1, in patients with bacterial sepsis and validated expansion of this cell subpopulation in publicly available transcriptomics data. Here, using published datasets, we show that the gene expression program associated with MS1 correlated with sepsis severity and was up-regulated in monocytes from patients with severe COVID-19. To examine the ontogeny and function of MS1 cells, we developed a cellular model for inducing CD14(+) MS1 monocytes from healthy bone marrow hematopoietic stem and progenitor cells (HSPCs). We found that plasma from patients with bacterial sepsis or COVID-19 induced myelopoiesis in HSPCs in vitro and expression of the MS1 gene program in monocytes and neutrophils that differentiated from these HSPCs. Furthermore, we found that plasma concentrations of IL-6, and to a lesser extent IL-10, correlated with increased myeloid cell output from HSPCs in vitro and enhanced expression of the MS1 gene program. We validated the requirement for these two cytokines to induce the MS1 gene program through CRISPR-Cas9 editing of their receptors in HSPCs. Using this cellular model system, we demonstrated that induced MS1 cells were broadly immunosuppressive and showed decreased responsiveness to stimulation with a synthetic RNA analog. Our in vitro study suggests a potential role for systemic cytokines in inducing myelopoiesis during severe bacterial or SARS-CoV-2 infection.

A recent estimate suggests that one in five deaths globally is associated with sepsis (1) . To date, no targeted treatment is available for this syndrome, likely because of substantial disease heterogeneity (2, 3) and our lack of insight into sepsis immunopathology (4) . These issues are highlighted by the current coronavirus disease 2019 (COVID- 19) pandemic, wherein many clinical manifestations of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection parallel those of bacterial sepsis (5) (6) (7) (8) . Sepsis is associated with profound alterations in the peripheral immune cell compartment, including a marked reduction in lymphocyte counts (9) (10) (11) and phenotypic alteration of myeloid cells (12) (13) (14) . Monocytes from patients with sepsis have decreased responsiveness to stimuli (14) (15) (16) and have lower expression of human leukocyte antigen-DR (HLA-DR) (17) (18) (19) (20) (21) (22) , which is characteristic of monocytic myeloid-derived suppressor cells (MDSCs). In line with these findings, we recently reported an expanded CD14 + monocyte state in patients with sepsis called MS1, which is reminiscent of MDSCs (23) . Previous studies performed in mice have reported opposing effects of MDSCs on sepsis outcomes (24) (25) (26) (27) , warranting further investigation into their ontogeny and function in humans and association with prognosis in patients with sepsis.

Recent reports have also noted the expansion of abnormal myeloid cells in the blood of patients with severe COVID-19 (7, (28) (29) (30) (31) . In particular, monocytes from patients with severe SARS-CoV-2 infection show increased expression of calprotectin (S100A8 and S100A9) and EN-RAGE (S100A12) and reduced expression of class II major histocompatibility complex (MHC-II) (7, 28, 30) . Furthermore, myeloid cells from severe cases of COVID-19 have a reduced interferon (IFN) response, which is associated with diminished capacity for viral control (7, 29) . These studies provide valuable data describing the phenotypic alteration of myeloid cells in severe COVID-19; however, they do not provide insight into the specific mechanisms through which these abnormal myeloid cells are produced. Identifying specific factors that result in the induction of these cells will facilitate the study of these cells in greater detail and further our understanding of their function during SARS-CoV-2 infection.

Here, we identified the gene expression programs associated with the MS1 cell state and correlated these programs with disease severity in patients with bacterial sepsis or COVID- 19 . We showed that the MS1 cell state was induced in CD34 + human hematopoietic stem and progenitor cells (HSPCs) in vitro after incubation with plasma from patients with bacterial sepsis or severe COVID-19. We generated MS1 monocytes and neutrophils from HSPCs in vitro and demonstrated that they were immunosuppressive; we then showed that the cytokines interleukin-6 (IL-6) and IL-10 contributed to MS1 gene program induction.

To characterize the gene expression programs associated with the MS1 cell state, we analyzed monocyte single-cell RNA sequencing (scRNA-seq) datasets from four cohorts including patients with sepsis and healthy controls from our earlier study (23) using consensus non-negative matrix factorization (cNMF; Fig. 1A and fig.  S1 , A to D) (32) . We found a gene expression program that included the MS1 marker genes RETN, ALOX5AP, and IL1R2 among the top 20 genes with highest loadings (fig. S1D ). This gene program was expressed more highly in patients with sepsis compared to healthy Reyes et al. (23) Reyes et al. (23) Pearson r

Lee et al. (29) Lee et al. (29) Arunachalam et al. (7) Arunachalam et al. (7) Lee et al. (29) Wilk et al. (31) Wilk et al. (31) Wilk et al. (31) Arunachalam et al. (7) Schulte-Schrepping et al. (28) Schulte-Schrepping et al. (28) Schulte-Schrepping et al. Forest plots indicate the effect size on patient survival (log 2 standardized mean difference) of inferred MS1 (left) or MHC-II (right) gene expression program usage in each dataset (patient cohort) from bulk gene expression deconvolution. Accession numbers for the data from each dataset are listed on the left. Blue boxes indicate the effect size in an individual study, with whiskers extending to the 95% confidence interval. Size of the box is proportional to the relative sample size of the study. Blue diamonds represent the summary effect size among the patient groups, determined by integrating the standardized mean differences across all studies. The width of the diamond corresponds to its 95% confidence interval. (D) Correlation matrix of the gene weights (z scores) for the monocyte gene expression programs across five scRNA-seq datasets for cohorts of patients with bacterial sepsis or COVID-19 is presented (table S7) (7, 23, 28, 29, 31) . Gene expression modules were derived in an unbiased manner from each dataset using consensus non-negative matrix factorization (cNMF (23) . Coexpression analysis within the MS1 program revealed that most genes were associated with S100A8 (Fig. 1B) , suggesting a role for S100A8 as a driver of the MS1 gene expression program. This gene and its partner S100A9 have been implicated in the development of MDSCs in cancer (33) and sepsis (24, 34) , indicating a similarity between monocytic MDSCs and the MS1 cell state. We hypothesized that expression of the MS1 gene program may be associated with worse outcomes in bacterial sepsis. In transcriptional data from our earlier study (23) , expression of the MS1 gene program correlated with sepsis severity in patients in the intensive care unit (ICU) but not in patients presenting to the hospital emergency department with milder disease (fig. S1H). Using gene expression deconvolution, we estimated the usage of the monocyte gene expression programs in 15 cohorts included in a recent meta-analysis examining sepsis mortality (fig. S1C and table S1) (35) . We found that expression of the MS1 gene program was negatively associated with patient survival (effect size = −0.32, FDR < 0.01), whereas expression of the MHC-II gene program had the opposite relationship (effect size = 0.53, FDR < 0.01; Fig. 1C ).

Gene expression signatures similar to that of MS1 have also been described recently in severe SARS-CoV-2 infection (28, 36) but have not been systematically analyzed across multiple cohorts. To determine whether MS1 cells are similarly expanded in severe COVID-19, we analyzed four independent COVID-19 scRNA-seq datasets (7, 28, 29, 31) independently and identified gene expression programs in CD14 + monocytes from each dataset using an unbiased cNMF method ( Fig. 1A and fig. S2 ). In each of the four COVID-19 datasets, we found gene expression programs corresponding to the MS1 and MHC-II modules, as evidenced by their strong correlation (Pearson r > 0.8) with gene expression programs from our sepsis datasets (Fig. 1D ). Sim ilar to the trends that we observed in bacterial sepsis, CD14 + monocytes from patients with severe SARS-CoV-2 infection or influenza A infection had higher and lower usage of the MS1 and MHC-II gene programs, respectively (FDR < 0.05; Fig. 1E ), and showed increased MS1 scores compared to healthy controls (P < 0.01; fig. S2 ). These findings suggest that the MS1 cell state is expanded in both bacterial sepsis and severe viral infection syndromes.

We previously demonstrated that MS1-like cells could be derived from immature progenitor cells through stimulation of total bone marrow mononuclear cells (BMMCs) with lipopolysaccharide or Pam3CSK4 (23) . Because BMMCs contain a heterogeneous mix of HSPCs and mature immune cells, potential paracrine interactions among the cell types confound identification of the precise factors that cause MS1 induction. Given this limitation, we sought a method for inducing the MS1 program directly in CD34 + HSPCs purified from healthy human bone marrow. We hypothesized that cytokines circulating in the blood of patients with sepsis might induce the differentiation of MS1 cells directly from HSPCs in vitro. Upon culturing HSPCs isolated from healthy human bone marrow with plasma from patients with urosepsis (urinary tract infection with organ dysfunction) or healthy controls, we found that sepsis plasma stimulated the production of monocytes and neutrophils more than healthy plasma did (P = 0.025 and 0.004 for CD34 − CD11b + CD14 + and CD34 − CD11b + CD15 + cells, respectively; Fig. 2A ). Single-cell analysis of the differentiated cell populations showed clear trajectories of myeloid differentiation ( Fig. 2B and fig. S3 , A to D). We observed that incubation of HSPCs with plasma from patients with urosepsis resulted in the emergence of CD14 + cells with high MS1 scores compared with control plasma (P < 0.01; Fig. 2B ). cNMF analysis of the scRNA-seq datasets generated in the plasma incubation experiments identified gene expression programs similar to the MS1 and MHC-II programs in patient peripheral blood mononuclear cells (PBMCs) (23) , as evidenced by strong correlations between their gene loadings (Pearson r = 0.73 and 0.78, respectively; Fig. 2C and fig. S3C ). The MS1 and MHC-II gene expression programs were also significantly up-regulated or down-regulated (P < 0.01), respectively, in CD14 + cells derived from HSPCs incubated with sepsis plasma in vitro (Fig. 2C ). These data support our hypothesis that cytokines circulating in the blood of patients with sepsis could induce the differentiation of MS1 cells from HSPCs in vitro.

We next hypothesized that similar effects might be observed when incubating HSPCs in vitro with plasma from patients with severe COVID-19. We performed the same experiments with heat-inactivated plasma from SARS-CoV-2-infected patients with mild to severe disease (CVD1-CVD4) and uninfected controls ( fig. S3 , E and F, and table S2). Plasma from COVID-19 patients who eventually died (CVD4) stimulated the production of monocytes more strongly than did plasma from nonhospitalized patients (CVD1) (FDR = 0.02; Fig. 2D ), although to a weaker extent than did plasma from patients with sepsis. Compared with plasma from CVD1 patients, plasma from CVD4 patients induced the production of CD14 + cells with higher MS1 scores (FDR < 0.01; Fig. 2E ) and caused increased and decreased expression of the MS1 and MHC-II gene programs, respectively (P < 0.01; fig. S3 , G and H). These findings suggested that circulating factors in the blood of patients with severe COVID-19 could stimulate the induction of MS1 cells from HSPCs in vitro.

Upon incubation of HSPCs with urosepsis plasma, we observed cells with high MS1 scores and a gene expression program similar to that of neutrophils produced from differentiating HSPCs ( Fig. 2B  and fig. S3C ). Neutrophils from the blood of critically ill patients with sepsis and patients with bacteremia expressed a gene program that included S100A8 among genes with the highest loadings ( fig.  S4 , A to C) and correlated with the MS1-like module in neutrophils generated from HSPCs treated with urosepsis plasma (Pearson r = 0.58; fig. S4D ). MS1 marker genes were also among the top differentially expressed genes in neutrophils from critically ill patients with sepsis or bacteremic patients compared to healthy individuals ( fig. S4E ). These findings are consistent with previous reports of MDSCs having both granulocytic and monocytic subtypes (37) and highlight the similarity between the MS1 cell state and MDSCs. Shown is the number of CD34 − CD11b + CD14 + (left) and C D 3 4 − C D 1 1 b + C D 1 5 + (right) myeloid cells produced after incubation of CD34 + HSPCs in vitro with control plasma or plasma from patients with urosepsis for 7 days. Six experiments were performed for each condition in (A) (three plasma donors with two technical replicates). P values were calculated using a two-tailed Wilcoxon rank sum test. (B) Shown are uniform manifold approximation and projection (UMAP) projections of scRNA-seq data from the experiment with HSPCs incubated with urosepsis plasma shown in (A). Colors indicate the plasma pool with which the CD34 + HSPCs were treated (left) or the MS1 gene expression score for each cell (right). Major immune cell types are labeled on the basis of expression of known marker genes. The experiment in (B) was performed on CD34 + HSPCs from two healthy bone marrow donors with two plasma donors for each condition; a total of 3039 and 5254 cells were profiled for the control plasma and urosepsis plasma treatment, respectively. (C) Gene weight correlations between the MS1 gene expression program (top) or MHC-II gene expression program (bottom) in experiments with HSPCs incubated with urosepsis plasma (x axis) and in PBMCs from patients with sepsis (y axis) are shown. Significance of the correlations (Pearson r) was calculated with a permutation test. Genes that were not detected (n.d.) in the HSPC-plasma incubation experiment but were among the top 30 genes for the corresponding gene program in the PBMC dataset are indicated. Insets show violin plots of gene program usage across different plasma treatment conditions. TPM, transcripts per million. (D) Shown is the number of CD34 − CD11b + CD14 + (left) and CD34 − CD11b + CD15 + (right) myeloid cells produced after 7 days of incubation of CD34 + HSPCs with plasma from patients that were not infected (Control, n = 10 patients) or infected with SARS-CoV-2 [CVD1 (n = 9), non-hospitalized; CVD2 (n = 14), hospitalized; CVD3 (n = 14), ICU; CVD4 (n = 10), deceased]. FDR values are shown when comparing plasma for each disease state to mild COVID-19 (CVD1) patient plasma (two-tailed Wilcoxon rank sum test, corrected for testing of multiple cohorts). (E) Shown are UMAP projections of scRNA-seq data from experiments incubating CD34 + HSPCs with COVID-19 plasma. Colors indicate the plasma pool with which the CD34 + HSPCs were treated (left) or the MS1 gene expression score for each cell (right). Major immune cell types are labeled on the basis of expression of known marker genes. The experiment in (E) was performed with HSPCs from two healthy bone marrow donors using pooled plasma from all donors in (D); a total of 4449, 4591, 3129, and 3711 cells were profiled after incubation of HSPCs with plasma from patients with mild to severe COVID-19, respectively. Inset shows violin plots of MS1 gene expression scores for CD14-expressing cells from each plasma treatment condition. Dashed line indicates the mean MS1 score in cells from the MS1 cluster in the PBMC dataset (23) . (F) Shown is a UMAP projection of the MS1 cell cluster differentiating from CD34 + HSPCs in vitro. Cells are colored by relative RNA velocity pseudotime (top) (70) , which is derived from the ratio of spliced and unspliced transcripts in the scRNA-seq dataset, and MS1 score (bottom). (G) A clustered heatmap of the top 30 genes in the MS1 module is presented. Columns indicate individual cells ordered by velocity pseudotime. Expression values are z score normalized for each gene. (H) Volcano plot shows differential gene expression analysis (exact test) between CD34 + HSPCs treated with control or urosepsis plasma for 24 hours. Genes with FDR of <0.1 are highlighted in red, and selected genes that are up-regulated early in the MS1 gene expression program are labeled. Four experiments were performed for each condition (two plasma donors with two technical replicates). (I) Quantitation of S100A8 intracellular staining of CD34 + HSPCs treated with control or urosepsis plasma for 24 hours is shown. Four experiments were performed with six donors for control and urosepsis plasma, with two technical replicates for each plasma sample. (J) Shown is the correlation between S100A8/9 concentrations in plasma and MS1 gene expression for plasma samples from patients with urosepsis and controls (n = 28). Line and shadow indicate linear regression fit and 95% confidence interval, respectively. Significance of the correlation (Pearson r) was calculated with a two-sided permutation test. 

LGALS1

: 

Up-regulation of S100A8 is associated with expression of the MS1 gene program Analysis of monocyte differentiation from HSPCs in vitro showed different pseudo-temporal dynamics of MS1 gene expression (Fig. 2F) . A number of genes including S100A8, VCAN, and MNDA were expressed during early differentiation and remained up-regulated (Fig. 2G ). Short-term stimulation (24 hours) of CD34 + HSPCs with sepsis plasma resulted in the up-regulation of these genes ( Fig. 2H and table S3) and an increase in the fraction of S100A8 + cells in the CD34 + population (Fig. 2I) . These genes were among the core genes expressed in the MS1 program (Fig. 1B) , and their early up-regulation was in line with our hypothesis that S100A8 was an important factor driving the induction and differentiation of MS1 cells from HSPCs. In addition, the plasma concentrations of the S100A8/9 dimer correlated significantly with the MS1 gene expression program in patients with sepsis (Pearson r = 0.51, P < 0.01; Fig. 2J ).

To determine which circulating factors in plasma induced the increased production of myeloid cells from HSPCs in vitro, we analyzed the concentrations of inflammatory cytokines implicated in cytokine storm (38) in plasma from patients with COVID-19 ( fig.  S5 ). We found that IL-6 concentrations correlated with a higher production of CD14 + cells from HSPCs in vitro (Fig. 3A) , suggesting that this cytokine may be involved in the induction of myelopoiesis.

To test this hypothesis, we added plasma from patients with urosepsis or COVID-19 or from corresponding healthy controls, to CD34 + HSPCs or CD34 + HSPCs lacking the IL-6 receptor and observed a marked reduction in CD14 + cell production in the latter group (Fig. 3B) . We also observed a weaker reduction in CD14 + cell production from HSPCs in vitro specifically for urosepsis plasma when HSPCs were differentiated in the presence of IL-6 neutralizing antibodies (Fig. 3C) . These findings suggested a role of circulating IL-6 in inducing myelopoiesis during severe infections.

To determine the specific cytokines that induced the differentiation of MS1 cells from HSPCs in vitro, we measured the concentrations of inflammatory cytokines in the plasma of patients with sepsis and controls from our previous study (23) . Several cytokines displayed increased concentrations in the plasma from patients with sepsis ( fig. S6 ), but both IL-6 and IL-10 concentrations specifically correlated with the MS1 gene expression program in monocytes (Fig. 3D ). This suggested the involvement of these two cytokines in MS1 cell induction and was consistent with their role in MDSC expansion in cancer (39, 40) . We found that short-term incubation of HSPCs with urosepsis plasma or IL-6, and variably with IL-10, but not with control plasma, resulted in phosphorylation of signal transducer and activator of transcription 3 (STAT3) (Fig. 3E) , the transcription factor downstream of both IL-6 and IL-10 (41, 42) . In addition, short-term treatment of HSPCs in vitro with recombinant IL-6 or IL-10 showed a dose-dependent up-regulation of the early MS1 genes S100A8, VCAN, and MNDA ( fig. S7A ) and an increase in S100A8 + progenitor cells ( fig. S7B ), similar to the effects observed with sepsis plasma.

To further test these dependencies, we used CRISPR-Cas9 to knock out the surface receptors for IL-6 (IL6RA and IL6ST) and IL-10 (IL10RA and IL10RB) in CD34 + HSPCs before incubation with urosepsis plasma in vitro ( fig. S7, C and D) . Knockout of the IL-6 receptor resulted in the production of CD14 + cells with lower MS1 scores from HSPCs (Fig. 3F ) and reduced expression of several, but not all, MS1 genes (Fig. 3G ). In addition, in these experiments, the CD14 + cells displayed increased surface expression of HLA-DR, consistent with the absence of the MS1 phenotype ( fig. S7E ). Similar trends were observed with knockout of the IL-10 receptor albeit with weaker effects.

To test whether these two cytokines were sufficient to induce the MS1 gene expression program, we differentiated CD34 + HSPCs in the presence of IL-6 or IL-10 and with or without granulocytemacrophage colony-stimulating factor (GM-CSF) and M-CSF, which are known growth factors that support the differentiation of HSPCs into monocytes ( fig. S8 ) (43) . We found that addition of GM-CSF and M-CSF, and to a lesser extent IL-6 or IL-10, increased the fraction of CD14 + cells produced by HSPCs, but IL-10 decreased the absolute number of cells produced. We found increased expression of several MS1 genes in CD14 + cells generated from HSPCs treated with IL-6 or IL-10 in the presence of GM-CSF and M-CSF (Fig. 3H ). cNMF analysis revealed gene expression programs resembling MS1 and MHC-II programs, both of which correlated significantly with the gene expression programs of PBMCs from patients with sepsis (Pearson r = 0.57 and 0.66, P < 0.01; Fig. 3I ). Expression of these programs showed opposite trends with IL-6 and IL-10 treatment (Fig. 3J) , consistent with the negative correlation observed in the PBMC dataset ( fig. S1G ).

Previous studies suggested that increased activation of inflammatory macrophages in infection sites is associated with disease severity in both bacterial and SARS-CoV-2 infections (44) (45) (46) . Because neither MS1 monocytes nor other peripheral immune cells express notable concentrations of IL-6 and IL-10 (23), we hypothesized that the increased concentrations of these cytokines in plasma were due to inflammatory activation of macrophages at sites of infection. We analyzed published scRNA-seq datasets from bronchoalveolar lavage fluid (BALF) samples and nasopharyngeal swabs from patients with SARS-CoV-2 infection and found significantly increased MS1 scores among monocytes and macrophages in patients with severe COVID-19 (P < 0.01; fig. S9, A and B) . We also found that subsets of these cells expressed IL6, IL10, or multiple chemokines (CCL3, CCL4, CXCL4, and CXCL8) and cytokines (TNF, IL1B, and IL1A) corresponding to inflammatory activation of macrophages. We observed that the expression of both IL6 and IL10, and the inflammatory macrophage gene program, positively correlated with expression of the MS1 program in each patient ( fig. S9, C to H) .

Given their similarity to MDSCs, we hypothesized that MS1 cells may suppress the activation of T cells (47) . We generated monocytes from CD34 + HSPCs either with GM-CSF and M-CSF alone (iMono) or with IL-6, IL-10, GM-CSF, and M-CSF (iMS1) and co-incubated the cells with allogeneic PBMCs activated with anti-CD3 and anti-CD28 antibodies. We found that coculture of T cells with iMS1 cells suppressed T cell proliferation, as evidenced by a reduction in the number of cell divisions undergone by both CD4 and CD8 T cells after 4 days of coculture (Fig. 4, A and B) . The suppression of T cell proliferation by iMS1 also depended on the ratio of monocytes added to the coculture (Fig. 4C) . Consistent with this finding, we found that expression of the MS1 gene module in monocytes correlated negatively with the fraction and absolute numbers of CD4 T cells, but not CD8 T cells, among total PBMCs from patients with sepsis compared to controls from our earlier study (Pearson r = −0.58 and −0.59, P < 0.01; Fig. 4D ) (23) . These results are consistent with previous studies showing that MDSCs in patients with sepsis have immunosuppressive effects on T cells (21, 48) . In published datasets comprising PBMC samples from patients with COVID-19 (7, 28, 29, 31) , we found that increased expression of the MS1 gene program in monocytes consistently correlated with lower CD8 T cell fractions in each of the four datasets that we analyzed (P < 0.05 in two of four datasets; Fig. 4D ) and less so with CD4 T cells (P < 0.05 in one dataset) (49) .

Sepsis is a systemic disease, so we tested the effects of co-incubating iMS1 cells with cell types from other organs that are commonly dysfunctional in sepsis ( fig. S10A) (50) . We found that co-incubation of iMS1 cells with tumor necrosis factor- (TNF)-activated human umbilical cord endothelial cells (HUVECs) resulted in lower expression of cytokines and chemokines (CXCL2, CXCL8, and IL6) and down-regulation of the TNF signaling pathway compared with co-incubation with iMono cells (fig. S10, B and C, and table S4 ). iMS1 cells exerted a similar albeit weaker effect on human renal epithelial cells, as evidenced by the down-regulation of a number of TNF signaling pathway genes (NFKBIA, CCL2, TNIP2, and RELB; fig. S10D and table S5). These effects were not observed in HUVECs or in human renal epithelial cells incubated with conditioned medium from iMS1 or iMono cells, suggesting a dependence on cell-cell contacts ( fig. S10E ). Co-incubation of human renal epithelial cells with iMS1 cells resulted in the up-regulation of MMP1, an important regulator of tissue remodeling and extracellular matrix homeostasis that has been previously shown to be up-regulated during the development of sepsis (51) . Similarly, PROS1, the gene encoding protein S, an important regulator of the clotting cascade, was also up-regulated in human renal epithelial cells co-incubated with iMS1 cells, suggesting a possible involvement of MS1 cells in sepsis-related coagulopathy (52, 53) .

We next generated CD14 + cells from HSPCs in vitro using plasma from patients with mild (CVD1) or severe (CVD4) COVID-19. As expected, the core MS1 gene S100A8 was the most significantly up-regulated gene in CD14 + cells generated from HSPCs using severe COVID-19 plasma (Fig. 4E and table S6 ). Next, we stimulated these CD14 + cell populations with high-molecular weight poly-(I:C), a synthetic RNA analog, or IFN-, an important mediator of antiviral responses. We observed a weaker induction of various cytokines and IFN-stimulated genes in CD14 + cells generated with plasma from patients with severe COVID-19 in response to poly-(I:C), whereas no such effect was observed in response to IFN- (Fig. 4, F and G) . Consistent with these observations, we found that expression of the MS1 gene program in CD14 + cells correlated negatively with expression of IFN response genes in monocytes from patients with severe COVID-19 from a published study (P < 0.01; Fig. 4H ) (28) . These results correlate with our findings in bacterial sepsis (23) , wherein MS1 cells showed a diminished response to pathogen-associated molecular patterns, and are consistent with a recent report demonstrating functional impairment of myeloid cells in patients with COVID-19 (7) .

Given the similarity between MS1 cells and monocytic MDSCs, we propose that the MS1 gene expression program could provide a precise definition for monocytic MDSCs in peripheral blood. Indeed, using MDSC data from a recent study (54) , we found that ~50% of genes expressed in monocytic MDSCs, but not in monocytes, overlapped with MS1-specific genes ( fig. S11 ). However, our single-cellbased MS1 gene signature contains many additional genes, and thus may more accurately define the MDSC subsets in sepsis. Our MS1 gene signature should be analyzed in MDSCs obtained from patients with other diseases (47, (55) (56) (57) . Given the multiple mechanisms underlying T cell suppression by MDSCs suggested in cancer studies (58) (59) (60) , a systematic analysis of the interaction between myeloid cells and T cells in severe infections is needed to determine the exact mechanism behind the lymphopenia observed in patients with bacterial sepsis and severe COVID-19.

In this study, we show that the systemic circulating cytokines IL-6 and to a lesser extent IL-10 induce differentiation of HSPCs in culture into myeloid cells that express the MS1 gene expression program. Our findings provide a potential explanation for the expansion of suppressive myeloid cells in bacterial sepsis and severe COVID-19. Our study reveals that induction of myelopoiesis and the MS1 gene expression program in HSPCs in vitro depends on systemic IL-6, a current therapeutic candidate for treating severe COVID-19 (61). Our results suggest that the role of IL-6 in infection is complex and more extensive than just the induction of an acute phase immune response (62) , consistent with the observed lack of an impact of IL-6 blockade in randomized clinical trials in patients with COVID-19 (63, 64) . Consistent with our experimental findings, a recent study has shown that treatment with tocilizumab, an anti-IL6-R antibody, resulted in decreased fractions of monocytic MDSCs in the blood of patients with severe COVID-19 (65) .

Our analysis supports a model in which cytokine production by inflammatory macrophages at sites of infection may result in increased concentrations of systemic IL-6 and IL-10, promoting the expansion of myeloid cells expressing the MS1 gene program. We propose that the expansion and recruitment of immunosuppressive MS1 cells could serve as a negative feedback response to dampen excessive inflammation and tissue damage at infection sites, which in certain situations may be clinically beneficial. This expansion of immunosuppressive MS1 cells could perhaps be sustained by a feedforward loop wherein an increase in systemic S100A8/9 protein released by MS1 cells could further stimulate Toll-like receptor 4 and lead to continued production of suppressive myeloid cells from progenitor cells (23, 66) .

We also demonstrate that MS1 cells generated from HSPCs in vitro have an anti-inflammatory effect on endothelial and epithelial cell lines, raising the possibility that MS1 cells exert their functional effects on nonimmune cell types during bacterial sepsis and severe COVID-19. This supports the hypothesis that the interaction of MS1 cells with other tissues affects pathogenesis of severe infections and suggests the need for profiling the function of myeloid cells from patients with severe infections and examining their impact on other cell types.

A key limitation of our study is the in vitro nature of the experiments. The role of cytokine production at sites of infection in stimulating HSPCs in the bone marrow must still be confirmed in animal models of sepsis and COVID-19. Furthermore, additional studies in vivo are needed to determine whether suppressive myeloid cells play a causal role in the development of or protection against organ dysfunction during sepsis and severe COVID-19.

Plasma samples from patients with bacterial sepsis or healthy controls were obtained from an existing cohort of patients as described previously (23) . COVID-19 plasma samples were obtained from patients entering the Massachusetts General Hospital (MGH) Emergency Department. This study was approved by the Partners Healthcare Institutional Review Board under protocol no. 2017P001681. Patients presenting to the MGH Emergency Department from March through May 2020 with respiratory distress suspected or known to be due to COVID-19 were enrolled. Inclusion criteria were age 18 years or older, clinical concern for COVID-19 upon Emergency Department presentation, and acute respiratory distress with at least one of the following: (i) tachypnea ≥ 22 breaths per minute, (ii) oxygen saturation ≤ 92% on room air, (iii) a requirement for supplemental oxygen, or (iv) a requirement for positive-pressure ventilation. Patients were categorized on the basis of disease outcomes as follows: CVD1, nonhospitalized (mild disease); CVD2, hospitalized without intensive care; CVD3, hospitalized and admitted to intensive care; and CVD4, hospitalized and eventually died; controls were SARS-CoV-2-negative patients. Blood samples were collected upon hospital presentation. None of the patients had a bloodstream infection, and only one patient in the CVD4 group had a possible coexisting bacterial infection (urinary tract). All of the patients with COVID-19 were adjudicated to have SARS-CoV-2 as a primary infection.

No prior sample size calculations were performed. Sample size was determined to be adequate on the basis of the degree and consistency of differences between groups. The number of cells, biological replicates, and technical replicates are indicated in each figure legend. Investigators were not blinded during experiments.

Patient cohorts, sample numbers, and sample types for the transcriptomic datasets analyzed in this study are detailed in table S7. Briefly, we analyzed seven scRNA-seq datasets, including those from PBMCs, BALF, and nasopharyngeal swab samples from patients with bacterial sepsis or COVID-19. To determine the association of monocyte gene programs with sepsis mortality, we analyzed a compilation of bulk transcriptomic datasets from a previous meta-analysis (35) .

Purified CD34 + bone marrow cells from healthy individuals were either purchased directly from STEMCELL Technologies or isolated from fresh human bone marrow. Bone marrow aspirates anticoagulated with EDTA were purchased from StemExpress and processed within 24 hours of isolation. To remove red blood cells (RBCs) from the bone marrow, 1× RBC lysis buffer (eBioscience) was added directly at a 10:1 ratio to the sample. After 5 min, the cells were centrifuged at 400g for 5 min and resuspended in 1X RBC lysis buffer to further clear the sample of RBCs. The cells were then centrifuged, resuspended in fluorescence-activated cell sorting (FACS) buffer [1X phosphate-buffered saline (PBS), 2.5% fetal bovine serum, and 2 mM EDTA; Invitrogen] and purified using human CD34 MicroBeads (Miltenyi Biotec). Isolated CD34 + cells were validated using flow cytometry (CD34-BV650, clone 561, BioLegend) to be of >90% purity. Cells were cryopreserved in CryoStor CS10 (StemCell Technologies) in aliquots of 200,000 cells each. The tubes were kept at −80°C overnight and then transferred to liquid nitrogen for longterm storage.

For each experiment, CD34 + cells were first thawed and rested for 48 hours in SFEM II (StemCell Technologies) with 75 nM StemRegenin 1 (StemCell Technologies), 3.5 nM UM171 (StemCell Technologies), stem cell factor (SCF), thrombopoietin (TPO), and FMS-related receptor tyrosine kinase 3 ligand (Flt3L) (40 ng/ml) (PeproTech), and 1X penicillin-streptomycin (Gibco). Cells were subsequently cultured in the same medium supplemented with either 20% plasma sterilized through a 0.2-m filter (Millipore) or 100 ng/ml of the cytokines IL-6, IL-10, GM-CSF, and M-CSF (PeproTech). iMS1 and iMono cells were isolated after 7 days of culture using human CD14 MicroBeads (Miltenyi Biotec). Isolated CD14 + cells were validated using flow cytometry [CD14-FITC (fluorescein isothiocyanate), clone M5E2, BioLegend] to be of >90% purity.

To assess the number of myeloid cells from plasma incubation experiments, cells were stained with the following panel: CD3-APC (allophycocyanin) (clone HIT3a), CD19-APC (clone HIB19), CD56-APC (clone 5.1H11), CD14-FITC (clone M5E2), CD15-AF700 (clone HI98), CD11b-PE-Cy7 (clone ICRF44), CD34-BV650 (clone 561), and CD38-PE (phycoerythrin)/Cy5 (clone HIT2) (BioLegend). After staining, cells were resuspended in FACS buffer with 2% CountBright beads (Invitrogen) to allow determination of absolute counts during analysis. Flow cytometry data were acquired on a CytoFLEX LX (Beckman Coulter) and analyzed using FlowJo v10.1.

scRNA-seq combined with cell hashing (67) was performed as previously described (23) . Briefly, cells from multiple culture conditions were labeled with hashtag oligo (HTO) antibodies (BioLegend) and loaded on the Chromium platform using the 3′ v3 profiling chemistry (10X Genomics). Libraries were sequenced to a depth of ~25,000 reads per cell on a NextSeq 550 (Illumina). The data were aligned to the GRCh38 reference genome using cellranger v3.1 (10X Genomics). Because of their technical incompatibility with droplet-based platforms, scRNA-seq of neutrophils was performed using plate-based Smart-Seq2 as previously described (68) .

Single-cell data analysis was performed using scanpy (69) with the same preprocessing and filtering parameters described in a prior publication (23) . To identify the major cell types in the differentiation experiments, we assessed the expression of the following marker genes: HSPCs, CD34 and CD38; monocytes, CD14 and LYZ; neutrophils, ELANE and MPO; and megakaryocytes, PPBP and PF4. MS1 scores were calculated for the top 30 genes from the MS1 module derived from the sepsis dataset using the "score_genes" in scanpy (ctrl_size = 50, n_bins = 25). RNA velocity analysis was performed using the scVelo package (70) using the default parameters.

cNMF analysis was performed as detailed in a previous publication (32) . Briefly, the top 3000 variable genes for each dataset were first selected to filter the gene expression matrix. NMF was then performed with k = 5 to 25 (10 iterations for each k). The number of modules (k) for downstream analysis was selected on the basis of biological interpretability of the modules and stability of the cNMF solution. To ensure that no modules from technical artifacts were analyzed, only gene programs with mean usage of >10 across all cells were included for further analysis.

The gene loading matrix from cNMF analysis of monocytes was used to construct a reference matrix for gene expression deconvolution. To reduce the number of genes in the reference matrix, only the top 1000 variable genes within the monocyte data were included. Sepsis datasets with survival annotation were obtained from a previously published meta-analysis (35) . Gene expression deconvolution was performed using CIBERSORT (71) with a no-sum-toone constraint and absolute scoring. The resulting score matrix was then used as an input to MetaIntegrator (72) , where the effect size of each gene module was visualized using forest plots.

Plasma samples from patients with bacterial sepsis were isolated by obtaining the top layer from Ficoll gradient separation of whole blood (diluted 1:1 with 1X PBS) and were centrifuged again at 1000g for 10 min to remove cell debris. Samples were immediately stored at −80°C. Before use in experiments, plasma samples were quickly thawed at 37°C. Plasma samples from patients with bacterial sepsis were not subjected to heat inactivation.

For patients with COVID-19, plasma samples were obtained as described above, but using whole blood diluted 1:2 in RPMI 1640. Before use in experiments, plasma samples were quickly thawed at 37°C and incubated for 1 hour at 53°C to inactivate viral particles. For the knockout and scRNA-seq experiments, plasma samples were pooled across each patient category.

Samples for neutrophil scRNA-seq were from patients enrolled in the Brigham and Women's Hospital (BWH); the criteria for patient recruitment for this cohort are described elsewhere (73, 74) . Control samples consisted of blood samples from age, gender, and ethnicitymatched healthy controls obtained from Research Blood Components (MA, USA).

To maintain viability of neutrophils for scRNA-seq, fresh blood was collected in EDTA Vacutainer tubes (BD Biosciences) and processed within 4 hours. To remove RBCs from the sample, 1X RBC lysis buffer (eBioscience) was added directly at a 10:1 ratio to the sample. After 5 min, the cells were centrifuged at 400g for 5 min and resuspended in 1X RBC lysis buffer to further clear the sample of RBCs. The remaining cells were stained with a general panel: 4′,6-diamidino-2-phenylindole, CD3-APC (clone HIT3a), CD19-APC (clone HIB19), CD56-APC (clone 5.1H11), CD14-FITC (clone M5E2), CD15-AF700 (clone HI98), and CD11b-PE-Cy7 (clone ICRF44) (BioLegend). Single neutrophils were sorted using an SH800 cell sorter (Sony) into 10 l of TCL buffer (Qiagen) with 1% -mercaptoethanol (BME; Sigma-Aldrich) in individual wells of a 96-well plate.

Bulk RNA-seq processing and data analysis Bulk RNA-seq was performed using Smart-Seq2 (75) with minor modifications, as described previously (76) , using 1000 cells as input. All RNA-seq libraries were sequenced with 38 × 38 paired-end reads using NextSeq (Illumina). RNA-seq libraries were sequenced to a depth of >2 million reads per sample. STAR was used to align sequencing reads to the UCSC hg19 transcriptome, and RSEM was used to generate an expression matrix for all samples. Both raw count and transcripts per million data were analyzed using edgeR and custom python scripts.

Intracellular staining with S100A8-PE (clone REA917, Miltenyi Biotec) was performed using the Cytofix/Cytoperm kit (BD Biosciences) following the manufacturer's protocol. For staining of phosphorylated STAT3 (Y705), cells were first fixed with 4% paraformaldehyde (Electron Microscopy Sciences) for 15 min at room temperature. The cells were then washed twice with 1X PBS and resuspended in 95% ice-cold methanol and left at −20°C overnight. The permeabilized cells were then stained with a pSTAT3-Y705 antibody (clone 13A3-1, BioLegend) for 30 min on ice. Flow cytometry data were acquired on a CytoFLEX LX (Beckman Coulter) and analyzed using FlowJo v10.1.

Cas9 protein, predesigned guide RNAs targeting S100A8, S100A9, IL6ST, IL6R, IL10RA, and IL10RB, and nontargeting guide RNAs (from GeCKO v2 library) were purchased from Integrated DNA Technologies. Ribonucleoprotein complexes were assembled by combining 2.1 l of 1X PBS, 1.2 l of 100 M guide RNA, and 1.7 l of Cas9 protein (10 g/ml) and incubating at room temperature for 15 min. The complexes were added to 50,000 to 100,000 CD34 + HSPCs resuspended in 20 l of P3 (Lonza) and electroporated (program code DZ-100) using the 4D-Nucleofector system (Lonza). After electroporation, the cells were immediately transferred to 500 l of HSPC media and rested for 48 hours. Knockout efficiency in CD34 + HSPCs was assessed after 48 hours via flow cytometry using the following panel: CD34-BV650 (clone 561), CD38-PE/Cy5 (clone HIT2), CD126-APC (clone UV4), and CD210-PECy7 (clone 3F9) (BioLegend).

Samples from patients with sepsis and controls were thawed and analyzed in parallel using the LEGENDplex Human Inflammation for IFN-2, IFN-, IL-1, IL-6, IL-10, TNF, and IL-18; the LEGENDplex Human Hematopoietic Stem Cell Panels for M-CSF and GM-CSF; and the LEGENDplex Human Vascular Inflammation Panel for S100A8/9 (BioLegend). Flow cytometry data were acquired on a CytoFLEX LX (Beckman Coulter) and analyzed using FlowJo v10.1.

Samples from patients with COVID-19 and controls were analyzed using a commercially available multiplexed proximity extension assay (Olink Proteomics) as detailed by Filbin et al. (77) . Normalized protein expression (NPX) values are calculated from the number of mapped counts on a NovaSeq run (Illumina) for each protein and represent a relative quantification scale (log 2 fold change over the mean of all proteins across all plasma samples).

T cell coculture and proliferation assay T cell coculture was performed as previously described (78) with minor modifications. Briefly, tissue culture plates were coated with purified anti-CD3 (5 g/ml) (clone HIT3a, BioLegend) at 4°C overnight and subsequently washed twice with 1X PBS. PBMCs from a healthy donor (Research Blood Components) were labeled with carboxyfluorescein diacetate succinimidyl ester (CFSE) (Invitrogen) following the manufacturer's protocol. PBMCs were resuspended in SFEM II (StemCell Technologies) with purified anti-CD28 (5 g/ ml) (clone CD28.2, BioLegend) and plated at a density of 1 million cells/ml. Isolated iMS1 or iMono cells were added at different ratios as indicated. The cells were left in culture for 3 to 4 days, with media replenished after 2 days. At the end of incubation, the cells were stained with CD3-AF700 (clone OKT3), CD4-APC (clone OKT4), and CD8a-PE (clone RPA-T8) (BioLegend) to determine the amount of CFSE dilution within the T cells.

HUVEC and human renal epithelial cell coculture and flow sorting Primary HUVECs and human renal epithelial cells were purchased from Lonza and cultured in EGM-2 and REGM, respectively. To improve cell viability, HUVECs were cultured in tissue culture plates precoated with Matrigel (Corning) diluted 1:100 in EBM-2 (Lonza). Cells were used for experiments within three to five passages. Before coculture, HUVECs and human renal epithelial cells (HREs) were treated with TNF (10 ng/ml) (PeproTech) for 24 hours. For the coculture experiments, CFSE-labeled iMono or iMS1 cells were added at a 1:1 ratio to a confluent monolayer of HUVECs or HREs and incubated for 2 hours. The cells were washed three times with 1X PBS and detached by adding 1X Accutase (Innovative Cell Technologies). The cells were transferred to FACS buffer after 15 min, and 1000 CFSE-negative cells were sorted into 10 l of TCL buffer (Qiagen) with 1% BME (Sigma) for bulk RNA-seq. Conditioned media were prepared by incubating iMS1 or iMono cells at 0.5 M cells/ml in EGM-2 or REGM for 24 hours overnight.

CVD1 and CVD4 plasma-generated CD14 + cells were isolated after 7 days of culture using human CD14 MicroBeads (Miltenyi Biotec). Isolated CD14 + cells were validated using flow cytometry (CD14-FITC, clone M5E2; BioLegend) to be of >90% purity. Cells were stimulated for 4 hours with IFN- (100 ng/ml) (PeproTech) or high-molecular weight poly-(I:C) (1 g/ml) (InvivoGen). After stimulation, cells were washed twice with FACS buffer and resuspended in 20 l of TCL buffer (Qiagen) with 1% BME (Sigma) for bulk RNA-seq.

Statistical analyses were performed with either GraphPad Prism 8 software or the SciPy statistics package. For most analyses, P values were calculated using two-tailed Wilcoxon rank sum tests and corrected for multiple testing using the Benjamini-Hochberg method. Differential gene expression analyses for bulk RNA-seq data were performed using edgeR exact tests.

stm.sciencemag.org/cgi/content/full/13/598/eabe9599/DC1 Figs. S1 to S11 Tables S1 to S7 Data files S1 to S4 View/request a protocol for this paper from Bio-protocol.

Global, regional, and national sepsis incidence and mortality

Presenting symptoms independently predict mortality in septic shock: Importance of a previously unmeasured confounder

The search for efficacious new therapies in sepsis needs to embrace heterogeneity

The immunopathology of sepsis and potential therapeutic targets

SARS-CoV-2 and viral sepsis: Observations and hypotheses

Systems biological assessment of immunity to mild versus severe COVID-19 infection in humans

Extrapulmonary manifestations of COVID-19

Persistent lymphopenia after diagnosis of sepsis predicts mortality

Sepsis-induced T cell immunoparalysis: The ins and outs of impaired T cell immunity

Impact of sepsis on CD4 T cell immunity

Sepsis triggers a late expansion of functionally impaired tissuevascular inflammatory monocytes during clinical recovery

Differential gene expression in circulating CD14 + monocytes indicates the prognosis of critically ill patients with sepsis

Human monocytes undergo functional re-programming during sepsis mediated by hypoxia-inducible factor-1

Monocyte deactivation in septic patients: Restoration by IFN-gamma treatment

Bacterial reprogramming of PBMCs impairs monocyte phagocytosis and modulates adaptive T cell responses

Human leucocyte antigen (HLA-DR) gene expression is reduced in sepsis and correlates with impaired TNF response: A diagnostic tool for immunosuppression

Proof of concept study of mass cytometry in septic shock patients reveals novel immune alterations

Downregulation of blood monocyte HLA-DR in ICU patients is also present in bone marrow cells

Immunosuppression in patients who die of sepsis and multiple organ failure

A high frequency of MDSCs in sepsis patients, with the granulocytic subtype dominating in gram-positive cases

Myeloid-derived suppressor cell function and epigenetic expression evolves over time after surgical sepsis

An immune-cell signature of bacterial sepsis

Intracellular S100A9 promotes myeloid-derived suppressor cells during late sepsis

S100-alarmin-induced innate immune programming protects newborn infants from sepsis

Myeloid-derived suppressor cells control microbial sepsis

Myeloid-derived suppressor cells evolve during sepsis and can enhance or attenuate the systemic inflammatory response

Immunophenotyping of COVID-19 and influenza highlights the role of type I interferons in development of severe COVID-19

Elevated calprotectin and abnormal myeloid cell subsets discriminate severe from mild COVID-19

A single-cell atlas of the peripheral immune response in patients with severe COVID-19

Identifying gene expression programs of cell-type identity and cellular activity with single-cell RNA-Seq

Inhibition of dendritic cell differentiation and accumulation of myeloid-derived suppressor cells in cancer is regulated by S100A9 protein

IL-10 induces an immune repressor pathway in sepsis by promoting S100A9 nuclear localization and MDSC development

A community approach to mortality prediction in sepsis via gene expression analysis

Single-cell profiling of Ebola virus infection in vivo reveals viral and host transcriptional dynamics

Myeloid-derived suppressor cells coming of age

HLH Across Speciality Collaboration, UK, COVID-19: Consider cytokine storm syndromes and immunosuppression

Ascites-derived IL-6 and IL-10 synergistically expand CD14+HLA-DR−/low myeloid-derived suppressor cells in ovarian cancer patients

Distinct myeloid suppressor cell subsets correlate with plasma IL-6 and IL-10 and reduced interferon-alpha signaling in CD4+ T cells from patients with GI malignancy

Cytokine response is determined by duration of receptor and signal transducers and activators of transcription 3 (STAT3) activation

Activation of STAT3 by IL-6 and IL-10 in primary human macrophages is differentially modulated by suppressor of cytokine signaling 3

Expansion and differentiation of CD14+CD16(−) and CD14+ +CD16+ human monocyte subsets from cord blood CD34+ hematopoietic progenitors

Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19

COVID-19 severity correlates with airway epithelium-immune cell interactions identified by single-cell analysis

Synergism of TNF- and IFN- triggers inflammatory cell death, tissue damage, and mortality in SARS-CoV-2 infection and cytokine shock syndromes

Tumor-induced tolerance and immune suppression depend on the C/EBPbeta transcription factor

Early expansion of circulating granulocytic myeloid-derived suppressor cells predicts development of nosocomial infections in patients with sepsis

UPenn COVID Processing Unit

Mechanisms and treatment of organ failure in sepsis

Role of plasma matrix-metalloproteases (MMPs) and their polymorphisms (SNPs) in sepsis development and outcome in ICU patients

Scientific and Standardization Committee on DIC, and the Scientific and Standardization Committee on Perioperative and Critical Care of the International Society on Thrombosis and Haemostasis, Diagnosis and management of sepsis-induced coagulopathy and disseminated intravascular coagulation

Dynamic evolution of coagulopathy in the first day of severe sepsis: Relationship with mortality and organ failure

Identification of monocyte-like precursors of granulocytes in cancer as a mechanism for accumulation of PMN-MDSCs

Activated human hepatic stellate cells induce myeloid derived suppressor cells from peripheral blood monocytes in a CD44-dependent fashion

A human promyelocytic-like population is responsible for the immune suppression mediated by myeloid-derived suppressor cells

Generation and functional characterization of MDSC-like cells

Arginase I in myeloid suppressor cells is induced by COX-2 in lung carcinoma

Regulatory myeloid cells paralyze T cells through cell-cell transfer of the metabolite methylglyoxal

Melanoma-educated CD14+ cells acquire a myeloid-derived suppressor cell phenotype through COX-2-dependent mechanisms

Focus shifts to antibody cocktails for COVID-19 cytokine storm

Interleukin-6: Designing specific therapeutics for a complex cytokine

RCT-TCZ-COVID-19 Study Group, Effect of tocilizumab vs standard care on clinical worsening in patients hospitalized with COVID-19 pneumonia: A randomized clincial trial

Profiling of immune dysfunction in COVID-19 patients allows early prediction of disease progression

Mrp8 and Mrp14 are endogenous activators of Toll-like receptor 4, promoting lethal, endotoxin-induced shock

Cell Hashing with barcoded antibodies enables multiplexing and doublet detection for single cell genomics

Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors

SCANPY: Large-scale single-cell gene expression data analysis

Generalizing RNA velocity to transient cell states through dynamical modeling

Robust enumeration of cell subsets from tissue expression profiles

Circulating mitochondrial DNA in patients in the ICU as a marker of mortality: Derivation and validation

Inflammasome-regulated cytokines are critical mediators of acute lung injury

Full-length RNA-seq from single cells using Smart-seq2

Multiplexed enrichment and genomic profiling of peripheral blood cells reveal subset-specific immune signatures

Plasma proteomics reveals tissue-specific cell death and mediators of cell

Methods to measure MDSC immune suppressive activity in vitro and in vivo

Count matrices are also deposited in the Broad Institute Single Cell Portal. Bulk RNA-seq data are included in this manuscript (data files S2 to S4). Raw sequencing data are available for download from the Broad Data Use Oversight System

Plasma from patients with bacterial sepsis or severe COVID-19 induces suppressive myeloid cell production from hematopoietic progenitors in vitro

is a consultant to and an equity holder in a company, 10X Genomics, whose products were