Efferocytosis of SARS-CoV-2-infected dying cells impairs macrophage anti-inflammatory functions and clearance of apoptotic cells

Macrophages engulf dying cells carrying viable SARS-CoV-2

SARS-CoV-2 infection causes cytopathic effects in human and primate epithelial cells (Chu et al., 2020; Zhu et al., 2020), possibly mediated by activation of multiple cell death pathways (Chan et al., 2020; Shaohua Li et al., 2020; Shufen Li et al., 2020; Mulay et al., 2021; Ren et al., 2020; Zhu et al., 2020). Notably, inhibition of apoptosis ameliorates cytokine expression and tissue damage in the lungs of SARS-CoV-2-infected mice, indicating that apoptosis in the lungs is pathogenic (Chu et al., 2021). Agreeing with previous reports (Chu et al., 2020; Ren et al., 2020; Zhu et al., 2020), infection with SARS-CoV-2 for 48 hours induced caspase-3 and caspase-8 activation in immortalized simian kidney epithelial Vero CCL81 (Figure 1A) and human lung epithelial Calu-3 cell lines (Figures 1A and B). We also found features consistent with the predominance of apoptotic cell death at 48 h post-infection by flow cytometric analysis of the pattern of annexin V and a permeability dye co-staining in infected cells (Figure 1C). Consistently, we observed low levels of cytolysis in these infected cells, as assessed by lactate dehydrogenase release (Figure 1D). Immunofluorescent labeling of cleaved caspase-3 in lung tissues of deceased COVID-19 patients confirmed the induction of apoptosis in the epithelia infected with replicant SARS-CoV-2 (Figure 1E, and Figure 1 – figure supplement 1A). Examination of lung tissues also revealed evidence of macrophages with internalized epithelial cells and SARS-CoV-2, suggesting the uptake of infected epithelial cells by macrophages (Figure 1F).

• Download figure
• Open in new tab

Figure 1 – figure supplement 1. Support data for Figure 1.

(A). Immunofluorescence images of tissue samples immunolabeled with secondary antibodies-only used for detection anti-cleaved Caspase-3 (red, anti-rabbit IgG conjugated to Alexa 594), anti-dsRNA (green, donkey anti-mouse IgG F(ab)’2-Alexa 488), and stained with Hoechst (nuclei, blue) in Figure 1E. Tissues were scanned by wide-field epifluorescence imaging. Scale bar: 20 μm, asterisks: erythrocytes.

(B) TCID₅₀ quantification of SARS-CoV-2 viral loads from detached, apoptotic Calu-3 and Vero CCL81 cells (CoV2-AC, 1.0 x 10⁶ cells) and their cell-free supernatants, collected 48 h post-infection and isolated as described in the Methods.

(C) Uptake of apoptotic Calu-3 and Vero CCL81 cells (labeled with CFSE) in response to UV irradiation (UV-AC) or infection with SARS-CoV-2 for 48 h (Cov2-AC) by THP-1-derived macrophages (labeled with CTFR). Macrophages were co-incubated with AC for 2 h, and internalization was assessed by flow cytometry. Representative gating strategy, plots, and percentage of engulfment (gated on Total cells/Single cells/Live cells/CTFR⁺/CTFR⁺ CFSE⁺ cells) are shown.

(D) Representative maximal projection of scanning confocal images showing the uptake of UV-AC or CoV2-AC (from Vero CCL81 cells, CFSE-labeled, green) by human monocyte-derived macrophages after 2 h of co-incubation. Samples were stained with Hoechst (nuclei, blue) and were immunolabeled with anti-Spike (SARS-CoV-2, red) and anti-CD11b (macrophage, white). Scale bar: 5 μm.

(E) Human monocyte-derived macrophages or THP-1-derived macrophages were unstimulated (Cntrl), co-incubated with apoptotic Vero CCL81 cells isolated from UV-irradiated (UV-AC) or SARS-CoV-2-infected (CoV2-AC) cultures, stimulated with the supernatants of infected apoptotic cells (CoV2-AC Sup), or infected with SARS-CoV-2 (CoV2), as indicated. The percentage of cellular viability was measured using a viability probe 24 h post-stimulation and analysis by flow cytometry (gated on Total cells/Single cells).

Boxes represent the mean of three biological replicates, and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by Student’s test (B-C) or ANOVA (E).

*, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups; ns: non-significant. Data shown are from one representative out of two (B, C) or three (E) experiments performed independently with similar results.

• Download figure
• Open in new tab

Figure 1. Macrophages engulf apoptotic cells carrying viable SARS-CoV-2

(A) Immunoblot of caspase-8 and caspase-3 cleavage in Calu-3 and Vero CCL81 cells in response to infection with SARS-CoV-2 for 48 h (CoV2, detected with rabbit anti-Spike), UV irradiation (UV) or freezing and thawing (F/T).

(B) Flow cytometric analysis of caspase-3 cleavage (Cl-caspase-3) in Calu-3 cells unstimulated (Cntrl) or infected with CoV2 for 48 h. The pan-caspase inhibitor z-VAD-FMK (zVAD; 20 µM) or vehicle (DMSO) were added to cell cultures after viral adsorption. Representative histogram, geometric mean fluorescence intensity (MFI), and frequency of Cl-caspase-3⁺ cells (gated on Total cells/Single cells) are shown. The black line on the histogram represents the control sample labeled with secondary antibody only.

(C) Flow cytometric analysis of apoptosis in Vero CCL81 cells in response to UV light irradiation or 48 h post-infection with CoV2, assessed by annexin V (FITC) and viability dye (Zombie NIR) co-staining. Representative histogram and frequency of each population (pre-gated on Total cells/Single cells) are shown.

(D) Cytotoxicity analysis by measurement of lactate dehydrogenase (LDH) release on the supernatants of Vero CCL81 cells in response to UV light irradiation or 48 h post-infection with CoV2.

(E, F) Representative histological findings in post-mortem lung tissue from COVID-19 patients obtained by ultrasound-guided minimally invasive autopsy. (E) Immunofluorescence images of tissue samples labeled with anti-Cl-caspase-3 (red), anti-dsRNA (green), and stained with Hoechst (nuclei, blue) for the detection of caspase-3 activation in infected epithelia (COVID-19, bottom panels) compared to samples from control patients (Cntrl, upper panels). Tissues were scanned by wide-field epifluorescence imaging. Scale bar: 20 μm, asterisks: erythrocytes. The mean fluorescence intensity (MFI) of Cl-caspase-3 of at least 300 cells are shown (n = 3 individuals per group). Technical control for secondary antibodies background is shown in Figure S1A. (F) Representative immunofluorescence image of a COVID-19 patient tissue sample labeled for the detection of efferocytosis in situ and superimposed with differential interference contrast. Samples were labeled with anti-CD68 (detecting macrophages, green), anti-cytokeratin18 (CK-18, detecting epithelial cells, white), anti-Spike (detecting SARS-CoV-2, red), and stained with Hoechst (nuclei, blue). Bottom panels show higher magnification of (a, b) CD68⁺ macrophage with CK-18 and Spike labeling in the cytosol and (c) CD68⁺ only macrophage. White arrows: CK-18⁺ cells; asterisks: erythrocytes. Tissues were scanned by wide-field epifluorescence imaging. Scale bar: 20 μm; scale bar for insets: 10 μm.

(G) Quantification by TCID₅₀ of SARS-CoV-2 viral loads in Vero CCL81 cells infected for 48 h. Viral loads were estimated for pre-sorted (Total cells) or annexin V⁺ (Ann⁺ cells) cells isolated by magnetic separation.

(H) Uptake of apoptotic Vero CCL81 cells (labeled with CTFR) in response to UV irradiation (UV-AC) or infection with SARS-CoV-2 for 48 h (Cov2-AC) by human monocyte-derived macrophages (labeled with CTV). Macrophages were co-incubated with AC for 2 h, and internalization was assessed by flow cytometry. Representative plot and percentage of engulfment (gated on Total cells/Single cells/Live cells/CTV⁺/ CTV⁺CTFR⁺ cells) are shown.

Boxes represent the mean of five (B, D), four (C, H), or three (G) biological replicates, and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA (B-D), Mann-Whitney test (E) or Student’s test (G, H). *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups. Data shown are from one representative out of two (A, C, G) or three (B, D, H) experiments performed independently with similar results.

We performed the Median Tissue Culture Infection Dose (TCID₅₀) assay and found that viral particles obtained from isolated apoptotic Vero and Calu-3 cells induced cytopathic effect to an equivalent extent as those released in the culture supernatants (Figure 1 – figure supplement 1B). Furthermore, we isolated annexin V-labeled apoptotic cells (AC) from SARS-CoV-2-infected cell culture and confirmed that they carry viable viral particles Figure 1G). Therefore, viable SARS-CoV-2 is retained within infected, apoptotic epithelial cells.

Exposure of phosphatidylserine (PtdSer) on the outer plasma membrane of cells undergoing regulated cell death (RCD), as observed by annexin V binding to infected dying cells (Figure 1C), is the most ubiquitous signal that triggers their phagocytosis (Fadok et al., 1992; Nagata, 2018). To determine whether phagocytes efficiently engulf SARS-CoV-2-infected apoptotic cells (CoV2-AC), we collected the loosely attached apoptotic cells from infected cultures and used flow cytometry to assess their uptake by primary human macrophages and macrophages differentiated from THP-1 monocytes. We observed that macrophages efficiently engulf CoV2-AC, similarly to the uptake of sterile, UV-irradiated apoptotic cells (UV-AC) (Figure 1H and Figure 1 – figure supplement 1C). Analysis by confocal microscopy and staining for Spike protein confirmed the presence of SARS-CoV-2 in engulfed cell corpses (Figure 1 – figure supplement 1D). Finally, either stimulation with infected or sterile dying cells did not robustly affect macrophages viability up to 24 h post-treatment (Figure 1 – figure supplement 1E).

These results demonstrate that macrophages engulf cell corpses carrying viable particles of SARS-CoV-2, offering a framework to investigate their effects on macrophage function.

Efferocytosis of SARS-CoV-2-infected dying cells impairs macrophage anti-inflammatory function

We next addressed the effect of the engulfment of SARS-CoV-2-infected apoptotic cells in macrophages. The uptake of sterile apoptotic cells often promotes macrophage functional polarization toward an anti-inflammatory and tissue repair phenotype (Doran et al., 2020). We first observed that stimulation of primary monocyte-derived macrophages with CoV2-AC, but not infection with SARS-CoV-2, reduced the expression of genes associated with alternative programming to tissue remodeling and secretion of immune-modulatory mediators such as CCL18, CD206 (also known as MRC1), MMP9, PPARG and CD163 (Figure. 2A). Furthermore, stimulation with UV-AC obtained from infected Vero cells increased both gene (Figure 2A) and surface protein (Figure 2B, Figure 2 – figure supplement 1A) expression of the anti-inflammatory marker CD206. However, upregulation of MRC1 at the gene or protein levels did not occur in response to stimulation with CoV2-AC (Figures 2A and B). We obtained similar results in THP-1-derived macrophages stimulated with CoV2-AC obtained from either Calu-3 or Vero cells infected with SARS-CoV-2 (Figure 2 – figure supplement 1B-D). Stimulation with the conditioned supernatants of infected dying cells (containing putative DAMP, cytokines produced by epithelial cells, extracellular vesicles, and released virions) or direct infection with SARS-CoV-2 did not induce MRC1 transcription (Figure 2A, Figure 2 – figure supplement 1B) or surface expression (Figure 2B, Figure 2 – figure supplement 1C and D) in macrophages. We also confirmed that modulation of CD206 expression occurred specifically in the macrophages that engulfed apoptotic cells and was not due to paracrine signaling (Figure 2C and Figure 2 – figure supplement 1E). To determine if viral viability was required to modulate CD206 expression by infected apoptotic cells, isolated CoV2-AC were exposed to UV irradiation for virus inactivation (Figure 2 – figure supplement 1F). Notably, macrophages treated with epithelial AC containing inactivated SARS-CoV-2 exhibited increased CD206 surface expression, similar to that of UV-AC (Figure 2D). This finding suggests that viable SARS-CoV-2 actively represses alternative programming carried in dying cells and is not primarily enforced by the type of cell death caused by the infection. Efferocytosis of sterile UV-AC in the presence of SARS-CoV-2 still induced higher expression of CD206 in macrophages (Fgure 2E). Therefore, the suppression of CD206 surface expression required delivery of viral particles within cell corpses to macrophages.

• Download figure
• Open in new tab

Figure 2 – figure supplement 1. Support data for Figure 2.

(A) Representative gating strategy used to determine the CD206 expression on the surface of macrophages in response to the different stimulus by flow cytometry (backgating corresponds to FMO control showed in Figure 2B).

(B) THP-1-derived macrophages were unstimulated (Cntrl), co-incubated with UV-AC, CoV2-AC, or CoV2-AC Sup obtained from Vero CCL81 cells, or infected with CoV2, as indicated. Expression of MRC1 in macrophages showed as fold change relative to Cntrl was assessed 24 h after stimulation. mRNA expression was determined by RT-qPCR and normalized to GAPDH. (C-D) THP-1-derived macrophages were incubated with UV-AC, Cov2-AC, or CoV2-Sup obtained from (C) Calu-3 or (D) Vero CCL81 infected cell cultures. CD206 expression on the cell surface 24h after stimulation was assessed by flow cytometry (gated on Total cells/Single cells/Live cells). Representative histograms and quantification of geometric mean fluorescence intensity (MFI) of CD206 are shown.

(E) Representative gating strategy used to determine the CD206 expression on the surface of CTFR-labeled macrophages stimulated with CTV-labeled AC in Figure 2C.

(F) Quantification by TCID₅₀ of the SARS-CoV-2 viral loads in 1.0 ×10⁶ AC isolated from Vero CCL81 infected cells, either intact (CoV2-AC) or UV-irradiated for 20 min for viral inactivation (CoV2-AC UV).

Boxes represent the mean of three biological replicates, and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA (B-D) *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups; ns: non-significant. Data shown are from one representative out of two experiments performed independently with similar results.

• Download figure
• Open in new tab

Figure. 2. Engulfment of SARS-CoV-2-infected dying cells impairs macrophage anti-inflammatory functions.

Macrophages were unstimulated (Cntrl), co-incubated with apoptotic Vero CCL81 cells (AC) isolated from UV-irradiated (UV-AC) or SARS-CoV-2-infected (CoV2-AC) cultures, stimulated with the supernatants of infected apoptotic cells (CoV2-AC Sup), or infected with CoV2 (CoV2), as indicated.

(A) Heatmap showing the expression of M2-marker genes (CCL18, MRC1, MMP9, PPARG, and CD163) in human monocyte-derived macrophages 24 h after stimulation, calculated as log2 fold change relative to Cntrl. mRNA expression was determined by RT-qPCR and normalized to GAPDH. Data represent the mean of biological triplicates.

(B) Flow cytometric analysis of CD206 expression on the cell surface of human monocyte-derived macrophages (gated on Total cells/Single cells/Live cells; gating strategy in Figure 2 – figure supplement 1A) 24 h after stimulation. Representative histograms and quantification of geometric mean fluorescence intensity (MFI) are shown. The black line on the histogram represents FMO control.

(C-E) Flow cytometric analysis of CD206 expression on the cell surface of THP1-derived macrophages 24 h after stimulation, showing quantification of CD206 MFI. (C) Macrophages (labeled with CTFR) were stimulated with UV-AC or CoV2-AC (labeled with CTV) and CD206 expression was measured on macrophages that engulfed (+) or not (-) AC. Gating strategy is depicted in Figure 2-figure supplement 1E. (D) Where indicated, isolated CoV2-AC were UV-irradiated for 20 min prior to co-incubation for viral inactivation. (E) The responses of macrophages to UV-AC and CoV2-AC were compared to infection with CoV2 in the presence of UV-AC (UV-AC + CoV2).

Boxes represent the mean of three biological replicates using cells from a single donor (A, B) or THP-1-derived cells (C-E), and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA. *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups. Data shown are from one representative out of at least two experiments performed independently with similar results.

Conversely, the uptake of CoV2-AC, but not UV-AC, significantly increased IL6 expression in monocyte- and THP-1-derived macrophages (Figure 3A and Figure 3 – figure supplement 1A, respectively). The uptake of CoV2-AC also induced robust secretion of inflammatory IL-6 and IL-1β, in both primary (Figure 3B) and THP-1-derived macrophages (Figure 3 – figure supplement 1B and C). We did not detect these cytokines in macrophages stimulated with the cell-free conditioned supernatants (Figure 3B and Figure 3 – figure supplement 1A-C), and therefore this effect was not due to immune mediators or viral particles released from infected epithelial cells. Importantly, assessment of IL-6 secretion in macrophages stimulated with sorted out populations of infected dying cells revealed that cytokine production occurred specifically in response to apoptotic cells, but not permeabilized dead cells (Figure 3C). In agreement, CoV2-AC obtained from epithelial cells infected in the presence of osmoprotectant glycine, which non-specifically prevents the leak of cytosolic content by lytic permeabilization (Evavold et al., 2017; Frank et al., 2000) still induced robust secretion of IL-6 by macrophages (Figure 3 – figure supplement 1D).

• Download figure
• Open in new tab

Figure 3 – figure supplement 1. Support data for Figure 3.

(A) THP-1-derived macrophages were incubated with UV-AC, Cov2-AC, or CoV2-Sup obtained from infected Vero CCL81 cells or infected with CoV2. Expression of IL6 in macrophages showed as fold change relative to Cntrl was assessed 24h after stimulation. mRNA expression was determined by RT-qPCR and normalized to GAPDH.

(B-C) THP-1-derived macrophages were incubated for 24 h with UV-AC, Cov2-C, or CoV2-Sup obtained from Calu-3 or Vero CCL81 infected cell cultures, as indicated. (B) IL-6 and (C) IL-1β in the culture supernatants were quantified by CBA.

(D) ELISA quantification of IL-6 in the culture supernatants of THP-1-derived macrophages incubated with UV-AC obtained from Vero CCL81 cells infected in the presence of glycine (Gly, 5 mM), where indicated, to inhibit release of cytosolic content due to permeabilization.

(E) Quantification by TCID₅₀ of the viral loads in 1.0 x 10⁶ AC isolated from Vero CCL81 cells similarly infected with either Coxsackie virus (Coxsackie-AC) or SARS-CoV-2 (CoV2-AC).

(F) ELISA quantification of IL-6 in the culture supernatants of THP-1-derived macrophages stimulated with Cov2-AC for 24 h. PtdSer on the cell surface was blocked to inhibit AC binding by incubation of CoV2-AC with Annexin V (Ann, 0.01 and 0.1 μg/mL) prior to addition to macrophage cultures.

Boxes represent the mean of three (A, B, C, F) or four (D, E) biological replicates, and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA (A-D, F) or Student’s test (E), *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups; ns: non-significant. Data shown are from one representative out of two experiments performed independently with similar results.

• Download figure
• Open in new tab

Figure. 3. Engulfment of SARS-CoV-2-infected dying cells causes inflammatory cytokine production

Macrophages were unstimulated (Cntrl), co-incubated with apoptotic Vero CCL81 cells (AC) isolated from UV-irradiated (UV-AC) or SARS-CoV-2-infected (CoV2-AC) cultures, stimulated with the supernatants of infected apoptotic cells (CoV2-AC Sup), or infected with CoV2, as indicated.

(A) IL6 expression in primary human monocyte-derived macrophages 24 h after stimulation showed as fold change relative to Cntrl. mRNA expression was determined by RT-qPCR and normalized to GAPDH.

(B) CBA quantification of IL-6 and IL-1β in the culture supernatants of the monocyte-derived macrophages 24 h after stimulation.

(C, D) ELISA quantification of IL-6 in the culture supernatants of THP1-derived macrophages 24 h after stimulation. (C) Macrophages were stimulated with early apoptotic AC (Annexin V⁺ Zombie^-) or late apoptotic AC (Annexin V⁺ Zombie⁺) sorted out by flow cytometry. (D) Where indicated, isolated CoV2-AC were UV-irradiated (UV) for 20 min or fixed with 2% paraformaldehyde (Fix) for 10 min prior to co-incubation.

(E-H) ELISA quantification of IL-6 in the culture supernatants of human primary monocyte-derived macrophages 24 h after stimulation. (E) IL-6 secretion in response to CoV2-AC was compared to apoptotic cells isolated from Vero CCL81 cells similarly infected with Coxsackievirus (Coxsackie-AC). (F) PtdSer on the surface of CoV2-AC was blocked by incubation with Annexin V (Ann, 0.1 μg/mL) prior to addition to macrophage cultures to inhibit AC binding. (G) CoV2-AC and macrophages were co-incubated in the presence of Cytochalasin D (Cyto D, 10 μM) to inhibit AC internalization. (H) Macrophages were treated with Remdesivir (Remd, 20 μM) following stimulation with CoV2-AC to inhibit viral RNA transcription.

Boxes represent the mean of three biological replicates using cells from a single donor (A, B, F, G), two donors (E), or THP-1-derived cells (C, D), and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA. *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups. Data shown are from one representative out of at least two experiments performed independently with similar results.

Further, the induction of inflammatory cytokines by SARS-CoV-2-loaded corpses required the presence of viable viral particles, as UV-treatment or paraformaldehyde-fixation of isolated CoV2-AC abrogated IL-6 secretion (Figure 3D). Notably, we did not observe robust production of IL-6 in response to dying epithelial cells infected with Coxsackievirus (Figure 3 – figure supplement 1E and Figure 3E). This result suggests that augmented pro-inflammatory cytokine production is not a universal response of macrophages to the uptake of cell corpses infected with positive single-strand viruses.

To confirm the requirement of recognition and binding of CoV2-AC to promote inflammatory skewing, we stimulated macrophages in the presence of annexin V to mask PtdSer exposed on the surface of the apoptotic cells. In the absence of PtdSer ligation, macrophages did not secrete IL-6 in response to CoV2-AC (Figure 3F and Figure 3 – figure supplement 1F). Stimulation of macrophages in the presence of actin depolymerization inducer Cytochalasin D, which allows corpse binding but not internalization, also attenuated IL-6 production (Figure 3G). Finally, we tested the effect of antiviral drugs targeting viral RNA transcription (RdRP-mediated RNA synthesis). We found that Remdesivir treatment of primary macrophages following stimulation with CoV2-AC partially reduced IL-6 secretion, suggesting that viral sgRNA expression is important to modulate the inflammatory response (Figure 3H).

Collectively, these findings support that sensing and engulfment of dying cells carrying viable SARS-CoV-2 switch macrophage anti-inflammatory, resolutive programming in response to efferocytosis toward an inflammatory phenotype. The exacerbated cytokine production observed in response to the efferocytosis of infected cell corpses by macrophages may contribute to the cytokine storm associated with COVID-19 hyperinflammatory syndrome.

Efferocytosis of SARS-CoV-2-infected dying cells suppresses continual clearance of apoptotic cells

Several receptors mediate efferocytosis through recognition of PtdSer on the surface of a dying cell, either by direct binding or through a bridging molecule (Boada-Romero et al., 2020; Penberthy and Ravichandran, 2016). We found that engulfment of CoV2-AC by primary macrophages reduced the transcription of such PtdSer receptors, including the scavenger receptors CD36 and SRA1, αVβ5 integrin (ITGB5), and T cell immunoglobulin mucin receptor 4 (TIM4) and MER proto-oncogene tyrosine kinase (MERTK) (Figure 4A). In THP-1-macrophages, we also observed reduced expression of SRA1, ITGB5, and TIM4 receptors in response to infection with SARS-CoV-2 (Figure 4 – figure supplement 1A). Phagocytes can ingest multiple corpses in subsequent rounds of efferocytosis (Miyanishi et al., 2007; Morioka et al., 2018; Park et al., 2011; Yurdagul et al., 2020). Previous in vivo work showed that macrophages must continually remove ACs to promote efficient repair of injury and prevent the accumulation of secondarily necrotic cells (Wang et al., 2017; Yurdagul et al., 2020). To determine if repression of efferocytic receptors affects additional uptake of dying cells, we treated macrophages with CoV2-AC and subsequently fed them with apoptotic human Jurkat cells (UV-Jurkat). We found that the engulfment of SARS-CoV-2-infected corpses suppressed the efferocytosis of other apoptotic cells (Figure 4B-E, Figure 4 – figure supplement 1B). As infection with SARS-CoV-2 reduced the expression of some efferocytic receptors in THP-1 macrophages (Figure 4 – figure supplement 1A), we also tested if it affected the uptake of apoptotic cells. Comparatively, infection with SARS-CoV-2 reduced AC clearance to a lower extent than prior uptake of CoV2-AC (Figure 4 – figure supplement 1C). Notably, treatment with Remdesivir did not improve the subsequent uptake by macrophages stimulated with CoV2-AC, arguing that viral RNA replication is not required for suppression of efferocytosis in response to CoV2-AC (Figure 4F). This result also suggests that different mechanisms lead to pro-inflammatory cytokine production and efferocytosis repression in response to the uptake of infected corpses.

• Download figure
• Open in new tab

Figure 4 – figure supplement 1. Support data for Figure 4.

(A) THP-1-derived macrophages were incubated with UV-AC or CoV2-AC obtained from infected Vero CCL81 cells or infected with CoV2. Expression of the efferocytic receptors CD36, SRA-I, ITGB5, TIMD4, and MERTK showed as fold change relative to Cntrl was assessed 24 h after stimulation. mRNA levels were determined by RT-qPCR and normalized to GAPDH.

(B) Representative gating strategy used to determine CTV-labeled macrophage populations for the two-step efferocytosis experiment showed in Figures 4D and E.

(C) THP-1-derived macrophages (CTV-labeled) were infected or not with SARS-CoV-2 (CoV2) for 18 h and subsequently were incubated with CTFR-labeled, UV-irradiated apoptotic Jurkat cells (UV-Jurkat) for 2 h. UV-Jurkat uptake is shown as the percentage of CTFR⁺ (gated on Total cells/Single cells/Live cells/ CTV⁺/ CTV⁺ CTFR⁺).

Boxes represent the mean of three biological replicates, and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA (A) or Student’s test (C). *, p<0.05, **, p<0.001, ***, p<0.0001 comparing the indicated groups; ns: non-significant. Data shown are from one representative out of at least two experiments performed independently with similar results.

• Download figure
• Open in new tab

Figure 4. Engulfment of SARS-CoV-2-infected dying cells suppresses continual efferocytosis by macrophages

(A) Expression of efferocytic receptors in human primary monocyte-derived macrophages. Macrophages were unstimulated (Cntrl), co-incubated for 24 h with apoptotic Vero CCL81 cells (AC) isolated from UV-irradiated (UV-AC) or SARS-CoV-2-infected (CoV2-AC) cultures, or infected with CoV2 (CoV2), as indicated. Expression of CD36, SRA1, ITGB5, TIMD4, and MERTK are showed as fold change relative to Cntrl. mRNA levels were determined by RT-qPCR and normalized to GAPDH.

(B-F) Flow cytometric analysis of two-step efferocytosis in vitro. THP-1-derived macrophages (CTV-labeled) were incubated with UV-AC or CoV2-AC isolated from Vero CCL81 cell cultures (labeled with pHRodo) for 18 h, and subsequently incubated with UV-irradiated apoptotic Jurkat cells (UV-Jurkat, CTFR-labeled) for 2 h. (B) Schematic representation of consecutive co-incubations of THP-1-derived macrophages with apoptotic cells prior to the assessment of cell corpse uptake by flow cytometry. (C) Representative maximal projection of scanning confocal images showing the uptake of AC (pHRodo) and UV-Jurkat (CTFR) by macrophages (CTV) after consecutive co-incubations. (D) Percentage of macrophages with internalized UV-AC or CoV2-AC (gated on Total cells/Single cells/Live cells/CTV⁺/CTV⁺ pHRodo⁺); gating strategy is displayed in Figure 4-figure supplement 1B). (E) Percentage of macrophages that engulfed UV-Jurkat within pHRodo⁺ populations (gated on Total cells/Single cells/Live cells/ CTV⁺/pHRodo⁺/ pHRodo⁺ CTFR⁺). (F) Percentage of macrophages that engulfed UV-Jurkat within pHRodo⁺ populations (gates as in (E)) wherein macrophages were treated with Remdesivir (Remd, 20 μM) following stimulation with CoV2-AC.

Boxes represent the mean of three (A) or four (D-F) biological replicates using cells from a single donor in (A) or THP-1-derived macrophages (D-F), and error bars are ± S.E.M. Each biological replicate is shown as a circle. Significance was calculated by ANOVA (A, F) or Student’s test (D, E). *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups. Data shown are from one representative out of three experiments performed independently with similar results.

• Download figure
• Open in new tab

Figure 5. Lung monocytes and macrophages of severe COVID-19 patients express reduced levels of efferocytic receptors.

(A-B) Representative histological findings in post-mortem lung tissue from control patients (Cntrl) and COVID-19 patients, obtained by ultrasound-guided minimally invasive autopsy. Tissue samples were immunolabeled with anti-S100A9 (phagocytes, green), (A) anti-CD36 (red), or (B) anti-MERTK (red) and stained with DAPI (nuclei, blue). Representative images show cropped details of lung tissues scanned by wide-field epifluorescence imaging. Scale bar: 10 µm. The mean fluorescence intensity (MFI) of CD36 and MERTK of at least 300 S1009⁺ cells are shown (n = 6 control and 4 COVID-19 patient samples). Dots represent an S1009⁺ cell (A and B), the crossing line represents the mean, and error bars are ± S.E.M. Significance was calculated by Mann-Whitney test (A and B), *, p<0.05, **, p<0.001, ***, p<0.0001, comparing the indicated groups.

(C) GSEA analysis for efferocytosis-related gene sets in early-infiltrating phagocytes (S1009⁺ CCL18^-) and in anti-inflammatory monocytes-derived macrophages (CD14⁺ S11009^- CCL18⁺) clusters from the bronchoalveolar lavage of mild (M) and severe (S) COVID-19 patients versus healthy individuals. The corrplot depicts the normalized enrichment score (NES), and p value for the gene sets indicated on the y axis.

(D) Enriched efferocytosis-related Gene Ontology (GO) terms in genes repressed in virus-positive versus virus-negative alveolar macrophages.

Thus, efferocytosis of SARS-CoV-2-infected dying cells affects the expression of efferocytic receptors and impairs the continual removal of apoptotic cells by macrophages.

Lung monocytes and macrophages of severe COVID-19 patients express reduced levels of efferocytic receptors

To gain insight into the contribution of dysfunctional efferocytosis in COVID-19 pathogenesis, we assessed the expression of efferocytic receptors by immunofluorescence in lung tissues obtained from autopsies of deceased COVID-19 patients. We found a reduction in the protein levels of CD36 in S1009⁺ infiltrating phagocytes in the lungs of COVID-19 patients compared to control tissues (Figure 5A). Phagocytes in the lungs of COVID-19 patients also expressed lower protein levels of MERTK (Figure 5B).

Using publicly available single-cell RNA sequencing (scRNA-seq) data from bronchoalveolar lavage (BAL) (Liao et al., 2020), we performed enrichment analysis using the genes differentially expressed (DEG) in mild and severe COVID-19 patients compared with healthy individuals. We targeted efferocytosis-related gene ontology (GO) annotated pathways and customized gene sets based on the literature (Boada-Romero et al., 2020; Penberthy and Ravichandran, 2016) (Supplementary File 1). Independent analysis of clusters previously assigned as macrophages (Liao et al., 2020) revealed significant repression of gene sets related to efferocytosis in severe patients (Figure 5C). While differences were evident between moderate and severe patients were observed in clusters expressing markers of early infiltrating phagocytes (S1009⁺ CCL18^-), they were more pronounced in mature, likely anti-inflammatory macrophages (S11009^- CCL18⁺). Notably, reassessment on scRNAseq by Ren et al. (2021) showed enrichment in pathways related to efferocytosis for genes repressed in alveolar macrophages which were positive for SARS-CoV-2 RNA compared with those which were negative to the virus (Figure 5D). This finding is consistent with the notion that uptake of infected cells reduces the capacity of macrophages to clear other dead cells. Altogether, these data support that macrophages in the lungs of severe COVID-19 patients may also fail in their efferocytic capacity.

Study approval

The procedures followed in the study were approved by the Research Ethics Committee of Hospital das Clínicas de Ribeirão Preto (CEP-FMRP/USP) and by the National Ethics Committee, Brazil (Comissão Nacional de Ética em Pesquisa (CONEP), protocols 30248420.9.0000.5440 and 39722020.9.0000.5440. Written informed consent was obtained from recruited donors.

Isolation of human peripheral blood mononuclear cells, monocyte seeding, and macrophage differentiation

Peripheral blood mononuclear cells (PBMC) from healthy donors were isolated from peripheral blood. The samples were centrifuged at 400 g for 10 min at room temperature (RT), the plasma was discarded, and the cellular portion was diluted in phosphate-buffered saline (PBS). PBMC were isolated by density–gradient centrifugation using Histopaque-1077 (Sigma-Aldrich). The collected PBMC fraction was treated with ACK lysis buffer to lyse erythrocytes and washed twice with PBS. For macrophages differentiation, monocytes were cultivated for 7 d in RPMI 1640 media supplemented with penicillin (10.000 U/mL, GIBCO), streptomycin (10.000 µg/mL, GIBCO), 1% glutamine (GIBCO), and 10% heat-inactivated human serum (Sigma Aldrich) at 37°C in a 5% CO₂ atmosphere.

Cell Lines

Calu-3 and Jurkat cells were originally from (ATCC^®) and purchased from Banco de Células do Rio de Janeiro (BCRJ). Vero CCL81 and THP-1 cells were from Dario Zamboni’s lab. No authentication method was used. Calu-3 and Vero CCL81 cells were maintained in DMEM with 4.5 g/L glucose (GIBCO) supplemented with penicillin (10,000 U/mL), streptomycin (10,000 μg/mL), and 10% heat-inactivated fetal bovine serum (FBS, GIBCO) (DMEMc) at 37°C in a 5% CO₂ atmosphere. Jurkat cells and THP-1 cells were maintained in RPMI1640 media supplemented with penicillin (10,000 U/mL), streptomycin (10,000 μg/mL), 1% glutamine, and 10% FBS (RPMIc) at 37°C in a 5% CO₂ atmosphere.

THP-1-derived macrophages were obtained by treatment with phorbol 12-myristate 13-acetate (PMA, 50 ng/mL, Sigma-Aldrich) for 24 h, followed by media replenishment with fresh RPMIc and other 24 h incubation period before stimulation.

Cell lines were routinely screened by PCR for mycoplasma contamination, and all cell lines tested negative.

Viral stock

The SARS-CoV-2 strain (Brazil/SPBR-02/2020) was isolated from the first Brazilian case of COVID-19. The viral stock was produced in a monolayer of Vero CCL81 cells expressing ACE2 (Rodrigues et al., 2020), maintained in serum-free DMEM at 37°C in a 5% CO₂ atmosphere. Coxsackievirus B5 (Coxsackie) was originally shared by Dr. Roger M. Loria (Virginia Commonwealth University, Richmond, Virginia, USA) (Gomes et al., 2009). The viral stock was produced in a monolayer of HeLa cells maintained in DMEM supplemented with 2% FBS at 37°C and 5% CO₂. Viral stocks were propagated under BSL3 conditions and stored at - 80°C. Viral loads were estimated by titration in Vero CCL81 cells seeded onto 96-well plates and standard limiting dilution to confirm the 50% tissue culture infectious dose.

Infection and UV-irradiation of epithelial cell lines

Cell monolayers were rinsed with PBS, replenished with a thin layer of serum-free media, and infection with a freshly thawed virus aliquot at a multiplicity of infection (MOI) of approximately 0.05 was carried out for 1 h at 37°C and 5% CO₂ for viral adsorption. The cells were then topped with fresh media and incubated for 48 h. Sterile apoptotic cells (UV-AC) were generated by exposure to UV-C radiation (Calu-3 cells at 500 mJ/cm² and Vero CCL81 cells at 350 mJ/cm²) in a UV crosslinker (Fisher Scientific), followed by 6h incubation in DMEMc at 37°C and 5% CO₂.

Cell death assays

Cells were collected 48 hours post-infection or 6 hours post-irradiation, spun at 300 g for 5 min, washed with PBS, and resuspended in the appropriate buffer before analysis.

For the analysis of caspase-8 and caspase-3 cleavage, 2.0 x 10⁶ Vero CCL81 and Calu-3 cells stimulated as described were lysed in RIPA buffer (10 mM Tris-HCl, pH 7.4, 1 mM EDTA, 150 mM NaCl, 1% Nonidet P-40, 1% deoxycholate, and 0.1% SDS) supplemented with protease inhibitor cocktail (Sigma). Precleared lysates were quantified by Bradford assay, boiled in Laemmli buffer, resolved by SDS-PAGE and transferred to PVDF membranes. Rabbit anti-cleaved caspase-3/caspase-3/cleaved Caspase-8/caspase-8 (all from Cell Signaling), rabbit anti-SARS-CoV-2 Spike (Abcam), HRP-conjugated goat anti-actin (Cell Signaling), and species-specific horseradish-conjugated secondary antibodies were used for antigen detection. For the analysis of cleaved caspase-3 by flow cytometry, infection of Calu-3 cells was carried out in the presence of Z-VAD-FMK (Selleckchem) at 20 μM, or mock treatment with vehicle (DMSO), added to cell cultures after virus adsorption. 10⁶ collected cells were washed with PBS, and fixation, permeabilization, and staining were carried out using BD Cytofix/Cytoperm™ Fixation/Permeabilization Kit (BD Biosciences), following the manufacturer’s instructions. Samples were incubated with rat serum (1:400) for 5 min at RT followed by immunolabelling of intracellular active caspase-3 was performed with anti-cleaved caspase-3 (Cell Signaling, 1:400) and secondary donkey anti-rabbit IgG conjugated to Alexa Fluor 488 (Thermo Fisher Scientific, 1:2,500). Sample acquisitions were performed in a FACSVerse^TM (BD Biosciences) flow cytometer and data processing was performed using FlowJo_V10 software (BD Biosciences).

Cell death profiling in infected or UV-irradiated Vero CCL81 cells was performed with a viability dye (Zombie NIR^TM Fixable Viability Dye, Biolegend, 1:400) and annexin V (FITC annexin V, BD Biosciences, 1:200) co-labeling. Labeling of 1.0 x 10⁶ cells with Zombie NIR was carried out for 10 min on ice, followed by annexin V labeling for 15 min on ice with cells resuspended annexin V binding buffer (0.01M HEPES, 0.14M NaCl, 2.5 mM CaCl2, pH 7.4). Sample acquisitions were performed in a in a BD FACSMelody^TM (BD Biosciences) cell sorter and data processing was performed using FlowJo_V10 software (BD Biosciences).

To determine LDH release, infection or UV-irradiation of 2.0 x 10⁵ Vero CCL18 cells were carried out in DMEMc without Phenol Red. LDH release was measured in the culture supernatants using the CytoTox 96 Non-Radioactive Cytotoxicity Assay (Promega), following the manufacturer’s instructions.

Isolation of apoptotic cells

Apoptotic cells were generated from SARS-CoV-2-infected (CoV2-AC) or UV-irradiated (UV-AC) Vero CCL81 and Calu-3 cells, or Coxsackievirus B5-infected Vero CCL81 cells (Coxsackie-AC). At the established time points, the supernatants were collected and temporarily stored for further processing. Then, a suspension containing CoV2-AC, Coxsackie-AC or UV-AC was obtained by carefully pipetting fresh media up and down the plate surface and collecting loosely attached cells. Following spun at 300 g for 5 min, pelleted AC were washed with PBS, resuspended in RPMIc, and counted prior to incubation with macrophages. Enrichment in apoptotic cells were confirmed by annexin V staining and flow cytometry.

To sort CoV2-AC populations, annexin V⁺ cells were pre-enriched using magnetic separation (Annexin V MicroBead Kit, Miltenyi Biotec) following the manufacturer’s instructions. Then cells were labeled with Zombie NIR and annexin V-FITC, as described above, and cell sorting were performed in a BD FACSMelody^TM (BD Biosciences) cell sorter. Isolated CoV2-AC were washed in PBS, resuspended in fresh RPMIc, and counted prior to incubation with macrophages

For virus inactivation, isolated CoV2-AC suspension was placed onto tissue culture hood surface with an open lid for 20 min with the UV-C lamp on. Virus inactivation was confirmed by TCID₅₀, as described below. Alternatively, isolated AC were fixed with 2% paraformaldehyde (PFA) in PBS for 10 min and washed twice with PBS before adding to macrophage cultures.

To prepare the cell-free conditioned supernatants of the infected AC (CoV2-AC Sup), collected supernatants were spun at 300 g for 5 min to remove cellular debris; the resulting supernatant was filtered using a 0.22 μm syringe filter and immediately used for macrophage stimulation. The volume of CoV2-AC Sup used to stimulate macrophages was normalized to the amount of collected AC.

TCID₅₀ quantification of viable viral particles

1.0 x 10⁶ isolated CoV2-AC and their equivalent purified conditioned supernatants, obtained as described above, were lysed by snap freeze-and-thaw. To estimate viral load in apoptotic cells isolated by annexin V labeling, all infected cells were collected, followed by magnetic separation with annexin V MicroBeads (Miltenyi Biotec) according to the manufacturer’s protocol. 1.0 x 10⁶ pre-sorted cells and annexin⁺ cells were lysed by freeze-and-thaw for quantification. Viral loads were estimated by titration in Vero CCL81 cells seeded onto 96-well plates and expressed as 50% tissue culture infectious dose (TCID₅₀). Quantification was performed with the Reed-Muench method and plotted as TCID₅₀ units per mL (Reed and Muench, 1938).

In vitro assays with THP-1- and PBMC-derived macrophages

1.0 x 10⁶ THP-1- and PBMC-derived macrophages seeded on 24-well tissue culture plates were rinsed with warm PBS before stimulation. Macrophage cultures were then topped with UV-AC, CoV2-AC, or Coxsackie-AC resuspended in the appropriate serum-containing macrophage media at a 1:1 ratio. Macrophage cultures were rinsed out 2 h after incubation to remove unbound AC and topped with fresh media at the indicated time points, cell supernatants were collected and macrophages were rinsed at least twice with warm PBS to remove cell debris before proceeding to further analysis. To block AC ligation and efferocytosis, CoV2-AC were pre-incubated with 0.1 or 0.01 µg/mL of annexin-V (BD Biosciences) and macrophage stimulation was carried out in RPMIc supplemented with Ca²⁺. To block AC internalization, macrophages were pre-treated for 30 min, and incubation was carried out in the presence of Cytochalasin D (10 µM, Sigma). Where indicated, macrophages were topped with fresh media containing Remdesivir (20 µM, BOC Sciences) after removal of unbound AC.

For stimulation with conditioned supernatants, CoV2-AC-Sup were added directly to macrophages cultures and topped with fresh media for volume adjustment. Infection with SARS-CoV-2 or Coxsackie virus was carried out with virus adsorption at MOI 0.1 in serum-free media for 1 h before topping with fresh serum-containing media.

For gene expression analysis by RT-qPCR, rinsed macrophage cultures were lysed in Trizol (Thermo Scientific) and immediately frozen at -80°C.

To estimate CD206 surface expression, rinsed macrophages were removed with a cell scraper, spun at 300 g for 5 min, washed with PBS, followed by incubation with rat serum (1:400) for 5 min at RT and cell surface staining for 15 minutes at 4°C. Cells were labeled with Zombie Violet™ Fixable Viability Kit (Biolegend) and anti-CD206-FITC antibody (clone 19.2, BD Biosciences). The cells were then fixed with 2% PFA for 10 min and washed before flow cytometric analyses in a FACSVerse (BD Biosciences) flow cytometer and data processing was performed using FlowJo_V10 software (BD Biosciences).

For immunofluorescence imaging of efferocytosis in vitro, human monocyte-derived macrophages seeded onto 13mm round coverslips were incubated with UV-AC or CoV2-AC labeled with 5 µM CellTrace™ CFSE dye (Thermo Scientific) for 2 h. After rinsing out non-engulfed AC, macrophages were fixed for 20 min with PBS 2% PFA at RT. Cells were then permeabilized with 0.1% Triton-X-100 in PBS for 10 min at room temperature RT. Next, cells were rinsed twice in PBS and incubated with PBS 1% bovine serum albumin (BSA) and 5 µg/mL normal donkey IgG (Jackson ImmunoResearch) for 45 min RT prior to overnight (ON) incubation at 4°C with primary antibodies: rabbit polyclonal anti-SARS-CoV-2 Spike antibody (Abcam, 1:1,000) and mouse mAb anti-rat CD11b (Abcam, 1:300) diluted in PBS containing 1% BSA. Next, cells were rinsed thoroughly in PBS and incubated for 30 min at RT with the secondary antibodies diluted in PBS: Alexa-594 goat-anti rabbit IgG (1:1,000; Thermo Scientific) and Alexa-647 donkey-anti mouse IgG (1:1,000; Thermo Scientific). Cell nuclei were stained with Hoechst 33342. Cells were then rinsed, and coverslips mounted on glass slides with Fluoromount-G (Invitrogen). Samples were imaged using a Zeiss LSM 780 laser scanning confocal microscope (Carl Zeiss), postprocessed for brightness and contrast adjustment with ImageJ image processing package (NIH) and layouts were built on Adobe Photoshop.

For the engulfment assays by flow cytometry, macrophages were labeled with 5 µM CellTrace™ Violet dye (CTV, Thermo Scientific), according to the manufacturer’s recommendations.

CellTrace™ Violet (CTV)-labeled macrophages were stimulated with UV-AC or CoV2-AC previously stained with 1 µM CellTrace™ Far Red dye (CTFR, Thermo Scientific), or 5 µM CellTrace™ CFSE dye, as described in the legend of the figures. Following the removal of unbound AC, macrophages were collected with a cell scraper, fixed with 2% PFA for 10 min, and washed before proceeding with the acquisition by flow cytometry. For the two-round efferocytosis assay, first-round UV-AC and CoV2-AC were labeled with a pH indicator fluorogenic intracellular probe (pHRodo Red AM Intracellular pH Indicator - Thermo Scientific). Second-round apoptotic Jurkat cells (UV-Jurkat) were generated by labeling with 1 µM CTFR and UV-C irradiation at 20 mJ/cm², followed by incubation at 37°C and 5% CO₂ for 6 h. CTV-labeled THP-1 cells were incubated with pHrodo-UV-AC or pHrodo-CoV2-AC for 18h at a 1:1 ratio. Cells were then rinsed with PBS, and CTV-labeled UV-Jurkat were added to the culture at a 1:1 ratio for 2h. Macrophage cultures were collected with a cell scraper, washed, fixed with 2% PFA for 10 min, and washed before proceeding with the acquisition by flow cytometry. Flow cytometric acquisition was performed in a FACSVerse (BD Biosciences) flow cytometer and data processing was performed using FlowJo_V10 software (BD Biosciences).

Cytokine quantification

IL-6 and IL-1β cytokine levels in macrophage supernatants were evaluated by ELISA assay (R&D Systems) or Cytometric Bead Array (CBA, BD Biosciences), as indicated in the legend of the figures, following the manufacturers’ recommendations.

RNA isolation and gene expression analyses

Total RNA extraction for human monocyte- or THP-1-derived macrophages stimulated in vitro was performed using Directzol^TM RNA miniPrep Kit (Zymo Research), following the manufacturer’s recommendations. Reverse transcriptase was performed using M-MLV reverse transcriptase (Thermo Scientific) according to the manufacturer’s recommendations.

Gene expression RT-qPCR was performed using SybrGreen Master Mix (Applied Biosystems). RT-qPCR was performed in fast mode, following the manufacturer’s recommendations. The evaluation of each gene expression was determined by the comparative CT method. Primer sequences used were:

Histological analysis

Ultrasound-guided minimally invasive autopsies for COVID-19 deceased patients were approved by the Research Ethics Committee of Hospital das Clínicas de Ribeirão Preto (CEP, protocol no. 4.089.567). Supplementary File 2 describes clinical, laboratory, and treatment records of COVID-19 patients. Non-neoplasic sections of lung parenchyma were obtained from lobectomies for lung cancer as a control group (Rodrigues et al., 2020). Paraffin sections from lungs were collected and processed as previously described (Rodrigues et al., 2020). Sections were rinsed in PBS, and the antigen retrieval was performed by incubation of samples with 0.1% trypsin (ThermoFisher Scientific) in PBS at 37°C for 15 min. Next, sections were rinsed and incubated with Image-iT FX Signal Enhancer (ThermoFisher Scientific) for 30 min at RT. Sections were then rinsed in PBS and incubated for 45 min at RT in PBS containing 0.5% BSA and 5µg/mL normal donkey IgG (Jackson ImmunoResearch).

For analysis of active caspase-3, lung sections were labeled with primary antibodies diluted in PBS and incubated overnight at 4°C: rabbit polyclonal anti-human cleaved caspase-3 (Asp175) (Cell Signaling; 1:1,000) and mouse anti-dsRNA (SCICONS; 1:1,000). Then, sections were rinsed thoroughly in PBS and incubated for 45 min at RT with the secondary antibodies: Alexa 594 goat anti-rabbit IgG (1:500; Thermo Fisher Scientific) and Alexa 488 goat anti-mouse IgG F(ab)’2 (1:500; Thermo Fisher Scientific) diluted in PBS. Cell nuclei were stained with Hoechst 33342. Cells were then rinsed, and coverslips mounted on glass slides with Fluoromount-G (Invitrogen). For in situ efferocytosis analysis, lung sections were labeled with primary antibodies mouse mAb anti-human CD68 (Abcam; 1:100), rabbit mAb anti-human cytokeratin 18 (GeneTex, 1:500) diluted in PBS and incubated overnight at 4°C. Then, sections were rinsed thoroughly in PBS and incubated for 45 min at RT with the secondary antibodies F(ab)’2 goat anti-mouse IgG conjugated to Alexa 488 (1:500; Thermo Scientific) and goat anti-rabbit IgG conjugated to Alexa 647 (1:500; Thermo Scientific). Next, sections were rinsed in PBS and incubated overnight at 4°C with rabbit polyclonal anti-SARS-CoV-2 Spike (Abcam; 1:1,000) previously conjugated to Alexa 568, according to the manufacturer’s instructions (Zenon Alexa Fluor 568 Rabbit IgG Labeling Kit, Thermo Scientific). Sections were then rinsed in PBS and incubated for 15 min at RT with Hoechst 33342 (Invitrogen), rinsed thoroughly, and mounted on glass slides with Fluoromount-G (Invitrogen). Samples from Cl-caspase-3 and in situ efferocytosis experiments were imaged using an ECLIPSE Ti2 Series microscope equipped with a DS-Qi2 camera (Nikon Instruments Inc.). Analysis was performed using NIS Elements (Nikon Instruments Inc.), wherein mean intensity fluorescence for Cl-caspase-3 was determined for randomly selected epithelial cells. Representative images were postprocessed for brightness and contrast adjustment with ImageJ image processing package (NIH) and layouts were built on Adobe Photoshop.

For CD36 or MERTK expression analysis, lung sections were labeled with the primary antibodies diluted in PBS and incubated overnight at 4°C: rabbit polyclonal anti-human CD36 (PA5-81996; 1:100 ThermoFisher Scientific) or rabbit mAb anti-human MERTK (Clone SR29-07; 1:150; ThermoFisher Scientific) and mouse mAb anti-human calprotectin (Clone MAC387; MA1-81381; 1:100; ThermoFisher Scientific). Then, sections were rinsed thoroughly in PBS and incubated for 45 min at RT with the secondary antibodies Alexa 488 donkey anti-rabbit IgG (1:500; Thermo Fisher Scientific) and Alexa 594 goat anti-mouse IgG (Thermo Scientific) diluted in PBS. Sections were then rinsed, and coverslips mounted with DAPI Fluoromount-G. Samples were imaged using an Olympus BX61 Fluorescence Motorized Slide Scanner

Microscope Pred VS120 (Olympus). Analysis was performed using ImageJ image processing package (NIH), wherein mean intensity fluorescence for CD36 and MERTK was determined for randomly selected S1009⁺ cells. Representative images were postprocessed for brightness and contrast adjustment with ImageJ image processing package (NIH) and layouts were built on Adobe Photoshop.

Re-analysis of scRNA-seq Data

We re-analyzed single-cell transcriptomic data from bronchoalveolar lavage fluid (BAL) from patients with varying severity of COVID-19 disease and their respective healthy controls (Liao et al., 2020). The dataset is publicly available at https://covid19-balf.cells.ucsc.edu/. Downloaded data was imported into R environment version v3.6.3. Differential expression analysis was conducted using FindMarkers function in Seurat using the Wilcoxon test to compare mild and severe COVID-19 patients with healthy individuals for each cluster previously identified by authors. Differentially expressed genes between mild/severe COVID-19 patients and controls for each cluster were identified considering genes expressed in at least 5% of cells and FDR < 0.05 and |avg_logFC| > 0.1. In addition, a differentially expressed gene list comparing infected (identified SARS-Cov-2 transcripts) against non-infected alveolar macrophages were obtained from supplemental materials from (Ren et al., 2021). Gene Set Enrichment Analysis (GSEA) was performed using ClusterProfiler and fgsea R packages (Korotkevich et al., 2021; Wu et al., 2021) for each differentially expressed gene or log2 fold-change lists from Liao et al. (2020) and Ren et al. (2021) dataset. A custom gene set was also utilized as a reference for GSEA, incorporating GO:0006911, GO:0043277 gene ontology terms, and additional efferocytic related genes from the literature (Penberthy and Ravichandran, 2016; Boada-Romero et al., 2020) (Table S1). Significantly enriched sets were identified considering terms with p-value < 0.05.

Quantification and statistical analyses

Please refer to the legend of the figures for description of samples, sample sizes (biological replicates) and experimental replicates. No statistical tests were used to estimate sample size. Data were plotted and analyzed with GraphPad Prism 8.4.2 software (GraphPad Prism Software Inc., San Diego, CA). The statistical tests used are listed in the legend of the figures. For the in vitro assays, we used Student’s t-test to compare two experimental groups or ANOVA and multiple comparison correction, to compare three or more experimental groups. For quantification of histological datasets, normality tests were performed, and samples with non-Gaussian distribution were analyzed applying the Mann–Whitney test.