Monocytes are phagocytic cells of the innate immune system implicated in immunopathogenesis of viral infection. Zheng et al., (2021) reported significant reduction of Mon frequency in PBMCs from COVID-19 patients [30], suggesting possible migration of Mon into the tissues. However, role of human Mon in immunopathogenesis of SARS-CoV-2 infection is poorly understood. In this study, we evaluated expression of myeloid and migration markers in both COVID-derived PBMCs and COVID-derived Mon. We observed that both PBMCs and Mon up-regulated expression of mRNAs of myeloid markers, including CD14, α-integrins (CD11b and CD11c), FCγ receptors (CD64A), and macrophage scavenger receptors (CD163), as a function of the severity of COVID-19 patients (Figs. 1A and C). Further, COVID-derived PBMCs and Mon up-regulated the expression of migration markers, including endothelial adhesion molecules, including L-selectin and Intercellular Adhesion Molecule 1 (ICAM1), and CC/CXC chemokine receptors (CCR1, CXCR1), as a function of the severity of COVID-19 patients (Figs. 1B and D). Together, the results show that up-regulation of myeloid and migration markers is consistent in both COVID-derived PBMCs and COVID-derived Mon, suggesting that monocytes could be the main subset into PBMCs that respond to SARS-CoV-2 infection in humans, suggesting a possible immunopathogenic role of Mon activation in SARS-CoV-2 infection.
Previous reports have shown that SARS-CoV-2 infection antagonizes IFN-dependent antiviral response [46][47] [48]. Additionally, Hadjadj et al., (2020) reported that PBMCs from COVID-19 patients did not express significantly high levels of any of the 3 types of IFN [46]. Here, we explore the mRNA expression of IFNs and IL27 subunits in both COVID-19 infection derived PBMCs and Mon. Transcriptomic analysis showed that neither COVID-derived PBMCs nor COVID-derived Monocytes express IFN-I, including IFNα1 and IFNβ1 (Figs. 2A-B and 3A-B, respectively). Neither IFN-II such as IFNγ (Figs. 2C and 3C, respectively), nor IFN-III, except for IFNλ1, which was expressed in COVID-derived PBMCs from critical COVID-19 patients (Fig. 2D), but not in COVID-derived Monocytes (Fig. 2D). Next, we investigated whether IL27 positively correlated with COVID-19 severity. The transcriptomic analysis shows that the expression of both IL27 subunits, IL27p28 and EBI3, increased with severity of COVID-19 in both PBMCs and Mon (Fig. 2E-F and 3E-F, respectively). Collectively, the data suggest that an increase in IL27 expression correlated with COVID-19 indicating that IL27 might be a marker of antiviral response to SARS-CoV-2 infection or, possible biomarker of COVID-19 progression.
3.3. Interleukin 27 signaling components are expressed in COVID-derived PBMCs and Monocytes
Considering that COVID-derived PBMCs and Monocytes show high levels of both IL27 subunits, we evaluated mRNA expression levels of IL27 signaling components in these cells. We observed that COVID-derived PBMCs and Monocytes showed significantly increased levels of receptor subunit Gp130, the receptor-associated kinase JAK2, STAT1 and STAT3, and the negative regulator SOCS3 (Fig. 4A and 5A, respectively), as a function of COVID-19 disease severity. As shown in the Fig. 4B, we observed a significant positive correlation between mRNA expression level of IL27p28 subunit and IL27 signaling components in COVID-derived PBMCs, including Gp130 (R = 0.2779, P = 0.0002), JAK2 (R = 0.4997, P < 0.0001), STAT1 (R = 0.3676, P < 0.0001), STAT3 (R = 0.3299, P < 0.0001), and SOCS3 (R = 0.4619, P < 0.0001).
Since we and others have reported that IL27 is essential for antiviral immunity [68][65], we proceeded to determine ISG mRNA levels in both COVID-derived PBMCs and monocytes. While the expression of Guanylate-binding Protein 1 (GBP1), Interferon-gamma Inducible Protein 16 (IFI16), Interferon Induced Protein with Tetratricopeptide Repeats 2 (IFIT2), Interferon Induced Transmembrane Protein 1 (IFITM1), and MX Dynamin Like GTPase 1 (MX1) was increased in COVID-19-derived PBMCs (Fig. 4C), the expression of APOBEC3A, GBP1, IFIT1, IFITM1, ISG15, MX1, and OAS1 was increased in COVID-19-derived monocytes (Fig. 5B). Similarly, we found that the expression pattern of STAT1-dependent cytokines and chemokines, including IL15, TNF-related Apoptosis-inducing Ligand (TRAIL), and B-cell Activating Factor (BAFF), CCL2, CCL7, CXCL10, and CXCL11 were increased with severity of COVID-19 both in PBMC and Mon, while CCL8 mRNA was significantly expressed in moderate but decreased in severe and critical COVID-19 patients; CXCL9 was significantly expressed in critical COVID-19 patients, but not in other groups (Fig. 4E). As shown in Fig. 4F, a moderate-high positive correlation between mRNA expression of STAT1 and their target genes, including ISGs such as GBP1 (R = 0.8809, P < 0.0001), IFITM1 (R = 0.7820, P < 0.0001), MX1 (R = 0.7156, P < 0.0001); cytokines such as BAFF (R = 0.6665, P < 0.0001) and TRAIL (R = 0.7639, P < 0.0001), and chemokines, including CCL2 (R = 0.4078, P < 0.0001), and CXCL10 (R = 0.4341, P < 0.0001) was observed in COVID-derived PBMCs. The data suggest that COVID-derived PBMCs and Monocytes activated IL27 signaling induce a robust STAT1-dependent pro-inflammatory and antiviral response dependent of ISGs expression.
3.4. Expression levels of TLRs, TLR adapter, as well as TLRs downstream signaling components, is associated with the severity of COVID-19.
TLRs activation drives transcriptional expression of both IL27p28 and EBI3 genes through the induction of IRF1 and NF-κB, respectively [17]. To determine whether any TLRs were positively correlated with COVID-19 severity, as IL27, we analyzed the dataset for TLR expression in patients with differing severities of COVID-19 (Figs. 6 and 7). We found that the expression of TLR2 was increased with severity of COVID-19 in both PBMCs and Mon (Figs. 6A and 7A, respectively), similar to that observed with IL27 subunits (Figs. 2E and F). We note that TLR1 was also increased with severity of CODIV-19 in Mon (Fig. 7A). As well, expression of TLR1, TLR4, TLR5, and TLR8 was significantly elevated in patients with severe and critical COVID-19 and the expression of TLR7 was increased in patients with moderate and critical (Fig. 6A). By contrast, expression of TLR4, TLR5, TLR6, and TLR7 was not altered in COVID-19-derived Mon (Fig. 7A). The expression of TLR3 did not show changes with the progression of COVID-19 in both PBMC and Mon (Figs. 6A and 7A). The data suggest a link of IL27 and certain TLRs with disease progression in patients with COVID-19.
MyD88 is important for proinflammatory cytokine production during SARS-CoV-1 infection [69]. To determine whether MyD88 or other TLR adapter TRIF (TIR-domain-containing adapter-inducing interferon-β), plays a role in SARS-CoV-2-induced inflammatory responses, we analyzed mRNA expression of these factors in PBMC and Mon from patients with differing severities of COVID-19. We found that the expression of MyD88 was increased with severity of COVID-19 in both PBMC and Mon (Figs. 6B and 7B). By contrast, the expression of TRIF was significantly increased only in PBMC from patients with severe and critical COVID-19 but not on Mon (Figs. 6B and 7B). Furthermore, as observed in Fig. 6B, mRNA levels of key mediators of TLR signaling, including IRAK2, and IRAK3 were increased in PBMC from patients with severe and critical COVID-19. Also, while IRAK2 was elevated in Mon from patients with moderate and severe COVID-19, IRAK3 was only increased in severe COVID-19 (Fig. 7B).
Although our transcriptomic analysis showed no expression of IFNs, we proceeded to evaluate the expression of IRFs, since they are induced by signaling via PRRs, including TLRs. We found that the expression pattern of IRF7 was increased with severity of COVID-19 in both PBMC and Mon (Figs. 6A and 7A). While the expression of IRF1 was significantly increased in patients with severe and critical COVID-19-derived PBMC (Fig. 6B), its expression in Monocytes was increased only in patients with severe COVID-19 (Fig. 7B), as was observed with IL27p28 and EBI3 (Figs. 2E and F, and 3E and 3F). IRF3 expression did not show any changes with the disease progression either in PBMC or in Mon of COVID-19 patients (Fig. 6A and 7B, respectively). The mRNA expression of various components of the NF-κB complex, including NF-κB1, NF-κB2, RELB, and IκBα, was significantly induced in PBMC from patients with severe and critical COVID-19 (Fig. 6C). In contrast, these components were increased only in Monocytes from patients with severe COVID-19 (Fig. 7C). Similarly, mRNA levels of NF-κB-target genes, including IL1β and IL10 were significantly increased independently of severity of COVID-19; CXCL1 and cyclooxygenase 2 (COX2) mRNA expression was significantly higher in PBMC from patients with severe and critical COVID-19 (Fig. 6D). TNFα mRNA expression level was significantly increased only in critical COVID-19-derived PBMC (Fig. 6D). In COVID-19-derived Mon, TNFα, IL1β, IL6, IL10, CXCL8/IL8, and COX2 were also significantly increased depending on the severity (Fig. 7D). Collectively, the results shown here suggest an association of TLRs and pathogenesis of COVID-19.
3.5. SARS-CoV-2 Spike protein modulates expression of TLR1/2-MyD88, downstream signaling components and induces pro-inflammatory and antiviral response IL27-dependent in human macrophages.
Both TLR2 and MyD88 expression have been associated with severity of COVID-19 disease since they induced production of pro-inflammatory cytokines by sensing structural SARS-CoV-2-PAMPs [30]. Further, activation of TLR1/2-MyD88 signaling induces a robust NF-κB-dependent pro-inflammatory response and a low but significant IL27-dependent antiviral response in both human and murine macrophages [17]. To determine whether TLR1/2-MyD88 induction by SARS-CoV-2 Spike protein (S-protein) plays a role in the establishment of pro-inflammatory and antiviral response IL27-dependent, in macrophages, we reanalyze a publicly available RNA-seq GSE173488 (GEO) [67]. As shown in Fig. 8A, transcriptomic analysis shows that human MDMs stimulated with SARS-CoV-2 S-protein significantly up-regulate mRNA expression of TLR1 and TLR2, as well as TLR adapters and mediators, including MyD88, TRIF, IRAK2, and IRAK3, as found in COVID-derived PBMCs and Mon. These results are consistent with significantly high expression levels of NF-κB-complex, including NF-κB1, NF-κB2, RELA, RELB, and IκBα (Fig. 8B), as well as genes targeted by NF-κB (TNFα, IL1β, IL6, IL10, CXCL1, CXCL8/IL8, and COX2; Fig. 8C)..Additionally, we found that IRF1 and IRF7 mRNA expression was significantly increased in MDMs stimulated with SARS-CoV-2 S-protein; in contrast, expression of IRF3 was downregulated (Fig. 8D). Interestingly, while MDMs stimulated with SARS-CoV-2 S-protein do not express the IFN-I (IFNα1 and IFNβ1), IFN-II (IFNγ), and IFN-III (IFNλ1) mRNAs, significantly increased levels of IL27p28 and EBI3 mRNA expression was observed (Fig. 8E). Also, we found that the expression of IL27 signaling components (Gp130, JAK1, STAT1, STAT3, and SOCS3) was significantly elevated in MDMs stimulated with SARS-CoV-2 S-protein (Fig. 8F). As expected, we found significantly increased expression of AVP-encoding ISG mRNAs, including GBP1, IFIT2, IFITM1, ISG15, MX1, OAS2, and Viperin, as well as STAT1-dependent cytokines such as IL7, IL15, and TRAIL in MDMs stimulated with SARS-CoV-2 S protein (Fig. 8G and H). As expected, we found significantly increased expression of AVP-encoding ISG mRNAs, including GBP1, IFIT2, IFITM1, ISG15, MX1, OAS2, and Viperin, as well as STAT1-dependent cytokines such as IL7, IL15, and TRAIL in MDMs that were stimulated with SARS-CoV-2 S protein (Fig. 8I). Overall, the data suggest that stimulation of human MDMs with SARS-CoV-2 S-protein induces TLR1/2-MyD88 signaling, which in turn promotes mRNA expression of both IL27 subunits, IL27p28 and EBI3, leading to a strong IL27-dependent pro-inflammatory and antiviral response, in IFN-independent manner.