RÉSUMÉ
Altered myeloid inflammation and lymphopenia are hallmarks of severe infections, including with SARS-CoV-2. Here, we identified a gene program, defined by correlation with EN-RAGE (S100A12) gene expression, which was up-regulated in airway and blood myeloid cells from COVID-19 patients. The EN-RAGE program was expressed in 7 cohorts and observed in patients with both COVID-19 and acute respiratory distress syndrome (ARDS) from other causes. This program was associated with greater clinical severity and predicted future mechanical ventilation and death. EN-RAGE+ myeloid cells express features consistent with suppressor cell functionality, with low HLA-DR and high PD-L1 surface expression and higher expression of T cell-suppressive genes. Sustained EN-RAGE signature expression in airway and blood myeloid cells correlated with clinical severity and increasing expression of T cell exhaustion markers, such as PD-1. IL-6 treatment of monocytes in vitro upregulated many of the severity-associated genes in the EN-RAGE gene program, along with potential mediators of T cell suppression, such as IL-10. Blockade of IL-6 signaling by tocilizumab in a placebo-controlled clinical trial led to a rapid normalization of ENRAGE and T cell gene expression. This identifies IL-6 as a key driver of myeloid dysregulation associated with worse clinical outcomes in COVID-19 patients and provides insights into shared pathophysiological mechanisms in non-COVID-19 ARDS.
RÉSUMÉ
Tissue- and organism-level biological processes often involve coordinated action of multiple distinct cell types. Current computational methods for the analysis of single-cell RNA-sequencing (scRNA-seq) data, however, are not designed to capture co-variation of cell states across samples, in part due to the low number of biological samples in most scRNA-seq datasets. Recent advances in sample multiplexing have enabled population-scale scRNA-seq measurements of tens to hundreds of samples. To take advantage of such datasets, here we introduce a computational approach called single-cell Interpretable Tensor Decomposition (scITD). This method extracts "multicellular gene expression patterns" that vary across different biological samples. These patterns capture how changes in one cell type are connected to changes in other cell types. The multicellular patterns can be further associated with known covariates (e.g., disease, treatment, or technical batch effects) and used to stratify heterogeneous samples. We first validated the performance of scITD using in vitro experimental data and simulations. We then applied scITD to scRNA-seq data on peripheral blood mononuclear cells (PBMCs) from 115 patients with systemic lupus erythematosus and 56 healthy controls. We recapitulated a well-established pan-cell-type signature of interferon-signaling that was associated with the presence of anti-dsDNA autoantibodies and a disease activity index. We further identified a novel multicellular pattern that appears to potentiate renal involvement for patients with anti-dsDNA autoantibodies. This pattern was characterized by an expansion of activated memory B cells along with helper T cell activation and is predicted to be mediated by an increase in ICOSLG-ICOS interaction between monocytes and helper T cells. Finally, we applied scITD to two PBMC datasets from patients with COVID-19 and identified reproducible multicellular patterns that stratify patients by disease severity. Overall, scITD is a flexible method for exploring co-variation of cell states in multi-sample single-cell datasets, which can yield new insights into complex non-cell-autonomous dependencies that define and stratify disease.
RÉSUMÉ
Type I interferon (IFN-I) neutralizing autoantibodies have been found in some critical COVID-19 patients; however, their prevalence and longitudinal dynamics across the disease severity scale, and functional effects on circulating leukocytes remain unknown. Here, in 284 COVID-19 patients, we found IFN-I autoantibodies in 19% of critical, 6% of severe and none of the moderate cases. Longitudinal profiling of over 600,000 peripheral blood mononuclear cells using multiplexed single-cell epitope and transcriptome sequencing from 54 COVID-19 patients, 15 non-COVID-19 patients and 11 non-hospitalized healthy controls, revealed a lack of IFN-I stimulated gene (ISG-I) response in myeloid cells from critical cases, including those producing anti-IFN-I autoantibodies. Moreover, surface protein analysis showed an inverse correlation of the inhibitory receptor LAIR-1 with ISG-I expression response early in the disease course. This aberrant ISG-I response in critical patients with and without IFN-I autoantibodies, supports a unifying model for disease pathogenesis involving ISG-I suppression via convergent mechanisms.