Abstract
A complete diagnostic autopsy is the gold-standard to gain insight into Coronavirus disease 2019 (COVID-19) pathogenesis. To delineate the in situ immune responses to SARS-CoV-2 viral infection, here we perform comprehensive high-dimensional transcriptional and spatial immune profiling in 22 COVID-19 decedents from Wuhan, China. We find TIM-3-mediated and PD-1-mediated immunosuppression as a hallmark of severe COVID-19, particularly in men, with PD-1+ cells being proximal rather than distal to TIM-3+ cells. Concurrently, lymphocytes are distal, while activated myeloid cells are proximal, to SARS-CoV-2 viral antigens, consistent with prevalent SARS-CoV-2 infection of myeloid cells in multiple organs. Finally, viral load positively correlates with specific immunosuppression and dendritic cell markers. In summary, our data show that SARS-CoV-2 viral infection induces lymphocyte suppression yet myeloid activation in severe COVID-19, so these two cell types likely have distinct functions in severe COVID-19 disease progression, and should be targeted differently for therapy.
Similar content being viewed by others
Introduction
Since the COVID-19 outbreak in January 2020, researchers have identified that immune dysregulation with lymphopenia and cytokine storm are hallmarks of severe disease1. However, much remains unknown about how this virus interacts with the human immune system. Most attempts to profile immune responses have been performed with peripheral blood samples2, but direct pathological evidence is critical to delineate the underlying pathophysiological mechanisms. Pathological reports are still limited by small sample sizes owing to onerous safety requirements and insufficient political support for performing autopsies3. Moreover, many prior studies have only utilized low throughput techniques like traditional immunohistochemistry (IHC) and PCR (Supplementary Table 1): these methods are inadequate for obtaining a complete picture of the heterogeneous immune landscape.
We performed a comprehensive, high-dimensional analysis of the immune signature in 22 decedents (53–88 years old) who had contracted the disease in Wuhan, China and succumbed to COVID-19 after experiencing a severe infection (Supplementary Table 2). We performed transcriptome analyses, Digital Spatial Profiling (DSP), and multiplex-IHC (mIHC) of autopsy tissue samples collected from patients admitted to the Huoshenshan Hospital of Wuhan, China in March 2020.
Here, using a combination of transcriptome analyses with DSP and mIHC, we report the high-dimensional immunopathological profiling of severe COVID-19 patients in multiple organs, including various immune activation-related markers (CD69, IFN-γ, granzyme B, etc.) and immune suppression and exhaustion related markers (TIM-3, PD-1, BTLA, NKG2A, etc.). We highlight the distinct patterns of severe immunosuppression in multiple organs and find TIM-3 and PD-1-mediated immunosuppression is a hallmark of severe COVID-19, and that the extent of immunosuppression correlates with male sex and advanced age. Particularly, activated myeloid cells, but not lymphocytes, are proximity to SARS-CoV-2 viral antigen, consistent with prevalent SARS-CoV-2 infection of myeloid cells in multiple organs. Our data show the COVID-19 infection landscape across multiple organs and demonstrate the viral load is positively associated with ACE2 expression as well as expression of markers of immunosuppression. In conclusion, we suggest that SARS-CoV-2 viral infection induces lymphocyte suppression yet myeloid activation in severe COVID-19. We should thus consider the two immune compartments separately when investigating disease progression.
Results
Severe COVID-19 infection is characterized by immunosuppression
To examine the associations between COVID-19 disease and patient characteristics, we performed bulk RNA-sequencing on lung samples from 11 COVID-19 decedents and 3 non-COVID-19 decedents. We detected clear differences in the transcriptomic profiles between COVID-19 decedents and non-COVID-19 decedents (Fig. 1a–b, Supplementary Fig. 1a–b). By gene set enrichment analysis (GSEA), we determined that the corresponding enriched pathways involved negative regulation of the immune response and reduced leukocyte activation in COVID-19 decedents (Fig. 1c). Interestingly, transcriptomic profiles of the lung tissues were more similar within than between the sexes (Fig. 1d, Supplementary Fig. 1c), suggesting that sex might influence the immune landscape in COVID-19 infection. This trend of sex differences has been further observed in the liver, kidney, spleen, and heart tissues (Fig. 1d). In particular, the immunosuppression signature was more pronounced in men compared to women (Fig. 1e–f) across all organs under consideration and in this subset of decedents translated to a shorter survival duration in men (Supplementary Fig. 1d).
TIM-3 and PD-1 mediate immunosuppression in men
To identify the sex-biased molecular mediators of immunosuppression in COVID-19 infection, we transcriptionally profiled the in situ immune environments of each sample (Fig. 2a). Corroborating the findings from bulk RNA-sequencing, we found that lung tissue from men expressed higher levels of inhibitory receptors, including HAVCR2 (TIM-3), PDCD1 (PD-1), ENTPD1 (CD39), BTLA and KIRs compared to women, but lower levels of immune activation-related markers, such as MKI67 (Ki-67), and CD69 (Fig. 2b–c). These trends were partially recapitulated in the heart and kidney but different from what was seen in the liver (Supplementary Fig. 2a–b). To validate the immunosuppressive signature at the protein level, we performed mIHC and similarly found that higher levels of TIM-3+ and PD-1+TIM-3+ cells were observed in the lung and kidney tissue of men compared to women but that the reverse was relationship observed in the liver (Fig. 2d, Supplementary Fig. 2c–d). In addition, we observed a positive association between lung TIM-3+ cells and age in men but not in women (Supplementary Fig. 2e).
Because clustering of multiple suppressive markers can compound immunosuppression4, we scrutinized the microenvironment surrounding these TIM-3+ cells by investigating their spatial relationship with other immunosuppressive markers. Interestingly, we observed that in the lung, liver, and kidney, the PD-1+ cells were significantly more likely to be proximal than distal from TIM-3+ cells (Fig. 2e–f). These findings suggest an association between TIM-3 and PD-1 upregulation and an unfavorable clinical outcome, and the potential use of anti-TIM-3 and/or anti-PD-1 therapy in combating COVID-19.
Lymphocytes are selectively suppressed in COVID-19
To identify the immune compartment affected by TIM-3 and PD-1-mediated immunosuppression, we separately compared CD3+- and CD68+ -rich regions (Fig. 3a, Supplementary Fig. 3a). In CD3+ -rich regions of the lung, we found a significantly higher expression of immunosuppression-associated genes and various immunosuppressive cytokines in men than in women (Fig. 3b–c). We observed a similar trend in the CD68+-rich region, although it did not statistically differ between the sexes (Supplementary Fig. 3b). We also found higher expression of genes associated with innate immune cells (CD83, NKG7 and CD68) and natural killer (NK) cell activation (GZMB, KLRK1, KLRF1, KLRC2 and NCR1) in the lung CD3+-rich regions of men compared to women (Supplementary Fig. 4a). Thereafter, to verify the gender biased immunosuppression in other organs, we compared the expression of genes associated with immune cells, immune activation, cytokine (Supplementary Fig. 4b–e) and immunosuppression (Supplementary Fig. 5a) from CD3+- and CD68+ -rich regions in the liver, kidney, heart and spleen. We observed similar trends in the CD3+- and CD68+-rich region in the heart, spleen and kidney, but not in the liver (Supplementary Fig. 5a). The difference in immunosuppression-associated gene expression such as PDCD1 in CD3+-rich regions was also more significant than the difference in CD68+-rich region in the heart, spleen and kidney, but not in the liver (Supplementary Fig. 5b–e).
We next asked if the immunosuppressive profile of COVID-19 decedents would manifest in the microenvironment surrounding SARS-CoV-2 viral antigen+ cells. To this end, we characterized the spatial relationship between SARS-CoV-2 nucleocapsid protein (NP) and surrounding cell types. We observed that activated GZMB+CD3+ cells did not exhibit any association with proximity to viral NP in the lung, kidney, heart, and spleen but were associated with proximity to viral NP in the liver (Fig. 3d–e). In lung and spleen tissue, NK cells and CD8+ cells were less likely to be proximal to viral NP whereas the reverse relationship was observed for the liver and heart (Fig. 3d–e). Together, these suggest that infected cells can escape detection by GZMB+CD3+ cells, NK cells, and CD8+ T-cells, particularly in the lung.
Notably, the myeloid population behaved differently as activated myeloid cells [CD80/83+ CD206+ macrophages and BDCA2+/BATF3+ dendritic cells (DCs)] were significantly associated with proximity to viral NP+ cells in the lung, liver, heart, and kidneys (Fig. 3f–g). Furthermore, the myeloid cells proximal to viral NP+ cells were more likely to be activated than those distal from viral NP+ cells across most organs (Fig. 3h). These findings indicate that severe COVID-19 might selectively suppress lymphocytes, particularly in the lung, while myeloid cells retain their capacity to mount an immune response.
ACE2 expression is associated to higher infection rate
Having established selective lymphocyte suppression in situ during severe COVID-19, we sought to identify additional correlates of viral load by stratifying the regions of interest based on the expression levels of SARS-CoV-2 ORF1ab. In doing so, we validated that PDCD1, BTLA, TIGIT and CD96 were highly expressed in virus-high regions, especially in the lungs, livers and spleens (Fig. 4a–b, Supplementary Fig. 6), lending support to our hypothesis that viral presence drives selective lymphocyte immunosuppression in situ. Reflecting its role in viral infection, ACE2 expression was higher in virus-high regions compared to virus-low regions of the lung, liver and spleen (Fig. 4b). This trend of expression is not significantly different in the heart and kidneys. We next characterized the distribution of viral antigens amongst various cell types and determined that in the lung, liver, and kidney and spleen, the predominant NP+ cells were DCs (BATF3/BDCA2+) and furthermore that macrophages (CD206+) constituted an appreciable fraction of NP+ cells only in the liver (Fig. 4c). We also discovered that BATF3/BDCA2+ DCs exhibited a higher tendency for infection over BATF3/BDCA2– cells in all organs under consideration (Fig. 4e) and that only in the liver, CD206+ macrophages were more likely to be infected than CD206– cells while in the lungs and heart, the converse was true (Fig. 4d). In the liver, BATF3/BDCA2+ DCs that were ACE2+ exhibited an even higher likelihood of infection than ACE2– BATF3/BDCA2+ DCs (Fig. 4d–e).
In sum, we characterized the COVID-19 infection landscape across multiple organs and found ACE2 expression to be correlated with infection susceptibility, and viral load to be positively associated with ACE2 expression as well as expression of markers of immunosuppression.
Discussion
COVID-19 infection elicits robust T-cell responses;5 however, poor clinical outcomes have been associated with lymphopenia1 and an increased number of non-functional T cells6. These findings expose the central role of the adaptive immune response in viral clearance. In this study, we suggest that beyond lymphopenia, selective lymphocyte immunosuppression drives disease pathology, particularly in men with severe COVID-19.
The higher morbidity and mortality rates described in men1, 2, 7, 8 follows a common theme in viral infections. Due to biological differences involving estrogen regulation of X-linked genes, incomplete X-inactivation, and modulation of interferon production by sex steroids9, 10, men may experience a higher case fatality rate than women. Based on the results of our study, we propose male-biased immunosuppression as an additional possible mechanism that accelerates disease severity. Specifically, we found in the tissue context that this mechanism operates via immune checkpoint markers on T cells and NK cells but not on DCs and macrophages.
Previous reports have described upregulated CD39, PD-1, BTLA or TIM-3 expression in T cells and NK cells in the peripheral blood of severe COVID-19 patients4, 11,12,13,14,15,16. Complementing these reports within the tissue context, we found in the organs autopsied that viral load correlated with immunosuppressive marker expression and that the immunosuppressive signature is more pronounced in men than women. We additionally identified a greater clustering of PD-1+ cells around TIM-3+ cells in men than women, which compounds the degree of inhibition of lymphocytes, as multiple marker expression parallels suppression severity4. Furthermore, severe COVID-19 has also been associated with non-functional T cells6 and a less expanded T-cell population17. As we detected a higher abundance of lymphocytes distant from rather than close to infected cells, the functionality of T cells is likely to be reined in by local virus-induced expression of immunosuppressive markers.
Taken together, these observations suggest consideration of the use of immunotherapy strategies such as co-blockade of PD-1 and TIM-3, especially in men, to reinvigorate the immune system against COVID-19 infection. One argument against the use of such checkpoint inhibitors is that it might upset the balance between viral clearance and immunopathology caused by the excessive T-cell response18. Several case reports involving cancer patients noticed severe COVID-19 development alongside treatment with checkpoint inhibitors18,19,20,21. However, with cancer as a co-morbidity, these patients are more susceptible to disease complications22 and therefore these reports remain inconclusive until larger studies are conducted17. Mounting evidence also suggests that the timing of administration of therapeutics might be a key determinant of immunopathological response23.
A common complication of severe COVID-19 is that of cytokine release syndrome, in which plasma levels of cytokines such as IL-2, IL-6, IFN-γ, and TNF-α are elevated24, 25. In the context of COVID-19, this phenomenon has been attributed to the myeloid rather than the lymphoid compartment14, also known as macrophage activation syndrome (MAS). In contrast to lymphocytes, we saw that DCs responded normally to infected cells and were activated in their vicinity. DCs also exhibited a higher tendency to be infected, and myeloid cells were the primary infected cell type among the immune compartments. This observation lends support to the hypothesis that myeloid cells home to the lung26, 27, get infected, and subsequently serve as a conduit for viral spread to other organs28, 29. Although viral replication within the cell might be limited, viral infected myeloid cells have the capacity to induce a great variety of pro-inflammatory cytokines and chemokines, which contribute to MAS30. Therefore, as a result of exposure to viral antigens, DCs and macrophages proximal to viral antigens have a higher tendency to be activated than distal DCs and macrophages. Ultimately, MAS corresponds to our model whereby such differential contributions originate from the greater proximity of myeloid cells than lymphocytes to infected cells.
Finally, myeloid cells such as monocytes, macrophages and DCs have been reported to express ACE230,31,32. As such, myeloid cells are indeed susceptible to direct SARS-CoV-2 infection. In particular, multiple IHC studies have clearly identified macrophages in the alveolar septum and lung cavity33,34,35,36,37,38. Positive staining for SARS-CoV-2 viral antigen on ACE2+ macrophages but not T cells were reported29, as well as flow cytometry analysis indicating SARS-CoV-2 spike protein interaction with ACE2+ macrophages, but not T cells32. It has been reported that monocytes undergo morphological and phenotypic alterations during COVID-19 infection39, 40, however the effects on DCs are less studied31. In this study, we observed that BATF3+BDCA2+ DCs were less abundant in women than in men. We also described an increased likelihood of ACE2+ cells to be infected over ACE2– cells. While ACE2 has been well-characterized as the receptor for the viral spike protein41, to our knowledge, ours is the first study to relate ACE2 expression to infection rate. Reports have suggested that the decreased COVID-19 severity in women might be related to the ACE2 locus being on the X chromosome42, however the interplay between DC abundance, ACE2 expression and the risk of cell infection and consequent disease severity remain to be further explored.
Despite our high-dimensional analysis, one limitation of our study is the sample size of 22 decedents, which given the rising global death toll, may now seem like a small number. Our samples originated from the main cohort of the limited number of decedents autopsied during the first outbreak of this pandemic in Wuhan, China. In the initial phase, autopsies were not favored due to onerous safety requirements and insufficient political support3. Hence, as a result of the sample size limitation, some statistical tests are under-powered, and we therefore sought to mitigate this shortcoming by profiling the samples through multiple modalities.
While the vast majority of studies that profiled COVID-19 immune responses have done so with peripheral blood samples2 and low throughput techniques (Supplementary Table 1) which only gather information from a single tissue slice, these are inadequate in obtaining a complete picture of the heterogenous immune landscape and lack the ability to delineate the underlying pathophysiological mechanisms of COVID-19. By comparison, our sample size is comparable to the current state of published COVID-19 autopsy reports, where only 2 multi-organ autopsies using IHC have included more than 10 patients (Supplementary Table 1). Finally, to our knowledge ours is the first study of the in situ immune microenvironment in different tissues during COVID-19 disease.
By combining transcriptome analyses with DSP and mIHC, we have achieved high-dimensional immunopathological profiling of COVID-19. From the resulting data, we propose that severe COVID-19 incites robust immune responses from myeloid cells while selectively suppressing lymphocytes, recapitulating in the tissue context COVID-19-associated observations such as MAS and lymphopenia. The distinct patterns of severe immunosuppression in multiple organs suggest that TIM-3 and PD-1-mediated immunosuppression is a hallmark of severe COVID-19, and that the extent of immunosuppression correlates with male sex and advanced age. At the protein level, PD-1+ cells can be found in close proximity to TIM-3+ cells. Activated myeloid cells, but not lymphocytes, were observed in proximity to SARS-CoV-2 viral antigen, thus providing direct pathological evidence for virus-induced myeloid cell activation that likely drives disease progression. Going forward, we propose that the two immune compartments should be considered separately when investigating disease progression.
Methods
Study cohort
The autopsy samples from 22 patients with severe COVID-19 who had undergone autopsy were obtained from the Huoshenshan Hospital, Wuhan, China. All safety precautions were in line with recently published guidelines43. All patients provided informed consent and voluntary donation certificates. The age of the patients ranged from 51 to 88 years, with a median age of 66 years. The three lung specimens of non-COVID-19 patients used in the bulk RNA-sequencing analysis were obtained from patients suffering from lung squamous cell carcinoma, pneumonia, and pulmonary bronchiectasis, respectively. The age of these patients ranged from 51 to 73 years, with a median age of 64.0 years. The clinical characteristics of all patients are provided in Supplementary Table 2. Bulk RNA-sequencing, Nanostring Digital Spatial Profiling and multiplex-immunohistochemistry were performed on lung, liver, heart, kidney, and salivary gland specimens from these patients. The specific uses of each autopsy sample used in each experiment are provided in Supplementary Table 4. All samples were anonymously coded in accordance with the Helsinki Declaration.
Study approval
This study was approved by the Ethics Committee of the First Affiliated Hospital of Army Medical University (KY2020027) and the First Affiliated Hospital of University of Science and Technology of China (2021KY193). Written informed consent was obtained from all patients. Family members’ consent for donation was obtained from all medicine autopsies on patients with COVID-19 infection. Brain and spinal cord samples were not taken due to the need to use a bone saw and the risk of excess aerosolization.
Bulk RNA sequencing
The raw RNA sequencing data were subjected to vigorous quality control (QC) procedures before downstream analysis. Raw fastq files were processed with Hisat2 software (v7.5.0). FastQC (v0.11.9) were also used to evaluate the sequencing quality. The clean data were compared using Hisat244. The sequence reads were aligned to the human genome hg19 (GRCh38), and the aligned reads were counted within exons using FeatureCounts (v2.0.1) with default parameters to generate the expression matrices of the raw read counts. For the raw data, the rowSums () function in the base package was used that filters out data with read counts < 2. The DESeq2 (v1.26.0) and edgeR (v3.28.1) package in Rstudio (https://rstudio.com/) was used to identify differentially expressed genes. All analyses were performed in R (https://www.r-project.org/;v3.6.1). Heat maps were drawn using Pheatmap(v1.0.12) and the GSEA and GO analysis of the enriched genes was performed using topGO (version2.38.1), clusterProfiler (version3.14.3), pathview (version1.26.0), AnnotationHub (version2.18.0), org.Hs.eg.db (version3.10.0), tidyverse (version1.3.0), GSEABase (version1.48.0), Rgraphviz (version2.30.0), showtext (version0.9), readr (version1.3.1), stats (version3.6.1), enrichplot (version1.6.1). Ggplot2 (version3.3.2), ggrepel (version0.8.2), RcolorBrewer (version1.12) were used drawing and color matching. TrimGalore (https://github.com/FelixKrueger/TrimGalore/issues/25) was used for data QC to remove joints in the pre-processing.
Digital spatial profiling (DSP)
NanoString spatial profiling was used for DSP analysis of three protein-level and 1800 transcriptomic-level immune-markers simultaneously on FFPE tissue slides45, 46. Following deparaffinization and antigen retrieval procedures, sections were simultaneously incubated overnight with fluorescent-labeled antibodies against CD3 (Origene), CD68 (Santa Cruz) and Pan-Cytokeratin (Novus Biologicals), smooth muscle actin (Invitrogen), CD20 (Novus Biologicals), or cytokeratin 8/18 (Novus Biologicals) to visualize morphological features in the regions of interest (ROI). After staining, the tissues were scanned using a GeoMx DSP instrument to generate digital fluorescent images and to select for individual ROIs inside a geometric section. To carry out high-resolution multiplex profiling, each ROI was assigned to CD3+ T-cell-rich or CD68+ macrophage-rich regions. UV light was directed through a programmable digital micro mirror device (DMD) or dual DMD (DDMD) to accurately illumine the ROI and cleave PC oligos from the selected region, which were collected by microcapillary tube inspiration and dispensed into a 96-well plate. Inside the microplate, the single-molecule counting nCounter System enabled the digital counting of released oligos47, 48.
Digital spatial profiling (DSP) analysis
Q3 Normalization – Individual counts were normalized against the 75th percentile of the signal from their own ROI49. We scale our AOIs so that they all have the same value for their Quartile 3 value. In particular, we (i) first, divide all the genes per AOI by their respective Q3 count, and (ii) second, multiply all the genes in all AOIs by a constant, defined as the geometric mean of Q3 counts for all AOIs. This approach generally performs well for CTA data and is often the preferred normalization. Normalize. Quantiles is used to correct the different batch.
Immunohistochemistry (IHC)
IHC was performed on FFPE tissue samples as previously described50,51,52. Lung, liver, heart, kidney, and spleen (4-μm thick) were labeled with antibodies targeting the SARS-CoV-2 NP (Supplementary Table 3). Appropriate controls were included. Images were captured using an IntelliSite Ultra-Fast Scanner (Philips, Eindhoven).
Multiplex immunohistochemistry (mIHC)
mIHC was performed using an Opal Multiplex fIHC kit (Akoya Biosciences, California). FFPE tissue sections were cut onto Bond Plus slides (Leica Biosystems, Richmond) and heated at 60 °C for 20 min. The tissue slides were then subjected to deparaffinization, rehydration and heat-induced epitope retrieval using a Leica Bond Max autostainer (Leica Biosystems, Melbourne) before endogenous peroxidase blocking (Leica Biosystems, Newcastle). Next, the slides were incubated with primary antibodies followed by incubation with polymeric HRP-conjugated secondary antibodies (Leica Biosystems, Newcastle) (Supplementary Table 3). Then, the samples were incubated with Opal fluorophore-conjugated tyramide signal amplification (TSA) (Akoya Biosciences, California) at a 1:100 dilution. The slides were rinsed with wash buffer (BOND Wash Solution 10X Concentrate) after each step. Following TSA deposition, the slides were again subjected to heat-induced epitope retrieval to strip the tissue-bound primary/secondary antibody complexes prior to further labeling. These steps were repeated until the samples were labeled with all six markers and spectral DAPI (Akoya Biosciences, California) at a 1:10 dilution. Finally, the slides were mounted in ProLong Diamond Anti-fade Mountant (Molecular Probes, Life Technologies, USA) and developed in the dark at room temperature for 24 h.
mIHC analysis
Multiple random ROIs were scored for each patient sample. Analyses were performed using Python(v3.8.2) on the cell-segmented output from the HALO software and all mIHC analysis plots were created using the matplotlib (v3.2.1) package. Statistical tests were performed using the scipy (v1.3.2), numpy (v1.18.5), pandas (v1.0.5) and statsmodels (v0.11.1) packages. Populations proximal or distal from “central cell” markers (TIM-3+ or SARS-CoV-2 NP+ in this paper) were determined by calculating a distance matrix between the centroids of each central cell to every other cell and identifying those cells <15 μm from a central cell (except the central cell itself) as proximal and all other cells as distal. ROIs from the same sample were combined by summation of proximal and distal cell counts. ROIs without central cells were disregarded. Odds ratios were aggregated to derive aggregated odds ratios and confidence intervals using Cochran-Mantel-Haenszel and Breslow-Day procedures as implemented by the Stratified Table class of the statsmodels module using default values. Comparison of cell percentages between sexes in Fig. 2d and Supplementary Fig. 2e were performed by averaging of cell percentages across all ROIs of each patient, followed by taking the average of those percentages in men and dividing by the average of those percentages in women, and transforming the resulting by log2. Odds ratios defined as (N(x)NP+/N(y)NP+)/(N(x)NP−/N(y)NP−), where (x) and (y) refer to two cell types compared to each other. When the odds ratio is close to 1, both cell types are as likely to be NP+ in Fig. 4d
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
The RNA-Seq and GeoMx data generated in this study have been deposited in the National Genomics Data Center database under the accession numbers HRA000974,OMIX488 and Gene Expression Omnibus under accession numbers GSE182917, GSE182920. The images can be queried in the database BioStudies under accession number S-BIAD170. Source data are provided with this paper.
Code availability
Code for the analysis of RNA-seq data can be accessed via: https://github.com/peiqi-sudo/Postmortem-high-dimensional-immune-profiling-of-severe-COVID-19-patients-reveals-distinct-patterns-o.
References
Chen, G. et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J. Clin. Investig. 130, 2620–2629 (2020).
Chen Z., John Wherry E. T cell responses in patients with COVID-19. Nat. Rev. Immunol. 20, 529–536 (2020).
Salerno, M. et al. No autopsies on COVID-19 deaths: a missed opportunity and the lockdown of science. J. Clin. Med. 9, 1472 (2020).
Xiong Q. et al. Single-cell immune profiling of matched disease and recovery phase blood of COVID-19 patients reveals anti-LAG3 based immune checkpoint therapy as a promising option. https://doi.org/10.21203/rs.3.rs-24576/v1 (2020).
Westmeier J. et al. Impaired Cytotoxic CD8(+) T cell response in elderly COVID-19 patients. mBio 11, e02243–20 (2020).
Zheng, H. Y. et al. Elevated exhaustion levels and reduced functional diversity of T cells in peripheral blood may predict severe progression in COVID-19 patients. Cell Mol. Immunol. 17, 541–543 (2020).
Meng, Y. et al. Sex-specific clinical characteristics and prognosis of coronavirus disease-19 infection in Wuhan, China: a retrospective study of 168 severe patients. PLOS Pathog. 16, e1008520 (2020).
Richardson, S. et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized With COVID-19 in the New York city area. JAMA 323, 2052–2059 (2020).
Scully, E. P., Haverfield, J., Ursin, R. L., Tannenbaum, C. & Klein, S. L. Considering how biological sex impacts immune responses and COVID-19 outcomes. Nat. Rev. Immunol. 20, 442–447 (2020).
Klein, S. L. & Flanagan, K. L. Sex differences in immune responses. Nat. Rev. Immunol. 16, 626–638 (2016).
Demaria O. et al. Identification of druggable inhibitory immune checkpoints on Natural Killer cells in COVID-19. Cell. Mol. Immunol. 17, 995–997 (2020).
Mathew D. et al. Deep immune profiling of COVID-19 patients reveals distinct immunotypes with therapeutic implications. Science 369, eabc8511 (2020).
Sekine T. et al. Robust T cell immunity in convalescent individuals with asymptomatic or mild COVID-19. bioRxiv 183, 158–168.e14. 2020.2006.2029.174888 (2020).
Wilk A. J. et al. A single-cell atlas of the peripheral immune response to severe COVID-19. medRxiv 26, 1070–1076. 2020.2004.2017.20069930 (2020).
De Biasi, S. et al. Marked T cell activation, senescence, exhaustion and skewing towards TH17 in patients with COVID-19 pneumonia. Nat. Commun. 11, 3434 (2020).
Schultheiss, C. et al. Next-generation sequencing of T and B cell receptor repertoires from COVID-19 patients showed signatures associated with severity of disease. Immunity 53, 442–455 e444 (2020).
Liao, M. et al. Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat. Med. 26, 842–844 (2020).
Kang, C. K. et al. Aberrant hyperactivation of cytotoxic T-cell as a potential determinant of COVID-19 severity. Int. J. Infect. Dis. 97, 313–321 (2020).
Xu, Z. et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respiratory Med. 8, 420–422 (2020).
Bonomi, L., Ghilardi, L., Arnoldi, E., Tondini, C. A. & Bettini, A. C. A rapid fatal evolution of coronavirus disease-19 in a patient with advanced lung cancer with a long-time response to nivolumab. J. Thorac. Oncol. 15, e83–e85 (2020).
Lovly C. M. et al. Rapidly fatal pneumonitis from immunotherapy and concurrent SARS-CoV-2 infection in a patient with newly diagnosed lung cancer. medRxiv https://doi.org/10.1101/2020.04.29.20085738 (2020).
Dai, M. et al. Patients with cancer appear more vulnerable to SARS-CoV-2: a multicenter study during the COVID-19 outbreak. Cancer Disco. 10, 783–791 (2020).
Park, A. & Iwasaki, A. Type I and Type III interferons - induction, signaling, evasion, and application to combat COVID-19. Cell Host Microbe 27, 870–878 (2020).
Huang, C. et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395, 497–506 (2020).
Wang, J., Jiang, M., Chen, X. & Montaner, L. J. Cytokine storm and leukocyte changes in mild versus severe SARS-CoV-2 infection: review of 3939 COVID-19 patients in China and emerging pathogenesis and therapy concepts. J. Leukoc. Biol. 108, 17–41 (2020).
Sanchez-Cerrillo I., et al. COVID-19 severity associates with pulmonary redistribution of CD1c+ DC and inflammatory transitional and nonclassical monocytes. J. Clin. Invest. 130, 6290–6300 (2020).
Peruzzi B. et al. Quantitative and qualitative alterations of circulating myeloid cells and plasmacytoid DC in SARS-CoV-2 infection. Immunology, 161, 345–353 (2020).
Liang Y. et al. Highlight of immune pathogenic response and hematopathologic effect in SARS-CoV, MERS-CoV, and SARS-Cov-2 infection. Front. Immunol. 11, 1022 (2020).
Chen Y. et al. The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) directly decimates human spleens and lymph nodes. medRxiv, https://doi.org/10.1101/2020.03.27.20045427 (2020).
Jafarzadeh, A., Chauhan, P., Saha, B., Jafarzadeh, S. & Nemati, M. Contribution of monocytes and macrophages to the local tissue inflammation and cytokine storm in COVID-19: Lessons from SARS and MERS, and potential therapeutic interventions. Life Sci. 257, 118102–118102 (2020).
Campana P. et al. Dendritic cells and SARS-CoV-2 infection: still an unclarified connection. Cells 9, 2046 (2020).
Wang C. et al. Alveolar macrophage dysfunction and cytokine storm in the pathogenesis of two severe COVID-19 patients. EBioMedicine 57, 102833 (2020).
Yao, X. H. et al. [A pathological report of three COVID-19 cases by minimal invasive autopsies]. Zhonghua Bing Li Xue Za Zhi 49, 411–417 (2020).
Barton L., Duval E., Stroberg E., Ghosh S., Mukhopadhyay S. COVID-19 Autopsies, Oklahoma, USA. Am. J Clin. Pathol. 153 https://doi.org/10.1093/ajcp/aqaa062 (2020).
Zhang Y., Gao Y., Qiao L., Wang W., Chen D. Inflammatory response cells during acute respiratory distress syndrome in patients with coronavirus disease 2019 (COVID-19). Ann. Intern. Med. 173, 402–404 (2020).
Buja, L. M. et al. The emerging spectrum of cardiopulmonary pathology of the coronavirus disease 2019 (COVID-19): report of 3 autopsies from Houston, Texas, and review of autopsy findings from other United States cities. Cardiovasc Pathol. 48, 107233 (2020).
Carsana L. et al. Pulmonary post-mortem findings in a series of COVID-19 cases from northern Italy: a two-centre descriptive study. Lancet Infect Dis https://doi.org/10.1016/S1473-3099(20)30434-5 (2020).
Adachi T. et al. Clinicopathologic and immunohistochemical findings from autopsy of patient with COVID-19, Japan. Emerg Infect Dis 26, 2157–2161 (2020).
Schulte-Schrepping, J. et al. Severe COVID-19 is marked by a dysregulated myeloid cell compartment. Cell 182, 1419–1440.e1423 (2020).
Zhang D. et al. COVID-19 infection induces readily detectable morphological and inflammation-related phenotypic changes in peripheral blood monocytes, the severity of which correlate with patient outcome. medRxiv, https://doi.org/10.1101/2020.03.24.20042655 (2020).
Walls, A. C. et al. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 181, 281–292 e286 (2020).
Martinez F. O., Combes T. W., Orsenigo F., Gordon S. Monocyte activation in systemic Covid-19 infection: assay and rationale. E. Bio. Medicine 59, 102964 (2020).
Hanley, B., Lucas, S. B., Youd, E., Swift, B. & Osborn, M. Autopsy in suspected COVID-19 cases. J. Clin. Pathol. 73, 239 (2020).
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667 (2016).
Gupta, S., Zugazagoitia, J., Martinez-Morilla, S., Fuhrman, K. & Rimm, D. L. Digital quantitative assessment of PD-L1 using digital spatial profiling. Lab. Investig. 100, 1311–1317 (2020).
Merritt, C. R. et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat. Biotechnol. 38, 586–599 (2020).
Beechem, J. M. High-plex spatially resolved RNA and protein detection using digital spatial profiling: a technology designed for immuno-oncology biomarker discovery and translational research. Methods Mol. Biol. 2055, 563–583 (2020).
Toki, M. I. et al. High-plex predictive marker discovery for melanoma immunotherapy-treated patients using digital spatial profiling. Clin. Cancer Res.: Off. J. Am. Assoc. Cancer Res. 25, 5503–5512 (2019).
Beasley G. M. et al. Dissecting the immune landscape of tumor draining lymph nodes in melanoma with high-plex spatially resolved protein detection. Cancer Immunol. Immunother. 70, 475–483 (2020).
Sun, H. et al. Human CD96 correlates to natural killer cell exhaustion and predicts the prognosis of human hepatocellular carcinoma. Hepatology 70, 168–183 (2019).
Sun, C. et al. Oncofetal gene SALL4 reactivation by hepatitis B virus counteracts miR-200c in PD-L1-induced T cell exhaustion. Nat. Commun. 9, 1241 (2018).
Sun, H. et al. Accumulation of tumor-infiltrating CD49a(+) NK cells correlates with poor prognosis for human hepatocellular carcinoma. Cancer Immunol. Res. 7, 1535–1546 (2019).
Acknowledgements
This work was supported by the National Natural Science Foundation of China (#92169118, #8202290021), National Key Research and Development Program of China (#2021YFC2300600), Anhui Provincial Natural Science Foundation (#2008085J35, 1908085MH281), and the Fundamental Research Funds for the Central Universities (WK9110000055).
Author information
Authors and Affiliations
Contributions
C.S. conceived and directed the study. H.B.W., Y.R., S.Q.X., H.L., W.W., Z.B.L., D.Y.Z., J.C., X.D.Z., D.P.J., X.C.F., L.Z., H.Z., Z.H.L., R.C., W.Q.L., C.F.W., S.Y.Z., J.W.Q., B.N. and Z.W. collect autopsy specimens. P.Q.H. and H.B.W. performed and interpreted the IHC data. Y.R. and Z.B. provide IHC technical support. P.Q.H. and H.B.W. performed biostatistical analysis. P.Q.H. analyzed Nanostring DSP profiling. C.S. drafted the manuscript and final approval of all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review information
Nature Communications thanks Xianwen Ren and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wu, H., He, P., Ren, Y. et al. Postmortem high-dimensional immune profiling of severe COVID-19 patients reveals distinct patterns of immunosuppression and immunoactivation. Nat Commun 13, 269 (2022). https://doi.org/10.1038/s41467-021-27723-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-021-27723-5
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.