HTNV infection promoted the expression of necroptosis-associated genes without triggering cell death.
Previous studies have demonstrated that programmed cell death (PCD), including apoptosis, necroptosis and pyroptosis, is a host defense against viral infection17. However, the relationship between cell death and HTNV, which causes a non-cytopathic effect (CPE), is still unclear. Interestingly, several death-associated pathways enriched according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis after RNA-sequencing both in HTNV and JEV infected cells (Fig. 1A). Heatmap showed that apoptosis and necroptosis-associated genes increased slightly post HTNV infection, while decreased after JEV infection (Fig. 1B). To confirm the expression level of cell death-associated genes post infection, we first collected RNA from primary macrophages after HTNV and JEV infection. qRT-PCR showed that RIPK3, a necroptosis-associate gene, presented an up-regulated trend and the expression of MLKL decreased after HTNV infection (Fig. 1C). We then investigated the dynamic expression of necroptosis-associated genes induced by HTNV infection. As shown in Fig. 1D, necroptosis-associated proteins (especially RIPK3) increased at the early stage of infection. Phosphorylated RIPK3 (Ser232) and phosphorylated MLKL (Ser345), the critical phosphorylation sites of necroptosis, did not change. Taken together, these results demonstrated that HTNV induced the expression of necroptosis-associated genes.
We then examined whether HTNV infection causes cell death. Obviously, ATP luminescence value showed no statistical difference after HTNV infection, and decreased at the early stage after JEV infection (Fig. 1E). The results of annexin V-PI staining showed that HTNV did not induce significant cell death in BMDMs. In contrast, the positive rate of annexin V-PI staining increased after JEV stimulation (Fig. 1F), which is considered to trigger cell death via activated RIPK315,18. We counted PI-positive cells under different conditions, and HTNV infection induced approximately 10% PI-positive cells with no significant difference compared to the uninfected group (Supplementary Fig. 1A). These observations indicated that HTNV infection could not trigger cell death and increased necroptosis-associated proteins might play noncanonical roles post-infection.
HTNV infection does not trigger cell death, and it is still unclear whether it can regulate death signaling. Flow cytometric assays showed that tumor necrosis factor α (TNF-α), smac mimetic and zVAD-fmk (TSZ) efficiently induced cell death in L929 cells compared with A549 cells (Supplementary Fig. 1B)19. The expression of MLKL and p-MLKL increased induced by TSZ (Fig. 1G). In addition, the ATP luminescence assay showed that HTNV infection could not repress the reduction in ATP content after TSZ stimulation for 6 h or 12 h (Fig. 1H). We then observed cell morphology under a microscope at 6 h after TSZ stimulation. HTNV could not prevent the cells from shrinking and becoming round (Fig. 1I). These results demonstrate that HTNV neither induced necroptosis nor antagonized cell death induced by TSZ.
RIPK3 regulated HTNV replication through JAK-STAT pathways.
To determine the functions of necroptosis-associated genes in virus infection, we first transfected primary macrophages with mouse Z-DNA Binding Protein 1 (ZBP1)/RIPK3/MLKL siRNA and we observed that silencing RIPK3 expression restricted nucleocapsid protein (NP) expression (Fig. 2A). Our previous results indicated that RIPK3 served as a regulatory protein for type I IFN signaling to promote JEV replication. We hypothesized that RIPK3 might assist in immune evasion post-HTNV infection. Consistent with our hypothesis, the level of HTNV decreased both at the RNA and protein levels in the RIPK3−/− BMDM (Fig. 2B). After transfection by RIPK3-Flag (pretreated with necrosulfonamide, an inhibitor of MLKL), the expression of NP was significantly increased (Fig. 2C).
To definitively confirm the function of RIPK3 in HTNV-infected cells, we performed RNA sequencing using infected BMDMs (wild type (WT) and RIPK3−/−, setting three replications) (Fig. 2D(i)). Based on the principal component analysis (PCA), the results showed the good fitness and predictive power of the model (Fig. 2D(ii)). In the volcano plot, differentially expressed genes (DEGs) were represented, of which, 133 were significantly downregulated (blue plots) and 636 were upregulated (orange plots), especially Mx1, Isg15 and Stat1(Fig. 2D(iii)). Furthermore, the IFN-associated PRRs signaling pathway and JAK-STAT pathway were enriched in KEGG analysis (Fig. 2E(i)). We found that the levels of STAT1 and STAT2 increased in the RIPK3−/− BMDMs. Knocking out RIPK3 significantly enhanced the expression of multiple inflammatory cytokines and ISGs (Fig. 2E(ii)), especially MX1, MX2, ISG15, interferon-induced protein with tetratricopeptide repeats 1 (IFIT1, also known as IFN-stimulated gene 56/ISG56), IFIT2 and IFIT3, which have been reported to be effective anti-hantaviral effectors9.
IFN-associated production lags behind HTNV replication at the early infection stage, and IFN responses activate at the late stage13. Primarily, we examined the dynamic expression of NP and IFN-associated signaling pathways post-HTNV infection. After infection with HTNV (MOI = 1), the viral level (RNA and protein) in BMDMs peaked at 48 hpi, together with the expression of PRRs (TLR3/TLR4) and IFN-α, as well as the downstream ISGs. However, the expression of p-JAK2 (Y1007, Y1008) and p-STAT1 (Ser727) increased at 12 hpi but gradually decreased from 24 to 36 hpi (Supplementary Fig. 2A, B).
To clarify the mechanism of RIPK3 in targeting IFN signaling, we used several reporter plasmids containing Luciferase activity. Interestingly, overexpression of RIPK3 in HEK293T cells downregulated STAT1-Luc and IRES-Luc expression but not IFN-Luc (Fig. 2F). Consistently, we detected IFN-α/β secretion by ELISA, and there were no differences between the WT and RIPK3−/− BMDMs after HTNV infection (Supplementary Fig. 3A). We then identified that HTNV-induced p-JAK2 (Y1007, Y1008) and p-STAT1 (Ser727) levels substantially increased in the RIPK3−/− BMDMs (Fig. 2G). As expected, the expression of STAT1 and p-STAT1 (Ser727) decreased after overexpressing RIPK3 (Fig. 2H). Of note, ISGs such as MX1, IFIT1 and ISG15 increased at 24 hpi in the RIPK3−/− BMDMs compared with those of the WT BMDMs (Fig. 2I), and ISGs were inhibited once RIPK3 was overexpressed (Fig. 2J). To determine whether RIPK3 regulates inflammation, we observed that the level of IL-6 increased while the expression of IL-10, an anti-inflammatory cytokine, reduced in the RIPK3−/− BMDMs, indicating that macrophages tended to exhibit pro-inflammatory M1 polarization (Supplementary Fig. 3B).
RIPK3 hindered anti-hantaviral ISG expression by interacting with STAT1.
Though RIPK3 could promote HTNV replication targeting JAK-STAT pathways, ambiguity remains regarding the mechanism by which they conduct it. We hypothesized that increased RIPK3 protein might interact with STAT1 and further inhibit its phosphorylation. As expected, we found that RIPK3 colocalized with STAT1 post HTNV infection (Fig. 3A). Additionally, the proportion of STAT1 which has translocated into the nucleus decreased after RIPK3 was over-expressed, indicating that RIPK3 acted as an inhibitor of the JAK-STAT pathway (Fig. 3B). Therefore, we carried out co-IP assays to determine whether RIPK3 interacts with STAT1. As expected, RIPK3-Flag physically interacted with exogenous STAT1-Myc (Fig. 3C, D). To investigate the interaction sites between RIPK3 and STAT1, we chalked out a rough sketch of the structural domains of RIPK3 and STAT1 (Fig. 3E). Mouse RIPK3 has 486aa with a kinase domain from 22 to 292aa and an RHIM motif from 440 to 461aa. Mouse STAT1 contains 750aa with an Src homolog 2 (SH2) domain from 573 to 670aa. Amino acid residues in the RIPK3 (Lys432, Thr338, Thr248 and Ser326, etc.) dynamically associate through hydrogen bond with residues of STAT1 (Thr719, Gln621, Asn530 and Arg619 etc.) (Fig. 3F). Based on the predicted interaction sites, we speculated that the interaction of RIPK3 and STAT1 might cap the phosphorylation site of STAT1, which suppress its further activation. And this noncanonical role of RIPK3 might be independent of the kinase activity and RHIM motif of RIPK3.
RIPK3-deficient mice showed clinical symptoms.
To assess the immune response of the RIPK3−/− mice to HTNV infection, we challenged RIPK3−/− mice with HTNV through intravenous injections on the first day and administered a booster on day six (Fig. 4A). All animals were monitored daily for clinical symptoms, body weight changes and temperature. Mice were sacrificed to analyze various laboratory parameters, viral loads, and pathology. Compared to the WT mice, RIPK3−/− mice showed clinical symptoms, such as ruffled fur and sluggish movement, as shown in the picture (Fig. 4B). The murine sepsis score (MSS) is used to assess the severity of sepsis in mice based on observational characteristics20. The results showed that murine sepsis symptoms alleviated in the RIPK3−/− mice (Fig. 4C). Additionally, the RIPK3−/− mice showed a slowing of weight gain (Fig. 4D). Notably, the core temperature of the RIPK3−/− mice indicated fever in the head, back and tail (Fig. 4E, supplementary Fig. 4A). Consistently, the body temperature of the RIPK3−/− mice were generally higher, while no fever was recorded in the WT mice, especially from 7 days to 9 days after the booster (Fig. 4F). We then performed a hematology test using whole blood collected from the RIPK3−/− and WT mice. The number of monocytes decreased in the RIPK3−/− mice, and neutrophils and lymphocytes increased post-infection, which indicated that there exited inflammation response in RIPK3−/− mice. The results showed that white blood cells (WBC) counts were higher in the RIPK3−/− group than in the WT mice. As thrombocytopenia is a common clinical feature in viral hemorrhagic fever, notably, blood platelet (PLT) counts in the infected WT mice decreased (≈ 5X10^10/L) but recovered in the RIPK3−/− mice (≈ 1.5X10^12/L) (Fig. 4G). We further collected serum to examine liver and kidney disorders. Compared with the WT mice, the serum levels of ALT, AST, and serum UA showed no elevation. These findings indicated that HTNV infection did not cause significant liver and kidney damage (Fig. 4H, supplementary Fig. 4B). All results mentioned above suggested that intravenous injections of HTNV into the RIPK3−/− mice, to some degree, alleviated the main pathological features of HTNV, which was possibly due to the activation of a robust immune response.
RIPK3-deficient mice showed HTNV clearance and an activated inflammatory response.
To evaluate the viral load in mice, we extracted HTNV RNA from the mouse serum to determine viral replication in the blood (Fig. 5A). The results showed that the viral copy number in the RIPK3−/− mice was lower than that in the WT mice (RIPK3−/− mice: 683.7 ± 305.5 copies per µL; WT mice: 1613 ± 273.1 copies per µL). Consistently, the expression of viral antigen in the major organs (including the lung, spleen and kidney) was significantly suppressed in the RIPK3−/− mice at protein levels (Fig. 5B).
Pathological changes in the HTNV-infected mice were evaluated by hematoxylin and eosin (H&E) staining and found to exist in the lungs of all infected RIPK3−/− mice, manifested as granulocyte and macrophage infiltration in the alveolus. The clinical signs caused by HTNV mainly include kidney injury characterized by acute tubulointerstitial nephritis involving infiltration of inflammatory cells21. However, pathological manifestations of the liver and kidney were not apparent in the HTNV-challenged RIPK3−/− mice (Fig. 5C). To further clarify macrophage infiltration into the lungs after HTNV infection, we performed IH staining with F4/80 (macrophages) and NP. The results showed that a group of immune cells positive for F4/80 infiltrated the alveolus after HTNV infection in the RIPK3−/− mice (Fig. 5D). We speculated that the infiltration of various inflammatory cells might cause lung damage post-infection in the RIPK3−/− mice. The red arrows in Fig. 5E showed the signs of emptying some lamellar bodies in alveolar type II cells. Importantly, increased macrophage infiltration and slight pulmonary fibrosis confirmed that RIPK3 could participate in the inflammatory response in vivo by regulating granulocyte and macrophage enrichment.
RIPK3-deficient mice increased the expression of ISGs by promoting the phosphorylation of STAT1.
To clarify the mechanism by which HTNV is cleared in the virus-challenged RIPK3−/− mice, we detected the expression level of STAT1 in organs. Consistent with the results in vitro, total STAT1 and p-STAT1 (Ser727) increased both in the lung and spleen with statistical analysis (Fig. 6A). We further found that the proportion of STAT1-positive cells was significantly higher in the lungs of the RIPK3−/− mice than in the lungs of the WT mice by IHE staining. Moreover, the proportion of STAT1 in the nucleus significantly increased in both NP-positive and NP-negative cells (Fig. 6B). We previously demonstrated that RIPK3 deficiency promoted HTNV clearance in BMDMs by inducing STAT1 phosphorylation and ISGs expression. However, whether these effects are macrophage-specific or universal in organs remains incompletely understood. To clarify this point, we collected RNA from the lung and spleen for qRT-PCR analysis. The results showed that HTNV infection significantly enhanced MX1, MX2, IFIT1 and ISG15 expression in the RIPK3−/− mice but did not affect TLR3 or TLR4, the PRRs that can be triggered by HTNV in previous reports (Fig. 6C, supplementary Fig. 4C)22. These above results indicated that RIPK3-deficient mice showed distinctive roles in clearing HTNV and promoting inflammation by regulating STAT1 activation.