ABSTRACT
During respiratory viral infections, the precise roles of monocytes and dendritic cells (DCs) in the nasopharynx in limiting infection and influencing disease severity are incompletely described. We studied circulating and nasopharyngeal monocytes and DCs in healthy individuals and in patients with mild respiratory infections (primarily influenza A virus, IAV). As compared to healthy controls (HCs), patients with acute IAV infection displayed reduced DC but increased intermediate monocytes frequencies in blood, and an accumulation of most monocyte and DC subsets in the nasopharynx. IAV patients had more mature monocytes and DCs in the nasopharynx, and higher levels of TNFα, IL-6 and IFNα in plasma and the nasopharynx. In blood, monocytes, the most frequent cellular source of TNFα during IAV infection, remained responsive to additional stimulation with TLR7/8L. Immune responses in older patients skewed towards increased monocytes rather than DCs suggesting a contributory role for monocytes in disease severity. In patients with other respiratory virus infections, we observed changes in monocyte and DC frequencies in the nasopharynx distinct from IAV patients, while differences in blood were more similar across patient groups. Together, our findings demonstrate tissue-specific and pathogen-specific patterns of monocyte and DC function during human respiratory viral infections and highlight the importance of comparative investigations in blood and the nasopharynx.
INTRODUCTION
Respiratory viral infections cause significant global disease burden with influenza, or flu, being responsible for a significant portion. An estimated 1 billion cases of influenza occur annually resulting in approximately 3-5 million severe cases and 290,000-650,000 deaths (1). In addition to seasonal epidemics caused by both influenza A or B virus (IAV and IBV, respectively), IAV can also cause pandemics. The majority of infections remain asymptomatic or develop mild to moderate respiratory disease, characterized by fever, nasal congestion, cough, and muscle aches. Severe disease mainly affects infants, pregnant women, and the elderly or immunocompromised, but can also occur in otherwise healthy individuals (2, 3). The determinants of disease severity are still incompletely understood but may include properties of the virus, environmental factors, genetics, and immune responses of the patient (4, 5). Other agents causing an influenza-like illness during influenza season include respiratory syncytial virus (RSV) and seasonal coronaviruses (OC43, HKU1, 229E, and NL63). Similar to IAV, coronaviruses are capable of causing pandemics, most notably the ongoing coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV- 2) (6).
IAV is primarily transmitted via inhalation of virus-containing aerosols or droplets, and mainly targets respiratory epithelial cells (7, 8), with the nasopharynx being the initial site of virus replication. IAV generally remains localized to the airways, despite signs of systemic inflammation (9). At the site of infection, resident innate immune cells including monocytes and dendritic cells (DCs) rapidly respond to the presence of virus by secreting cytokines, interferons, and chemokines to limit viral spread and recruit immune cells (10, 11). Monocytes and DCs shape the specificity and strength of the subsequent adaptive responses (12, 13). In blood, three subsets of monocytes are found: the CD14+CD16– classical monocytes (CMs), the most frequent subset at steady state, and the further differentiated CD14+CD16+ intermediate monocytes (IMs) and CD14–CD16+ nonclassical monocytes (NCMs) (14–16). Blood IMs expand rapidly in response to inflammation, infection, or vaccination (3, 12, 17). In addition, monocytes extravasate to tissue where they play an important role in innate immune protection. Monocytes also secrete TNFα, a major regulator of innate immune function that is central to the cytokine storm associated with IAV infection (18). Of the DCs, the CD1c+ myeloid DCs (MDCs) excel at activating naïve T cells (19); the CD141+ MDCs can cross-present antigens via MHC-I (20); and the CD123+ plasmacytoid DCs (PDCs) mediate type I IFN responses (21). Monocytes and DCs vary in distribution and function, depending on the anatomical compartment (22, 23). Moreover, monocytes and DCs are susceptible to IAV infection in vitro, and the cytopathic nature of the virus may impair their antigen processing and presenting functions (24–26), delaying recovery and normalization of immune cell distribution and function (27, 28).
Studies have shown that monocytes and DCs are recruited to the nasopharynx following infection with 2009 H1N1pdm IAV strains (3), and in individuals hospitalized with severe influenza infections (29, 30). Disease severity in hospitalized patients has been shown to correlate with (i) monocyte recruitment and increased levels of MCP3, IFNα-2 and IL-10 in the nasal compartment (3) and (ii) strong TNF-producing monocytic responses in blood (2) and inflammatory, neutrophil-dominant patterns (31). Despite accounting for a comparatively greater burden of disease, immune responses during mild seasonal influenza infections remain less studied. Therefore, the roles played by monocytes and DCs in contributing to or mitigating mild influenza disease are largely unknown. Additionally, while the response of blood monocytes and DCs to IAV has been studied well in vitro and in animal models (24, 25, 27, 32- 35), few studies compare responses between blood and the nasopharynx in human infections (31, 32, 36). Immune cell behaviour in the nasopharynx during IBV and RSV infections has not been studied in great detail but evidence of DC mobilization to the nasal cavity has been reported (29). Studies on immune responses to mild SARS-CoV-2 infection have also primarily focused on blood and rarely the upper airways (37).
Here, we determined monocyte and DC subset distribution, maturation and function in both blood and, for the first time, the nasopharynx, in patients with mild to moderate seasonal influenza and influenza- like infections. The methods described in this study allowed us to investigate airway immunity in a larger cohort of patients with SARS-CoV-2 infection (38, 39), showing that methods to study the immune responses in the nasopharynx during acute disease are essential tools as we face the possibility of future pandemics. Comparing the dynamics of systemic and nasopharyngeal immune function will add to our understanding of the roles of monocytes and DCs in shaping the nature and magnitude of inflammation and subsequently disease severity during respiratory viral infections.
RESULTS
Study subject characteristics
During three consecutive influenza seasons (2016–2018), 84 adults with symptoms of influenza-like illness (ILI) were included in the study. Blood, nasal swabs and nasopharyngeal aspirates were collected (Figure 1A). IAV infection was confirmed by PCR in 40 patients while 44 patients had other infections (IBV: 10, RSV: 6, other viruses: 2, bacteria: 3 and unknown aetiology: 23) despite presenting with similar symptoms (Figure 1B), consistent with inclusion based on ILI. Patients with bacterial infection, other viral infection or infections of unknown etiology were excluded from further analysis leaving 56 patients with confirmed IAV, IBV or RSV infection. The severity of disease in patients was categorized as “mild” or “moderate” (detailed description in Methods). Of the 56 IAV/IBV/RSV patients, 8 had moderately severe disease and 21 were hospitalized (Table 1). During the COVID-19 pandemic in 2020, we included 8 adults with PCR+ SARS-CoV-2 infection (7 with mild and 1 with moderate disease) for this study based on the availability of paired blood and NPA samples from the time of inclusion. Sixteen healthy controls (HCs) were included and sampled identically as patients.
IAV patients had a mean age of 59 (range: 20-98 years) and on average sought medical attention 3.3 days after the onset of symptoms (SD: 2.1 days) (Table 1). HCs had slightly lower mean age of 46 (range: 28-59 years). IAV patients had a mean Charlson co-morbidity index (CCI) of 2.46 (SD: 2.22), and 36 IAV patients had at least one underlying comorbidity (chronic heart/lung diseases, reduced lung function, kidney insufficiency, diabetes mellitus, asplenia/hyposplenia and malignancies).
Human IAV infection is characterized by an influx of CD11c+ cells into the nasopharynx
Blood from IAV patients yielded significantly fewer PBMCs/mL compared to HCs (Figure 1C). In contrast, 3-fold higher cell numbers were recovered from the nasopharynx of IAV patients compared to HCs (median 0.77 vs. 0.25 x 106 cells). In fact, 43% of IAV patients had more than 1 x 106 cells recovered from their NPA sample. Viability of PBMCs and NPA cells was variable across individuals, with no statistically significant differences between patients or HCs (Figures 1D). We determined the immune cell distribution in blood and the nasopharynx by flow cytometry (Supplementary Figure 1) on matched PBMC and NPA samples with a minimum 105 cells and ≥70% viability (n=22 IAV patients and n=16 HCs) to obtain high quality data (stained with identical panels and clones of antibodies). The frequencies of live CD45+ immune cells were increased in the NPA of patients as compared to HCs but remained similar in blood between the groups (data not shown). Among the immune cells, we found significantly higher frequencies of lineage (CD3, CD19, CD20, CD56, CD66abce) negative, HLA-DR+ cells, the compartment where monocytes and myeloid DCs can be identified, in both blood (p<0.0001), and the NPA (p<0.01) of IAV patients as compared to HCs (Figure 1E). Therefore, our data show that acute IAV infection results in an influx of monocytes and DCs to the nasopharynx.
Increased frequencies of classical and intermediate monocytes in the nasopharynx during IAV infection
To identify which monocyte subsets contributed to the changes observed during IAV infection, we analyzed the distribution of the different monocyte subsets (Figure 2A). As expected, CMs were the most frequent monocytes in blood in both patients and HCs, and remained comparable. However, in the nasopharynx of IAV patients as compared to HCs, blood CM frequencies were significantly increased (Figure 2B). Strikingly, IM frequencies were significantly elevated, in both blood and NPA of IAV patients (Figure 2C), while blood NCM appeared to be lower in patients compared to HCs (Figure 2D). Comparing frequencies of monocytes in blood and NPA in the same individual (Figure 2E) illustrates that IAV infection disrupts the pattern of monocyte distribution in the two anatomical compartments seen in healthy individuals- primarily in CMs and IMs. The frequencies of IMs in blood and NPA, and of the CMs in the NPA (solid lines) displayed greatest fluctuation from HCs (dashed line) at fewer days with symptoms (Figures 2F-G). This suggests that changes in the monocyte compartment may occur earlier in the course of infection but these findings should be interpreted with caution since they are based on a single acute phase sample from each patient. In future studies, longitudinal sampling of the same patient would allow for mapping kinetics of monocyte redistribution over the course of acute IAV infection and convalescence. We also compared the frequency of IMs in blood and nasopharynx with the age of IAV patients and found a negative correlation in blood (R=–0.55, p=0.0008) (Figure 2H) but a positive correlation in the nasopharynx (R=0.46, p=0.045) (Figure 2I). In HCs, age and IM frequencies were not significantly correlated in blood or NPA (data not shown). A subset of IAV patients (n=11) returned for sampling during convalescence (≥ 4 weeks after initial sampling). We observed that frequencies of CMs (in blood) and IMs (blood and NPA) in convalescent individuals were reduced and closer to values seen in HCs (Figure 2J). Collectively, we found that the increased immune cell presence in the nasopharynx during acute IAV infection could be, to a large extent, attributed to increased frequencies of IMs as well as CMs which normalized during convalescence. Acute IAV infection resulted in altered monocyte distribution, in particular at the site of infection and was more pronounced in older patients.
Dendritic cell subsets infiltrate the nasopharynx during IAV infection
DCs may also account for redistribution in the myeloid cell compartment. In HCs, the nasopharynx was virtually devoid of DCs, and CD1c+ MDCs, CD141+ MDCs and plasmacytoid DCs (PDCs) were only identified in a subset of HCs (Supplementary figure 1 and Figure 3A). In contrast, during acute IAV infection, all DC subsets were significantly increased in the nasopharynx of patients (Figure 3B-D). Meanwhile in blood, DC frequencies decreased significantly in IAV patients compared to HCs (Figures 3B-D). This reciprocal pattern is further illustrated by comparing the DC frequencies in blood and NPA in the same individual (Figure 3E). Looking over days of symptoms, the frequencies of DCs in blood remained lower in IAV patients (solid lines) as compared to HCs (dashed lines) throughout the sampling window (0-9 days with symptoms), with the most striking reduction evident at the earliest times (Figures 3F-G). Conversely in the nasopharynx, DC frequencies in patients fluctuated over time but remained higher than in controls. In contrast to IMs, we found an inverse correlation between age and the frequency of CD1c+ MDCs in NPA (R=–0.558, p=0.016) in IAV patients (Figure 3H). During convalescence, CD1c+ MDC frequencies were increased in blood and lowered in nasopharynx (Figure 3I). Taken together, we showed that the DC compartment also underwent substantial changes in both blood and the nasopharynx during acute IAV infection, potentially suggesting that DCs also received signals to induce migration and/or maturation.
Monocytes and DCs recruited to the human nasopharynx during IAV infection are more mature
We next analyzed the maturation status of DCs and monocytes in blood and nasopharynx samples from IAV patients and HCs (Figure 4A-D). Cells from HCs had low and comparable levels of surface HLA- DR in both blood and the nasopharynx (Figure 4A). In contrast, in IAV patients, monocytes and DCs in the nasopharynx expressed higher levels of HLA-DR than those in blood (Figure 4A and C). We also found that in IAV patients, nasopharyngeal CD1c+ MDCs and CD141+ MDCs expressed more CD86 than cells in blood (Figure 4D), while nasopharyngeal PDCs showed significant upregulation of CD83 during IAV infection as compared to blood PDCs (Figure 4E). When we compared maturation of CM and MDCs (i.e., CD86 expression) with viral RNA load (i.e., Ct values- number of cycles required to amplify viral RNA), we observed a significant inverse correlation (R=–0.696, p=0.0019), implying that higher viral RNA loads (low Ct value) were associated with increased maturation of CMs (Figure 4F). We also observed a positive correlation between CD86 expression on nasopharyngeal CMs, and both nasopharyngeal CD1c+ MDCs (R=0.735, p=0.0012) and CD141+ MDCs (R=0.832, p=0.0001) (Figures 4G and H, respectively). This indicates that in patients with greater viral RNA load who also have more mature CMs in the nasopharynx, there is a greater likelihood of finding mature MDCs in the nasopharynx as well.
Elevated cytokine levels observed in plasma and the nasopharyngeal secretions correlate with increased monocyte frequencies in respective compartment
Pronounced cytokinemia is a hallmark of severe influenza disease (18). In order to characterize the degree of inflammation in IAV patients, we measured local and systemic cytokine levels. In agreement with earlier reports (2, 3, 29, 30), we observed elevated levels of TNFα, IL-6 and IFNα in nasopharyngeal secretions (Figures 5A-C) as well as in plasma of IAV patients as compared to HCs (Figures 5D-F). We also observed elevated levels of plasma IL-10, IL-15 and IL-18 in IAV patients as compared to HCs (data not shown). We compared soluble TNFα levels against frequencies of monocytes and DCs (i.e., potential cellular sources) at the respective anatomical sites; and found a positive correlation between soluble TNFα and CM frequency, both in blood (Figure 5G); and in the nasopharynx of IAV patients (Figure 5H). However, such correlation was not observed for DCs in blood or NPA (data not shown). Interestingly, we also observed positive associations between age and level of TNFα in plasma and NPA in IAV patients (Figure 5I). Furthermore, we observed that plasma IFNα levels were positively correlated with Ct value (suggesting inverse correlation with viral RNA load) (Figure 5J). Patients with higher viral RNA loads had lower amounts of IFNα in circulation, which may suggest incomplete protection from infection. For a subset of patients, we assessed chemokines (CCL2, CCL3 and CCL7) in circulation and at the site of infection. Interestingly, in blood, elevated plasma levels of CCL2 correlated positively with IM frequencies, which were significantly elevated in IAV patients, indication a role for CCL2 in the changes to the IMs during infection (Figure 5K). Moreover, during convalescence, TNFα and IL-6 were also reduced both in the nasopharynx and in blood suggesting ablation of both local and systemic inflammation (Figures 5L).
Monocytes are a potent source of systemic TNFα during IAV infection
The localization of mature monocytes and DCs in the nasopharynx during IAV infection, and increased cytokine levels in both compartments led us to question whether the cells in circulation were directly involved in inflammation during ongoing infection, or primarily provided a cache of differentiated cells that can migrate to the site of infection. Limited by the number of viable cells obtained from the nasopharynx, we tested the functional response of circulating monocytes and DCs to in vitro stimulation with a TLR7/8 agonist which mimics ssRNA and then quantified the frequency of TNF-producing cells in each monocyte and DC subset (Figures 6A and B). We observed that blood monocytes and DCs from IAV patients produced TNFα spontaneously, in the absence of any external stimulus, with most of the cytokine coming from monocytes (CMs > IMs > NCMs) (Figure 6C). In contrast, blood monocytes and DCs from HCs only produced TNFα upon stimulation with TLR7/8L (Figure 6D). Importantly, while many of the IAV patients had cells producing TNFα spontaneously, monocytes and DCs remained responsive and had the potential for further increased frequency of TNF-producing cells in response to TLR7/8L stimulation (Figure 6D). In summary, our data suggest that during IAV infection, mature monocytes and DCs accumulate in the nasopharynx, and blood monocytes and DCs function as a general source of TNFα, potentially contributing to the systemic inflammatory effects accompanying influenza infections.
Nasopharyngeal aspirates allow assessment of in situ immune responses to mild infections with other respiratory viruses including SARS-CoV-2
In order to determine whether our findings were pathogen-specific to IAV or a reflection of immune responses to respiratory viral infections in general, we analyzed samples from patients with confirmed IBV, RSV or SARS-CoV-2 infections. IBV and RSV patients were closer to IAV patients in age (mean: 53 and 61 years, respectively) and days with symptoms (mean: 5.0 and 6.5 days respectively). In contrast, patients with mild SARS-CoV-2 were younger than IAV patients (mean: 44 years and 59 years respectively) and sought medical attention later: mean 13.2 days with symptoms (SARS-CoV-2) as compared to 3.3 days (IAV) (Table 2). Similar to IAV patients, patients with SARS-CoV-2 and RSV had a higher yield of cells in the nasopharynx than HCs (Figure 7A). The expansion in monocyte and DC frequencies in blood we saw in IAV patients, was not observed in other infections. Interestingly, in the nasopharynx, monocytes and DCs were significantly elevated during IAV infection, but not during IBV infection (Figure 7B). Within the monocyte compartment, we noted differences between the different pathogens. IMs were not increased during IBV or SARS-CoV-2 infection but were increased during RSV infection in both blood and the nasopharynx (Figure 7C). The DC compartment too showed differences between the groups in NPA. In blood, all DC subsets were decreased in all groups (Figure 7D). CD1c+ MDCs were increased in the nasopharynx only during IAV and SARS-CoV-2 infections. Interestingly, CD141+ MDCs were significantly increased in the nasopharynx only during IAV infection. PDCs in NPA were increased in all groups.
In plasma, individuals with SARS-CoV-2 infection had elevated levels of TNFα in comparison with HCs and IAV patients (Figure 7E). SARS-CoV-2 infection was associated with elevated levels of IL-6 compared to HCs, but not compared to the IAV and IBV groups. In NPA, only the IAV and IBV groups had increased levels of IL-6. Despite having increased frequencies of PDCs in NPA, nasopharyngeal levels of IFNα were not elevated in SARS-CoV-2 patients. Together, these data show different patterns of monocyte and DC engagement in the nasopharynx and in blood, and also between IAV, IBV, RSV and SARS-CoV-2, suggesting a requirement for further scrutiny.
DISCUSSION
In this study, we mapped monocyte and DC distribution and function in blood and in the nasopharynx (the initial site of infection) of patients with ongoing IAV, IBV, RSV or SARS-CoV-2 infections and healthy controls. Several studies have examined immune responses during severe flu, in particular, focusing on patients sampled during or immediately following the 2009 H1N1 pandemic (2, 3, 31, 40), or patients hospitalized with severe respiratory symptoms (2, 29–31, 40); and more recently, the immune dysregulation during severe COVID-19 (41–44). Here, we elucidated the innate myeloid cell composition and responses in a cohort with relatively mild symptoms, advanced age and underlying comorbidities, typical of seasonal influenza. Finally, we showed that our methods can be adapted to study immune responses in other viral infections, including the current COVID-19 pandemic.
Flow cytometric characterization revealed an influx of monocytes and DCs to the nasopharynx during infection, in line with previous reports in more severe IAV patients (3, 29, 30, 45). Interestingly, age appeared to differentially skew the IM response in IAV patients. Despite having higher baseline frequencies of circulating CD16-expressing monocytes (46), older patients have fewer IMs in blood during acute IAV infection and concurrently more IMs in the nasopharynx. The increase in IMs suggests a more inflammatory milieu at the site of infection in older patients, perhaps contributing to disease severity as previously suggested in pandemic influenza (2–5). Despite the general increase in IM frequencies, CMs remained the most frequent monocyte/DC subset in either compartment implying a functional role for CMs during IAV infection. Similar to studies in paediatric patients (29, 30), we also found that over the course of the illness, DCs appeared to migrate from blood to infiltrate the nasopharynx. However, older patients displayed a weaker recruitment of CD1c+ MDCs. In vitro studies have shown that monocytes can differentiate into type I IFN producing monocyte-derived DCs in response to IAV exposure (47). It is therefore possible that CMs recruited to the nasopharynx in older individuals preferentially differentiate into IMs rather than DCs. Moreover, CD1c+ DC function is critical for clearance of IAV infection (35). Therefore, diminished recruitment of CD1c+ MDCs and or reduced in situ DC differentiation (and therefore delayed or attenuated CD8+ T cell responses) may also contribute more severe disease in older patients.
During acute disease, monocytes present in/recruited to the nasopharynx likely contribute to sustained DC recruitment by secreting TNFα, CCL2, CCL3 and CCL7 locally (3, 29, 40). Due to the limited availability of cells in the NPA samples, demonstrating cytokine secretion from nasopharyngeal cells was not possible. However, we noted that in the nasopharynx, CMs, CD1c+ MDCs and CD141+ MDCs all appeared to mature by upregulating HLA-DR and/or CD86 expression as compared to HCs, suggesting local activation, likely in response to the inflammatory milieu and/or virus exposure. CD86, a critical costimulatory molecule and ligand for CD28 and CTLA-4, is upregulated on monocytes and DCs upon virus infection in vitro and in vivo (48). Lower expression of CD86 on circulating myeloid cells as compared to nasopharyngeal cells during IAV infection suggests migration of mature cells to the nasopharynx, and/or in situ differentiation of monocytes to induced-DCs (that may not upregulate maturation markers) in blood during acute IAV infection as previously speculated (47). A non-human primate model of chronic simian immunodeficiency virus (SIV) infection revealed sustained mobilization of MDCs from the bone marrow via blood to the intestinal mucosa, where MDCs remained activated and eventually underwent apoptosis (49). In IAV infection, monocytes and DCs are likely recruited to the nasopharynx for a shorter period of time, where they face a similar fate. Some DCs likely traffic antigens to the lymph node to support adaptive responses (21, 35).
Nasopharyngeal CMs were also more mature in patients with higher viral RNA loads and their recruitment and/or differentiation may be aided by local cytokine production (TNFα and IL-6). Key innate inflammatory cytokines, including TNFα, IL-6, IFNα, IL-10, IL-15, and IL-18, were all significantly elevated early during infection (days 1-5) compared to healthy controls. TNFα levels strongly correlated with the presence of CMs in both blood and nasopharynx; despite the significant expansion of IMs in blood, supporting previous studies (2, 14, 31, 50). We also demonstrated that be circulating CMs was the most frequent cellular source of TNFα among blood monocytes and DCs, while the TNFα at the site of infection may come from both CMs and the IMs (both subsets found in increased frequencies). Ex vivo, at steady state, human blood CMs and IMs have comparable TNFα secretion irrespective of the nature of stimulation (14). IAV infection, therefore, may skew monocyte subset function, in addition to monocyte distribution and maturation. The broad range of TNFα, IL-6 and IFNα responses seem among the IAV patients suggest that seasonal IAV strains are associated with milder cytokine responses than those observed during infections with 2009 H1N1 pandemic influenza or highly pathogenic avian H5N1 strains (3, 10, 31, 51). Within the setting of seasonal influenza however, older patients had more IMs in the NPA and displayed stronger TNFα responses, locally and systemically.
During IAV infection, the significant expansion of IMs in blood and the nasopharynx could be mediated by the presence of CCL2, a chemokine critical for monocyte chemotaxis (52). Upon CCL2 ligation, CCR2 facilitates CM and IM adhesion and transmigration from blood to tissues (53). Therefore, in IAV patients, CCL2 may contribute to recruitment of monocytes and DCs from circulation to the nasopharynx to aid inflammation (3, 50). Furthermore, elevated IFNα in circulation correlated with reduced viral RNA load, although we saw no evidence of this locally at the site of infection; suggesting that perhaps lower systemic IFN responses can indicate impaired control of viral replication in the nasopharynx, as previously reported (31). This finding also supports previous studies showing that attenuated RIG-I signalling impairs IFN responses in the elderly (54) and very young (55), both groups at high risk of influenza-associated mortality.
During convalescence, monocyte and DC frequencies, and cytokine levels normalized in blood and the nasopharynx to levels seen in healthy individuals. The local and systemic inflammation observed, therefore, appears limited to the acute phase of disease and wanes during convalescence. Unlike patients infected with the pandemic 2009 H1N1 IAV strain, where long lasting immune perturbations were documented (28); the repercussions of relatively mild seasonal IAV infections appear short-lived. This is also in contrast to in vivo studies which demonstrated prolonged DC activation and disruption of DC subset reconstitution in the lungs following IAV infection (56).
In April-May 2020, during the COVID-19 pandemic, the approach and methods described here allowed us to implement this study plan in a rapidly evolving pandemic. Although the patient cohorts were not identical, this endeavour proved the feasibility of using nasopharyngeal aspiration to assess local immune responses during respiratory viral infections. We observed that patients with IBV, RSV or SARS-CoV-2 also displayed an influx of cells to the nasopharynx during infection. Despite annual reoccurrence of IBV and RSV infections, detailed investigations of the innate cells at the site of infection in patients are quite rare (29, 30), and the technique we describe allows for minimally invasive longitudinal measurement of cellular and molecular aspects of infection. Patients with mild SARS-CoV-2 infection have not been studied as extensively as patients with severe or fatal COVID-19. The cell influx to the nasopharynx in COVID-19 patients did not appear to be due to monocytes but rather due to CD1c+ MDCs and PDCs, which were reduced in blood, in line with previous reports (41, 57). Also in line with previous studies (42–44), we found elevated levels of IL-6 in plasma, and also of TNFα. As different SARS-CoV-2 variants emerge, with concurrently increasing global vaccination rates, incomplete neutralization and/or protection and herd immunity may result in recurring, mild SARS-CoV- 2 infections. Understanding site-of-infection responses to mild SARS-CoV-2 infections is therefore increasingly relevant, and can be accomplished by the methods we describe in our study.
The present study augments our current understanding of the role of monocyte and DC subsets during human respiratory viral infections by highlighting unique dynamics and location-specific functions for each subset over the course of infection. Acute human IAV infection is characterized by an expansion of IMs and monocyte-mediated cytokinemia driving monocyte and DC migration from blood to the nasopharynx, the initial site of IAV infection. Older patients with comorbidities display a recruitment pattern skewed towards IMs and increased TNF rather than CD1c+ MDCs which may contribute to more severe disease and longer duration of hospitalization observed in these patients. IBV, RSV and mild SARS-CoV-2 infections elicit different patterns on immune cell recruitment to the nasopharynx as compared to IAV infections. We also illustrate the value of comparative high-resolution studies of immune cells in blood and at the site of infection, in order to fully understand their individual contributions to disease and how they orchestrate inflammation synergistically. Similar studies, carried out longitudinally across tissues, will aid resolution of these findings, and allow multivariate modelling of biomarkers of disease severity. Therapeutic approaches which allow selective modulation of monocyte and DC redistribution, maturation and cytokine/chemokine function may hold the key to reducing influenza-associated disease burden and mortality in the future.
METHODS
Study subjects
Patients seeking medical attention for influenza-like illness at the Emergency Department at the Karolinska University Hospital in Stockholm, Sweden during three consecutive influenza seasons (January - March) of 2016-2018 were recruited to the study following informed consent. The inclusion criteria for enrolment of patients were (1) age >18 years, (2) no known immunodeficiency, (3) not taking antibiotics, immunomodulatory or anti-inflammatory medication at time of inclusion, (4) presenting with fever and at least one of the following symptoms of influenza-like illness: cough, nasal congestion, headache or muscle ache. Convalescent samples were collected after at least 4 weeks, ensuring absence of respiratory symptoms in the prior week. Healthy volunteers were recruited and sampled similarly outside of influenza season. During the COVID-19 pandemic, additional patients were included between April - May 2020 at the Infectious Diseases ward at the Karolinska University Hospital or the Haga Outpatient Clinic (Haga Närakut) in Stockholm, Sweden, as well as mild/asymptomatic household contacts of patients with confirmed COVID-19 were screened by PCR and enrolled if positive. Clinical data were obtained from the patients and medical records (Table 1) and are extensively discussed in a previous publication (58). Total burden of comorbidities was assessed using the CCI (59).
The severity of disease was categorized using the respiratory domain of the sequential organ failure assessment score (SOFA) (60). In the absence of arterial partial pressure of oxygen (PaO2), peripheral transcutaneous haemoglobin saturation (SpO2) was used instead to calculate a modified SOFA score (mSOFA) (61). Fraction of inspired oxygen (FiO2) estimation based on O2 flow was done in accordance with the Swedish Intensive Care register definition (62). Mild disease was defined as PaO2/FiO2 (PFI) >53 kPa (>400 mmHg) or SpO2/FiO2 (SFI) >400. Moderate disease was defined as PFI >27-53 kPa (>200-400 mmHg) or SFI 235-400. The disease severity score is more extensively described in a previous publication (38).
Sample collection
Blood, nasal swabs and nasopharyngeal aspirates (NPA) were obtained from all patients (acute and convalescent phase samples) and healthy controls (Figure 1A). Briefly, up to 30mL venous blood was collected in Vacutainer® tubes containing EDTA, for blood counts and PBMC isolation. Nasopharyngeal swabs (Sigma Virocult®) were collected for diagnostic qPCR. NPA samples were collected into a vacuum trap by inserting a thin catheter through the naris, deep into the nasopharynx and applying gentle suction for 1-3 minutes. The vacuum trap and tubing were rinsed out with 3mL sterile PBS. All samples were processed within 2 hours of sampling.
Diagnostic tests to determine etiology of infection
Nasal swab samples were analyzed for influenza A virus (IAV), influenza B virus (IBV) and respiratory syncytial virus (RSV) by real-time PCR using the commercial Simplexa systemTM (63), as well as bacteria (by culture methods). The tests were performed at the Department of Clinical Microbiology, Karolinska University Laboratory as part of routine diagnostics for respiratory viral infections. Cycle threshold (Ct) values from the SimplexaTM assay were considered as semi-quantitative measures of virus levels in statistical analyses. Bacterial culture results were retrieved wherever available from the Department of Clinical Microbiology at the Karolinska University Hospital. Convalescent individuals and healthy controls were confirmed IAV– by qPCR. SARS-CoV-2 infection was diagnosed similarly using the GeneXpert SARS-CoV-2 detection system (Cepheid). Supplementary data on bacterial cultures were retrieved from the microbiology lab/clinical records.
Isolation of cells from blood and nasopharyngeal aspirates
Blood samples were centrifuged at 800g/10 min/room temperature (RT) and plasma was frozen at – 20°C. The blood volume was reconstituted with sterile PBS and PBMCs were obtained by density- gradient centrifugation using Ficoll-Paque Plus (GE Healthcare) after centrifugation at 900g/25 min/RT (without brake). NPA samples were centrifuged at 400g/5 min/RT and the supernatant was frozen at – 20°C. The cells were washed with 5mL sterile PBS to remove mucus, by filtering through a 70μm cell strainer followed by centrifugation for 400g/5 min/RT. Cell counts and viability were assessed using Trypan Blue (Sigma) exclusion and an automated Countess cell counter (Invitrogen) or manual counting using a light microscope (when the counter was unable to distinguish live and dead NPA cells). Cells were used fresh for flow cytometry (all NPA cells when samples contained >105 cells and had ≥70% viability; and at least 1x106 PBMCs) or frozen in RNALater (<105 NPA cells) (ThermoFisher). Excess PBMCs were cryopreserved in FBS (Gibco) + 10% DMSO (Sigma) and stored in liquid nitrogen.
Flow cytometry analysis
Phenotypic analysis was performed on NPA cells and PBMCs using Live/Dead Blue; lineage markers CD3 (SK7; BD), CD19 (HIB19; Biolegend), CD20 (L27, BD), CD45 (HI30; BD), CD56 (HCD56, BD) and CD66abce (TET2, Miltenyi Biotec); HLA-DR (TU36, Life technologies), CD14 (M5E2, BD), CD16 (3GE, Biolegend), CD11c (B-Ly6, BD), CD1c (AD5-8E7; Miltenyi), CD141 (AD5-14H12; Miltenyi), CD123 (7G3; BD); maturation markers CD83 (HB15e, Biolegend) and CD86 (2331; BD); adhesion marker CD62L (SK11, BD); and migration markers CCR2 (K036C2, Biolegend) and CCR7 (150503: BD); and fixed with 1% paraformaldehyde for flow cytometry. Samples were acquired on an LSRFortessa flow cytometer (BD Biosciences) and analyzed using FlowJo software v10 (Tree Star).
TNF-α release assay
TNFα secretion from fresh PBMCs in response to TLR stimulation was assessed using the TNF-α Secretion Assay-Detection Kit (PE) (130-091-268, Miltenyi Biotec) according to the manufacturer’s instructions. TNFα secretion was measured over 2 hours of incubation at 37°C with shaking (200 rpm) in the presence or absence of 1 μg/mL 3M-019 (7/8L) (Invivogen). Briefly, the capture antibody immobilized all secreted TNFα on the surface of the cell, and the detection antibody (PE-labelled) was incorporated into the above panel of antibodies to determine the cellular source of TNFα. Secretion of TNFα was quantified using flow cytometry.
ELISA and Luminex
Plasma samples and NPA supernatants were assayed for soluble markers using ELISA and Luminex. Human IFNα All Subtype ELISA was performed according to the manufacturer’s instructions (PBL Assay Science). TNFα and IL-6 ELISAs were performed using DuoSet® kits (R&D Systems). Luminex assays were performed using custom-designed 9-plex (IFNγ, IL-8, IL-10, IL-18, CCL2, CCL3, CCL7, IL-1β and IL-12p70) kit (R&D Systems) and analyzed on the Bio-Plex® 200 instrument (Bio-Rad).
Statistics
Data were analyzed using GraphPad Prism version 6.0 (GraphPad Software) and JMP®, version 14.2. (SAS Institute Inc., Cary, NC, 1989-2019). Differences between frequencies in IAV patients and HCs were assessed using nonparametric tests after assessing normality, using the Mann-Whitney test (at 95% confidence intervals). For comparisons between exposure conditions, two-way ANOVA with Sidak’s multiple comparisons test was applied. For comparisons between acute and convalescent phase data, Wilcoxon matched-pairs signed rank test was used. Bivariate and multivariate linear regression analysis was performed using JMP®, choosing Spearman’s rank correlation coefficient for nonparametric analyses. Differences between HCs or IAV patients and other patient groups were assessed using nonparametric tests after assessing normality, using the Kruskal-Wallis test with Dunn’s multiple comparisons test (at 95% confidence intervals). Data were considered significant at p<0.05.
Study approval
Informed consent was obtained from all patients and volunteers following verbal and written information. The study was approved by the Swedish Ethical Review Authority (No. 2015/1949-31/4) and performed according to the Declaration of Helsinki.
ADDITIONAL SUPPLEMENTARY MATERIAL
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
AUTHOR CONTRIBUTIONS
SV, NJ, AF and AS-S planned the study. SV, MY, SF-J and SL performed experiments. SF-J, BÖ, KS, AF and NJ included and sampled patients, and collected clinical data. JA provided relevant pseudonymized patient clinical data. SV, MY, SF-J, BÖ, SL, MA and AS-S analyzed data. SV and AS- S prepared figures and wrote the manuscript. All co-authors edited the manuscript.
DISCLOSURES
AS-S is a consultant to Astra-Zeneca on studies not related to the present study.
ACKNOWLEDGMENTS
We thank the patients and healthy volunteers who have contributed clinical material to this study. We would like to thank Adeline Mawa, Kevin Fallahi, Marcus Nordin, Roosa Vaitiniemi, Eric Åhlberg and the research nurses at the Emergency and Infectious Diseases Departments for technical assistance. This work was supported by grants to AS-S from the Swedish Research Council (VR), the Swedish Heart- Lung Foundation, the Bill and Melinda Gates Foundation, the Swedish Childhood Cancer Fund and Karolinska Institutet.
Footnotes
As part of submission to a journal, minor edits were made to Figures 1-7, and the text. Edits include updated statistics information, and incorporation of figures into running text in sequence.