ABSTRACT
Background Innate lymphoid cells (ILCs) are involved in immunity and homeostasis but, except for natural killer (NK) cells, their role in human viral infections is not well known. Puumala virus (PUUV) is a hantavirus that causes the acute zoonotic disease hemorrhagic fever with renal syndrome (HFRS). HFRS is characterized by strong systemic inflammation and NK cells are highly activated in HFRS, suggesting that also other ILCs might be responding to infection.
Methods Here we phenotypically analyzed peripheral ILCs in acute and convalescent PUUV-infected HFRS patients. Additionally, plasma levels of soluble factors and viral load were analyzed.
Findings Overall, the frequencies of NK cells and naïve ILCs were reduced while the frequency of ILC2, in particular the ILC2-lineage committed c-Kitlo ILC2 subset, was increased during acute HFRS. Interestingly, we observed a negative correlation between viral load and frequencies of both NK and non-NK ILCs in acute HFRS. Phenotypically, ILCs displayed an activated profile with increased proliferation, and showed altered expression of several homing markers during acute HFRS. In line with the observation of activated ILCs, plasma levels of inflammatory proteins, including the ILC-associated cytokines interleukin (IL)-13, IL-23, IL-25, IL-33, and thymic stromal lymphopoietin (TSLP), were elevated during acute HFRS.
Interpretation These findings indicate a general involvement of ILCs in response to human hantavirus infection. Further, this constitutes the first comprehensive study of ILCs in a hantavirus-caused disease, aiding in further understanding the role of these cells in disease pathogenesis and in human viral infections in general.
Funding A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.
Evidence before this study Innate lymphoid cells (ILCs) include a broad range of innate cell subsets involved, among other functions, in early response to infections. Natural killer (NK) cells have been broadly characterized in human viral infections, but much less is known in such context about other more recently discovered ILCs. Puumala virus is one of the causative agents of hemorrhagic fever with renal syndrome (HFRS), an acute zoonotic disease characterized by systemic inflammation. No current treatment or vaccines are available to date for hantavirus-caused diseases and pathogenesis is still not fully understood. A PubMed search up until April 2022 using a combination of the terms virus, Puumala virus, hantavirus, HFRS, ILCs, and NK cells shows that while studies have been performed characterizing NK cells in HFRS, no studies are available on the rest of ILCs in hantavirus-caused diseases, which in addition have been studied in only a handful of human viral infections.
Added value of this study In this study we thoroughly characterized the ILC landscape in circulation and its milieu in acute and convalescent Puumala-infected HFRS patients. We found that the frequency of NK cells and other ILC subsets was altered in the acute HFRS patients as compared to convalescent HFRS patients and control individuals, with a decrease in NK cells and naïve ILCs, and an increase in ILC2. In particular the ILC2-lineage committed c-Kitlo ILC2 subset was found increased. Total ILCs showed increased levels of activation and proliferation, and signs of altered migration patterns in acute HFRS patients. The finding of elevated levels of several soluble inflammatory proteins associated with ILCs in acute HFRS patients further suggests an implication of ILCs in HFRS. Moreover, an association between the frequency of NK cells and non-NK ILCs and plasma viral load suggest a possible connection between viremia and the activity of these cell types during the course of disease.
Implications of all the available evidence Our research shows that all circulating ILCs, including non-NK ILCs, are activated, proliferating, and correlate with viral load in acute HFRS patients, indicating a general involvement in human hantavirus infection and disease. Being the first comprehensive study on ILCs in hantavirus infections, further research will give relevant insight into the roles of ILCs in disease pathogenesis and protection.
INTRODUCTION
Innate lymphoid cells (ILCs) are a group of innate immune cells that play important roles in the modulation of immune and inflammatory responses.1 Naïve ILCs (nILC) constitute an immature subset 2,3 that can home from peripheral blood to tissues where they give rise to the mature ILC subsets 4. Mature ILCs are classified in five main subsets based on the transcription factors they express and the cytokines they produce: natural killer (NK) cells, ILC1, ILC2, ILC3, and lymphoid tissue inducer (LTi) cells.5,6 NK cells share the same features as ILC1, but in addition can, as opposed to the other ILCs, kill virus-infected cells.6 ILC1 and ILC3 are mainly found in mucosal tissues while NK cells and ILC2 are found both in tissues and in peripheral blood.2,4 Furthermore, two functionally distinct subsets of ILC2 can be found in peripheral blood: c-Kitlo and c-Kithi ILC2, the first being more committed to the ILC2 lineage.7
NK cells have been extensively described in different viral infections.8–10 On the other hand, non-NK ILCs (hereinafter referred to as ILCs) were more recently discovered,11,12 and thus less is known regarding their role in viral infections. Due to their location in mucosal tissue, ILCs are on the first line of defense and hence potentially essential in the early phases of viral infections.13 Their enrichment in lungs suggests an important role for ILCs in respiratory viral infections.13,14 There are conflicting data regarding ILC2 and virus infections: they have been suggested to promote tissue repair and protect the lungs of influenza-infected mice,15,16 while other studies reported that ILC2 induce airway hyperreactivity in influenza-infected 17–19 and respiratory syncytial virus-infected 20 mice. Furthermore, ILCs have recently been investigated in the context of a few human viral diseases. In HIV-1 infected individuals, levels of ILCs were found to be reduced in ileum and colon, 21 as well as in circulation, and to negatively correlate with viral load.22 Total peripheral ILC levels were also found to be decreased in SARS-CoV-2-infected coronavirus disease-19 (COVID-19) patients, with ILC2 decreased in severe but not in moderate patients.23–25 Moreover, in infants with respiratory syncytial virus bronchiolitis, elevated levels of ILC2 in the airways were found to associate with disease severity.26 Overall, these studies show an effect of viral infections on the ILC landscape.
Hantaviruses (genus Orthohantavirus, family Hantaviridae, order Bunyavirales) are RNA viruses that can cause zoonotic diseases in humans 27. The reservoirs of disease-causing hantaviruses are rodents, and humans are normally infected via inhalation of hantavirus-infected rodent excreta.28,29 Hantaviruses can cause hemorrhagic fever with renal syndrome (HFRS) in Eurasia and hantavirus pulmonary syndrome (HPS) in the Americas, with up to 10% and around 35% case fatality rate, respectively.28,30,31 Puumala virus (PUUV) is endemic in Europe, where it is the most common causative agent of HFRS.30,32 Hantavirus-infected patients initially develop non-specific symptoms such as high fever and headache, and many eventually present kidney dysfunction (primarily in HFRS), lung dysfunction (primarily in HPS), and gastrointestinal symptoms, including abdominal pain, diarrhea, vomiting, and gastrointestinal bleeding.31,33,34 As of today, no specific treatment nor United States Food and Drug Administration / European Medicines Agency-approved vaccines are available for hantavirus-caused diseases.35 Hantaviruses trigger immunopathogenic responses that likely contribute to hyperinflammation 36–43 and vascular leakage in patients, but the exact mechanisms leading to these events remain to be understood.32,44,45 Peripheral NK cells show signs of strong activation and proliferation in the acute phase of HFRS.46–49 B and T cells are highly expanded in the circulation of HFRS and HPS patients, with concomitant elevated levels of CD8+ T cells in the respiratory airways, 36,50–54. Mucosa-associated invariant T (MAIT) cells were recently shown to be reduced but highly activated in circulation during PUUV-caused HFRS.42 Levels of neutrophils are highly increased and also strongly activated in hantavirus-infected patients.55–57 Mononuclear phagocytes are susceptible to infection with hantaviruses, and their activation and redistribution from circulation towards the airways and kidneys in HFRS patients have been reported.58–61 Combined, these reports show that hantaviruses trigger strong immune cell responses, indicating a possible involvement of ILCs as well.
Here we performed a detailed characterization of peripheral blood ILCs and NK cells, as well as of their cytokine and chemokine milieu, in PUUV-infected HFRS patients. We observed increased plasma levels of inflammatory proteins, including ILC-associated cytokines. We showed that NK cell frequencies are reduced during acute HFRS but recover during convalescence. While total ILC frequencies did not change, we report increased frequency of ILC2 and a concomitant decreased frequency of nILC during acute HFRS. In particular, the ILC2-lineage committed c-Kitlo ILC2 subset was increased during acute HFRS. Furthermore, NK cells and ILCs displayed an activated phenotype and ongoing proliferation during acute HFRS. Interestingly, we observed a negative correlation between viral load and the frequencies of both NK cells and ILCs in acute HFRS, suggesting a potential direct or indirect influence of hantaviruses on the ILC landscape in HFRS patients.
MATERIAL AND METHODS
Patient samples
17 patients with serologically confirmed acute HFRS were included in the study. The patients were diagnosed with PUUV infection at the Tampere University Hospital, Finland during 2002–2007. Whole blood samples were collected, and peripheral blood mononuclear cells (PBMCs) were isolated as previously described and stored in liquid nitrogen at -150°C, while plasma was stored at -80°C, until further use. 61
The study was approved by the Ethics Committee of Tampere University Hospital (ethical permit nr. R04180) and all subjects gave written informed consent. The samples were collected in the acute phase of disease (5-8 days after onset of symptoms), during the early convalescence phase (20-27 days after onset of symptoms), and at a late convalescence phase at 180 or 360 days after onset of symptoms. Patients were stratified as having either mild or severe HFRS based on a scoring system adapted from the sequential organ failure assessment scoring system, where the maximum levels of creatinine (4 = > 440, 3 = 300-440, 2 = 171-299, 1 = 110-170, and 0 = < 110 µmol/l), minimum level of platelets (4 = < 20, 3 = 20-49, 2 = 50-99, 1 = 100-150, and 0 = > 150 x 103/µl), and minimum mean arterial blood pressure (1 = < 70 and 0 = ≥ 70 mmHg) were ranked. A total score of ≥ 5 was considered severe and < 5 mild. 61,62
As controls, PBMCs were obtained from buffy coats from 10 healthy blood donors from the Blood Transfusion Clinic at the Karolinska University Hospital Huddinge, Stockholm, Sweden (ethical permit nr. 2020-02604) and stored in liquid nitrogen until further use. PBMCs were isolated from the buffy coats by density centrifugation using Lymphoprep (StemCell Technologies), according to manufacturer’s guidelines.
Viral load in plasma
RNA was isolated from 140 µl of patient EDTA-plasma using column-based RNA isolation kit following the manufacturer’s instructions (Viral RNA mini kit, Qiagen). Isolated RNA was subjected to PUUV S RNA RT-qPCR analysis based on a previously described protocol,63 with TaqMan fast virus 1-step master mix (Thermo Scientific) using AriaMx instrumentation (Agilent).
Flow cytometry analysis
PBMC samples were thawed in RPMI (Cytiva) complete media [L-glutamine (ThermoFisher Scientific), FCS (Sigma-Aldrich), penicillin/streptomycin (Cytiva)] with DNase (Roche) and counted. 3-4 million cells per sample were stained. Briefly, cells were incubated with LIVE/DEAD Fixable Green Dead Cell Stain Kit (ThermoFischer Scientific) -used as a viability marker- and fluorochrome-conjugated antibodies directed against surface markers (Suppl Table 1) for 20 min at room temperature in the dark followed by 2 washes with flow cytometry buffer (2 mM EDTA in PBS). Cells were then fixed and permeabilized using FACS Lysing solution and FACS Permeabilizing solution (BD Biosciences), and next incubated with intracellular staining antibodies (Suppl Table 1) for 30 min at 4°C in the dark. Cells were then washed and resuspended in flow cytometry buffer. Samples were acquired with a BD LSR Fortessa (BD Biosciences) flow cytometer. Flow cytometric analysis was performed using FlowJo version 10.7.2 (TreeStar, Ashland). To ensure unbiased manual gating, a blinded analysis was implemented, whereby all FCS3.0 files were renamed and coded by one person and blindly analyzed by another person. All samples were compensated electronically, and gatings were based on fluorescent-minus-one (FMO) or negative controls. After all gatings were performed, samples were decoded and data analysis was performed.
ILCs were defined as live (DCM-) CD1a- CD14- CD19- CD34- CD123- BDCA2- FcεR- TCRα/β- TCRγδ- CD3- CD45+ CD127hi. NK cells were defined as DCM- CD1a- CD14- CD19- CD34- CD123- BDCA2- FcεR- TCRα/β- TCRγδ-CD3- CD45+CD56+ CD127lo/hi. For a detailed gating strategy see Suppl Fig 1.
Multiplex immunoassay
Plasma levels of IL-5, IL-6, IL-7, IL-10, IL-13, IL-15, IL-17A, IL-18, IL-23, IL-25, IL-33, IFN-γ, TNF, CCL20, CCL27, CCL28, GM-CSF, TSLP, and granzyme A were measured in plasma diluted 1:2, using a custom made Magnetic Luminex Screening assay (R&D Systems, Minneapolis) and analysed in a Magpix instrument (Luminex), according to manufacturer’s guidelines. For graphing purposes, values not calculated by the instrument’s software - because they were out of range in the lower standard range-were replaced by half of the value of the lowest detectable value for the given soluble factor measured.
Statistical analysis
Statistical analyses were performed using GraphPad Prism software v.9.2 for MacOSX (GraphPad Software). Statistical differences between healthy controls and HFRS patients were analyzed with Kruskal-Wallis test followed by Dunn’s multiple comparison post hoc test. Paired comparisons between HFRS in the different disease stages (acute, early convalescence, and late convalescence) were performed using Wilcoxon test. p-values < 0.05 were considered statistically significant. Spearman’s rank correlation coefficient was used for assessing correlations. Spearman’s correlation matrixes were generated with R (v.4.1.1; R Core Team, 2020) using package corrplot (v.0.9). Principal component analysis (PCA) was performed in R (v.4.1.1) using packages Factoextra (v.1.0.7), FactoMineR (v.2.4), RColorBrewer (v.1.1-2), and ggplot2 (v.3.3.5). Data was normalised in R using the scale argument within the PCA function. Where data was missing, the values were imputed using package missMDA (v.1.18).
Role of the funding sources
The funders of this study had no role in the study design, data collection, data analysis, data interpretation, or writing of the report.
RESULTS
Study design and patient characteristics
A total of 17 PUUV-infected hospitalized HFRS patients and 10 healthy controls were included in the study (Table 1). Peripheral blood samples were obtained at the acute (5-8 days after onset of symptoms), early convalescent (20-27 days after onset of symptoms), and late convalescent (180 or 360 days after onset of symptoms) phase of disease (Fig 1a). Hospitalized HFRS patients showed a typical clinical presentation during the acute phase, with thrombocytopenia, elevated C-reactive protein (CRP) and creatinine plasma levels, and viral load (Table 1, Fig 1b, and Suppl Table 2). For most of the patients all parameters normalized to levels within the normal range during the convalescent phase of disease (Fig 1b). The severity of the patients was assessed with a scoring system based on platelet counts, creatinine values, and mean arterial blood pressure values, as previously describe. 61,62 Two patients scored as severe, while all other scored as mild (Suppl Table 2).
HFRS patients present a strong inflammatory response during the acute phase of disease
Hantavirus disease is characterized by strong systemic inflammatory responses.37,40–43,64–67 Using a multiplex immunoassay, we assessed the plasma levels of 19 cytokines in HFRS patients during the acute and convalescent phase of disease. Except for interleukin (IL)-5, IL-7, and IL-17A (Suppl Fig 1a), we observed differences in levels of all proteins between acute and later stages of HFRS (Fig 2a and b). As previously reported,37,40–43,64–67 the plasma levels of tumor necrosis factor (TNF), IL-6, granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-10, interferon gamma (IFN-γ), IL-15, IL-18, and granzyme A (GrzA) were all significantly higher in acute HFRS as compared to the convalescent phases (Fig 2b). Further, we observed significantly higher levels of the type 2-associated cytokines IL-13, IL-25 (also called IL-17E), IL-33, and thymic stromal lymphopoietin (TSLP), related to ILC2 and T helper 2 cell activity and involved in functions such as tissue repair.68,69 We also observed significantly higher levels of the type 3-associated cytokine IL-23 in the acute phase of HFRS, involved in the activation of immune cells such as ILC3 and T helper 17 cells which have a role in protection from tissue damage (Fig 2b).5,70 Moreover, we observed that levels of the chemokines CCL20 and CCL27 were significantly higher in acute samples, while the level of CCL28 was significantly decreased (Fig 2b).
Principal component analysis (PCA) showed that samples from the acute phase separated from samples from the early and late convalescent phase. This separation was mainly driven by IL-10, GrzA, IFN-γ, TNF, IL-18, CCL27 and IL-33 (Fig 2c). Interestingly, patient 3, one of the two most severely ill patients in the cohort, deviated from the rest of the patients, both in the acute and early convalescent phase (Fig 2c). This patient showed the highest levels of several soluble factors and presented higher levels of many of them in the early convalescent phase than in the acute phase (Fig 2a), suggesting a longer than usual acute phase.
Additionally, we observed significant positive correlations between the type 2- associated cytokines IL-25 and IL-13 and between IL-25 and TSLP during the acute phase of HFRS (Suppl Fig 1b and Fig 2d and e). Further, out of the 15 acute HFRS patients, 13 were positive for PUUV S RNA in blood and a positive correlation was seen as well between viral load and IFN-γ (Suppl Fig 1b and Fig 2f).
Altogether, these results showed that PUUV-infected HFRS patients display a strong inflammatory response, including elevated levels of 16 cytokines many of which are known to be produced by or involved in the activation of ILCs and NK cells.
Peripheral NK cells are activated but decreased in frequency during acute HFRS
We next characterized the ILC and NK cell compartments in PBMCs from HFRS patients. For the identification and analysis of ILCs and NK cells, we used 18-parameter flow cytometry and a modification of a well-established gating strategy (Suppl Fig 2).71
We observed decreased frequencies of total CD56+ NK cells in peripheral blood in the acute and early convalescent phase of disease, which normalized in late convalescence (Fig 3a). There was a decreased frequency of CD56dim NK cells, with a concomitant increase of the smaller population of CD56bright NK cells, during the acute phase of HFRS (Fig 3b and c). High frequencies of CD69+ and HLA-DR+ total CD56+ NK cells were detected in the acute phase of HFRS, indicating NK cell activation (Fig 3d). As expected from the general NK cell activation, frequencies of NKp44+ and NKG2A+ NK cells were also increased in the acute phase of HFRS (Fig 3d). Furthermore, the frequency of Ki-67+ NK cells was increased, showing that the NK cells proliferated during acute HFRS (Fig 3d). Additionally, when analyzing for homing receptors, we observed a decreased frequency of α4β7+ NK cells and an increased frequency of CCR6+ and CCR10+ NK cells in the acute phase of HFRS (Fig 3d), suggesting an effect on NK cell migration. The frequency of CD45RA+ and CD161+ NK cells was decreased in HFRS patients as compared to healthy controls, with a tendency to recovery in the convalescent phase of HFRS (Fig 3d). Analysis of surface markers in the CD56bright and CD56dim NK cell subsets showed similar results as observed for total NK cells (Suppl Fig 3). A PCA revealed a separation of acute HFRS from convalescent HFRS patients and healthy controls based on the level of expression of the different analyzed surface markers in total NK cells (Fig 3e).
Next, we examined correlations between frequencies of NK cells and both soluble plasma proteins and clinical parameters (Suppl Fig 4a-d). IL-10 plasma levels positively correlated with the frequencies of both activated (CD69+) and proliferating (Ki-67+) NK cells (Fig 3f and g). Granzyme A levels correlated positively with the frequency of CD69+ NK cells (Fig 3h), and negatively with the frequency of CD161+ NK cells (Fig 3i). Interestingly, the frequency of activated (CD69+) NK cells correlated positively with viral load while the frequency of total NK cells showed a negative correlation with viral load during the acute phase of HFRS (Fig 3j and k).
Peripheral ILCs are activated and proliferate during acute HFRS
We next characterized the peripheral ILC responses in the HFRS patients. No significant difference in total ILC frequency was observed between the patients and the healthy controls (Fig 4a). Interestingly, as for NK cells, we observed a negative correlation between viral load and the frequency of ILCs in acute HFRS, showing reduced frequencies of peripheral ILCs in patients with higher viral loads (Fig 4b). Furthermore, we observed increased frequencies of activated (CD69+) and proliferating (Ki-67+) ILCs during the acute phase of HFRS (Suppl Fig 5a), while no differences were observed in the frequencies of NKp44, HLA-DR, and CD45RA-expressing ILCs (Suppl Fig 5a). When assessing expression of homing markers, we found a decreased frequency of α4β7+ ILCs in acute HFRS, but no significant difference in frequencies of ILCs expressing the chemokine receptors CCR6 and CCR10 (Suppl Fig 5a).
Next, we explored specific ILC subsets. CD117neg ILCs in peripheral blood have been shown to make up a heterogenous population with yet undefined functions.2 We therefore decided to focus our analysis on the more well-defined nILC and ILC2. The composition of these ILC subsets changed over time in the HFRS patients (Fig 4c). We observed an increase in ILC2 frequency and a decreased frequency of naïve ILC (nILC) in the acute phase of HFRS (Fig 4c and d).
Peripheral c-Kitlo ILC2 are increased in frequency during HFRS
Next, we characterized the phenotype of the ILC subsets. Increased frequencies of activated (CD69+) and proliferating (Ki-67+) nILC were observed in acute HFRS (Fig 4e). Moreover, similar to NK cells (Fig 3d), a decreased frequency of α4β7+ nILC was observed in acute HFRS (Fig 4e). No differences were observed in the frequencies of nILCs expressing HLA-DR, NKp44, CCR6, and CCR10 in HFRS patients as compared to healthy controls (Fig 4e). PCA based on the frequency of expression of surface markers in nILC showed a separation between the acute HFRS patients and healthy controls (Fig 4f).
The ILC2 population showed a similar phenotypic pattern as the nILC, with increased frequency of activated (CD69+) and proliferating (Ki-67+) cells during acute HFRS (Fig 4g). Additionally, a significantly decreased frequency of CCR6+ ILC2 was observed in the acute phase of disease (Fig 4g), while no differences were observed in the frequencies of ILC2 expressing NKp44, HLA-DR, CCR10, and α4β7 in HFRS as compared to healthy controls (Fig 4g). In line with these findings, a PCA based on the frequency of expression of surface markers on ILC2s showed no clear separation between the acute and convalescent HFRS, but a separation of the acute and control samples was observed (Fig 4h). Further, when analyzing for possible correlations to soluble proteins (Suppl Fig 6), we observed a positive correlation between the plasma levels of IL-10 and the frequency of CD69+ ILC2 (Fig 4i) and between plasma levels of the CCR10 ligand CCL27 and CCR10+ ILC2 (Fig 4j), as well as a negative correlation between plasma levels of TSLP and frequency of CCR10+ ILC2 in acute HFRS patients (Fig 4k). The first two correlations were also observed for total ILCs, as well as a positive correlation between plasma levels of IL-7 and Ki-67+ ILC (Suppl Fig 5b).
Having observed an increased frequency of ILC2 (Fig 4d), a decreased frequency of CCR6+ ILC2 (Fig 4g), and increased plasma levels of type 2 cytokines in acute HFRS (Fig 2), we next assessed whether there were changes in the ILC2 subsets in HFRS patients. Indeed, we observed increased frequencies of c-Kitlo ILC2 and, concomitantly, decreased frequencies of c-Kithi ILC2 during the acute phase of HFRS, as compared to convalescent HFRS patients and healthy controls (Fig 5a and b). Moreover, c-Kithi ILC2 showed higher frequency of CCR6 expression than c-Kitlo ILC2 both in patients and healthy controls (Fig 5c) and, aligning with the relative depletion of CCR6+ ILC2 in acute HFRS, (Fig 4g) we observed a lower frequency of CCR6+ c-Kithi ILC2 in the acute phase of HFRS as compared to the convalescent phase (Fig 5c).
DISCUSSION
Here we provide a detailed characterization of total ILCs, including both NK cells and non-NK ILCs, in blood samples from PUUV-infected HFRS patients. We reveal that total ILCs are activated and show an altered composition correlating to viral load.
Whereas NK cells have been extensively described in several human viral infections 8–10 including hantavirus infections, 46–49,72,73 ILCs remain understudied in human viral infections. Recent studies reveal dysregulated frequencies and phenotypes of ILCs in HIV-1,21,22 SARS-CoV2,23–25 rhinovirus,74 and respiratory syncytial virus (RSV) infection.26 Here we characterized ILCs and NK cell responses in the context of human hantavirus infection.
HFRS is characterized by a strong general immune activation, including hyperinflammation. 30,75 Some typical laboratory features are leukocytosis, thrombocytopenia, high CRP and hematocrit levels, and due to acute kidney injury, increased serum creatinine, proteinuria, and haematuria.28,31,76–78 Viral load peaks between the first 3 to 5 days during the febrile phase of disease, and viral load is normally not detectable in the convalescent phase of disease 31,32,79,80. In line with this, the HFRS patients in our study presented with typical laboratory findings together with a strong inflammatory response in the acute phase of disease.37,40–43,64–67 We observed an increase of several type 2-associated cytokines such as IL-13, IL-25, IL-33, and TSLP in acute PUUV-infected HFRS patients. The alarmins IL-25, IL-33, and TSLP are known activators of ILC2, which upon activation can secrete IL-13.15,81 These alarmins are upregulated upon infection or damage mainly of epithelial, endothelial, and stromal cells.82–84 Given that the main target cell of hantaviruses are endothelial cells,30 it is likely that hantavirus induce the secretion of the mentioned alarmins by these cells in HFRS patients. Interestingly, IL-33 has previously been found to be also elevated in plasma of Hantaan virus-infected HFRS patients and to positively correlate with disease severity.85
As recently shown by Resman Rus and colleagues,48 we observed transiently reduced NK cell frequencies in peripheral blood of PUUV-infected HFRS patients, returning to normal levels during convalescence. Moreover, as earlier shown,47 we observed that remaining circulatory NK cells were highly proliferating, with approximately half of them expressing Ki-67. During the acute phase of HFRS, NK cells also showed increased expression of several activation markers and altered expression of chemokine receptors associated with migration to tissues, such as lung (CCR6) and intestine (α4β7).86–89 This suggests migration of NK cells to tissues, which could be the cause of the observed decreased frequencies of NK cells in peripheral blood during acute HFRS. Interestingly, viral load correlated negatively with the frequency of NK cells in peripheral blood but positively with the frequency of activated NK cells, suggesting that active viral replication impacts, in an unknown manner, the activation of NK cells during HFRS. We have previously shown that hantaviruses have strong anti-apoptotic properties that potentially protect infected cells from cytotoxic lymphocyte mediated killing and may trigger NK cell mediated bystander killing of uninfected cells.72,73,90,91 Whether this holds true in vivo remains to be investigated, but the strong activation of NK cells we observe in HFRS patients could be a sign of continuous attempts of NK cells to kill hantavirus-infected cells. Of further interest, we observed that levels of IFN-γ also positively correlated with viral load in HFRS patients. IFN-γ is a cytokine with strong antiviral effects produced by a wide array of innate and adaptive lymphocytes, including ILCs, NK cells, NKT cells, and T cells.92–95 The main targets of hantavirus infection, such as endothelial cells and monocytes,30,59 are not known to express IFN-γ, so this correlation further suggests that active viral replication in infected cells indirectly leads to enhanced IFN-γ production by immune cells, such as activated NK cells.
We observed a decreased frequency of CD161+ NK cells in acute HFRS patients. The role of CD161 in NK cells is not fully understood, but it has been suggested to mark pro-inflammatory NK cells with a high ability to respond to innate cytokines.96 There are previous reports of modulation of CD161 expression in NK cells in other viral infections, such as a reduced CD161 expression in NK cells in acute hepatitis C virus infection which associated with enhanced viral clearance,97 reduced frequency of CD161+ NK cells in Chikungunya-infected patients,98 and depletion of CD161+ NK cells in cytomegalovirus infection.96 Paralleling the findings on NK cells, downregulation of CD161 expression on MAIT cells has been described in HFRS,42 as well as in HIV-infected individuals.99 Interestingly, we also observed a negative correlation between the plasma levels of granzyme A and the frequency of CD161+ NK cells in acute HFRS. Even though NK cells are not the sole source of granzyme A in plasma, this observation could suggest that the decrease in frequency of CD161+ NK cells - and thus increase of a CD161- NK population - is related to higher NK cell degranulation and cytotoxicity for viral control in HFRS patients. Along this line, we also observed a positive correlation between activated NK cells and the levels of IL-10, a cytokine that has been reported to enhance the effector functions of NK cells.100
The frequency of ILCs is reduced in the circulation of acute HIV-infected individuals, and correlates negatively with viral load.22 In COVID-19 we also recently reported a decrease in peripheral ILC frequencies and numbers.23 In contrast, here we did not observe a significant change in the frequency of ILCs in HFRS as compared to healthy controls. Interestingly though, and in line with the previous report for HIV infection,22 frequencies of peripheral ILCs in acute HFRS showed a strong negative correlation with viral load. Furthermore, as described for moderate COVID-19 patients,23 acute HFRS patients presented increased frequencies of ILC2 with a concomitant reduction of nILC frequencies. In contrast to our findings in COVID-19 patients,23 we observed increased frequencies of Ki-67-expressing ILC2 and nILC in HFRS patients. Moreover, the decreased levels of α4β7+ nILCs in acute HFRS suggest that the reduced frequency of peripheral blood nILCs could be due to their migration to tissues 4. Alternatively, nILC might differentiate into mature ILC subsets, such as ILC2s, in the circulation. Of interest, acute HFRS patients presented increased frequencies of c-Kitlo ILC2. Two ILC2 subsets have been defined, differing in their surface expression of c-Kit and their functionality. c-Kitlo ILC2 are more mature and ILC2-lineage committed, while c-Kithi ILC2 show plasticity towards an ILC3 phenotype and functionality.7 Moreover, c-Kitlo ILC2 express less CCR6 compared to c-Kithi ILC2.7 In line with this, acute HFRS patients showed decreased frequency of CCR6+ c-Kithi ILC2 as compared to convalescent patients and healthy controls. This suggests a skewing of c-Kithi ILC2 towards more ILC2-commited cells in HFRS, possibly explaining the increase in c-Kitlo ILC2 frequencies in acute HFRS. Alternatively, this decrease in CCR6+ c-Kithi ILC2 levels could also be due to migration of these cells to lungs, where ILC2 have been shown to play an important role in lung tissue repair during influenza infection in mice.15,16
Here we have characterized peripheral ILCs in PUUV-infected HFRS patients. Future characterization of ILCs in tissue samples, such as lung and intestines, can add important knowledge regarding circulating ILC infiltration and local ILC responses.
In conclusion, this study provides the first comprehensive characterization of total circulating ILCs in hantavirus-infected patients. We report an overall activated and proliferating ILC profile in these patients, with a particular increased frequency of the ILC2 subset, and a skewing towards the ILC2-lineage committed c-Kitlo ILC2 in acute HRFS. Additionally, we show that NK cells are reduced in frequencies and confirm that remaining circulating NK cells are highly activated and proliferating in acute HFRS. Moreover, we report a negative correlation between viral load and the frequencies of both NK cells and ILCs in acute HFRS, suggesting a potential influence of viral replication on these cells during the acute phase of hantavirus-caused disease.
Data Availability
The authors declare that the data supporting the findings of this study are available from the corresponding authors upon reasonable request.
CONTRIBUTORS
All authors read and approved the final version of the manuscript.
MG: data curation, funding acquisition, formal analysis, investigation, methodology, visualization, writing-original draft, writing-review & editing;
ACG: data curation, formal analysis, investigation, methodology, visualization, writing-review & editing;
JT: methodology;
KM: methodology;
AV: resources, funding acquisition;
SM: resources, funding acquisition;
JuM: resources, funding acquisition;
ASS: resources;
TS: resources; funding acquisition;
JeM: conceptualization, investigation, visualization, writing-review & editing, supervision;
JK: conceptualization, funding acquisition, investigation, visualization, writing-review & editing, supervision.
DECLARATION OF INTERESTS
The authors have no competing financial interests or conflicts.
DATA SHARING STATEMENT
The authors declare that the data supporting the findings of this study are available from the corresponding authors upon reasonable request.
ACKNOWLEDGMENTS
We thank the patients and volunteers who have contributed with clinical material to this study. We also thank Ms. Sanna Mäki for expert technical assistance. This study was supported by grants from the Swedish Research Council (projects K2015-56X-22774-01-3 and 2018-02646 to JK), Karolinska Institutet Research Foundation (project 2020-01469 to MG), the Academy of Finland (project 321809 to TS), Magnus Ehrnrooth Foundation (to AV), and Sigrid Jusélius Foundation (to AV and JuM), and the Competitive State Research Financing of the Responsibility Area of Tampere University Hospital (9AA050 and 9AB046 to JuM and 9AA052 to SM).
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.
- 39.
- 40.↵
- 41.
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.
- 52.
- 53.
- 54.↵
- 55.↵
- 56.
- 57.↵
- 58.↵
- 59.↵
- 60.
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.
- 66.
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.
- 84.↵
- 85.↵
- 86.↵
- 87.
- 88.
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.
- 94.
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵