- Open Access
IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation
Genome Medicine volume 13, Article number: 64 (2021)
Immunosuppressive and anti-cytokine treatment may have a protective effect for patients with COVID-19. Understanding the immune cell states shared between COVID-19 and other inflammatory diseases with established therapies may help nominate immunomodulatory therapies.
To identify cellular phenotypes that may be shared across tissues affected by disparate inflammatory diseases, we developed a meta-analysis and integration pipeline that models and removes the effects of technology, tissue of origin, and donor that confound cell-type identification. Using this approach, we integrated > 300,000 single-cell transcriptomic profiles from COVID-19-affected lungs and tissues from healthy subjects and patients with five inflammatory diseases: rheumatoid arthritis (RA), Crohn’s disease (CD), ulcerative colitis (UC), systemic lupus erythematosus (SLE), and interstitial lung disease. We tested the association of shared immune states with severe/inflamed status compared to healthy control using mixed-effects modeling. To define environmental factors within these tissues that shape shared macrophage phenotypes, we stimulated human blood-derived macrophages with defined combinations of inflammatory factors, emphasizing in particular antiviral interferons IFN-beta (IFN-β) and IFN-gamma (IFN-γ), and pro-inflammatory cytokines such as TNF.
We built an immune cell reference consisting of > 300,000 single-cell profiles from 125 healthy or disease-affected donors from COVID-19 and five inflammatory diseases. We observed a CXCL10+ CCL2+ inflammatory macrophage state that is shared and strikingly abundant in severe COVID-19 bronchoalveolar lavage samples, inflamed RA synovium, inflamed CD ileum, and UC colon. These cells exhibited a distinct arrangement of pro-inflammatory and interferon response genes, including elevated levels of CXCL10, CXCL9, CCL2, CCL3, GBP1, STAT1, and IL1B. Further, we found this macrophage phenotype is induced upon co-stimulation by IFN-γ and TNF-α.
Our integrative analysis identified immune cell states shared across inflamed tissues affected by inflammatory diseases and COVID-19. Our study supports a key role for IFN-γ together with TNF-α in driving an abundant inflammatory macrophage phenotype in severe COVID-19-affected lungs, as well as inflamed RA synovium, CD ileum, and UC colon, which may be targeted by existing immunomodulatory therapies.
Tissue inflammation is a unifying feature across disparate diseases. While tissue- and disease-specific factors shape distinct inflammatory microenvironments, seemingly unrelated diseases can respond to the same therapy. For example, anti-tumor necrosis factor (TNF) therapies have revolutionized treatment for joint inflammation in autoimmune rheumatoid arthritis (RA) , while patients with intestinal inflammation due to Crohn’s disease (CD) and ulcerative colitis (UC), collectively known as inflammatory bowel disease (IBD), also respond to anti-TNF medications . Here, we posit that the deconstruction of tissues to the level of individually characterized cells and subsequent integration of these cells from various types of inflamed tissues could provide a platform to identify shared pathologic features across diseases and provide rationale for repurposing medications in outwardly dissimilar conditions.
Recent studies have detailed features of local tissue inflammation and immune dysfunction in COVID-19 and related diseases caused by SARS and MERS coronaviruses . Consensus is building that extensive unchecked inflammation involving so-called “cytokine storm” is a driver of severe late-stage disease. A single-cell study of bronchoalveolar lavage fluid (BALF) in intubated COVID-19 patients identified two inflammatory macrophage subsets—one characterized by CCL2, CCL3, and CXCL10 expression and a second by FCN1 and S100A8—as potential mediators of pathology in this late-stage disease . The presence of these macrophage subsets in the lung correlated with elevated circulating cytokines and extensive damage to the lung and vascular tissue. Reports looking at peripheral blood from large numbers of COVID-19 patients have consistently documented lymphopenia (reduced lymphocyte frequency) paired with increased levels of CD14+ monocytes and inflammatory cytokines, such as IL1B, TNF-α, IFN-α, and IFN-γ [5,6,7]. These factors are ineffective in lowering viral load while possibly contributing to cytokine release syndrome (CRS) . Together, these studies indicate the importance of uncovering the full extent of cell states present in COVID-19 patients including within affected tissues, and in particular among macrophages. Further, the extent to which these cell states are shared between COVID-19 and other inflammatory diseases and their disease association may further clarify disease mechanisms and precisely define therapeutic targets.
Macrophages are pervasive throughout the body and pivotal to tissue homeostasis, where they tailor their function to the parenchymal functions of each tissue type. In inflammation, tissue-resident macrophages and infiltrating monocytes are activated not only by factors from the unique tissue microenvironment, but also by disease-associating factors such as byproducts of deregulated tissue homeostasis, tissue damage, gene expression differences due to genetic variants, immune reactions, and in some cases, infecting pathogens. The unprecedented plasticity and robust reactivity of macrophages and monocytes generates a spectrum of phenotypes yet to be fully defined in human disease that mediate clearance of noxious elements, but in some cases, such as in cytokine storm, aggravate disease pathology. These phenotypes include a range of pro-inflammatory and anti-microbial states that secrete key cytokines (e.g., TNF and IL-1B) and chemokines (e.g., CXCL10 and CXCL11) and other functional states geared towards debris clearance, dampening inflammation, and tissue reconstruction, as well as a variety of intermediate states [8,9,10,11]. Meta-analysis of reactive macrophage phenotypes in inflamed tissues across diseases may further refine our understanding of the complexity of human macrophage functions, identifying subsets potentially shared across immune disorders, and thereby providing a promising route towards repurposing therapeutic strategies.
Single-cell RNA-seq (scRNA-seq) has provided an opportunity to interrogate inflamed tissues and identify expanded and potentially pathogenic immune cell types . We recently defined a distinct CD14+ IL1B+ pro-inflammatory macrophage population that is markedly expanded in RA compared to osteoarthritis (OA), a non-inflammatory disease [13, 14]. Likewise, scRNA-seq studies on inflamed colonic tissues have identified inflammatory macrophage and fibroblast phenotypes with high levels of Oncostatin M (OSM) signaling factors that are associated with resistance to anti-TNF therapies . Only very recently, developments in computational methods have made it possible to meta-analyze an expansive number of cells across various tissue states, while mitigating experimental and cohort-specific artifacts [16,17,18,19,20,21,22], therein assessing shared and distinct cell states in disparate inflamed tissues.
To define the key shared immune cell compartments between inflammatory diseases with COVID-19, we meta-analyzed and integrated tissue-level single-cell profiles from five inflammatory diseases and COVID-19. We created an immune cell reference consisting of 307,084 single-cell profiles from 125 donor samples from RA synovium, systemic lupus erythematosus (SLE) kidney, UC colon, CD ileum, interstitial lung disease, and COVID-19 BALF. This single-cell reference represents comprehensive immune cell types from different disease tissues with different inflammation levels, which can be used to investigate inflammatory diseases and their connections with COVID-19 in terms of immune cell responses. Using our meta-dataset reference, we identified major immune cell lineages including macrophages, dendritic cells, T cells, B cells, NK cells, plasma cells, mast cells, and cycling lymphocytes. Among these, we found two inflammatory CXCL10+ CCL2+ and FCN1+ macrophage states that are shared between COVID-19 and several of the inflammatory diseases we analyzed. To understand the factors driving these phenotypes, we stimulated human blood-derived macrophages with eight different combinations of inflammatory disease-associated cytokines and tissue-associating stromal cells. We demonstrated that the CXCL10+ CCL2+ macrophages from severe COVID-19 lungs share a transcriptional phenotype with macrophages stimulated by TNF-α plus IFN-γ. Notably, the other two conditions wherein these macrophages are most abundant are RA and CD. As patients with RA and CD show response to anti-TNF therapies, this finding supports the approach of identifying shared cellular states in unrelated inflamed tissues to define shared responses to medications. Furthermore, janus kinase (JAK) inhibitors have also proved effective in RA, presumably in large part through targeting IFN-γ responses [8, 23, 24]. Our data collectively support the potential efficacy of JAK inhibitors and anti-TNF therapies in inflammatory macrophage responses in COVID-19 due to cellular phenotype associations with select inflammatory tissue diseases already proven to respond to these medications.
Integration of scRNA-seq profiles from multiple datasets
scRNA-seq data collection, remapping, and aggregation
To build a multi-tissue immune cell reference, we obtained the raw FASTQ files and raw count matrices from the following publicly available scRNA-seq datasets: RA synovial cells from dbGaP (Zhang, et al, 2019; phs001457.v1.p1)  and dbGaP (Stephenson, et al, 2018; phs001529.v1.p1) , SLE kidney cells from dbGaP (Arazi, et al, 2019; phs001457.v1.p1) , UC colon cells from Single Cell Portal (Smillie, et al, 2019; SCP259) , CD ileum cells from GEO (Martin, et al, 2019; GSE134809) , interstitial and pulmonary lung disease from GEO (Reyfman, et al, 2019; GSE122960) , and COVID-19 and healthy BALF cells from GEO (Liao, et al, 2020; GSE145926) . We also use the datasets from Grant et al. (GSE155249)  and Xue et al. (GSE47189)  as additional validations.
For the FASTQs that we obtained, we used Kallisto  to map the raw reads to the same kallisto index generated from GRCh38 Ensembl v100 FASTA files. We pseudo-aligned FASTQ files to this reference, corrected barcodes, sorted BUS files, and counted unique molecular identifiers (UMIs) to generate UMI-count matrices. We aggregated all the cell barcodes from 125 donor samples into one matrix. We performed consistent QC to remove the cells that expressed fewer than 500 genes or with more than 20% of the number of UMIs mapping to the mitochondrial genes, resulting in 307,084 cells in total. The number of donor samples and cells that passed QC for each tissue source, disease status, technology, and clinical data are shown in Additional file 1: Table S1.
Normalization, scaling, and feature selection
We aggregated all samples on the overlapped 17,054 genes. We then normalized each cell to 10,000 reads and log-transformed the normalized data. We then selected the top 1,000 most highly variable genes based on dispersion within each donor sample and combined these genes to form a variable gene set. Based on the pooled highly variable genes, we then scaled the aggregated data matrix to have mean 0 and variance 1. We normalized the expression matrix using the L2 norm.
Dimensionality reduction and batch effect correction
To minimize the effect from multiple datasets with different cell numbers during an unbiased scRNA-seq data integration, we performed weighted principal component analysis (PCA) and used the first 20 weighted PCs for follow-up analysis. The summation of the weights for cells from each separate single-cell dataset is equal so that each dataset contributed equally to the analysis. For all cell-type integration, we corrected batch effects on three different levels (sequencing technology, tissue source, and donor sample) simultaneously using Harmony . We use default parameters and also specified theta = 2 for each batch variable, max.iter.cluster = 30, and max.iter.harmony = 20. For Harmony batch correction, we use the same weights from the weighted PCA. For macrophage only integration, we corrected the effect from donors for the 10X data, and dataset for the CEL-seq2 data since each donor generated from CEL-seq2 data only has less than 100 cells. As outputs, we obtained batch-corrected PC embeddings where the effects from different single-cell datasets and donors are removed in low-dimensional PC space.
Quantitative evaluation of batch correction and dataset integration
Variance explained from different sources: To quantitatively measure the mixture of batch effects after correction, we estimated the sources of variance explained from gene expression on the first ten principal component embeddings. We show the proportion of variance explained from the original pre-defined immune cell type, tissue origin, technology, and donor sample. We used the R package limma  to fit the model and ANOVA to compute the percentage of variance explained:
LISI score: Meanwhile, we used a LISI (local inverse Simpson’s Index) metric to measure the mixture levels of batch labels based on local neighbors chosen at a specific perplexity [16, 22]. Specifically, we built Gaussian distribution of neighborhoods and computed these local distributions of batch probabilities p(b) using perplexity 30 on the first 20 principal components. B is the number of batches. Then, we calculated the inverse Simpson’s index:
An iLISI (integration LISI) score ranges from 1.0, which denotes no mixing, to B (the maximum score is the total number of levels in the categorical batch variable) where higher scores indicate better mixing of batches. Here batch can be tissue source, donor sample, and sequencing technology. We also calculated the cLISI (cell-type LISI), which measures integration accuracy of pre-defined cell-type annotations instead of using the same formulation. An accurate embedding has a cLISI close to 1 for every cell neighborhood, reflecting separation of distinct cell types.
We then applied unbiased graph-based clustering (Louvain ) on the top 20 batch-corrected PCs at various resolution levels (0.2, 0.4, 0.6, 0.8, 1.0). We use 0.4 as the final resolution value to gain the biological interpretations that make most sense. Then, we furthermore performed dimensionality reduction using UMAP .
Pseudo-bulk differential expression analysis
To identify robust single-cell cluster marker genes that are shared between diseases, we performed pseudo-bulk analysis by summing the raw UMI counts for each gene across cells from the same donor sample, tissue source, and cluster assignment. We modeled raw count as a negative binomial (NB) distribution and fitted a generalized linear model (GLM) for each gene accounting for tissue, sample, and nUMI using DESeq2 . We also computed AUC and P using the Wilcoxon rank-sum test by comparing pseudo-bulk samples from one cluster to the others. We use several criteria to decide statistically significant marker genes: (1) GLM-β, (2) fold change, (3) AUC, and (4) Wilcoxon rank-sum test and Bonferroni-corrected P (threshold 10−5, 0.05/5,000 tested highly variable genes). We tested all genes that were detected in more than 100 cells with non-zero UMI counts.
Identification of major immune cell-type clusters
We carefully annotated each identified immune cell-type cluster in two ways. First, we mapped the original published annotation labels [4, 13, 15, 26, 27] to our UMAP embeddings when applicable. We are able to reproduce the original cell-type subsets in our cross-disease integrative analysis. Second, we annotated the identified clusters using cell-type lineage marker genes: T cells (CD3D), NK cells (NCAM1), B cells (MS4A1), plasma cells (MZB1), macrophages (FCGR3A/CD14), dendritic cells (DCs, CD1C), mast cells (TPSAB1), and cycling cells (MKI67).
Cell culture for human blood-derived macrophages and synovial fibroblasts
We obtained human leukocyte-enriched whole blood samples from 4 healthy blood donors from the New York Blood Center and purified peripheral blood mononuclear cells (PBMC) from each using Ficoll gradient centrifugation. We isolated CD14+ monocytes from each sample using human CD14 microbeads (Miltenyi Biotec) and differentiated these cells into blood-derived macrophages for 1 day at 37 °C in Macrophage-Colony Stimulating Factor (M-CSF); 10 ng/mL) (PeproTech) and RPMI 1640 medium (Corning) supplemented with 10% defined fetal bovine serum (FBS) (HyClone), 1% penicillin-streptomycin (Thermo Fisher Scientific), and 1% l-glutamine (Thermo Fisher Scientific) in a 6-well plate at a concentration of 1.2 million cells/mL.
In parallel, we obtained human synovial fibroblasts derived from deidentified synovial tissues from RA patients undergoing arthroplasty (HSS IRB 14033). Two unique fibroblast lines were used, each paired with two distinct blood-derived macrophage donor samples. We cultured fibroblasts in alpha minimum essential medium (aMEM) (Gibco) supplemented with 10% premium FBS (R&D Systems Inc), 1% penicillin-streptomycin (Thermo Fisher Scientific), and 1% l-glutamine (Thermo Fisher Scientific) for 4 to 6 passages. To create each transwell, we seeded the mesh of polyester chambers with 0.4μm pores (Corning) with either 200,000 synovial fibroblasts or without fibroblasts for 1 day at 37 °C.
The following day, we suspended each transwell—3 with fibroblasts and 6 without fibroblasts per donor—above one well of cultured macrophages. Those transwells with fibroblasts had a fibroblast-to-macrophage ratio of 1:15. In total, we created 9 wells per donor. Next, we added IFN-β (200 pg/mL), IL-4 (20 ng/ mL), TNF-α (20 ng/mL), and/or IFN-γ (5 ng/mL) to each transwell and underlying plate per donor. All plates were incubated at 37 °C for 19 h.
RNA library preparation and sequencing
We applied a modified version of the staining protocol from CITE-seq, using only Totalseq™-A Hashing antibodies from Biolegend . We harvested macrophages from each well and aliquoted one fifth of the cells, ~ 750,000 cells per condition, for staining in subsequent steps. We washed the cells in filtered labeling buffer (PBS with 1% BSA) and resuspended in 50 μL of labeling buffer with Human TruStain FcX™ (Biolegend Cat #422302, 5 μL per stain) for 10 min at 4 °C. Next, we added 50 μL of labeling buffer for a final concentration of 1.6 ng/μL of a total-seq hashtag (1, 2, 4–9, or 12) per condition per donor for 25 min at 4 °C. Next, we washed all samples in 2 mL, 1 mL, and 1 mL of labeling buffer, sequentially. We counted the remaining cells using a cellometer (Nexcelom Cellometer Auto 1000) and aliquoted the equivalent of 60,000 cells from each condition into one Eppendorf tube per donor. From here, we filtered through a 40-μm mesh and resuspended in PBS with 0.04% BSA to a concentration of 643.7 cells/μl. We followed the Chromium Single Cell 3′ v3 kit (10x Genomics) processing instructions and super-loaded 30,000 cells per lane. We used one lane per donor, with 9 conditions multiplexed per donor sample. After cDNA generation, samples were shipped to the Brigham and Women’s Hospital Single Cell Genomics Core for cDNA amplification and sequencing. Pairs of libraries were pooled and sequenced per lane on an Illumina NovaSeq S2 with paired-end 150 base-pair reads.
Processing FASTQ reads into gene expression matrices and cell hashing
We quantified mRNA and antibody UMI counts, respectively. Cellranger v3.1.0 was used to process the raw BCL files and produce a final gene by cell barcode UMI count matrix. First, raw BCL files were demultiplexed using cellranger mkfastq to generate FASTQ files with default parameters. Then, these FASTQ files were aligned to the GRCh38 human reference genome. Gene/antibody reads were quantified simultaneously using cellranger count. Cell barcodes and UMIs were extracted for gene/hashtag antibodies for each run.
For quality control of the cells, we first performed mRNA-level cell QC and then hashtag-level QC. For the mRNA-level QC, we removed the cells that expressed fewer than 1,000 genes or more than 10% of UMIs mapping to the mitochondrial genes. For the hashtag QC, we removed the cells whose proportion of UMIs for the most abundant hashing antibody is less than 90%, and removed the cells whose ratio of the second most-abundant and first most-abundant antibody is greater than 0.10. After filtering, each cell was assigned a hashing antibody and donor sample on the most abundant hashing antibody barcode. After QC, we obtained 9,399, 8,775, 4,622, and 3,027 cells for the 4 donor samples. We then normalized UMI counts from each cell based on the total number of UMIs and log-transformed the normalized counts.
Linear modeling for experimental stimulation-specific genes from cell culture single-cell profiles
To more accurately identify gene signatures that are specific to each of the eight stimulatory conditions, we used linear models to test each gene for differential normalized gene expression across contrasts of interest. Specifically, we fit the following models:
where stim is a categorical variable that represents eight stimuli and an untreated status, 1 ∣ sample is the random effect of the 4 replicated donor samples, and nUMI (number of unique molecular identifiers) represents the technical cell-level fixed effect. We obtained the fold change, T and P value, and Bonferroni-corrected P value as measurements for each tested gene signature for each applied condition. We then generated a list of differentially expressed genes whose fold change is greater than 2 and P is smaller than the Bonferroni correction threshold 10−7 (0.05/7,000 highly variable genes × 9 conditions) for each stimulatory condition.
Testing integrative macrophage clusters for association with severe/inflamed status
We tested the association of each macrophage cluster with severe/inflamed status compared to healthy with MASC (mixed-effects modeling of associations of single cells) . We fit a logistic regression model for each identified cluster within one tissue and set the nUMIs and percent MT (% MT) content as cell-level fixed effects, and donor sample as a random effect:
where Yi,c is the odds of cell i in cluster c, βcase is the effect log (odds ratio) for case (severe COVID-19)-control (healthy) status, βtech1 is a vector of technical cell-level (nUMIs) covariate, βtech2 is a vector of technical cell-level (% mitochondrial genes) covariate, Xi is the values for cell i in technology as appropriate, and (φd| d ) is the random effect of donor d. Thus, we used this logistic regression model to test for differentially abundant macrophage clusters associated with severe COVID-19 by correcting for the technical cell-level and donor-level covariates. Similarly, we also tested for differentially abundant macrophage clusters associated with inflamed CD compared to non-inflamed CD, RA compared to OA, and inflamed UC compared to healthy colon, accounting for technical cell-level and donor-specific covariates. We generated log likelihood-ratio test MASC P values and odds ratios for each tested cluster and used Bonferroni correction to report the macrophage clusters that are statistically significantly more abundant in severe/inflamed samples compared to healthy or non-inflamed controls.
Gene score calculation
We calculated a CXCL10+ CCL2+ gene score for each single-cell profile from an external single-cell RNA-seq dataset from severe COVID-19 BALF . The gene score was calculated as the sum of counts for CXCL10+ CCL2+ genes (n = ~ 70) as a percent of total gene counts for each cell.
Pathway enrichment analysis
For pathway gene set enrichment, we use the msigdbr R package on 4872 genesets including C5 (Gene Ontology), C7 (immunologic signature), and H (Hallmarks) from MSigDB  to calculate enriched pathways of macrophage states for each disease tissue.
For all the analysis and plots, sample sizes and measures of center and confidence intervals (mean ± SD or SEM), and statistical significance are presented in the figures, figure legends, and in the text. Results were considered statistically significant when P < 0.05 by Bonferroni correction as is indicated in figure legends and text.
A reference of > 300,000 immune single-cell profiles across inflammatory diseases and COVID-19
To compare hematopoietic cells across inflammatory diseases and COVID-19 in an unbiased fashion, we aggregated 307,084 single-cell RNA-seq profiles from 125 healthy or inflammatory disease-affected tissues spanning six disorders: (1) colon from healthy individuals and patients with inflamed or non-inflamed UC ; (2) terminal ileum from patients with inflamed or non-inflamed CD ; (3) synovium from patients with RA or OA [13, 25]; (4) kidney from patients with SLE or healthy controls , (5) lung from patients with interstitial lung disease , and (6) BALF from healthy individuals and those with mild or severe COVID-19  (Fig. 1a, b, Additional file 2: Figure S1a, Additional file 1: Table S1). We developed a pipeline for multi-tissue integration and disease association at the single-cell level (Fig. 1a, “Methods”). Where feasible, we obtained raw reads and re-mapped them to the GRCh38 genome assembly. We then aggregated raw counts for 17,054 shared genes across studies into a single matrix, performed consistent quality control (QC), library size normalization, and principal component analysis  (PCA) (“Methods”). To account for different cell numbers from different datasets, we performed weighted PCA, assigning higher weights to cells from datasets with a relatively small number of cells and vice versa. In the integrated PCA embedding, we modeled and removed the effects of technology, tissue, and donor with Harmony  to identify shared cell states across studies and diseases (“Methods”). Before Harmony, cells grouped primarily based on tissue source (Additional file 2: Figure S1b). After Harmony, < 1% of the variation explained by PC1 and PC2 was attributable to tissue source and sample, while > 60% was attributable to previously defined cell types (Fig. 1c). Importantly, rare pathogenic cell types within tissue, such as germinal center B cells in inflamed UC colon and age-associated B cells in RA synovium, were identifiable in the integrated space (Additional file 2: Figure S1c). We confirmed the degree of cross sample, tissue, technology, and cell-type mixing with an independent measure of single-cell integration: LISI [16, 22] (Local Inverse Simpson’s Index). An increased iLISI (integration LISI) score after batch correction compared to before batch correction indicates a better mixing of batches after correction (Fig. 1d and Additional file 2: Figure S2a).
In this integrated space, we performed graph-based clustering  and visualization with UMAP (Uniform Manifold Approximation and Projection) . We identified 9 major cell-type clusters (Fig. 1e) present in all six tissues (Fig. 1f) and diseases (Additional file 2: Figure S2b). We labeled the clusters with canonical markers (Fig. 1g, Additional file 3: Table S2): CD3D+ T cells, NCAM1+ NK cells, MS4A1+ B cells, MZB1+ plasma cells, FCGR3A+/CD14+ macrophages, CD1C+ dendritic cells (DCs), TPSAB1+ mast cells, and MKI67+ cycling T and B cells.
While the proportion of these immune populations differed substantially among tissues, macrophages represented a major component in each tissue (Additional file 2: Figure S2c). For example, samples obtained from lung tissues and BALF, whether from healthy controls or patients with ILD and COVID-19, contained the highest proportion of macrophages (74.8% of total hematopoietic cells) (Fig. 1f, Additional file 2: Figure S2c). In contrast, while RA synovium, SLE kidney, and CD ileum contained 9.4% macrophages, T lymphocytes comprised the majority of cells in these tissues (55.7%). The UC colon samples contained 8.3% macrophages, but had a distinctively high abundance of plasma cells (42.4%) (Additional file 2: Figure S2c).
Identification of shared inflammatory macrophage states across inflammatory disease tissues and COVID-19 lungs
To resolve the heterogeneity within the macrophage compartment, we analyzed 74,373 macrophages from 108 donors and performed weighted PCA and fine clustering analysis to define shared and distinct states across diseases (Fig. 2a, Additional file 2: Figure S3a, Additional file 4: Table S3). We identified four shared macrophage states defined by different marker sets: (1) CXCL10+ CCL2+ cells, (2) FCN1+ cells, (3) MRC1+ FABP4+ cells, and (4) C1QA+ cells (Fig. 2a, b, Additional file 2: Figure S3b). The CXCL10+ CCL2+ cells and the FCN1+ cells expressed classic inflammatory genes  including IL1B, S100A8, CCL3, CXCL11, STAT1, IFNGR1, and NFKB1 (Fig. 2b, c). A higher proportion of inflammatory macrophages in severe COVID-19 expressed these inflammation-associated genes compared to healthy BALF (Additional file 2: Figure S3c). We detected the gene signature for the CXCL10+ CCL2+ inflammatory macrophage state in a higher proportion of macrophages from severe COVID-19 BALF than from other inflamed tissues (Fig. 2c).
Liao et al.  previously identified CXCL10+ CCL2+ and FCN1+ populations as inflammatory states in the COVID-19 BALF samples used in this integrated analysis. In our multi-disease clustering, the inflammatory macrophages from inflamed RA synovium and UC and CD intestinal tissue largely mapped to the same two inflammatory macrophages seen in severe COVID-19 (Fig. 2d, Additional file 2: Figure S3d-e). In most tissue types, we found all four states represented in all six tissues, and we quantified this overlap with LISI and estimated the variance explained in the PC space (Additional file 2: Figure S3f, g). Strikingly, we observed that the FCN1+ inflammatory macrophage state dominated in SLE kidney, with few in the CXCL10+ CCL2+ macrophages (Fig. 2d), suggesting that our integrative analysis was effective in identifying both shared inflammatory states while maintaining distinct patterns in a subset of tissues.
To comprehensively define markers for the two inflammatory tissue macrophage states shared across COVID-19, RA, UC, and CD, we performed a pseudo-bulk differential expression analysis (“Methods,” Additional file 5: Table S4, fold change > 2, AUC > 0.6, Bonferroni-adjusted P < 10−5). The CXCL10+ CCL2+ inflammatory macrophages displayed significantly higher expression of CXCL10, CXCL11, CCL2, CCL3, GBP1, and IDO1 in severe COVID-19, inflamed RA, and CD compared to the FCN1+ macrophages (Fig. 2e). In contrast, the FCN1+ macrophages displayed high expression of FCN1 (Ficolin-1) and a series of alarmins such as S100A8 and S100A9 in most of the inflamed tissues (Fig. 2e). Both inflammatory macrophage states showed high expression of transcription factors that promote a pro-inflammatory macrophage phenotype, STAT1 and IRF1, in inflamed RA, UC, CD, and COVID-19 BALF relative to healthy or non-inflamed tissues (Fig. 2e). Within the CXCL10+ CCL2+ state, there was notable heterogeneity across cells in terms of IL1B expression indicating the possibility of further delineation of this macrophage state (Additional file 2: Figure S4a-b). Moreover, the effect size of all genes in CXCL10+ CCL2+ and FCN1+ subsets compared with the MRC1+ FABP4+ macrophages for each tissue further highlighted a similar set of inflammatory genes with greatest fold changes across all diseases for each subset (Additional file 2: Figure S5).
As validation, we assessed the macrophage phenotypes found in a recent analysis of single cells from severe COVID-19 BALF . Notably, we observed a significant correlation between the cross-disease shared CXCL10+ CCL2+ macrophages and two monocyte-derived alveolar macrophage (MoAM) inflammatory phenotypes from this independent severe COVID-19 cohort (wherein they were referred to as MoAM1 and MoAM2)  (Additional file 2: Figure S6a-d). We further examined CXCL10+ CCL2+ macrophage-associated genes with CD14+ cells from inflamed (leukocyte-rich) RA, non-inflamed (leukocyte-poor) RA, and OA ; we observed significant enrichment of CXCL10+ CCL2+ state-specific genes (CXCL10, CXCL9, CCL3, GBP1, and IDO1), FCN1+ state-specific genes (FCN1, S100A9, CD300E, IFITM3, and CFP), and genes (IRF1, BCL2A1, and STAT1) associated with both states in the macrophages from inflamed RA compared to non-inflamed RA and OA (Additional file 2: Figure S6e). By integrating macrophages across multiple inflamed tissues, we show that inflammatory subsets identified in COVID-19 may share common phenotypes with macrophages from other inflammatory conditions.
To elucidate cell states that were phenotypically associated, we tested the association of each state with severe COVID-19 compared to healthy BALF using a logistic regression model accounting for technical cell-level and donor-specific effects  (“Methods”). We observed the CXCL10+ CCL2+ and FCN1+ states are abundant in severe COVID-19 compared to healthy BALF (Fig. 2f). The CXCL10+ CCL2+ inflammatory state was also expanded in inflamed CD compared to non-inflamed CD, RA compared to non-inflammatory OA, and inflamed UC compared to healthy colon, respectively (Fig. 2f). We indeed observed significant enrichment of the TNF-alpha signaling via nuclear factor-κB (NF-kB) pathway and the response to interferon gamma pathway in the CXCL10+ CCL2+ cells from examined inflamed tissues (Fig. 2g). Consistent with this result, we also observed reduced frequencies of MRC1+ FABP4+ macrophages in each inflamed tissue (Fig. 2f). Taken together, these results indicate that the shared CXCL10+ CCL2+ inflammatory macrophage phenotype is expanded in inflamed tissues and severe COVID-19 BALF.
Tissue inflammatory conditions that drive distinct macrophage phenotypes
To define the factors that shape disease-associated macrophage states in affected tissues, we generated human blood-derived macrophages from four donors and activated them with eight defined mixtures of inflammatory factors, focusing particularly on the effects of antiviral interferons (IFN-β and IFN-γ) and pro-inflammatory cytokines such as TNF that mediate CRS and tissue pathology in RA and IBD  (Fig. 3a, Additional file 2: Figure S7a, “Methods”). Co-cultured fibroblasts were a component in some conditions to generate factors produced by resident stroma. To reduce confounding batch effects during scRNA-seq barcode labeling, we used a single-cell antibody-based hashing strategy  to multiplex samples from different stimulatory conditions in one sequencing run (Additional file 6: Table S5, Additional file 7: Table S6). We obtained 25,823 post-QC cells after applying 10X Genomics droplet-based single-cell assay (Additional file 2: Figure S7b-d, “Methods”). In the UMAP space, a strong response to IFN-γ drove much of the observed variation; cells treated with IFN-γ clustered well apart from all other conditions (Fig. 3b). All conditions containing IFN-γ (Type II interferon) resulted in macrophages with high expression levels of the transcription factor STAT1, interferon-stimulated genes CXCL9 and CXCL10, and inflammatory receptors such as FCGR1A  (Fig. 3c). Consistent with well-established effects, macrophages stimulated by TNF induced MMP9, IL1B, and PLAUR expression while IL-4 stimulation increased expression of CCL23, MRC1, and LIPA (Fig. 3c).
Using linear models, we identified the genes with the greatest changes in expression after each stimulation and estimated the effect sizes (“Methods”). We found that 403 genes (fold change > 2, FDR < 0.05) were significantly enriched in the TNF-α and IFN-γ stimulation compared to untreated macrophages. All conditions with IFN-γ resulted in similar effect sizes for induction of CCL2, CXCL9, CXCL10, SLAMF7, and STAT1 expression—indicating a robust IFN-γ driven macrophage signature (Fig. 3d left, Additional file 2: Figure S7e). This included robust induction by IFN-γ in macrophages co-treated with TNF (Fig. 3e left). Collectively, the TNF-driven gene expression patterns appeared more modifiable by co-stimulatory factors than IFN-γ. For example, co-cultured fibroblasts further increased TNF-induced MMP9, PLAUR, and VCAN expression, while co-stimulating with IFN-γ repressed TNF induction of these genes (Fig. 3d right). Nonetheless, a portion of the TNF effect was well preserved in TNF plus IFN-γ co-stimulated cells, including genes such as CCL2, CCL3, IL1B, and NFKBIA (Fig. 3e right). TNF-α and IFN-γ ultimately generated a macrophage phenotype with increased expression of NF-kB targets such as NFKBIA, IL1B, and HLA-DRA together with STAT1 targets such as CXCL9 and CXCL10, and GBP1 and GBP5 (Fig. 3d, e).
Identification of an IFN-γ and TNF-α synergistically driven inflammatory macrophage phenotype expanded in severe COVID-19 lungs and other inflamed disease tissues
Our cross-tissue integrative analysis revealed two shared inflammatory macrophage states (Fig. 2). To further understand these cell states and the in vivo inflammatory tissue factors driving them, we integrated the single-cell transcriptomes of both the tissue macrophages and our experimentally stimulated macrophages. After combining and correcting for tissue, technology, and donor effects, we identified 7 distinct macrophage clusters (Fig. 4a). We evaluated the robustness of the clustering and observed that our clusters were stable to the choice of the variable genes used in the analysis (Additional file 2: Figure S8a). The tissue CXCL10+ CCL2+ inflammatory macrophages from UC colon, CD ileum, RA synovium, and COVID-19 BALF were transcriptionally most similar to macrophages stimulated by the combination of TNF-α plus IFN-γ in cluster 1 (Fig. 4b, c, Additional file 2: Figure S8b-c). The blood-derived macrophages in cluster 1 included macrophages stimulated by four different conditions all including IFN-γ, of which the most abundant population (37.5%) were macrophages stimulated by TNF-α with IFN-γ (Fig. 4c, d). Comparing our results to a previously reported macrophage spectrum with 28 unique stimulatory conditions , we observed the highest expression of cluster 1-associated genes in their macrophages exposed to conditions including both TNF and IFN-γ (Additional file 2: Figure S9a).
We further identified a principal component (PC1) that captures a gradient from the FCN1+ state to the CXCL10+ CCL2+ state by applying PCA analysis to the tissue-level inflammatory macrophages (Fig. 4e), suggesting a potential continuum between the inflammatory FCN1+ and CXCL10+ CCL2+ states. Aligning cells from separate tissues along PC1, we found that the majority of inflammatory macrophages in RA, UC, and CD align more closely with the FCN1+ state (Additional file 2: Figure S9b). In severe COVID-19, we observed a shift in cell frequency between the FCN1+ and CXCL10+ CCL2+ macrophages (Wilcoxon rank-sum test P = 1.4e−07, Fig. 4f). Furthermore, we mapped the experimentally stimulated blood-derived macrophages to PC1 based on the top 50 genes with the largest and smallest PC1 gene loadings. Strikingly, the TNF-α stimulated macrophages (mean − 0.27) map to the left of the FCN1+ tissue macrophages (mean − 0.14), while the IFN-γ (mean 0.10), and TNF-α and IFN-γ (mean 0.23), stimulated macrophages map to the right of the CXCL10+ CCL2+ tissue macrophages (− 0.03) (Fig. 4f). This suggests the importance of IFN-γ stimulation in order to drive a phenotype most similar to the CXCL10+ CCL2+ state, with the addition of TNF stimulation resulting in further pushing of the macrophage phenotype along the PC1 trajectory. We observed higher expression levels of PC1-associated genes, for example CXCL10, STAT1, CCL2, CCL3, NFKBIA, and GBP1, in CXCL10+ CCL2+ severe COVID-19 compared to FCN1+ cells, and higher induced expression levels of these same genes in TNF-α and IFN-γ stimulation compared to TNF-α stimulation alone (Fig. 4g). Taken together, these results suggest we are able to recapitulate the gradient observed in vivo across multiple diseases by stimulating macrophages ex vivo with synergistic combinations of IFN-γ and TNF-α.
Our study demonstrates the power of a multi-disease reference dataset to interpret cellular phenotypes and tissue states, while placing them into a broader context that may provide insights into disease etiology and rationale for repurposing medications. Such meta-datasets can increase the resolution of cell states and aid understanding of shared cellular states found in less well-understood diseases such as COVID-19. Amassing diverse tissues from > 120 donors with a wide range of diseases, we built a human tissue inflammation single-cell reference. Applying powerful computational strategies, we integrated > 300,000 single-cell transcriptomes and corrected for factors that interfere with resolving cell-intrinsic expression patterns. In particular, we have identified a CXCL10+ CCL2+ inflammatory macrophage phenotype shared between tissues affected in autoimmune disease (RA), inflammatory diseases (CD and UC), and infectious disease (COVID-19). We observed that the abundance of this population is associated with inflammation and disease severity. With integrated analysis of an ex vivo dataset, we elucidated its potential cytokine drivers: IFN-γ together with TNF-α.
Macrophages are ideal biologic indicators for the in vivo state of a tissue due to their dynamic nature, robust responses to local factors, and widespread presence in most tissues. Through our cross-disease analysis, we defined two inflammatory macrophage states that can be found in selected groups of seemingly unrelated tissues and diseases. Most notably, the CXCL10+ CCL2+ inflammatory macrophages predominate in the bronchoalveolar lavage of patients with severe COVID-19, and are also detected in synovial tissue affected by RA and inflamed intestine from patients with IBD. These cells are distinguished by high levels of CXCL10 and CXCL11, STAT1, IFNGR1, and IFNGR2, as well as CCL2 and CCL3, NFKB1, TGFB1, and IL1B. This gene expression pattern of the JAK/STAT and NF-kB-dependent cytokines implicates induction by an intriguing combination of both the IFN-induced JAK/STAT and TNF-induced NF-kB pathways and, in conjunction, the overall transcriptome program most closely aligns with macrophages stimulated by IFN-γ plus TNF-α. As both JAK inhibitors and anti-TNF medications have outstanding efficacy in treating RA and anti-TNFs are the most common medications treating inflammatory bowel disease, including Crohn’s Disease , these therapies may target the inflammatory macrophages in severe COVID-19 lung during the phase involving cytokine release syndrome .
Infection with SARS-CoV2 triggers local immune response and inflammation in the lung compartment, recruiting macrophages that release and respond to inflammatory cytokines and chemokines . This response may change with disease progression, in particular during the transition towards the cytokine storm associated with severe disease. Intriguingly, our cross-disease tissue study strongly suggests that IFN-γ is an essential component in the inflammatory macrophage phenotype in severe COVID-19. Most studies on interferons and coronaviruses have focused on Type I interferons, such as IFN-β, due to their robust capacity to interfere with viral replication . Indeed, ongoing research into the administration of recombinant IFN-β has shown promise in reducing the risk of severe COVID-19 disease . However, other studies have indicated that targeting IFN-γ may be an effective treatment for cytokine storm, a driver of severe disease in COVID-19 patients [46, 47]. Additionally, several studies have indicated that targeting IFN-γ using JAK inhibitors such as ruxolitinib, baricitinib, and tofacitinib offers effective therapeutic effects in treating severe COVID-19 patients [43, 48,49,50,51]. Clinical trials of Type II interferon inhibitors in COVID-19 are under way (NCT04337359, NCT04359290, and NCT04348695) . Recent research has also identified that the synergism of TNF-α and IFN-γ can trigger inflammatory cell death, tissue damage, and mortality in SARS-CoV-2 infection , and shown increased levels of IFN-γ, TNF-α, CXCL10, and CCL2 in the serum of severe COVID-19 patients . In agreement with these studies, our findings indicate that IFN-γ is an important mediator together with TNF-α of severe disease, in part through activating the inflammatory CXCL10+ CCL2+ macrophage subset. We hypothesize that anti-Type II interferon (like JAK inhibitors) and anti-TNF combinatorial treatment might prove effective at inhibiting the cytokine storm driving acute respiratory distress syndrome in patients with severe COVID-19. We are aware of the limited number of longitudinal BALFs from COVID-19 patients involved in our across-tissue study due to the current crisis situation, so we expect to replicate our findings in a broader generalization of COVID-19 patients in the future. Of course, the presence of an IFN-γ and TNF phenotype is an association that may not be causal. Whether targeting these cytokines is reasonable or not will depend on additional clinical investigation.
In this study, we built a single-cell immune reference from multiple inflamed disease tissues and identified two inflammatory macrophage states, CXCL10+ CCL2+ and FCN1+ inflammatory macrophages, that were shared between COVID-19 and inflammatory diseases such as RA, CD, and UC. We demonstrated that the CXCL10+ CCL2+ macrophages are transcriptionally similar to human blood-derived macrophages stimulated by IFN-γ and TNF-α and were expanded in severe COVID-19 lungs and inflamed RA, CD, and UC tissues. This finding indicates that Type II interferon and TNF responses may be involved in late-stage cytokine storm-driven severe COVID-19 and inhibiting these responses in the inflammatory macrophages may be a promising treatment. Our cross-tissue single-cell integrative strategy along with our disease association analysis provides a proof-of-principle that identifying shared pathogenic features across human inflamed tissues and COVID-19 lungs has the potential to guide drug repurposing.
Availability of data and materials
The single-cell RNA-seq data for blood-derived macrophages are available in the Gene Expression Omnibus database with accession number GSE168710, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168710 . Source code repository to reproduce analyses is located at https://github.com/immunogenomics/inflamedtissue_covid19_reference .
The publicly available datasets analyzed during the study are available from the GEO repository:
Single Cell Portal:
McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. N Engl J Med. 2011;365(23):2205–19. https://doi.org/10.1056/NEJMra1004965.
Neurath MF. Cytokines in inflammatory bowel disease. Nat Rev Immunol. 2014;14(5):329–42. https://doi.org/10.1038/nri3661.
Liu J, Zheng X, Tong Q, Li W, Wang B, Sutter K, Trilling M, Lu M, Dittmer U, Yang D. Overlapping and discrete aspects of the pathology and pathogenesis of the emerging human pathogenic coronaviruses SARS-CoV, MERS-CoV, and 2019-nCoV. J Med Virol. 2020;92(5):491–4. https://doi.org/10.1002/jmv.25709.
Liao M, Liu Y, Yuan J, Wen Y, Xu G, Zhao J, Cheng L, Li J, Wang X, Wang F, Liu L, Amit I, Zhang S, Zhang Z. Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat Med. 2020;26(6):842–4. https://doi.org/10.1038/s41591-020-0901-9.
Wen W, Su W, Tang H, le W, Zhang X, Zheng Y, Liu X, Xie L, Li J, Ye J, Dong L, Cui X, Miao Y, Wang D, Dong J, Xiao C, Chen W, Wang H. Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing. Cell Discov. 2020;6(1):31. https://doi.org/10.1038/s41421-020-0168-9.
Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, Zhang L, Fan G, Xu J, Gu X, Cheng Z, Yu T, Xia J, Wei Y, Wu W, Xie X, Yin W, Li H, Liu M, Xiao Y, Gao H, Guo L, Xie J, Wang G, Jiang R, Gao Z, Jin Q, Wang J, Cao B. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497–506. https://doi.org/10.1016/S0140-6736(20)30183-5.
Lucas C, et al. Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature. 2020;584(7821):463–9. https://doi.org/10.1038/s41586-020-2588-y.
He W, Kapate N, Shields CW 4th, Mitragotri S. Drug delivery to macrophages: a review of targeting drugs and drug carriers to macrophages for inflammatory diseases. Adv Drug Deliv Rev. 2019;165-166:15–40. https://doi.org/10.1016/j.addr.2019.12.001.
Kinne RW, Bräuer R, Stuhlmüller B, Palombo-Kinne E, Burmester GR. Macrophages in rheumatoid arthritis. Arthritis Res. 2000;2(3):189–202. https://doi.org/10.1186/ar86.
Ma W-T, Gao F, Gu K, Chen D-K. The role of monocytes and macrophages in autoimmune diseases: a comprehensive review. Front Immunol. 2019;10:1140. https://doi.org/10.3389/fimmu.2019.01140.
Xue J, Schmidt SV, Sander J, Draffehn A, Krebs W, Quester I, de Nardo D, Gohel TD, Emde M, Schmidleithner L, Ganesan H, Nino-Castro A, Mallmann MR, Labzin L, Theis H, Kraut M, Beyer M, Latz E, Freeman TC, Ulas T, Schultze JL. Transcriptome-based network analysis reveals a spectrum model of human macrophage activation. Immunity. 2014;40(2):274–88. https://doi.org/10.1016/j.immuni.2014.01.006.
Papalexi E, Satija R. Single-cell RNA sequencing to explore immune cell heterogeneity. Nat Rev Immunol. 2018;18(1):35–45. https://doi.org/10.1038/nri.2017.76.
Zhang F, et al. Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry. Nat Immunol. 2019;20(7):928–42. https://doi.org/10.1038/s41590-019-0378-1.
Kuo D, Ding J, Cohn IS, Zhang F, Wei K, Rao DA, Rozo C, Sokhi UK, Shanaj S, Oliver DJ, Echeverria AP, DiCarlo EF, Brenner MB, Bykerk VP, Goodman SM, Raychaudhuri S, Rätsch G, Ivashkiv LB, Donlin LT. HBEGF+ macrophages in rheumatoid arthritis induce fibroblast invasiveness. Sci Transl Med. 2019;11(491):eaau8587. https://doi.org/10.1126/scitranslmed.aau8587.
Smillie CS, et al. Intra- and inter-cellular rewiring of the human colon during ulcerative colitis. Cell. 2019;178:714–730.e22.
Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, Baglaenko Y, Brenner M, Loh PR, Raychaudhuri S. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods. 2019;16(12):1289–96. https://doi.org/10.1038/s41592-019-0619-0.
Stuart T, Satija R. Integrative single-cell analysis. Nat Rev Genet. 2019;20(5):257–72. https://doi.org/10.1038/s41576-019-0093-7.
Hie B, Bryson B, Berger B. Efficient integration of heterogeneous single-cell transcriptomes using Scanorama. Nat Biotechnol. 2019;37(6):685–91. https://doi.org/10.1038/s41587-019-0113-3.
Haghverdi L, Lun ATL, Morgan MD, Marioni JC. Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors. Nat Biotechnol. 2018;36(5):421–7. https://doi.org/10.1038/nbt.4091.
Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36(5):411–20. https://doi.org/10.1038/nbt.4096.
Polański K, Young MD, Miao Z, Meyer KB, Teichmann SA, Park JE. BBKNN: fast batch alignment of single cell transcriptomes. Bioinformatics. 2020;36(3):964–5. https://doi.org/10.1093/bioinformatics/btz625.
Tran HTN, Ang KS, Chevrier M, Zhang X, Lee NYS, Goh M, Chen J. A benchmark of batch-effect correction methods for single-cell RNA sequencing data. Genome Biol. 2020;21(1):12. https://doi.org/10.1186/s13059-019-1850-9.
Ivashkiv LB. IFNγ: signalling, epigenetics and roles in immunity, metabolism, disease and cancer immunotherapy. Nat Rev Immunol. 2018;18(9):545–58. https://doi.org/10.1038/s41577-018-0029-z.
Barrat FJ, Crow MK, Ivashkiv LB. Interferon target-gene expression and epigenomic signatures in health and disease. Nat Immunol. 2019;20(12):1574–83. https://doi.org/10.1038/s41590-019-0466-2.
Stephenson W, Donlin LT, Butler A, Rozo C, Bracken B, Rashidfarrokhi A, Goodman SM, Ivashkiv LB, Bykerk VP, Orange DE, Darnell RB, Swerdlow HP, Satija R. Single-cell RNA-seq of rheumatoid arthritis synovial tissue using low-cost microfluidic instrumentation. Nat Commun. 2018;9(1):791. https://doi.org/10.1038/s41467-017-02659-x.
Arazi A, et al. The immune cell landscape in kidneys of patients with lupus nephritis. Nat Immunol. 2019;20(7):902–14. https://doi.org/10.1038/s41590-019-0398-x.
Martin JC, et al. Single-cell analysis of Crohn’s disease lesions identifies a pathogenic cellular module associated with resistance to anti-TNF therapy. Cell. 2019;178:1493–1508.e20.
Reyfman PA, Walter JM, Joshi N, Anekalla KR, McQuattie-Pimentel AC, Chiu S, Fernandez R, Akbarpour M, Chen CI, Ren Z, Verma R, Abdala-Valencia H, Nam K, Chi M, Han SH, Gonzalez-Gonzalez FJ, Soberanes S, Watanabe S, Williams KJN, Flozak AS, Nicholson TT, Morgan VK, Winter DR, Hinchcliff M, Hrusch CL, Guzy RD, Bonham CA, Sperling AI, Bag R, Hamanaka RB, Mutlu GM, Yeldandi AV, Marshall SA, Shilatifard A, Amaral LAN, Perlman H, Sznajder JI, Argento AC, Gillespie CT, Dematte J, Jain M, Singer BD, Ridge KM, Lam AP, Bharat A, Bhorade SM, Gottardi CJ, Budinger GRS, Misharin AV. Single-cell transcriptomic analysis of human lung provides insights into the pathobiology of pulmonary fibrosis. Am J Respir Crit Care Med. 2018;199(12):1517–36. https://doi.org/10.1164/rccm.201712-2410OC.
Grant RA, et al. Circuits between infected macrophages and T cells in SARS-CoV-2 pneumonia. Nature. 2021;590(7847):635–41. https://doi.org/10.1038/s41586-020-03148-w.
Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34(5):525–7. https://doi.org/10.1038/nbt.3519.
Ritchie ME, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Blondel VD, Guillaume J-L, Lambiotte R, Lefebvre E. Fast unfolding of communities in large networks. arXiv [physics.soc-ph]. 2008;(10):P10008. https://arxiv.org/abs/0803.0476.
McInnes L, Healy J, Melville J. Umap: Uniform manifold approximation and projection for dimension reduction. arXiv preprint arXiv:1802.03426. 2018. https://arxiv.org/abs/1802.03426.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.
Stoeckius M, Hafemeister C, Stephenson W, Houck-Loomis B, Chattopadhyay PK, Swerdlow H, Satija R, Smibert P. Simultaneous epitope and transcriptome measurement in single cells. Nat Methods. 2017;14(9):865–8. https://doi.org/10.1038/nmeth.4380.
Fonseka CY, Rao DA, Teslovich NC, Korsunsky I, Hannes SK, Slowikowski K, Gurish MF, Donlin LT, Lederer JA, Weinblatt ME, Massarotti EM, Coblyn JS, Helfgott SM, Todd DJ, Bykerk VP, Karlson EW, Ermann J, Lee YC, Brenner MB, Raychaudhuri S. Mixed-effects association of single cells identifies an expanded effector CD4+ T cell subset in rheumatoid arthritis. Sci Transl Med. 2018;10(463):eaaq0305. https://doi.org/10.1126/scitranslmed.aaq0305.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25. https://doi.org/10.1016/j.cels.2015.12.004.
Raychaudhuri S, Stuart JM, Altman RB. Principal components analysis to summarize microarray experiments: application to sporulation time series. Pac Symp Biocomput. 2000:455–66. https://pubmed.ncbi.nlm.nih.gov/10902193/.
Blanco-Melo D, et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell. 2020;181:1036–1045.e9.
Robinson PC, Liew DFL, Liew JW, Monaco C, Richards D, Shivakumar S, Tanner HL, Feldmann M. The Potential for Repurposing Anti-TNF as a Therapy for the Treatment of COVID-19. Med. 2020;1(1):90–102. https://doi.org/10.1016/j.medj.2020.11.005.
Stoeckius M, Zheng S, Houck-Loomis B, Hao S, Yeung BZ, Mauck WM III, Smibert P, Satija R. Cell Hashing with barcoded antibodies enables multiplexing and doublet detection for single cell genomics. Genome Biol. 2018;19(1):224. https://doi.org/10.1186/s13059-018-1603-1.
Dallagi A, Girouard J, Hamelin-Morrissette J, Dadzie R, Laurent L, Vaillancourt C, Lafond J, Carrier C, Reyes-Moreno C. The activating effect of IFN-γ on monocytes/macrophages is regulated by the LIF-trophoblast-IL-10 axis via Stat1 inhibition and Stat3 activation. Cell Mol Immunol. 2015;12(3):326–41. https://doi.org/10.1038/cmi.2014.50.
Luo W, Li YX, Jiang LJ, Chen Q, Wang T, Ye DW. Targeting JAK-STAT signaling to control cytokine release syndrome in COVID-19. Trends Pharmacol Sci. 2020;41(8):531–43. https://doi.org/10.1016/j.tips.2020.06.007.
Wang BX, Fish EN. Global virus outbreaks: Interferons as 1st responders. Semin Immunol. 2019;43:101300. https://doi.org/10.1016/j.smim.2019.101300.
Davoudi-Monfared E, Rahmani H, Khalili H, Hajiabdolbaghi M, Salehi M, Abbasian L, Kazemzadeh H, Yekaninejad MS. Efficacy and safety of interferon β-1a in treatment of severe COVID-19: A randomized clinical trial. Antimicrobial Agents and Chemotherapy. 2020. https://aac.asm.org/content/64/9/e01061-20.
Nile SH, Nile A, Qiu J, Li L, Jia X, Kai G. COVID-19: pathogenesis, cytokine storm and therapeutic potential of interferons. Cytokine Growth Factor Rev. 2020;53:66–70. https://doi.org/10.1016/j.cytogfr.2020.05.002.
Ye Q, Wang B, Mao J. The pathogenesis and treatment of the `cytokine storm’ in COVID-19. J Inf Secur. 2020;80:607–13.
Cao Y, et al. Ruxolitinib in treatment of severe coronavirus disease 2019 (COVID-19): A multicenter, single-blind, randomized controlled trial. J Allergy Clin Immunol. 2020;146:137–146.e3.
Ahmed A, Merrill SA, Alsawah F, Bockenstedt P, Campagnaro E, Devata S, Gitlin SD, Kaminski M, Cusick A, Phillips T, Sood S, Talpaz M, Quiery A, Boonstra PS, Wilcox RA. Ruxolitinib in adult patients with secondary haemophagocytic lymphohistiocytosis: an open-label, single-centre, pilot trial. Lancet Haematol. 2019;6(12):e630–7. https://doi.org/10.1016/S2352-3026(19)30156-5.
Zizzo G, Cohen PL. Imperfect storm: is interleukin-33 the Achilles heel of COVID-19? Lancet Rheumatol. 2020;2(12):e779–90. https://doi.org/10.1016/S2665-9913(20)30340-4.
Kalil AC, Patterson TF, Mehta AK. Baricitinib plus remdesivir for hospitalized adults with COVID-19. N Engl J Med. 2021;384(9):795–807. https://doi.org/10.1056/NEJMoa2031994.
Karki R, Sharma BR, Tuladhar S, Williams EP, Zalduondo L, Samir P, Zheng M, Sundaram B, Banoth B, Malireddi RKS, Schreiner P, Neale G, Vogel P, Webby R, Jonsson CB, Kanneganti TD. Synergism of TNF-α and IFN-γ triggers inflammatory cell death, tissue damage, and mortality in SARS-CoV-2 infection and cytokine shock syndromes. Cell. 2021;184(1):149–68.
Garcia-Beltran WF, et al. COVID-19-neutralizing antibodies predict disease severity and survival. Cell. 2021;184:476–488.e11.
Zhang F, Mears JR, Shakib L, Beynor JI, Shanaj S, Korsunsky I, Nathan A, Accelerating Medicines Partnership Rheumatoid Arthritis and Systemic Lupus Erythematosus (AMP RA/SLE) Consortium, Donlin LT, Raychaudhuri S. IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation. GSE168710, Gene Expression Omnibus, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168710 (2021).
Zhang F, Mears JR, Shakib L, Beynor JI, Shanaj S, Korsunsky I, Nathan A, Accelerating Medicines Partnership Rheumatoid Arthritis and Systemic Lupus Erythematosus (AMP RA/SLE) Consortium, Donlin LT, Raychaudhuri S. IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation. Github, https://github.com/immunogenomics/inflamedtissue_covid19_reference (2021).
We thank the Brigham and Women’s Hospital Single Cell Genomics Core for assistance in the single-cell hashing experiment. We thank members of the Raychaudhuri Laboratory for discussions.
Accelerating Medicines Partnership Rheumatoid Arthritis & Systemic Lupus Erythematosus (AMP RA/SLE) Consortium:
Jennifer Albrecht9, Jennifer H. Anolik9, William Apruzzese5, Brendan F. Boyce9, Christopher D. Buckley10, David L. Boyle11, Michael B. Brenner5, S. Louis Bridges Jr12, Jane H. Buckner13, Vivian P. Bykerk7, Edward DiCarlo14, James Dolan15, Andrew Filer10, Thomas M. Eisenhaure4, Gary S. Firestein10, Susan M. Goodman7, Ellen M. Gravallese5, Peter K. Gregersen16, Joel M. Guthridge17, Nir Hacohen4, V. Michael Holers18, Laura B. Hughes12, Lionel B. Ivashkiv19,20, Eddie A. James13, Judith A. James17, A. Helena Jonsson5, Josh Keegan15, Stephen Kelly21, Yvonne C. Lee22, James A. Lederer15, David J. Lieb4, Arthur M. Mandelin II22, Mandy J. McGeachy23, Michael A. McNamara7, Nida Meednu9, Larry Moreland23, Jennifer P. Nguyen15, Akiko Noma4, Dana E. Orange24, Harris Perlman22, Costantino Pitzalis25, Javier Rangel-Moreno9, Deepak A. Rao5, Mina Ohani-Pichavant26,27, Christopher Ritchlin9, William H. Robinson26,27, Karen Salomon-Escoto28, Anupamaa Seshadri15, Jennifer Seifert18, Darren Tabechian9, Jason D. Turner10, Paul J. Utz26,27, Kevin Wei5.
9Division of Allergy, Immunology and Rheumatology, Department of Medicine, University of Rochester Medical Center, Rochester, NY, USA.
10Rheumatology Research Group, Institute for Inflammation and Aging, NIHR Birmingham Biomedical Research Center and Clinical Research Facility, University of Birmingham, Queen Elizabeth Hospital, Birmingham, UK.
11Department of Medicine, Division of Rheumatology, Allergy and Immunology, University of California, San Diego, La Jolla, CA, USA.
12Division of Clinical Immunology and Rheumatology, Department of Medicine, Translational Research University of Alabama at Birmingham, Birmingham, AL, USA.
13Translational Research Program, Benaroya Research Institute at Virginia Mason, Seattle, WA, USA.
14Department of Pathology and Laboratory Medicine, Hospital for Special Surgery, New York, NY, USA.
15Department of Surgery, Brigham and Women’s Hospital and Harvard Medical School,
Feinstein Boston, MA, USA.
16Feinstein Institute for Medical Research, Northwell Health, Manhasset, NY, USA.
17Department of Arthritis & Clinical Immunology, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA.
18Division of Rheumatology, University of Colorado School of Medicine, Aurora, CO, USA.
19Graduate Program in Immunology and Microbial Pathogenesis, Weill Cornell Graduate School of Medical Sciences, New York, NY, USA.
20David Z. Rosensweig Genomics Research Center, Hospital for Special Surgery, New York, NY, USA.
21Department of Rheumatology, Barts Health NHS Trust, London, UK.
22Division of Rheumatology, Department of Medicine, Northwestern University Feinberg School of Medicine, Chicago, IL, USA.
23Division of Rheumatology and Clinical Immunology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA.
24The Rockefeller University, New York, NY, USA.
25Centre for Experimental Medicine & Rheumatology, William Harvey Research Institute, Queen Mary University of London, London, UK.
26Division of Immunology and Rheumatology, Department of Medicine, Stanford University School of Medicine, Palo Alto, CA, USA.
27Immunity, Transplantation, and Infection, Stanford University School of Medicine, Stanford, CA, USA.
28Division of Rheumatology, Department of Medicine, University of Massachusetts Medical School, Worcester, MA, USA.
This work is supported in part by funding from the National Institutes of Health (NIH) Grants UH2AR067677, U01HG009379, and R01AR063759 (to S.R.) and NIH R01AI148435, UH2 AR067691, Carson Family Trust, and Leon Lowenstein Foundation (to L.T.D.).
Ethics approval and consent to participate
Healthy blood samples were purchased from the New York Blood Center (NYBC), provided by volunteer donors who consented for the blood to be used in biomedical research and other uses at the discretion of NYBC. The samples are deidentified by the NYBC, the research study investigators had no access to identifiable private information. As per the NIH guidelines, this does not constitute Human Subjects research. For the stimulated blood-derived macrophage experiment, co-cultures with synovial fibroblast involved synovial fibroblast lines generated from patients with RA undergoing arthroplasty (HSS IRB 14-033). Patients provided informed consent and all appropriate measures were taken for compliance with the Helsinki Declaration.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Basic information and demography of multiple single-cell datasets.
Overall integration of immune cells from multiple scRNA-seq datasets. Figure S2. Quantification of the performance of all cell type multi-disease tissue integration. Figure S3. Tissue-level macrophage integrative analysis of multiple scRNA-seq datasets. Figure S4. Heterogeneity of shared inflammatory macrophages from multiple tissues. Figure S5. Single-cell differential gene expression analysis of comparing inflammatory macrophages with non-inflammatory macrophages within each individual tissue source. Figure S6. Examination of the CXCL10+ CCL2+ macrophage marker genes in additional diseased cohort studies. Figure S7. Experimental design and quality control of human blood-derived macrophages stimulated by different conditions. Figure S8. Integrative analysis of tissue-level macrophages and human blood-derived macrophages. Figure S9. Assessment of previously reported stimulated macrophage spectrum analysis and alignment of macrophages from different disease tissues to a trajectory.
Cell type marker genes and statistics.
Number of cells per cluster, per disease and tissue for macrophage integration analysis.
Macrophage cluster marker genes and relative statistics.
Hashtag antibodies for the 10X single-cell cell hashing experiment.
Details for the 10X single-cell cell hashing experiment.
About this article
Cite this article
Zhang, F., Mears, J.R., Shakib, L. et al. IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation. Genome Med 13, 64 (2021). https://doi.org/10.1186/s13073-021-00881-3
- Single-cell transcriptomics
- Single-cell multi-disease tissue integration
- Inflammatory diseases
- Macrophage stimulation
- Macrophage heterogeneity