Abstract
Rhinoviruses (RV) and inhaled allergens, such as house dust mite (HDM) are the major agents responsible for asthma onset, its life-threatening exacerbations and progression to severe disease. The role of severe acute respiratory syndrome coronavirus (SARS-CoV-2) in exacerbations of asthma or the influence of preexisting viral or allergic airway inflammation on the development of coronavirus disease 2019 (COVID-19) is largely unknown. To address this, we compared molecular mechanisms of HDM, RV and SARS-CoV-2 interactions in experimental RV infection in patients with asthma and healthy individuals. RV infection was sensed via retinoic acid-inducible gene I (RIG-I) helicase, but not via NLR family pyrin domain containing 3 (NLRP3), which led to subsequent apoptosis-associated speck like protein containing a caspase recruitment domain (ASC) recruitment, oligomerization and RIG-I inflammasome activation. This phenomenon was augmented in bronchial epithelium in patients with asthma, especially upon pre-exposure to HDM, which itself induced a priming step, pro-IL-1β release and early inhibition of RIG-I/TANK binding kinase 1/IκB kinase ε /type I/III interferons (RIG-I/TBK1/IKKε/IFN-I/III) responses. Excessive activation of RIG-I inflammasomes was partially responsible for the alteration and persistence of type I/III IFN responses, prolonged viral clearance and unresolved inflammation in asthma. RV/HDM-induced sustained IFN I/III responses initially restricted SARS- CoV-2 replication in epithelium of patients with asthma, but even this limited infection with SARS-CoV-2 augmented RIG-I inflammasome activation. Timely inhibition of the epithelial RIG-I inflammasome and reduction of IL-1β signaling may lead to more efficient viral clearance and lower the burden of RV and SARS- CoV-2 infection.
1. Introduction
Asthma is one of the most common chronic inflammatory lung diseases affecting more than 5% of the global population1. Its pathogenesis and clinical presentation is complex2, with a common feature of susceptibility to exacerbations leading to loss of disease control, hospitalizations, and in some cases, progressive loss of lung function3, 4. Exacerbations of asthma are most often caused by common respiratory viruses5, 6, with rhinoviruses (RV) responsible for up to 80% of asthma attacks5. RVs that have been initially considered as benign viruses, now are also linked to the early-life development of asthma, severe bronchiolitis in infants and fatal pneumonia in elderly and immunocompromised patients7–10. Likewise, human coronaviruses have not been strongly linked with asthma pathology11. However, the current pandemic of severe acute respiratory syndrome coronavirus (SARS-CoV-2) has been challenging this view, resulting in contradictory observations of asthma being considered a risk factor for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19) severity12–14 or constituting a protection from the disease15, 16. Another important factor for asthma development and exacerbations is exposure to inhaled allergens. House dust mite (HDM) is the most significant source of perennial allergens worldwide. HDM sensitization is found in around 50%-85% of patients with asthma, and HDM exposure correlates with asthma severity17, 18. There are strong epidemiological links between RV infections, allergen exposure and sensitization on the risk of asthma development and the rates of exacerbations9, 19. Children with early life RV-induced wheezing and aeroallergen sensitization have an extremely high incidence of asthma in later years9. Combination of virus detection in the airways with the high allergen exposure markedly increases the risk of hospital admission20. In line with this, HDM immunotherapy significantly reduces risk of asthma exacerbations21. It has been also recently suggested that allergen exposure might influence SARS-CoV-2 infection patterns in the general population22, 23. However, the underlying mechanisms of these noxious, reciprocal allergen-virus effects in asthma are incompletely understood24.
The host response to the RV infection encompasses its RNA recognition by the endosomal toll-like receptor 3 (TLR) 3, TLR7/8 and cytoplasmic RNA helicases: retinoic acid-inducible gene I (RIG-I) and melanoma- differentiation-associated gene 5 (MDA5)25–27, whereas its capsid might interact with the cell surface TLR2 and initiate myeloid differentiation primary response 88 (MyD88)-dependent nuclear factor ’kappa-light-chain- enhancer’ of activated B cells (NF-kB) activation26, 28. RIG-I in its monomeric form binds to the 5’ end of viral RNA and undergoes conformational changes, oligomerizing and exposing the caspase activation and recruitment domains (CARD), enabling further CARD-CARD interactions with mitochondrial antiviral signaling protein (MAVS)29. MAVS recruits tumor necrosis factor receptor-associated factor 3 (TRAF3) to activate TRAF family member-associated nuclear factor kappa B activator (TANK)-binding kinase (TBK)-1 and IκB kinase ε (IKKε) complex. TBK1 complex mediates phosphorylation of interferon (IFN) response factors and subsequent induction of type I and type III (I/III) IFNs30–32. Interferons further signal via their respective receptors which leads to the broad expression of interferon-stimulated genes (ISGs)30–33. Epithelial antiviral response should be sufficient to clear RV infection in healthy airways34. We and others demonstrated several alternations in RV- induced type I/III IFN responses and other antiviral mechanisms in asthma35–39. Despite increased understanding of IFN signaling, no preventive methods or treatments targeting these pathways are available for virus-induced asthma exacerbations, suggesting greater complexity than previously anticipated40–42.
RV infection also leads to the release of proinflammatory cytokines, chemokines and growth factors via activation of MyD88/NF-kB pathway26, 27, 43. This response is likely beneficial in the early infection phase and necessary for effective viral clearance. However, when it is sustained and excessive, it might lead to tissue damage and unresolved inflammation44. HDM also induces expression of proinflammatory mediators in the airway epithelium24 via activation of MyD88/NF-kB and other transcriptional pathways, leading to an increase in expression of proinflammatory proteins45–47. Expression and release of mature IL-1β needs to be tightly regulated, by the transcriptional activation of pro-IL-1β, called priming, followed by an activation of supramolecular complexes called inflammasomes and release of mature, active forms of IL-1β and/or proinflammatory cell death called pyroptosis48. Inflammasomes are composed of at least a sensor protein and caspase-1 and often the adaptor protein apoptosis-associated speck-like protein containing a CARD (ASC)49. ASC is recruited by receptors not containing CARD domains, like NLR family pyrin domain containing 3 (NLRP3), to start oligomerization and further recruitment and activation of pro-caspase 1, via CARD-CARD interactions49. In the case of RIG-I (containing two CARD domains), ASC is recruited to couple and enhance caspase-1 recruitment and mature IL-1β50 production, whereas MAVS interaction with CARD9 is responsible for NF-kB activation and transcriptional priming of pro-IL-1β50. Activation of NLRP3 as well as RIG-I inflammasome has been demonstrated in the context of some respiratory RNA viruses, including RV51–53, influenza A (IAV)52, 54, 55, SARS-CoV-156, 57 and most recently SARS-CoV-258, 59 with major differences depending on the hematopoetic50, 51 or epithelial cellular origin52 and their differentiation state53. It remains unclear whether activation of airway epithelial inflammasomes is necessary to clear infection or in contrast, whether it initiates mucosal hyperinflammation delaying virus clearance60, 61. Additionally, an involvement of NLRP3 priming and/or inflammasome activation in HDM-models of asthma and in severe asthma in humans has been demonstrated, however data are conflicting and remain poorly understood46, 62–64. Finally, little is known about the airway epithelial response in health or during the preexisting HDM-induced inflammation in asthma and combined infection with RV and SARS-CoV-2.
Therefore, in the current study, we investigated the effects of HDM, RV and SARS-CoV-2 on differentiated primary human bronchial epithelium in vitro and in experimental in vivo RV infection in healthy subjects and in patients with asthma. We identified that RV infection and replication activated the RIG-I/ASC, but not the NLRP3, inflammasome, which was further augmented in the presence of HDM. Activation of the RIG-I inflammasome was stronger in patients with asthma and was partially responsible for disturbed RIG-I-dependent IFN signaling, enhanced and sustained inflammation and infection in those patients, especially upon HDM pre-exposure. Finally, we observed that preexisting RV infection and induction of IFNs in the epithelium of patients with asthma decreased SARS-CoV-2 infection, but the presence of both viruses and HDM further increased RIG-I inflammasome activation and release of proinflammatory mediators.
2. Results
2.1 House dust mite increased rhinovirus-induced inflammasome activation in fully differentiated human bronchial epithelium
First, we aimed to investigate mechanisms of the release of mature IL-1β in differentiated primary human bronchial epithelial cells (HBECs) upon rhinovirus A16 (RV-A16) infection with and without exposure to HDM in a dose and time-dependent manner (Fig. 1A-C, S1A-D). RV-A16 infection at a multiplicity of infection (MOI) 0.1, but not UV-inactivated RV-A16 (UV-RV-A16) or HDM alone, significantly induced secretion of mature IL- 1β in patients with asthma 24h after infection (Fig. 1B-C, S1C-D). This increase was accompanied by a higher virus replication at this time point (Fig. S1E). Interestingly, RV-A16-induced IL-1β secretion was even further increased by HDM in patients with asthma or significantly induced in controls (Fig. 1B-C, S1B). Still, this effect was significantly more pronounced in patients with asthma (Fig. 1B-C). Next, we demonstrated that mature IL-1β came from inflammasome activation and formation of ASC specks in HBECs from control individuals and patients with asthma infected with RV-A16 (Fig. 1D, S1F). Complementary to IL-1β secretion, HDM prestimulation increased RV-A16-induced ASC specks formation in HBECs from controls and patients with asthma (Fig. S1G). We did not observe ASC-speck formation after HDM or UV-RV-A16 alone. Once we demonstrated that RV-A16-induced inflammasome activation is enhanced in patients with asthma, we also looked at the baseline status of IL-1β expression and the efficiency of the priming step in both groups. We observed higher expression of pro-IL-1β protein in HBECs from patients with asthma at baseline in vitro (Fig. 1C, S1H), as well as a moderate upregulation of IL-1β concentrations in bronchoalveolar lavage (BAL) fluid in vivo (Fig. S1I). RV-A16, UV-RV-A16 and HDM+RV-A16 stimulation further increased expression of pro-IL-1β mRNA and protein, especially in epithelium in asthma (Fig. 1C, E, F, S1H), suggesting combined effects of HDM and RV-A16-replication independent and dependent mechanisms on pro-IL-1β expression. HDM stimulation, RV-A16 infection alone, or combined with HDM stimulation did not affect protein expression of ASC or pro- caspase-1 (Fig. 1C and S1 J, K). HDM alone increased pro-IL-1β protein expression and release of non-cleaved, pro-IL-1β in control individuals and patients with asthma (Fig. 1C, F). Next, we observed significant inhibition of HDM+RV-A16-induced inflammasome activation upon caspase-1 inhibitor (YVAD) treatment (Fig. 1G), while it did not affect either inflammasome priming (Fig. 1H), nor RV-A16 infection (Fig. S1L). Finally, to investigate whether active RV-A16 infection is necessary for inflammasome activation in HBECs, we blocked RV-A16 entry to the cells using monoclonal antibodies, blocking its receptor ICAM-1. Indeed, significantly diminished RV- A16 infection (Fig. S1M) led to decreased mature IL-1β secretion (Fig. 1I), with no effect on the priming step (Fig. 1J). In summary, we demonstrated here that RV-A16 infection led to the formation and activation of inflammasome in the differentiated primary HBECs, which was enhanced by the exposure to HDM, especially in patients with asthma, who already had higher pro-IL-1β expression at baseline.
2.2 Rhinovirus infection activated the RIG-I ASC inflammasome, but not NLRP3 or MDA5 inflammasomes in human airway epithelium
Next, we investigated which of the pattern recognition receptors (PRR) expressed in human epithelium is a sensor and activator of inflammasome assembly. We found extremely low expression of NLRP3 mRNA in fully differentiated, mature bronchial epithelial cells from patients with asthma and control individuals at baseline (Fig. S2A) or after RV-A16 infection (Fig. S2B). In line with that, we did not detect the expression of NLRP3 protein in differentiated HBECs from controls and patients with asthma at baseline, after RV-A16 infection, HDM stimulation or both (Fig. 2A, B). Lastly, using a specific NLRP3 inflammasome inhibitor (MCC950), we demonstrated, that it did not affect the secretion of mature IL-1β in epithelium from patients with asthma or control individuals upon RV-A16 infection with or without HDM (Fig. 2C). Therefore, since we found that NLRP3 priming or inflammasome activation is not involved in the epithelial response to infection with RV-A16 or HDM stimulation, we analyzed another PRR – RIG-I (DDX58), which was shown previously to form inflammasome in response to RNA viruses either in hematopoietic cells50, 51 or in epithelium52. In contrast to NLRP3, DDX58 (RIG-I) mRNA was expressed in HBECs from both controls and patients with asthma at baseline (Fig. S2C). We observed further increases in DDX58 (RIG-I) mRNA (Fig. 2D, S2D) and RIG-I protein (Fig. 2E, F) expression in HBECs from control individuals and patients with asthma upon RV-A16 infection in the presence or absence of HDM stimulation. Interestingly, RV-A16-induced upregulation of DDX58 (RIG-I) was more enhanced in asthma as compared to control (Fig. 2D). We observed formation of RIG-I speck-like structures upon RV-A16 infection in the confocal microscopy, which further increased in the presence of HDM (Fig. 2G), suggesting RIG-I inflammasome formation. Indeed, coprecipitation of ASC with RIG-I, confirmed RIG-I binding to ASC upon RV-A16 infection and HDM pre-stimulation (Fig. 2H). Notably, MDA5 was not bound to ASC upon RV-A16 infection (Fig. S2E). Since ASC specs and assembled inflammasome complexes are often released from the cells together with mature IL-1β65, we analyzed RIG-I protein expression also in the supernatants of the cells. Accordingly, we found that activation of RIG-I inflammasome led to increased release of RIG-I protein from the epithelial cells (Fig. 2I). Thus, we showed here that infection of human bronchial epithelium with RV-A16 led to RIG-I inflammasome activation, which was strongly enhanced in the epithelium of patients with asthma, especially in the presence of HDM exposure. NLRP3 and MDA5 inflammasomes were not involved in the responses to RV-A16 infection or HDM exposure in differentiated human primary bronchial epithelial cells.
2.3 RIG-I inflammasome- and IL-1β-mediated immune responses are augmented in asthma at an early time point after rhinovirus infection
To evaluate whether the enhanced epithelial RIG-I inflammasome activation in response to RV infection might have further implications on the overall inflammatory responses at the bronchus barrier sites in asthma, we investigated inflammasome- and IL-1β-mediated immune responses. Using targeted proteomics, we found that in addition to IL-1β, also IL-18, IL-1α, tumor necrosis factor (TNF) and TNF-related activation-induced cytokine (TRANCE) were released 24h after RV-A16 infection in HBECs from both control individuals and patients with asthma (Fig. 3A). Monocyte chemoattractant protein (MCP)-1 was significantly upregulated only in asthma, while MCP-3 was significantly upregulated only in controls, but overall, we observed rather similar changes in the expression of the analyzed proinflammatory proteins in HBECs from both patients with asthma and control individuals 24h after infection. Therefore, to evaluate if the enhanced IL-1β release and stronger RV-A16 infection in patients with asthma might influence an even broader range of responses, we analyzed publicly available next-generation sequencing (NGS) raw data of RV-A16-infected HBECs66. An unbiased analysis of pathways and ontologies, revealed upregulated interferon signaling and innate immune responses to RNA viral infections (Fig. 3B, Suppl. Table S1). We also noted a significant enrichment of inflammasome- mediated immune responses, in both – control and asthma samples, but it ranked slightly higher in patients with asthma (Fig. 3B). Accordingly, epithelium from both patients with asthma and control individuals showed increased inflammasome-mediated immune responses after RV-A16 infection (Fig. 3C, D), but the inflammasome-related molecules, cytokines and chemokines, such as CASP1 (caspase-1), IL6, NLR family CARD domain containing 5 (NLRC5), CXCL1 and others were significantly more upregulated in asthma (Fig. 3E). In agreement with the NGS results, we also observed that infection with RV-A16 increased expression of NLRC5 and CASP1 (caspase-1), which was significantly more enhanced in patients with asthma than control individuals (Fig. 3F). In summary, we demonstrated here that early after RV-A16 infection, other inflammasome- and IL-1β-mediated immune responses are boosted in both groups, controls and patients with asthma, but in case of several inflammasome-related molecules, cytokines and chemokines this increase is much more pronounced in patients with asthma.
2.4 Sustained bronchial RIG-I inflammasome activation and inflammasome-mediated immune responses in asthma after rhinovirus infection in vivo in humans
Next, we analyzed samples from a previously reported67 controlled, experimental RV-A16 infection study of patients with asthma and control individuals, including BAL fluid, bronchial biopsies and bronchial brushings taken two weeks before infection (d-14, baseline) and 4 days after infection (Fig. 4A). Both groups were seronegative for anti-RV-A16 antibodies prior to infection and only individuals without recent natural respiratory infection underwent the experimental infection, as described previously67. Quite strikingly, we found a strong upregulation of the same transcriptome profiles of inflammasome-mediated immune responses in bronchial brushings from patients with asthma, in a sharp contrast to the downregulation of similar genes in control individuals at the same time (4 days after infection) (Fig. 4B, C, S3A, B). We further validated the expression of IL-1β, caspase-1, and RIG-I proteins in bronchial biopsies from the same patients. In line with the results demonstrated above, we found higher expression of IL-1β in the epithelial area of the bronchial biopsies of patients with asthma at baseline as compared to healthy controls (Fig. 4D). In accordance with its gene expression in bronchial brushings, in healthy controls, we noted a decrease in IL-1β expression following RV-A16 infection (Fig. 4D). Additionally, we assessed concentrations of mature IL-1β protein secreted into the bronchoalveolar lavage fluid. In agreement with the epithelial mRNA and protein expression of IL-1β, we found significantly decreased IL-1β protein concentration in BAL fluid from control individuals 4 days after infection, whereas in patients with asthma IL-1β protein concentrations in BAL fluid tended to be increased, though this increase was not statistically significant (Fig. 4E). Likewise, we also noted downregulation of epithelial expression of caspase-1 in bronchial biopsies from healthy controls 4 days after infection, whereas it did not change in patients with asthma (Fig. 4F, S3C). Here we also noticed that in vivo RIG-I protein expression in bronchial biopsies at baseline was higher in healthy controls than in patients with asthma, confirming our above-demonstrated in vitro data (Fig. S3D). After RV-A16 infection, there was a trend towards an increase in RIG-I expression in bronchial biopsies from patients with asthma, and towards a decrease in controls, but those trends were not statistically significant. Next, we investigated the dynamics and relationship of inflammasome-mediated immune responses after in vitro, and in vivo RV-A16 infection in patients with asthma and control subjects. We observed a strong negative correlation between gene expression of inflammasome-mediated immune responses after infection at 24h in vitro and at 4 days in vivo in control individuals (Fig. 4G, left panel). Genes upregulated in vitro at 24h post infection, were downregulated or not changed anymore at 4 days after in vivo infection. In contrast, in patients with asthma there was a strong positive correlation between the expression of similar genes at 24h in vitro and 4 days in vivo (Fig. 4G, right panel) showing that genes already strongly upregulated in asthma at an early time point after infection, still stayed significantly upregulated in vivo 4 days after infection. Altogether, combining the data from in vivo and in vitro approaches in humans, we demonstrated that the RIG-I inflammasome is indeed activated in human bronchial epithelium after in vivo infection with RV-A16. Moreover, our data revealed that RIG-I inflammasome-mediated inflammation in healthy individuals was either being actively suppressed or already resolved 4 days after infection, whereas it was still ongoing in asthma.
2.5 Sustained rhinovirus infection and activation of antiviral responses in bronchial epithelium of patients with asthma in vivo
Since the major function of RIG-I is recognition of RNA viruses30, we also analyzed the status of antiviral genes and proteins involved in in vivo and in vitro responses to RV-A16 infection. In line with inflammasome- mediated immune responses, the majority of genes encoding antiviral pathways were still upregulated 4 days after in vivo RV-A16 infection in patients with asthma while they were either downregulated or not changed in healthy controls (Fig. 5A, B, Supplementary Table S2). These data suggest less effective resolution of RV- A16 infection and delayed clearance of the virus in asthma. Indeed, RV-A16 load in the bronchoalveolar lavage fluid in asthma was around 100-fold higher than in controls and the peak nasal lavage virus load was 25-fold higher in patients with asthma than in healthy controls67, though this difference was not statistically significant (Fig. 5C, D). In vitro at 24h post infection, antiviral gene (Fig. 5E, F, Supplementary Table S2) and protein responses (Fig. 5G) were enhanced in HBECs from both patients with asthma and control individuals. Interestingly, virus load was already slightly higher in epithelium from patients with asthma at this time-point (Fig. S1E). Therefore, as in the case of inflammasome mediated-immune responses, due to different dynamics of infection and viral clearance, we found a significant negative correlation of antiviral gene expression in HBECs from controls 24h after in vitro, and in bronchial brushings 4 days after in vivo RV-A16 infection (Fig. 5H, left panel). In contrast, in asthma, the delayed viral clearance and sustained high viral load resulted in a very strong positive correlation between in vitro upregulation of antiviral genes at an early time point and their ongoing upregulation observed in vivo at 4 days post infection (Fig. 5H, right panel). Thus, we showed here that bronchial epithelium from healthy individuals can efficiently respond to RV-A16 infection which leads to rapid virus clearance and subsequent resolution of antiviral responses. In contrast, in asthma, the lack of resolution of antiviral responses and delayed virus clearance suggest that there is an ongoing process in epithelium, which impairs the effectiveness of RIG-I-induced antiviral mechanisms.
2.6 Activation of the RIG-I inflammasome impaired RIG-I dependent interferon signaling in bronchial epithelium of patients with asthma
Having demonstrated excessive RIG-I inflammasome activation and IL-1β secretion in epithelium of patients with asthma in response to RV-A16, we hypothesized that this could result in persistent, but less efficient anti- RV-A16 response in asthma. Therefore, we further studied whether RIG-I activation of MAVS/TBK1/IKKε and downstream interferon signaling inhibits RIG-I inflammasome activation and conversely if formation of RIG-I inflammasome inhibits interferon signaling. First, we used BX795, a chemical inhibitor of TBK1 and IKKε. As expected, it blocked expression of IFNL2/3 (IFN-λ) in HBECs of patients with asthma (Fig. 6A) and subsequently decreased expression of DDX58 (RIG-I) (Fig. 6B). BX795 treatment also reduced expression of interferon- stimulated chemokines: CXCL10, CXCL11 and CCL3 (Fig. 6C). It also led to a trend to increased RV-A16 infection (Fig. 6D), as well as to significantly augmented inflammasome priming (Fig. 6E) and activation (Fig. 6F). Next, we blocked IL-1β processing by RIG-I inflammasome with the use of the caspase-1 inhibitor YVAD and investigated interferon signaling and RV-A16 infection in bronchial epithelium. Inhibition of RIG-I inflammasome activation and subsequent IL-1β signaling indeed led to trends of increasing expression of IFNB (IFN-β) (Fig. 6G) and DDX58 (RIG-I) mRNA (Fig. 6H) and it further increased the production of CXCL10, CXCL11, CCL3, and CCL4 (Fig. 6I). However, it did not change the infection at the same early timepoint (Fig. S1L). In summary, these data suggest that increased RV-A16-dependent RIG-I inflammasome activation in bronchial epithelium disturbed the effectiveness of RIG-I dependent anti-RV-A16 responses in asthma. Therefore, a timely blockade of the excessive RIG-I inflammasome activation and IL-1β signaling may lead to more efficient viral clearance and lower burden of infection.
2.7 House dust mite impaired interferon responses in rhinovirus-infected bronchial epithelium of patients with asthma
Having demonstrated that HDM increased RV-A16-induced RIG-I inflammasome activation, we continued to explore the effect of HDM pre-exposure on the timing and strength of antiviral responses. We found that HDM pre-treatment decreased RV-A16-induced mRNA expression of IFNB (IFN-β) and DDX58 (RIG-I) only in HBECS from patients with asthma 24h, but not 6h after RV-A16 infection in vitro (Fig. 7A, B, S4A, B). Accordingly, when we analyzed protein expression and enriched biological pathways by targeted proteomics in the same conditions, we observed decreased cell- and IFNs- mediated antiviral responses in HBECs from control individuals and patients with asthma at 24h post infection (Fig. 7C). In addition, HDM in the presence of RV- A16 stimulated the release of epithelial to mesenchymal transition factors such as interleukin 15 receptor subunit alpha (IL-15RA)68, 69, artemin (ARTN)70, tolerance inducing TNF receptor superfamily member 9 (TNFRSF9)71 (also called 4-1BB and CD137), extracellular newly identified RAGE-binding protein (EN-RAGE) alarmin72 in both studied groups. However, only in asthma, complementary with all our previous data, HDM simultaneously increased the activation status and release of proinflammatory and pro-remodeling proteins, such as IL-1α, signaling lymphocytic activation molecule family member 1 (SLAMF1)73, cluster of differentiation (CD) 4074, 75 and TRANCE (RANKL)76, 77 (Fig. 7C, Supplementary Table S3). Importantly, HDM pre-stimulation had a slightly additive effect to the antiviral RIG-I pathway inhibitor (BX795) and further reduced BX795-decreased protein expression of the ISGs: CXCL10, CXCL11, CCL3 and CCL4 (Fig. 7D). It all suggests that pre-exposure to HDM, before RV-A16 infection decreases IFN type I response in a non-specific way in patients with asthma and control individuals. However, only in patients with asthma, HDM exposure contributes further to the enhanced inflammasome-mediated and other immune responses, and associated impairment of the effectiveness of antiviral responses.
2.8 Pre-existing rhinovirus infection attenuated SARS-CoV-2 infection, but augmented epithelial inflammation in asthma
Finally, facing the current pandemic and noting the contradictory results about asthma as a risk factor for COVID-19 in different populations12–16, we investigated if RV-A16-induced RIG-I inflammasome activation and HDM-mediated decrease of IFN responses may affect SARS-CoV-2 infection. We first treated primary HBECs from healthy controls and patients with asthma with or without HDM, next after 24h we infected them with RV-A16, and after a further 24h we infected them with SARS-CoV-2 for 48h (Fig. 8A). We confirmed infection with SARS-CoV-2 by the detection of its nucleocapsid protein N (Fig. 8B) and the increase of SARS-CoV-2 viral RNA (Fig. 8C). In patients with asthma, but not in healthy controls, we observed lower infection with SARS- CoV-2 in samples pre-infected with RV-A16 (Fig. 8C). Individual samples with high RV-A16 virus loads (Fig. 8D) had lower SARS-CoV-2 infection, and vice versa. This effect was diminished upon HDM pre-stimulation in samples from patients with asthma, suggesting that HDM pre-stimulation tended to increase SARS-CoV-2 infection in RV-A16+SARS-CoV-2 infected epithelium in asthma (Fig. 8C). To understand these different patterns of infection, we analyzed the expression of DDX58 (RIG-I), type I/III IFNs and IL1B (IL-1β). In line with our previous data, RV-A16 infection increased the expression of DDX58 (RIG-I) (Fig. 8E, S5A), IFNB (IFN-β) (Fig. 8F, S5B), IFNL1 (IFN-λ1) (Fig. 8G, S5C) and IFIH1 (MDA5) (Fig. 8H, S5D). However, we noted that infection with SARS-CoV-2 alone did not induce expression of DDX58 (RIG-I), IFIH1 (MDA5), IFNB (IFN-β), IFNL1 (IFN-λ1) (Fig. 8 E-H, S5 A-D) or secretion of proinflammatory proteins (Fig. 8 I-J, Suppl. Table S4) at this timepoint. In contrast, HDM prestimulation increased the release of IL-18 in SARS-CoV-2-infected HBECs from patients with asthma (Fig. 8 I-J). Importantly though, we observed that infection with SARS-CoV-2 on top of the infection with RV- A16 resulted in the tendency to increase IL-18 (Fig. 8 I-J, Suppl. Table S4) and IL1B (IL-1β) (Fig. 8K, S5E). Notably, HDM pre-stimulation resulted in an amplified secretion of proinflammatory proteins, such as TRAIL, IL-15, CXCL9, IL-17C and CCL8 only in epithelium of patients with asthma infected with both viruses, whereas these conditions induced IL-18 secretion both in patients with asthma and control individuals (Fig. 8 I-J, Suppl. Table S4). In summary, we found here that pre-existing RV-A16 infection restricted SARS-CoV-2 replication in asthma, but not in controls, which was in line with the sustained type I/III IFNs at this time-point, induced by RV-A16 in asthma. SARS-CoV-2 infection alone did not induce DDX58 (RIG-I), IFIH1 (MDA5), IFNB (IFN-β), IFNL1 (IFN-λ1) at 48h after infection. However, especially in the presence of HDM, we observed enhanced proinflammatory responses in patients with asthma after SARS-CoV-2 co-infection with RV-A16, in spite of the reduced viral load of SARS-CoV-2.
Discussion
In the current study, we investigated mechanisms, interactions and dynamics of inflammasomes, antiviral and proinflammatory responses in the airways of healthy controls and patients with asthma upon combined exposures with house dust mite, rhinovirus and SARS-CoV-2. We analyzed the direct in vivo responses in humans, paired with the experiments in patients’ primary cells. We showed, that recognition of replicating RV in bronchial epithelial cells by RIG-I helicase led to ASC recruitment, ASC oligomerization, activation of caspase- 1, processing and release of mature IL-1β, independently of NLRP3 inflammasome activation. The same mechanisms occurred both in healthy individuals and in patients with asthma, but they were significantly more pronounced in asthma, especially when preceded with HDM exposure. Overactivation of epithelial RIG-I inflammasome in asthma compromised the dynamics of RIG-I-dependent type I/III IFNs and ISG responses, leading to less effective virus clearance, as well as to sustained inflammasome-, and IFN-dependent airway inflammation. In addition to enhancing RV-dependent RIG-I inflammasome activation, HDM also partly inhibited early type I/III IFN responses. Interestingly, this sustained RV-induced antiviral response in HBECS from patients with asthma led to restrained SARS-CoV-2 replication, but augmented inflammasome activation and proinflammatory responses.
There are four main inflammasomes described to date to be involved in innate antiviral immunity against RNA viruses – the NLRP3, RIG-I and in some cases MDA5 and absent in melanoma 2 (AIM2) inflammasomes 50, 55, 78, 79. They are activated by several stimuli involved in the viral infection, such as viral nucleic acids, viroporins, RNA- modulating proteins, reactive oxygen species (ROS) and others61. Here, we found, that in vivo in humans and in fully differentiated primary human bronchial epithelium, infection with RV, a single stranded RNA virus, leads to increased priming of pro-IL-1β in a replication independent and dependent manner, and to assembly of RIG-I/ASC inflammasome, in a replication dependent manner. We did not see involvement of NLRP3 inflammasome upon RV infection or HDM exposure or even significant expression of NLRP3 in airway epithelium at baseline in any of our in vivo, in vitro or data mining approaches66, 80. We also did not see MDA5 forming inflammasome. Other groups observed that infection of human peripheral blood mononuclear cells, human macrophages and mouse bone-marrow derived cells with other single stranded RNA viruses, vesicular stomatitis virus (VSV) and IAV, activates RIG-I/MAVS-dependent pro-IL-1β transcription and RIG-I/ASC- dependent, but NLRP3-independent, inflammasome activation and mature IL-1β and IL-18 production50, 81. In contrast, in undifferentiated, submerged cultures of primary human airway epithelial cells infection with IAV revealed RIG-I and NLRP3-inflammasome-dependent mature IL-1β release52, while infection with RV in a similar model led to activation of NLRP3/NLRC5/ASC complexes and mature IL-1β release53. These differences might come from i) the undifferentiated state of the cells in the previous studies, as in non-differentiated epithelium, lacking ciliated cells, viral infection might engage different pathways than in human airways in vivo82, 83; ii) from different expression of inflammasome components in undifferentiated and differentiated mature epithelium lining human airways, as we showed previously84; as well as iii) from the differences in the virus strains and serotypes85. Importantly, our in vitro and in vivo data consistently showed the same results that RIG-I is engaged as inflammasome in bronchial epithelium upon RV infection, which constitutes an important early time-point event, triggering subsequent airway inflammation. It is also possible, that in vivo in humans both inflammasomes are engaged in different cellular compartments: RIG-I/ASC in airway epithelium and NLRP3 in the infiltrating inflammatory cells in the airways. Indeed, RV infection in mice leads to partly macrophage-derived, NLRP3 inflammasome-dependent airway inflammation51. However, neither depletion of macrophages nor NLRP3 knockout leads to the complete blockade of mature caspase-1 and IL-1β processing upon RV infection, underlining that RV might induce other inflammasomes in airway bronchial epithelium51.
An appropriate balance between activation of RIG-I epithelial inflammasome and subsequent IL-1β/IL-1 receptor (IL1R) signaling with RIG-I-dependent type I/III IFN responses should lead to the limitation of viral replication, efficient virus clearance and timely resolution of airway inflammation52. Indeed, we observed here that in the bronchial epithelium of healthy subjects at early time points during RV infection there was an activation of RIG-I inflammasome and inflammasome-mediated immune responses, together with efficient type I/III IFN and ISG-responses. Importantly all of these responses were actively inhibited or went back to the pre-infection state, already 4 days after in vivo infection. In contrast, in epithelium of patients with asthma, there was enhanced RIG-I inflammasome activation accompanied by augmented inflammasome/IL1R- mediated proinflammatory responses starting early after infection and still non-resolved in vivo 4 days after infection. Overactivation of epithelial RIG-I inflammasome and subsequent increases in mature IL-1β release might be at least partially responsible for the delayed and sustained type I/III IFN/ISG responses. We demonstrated this here by blocking caspase-1 with YVAD which led to an increase in IFN-β (IFNB) and RIG-I (DDX58) mRNA together with IFN-responsive chemokines such as CXCL10, CXCL11, CCL3, and CCL4. Our findings are in line with early observations showing that IL-1β is able to attenuate transcription and translation of type I IFNs and excessive IFNα/β-induced effects via proteasome-dependent mechanisms or by induction of prostaglandin E286–88. Interestingly, we also noted here that this cross-talk between IL-1β and type I IFNs is reciprocal, as blocking phosphorylation of TBK1 and IKKε by their inhibitor BX795 and thus reducing RIG-I- induced type I interferons, significantly increased pro-IL-1β transcription and its processing by RIG-I inflammasome in airway epithelium. This was previously elegantly demonstrated in the case of NLRP3 inflammasome89 and in transcriptional IL-1β regulation in macrophages and other myeloid cells90, but it has not been demonstrated for RIG-I inflammasome or considered to be important at the epithelial barrier.
There is a long-lasting discussion about the underlying origins and mechanisms of the frequent viral infections and exacerbations in patients with asthma, which gained a special weight during the current SARS-CoV-2 pandemic16, 91. An impairment in epithelial type I/III IFN responses has been demonstrated in patients with uncontrolled and/or severe asthma, while it is less or not detectable in mild or well-controlled asthma92–94. More recently, it was suggested that this impairment might result from the unfitting grade (too high or too low), incorrect timing of initiation (delayed)95 or resolution of the antiviral and inflammatory response (sustained) in bronchi of patients with asthma96, 97. Treatment with nebulized IFN-β initiated early after development of common cold symptoms, did not provide improvement in comparison to placebo in the whole studied group including a majority of patients with mild asthma42, however in those with moderate/severe asthma, improvements in symptoms and lung function were observed42, 98. Indeed, adding INF-β in vitro to monocyte-derived macrophages and epithelial cells prior to, but not after infection with IAV drastically reduced the number of nucleoprotein-1 positive cells97. These data combined suggest that type I/III IFNs might be extremely important in the antiviral response in the very early phase of infection, whereas they might not be relevant or could be even detrimental in the latter stages. Interestingly, the phenomenon of the frequent and severe exacerbations appears in both atopic, type-2 asthma, as well as in non-atopic, type-2 low asthma, suggesting that there might be a common mechanistic link to it 44, 99. Here, we showed that RV infection in asthma leads to the simultaneous overactivation of RIG-I inflammasome in their airway epithelium and impaired early RIG-I-dependent type I/III IFNs, which overall led to sustained inflammation, but insufficient viral clearance as compared to healthy controls. In patients with asthma, partially due to the higher pro-IL-1β expression, but lower RIG-I expression at baseline, RV infection led to the strong stimulation of RIG-I transcription and translation, and overactivation of both RIG-I-induced pathways: inflammasome and type I IFN responses. Since RIG-I inflammasome is also induced in healthy individuals in the early phase of RV infection, it might provide a physiological mechanism to fight an infection, as shown in ferrets, where the N protein of H5N1 influenza inhibiting this mechanism, resulted in higher mortality52. However, its overactivation in patients with asthma, especially at the very early stage of the infection, might contribute rather to the damage of the infected ciliated epithelium than to the effective antiviral response. Such damage might occur via pyroptosis100 and lead to the release of active IL-1β and RIG-I/ASC/caspase-1 inflammasome complexes to the extracellular space and subsequent barrier impairment2, 101, recruitment of macrophages and neutrophils101, and as we also showed here, to the inhibition of a timely and functional antiviral response and delayed in vivo viral clearance.
Presence of other airway barrier-damaging and/or activating factors, such as exposure to allergens in addition to the viral infection, worsens the clinical outcomes, leads to a more severe exacerbation, hospitalization or respiratory failure11. Mechanistically, it might be connected with the multiplication of pathways activated in airway epithelium and in cells infiltrating the airways, and/or with the additive effects of different triggers on the same pathway24, 91. HDM activates airway epithelium via, among others, TLR2/4, C-type lectins, and PARs in an allergen-non-specific way to initiate allergen sensitization, but also to perpetuate already developed allergic and probably non-allergic airway inflammation in the absence of the sufficiently developed inhibitory signals45–47. It is known that one of the strongest effects of HDM-specific immunotherapy and anti- immunoglobulin E (IgE) treatment is the reduction of the rate of asthma exacerbations. It might be at least partially connected to the fact that in sensitized individuals, allergen binding to specific IgE, captured by the Fc epsilon receptor (FcεRI) receptor on the surface of plasmacytoid dendritic cells (DCs) may lead to the additional impairment of antiviral responses102. However, it is not well understood if and how HDM-induced signaling in the airway epithelium interferes with the effectiveness of early antiviral response. Here, we found that HDM contributes significantly to RV-induced pathologic responses in human airway epithelium in patients with asthma by i) enhancement of priming and non-mature pro-IL-1β release, ii) overactivation of RIG-I inflammasome and subsequent release of mature caspase-1, mature IL-1β, and RIG-I , iii) inhibition of type I/III IFNs and ISG-responses and iv) activation of extra proinflammatory and pro-remodeling proteins, such as IL-1α, SLAMF173, CD4074, 75 and TRANCE (RANKL)76, 77. In the presence of both HDM and RV, we noted that increased usage of RIG-I protein engaged in inflammasome formation and its subsequent expulsion is paired with HDM-dependent inhibition of RIG-I, type I/III IFNs and several ISGs on the mRNA and protein level. At an early stage of RV infection, it might explain disturbance of the antiviral response dynamics observed in asthma by us and others103. Due to this functional reduction of RIG-I availability, the following type I/III IFN-response is less effective, not able to quickly and efficiently clear the infection and thus it is sustained up to later time points, together with the enhanced inflammasome-related proinflammatory response, as we demonstrated here in vivo and in vitro.
Co-infections with two or more respiratory viruses occur often and likely acts as additional factors increasing airway epithelial damage. Patients with asthma are at a greater risk of developing respiratory failure as shown in the case of H1N1 influenza infection104. Thus, it has been somewhat surprising that so far during the current SARS-CoV-2 pandemic, epidemiological cohorts of COVID-19 patients from different geographical locations have resulted in partially contradictory observations that asthma is (USA, United Kingdom, Australia) or is not (Europe, China) a risk factor for SARS-CoV-2 infection and/or severity of COVID-1915, 105, 106. We and others demonstrated that the expression of angiotensin-converting enzyme 2 (ACE2), the main SARS-CoV-2 receptor and other plausible points of entry, are not changed in patients with asthma, even though different types of airway inflammation or inhaled steroids might modulate their expression107, 108 and as such it seems unlikely that it would be the main reason for observed discrepancies. Different geographical locations might represent variable levels and quality of environmental exposures such as viruses or allergens, which may interfere with the rate of SARS-CoV-2 infections23, 109 and COVID-19 severity. SARS-CoV-2, an enveloped, positive-sense, ssRNA virus110 has been shown to be sensed, depending on the cell type, by MDA5111, RIG-I111 and NLRP3112, 113. However, due to several evasion properties and encoding by non-structural (Nsp) and accessory proteins, such as Nsp1,6,12,13114, various open read frames (ORFs)115, protein M116, protein N117 and others, which antagonize interferon pathways on many levels118 , induction of IFNs by SARS-CoV-2 itself is reduced or delayed114, 119 with the augmented proinflammatory mediator release. Several in vitro and in vivo animal studies revealed that SARS-CoV-2, similarly to SARS-CoV-1, is very sensitive to pre-treatment with type I/III IFNs120–122, which inhibit its replication. Patients with disrupted IFN gene expression and production or patients with autoantibodies against type I IFNs have been shown to be at a greater risk of severe COVID-19. In accordance with these studies, we observed here that infection of epithelium from patients with asthma at the moment of heightened and sustained IFNs type I/III response induced by RV, results in reduced SARS-CoV- 2 viral load, in contrast to the healthy epithelium, where the RV-induced IFN-I/III response has been already actively resolved. However, even if restricted, SARS-CoV-2 infection in combination with RV and HDM led to an increase in activation of epithelial inflammasome and release of higher amounts of IL-18 and other proinflammatory cytokines. In context of timing and possible clinical relevance it may mean that patients with asthma with pre-existing RV-infection might have a slightly restricted SARS-CoV-2 infection at first, but due to excessive inflammasome-related damage and proinflammatory signaling, together with SARS-CoV-2-induced inhibition of type I/III IFNs, they may in fact succumb eventually to more severe COVID-19. Importantly, in the presence of HDM, these potentially adverse effects are further heightened, meaning that HDM reduces RV- induced IFN response, which leads to higher SARS-CoV-2 replication and it enhances RIG-I inflammasome activation, inflammation and tissue damage.
All in all, we showed here in vivo and in vitro that the lack of balance between activation of RIG-I inflammasome and the RIG-I-IFNs-axis in response to a common respiratory virus, is an important driving factor of epithelial damage, lack of viral clearance and sustained airway inflammation in patients with asthma. Timely targeting of this abnormal response by the yet-to-be developed early therapeutics or even prophylactic approaches might provide in the future a beneficial strategy to prevent RV-induced exacerbations of asthma and potentially severe COVID-19.
3. Material and Methods
Reagents
House dust mite (HDM) extract (HDM) (Allergopharma, Reinbek, Germany) and house dust mite extract B (HDM B) (Citeq, York, UK) were diluted in sterile 0.9% saline (NaCl) and stored in -20°C. The concentration of the HDM used for the experiments was calculated according to the total protein content. Detailed description of the HDM extracts including protein, main allergens and endotoxin content is presented in the Supplementary Table S5.
All commercially available antibodies and reagents used in the manuscript are described in the Supplementary Table S6.
Mouse IgG2a monoclonal anti-human ICAM-1 antibody (antibody R6.5) was produced from the hybridoma cells (ATCC HB-9580, mouse hybridoma). The hybridoma cells were expanded in RPMI 1640 medium supplemented with 10% (v/v) IgG depleted fetal bovine serum, 2 mM L-glutamine, 1.25 g/L D-(+)-Glucose, 1 mM Sodium Pyruvate, 10 mM HEPES, 100 U/mL Penicillin-Streptomycin at 37°C in 5% CO2. The IgG in culture media were affinity purified from the cell culture supernatant using a 1 mL HiTrap™ Protein G HP column (GE Healthcare, 29-0485-81). Eluted fractions were immediately neutralized and buffer exchanged into PBS using dialysis. The antibody was then filtered through a 0.22 µm filter and stored at 4°C.
Viruses
Rhinovirus A16 (RV-A16) was purchased from Virapur (San Diego, USA). UV-light inactivated RV-A16 (UV-RV- A16) was used as a control after an exposure to UV-light of the 254 nm at a 2 cm distance for 60 min. Cells were infected with RV-A16 or UV-RV-A16 at the multiplicity of infection (MOI) 0.1, 0.01 and 0.001 as determined by plaque assay. Briefly, HeLa cells were infected with the virus serial dilutions from 10-2 to 10-8 in duplicates. Seven days after infection, cells were fixed with formaldehyde solution and stained with 1% crystal violet in 20% ethanol and dH20. Visible plaques were counted under a microscope.
SARS-CoV-2 viral genome was generated from the synthetic DNA fragments produced by GenScript (Piscataway, USA) using the in-yeast transformation-associated recombination (TAR) cloning method, as previously described123. In-vitro transcription was performed for the linearized Yeast artificial chromosome (YAC), containing the cDNA of the SARS-CoV-2 genome, as well as a PCR amplified SARS-CoV-2 N gene using theT7 RiboMAX Large Scale RNA production system (Promega, Madison, USA) with m7G(5ʹ)ppp(5ʹ)G cap provided as described previously124. Transcribed, capped mRNA was subsequently electroporated into baby hamster kidney cells (BHK-21) expressing SARS-CoV N protein. Co-culture of electroporated BHK-21 cells with susceptible Vero E6 cells produced passage 0 of SARS-CoV-2 virus. Passage 0 was used to infect Vero E6 cells to generate passage 1 working stocks, which were used for all experiments. Titers were determined using standard plaque assay, as described previously123.
Patients groups and experimental in vivo RV-A16 infection in humans
Experimental in vivo rhinovirus infection in 11 control individuals and 28 patients with asthma was performed as reported previously67. Briefly, non-smoking, non-atopic control individuals, and non-smoking mild/moderate patients with asthma without any recent viral illness and without serum neutralizing antibodies towards RV-A16, who passed inclusion criteria, underwent infection on day 0 with RV-A16 at the dose of 100 TCID50. Bronchial brushings, bronchial biopsies and bronchoalveolar lavage (BAL) fluid were collected around 2 weeks before and at 4 days after RV-A16 infection. Additionally, nasal lavage (NL) samples at the peak of RV-A16 infection were collected to assess RV-A16 infection rates. Only subjects who had sufficient remaining samples to be analyzed in this study and/or subjects who had successful infection in the lungs, as assessed by viral RNA copies by qPCR, were included in the BAL, NL, and biopsies analyses (n=9 healthy control, n=19 patients with asthma), and bronchial brushing microarray analysis (n=7 healthy controls, n=17 patients with asthma). The study received ethical approval from the St. Mary’s Hospital Research Ethics Committee (09/H0712/59). All participants gave written, informed consent. The clinical characteristics of the 9 control and 19 asthma study participants who had sufficient remaining samples to be analyzed in this study is presented in Supplementary Table S7.
Control individuals and patients with asthma were enrolled in the ALL-MED Medical Research Institute, Wroclaw, Poland; the Pulmonary Division, University Hospital of Zurich, Switzerland (cohort SIBRO)125, or at the University Hospital, Jagiellonian University Medical College, Cracow, Poland (cohort A)126, as described previously. Briefly, bronchoscopy with epithelial cells brushings and BAL fluid collection was performed. The study was granted ethical permission from Switzerland and Poland (KEK-ZH-Nr. 20212-0043 – Kantonale Ethik- Kommission Zürich; KB-70/2013 and KB-567/2014 – Bioethical Committee, Wroclaw Medical University) or the Jagiellonian University Bioethics Committee (KBET/68/B/2008 and KBET/209/B/2011). Asthma diagnosis and severity were assessed according to the GINA guidelines3. All participants gave written, informed consent. Clinical characteristics of the study participants is presented in the Supplementary Table S7.
Primary Human Bronchial Epithelial cells (HBECs) were obtained from the above-listed cohorts or from the doctor-diagnosed asthma and control individuals from two independent commercial sources: Lonza (Basel, Switzerland), and Epithelix (Plan-les-Ouates, Switzerland). Characteristics of the HBECs used in the manuscript are presented in Supplementary Table S8.
Air-liquid interface (ALI) cultures of bronchial epithelium from control individuals and patients with asthma
HBECs from the control subjects and patients with asthma were cultured and differentiated in the air-liquid interface (ALI) conditions as described previously, with minor alterations of the previous protocol126. Briefly, cells from passage 2 were grown in 20 mL of bronchial epithelial basal medium (Lonza, Basel, Switzerland) supplemented with the SingleQuot Kit (Lonza, Basel, Switzerland) placed in 150cm2 T-flask in humidified incubator at 37°C with 5% CO2 for maximum 10 days, or until 80%-90% confluency. Next, cells were trypsinized (ThermoFisher Scientific, Waltham, USA) and seeded at a density of 1.5x105 cells/well on the 6.5-mm-diameter polyester membranes with the 0.4 μm pore size and growth area of 0.33 cm2 (Costar, Corning, NY, USA; Oxyphen, Wetzikon, Switzerland) in 24-well culture plates. Bronchial Epithelial Growth Medium (BEGM) (Lonza, Basel, Switzerland) supplemented with the SingleQuot kit (Lonza, Basel, Switzerland), with an exception of the retinoic acid (ATRA, Sigma-Aldrich, St. Louis, USA) and triiodothyronine, was mixed in the 1:1 ratio with the Dulbecco modified Eagle medium (DMEM, Gibco, Thermofisher Scientific, Waltham, USA). Fresh all-trans ATRA (Sigma-Aldrich, Merck, Kenilworth, USA) was supplemented at a concentration of 15 ng/mL. Cells were grown submerged for 3-5 days in the apical medium and were in contact with the basolateral medium. After they obtained a full confluence, the apical medium was removed and cells were kept in the air- liquid interface (ALI) cultures for at least 21 days. BEGM/DMEM/ATRA medium was maintained only basolaterally to differentiate the HBECs. During the cell culture process, medium was exchanged every 2-3 days and, periodically, excess of produced mucus was removed from the wells. All experiments were performed on the fully differentiated HBECs from the same passage, between 21 and 28 days of ALI culture (Fig. S1A).
House dust mite stimulation and rhinovirus A16 in vitro infection model in the primary HBECs
House dust mite (HDM) stimulation, followed by rhinovirus A16 (RV-A16) infection experiments were performed in the OptiMEM medium (LifeTechnologies, ThermoFisher Scientific, Waltham, USA). ALI- differentiated HBECs from control individuals and patients with asthma were treated apically with the HDM extract (Allergopharma, Reinbek, Germany) at a dose of 200 μg/mL of the total protein in 200 μl OptiMEM on the apical side, and 600 μl of clear OptiMEM on the basolateral side (Fig. S1A), in the humidified incubator at 37°C with 5% CO2. After 24h of HDM stimulation cells, were apically infected with RV-A16 at the MOI of 0.1 or as otherwise specified, or stimulated with UV-RV-A16 at the same MOI, and cultured in the humidified incubator at 34.5°C with 5% CO2 for the next 24h (Fig. S1A). Next, cell supernatants (apical and basolateral), RNA, and protein cellular lysates were collected and stored in -80°C. Some cells were fixed with 4% PFA (Fluka/Sigma Aldrich Buch, Switzerland) and were stored wet at 4°C for 1-2 weeks before the subsequent confocal analyses. All doses and time-points used for the final experiments were based on the preliminary dose-dependent and time-course experiments. Briefly, two different HDM extracts: main HDM extract used in the manuscript (Allergopharma, Reinbek, Germany) at the dose of 200 μg/mL, and HDM extract B (Citeq, York, UK) at the dose of 200 and 100 μg/ml were investigated (Fig. S1B, S4A-B). RV-A16-infection of HBECs from patients with asthma was performed in 6h and 24h time-points (Fig. S1C). Lastly, RV-A16 infection at the MOI 0.001, MOI 0.01, and MOI 0.1 was investigated (Fig. S1D). Based on the secretion of the mature IL-1β, HDM extract from the Allergopharma, Reinbek, Germany at the dose of 200 μg/mL, RV-A16 MOI 0.1 and the 24h time-point were chosen, and are presented through the manuscript, if not mentioned differently. For the experiments with inhibitors, 40 μM of the caspase-1 inhibitor: ac-YVAD-cmk (Acetyl-tyrosine-valine-alanine- aspartate-chloromethyl ketone, Invivogen, San Diego, USA), 1μM of the IKKε/TBK1 inhibitor: BX795 (N-[3-[[5- iodo-4-[[3-[(2-thienylcarbonyl)amino]propyl]amino]-2-pyrimidinyl]amino]phenyl]-1-Pyrrolidinecarboxamide hydrochloride, Sigma Aldrich, Merck, Kenilworth, USA), or 1μM of the NLRP3 inflammasome inhibitor: MCC950 (C20H23N2NaO5S, Avistron, Bude, UK) or appropriate vehicle controls were used apically and basolaterally, 24h prior RV-A16 infection. To block ICAM-1, a receptor responsible for RV-A16 infection of HBECs, anti-ICAM-1 antibodies were added to the apical and basolateral compartment, 3h prior RV-A16 infection at the dose of 10ug/mL (Fig. S1A).
Rhinovirus and severe acute respiratory syndrome coronavirus-2 (SARS-COV-2) in vitro co-infection model in primary HBECs from control individuals and patients with asthma
The ALI-differentiated MucilAir cultures (Epithelix, Plan-les-Ouates, Switzerland) from primary human bronchial epithelium obtained from 4 control individuals and 5 patients with asthma (Supplementary Table S8) were cultured for 7 days in the MucilAir Medium (Epithelix, Plan-les-Ouates, Switzerland) in ALI conditions, with basolateral medium changed every other day. At the day of the experiment, performed in a biosafety level 3 (BSL3) laboratory, cells were washed with warm PBS to remove an excess of mucus. The experiment was performed in the OptiMEM medium (LifeTechnologies, ThermoFisher Scientific, Waltham, USA) in the volume of 250 μl on the apical, and 600 μl on the basolateral side. Through the whole experiment cells were kept in the humidified incubator at 37°C with 5% CO2. First, HBECs were stimulated apically with 200 μg/mL of protein content of HDM extract or vehicle. 24h after HDM stimulation, cells were apically infected with/without RV-A16 at the MOI of 0.1. After next 24h, HBECs were apically infected with/without SARS-CoV- 2 at the MOI of 0.1. Finally, 48h after SAR-CoV-2 infection experiment was harvested (Fig. 8A). In order to inactivate SARS-CoV-2, all collected supernatants were treated with 65°C for 30 min. Cells were fixed in 4% PFA for at least 20 min. Inactivated supernatants were frozen in -80°C until further analyses. For RNA analyses, insert with the fixed cells were preserved in RNAlater (Qiagen, Hilden, Germany), left overnight in 4°C, and stored in -20°C in the new, dry tube. In order to perform confocal staining, inserts with the fixed cells were snap frozen in the Clear Frozen Section Compound (FSC22, Leica, Wetzlar, Germany).
THP-1 cell culture
THP-1-XBlue cells (Invivogen, San Diego, USA) were defrosted in 32 mL of RPMI-1640 medium (Sigma-Aldrich, St. Louis, USA) supplemented with the Penicillin/Streptomycin/Kanamycin, MEM vitamins, Na-Pyruvate/MEM Non-essential Amino Acid Solution and heat-inactivated FCS (cRPMI medium) in the 75cm2 T-flask, and cultured for 1 day in the humidified incubator at 37°C with 5% CO2. In the following day, cells were counted, checked for viability (98%) and transferred to the 12-well cell cultures plate (0.5 mio cells/well in 1 mL of cRPMI medium). Next day, cells were stimulated with LPS (100 ng/mL, Invivogen, San Diego, USA) or vehicle for 4h followed by 2 mM ATP or vehicle (Invivogen, San Diego, USA) for 20 min. Cytospins (250xg, 3 mins, Shandon Cytospin 2, Marshall Scientific, Hampton, USA) were prepared, and cells were immediately fixed with 4% PFA (Fluka/Sigma Aldrich, Buchs, Switzerland), and stored in wet chamber before the confocal staining.
Immunoassays
ELISA and MSD multiplex
Secreted IL-1β in majority of in vitro experiments was measured with the ELISA kit (R&D Systems, McKinley Place, USA), according to the manufacturers instruction. Detection limit for the kit is 3.91 pg/mL. IL-1β in BAL fluid from control individuals and patients with asthma was measured using mesoscale discovery platform (MSD, Kenilworth, USA) kits as described previously125. Additionally, BAL samples from controls and patients with asthma before and after experimental intranasal RV-A16 infection in vivo were analyzed with V-PLEX human IL-1β Kit (MSD, Kenilworth, USA), according to the manufacturer’s instructions. V-PLEX IL-1β is presented as log10 arbitrary units (au) corresponding to percentage of total protein concentration measured by BCA (ThermoFisher Scientific, Waltham, USA) and multiplied by factor 1000000.
Proximity extension assay (PEA) targeted proteomics
Protein expression in the apical compartments of the HBECs were measured using the proximity extension assay targeted proteomics technology (Olink, Stockholm, Sweden). Targeted 96-proteins Inflammation (control and asthma), Immune Response and Immuno-Oncology (asthma) panels were used according to the manufacturer’s instructions with the suggested adaptations to the cell cultures conditions. Apical compartments from RV-A16+SARS-CoV-2 model were analysed with the human Target 48 Cytokine Panel (Olink, Uppsala, Sweden). Final results for 96-plex assay are reported as Normalized Protein eXpression (NPX), an arbitrary unit in log2-scale. Results from the Target 48 Cytokine Panel are in pg/mL.
Western Blotting
Western Blotting experiments from the cell lysates and the apical supernatants were performed as previously described126, 127. Briefly, cells were lysed in RIPA Lysis and Extraction buffer (ThermoFisher Scientific, Waltham, USA) supplemented with the protease inhibitor (Roche, Merck, ThermoFisher Scientific, Waltham, USA) for 15 minutes on ice, centrifuged (15 min, full speed, 4°C), and debris-free protein lysates were frozen in -80°C for further analyses. Protein concentration was assessed with the BCA kit (ThermoFisher Scientific, Waltham, USA), according to the manufacturer’s instructions. Protein from the apical supernatants was precipitated with 1 volume of methanol (Fisher Scientific, Reinach, Switzerland) and ¼ volume of chloroform (Merck Millipore, Burlington, USA) as described previously127. Equal amounts of cell lysate proteins (10-20μg) were loaded on the 4-20% Mini-PROTEAN TGX Gel (Bio-Rad, Hercules, USA) or 4-20% SuperPAGE gel (GenScript, Leiden, Netherlands) in Tris/Glycine/SDS buffer (Bio-Rad Lab, Hercules, USA) or MOPS buffer (GenScript, Leiden, Netherlands). After electrophoresis, the proteins were transferred to the nitrocellulose membranes (Bio-Rad, Hercules, USA or Advansta, San Jose USA) using the Trans-Blot Turbo Blotting System (Bio-Rad, Hercules, USA) or eBlot L1 Protein Transfer System (GenScript, Leiden, Nederlands). The membranes were blocked with 5% nonfat milk in 0.1% Tween20 PBS (PBST) for 1h, washed with 10x PBST, and incubated with the primary antibodies for overnight in 4°C. Dilutions of primary antibodies used for the analyses of the cell lysates: 1:100 IL-1β (R&D Systems, McKinley Place, USA), 1:200 RIG-I (Santa Cruz Biotechnology, Dallas, USA), 1:1000 ASC (Santa Cruz Biotechnology, Dallas, USA), 1:1000 caspase-1 (Cell Signaling, Danvers, USA), and 1:1000 NLRP3 (Adipogen, San Diego, USA). The membranes were subsequently washed in 10xPBST and incubated with an appropriate secondary antibody conjugated with the horseradish peroxidase (HRP) (1:10,000 dilution) (Jackson ImmunoResearch, West Grove, USA; Santa Cruz Biotechnology, Dallas, USA) for 1h at room temperature. β-actin expression was determined with HRP-conjugated antibodies in 1:25,000 dilution (Abcam, Cambridge, UK). Protein samples precipitated from the apical supernatants were analyzed with the use of goat anti-IL-1β antibodies (1:1000, R&D Systems, Minneapolis, USA), rabbit anti-caspase-1 (1:1000, Cell Signaling Technology, Danvers, USA), and HRP conjugated anti-goat (1:10,000, Santa Cruz Biotechnology, Santa Cruz, USA), and anti-rabbit (1:10,000, Jackson ImmunoResearch, West Grove, USA) antibodies. After washing with 10xPBST, the blots were developed with the SuperSignal West Femto Kit (ThermoFisher Scientific, Waltham USA) or WesternBright Quantum/Sirius HRP substrate (Advansta, San Jose, USA) and visualized on the Luminescent Image Analyzed LAS-1000 (Fujifilm, Tokyo, Japan) or the Fusion FX (Vilber, Collegien, France). To assess more than one protein, the membranes were stripped with the Restore PLUS Western Blot Stripping Buffer (ThermoFisher Scientific, Waltham, USA). Quantification of the protein expression was performed in Fiji Software. 128 Briefly, an area of the peak of the protein of interest was measured in triplicates, and average was used to calculate the ratio between expression of the protein of interest and β-actin (protein/β-actin). Protein/β-actin ratio from HBECs from control individuals and patients with asthma upon HDM stimulation, RV-A16 infection, or both was further normalized to the vehicle control condition from control individuals, and log transformed with Y=log(Y) function.
Co-immunoprecipitation
For co-immunoprecipitation cells were lysed with the Lysing Buffer (1μM DTT + 10% Triton X100 in ddH20 supplemented with the protease inhibitor (Roche, Merck, ThermoFisher Scientific, Waltham, USA) for 15 minutes on ice, centrifuged (15 min, full speed, 4°C), and the debris-free protein lysates were frozen in -80°C for further analyses. Protein concentration was assessed with the BCA kit (ThermoFisher Scientific, Waltham, USA) according to the manufacturer’s recommendation. 100 μg of the pooled protein was pre-cleared with 100 μl of Protein A beads (Bio-Rad, Hercules, USA), magnetized, and incubated with 10 μg of anti-ASC antibodies (Santa Cruz Biotechnology, Santa Cruz, USA) overnight at 4°C, followed by ASC immunoprecipitation with 100 μl of Protein A beads (Bio-Rad, Hercules, USA) for 2h in RT. Co-IP samples and input (protein not bound to the beads) were collected and together with the cell lysates were further analyzed with the Western Blot protocol with the use of mouse anti-human RIG-I antibodies (1:200, Santa Cruz Biotechnology, Santa Cruz, USA) or rabbit anti-human MDA5 antibodies (1:1000, Abcam, Cambridge, UK) and HRP conjugated anti-mouse or anti-rabbit antibodies (1:10,000, Jackson ImmunoResearch, West Grove, USA).
mRNA isolation and RT-PCR
HDM and RV-A16 model
Cells were lysed on ice with the RLT buffer (Qiagen, Hilden, Germany) supplemented with β-mercaptoetanol (Sigma-Aldrich, St. Louis, USA), and stored at -80°C until further analyses. mRNA isolation was performed with the RNeasy Micro Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Quality and quantity of isolated RNA was assessed by the Nanodrop 2000 (ThermoFisher Scientific, Waltham, USA). Reverse transcription was performed with RevertAid RT kit (ThermoFisher Scientific, Waltham, USA) with random hexamers, according to the manufacturer’s recommendations. Gene expression (5-10 ng of cDNA/well) was assessed by RT-PCR using the SYBR Green/ROX qPCR Master Mix (ThermoFisher Scientific, Waltham, USA), performed on the QuantStudio 7 Flex Real-Time PCR System (ThermoFisher Scientific, Waltham, USA). The sequences of used primers (5μM) are summarized in the Supplementary Table S9. Gene expression was normalized to the housekeeping gene, elongation factor 1α (EEF1A), and presented as a relative quantification calculated with the ΔΔCt formula, as described previously46. Depending on the analyses, data were calibrated according to the vehicle condition from HBECs from control individuals, or vehicle condition calculated separately for control individuals and patients with asthma. Data are presented as 2-ΔΔCt values, or percentage change normalized to the specified condition.
HDM, RV-A16 and SARS-CoV-2 model
Samples preserved in the RNAlater, as described above, were immersed in the increasing concentrations of ethanol (30% up to 100%, increasing every 10%). After this initial step, RNA was isolated with use of RecoverAll kit (Invitrogen, Waltham, USA) according to the manufacturer’s protocol. Isolated mRNA was concentrated with the use of SpeedVac (DNA Speed Vac, DNA110, Savant, Hayanis, USA) for 1h in the ambient temperature. Sample concentration did not affect its quality, as measured with use of NanoDrop OneC (ThermoFisher Scientific, Waltham, USA). Reverse transcription was performed with use of the SuperScript IV VILO Master Mix (Thermofisher Scientific, Waltham USA), according to the manufacturer’s recommendations. Gene expression (5ng of cDNA/well) was assessed by RT-PCR using i) SYBR Green PCR Master Mix (ThermoFisher Scientific, Waltham, USA) for DDX58, IFNB, IFNL1, IL1B, MDA5, and ii) TaqMan assays for RV-A16 and SARS- CoV-2 Protein N, Protein S, ORF1AB detection (ThermoFisher Scientific, Waltham, USA) and was performed on the QuantStudio 7 Flex Real-Time PCR System (ThermoFisher Scientific, Waltham, USA). The sequences of used primers are summarized in the Supplementary Table S9. Gene expression was normalized to the elongation factor 1α (EEF1A) (SYBR Green analyses), or RNase P (TaqMan assays), and presented as a relative quantification calculated with -ΔΔCt formula, compared to the vehicle condition separately for controls and patients with asthma, as described previously46. SARS-CoV-2 viral RNA in presented as 2-ΔΔCt values averaged from the expression of N protein, S protein and ORF1AB in each condition.
Rhinovirus infection quantification in the nasal and bronchoalveolar lavage fluid in humans
Rhinovirus infection in control individuals and patients with asthma after experimental RV-A16 infection in vivo was performed in the nasal lavages (peak of infection) and BAL fluid (4 days post infection) with use of qPCR, as previously described. 67 Results are presented as log10 of viral RNA copies in 1 mL.
Confocal microscopy
Cells were fixed on the inserts with 4% paraformaldehyde (Fluka/Sigma Aldrich, Buch, Switzerland) for 10 minutes, permeabilized with detergent (PBS + 0.1% TritonX100 + 0.02% SDS) for 5 minutes and blocked with 10% goat serum (Dako, Agilent, Santa Clara, USA) in 1% BSA/PBS for 60 min at room temperature (RT). All antibodies were diluted in 4% goat serum + 1% BSA/PBS, and cells were stained from apical and basolateral sides with 100 μl of antibodies working solution. Samples were stained for ASC (2μg/mL, mouse anti-ASC, Santa Cruz Biotechnology, Santa Cruz, USA), IL-1β (10μg/mL, mouse anti- IL-1β, Abcam, Cambridge, UK), and RIG-I (2μg/mL, mouse anti-RIG-I, Santa Cruz Biotechnology, Santa Cruz, USA) for 60 minutes at RT. Proper mouse isotype controls, in the corresponding concentrations were used to control for unspecific binding. (Dako, Agilent, Santa Clara, USA). Subsequently, samples were incubated with the goat anti-mouse Alexa Fluor 488 (for ASC), and the goat anti-mouse Alexa Fluor 546 (for IL-1β and RIG-I) secondary antibodies at the concentrations of 1:2000 (Invitrogen, Waltham, USA) for 60 minutes at RT. Samples were mounted in the ProLong Gold mounting medium containing DAPI (Life Technologies, Carlsbad, USA) according to the manufacturer’s instructions, analyzed with a Zeiss LSM780 confocal microscope (Zeiss, Oberkochen, Germany) and Zen Software (Zeiss, Oberkochen, Germany). All pictures were taken at the 40x magnification and are presented as maximal projection (orthogonal projection) from z-stacks (3-22 for ASC, 4 for IL-1β and RIG-I). ¼ of the original picture is shown on the figures, with appropriate scale bar and further magnification of the area of interest. Additionally, to quantify ASC speck formation specks from 3-5 pictures per condition were counted in duplicates and averaged.
Differentiated HBECs upon HDM+ RV-A16 stimulation and THP-1 cells stimulated with LPS+ATP were used for NLRP3 and Occludin visualization, whereas HBECs from RV-A16+SARS-CoV-2 model were used for ACE2 and SARS-CoV-2 Protein N staining. Following fixation in 4% (wt/vol) PFA in PBS (Fluka, Fluka/Sigma Aldrich, Buchs, Switzerland), THP-1 cytospins were lined with the PAP pen (Sigma Aldrich, St. Louis, USA), HBECs were stained on the insert, whereas ¼ of the ALI insert from RV-A16+SARS-CoV-2 model were prepared for cryosections by freezing in Clear Frozen Section Compound (FSC22, Leica, Wetzlar, Germany), cut at 6 µm in a cryostat (LEICA CM3050S, Leica Microsystems, Wetzlar, Germany) and mounted on SuperFrost PlusTM glass slides (Menzel, ThermoFisher, Waltham, USA). Samples were incubated in the blocking solution (10% normal goat serum, 1% bovine serum albumin and 0.2% TritonX-100 in PBS) (Dako, Agilent, Santa Clara, USA) for 1h at RT. Primary antibodies for NLRP3 (5 µg/mL, mouse anti-NLRP3, Adipogen, San Diego, USA), Occludin (2.5 µg/mL, mouse anti-Occludin, ThermoFisher, Waltham, USA), ACE-2 (2 µg/mL, rabbit anti-ACE2, Abcam, Cambridge, UK), and SARS-CoV-2 Protein N (1 µg/mL, mouse anti-Protein N, ThermoFisher, Waltham, USA) diluted in blocking solution (1:1 in PBS) and incubated at 4°C overnight. Proper isotype controls, in the concentrations corresponding to the antibodies, were used (Dako, Agilent, Santa Clara, USA). Following three washing steps in 0.05% Tween20 in PBS, secondary antibodies with DAPI (1 µg/mL, Sigma Aldrich, St. Louis, USA) in diluted blocking solution (1:1 in PBS) were applied for 2h at RT in the dark. Sections have been washed three times in 0.05%Tween20 in PBS before mounting with Fluoromount (Sigma Aldrich, St. Louis, USA). Image acquisition was performed with Zeiss LSM780 (Zeiss, Oberkochen, Germany), by using 40x objective and ZEN software (Zeiss, Oberkochen, Germany). ImageJ/Fiji128 was used for tale scan stitching and image analysis.
Bronchial biopsies were collected from the control individuals and patients with asthma at baseline and 4 days after in vivo RV-A16 infection. Biopsies were fixed and embedded in the paraffin blocks, sections were cut, and placed on the glass slides as described previously96. Prepared slides were baked 30 minutes in 65°C, followed by the deparaffinization with xylol (2x10 min), graded isopropanol (2x3 mins 100%, 2x3 mins 96% and 3 mins 70%), and rehydration (2x5 mins in H202). Samples were boiled in the sodium citrate buffer (10mM sodium citrate with 0.05% Tween 20 in PBS at pH6) in the pressure cooker for 4 mins, as described previously129. Samples were permeabilized and blocked with the Perm/Block Buffer (1%BSA+0.2%TritonX100+10% goat serum in PBS) for 25 minutes in RT. All antibodies were diluted in 4% goat serum + 0.05% Tween20 in PBS, and 50 μl of antibodies per sample were used. Primary antibodies for IL-1β (10 μg/mL, mouse anti-IL1β, Abcam, Cambridge, UK), caspase-1 (1:50, rabbit anti-caspase 1, Cell Signaling, Danvers, USA), and RIG-I (1:250, mouse anti RIG-I, Santa Cruz Biotechnology, Santa Cruz, USA) were incubated in the wet chamber overnight at 4°C. Appropriate mouse and rabbit isotype controls, in the concentrations corresponding to the antibodies, were used (Dako, Agilent, Santa Clara, USA). Subsequently, samples were incubated with the goat anti-rabbit Alexa Fluor 488 (for caspase-1), and the goat anti-mouse Alexa Fluor 546 (for IL-1β and RIG-I) secondary antibodies for 60 minutes in RT in the concentration of 1:1000 (Invitrogen, Waltham, USA). After 3 minutes’ incubation in 1% PFA in room temperature, samples were mounted in the ProLong Gold mounting medium containing DAPI (Life Technologies, Carlsbad, USA) according to the manufacturer’s instructions, analyzed with a Zeiss LSM780 (Zeiss, Oberkochen, Germany) and Zen Software (Zeiss, Oberkochen, Germany). All pictures were taken at the 40x magnification and are presented as the maximal projection (orthogonal projection) from 4 z-stacks, with appropriate scale bar. For the quantification, 10 equal squares from epithelial areas of the tissue (assessed using Haematoxylin and Eosin staining of the adjacent slice) from stained and isotype control samples were measured for signal intensity and averaged. Isotype control signal was further subtracted from the intensity of the stained samples, and values from 10 squares per sample were presented as the mean fluorescent intensity (MFI) of the protein expression.
Transcriptome analyses
Next generation sequencing (NGS) from the differentiated human bronchial epithelial cells (HBECs) from control individuals and patients with asthma (cohort A) was performed as previously described107. Briefly, total RNA was isolated with a RNeasy Plus Micro Kit (Qiagen, Hilden, Germany). Library was prepared with the TruSeq Stranded mRNA Sample Prep Kit (Illumina, San Diego, USA), and sequenced on the Illumina HiSeq 4000 platform. Description of the study subjects is presented in the Supplementary Table S7.
HBECs from 6 control individuals and 6 patients with asthma, infected with RV-A16 in the MOI 10 for 24h were harvested and sequenced with the use of Illumina HiSeq 2000 platform, as described previously66. The mRNA expression data are publicly available at the Gene Expression Omnibus platform (https://www.ncbi.nlm.nih.gov/) under the accession number: GSE6114166.
Bronchial brushings from control individuals and patients with asthma before and after experimental RV-A16 infection in vivo were analyzed by Affymetrix HuGene 1.0 array according to the manufacturer’s instructions and Transcriptome Analysis Console v4.0 (Santa Clara, United States).
Data analysis
Distribution (normality) of the data was assessed with Shapiro-Wilk test. One-way ANOVA (Kruskal-Wallis test), RM one-way ANOVA (Friedman test) or mixed-effects model tests were performed for more than three groups comparisons depending on the data relation (paired/not-paired) and distribution (normal/not- normal). Two-tailed paired/not-paired t-test or Wilcoxon/U-Mann-Whitney tests were performed for two groups comparisons depending on the data relation and distribution. The data are presented as the mean ± SEM, with the number of samples in each experiment indicated in the figure description. IL-1β expression in BAL fluid from cohort SIBRO was analyzed with the Welsh’s test. Correlations regarding inflammasome and antiviral signatures between in vitro and in vivo RV-A16-infection models were calculated with Pearson’s correlation coefficient, or Spearman’s rank correlation coefficient (depending on data distribution) between genes significant both in vivo and in vitro (red) and gene significant in at least one type of experiment (black). All differences were considered significant when p≤0.05 and defined as *p-value≤0.05, **p-value≤0.01, ***p- value≤0.001, ****p-value≤0.0001. Statistical analysis was performed with the Prism 9 software (Redmond, USA).
Transcriptome data were processed with the workflow available at https://github.com/uzh/ezRun, with the significance threshold for differentially expressed genes set to p-value<.05 calculated for the entire gene lists in each project using the edgeR R package130. Microarray data was analyzed by the following Bioconductor microarray analysis workflow https://www.bioconductor.org/packages/release/workflows/vignettes/arrays/inst/doc/arrays.html.
Differentially expressed probe was identified by the limma R package with empirical Bayes estimation. Threshold for significance for transcriptome data presented on the figures are as follows: p-value: *<0.05; **<0.005; ***<0.0005, ****<0.00005. Heatmaps display normalized gene expression across the gene in the groups (row normalization). Additionally, top 100 genes upregulated after RV-A16 infection in the HBECs from control individuals and patients with asthma from GSE6114166 were analyzed for the enriched pathways using Metacore software version 20.3.70200 (Thomson Reuters, Toronto, Canada) (Supplementary Table S1). Inflammasome-mediated immune responses and antiviral responses gene sets were curated from GSEA and MSigDB Database (Broad Institute, Massachusetts Institute of Technology, and Reagent of the University of California, USA). Full sets of analyzed genes are described in the Supplementary Table S10.
Proximity Extension Assay (PEA) normalizes protein expression (NPX) data were analyzed with the use of the internal Shiny App Olink data analysis toolkit http://46.14.201.130:3838/OlinkR_shinyApp. The statistical comparison of protein expression between groups was performed with the Bioconductor limma package131. The fold change and p-value were estimated by fitting a linear model for each protein. Proteins with p- value<0.05 were considered significant. Additionally, for Target 96 Inflammation panel data are presented as: i) heatmaps of curated signatures of inflammasome-mediated immune responses and antiviral responses (Supplementary Table S10) and ii) protein interactions and pathways analysis prepared using the STRING (version 11.0)132, and further processed with the Cytoscape software (version 3.8.2)133 (Supplementary Table S3). The list of all proteins available for PEA measurements at the moment of the current analysis, that were used as a background reference for STRING analyses for targeted proteomics data is presented in Supplementary Table S11.
4. Data availability
Transcriptome data from bronchial brushings from control individuals and patients with asthma experimentally infected with RV-A16 has been submitted to the NCBI GEO: GSE185658 and will be publicly available at the time of publication. All other data are included in the Online Supplement or are available from the corresponding author upon request.
5. Code availability
The codes for transcriptome data analysis are available here: https://github.com/uzh/ezRun (NGS), https://github.com/ge11232002/p1688-Ula (microarray). Code for Proximity Extension Assay data analysis is available from the corresponding author upon request. All codes will be publicly available at the time of publication.
Data Availability
Transcriptome data from bronchial brushings from control individuals and patients with asthma experimentally infected with RV-A16 has been submitted to the NCBI GEO: GSE185658 and will be publicly available at the time of publication. All other data are included in the Online Supplement or are available from the corresponding author upon request. The codes for transcriptome data analysis are available here: https://github.com/uzh/ezRun (NGS), https://github.com/ge11232002/p1688-Ula (microarray). Code for Proximity Extension Assay data analysis is available from the corresponding author upon request. All codes will be publicly available at the time of publication.
8. Authors contributions
Concept design: MiSo. Data collection: UR, AE, MiSo, NS, AH, PaWe, SS, AD, PaWa, BR, JRC, DZ, MH, JB, MaSa, LOM, DJJ, MRE, TK. Analysis/interpretation of data: UR, AE, MiSo, GT, CA. Write/intellectual contribution: UR, AE, MiSo, CA, MM, TV, SJ. All authors approved the final version of the manuscript.
9. Ethics declaration
CA reports research grants from Allergopharma, Idorsia, Swiss National Science Foundation, Christine Kühne- Center for Allergy Research and Education, European Commission’s Horison’s 2020 Framework Programme “Cure”, Novartis Research Institutes, Astrazeneca, SciBase, Stanford University SEAN Parker Asthma and Allergy Center; advisory board of Sanofi/Regeneron, GSK and Novartis, consulting fees from Novartis; Editor- in-Chief Allergy, Co-Chair EAACI Environmental Science in Allergic Diseases and Asthma Guidelines. AE reports National Science Centre Grant No. 2020/37/N/NZ5/04144, National Centre for Research and Development No. STRATEGMED2/269807/14/NCBR/2015, National Centre for Research and Development (POLTUR3/MT- REMOD/2/2019). DJJ reports advisory board and speaker fees from AstraZeneca, GSK and Sanofi and research grants from AstraZeneca. SLJ reports grants/contracts from European Research Council ERC FP7 grant number 233015, Chair from Asthma UK CH11SJ, Medical Research Council Centre grant number G1000758, NIHR Biomedical Research Centre grant number P26095, Predicta FP7 Collaborative Project grant number 260895, NIHR Emeritus NIHR Senior Investigator; consulting fees from Lallemand Pharma, Bioforce, resTORbio, Gerson Lehrman Group, Boehringer Ingelheim, Novartis, Bayer, Myelo Therapeutics GmbH; patents issued/licensed: Wark PA, Johnston SL, Holgate ST, Davies DE. Anti-virus therapy for respiratory diseases. UK patent application No. GB 0405634.7, 12March 2004. Wark PA, Johnston SL, Holgate ST, Davies DE. Interferon-Beta for Anti-Virus Therapy for Respiratory Diseases. International Patent Application No. PCT/GB05/50031, 12 March 2004. Davies DE, Wark PA, Holgate ST, JohnstonSL. Interferon Lambda therapy for the treatment of respiratory disease. UK patent application No. 6779645.9, granted15th August 2012; Participation on a data safety monitory board or advisory board: Enanta Chair of DSMB, Virtus Respiratory Research Board membership. MM reports personal payments from Astra Zeneca, GSK, Berlin-Chemie/Menarini, Lek-AM, Takeda, Celon and support for attending meetings from Astra Zeneca, GSK, Berlin-Biochemie/Menarini. UR reports board secretary position of Working Group of Genomics and Proteomics of the European Academy of Allergy and Clinical Immunology (EAACI). JRC reports Pre-doctoral grant FPI from Universidad CEU San Pablo, Swiss European Mobility Program grant from University of Zurich, EAACI Mid-term Fellowship. MiSo reports research grants from Swiss National Science Foundation, GSK, Novartis and speakers fee from AstraZeneca and board secretary position of the Basic and Clinical Immunology Section of the European Academy of Allergy and Clinical Immunology (EAACI). SS reports funding from National Center of Competence in Research (NCCR) on RNA and Disease to VT (https://nccr-rna-and-disease.ch/). VT reports grant from Swiss National Science Foundation. All other authors report no conflict of interest regarding this work.
Supplementary Tables
7. Acknowledgements
The Authors are grateful to all patients, clinicians and research staff involved in securing and processing patients samples. We would like to thank David Mirer for rhinovirus plaque assay analyses. HDM extract was a gift from Allergopharma AG. This work was supported by the Swiss National Science Foundation (SNSF) grant nr 310030_189334/1, European Respiratory Society Long-term research fellowships, European Academy of Allergy and Clinical Immunology long-term research fellowship and GSK research grant (to MS), and the SNF grant nr 320030_176190 and GSK (to CA). The in vivo work was supported by the European Research Council (ERC FP7 grant number 233015); a Chair from Asthma UK (CH11SJ); the Medical Research Council Centre (grant number G1000758); National Institute of Health Research (NIHR) Biomedical Research Centre (grant number P26095); Predicta FP7 Collaborative Project (grant number 260895); and the NIHR Biomedical Research Centre at Imperial College London. SLJ is an Emeritus NIHR Senior Investigator. The in vivo transcriptome work was supported by GSK. UR, AE and MM were supported from funds from the Leading National Research Centre (KNOW) in Bialystok, Poland. JRC was supported by an FPI-CEU predoctoral fellowship and the Swiss-European Mobility Programme.
Footnotes
↵& Senior Co-Authorship
Abbreviations ACE2, angiotensin-converting enzyme 2; AIM2, absent in melanoma 2; ARTN, artemin; ASC, apoptosis-associated speck-like protein containing CARD; BAL, bronchoalveolar lavage; CARD, caspase-recruitment domain; CCL, C-C motif chemokine; CDC, carboxy-terminal domain; CLRs, C-type lectin receptors; COVID-19, coronavirus disease 2019; CXCL, C-X-C motif chemokine; dsRNA, double-stranded RNA; EN-RAGE, extracellular newly identified RAGE-binding protein; FCεRI, Fc epsilon receptor; HBECs, human bronchial epithelial cells; IAV, influenza A virus; ICAM-1, intracellular adhesin molecule-1; IFN, interferon; IgE, immunoglobulin E; IKKε, inhibitor-κB kinase epsilon; IL, interleukin; IL-1R, IL-1 receptor; ISG, interferon-stimulated genes; MAPK, mitogen-activated protein kinase; MAVS, mitochondrial antiviral-signaling protein; MCP, monocyte chemoattractant protein; MDA5, melanoma differentiation-associated protein 5; MOI, multiplicity of infection; MX1, MX dynamin like GTPase 1; MyD88, MyD88 innate immune signal transduction adaptor; NL, nasal lavage; NFκB, nuclear factor kappa-light-chain-enhancer of activated B cells; NGS, next generation sequencing; NLRC5, NOD- like receptor family CARD domain containing 5; NLRP3, NLR family pyrin domain containing 3; Nsp, non-structural proteins; OAS1, 2’- 5’-Oligoadenylate Synthetase 1; ORF, open read frames; PARs, protease activating receptors; RV, rhinovirus; RV-A16, rhinovirus A16; PRR, patter recognition receptor; RIG-I, retinoic acid-inducible gene I; ROS, reactive oxygen species; SARS-CoV-1, severe acute respiratory syndrome virus 1; SARS-CoV-2, severe acute respiratory syndrome virus 2; SLAMF1, signaling lymphocytic activation molecule precursor; ssRNA, single-stranded RNA; TBK1, TANK binding kinase 1; TLR, toll-like receptor; TNF, tumor necrosis factor; TNFRSF9, TNF receptor superfamily member 9; TRAF, TNF receptor associated factors; TRANCE, TNF-related activation-induced cytokine; UV-RV-A16, UV-light-inactivated rhinovirus A16; VSV, vescicular stomatitis virus.
References
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.
- 37.
- 38.
- 39.↵
- 40.↵
- 41.
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵
- 120.↵
- 121.
- 122.↵
- 123.↵
- 124.↵
- 125.↵
- 126.↵
- 127.↵
- 128.↵
- 129.↵
- 130.↵
- 131.↵
- 132.↵
- 133.↵