Human influenza virus infection elicits distinct patterns of monocyte and dendritic cell mobilization in blood and the nasopharynx
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 controls (HCs) and in patients with mild to moderate infections (primarily influenza A virus [IAV]). As compared to 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 than HCs. In blood, monocytes were the most frequent cellular source of TNFα during IAV infection and remained responsive to additional stimulation with TLR7/8L. Immune responses in older patients skewed towards increased monocyte frequencies 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 infection groups. Using SomaScan, a high-throughput aptamer-based assay to study proteomic changes between patients and HCs, we found differential expression of innate immunity-related proteins in plasma and nasopharyngeal secretions of IAV and SARS-CoV-2 patients. 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.
Editor's evaluation
This study presents a valuable evaluation of the distribution of monocytes and dendritic cells in the blood and nasopharyngeal aspirates of patients with mild respiratory tract infections. They report solid differences between monocytes and dendritic cells and variation with patient age, showing that during influenza A virus infection, the classical and intermediate monocyte numbers increased both, in blood and nasopharyngeal aspirates while DCs increased in the nasopharyngeal aspirates only. The work will be of broad interest to immunologists, lung biologists and infection disease community.
https://doi.org/10.7554/eLife.77345.sa0Introduction
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 (WHO, 2018). 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 (Cole et al., 2017; Oshansky et al., 2014). 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 (Kuiken and Taubenberger, 2008; Kuiken et al., 2010). Other agents causing an influenza-like illness (ILI) 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) (Zhou et al., 2020).
IAV is primarily transmitted via inhalation of virus-containing aerosols or droplets, and mainly targets respiratory epithelial cells (Tellier, 2006; Sanders et al., 2011), with the nasopharynx being the initial site of virus replication. IAV generally remains localized to the airways, despite signs of systemic inflammation (Stramer et al., 2012). 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 (Sakabe et al., 2011; Wang et al., 2012). Monocytes and DCs shape the specificity and strength of the subsequent adaptive responses (Kwissa et al., 2014; Spörri and Reis e Sousa, 2005). In blood, three subsets of monocytes are found: the CD14+CD16– classical monocytes (CMs), the most frequent subset at steady state; the terminally differentiated CD14–CD16+ nonclassical monocytes (NCMs); and CD14+CD16+ intermediate monocytes (IMs) that represent CMs transitioning into NCMs (Boyette et al., 2017; Jakubzick et al., 2013; Patel et al., 2017). Blood IMs are readily identifiable at steady state and expand rapidly in response to inflammation, infection, or vaccination (Oshansky et al., 2014; Kwissa et al., 2014; Fingerle et al., 1993). 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 (Tisoncik et al., 2012). Conventionally, DCs were grouped as the CD1c+ myeloid DCs (MDCs) that excel at activating naïve T cells Hémont et al., 2013; the CD141+ MDCs that can efficiently can cross-present antigens via MHC-I Jongbloed et al., 2010; and the CD123+ plasmacytoid DCs (PDCs) that mediate type I IFN responses (Cella et al., 1999). More recent single-cell sequencing and unbiased clustering revealed six subtypes of DCs and four subtypes of monocytes (Villani et al., 2017). In this manuscript, we focus on cDC1 (previously CD141+ MDCs) and cDC2 (previously CD1c+ MDCs). We also speculate the roles for DC3s (inflammatory CD14-expressing CD1c+ MDCs, which seemingly arise in a GM-CSF-dependent manner) (Bourdely et al., 2020) and pre-DCs (AXL+SIGLEC6+CD11clowCD123+DCs that do not produce IFNα upon TLR9 stimulation but produce IL-12p70 at levels comparable to cDCs) (Villani et al., 2017; See et al., 2017). Monocytes and DCs vary in distribution and function, depending on the anatomical compartment (Baharom et al., 2016; Alcántara-Hernández et al., 2017). 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 (Smed-Sörensen et al., 2012; Diao et al., 2014; Baharom et al., 2015), delaying recovery and normalization of immune cell distribution and function (Jochems et al., 2018; Lichtner et al., 2011).
Studies have shown that monocytes and DCs are recruited to the nasopharynx following infection with 2009 H1N1pdm IAV strains (Oshansky et al., 2014), and in individuals hospitalized with severe influenza infections (Gill et al., 2008; Gill et al., 2005). Disease severity in hospitalized patients has been shown to correlate with (1) monocyte recruitment and increased levels of MCP3, IFNα-2, and IL-10 in the nasal compartment (Oshansky et al., 2014) and (2) strong TNF-producing monocytic responses in blood (Cole et al., 2017) and inflammatory, neutrophil-dominant patterns (Dunning et al., 2018). 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 (Smed-Sörensen et al., 2012; Diao et al., 2014; Jochems et al., 2018; Vangeti et al., 2019; McGill et al., 2008; Mount and Belz, 2010; GeurtsvanKessel et al., 2008), few studies compare responses between blood and the nasopharynx in human infections (Dunning et al., 2018; Vangeti et al., 2019; Segura et al., 2013). 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 (Gill et al., 2008). Studies on immune responses to mild SARS-CoV-2 infection have also primarily focused on blood and rarely the upper airways (Mick et al., 2020).
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 (Falck-Jones et al., 2021; Cagigi et al., 2021), 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 five consecutive influenza seasons (2016–2020), 121 adults with symptoms of ILI were included in the study. Blood, nasal swabs, and nasopharyngeal aspirates were collected (Figure 1A). IAV infection was confirmed by PCR in 70 patients while 51 patients had other infections (IBV: 10, RSV: 6, SARS-CoV-2: 35) despite presenting with similar symptoms (Figure 1B), consistent with inclusion based on ILI. No annual pattern of viral etiology was observed across the 5 years, except for SARS-CoV-2 which emerged in late 2019 and was only seen in patients from 2020 (Figure 1C). The severity of disease in patients was categorized as ‘mild’ or ‘moderate’ (detailed description in Methods). Of the 121 patients, 35 had moderately severe disease and 64 were hospitalized (Table 1). Of the SARS-CoV-2 patients, 15 were characterized as mild and 20 had moderate disease. Twenty healthy controls (HCs) were included and sampled identically.
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data1-v2.xlsx
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data2-v2.xlsx
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data3-v2.xlsx
-
Figure 1—source data 4
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data4-v2.xlsx
-
Figure 1—source data 5
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data5-v2.xlsx
-
Figure 1—source data 6
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig1-data6-v2.xlsx
IAV patients had a median age of 59 (range: 20–98 years) and sought medical attention after a median of 4 days following onset of symptoms (IQR: 2–5) (Table 1). HCs had a median age of 44 (range: 22–59 years). IAV patients had a median Charlson comorbidity index (CCI) of 1.5 (IQR: 0–4), 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 1D). In contrast, threefold higher cell numbers were recovered from the nasopharynx of IAV patients compared to HCs (median 0.77 vs. 0.25 × 106 cells). In fact, 43% of IAV patients had more than 1×106 cells recovered from their NPA sample. Viability of PBMCs and NPA cells was variable across individuals, with no statistically significant differences between patient or HC groups (Figure 1E). We determined the immune cell distribution in blood and the nasopharynx by flow cytometry (Figure 1—figure supplement 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 MDCs 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 CM and IM in the nasopharynx during IAV infection
To identify which monocyte and DC subsets contributed to the changes observed in the myeloid cell compartment during IAV infection, we compared the distribution of the different cell subsets. 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 2A). Strikingly, IM frequencies were significantly elevated, in both blood and NPA of IAV patients (Figure 2B), while blood NCM appeared to be lower in patients compared to HCs (Figure 2C). In HCs, the nasopharynx was virtually devoid of DCs– cDC2s, cDC1s, PDCs, and CD14+CD1c+ monocyte-derived DCs (mo-DCs) were only identified in a subset of HCs (Figure 1—figure supplement 1 and Figure 2D). In contrast, during acute IAV infection, all DC subsets were significantly decreased in the blood of patients (Figure 2E–H). Meanwhile in the nasopharynx, cDC2, cDC1, and PDC frequencies were significantly higher in IAV patients compared to HCs (Figure 2E–H). Additionally, we looked from the rare cDC1s with low CD11c expression (Figure 2—figure supplement 1). Despite very low numbers, the same trend of significant reduction in blood and increased presence in the nasopharynx was observed for these cells as well.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data1-v2.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data2-v2.xlsx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data3-v2.xlsx
-
Figure 2—source data 4
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data4-v2.xlsx
-
Figure 2—source data 5
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data5-v2.xlsx
-
Figure 2—source data 6
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data6-v2.xlsx
-
Figure 2—source data 7
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data7-v2.xlsx
-
Figure 2—source data 8
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data8-v2.xlsx
-
Figure 2—source data 9
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data9-v2.xlsx
-
Figure 2—source data 10
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data10-v2.xlsx
-
Figure 2—source data 11
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data11-v2.xlsx
-
Figure 2—source data 12
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data12-v2.xlsx
-
Figure 2—source data 13
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data13-v2.xlsx
-
Figure 2—source data 14
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data14-v2.xlsx
-
Figure 2—source data 15
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data15-v2.xlsx
-
Figure 2—source data 16
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data16-v2.xlsx
-
Figure 2—source data 17
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data17-v2.xlsx
-
Figure 2—source data 18
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data18-v2.xlsx
-
Figure 2—source data 19
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data19-v2.xlsx
-
Figure 2—source data 20
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data20-v2.xlsx
-
Figure 2—source data 21
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig2-data21-v2.xlsx
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). In contrast to IMs, we found an inverse correlation between age and the frequency of cDC2s in NPA (R=–0.558, p=0.016) in IAV patients (Figure 2K).
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), IMs (blood and NPA), and cDC2s (blood) in convalescent individuals were reduced and closer to values seen in HCs (Figure 2L). 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.
Monocytes and DCs recruited to the human nasopharynx during IAV infection are mature
We next analysed the maturation status of DCs and monocytes in blood and nasopharynx samples from IAV patients and HCs (Figure 3A–D). Cells from HCs had low and comparable levels of surface HLA-DR in both blood and the nasopharynx (Figure 3A). In contrast, in IAV patients, monocytes and DCs in the nasopharynx expressed higher levels of HLA-DR than those in blood (Figure 3A and C). We also found that in IAV patients, nasopharyngeal cDC2s, and cDC1s expressed more CD86 than cells in blood (Figure 3D), while nasopharyngeal PDCs showed significant upregulation of CD83 during IAV infection as compared to blood PDCs (Figure 3E). When we compared maturation of CM and cDCs (i.e., CD86 expression) with viral RNA load (i.e., cycle threshold [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 3F). We also observed a positive correlation between CD86 expression on nasopharyngeal CMs, and both nasopharyngeal cDC2s (R=0.735, p=0.0012) and cDC1s (R=0.832, p=0.0001) (Figure 3G 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.
-
Figure 3—source data 1
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data1-v2.xlsx
-
Figure 3—source data 2
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data2-v2.xlsx
-
Figure 3—source data 3
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data3-v2.xlsx
-
Figure 3—source data 4
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data4-v2.xlsx
-
Figure 3—source data 5
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data5-v2.xlsx
-
Figure 3—source data 6
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data6-v2.xlsx
-
Figure 3—source data 7
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data7-v2.xlsx
-
Figure 3—source data 8
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig3-data8-v2.xlsx
Elevated cytokine levels observed in plasma and nasopharyngeal secretions correlate with increased monocyte frequencies in the respective compartment
Cytokinemia is a hallmark of severe influenza disease (Tisoncik et al., 2012). In order to characterize the degree of inflammation in IAV patients, we measured local and systemic cytokine levels. In agreement with earlier reports (Cole et al., 2017; Oshansky et al., 2014; Gill et al., 2008; Gill et al., 2005), we observed elevated levels of TNFα, IL-6, and IFNα in nasopharyngeal secretions (Figure 4A–C) as well as in plasma of IAV patients as compared to HCs (Figure 4D–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 4G) and in the nasopharynx of IAV patients (Figure 4H). 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 4I). Furthermore, we observed that plasma IFNα levels were positively correlated with Ct value (suggesting inverse correlation with viral RNA load) (Figure 4J). 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 (R=0.639, p=0.048), which were significantly elevated in IAV patients, indicating a role for CCL2 in the changes to the IMs during infection (Figure 4K). 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 (Figure 4L).
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data1-v2.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data2-v2.xlsx
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data3-v2.xlsx
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data4-v2.xlsx
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data5-v2.xlsx
-
Figure 4—source data 6
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data6-v2.xlsx
-
Figure 4—source data 7
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data7-v2.xlsx
-
Figure 4—source data 8
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data8-v2.xlsx
-
Figure 4—source data 9
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data9-v2.xlsx
-
Figure 4—source data 10
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data10-v2.xlsx
-
Figure 4—source data 11
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data11-v2.xlsx
-
Figure 4—source data 12
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig4-data12-v2.xlsx
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 subset and in cDCs (Figure 5A 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 5C). In contrast, blood monocytes and DCs from HCs only produced TNFα upon stimulation with TLR7/8L (Figure 5D). 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 5D). 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 specific to IAV or a reflection of immune responses to respiratory viral infections in general, we analysed samples from patients with confirmed IBV, RSV, or SARS-CoV-2 infections. All patient groups were relatively close in age (Table 1). SARS-CoV-2 patients tended to seek medical attention later: median 12 days with symptoms (SARS-CoV-2) as compared to 4 days (IAV), 5 days (IBV), and 6 days (RSV) (Table 1). Similar to IAV patients, patients with SARS-CoV-2 and RSV had a higher yield of cells in the nasopharynx than HCs (Figure 6A). The expansion in monocyte (IM) 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 6B). 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 6C). The DC compartment too showed differences between the groups in NPA. In blood, all DC subsets were decreased in all groups (Figure 6D). cDC2s were increased in the nasopharynx only during IAV and SARS-CoV-2 infections. Interestingly, cDC1s were significantly increased in the nasopharynx only during IAV infection. PDCs in NPA were increased in all groups.
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data1-v2.xlsx
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data2-v2.xlsx
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data3-v2.xlsx
-
Figure 6—source data 4
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data4-v2.xlsx
-
Figure 6—source data 5
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data5-v2.xlsx
-
Figure 6—source data 6
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data6-v2.xlsx
-
Figure 6—source data 7
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data7-v2.xlsx
-
Figure 6—source data 8
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data8-v2.xlsx
-
Figure 6—source data 9
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data9-v2.xlsx
-
Figure 6—source data 10
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data10-v2.xlsx
-
Figure 6—source data 11
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data11-v2.xlsx
-
Figure 6—source data 12
- https://cdn.elifesciences.org/articles/77345/elife-77345-fig6-data12-v2.xlsx
In plasma, individuals with SARS-CoV-2 infection had elevated levels of TNFα in comparison with HCs and IAV patients (Figure 6E). 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.
Innate immunity-related proteins are differentially expressed in plasma and nasopharyngeal secretions of IAV and SARS-CoV-2 infected patients
To further emphasize the distinctive aspects of immune responses between anatomical compartments, and across infections, we utilized SomaScan, an aptamer-based proteomics platform assay to provide a broad assessment of the human proteome (n=7288 analyses targeting 6595 unique human protein measurements) between acute (IAV and SARS-CoV-2) patients and HCs. Based on PCA, plasma and nasopharyngeal proteome composition displayed an overlap regardless of infection or healthy status, separating mainly on the basis of biological site (Figure 7A). In IAV patients compared to HCs, more proteins were differentially expressed in blood (n=548) than in the nasopharynx (n=281) (Figure 7B), with an overlap of only 8.4% between anatomical sites (n=60 proteins) (Figure 7C–D). Interestingly, in SARS-CoV-2 patients compared to HCs, the nasopharynx (n=1868) had considerably more differentially expressed proteins than blood (n=693), over 1500 of which were unique to the nasopharynx. Also, in both anatomical compartments, many differentially expressed proteins were unique to IAV or SARS-CoV-2 (Figure 7D, right panels). Pathway analysis was performed to identify significantly enriched GO biological processes and pathways across the four t test models. The number and type of enriched pathways varied across models (Figure 7E and Supplementary file 2). The largest cluster, primarily mapped to IAV (plasma and NPA), and to a lesser extent, to SARS-CoV-2 NPA, where proteins mapped to host immune responses and innate immunity, supporting the data described elsewhere in this study. Notably, pathways containing TNF, IL-6, ISG15, IL-18R, CCL7, CXCL10 (IP-10), CXCL11, GZMB, SEMA4A, S100A8, S100A9 were all associated with the nasopharynx in IAV patients, with IL-6 pathways showing significant association with nasopharynx in SARS-CoV-2 patients. Both NPA models had the greatest number of significantly enriched pathways (IAV = 51 and SARS-CoV-2=64) compared to the plasma models (IAV = 12 and SARS-CoV-2=19).
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 in HCs. Several studies have examined immune responses during severe influenza, in particular, focusing on patients sampled during or immediately following the 2009 H1N1 pandemic (Cole et al., 2017; Oshansky et al., 2014; Dunning et al., 2018; Marion et al., 2016), or patients hospitalized with severe respiratory symptoms (Cole et al., 2017; Gill et al., 2008; Gill et al., 2005; Dunning et al., 2018; Marion et al., 2016); and more recently, the immune dysregulation during severe COVID-19 (Zingaropoli et al., 2021; Del Valle et al., 2020; Merad and Martin, 2020; Mudd et al., 2020). 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 (Oshansky et al., 2014; Gill et al., 2008; Gill et al., 2005; Herold et al., 2006). Interestingly, age appeared to differentially skew the CD14+CD16+ IM response in IAV patients. Despite having higher baseline frequencies of circulating CD16-expressing monocytes (The Milieu Intérieur Consortium et al., 2018), older patients have fewer CD14+CD16+ 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 (Cole et al., 2017; Oshansky et al., 2014; Kuiken and Taubenberger, 2008; Kuiken et al., 2010). Despite the general increase in CD14+CD16+ IM frequencies, CD14+CD16– 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 (Gill et al., 2008; Gill et al., 2005), 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 mo-DCs in response to IAV exposure (Cao et al., 2012). It is therefore possible that CMs recruited to the nasopharynx in older individuals preferentially differentiate into IMs rather than mo-DCs. Moreover, cDC2 function is critical for clearance of IAV infection (GeurtsvanKessel et al., 2008). Therefore, diminished recruitment of cDC2s 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.
At the time this study was conceived, DC3s were not yet identified. Guided by a canonical approach and limited by nasopharyngeal cell recovery, we were unable to capture the true heterogeneity of monocyte and DC subsets within the human airways in this study. From later reports establishing DC3 as a distinct DC subset arising from a specific progenitor (Bourdely et al., 2020; Cytlak et al., 2020; Dutertre et al., 2019), especially those describing dysregulated immune responses during SARS-CoV-2 infection (Winheim et al., 2021; Kvedaraite et al., 2021; Affandi et al., 2021), it appears that DC3 may be depleted during severe infection. As the pattern of overall DC depletion from circulation is shared between IAV and SARS-CoV-2 infections, it is plausible that DC3 may be present in the nasopharynx during acute IAV infection. As our comprehension of monocyte and DC heterogeneity expands, it is vital that studies like the present are carried out to accurately map the fate of newly discovered subsets in the airways.
During acute disease, monocytes present in/recruited to the nasopharynx likely contribute to sustained DC recruitment by secreting TNFα, CCL2, CCL3, and CCL7 locally (Oshansky et al., 2014; Gill et al., 2008; Marion et al., 2016). 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, cDC2s, and cDC2s all appeared to mature by upregulating HLA-DR and/or CD86 expression as compared to cells in 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 (Hubo et al., 2013). 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 (Cao et al., 2012). A non-human primate model of chronic simian immunodeficiency virus infection revealed sustained mobilization of MDCs from the bone marrow via blood to the intestinal mucosa, where MDCs remained activated and eventually underwent apoptosis (Wijewardana et al., 2013). 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 (Cella et al., 1999; GeurtsvanKessel et al., 2008).
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 HCs. 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 (Cole et al., 2017; Boyette et al., 2017; Dunning et al., 2018; Sprenger et al., 1996). We also demonstrated that the circulating CMs were 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 (Boyette et al., 2017). 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 seen 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 (Oshansky et al., 2014; Sakabe et al., 2011; Dunning et al., 2018; Gerlach et al., 2013). 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 CD14+CD16+ IMs in blood and the nasopharynx could be mediated by the presence of CCL2, a chemokine critical for monocyte chemotaxis (Dean et al., 2008). Upon CCL2 ligation, CCR2 facilitates CM and IM adhesion and transmigration from blood to tissues (Tsou et al., 2007). Therefore, in IAV patients, CCL2 may contribute to recruitment of monocytes and DCs from circulation to the nasopharynx to aid inflammation (Oshansky et al., 2014; Sprenger et al., 1996). 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 (Dunning et al., 2018). This finding also supports previous studies showing that attenuated RIG-I signalling impairs IFN responses in the elderly (Molony et al., 2017) and very young (Marr et al., 2014), 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 during the pandemic years (2009–2010), where long-lasting immune perturbations were documented Lichtner et al., 2011; 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 in mice (Strickland et al., 2014).
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 (Gill et al., 2008; Gill et al., 2005), 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 cDC2s and PDCs, which were reduced in blood, in line with previous reports (Zingaropoli et al., 2021; Peruzzi et al., 2020). Also in line with previous studies (Del Valle et al., 2020; Merad and Martin, 2020; Mudd et al., 2020), 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.
Finally, to further highlight the importance of studies like ours, we employed a high sensitivity, high-throughput proteomics platform to quantify differences in the nasopharyngeal proteome as compared to blood, in both IAV and SARS-CoV-2 infections. Additionally, this approach allows identification of proteins that cannot be identified by conventional assays that may be limited in breadth of proteome coverage or by limits of detection. We found differentially expressed proteins unique to blood or nasopharynx, and also to IAV or SARS-CoV-2. Preliminary pathway analysis reveals sizeable changes within the nasopharyngeal proteome during infection, especially during SARS-CoV-2 infection. In IAV patients, in line with previous reports (Marion et al., 2016) several innate immune pathway proteins were differentially enriched in the nasopharyngeal proteome, including TNF and IL-6, as seen from ELISA data.
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 modeling 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 ILI 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 ILI: cough, nasal congestion, headache, or muscle ache. For the current study, only data from patients who had confirmed IAV, IBV, RSV, or (later) SARS-CoV-2 infections were used. 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 and 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 (Falck-Jones et al., 2021). Total burden of comorbidities was assessed using the CCI (Charlson et al., 1987).
The severity of disease was categorized using the respiratory domain of the sequential organ failure assessment (SOFA) score (Vincent et al., 1996). 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) (Grissom et al., 2010). Fraction of inspired oxygen (FiO2) estimation based on O2 flow was done in accordance with the Swedish Intensive Care register definition (Intensivvårdsregistret, 2018). 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 (Falck-Jones et al., 2021).
Sample collection
Blood, nasal swabs, and nasopharyngeal aspirates (NPA) were obtained from all patients (acute and convalescent phase samples) and HCs (Figure 1A). Briefly, up to 30 mL 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 min. The vacuum trap and tubing were rinsed out with 3 mL sterile PBS. All samples were processed within 2 hr of sampling.
Diagnostic tests to determine etiology of infection
Nasal swab samples were analysed for IAV, IBV, and RSV by real-time PCR using the commercial Simplexa system (Svensson et al., 2014), 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. Ct values from the Simplexa 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 HCs 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 800 g/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 900 g/25 min/RT (without brake). NPA samples were centrifuged at 400 g/5 min/RT and the supernatant was frozen at –20°C. The cells were washed with 5 mL sterile PBS to remove mucus, by filtering through a 70 μm cell strainer followed by centrifugation for 400 g/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 1×106 PBMCs) or frozen in RNALater (<105 NPA cells) (Thermo Fisher). 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 analysed 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 hr 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 analysed on the Bio-Plex 200 instrument (Bio-Rad).
SomaScan assay
Proteomic profiling of NPA and EDTA plasma samples was conducted using SomaScan v4.1 (SomaLogic). The SomaScan proteomic platform has been described previously (Rohloff et al., 2014). Briefly, the platform utilizes DNA-based binding reagents (slow off-rate modified aptamers) to quantify by fluorescence the availability of binding epitopes (which represents protein abundance, shape, and charge) for over 7000 protein targets with high specificity and limits of detection comparable to antibody-based assays. The assay captures both high- and low-abundance proteins over a dynamic range of detection of approximately 10 logs (Williams et al., 2022). Normalization of SomaScan data was performed on plasma and NPA samples separately for all statistical analysis except PCA visualization. This was done because of the dramatically different signal distributions observed for the two anatomical compartments. Median scaling of samples to the global population median was performed. Additional normalization methods were evaluated and resulted in similar outcomes and performance (data not shown). Only samples derived from patients with a disease severity ≤ 4 (moderate severity) were evaluated. Outlier samples were detected and removed based on correlation with other samples. All pairwise correlations for a given sample with other samples in the same tissue type were calculated. If the median pairwise correlation for a given sample fell below 3 robust standard deviations (MAD) of the median pairwise correlation for all samples that sample was considered an outlier and removed from t test and pathway analysis (see below).
Statistics
Data were analysed using GraphPad Prism v6.0 (GraphPad Software), JMP, v14.2 (SAS Institute Inc, Cary, NC, 1989–2019) and R (v4.0.3). 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. Principal components analysis was performed to visualize the global impact of NPA and plasma tissue type and the impact of either IAV or SARS-CoV-2 infection. Non-parametric ellipses cover 95% of samples within a group; t tests were performed comparing the acute phase samples with HCs. Test populations were stratified by tissue type and infectious agent. Differentially expressed proteins were selected based on a p<0.05 and absolute fold-change >1.5 criteria. In cases where individual SOMAmers may interrogate more than one gene product (EnsembleGene) the different genes were separated for Venn diagrams and pathway analysis. Final gene lists are made up of unique gene lists because of the one-to-many and many-to-one relationship between SOMAmer reagents and protein/gene products.
Pathway analysis
Analysis of enriched pathways was performed for the four t test comparisons representing IAV and SARS-CoV-2 from either NPA or plasma specimens. Differentially expressed proteins detected by t test (p<0.05) and absolute fold-change >1.5 were used as input for the g:profiler pathway enrichment method using gprofiler2 v0.2.1 (Kolberg et al., 2020; Raudvere et al., 2019) in R. Annotated gene sets representing Gene Ontology biological processes, Hallmark pathways, C6 Oncogenic signature gene sets, and Canonical Pathways (incl. Biocarta, KEGG, PID, Reactome, and Wikipathways) were downloaded from the Molecular Signatures Database v7.5.1 for use as target pathways (Liberzon et al., 2011). Due to the semi-targeted nature of the SomaScan protein quantification method, a unique list of all proteins (ensembleGeneSymbol) interrogated on the platform (n=6370) were used as a statistical null background. An FDR threshold of q<0.05 was used to select enriched gene sets for further analysis. The Cytoscape Enrichment Map plugin v3.3.4 (Merico et al., 2010) was used for clustering and visualization of significantly enriched gene sets detected for each cell line xenograft model. Default clustering parameters were applied. Singleton gene sets were removed from visualization if only a single t test comparison contributed to enrichment and the gene set was not connected to a larger functional cluster.
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.
Appendix 1
Data availability
Source Data files have been provided for figures 1-7 and Appendix 2_Table 2 as separate excel files.
References
-
Cd169 defines activated CD14+ monocytes with enhanced CD8+ T cell activation capacityFrontiers in Immunology 12:697840.https://doi.org/10.3389/fimmu.2021.697840
-
A new method of classifying prognostic comorbidity in longitudinal studies: development and validationJournal of Chronic Diseases 40:373–383.https://doi.org/10.1016/0021-9681(87)90171-8
-
Functional monocytic myeloid-derived suppressor cells increase in blood but not airways and predict COVID-19 severityThe Journal of Clinical Investigation 131:e144734.https://doi.org/10.1172/JCI144734
-
The novel subset of CD14+/CD16+ blood monocytes is expanded in sepsis patientsBlood 82:3170–3176.
-
Clearance of influenza virus from the lung depends on migratory langerin+cd11b- but not plasmacytoid dendritic cellsThe Journal of Experimental Medicine 205:1621–1634.https://doi.org/10.1084/jem.20071365
-
Mobilization of plasmacytoid and myeloid dendritic cells to mucosal sites in children with respiratory syncytial virus and other viral respiratory infectionsThe Journal of Infectious Diseases 191:1105–1115.https://doi.org/10.1086/428589
-
Differential recruitment of dendritic cells and monocytes to respiratory mucosal sites in children with influenza virus or respiratory syncytial virus infectionThe Journal of Infectious Diseases 198:1667–1676.https://doi.org/10.1086/593018
-
A modified sequential organ failure assessment score for critical care triageDisaster Medicine and Public Health Preparedness 4:277–284.https://doi.org/10.1001/dmp.2010.40
-
Human blood MDC subsets exhibit distinct TLR repertoire and responsivenessJournal of Leukocyte Biology 93:599–609.https://doi.org/10.1189/jlb.0912452
-
Human CD141+ (BDCA-3) + dendritic cells (DCs) represent a unique myeloid DC subset that cross-presents necrotic cell antigensThe Journal of Experimental Medicine 207:1247–1260.https://doi.org/10.1084/jem.20092140
-
Pathology of human influenza revisitedVaccine 26 Suppl 4:D59–D66.https://doi.org/10.1016/j.vaccine.2008.07.025
-
Comparative pathology of select agent influenza A virus infectionsVeterinary Pathology 47:893–914.https://doi.org/10.1177/0300985810378651
-
Molecular signatures database (msigdb) 3.0Bioinformatics 27:1739–1740.https://doi.org/10.1093/bioinformatics/btr260
-
Protective influenza-specific CD8 T cell responses require interactions with dendritic cells in the lungsThe Journal of Experimental Medicine 205:1635–1646.https://doi.org/10.1084/jem.20080314
-
Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophagesNature Reviews Immunology 20:355–362.https://doi.org/10.1038/s41577-020-0331-4
-
Mouse models of viral infection: influenza infection in the lungMethods in Molecular Biology 595:299–318.https://doi.org/10.1007/978-1-60761-421-0_20
-
Mucosal immune responses predict clinical outcomes during influenza infection independently of age and viral loadAmerican Journal of Respiratory and Critical Care Medicine 189:449–462.https://doi.org/10.1164/rccm.201309-1616OC
-
The fate and lifespan of human monocyte subsets in steady state and systemic inflammationThe Journal of Experimental Medicine 214:1913–1923.https://doi.org/10.1084/jem.20170355
-
G: profiler: a web server for functional enrichment analysis and conversions of gene Lists (2019 update)Nucleic Acids Research 47:W191–W198.https://doi.org/10.1093/nar/gkz369
-
Cytokine production by primary human macrophages infected with highly pathogenic H5N1 or pandemic H1N1 2009 influenza virusesThe Journal of General Virology 92:1428–1434.https://doi.org/10.1099/vir.0.030346-0
-
Respiratory epithelial cells in innate immunity to influenza virus infectionCell and Tissue Research 343:13–21.https://doi.org/10.1007/s00441-010-1043-z
-
Similar antigen cross-presentation capacity and phagocytic functions in all freshly isolated human lymphoid organ-resident dendritic cellsThe Journal of Experimental Medicine 210:1035–1047.https://doi.org/10.1084/jem.20121103
-
Selective induction of monocyte and not neutrophil-attracting chemokines after influenza A virus infectionThe Journal of Experimental Medicine 184:1191–1196.https://doi.org/10.1084/jem.184.3.1191
-
Sensitive detection assays for influenza RNA do not reveal viremia in US blood donorsThe Journal of Infectious Diseases 205:886–894.https://doi.org/10.1093/infdis/jir863
-
Performance of the simplexaperformance of the simplexaScandinavian Journal of Infectious Diseases 46:825–831.https://doi.org/10.3109/00365548.2014.946444
-
Review of aerosol transmission of influenza A virusEmerging Infectious Diseases 12:1657–1662.https://doi.org/10.3201/eid1211.060426
-
Into the eye of the cytokine stormMicrobiology and Molecular Biology Reviews 76:16–32.https://doi.org/10.1128/MMBR.05015-11
-
Critical roles for CCR2 and MCP-3 in monocyte mobilization from bone marrow and recruitment to inflammatory sitesThe Journal of Clinical Investigation 117:902–909.https://doi.org/10.1172/JCI29919
-
The SOFA (sepsis-related organ failure assessment) score to describe organ dysfunction/failureIntensive Care Medicine 22:707–710.https://doi.org/10.1007/BF01709751
-
A proteomic surrogate for cardiovascular outcomes that is sensitive to multiple mechanisms of change in riskScience Translational Medicine 14:eabj9625.https://doi.org/10.1126/scitranslmed.abj9625
Article and author information
Author details
Funding
Vetenskapsrådet
- Anna Smed-Sörensen
Hjärt-Lungfonden
- Anna Smed-Sörensen
Bill & Melinda Gates Foundation
- Anna Smed-Sörensen
Barncancerfonden
- Anna Smed-Sörensen
Karolinska Institutet
- Anna Smed-Sörensen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
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 & Melinda Gates Foundation, the Swedish Childhood Cancer Fund, and Karolinska Institutet.
Ethics
Human subjects: 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.
Copyright
© 2023, Vangeti et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,880
- views
-
- 231
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
- Neuroscience
Acute retinal ischemia and ischemia-reperfusion injury are the primary causes of retinal neural cell death and vision loss in retinal artery occlusion (RAO). The absence of an accurate mouse model for simulating the retinal ischemic process has hindered progress in developing neuroprotective agents for RAO. We developed a unilateral pterygopalatine ophthalmic artery occlusion (UPOAO) mouse model using silicone wire embolization combined with carotid artery ligation. The survival of retinal ganglion cells and visual function were evaluated to determine the duration of ischemia. Immunofluorescence staining, optical coherence tomography, and haematoxylin and eosin staining were utilized to assess changes in major neural cell classes and retinal structure degeneration at two reperfusion durations. Transcriptomics was employed to investigate alterations in the pathological process of UPOAO following ischemia and reperfusion, highlighting transcriptomic differences between UPOAO and other retinal ischemia-reperfusion models. The UPOAO model successfully replicated the acute interruption of retinal blood supply observed in RAO. 60 min of Ischemia led to significant loss of major retinal neural cells and visual function impairment. Notable thinning of the inner retinal layer, especially the ganglion cell layer, was evident post-UPOAO. Temporal transcriptome analysis revealed various pathophysiological processes related to immune cell migration, oxidative stress, and immune inflammation during the non-reperfusion and reperfusion periods. A pronounced increase in microglia within the retina and peripheral leukocytes accessing the retina was observed during reperfusion periods. Comparison of differentially expressed genes (DEGs) between the UPOAO and high intraocular pressure models revealed specific enrichments in lipid and steroid metabolism-related genes in the UPOAO model. The UPOAO model emerges as a novel tool for screening pathogenic genes and promoting further therapeutic research in RAO.
-
- Immunology and Inflammation
The gram-negative bacterium lipopolysaccharide (LPS) is frequently administered to generate models of systemic inflammation. However, there are several side effects and no effective treatment for LPS-induced systemic inflammation. PEGylated PDZ peptide based on zonula occludens-1 (ZO-1) was analyzed for its effects on systemic inflammation induced by LPS. PDZ peptide administration led to the restoration of tissue injuries (kidney, liver, and lung) and prevented alterations in biochemical plasma markers. The production of pro-inflammatory cytokines was significantly decreased in the plasma and lung BALF in the PDZ-administered mice. Flow cytometry analysis revealed the PDZ peptide significantly inhibited inflammation, mainly by decreasing the population of M1 macrophages, and neutrophils (immature and mature), and increasing M2 macrophages. Using RNA sequencing analysis, the expression levels of the NF-κB-related proteins were lower in PDZ-treated cells than in LPS-treated cells. In addition, wild-type PDZ peptide significantly increased mitochondrial membrane integrity and decreased LPS-induced mitochondria fission. Interestingly, PDZ peptide dramatically could reduce LPS-induced NF-κB signaling, ROS production, and the expression of M1 macrophage marker proteins, but increased the expression of M2 macrophage marker proteins. These results indicated that PEGylated PDZ peptide inhibits LPS-induced systemic inflammation, reducing tissue injuries and reestablishing homeostasis, and may be a therapeutic candidate against systemic inflammation.