A paucigranulocytic asthma host environment promotes the emergence of virulent influenza viral variants
Abstract
Influenza virus has a high mutation rate, such that within one host different viral variants can emerge. Evidence suggests that influenza virus variants are more prevalent in pregnant and/or obese individuals due to their impaired interferon response. We have recently shown that the non-allergic, paucigranulocytic subtype of asthma is associated with impaired type I interferon production. Here, we seek to address if this is associated with an increased emergence of influenza virus variants. Compared to controls, mice with paucigranulocytic asthma had increased disease severity and an increased emergence of influenza virus variants. Specifically, PB1 mutations exclusively detected in asthmatic mice were associated with increased polymerase activity. Furthermore, asthmatic host-derived virus led to increased disease severity in wild-type mice. Taken together, these data suggest that at least a subset of patients with asthma may be more susceptible to severe influenza and may be a possible source of new influenza virus variants.
Introduction
Similar to many other RNA viruses, the influenza virus RNA polymerase has poor proofreading activity, resulting in a high rate of viral mutations (Van den Hoecke et al., 2015). Accordingly, even if an influenza virus infection occurs with a defined molecular clone, there will be rapid generation of viral variants that possess one or more mutations relative to the original virus strain. Many viral mutants are of reduced or unaltered virulence compared to the wild-type (WT) parental viral strain (Van den Hoecke et al., 2015). However, viral variants that display increased virulence can also emerge. Modelling studies suggest that there are a higher number of influenza virus variants in individuals that suffer from severe disease (Reperant et al., 2014). This may reflect the fact that mutations that are of increased virulence but of reduced fitness are able to replicate in immunocompromised individuals and not in individuals that have a fully competent immune response. Alternatively, in individuals with severe disease there may be additional rounds of virus replication, which represent additional opportunities for viral mutagenesis (Reperant et al., 2014). Consistent with these observations, pregnant mice infected with influenza virus have an increased rate of viral mutations compared to their non-pregnant counterparts (Engels et al., 2017). Similarly, the impaired interferon response of obese mice is thought to drive the emergence of virulent influenza virus variants (Honce et al., 2020). These data raise the question whether asthma, a host co-morbidity frequently associated with an impaired interferon response (Wark et al., 2005; Contoli et al., 2006), also facilitates the emergence of influenza virus variants.
The impact of pre-existing asthma on the severity of influenza virus infection is controversial. During the 2009 H1N1 pandemic asthma was identified as the single most common underlying medical condition in individuals hospitalised with influenza (Nguyen-Van-Tam et al., 2010; Morris et al., 2012). Individuals with asthma are known to have a greater risk of (i) admission to the hospital (Neuzil et al., 2000; O'Riordan et al., 2010; Van Kerkhove et al., 2011; Miller et al., 2008), (ii) being treated in the intensive care unit (Bueving et al., 2004; Libster et al., 2010), (iii) requiring a greater number of outpatient visits (Morris et al., 2012; Neuzil et al., 2000; Miller et al., 2008), and (iv) developing pneumonia (O'Riordan et al., 2010; Dawood et al., 2011) upon influenza virus infection. Whilst some of these observations may be attributable to the fact that influenza virus can exacerbate the symptoms of asthma (Obuchi et al., 2013), Dawood and colleagues demonstrated that the majority of patients with asthma were hospitalised for influenza in the absence of an acute asthma attack (Dawood et al., 2011). However, a global analysis of pH1N1 cases suggests that once hospitalised, influenza patients with asthma have a greater survival rate than patients with other conditions (Van Kerkhove et al., 2011). Consistent with these findings, studies using a mouse model of disease actually suggest that mice with the cardinal features of asthma display reduced disease severity following infection with influenza virus (Furuya et al., 2015; Samarasinghe et al., 2014; Ishikawa et al., 2012).
The controversies in the literature regarding the interactions between asthma and influenza appear hard to reconcile. However, it is important to recognise that asthma is a heterogenous disease, consisting of many different phenotypes (Svenningsen and Nair, 2017). It is possible that different subtypes of asthma display differing degrees of synergism with influenza virus. For example, the aforementioned mouse studies used either ovalbumin (OVA), fungal allergens, or house dust mite for sensitisation (Furuya et al., 2015; Samarasinghe et al., 2014; Ishikawa et al., 2012), all of which induce an eosinophilic asthma which was characterised by high levels of IgE and eosinophilia. Whilst approximately 40% of adults suffer from this allergic ‘eosinophilic asthma’, a large percentage of asthma patients suffer from non-allergic asthma (i.e. where asthma is observed in the absence of eosinophilic asthma). One such subtype is paucigranulocytic asthma in which airway remodelling occurs in the absence of eosinophilia or neutrophilia.
Here, we sought to determine the role of paucigranulocytic asthma, using an established experimental mouse model (Arikkatt et al., 2017), in the generation of influenza virus variants. We chose to focus on paucigranulocytic asthma for this study as (i) the relationship between this specific phenotype of non-allergic asthma and severe influenza has yet to be studied, (ii) in some populations paucigranulocytic asthma can affect up to 40% of asthmatic patients (Schleich et al., 2013), and (iii) paucigranulocytic asthma is associated with an impaired interferon response (Arikkatt et al., 2017). We show that mice with a paucigranulocytic asthma phenotype (henceforth referred to as ‘asthmatic mice’) have increased intra-host influenza viral diversity compared to control mice. Importantly, two of the viral mutations that emerged in asthmatic mice were associated with increased viral replication and polymerase activity. Together, these data provide strong evidence that host–viral interactions are bidirectional and that host co-morbidities can influence the evolution of influenza virus.
Results
Pneumonia virus of mice (PVM)-infection of RAGE-deficient mice induces the cardinal features of paucigranulocytic asthma
To investigate the effect of asthma on the severity of influenza we employed a previously established model of paucigranulocytic asthma (Arikkatt et al., 2017). In this model, an early-life PVM infection followed by reinfection with PVM later in life in mice deficient in the receptor for advanced glycation endproducts (RAGE), encoded by the Ager gene, promotes airway remodelling characteristic of asthma (Arikkatt et al., 2017). To confirm this phenotype, the number of mucus-secreting cells and airway smooth muscle mass (ASM) was assessed in PVM-infected RAGE-deficient mice and mock-infected RAGE-deficient mice. Importantly, RAGE-deficient mice that received only an early life PVM infection were not included in these experiments as we have previously shown that an early life PVM infection, followed by a later life infection with influenza A virus (IAV), also results in an asthma-like phenotype (Arikkatt et al., 2017). Consistent with our previous findings (Arikkatt et al., 2017), an early and later life PVM infection in RAGE-deficient mice resulted in a significant increase in the number of mucus secreting cells and ASM mass (Figure 1A). To confirm that the observed structural changes were not associated with a marked pro-inflammatory phenotype or strong interferon response, various cytokines were assessed in the lungs of asthmatic and non-asthmatic mice. Consistent with our previous findings (Arikkatt et al., 2017) asthmatic mice had significantly lower pulmonary levels of IFN-γ, CXCL1, TNF-α, CCL2, IL12, IL1b, IL10, and IFN-α compared to non-asthmatic control mice (Figure 1B). These data indicate that this is an appropriate mouse model to study the effects of asthma, and an impaired interferon response, on the emergence of influenza virus variants.

Pneumonia virus of mice (PVM)-infection (at 7 days of age and reinfection at 35 days of age) of RAGE-deficient mice induces the cardinal features of asthma.
(A) (Left) Percentage of airway epithelial cells (AECs) producing mucus was quantified by immunohistochemistry at 7 days post-reinfection (i.e. 42 days of age) from uninfected RAGE-deficient mice (non-asthmatic) and infected RAGE-deficient mice (asthmatic). (Right) Lung sections from non-asthmatic and asthmatic mice at 42 days of age were stained for smooth muscle actin by immunohistochemistry and the airway smooth muscle (ASM) area was calculated relative to the basement membrane (BM) length of small airways. (B) Cytokine levels in the lung of RAGE-deficient mice 7 days post-PVM or mock infection (i.e. 42 days of age). Statistical analysis was performed as described in Materials and methods with *p≤0.05. Data are pooled from at least two independent experiments and shown as mean ± SEM.
Mice with the features of paucigranulocytic asthma had increased IAV severity
To the best of our knowledge the association between non-allergic asthma and influenza has not been investigated. Accordingly, we sought to investigate the severity of IAV infection in the non-allergic model of asthma described above. Asthmatic mice infected with IAV lost significantly more body weight over time compared to non-asthmatic control mice (Figure 2A). This was associated with increased IAV replication in the lungs of asthmatic mice at 6 days post-infection (d.p.i)(Figure 2B). While there was an increase in the histopathology score of mock infected asthmatic mice compared to non-asthmatic mice, there was no observable difference following IAV infection (Figure 2C). Furthermore, asthmatic mice showed a significantly impaired type I interferon response, with a downregulation of both IFN-α and interferon inducible MX1 following IAV infection (Figure 2D).

RAGE-deficient asthmatic mice experience more severe influenza than non-asthmatic mice.
RAGE-deficient non-asthmatic and asthmatic mice were infected with 100 (PFU) of H1N1/Auck/09. Mock control groups received PBS only. (A) Percentage weight loss of influenza A virus (IAV) or mock infected mice. Weights are displayed as percentage of a mouse’s weight at the time of infection. Each data point represents mean ± SEM, with at least 17 mice per group. (B) Viral titres present in lung homogenate at 4 and 6 days post-infection (d.p.i.) in IAV and mock-infected mice. (C) Histopathology score of lung sections for vascular changes, bronchitis, interstitial inflammation, alveolar inflammation, pneumocyte hypertrophy, and pleuritis combined. (D) Cytokines in lung homogenate at 4 and 6 d.p.i. Data are normalised to Gapdh expression and fold change was calculated using the ΔΔCt method, expressed relative to mock infected mice. Statistical analysis was performed as described in Materials and methods with *p≤0.05, comparing non-asthmatic IAV to asthmatic IAV. Data are pooled from at least two independent experiments and shown as mean ± SEM (A–D), with a data point representing one mouse (B–D).
An asthmatic host environment led to increased non-synonymous IAV mutations
Having established that the asthma phenotype employed herein was associated with increased IAV severity we next sought to determine if asthmatic mice had an increased number of IAV variants. To do so, we used deep sequencing to analyse viral sequences isolated from the lungs of infected mice at 4 and 6 d.p.i. Mutations were only included in the analysis if they constituted >1% of the total viral population and were detected in a minimum of 1000 reads following quality control. Shannon’s entropy was determined for each of the eight influenza genome segments in both 4 and 6 d.p.i. samples (Figure 3). To account for any bias from difference in sequencing coverage, the read count for each segment was analysed. No significant difference was observed in sequence coverage for any of the influenza virus genes between the two treatment groups (Figure 3—figure supplement 1). Asthmatic host-derived samples showed significantly more within-host viral diversity of the PB1 segment at 6 d.p.i. (Figure 3B). To determine if the mutations were associated with a specific section of the PB1 gene, their frequency was visualised against the nucleotide position, where no discernible trend was observed (Figure 4).

Asthmatic hosts produced more diverse viral variants in PB1.
Influenza viral RNA isolated from the lung tissue of asthmatic (n = 27) and non-asthmatic (n = 25) mice was analysed for viral variants at 4 and 6 days post-infection (d.p.i). (A–H) Measurement of alpha diversity presented as Shannon entropy of viral mutations at a frequency above 1% by genome segment at the nucleotide level. Statistical analysis was performed as described in Materials and methods with *p≤0.05. Data are shown as mean ± SEM.

Asthmatic hosts produced more single nucleotide variants (SNVs) in the PB1 genome segment.
Influenza viral RNA isolated from the lung tissue of asthmatic (n = 27) and non-asthmatic (n = 25) mice was analysed for viral variants. Each data point represents the number of viral variants detected at that nucleotide position within the viral genome PB1 segment. bp: base pair.
To determine if any of the mutations detected in the asthmatic mice could have contributed to the increased IAV severity observed, we focused on unique (occurring specific to either non-asthmatic or asthmatic hosts) non-synonymous mutations within the PB1 segment, as these were the most likely to be of functional significance. This produced a list of 404 unique mutations across all groups (Supplementary file 2). We then refined our analysis to only include non-synonymous mutations that were detected in two or more mice. Overall, there was an increase in unique mutations detected in asthmatic compared to non-asthmatic host samples, with 32 and 5, respectively (Table 1). Interestingly, the bulk of the mutations were focused around the priming loop, thumb, and C-terminal extension of PB1.
Frequency of exclusive non-synonymous viral mutations detected in PB1 at 4 and 6 days post-influenza A virus (IAV) inoculation and their frequency in clinical samples.
Non-asthmatic | Asthmatic | Clinical samples (n = 51 per group) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Day 4 (n* = 13) | Day 6 (n = 12) | Day 4 (n = 13) | Day 6 (n = 14) | Non-asthmatic | Asthmatic | p-Value‡ | ||||
M1I | 50%† | 0% | 0% | |||||||
M1L | 17% | 0% | 0% | |||||||
D2H | 14% | 0% | 0% | |||||||
N77I | 14% | 0% | 0% | |||||||
P74A | 14% | 0% | 0% | |||||||
E75D | 29% | 0% | 0% | |||||||
E75V | 14% | 0% | 0% | |||||||
D76A | 14% | 0% | 0% | |||||||
D76H | 14% | 0% | 0% | |||||||
N77V | 14% | 0% | 0% | |||||||
E78G | 14% | 0% | 0% | |||||||
S80H | 14% | 0% | 0% | |||||||
I205M | 14% | 0% | 0% | |||||||
I248V | 14% | 0% | 0% | |||||||
M616I | 14% | 0% | 0% | |||||||
R621P | 14% | 0% | 0% | |||||||
C625S | 14% | 0% | 0% | |||||||
A652P | 14% | 0% | 0% | |||||||
S654N | 14% | 4% | 18% | 0.0514 | ||||||
Y657H | 14% | 0% | 0% | |||||||
N671H | 17% | R670S | 14% | 0% | 0% | |||||
R672P | 15% | 0% | 0% | |||||||
S673A | 17% | 0% | 0% | |||||||
I682F | 15% | 0% | 0% | |||||||
D685E | 15% | 0% | 0% | |||||||
D685H | 14% | 0% | 0% | |||||||
M688I | 14% | 0% | 0% | |||||||
Y689D | 21% | 0% | 0% | |||||||
Y689R | 14% | 0% | 0% | |||||||
Q690S | 15% | 0% | 0% | |||||||
K691R | 15% | 4% | 16% | 0.092 | ||||||
C692S | 14% | 0% | 0% | |||||||
C693N | 14% | 0% | 0% | |||||||
N694R | 14% | 0% | 0% | |||||||
L695I | 15% | 0% | 0% | |||||||
R727T | 14% | 0% | 0% |
-
*n = Total number of mouse samples.
†Percentage of mice with identified mutation.
-
‡Statistical comparisons were made using a two-sided fisher’s exact test.
We next sought to determine if any of these IAV mutations detected in mouse samples had also been reported in pH1N1 IAV patients with asthma. To do so, post-2009 H1N1 viral sequences from patients that listed ‘asthma’ as a medical co-morbidity (n = 51) were obtained from the Influenza Research Database (IRD). This was compared to viral sequences obtained from 51 age and sex-matched pH1N1 patients on the IRD with no medical co-morbidities. Strikingly, of the unique mouse single nucleotide variants, two that were unique to asthmatic mouse samples were also detected in clinical samples at a higher frequency in asthmatic patient samples (Table 1). Interestingly, of the mutations identified in the clinical consensus sequence samples, the K691R mutation reached a high enough frequency in two asthmatic mouse samples to alter the consensus sequence.
Asthmatic host-derived samples showed increased replication kinetics in vitro
Having established that an asthmatic host generated more viral variants, specifically within PB1, we next sought to assess if these mutations influenced the replicative ability of asthmatic and non-asthmatic mouse-derived IAV samples in vitro. To do so, we infected Madin–Darby canine kidney (MDCK) cells, human alveolar basal epithelial cells (A549 cells), and Immortalized Mouse Mammary Epithelial Cells (iMMEC) at a multiplicity of infection (MOI) of 0.01 with either lung homogenates derived from three asthmatic mice or three non-asthmatic mice. The viruses in the asthmatic lung homogenates replicated to significantly higher titres than those from non-asthmatic lung homogenates at 24 hr post-infection in MDCK cells (Figure 5A). A similar trend was observed in mouse and human epithelial cells (Figure 5B and C) but with a delayed rate of replication.

Increased replication of influenza A virus (IAV) was observed in samples derived from asthmatic hosts in vitro.
IAV isolated from asthmatic hosts replicates more quickly than non-asthmatic hosts in (A) Madin–Darby canine kidney (MDCK) cells, (B) Immortalized Mouse Mammary Epithelial Cells (iMMEC), and (C) adenocarcinomic human alveolar basal epithelial cells (A549 cells). Each cell type was infected with the indicated viruses at a multiplicity of infection (MOI) of 0.01 (n = 3 wells/virus/time point). Mock infected data not shown. Data were analysed using a one-way ANOVA with Holm–Sidak's multiple comparisons test and are represented as mean ± SEM. *p<0.05. Dashed line indicates the detection limit of the assay.
Interestingly, comparing two of the viral samples pre and post in vitro replication in MDCKs we did not detect a change in the consensus sequence of the daughter population after in vitro replication compared to the consensus sequence of the parent. In terms of minority viral populations, in one daughter population no viral variants (with non-synonymous mutations) were detected that constituted between 5 and 49% of the viral population. In the second daughter population, one non-synonymous mutation (PB1 I205M) was detected that constituted 28% of the viral population. However, this mutation was present in the consensus sequence of the parent isolate, albeit at a higher percentage (95%). Together, these data suggest that the increased replication of viral isolates from asthmatic mice was not associated with a large number of additional in vitro adaptations.
PB1 C693Y and P651Q may increase the stability of the influenza RNA polymerase protein dimer
We next sought to elucidate if any of the observed mutations identified in the asthmatic group could be responsible for increased viral replication. To do so, four mutations were selected for structural modelling; two that appeared in multiple mouse samples as well as clinical samples (S694N and K691R), and two that were found in areas likely to affect function (C693Y and P651Q) (Supplementary file 2).
None of the mutations were predicted to have a major destabilising effect and all mutations maintained the secondary structures within PB1. The Foldx results predicted S654N (ΔΔG −0.07 kcal/mol) and K691R (ΔΔG 0.43 kcal/mol) to be neutral and C693Y (ΔΔG −1.86 kcal/mol) and P651Q (ΔΔG −1.78 kcal/mol) to be stabilising as defined by Studer et al., 2014. As the C-terminus extension interacts with PB2, Foldx was used to highlight potential interactions between the C-terminal mutated residues and residues of PB2. The structural modelling of C693Y identified a side chain–side chain interaction in the mutant between PB1 C693Y and PB2 K32 in the third α-helix (Figure 6A). In contrast, PB1 WT C693 was only predicted to interact with resides within the first PB1 α-helix. Together, these data suggest that the C693Y mutation identified exclusively in asthmatic hosts could alter the interaction between PB1 and PB2 within the polymerase heterotrimer. No close contacts were identified between K691R and residues of PB2 (Figure 6A). S654N and P651Q reside within the conserved anti-parallel β-hairpin with P651Q sitting in the tip of the β-hairpin (Figure 6B). Both mutations within the β-hairpin changed the orientation of the residue, which could suggest a functional change.

Crystal structures of the mutations within PB1.
(A) Mutations C693Y and K691R. The molecule is displayed in surface view with cartoon representation of the protein backbone (PB1 magenta and PB2 sky blue). The α-helices are labelled and PB1 C693Y and K691R and PB2 K32 are shown as sticks in the respective colours of the subunits with the wild-type (WT) residues in pale yellow. (B) Mutations P651Q and S654N within the conserved β-hairpin (cartoon representation in pale green). The sites of mutations are shown as sticks with the WT residues in pale yellow.
PB1 C693Y and P651Q increase viral polymerase activity
We next sought to investigate the individual effect of the PB1 C693Y and P651Q mutations on viral polymerase activity using a minigenome assay, given their predicted increase in stability. Both the PB1 C693Y and P651Q mutations within the PB1 gene increased polymerase activity when compared to the WT virus polymerase (Figure 7). Together these data suggest that the C693Y and P651Q mutations identified exclusively in asthmatic hosts may alter the interaction between PB1 and PB2 in a stabilising manner, resulting in increased polymerase activity.

The PB1 C694Y and P651Q mutations showed increased polymerase activity compared to the wild type (WT).
Influenza RNA polymerase-driven activity is increased compared to WT in HEK 293 T cells transfected with either the C693Y or P651Q mutation in PB1. Data are pooled from four independent experiments (with five biological replicates per group, per experiment). Results are presented as a fold change relative to WT. Statistical significance was determined as described in the Materials and methods. Data shown as mean ± SEM. *p<0.05.
Asthmatic host-derived IAV exhibits greater disease severity and viral diversity upon infection of a WT host
Vulnerable patient groups may be a source for the emergence of new viral strains (Honce et al., 2020). As such, we next sought to determine if asthmatic host-derived IAV would induce greater disease severity in WT mice. We inoculated WT mice with either asthmatic or non-asthmatic host-derived IAV and monitored them for 4 days, at which point they were euthanised. Mice infected with asthmatic host-derived IAV had significantly more weight loss (Figure 8A), lower blood oxygen saturation (Figure 8B), and higher viral titres (Figure 8C) compared to mice infected with non-asthmatic host-derived IAV. There was also a non-significant trend towards increased viral mRNA present (Figure 8D). Mice infected with asthmatic host-derived IAV had a significant upregulation in various anti-viral and pro-inflammatory cytokines (Figure 8E).

Mice infected with asthmatic host-derived virus experience more severe influenza than mice infected with non-asthmatic host-derived virus.
C57BL/6 mice were infected with 2000 PFU of either asthmatic or non-asthmatic host-derived influenza A virus (IAV). (A) Percentage weight loss of mice following IAV infection. Weights are displayed as percentage of a mouse’s weight at the time of infection. Each data point represents mean ± SEM (n = 6 per group). (B) Percentage blood oxygen saturation of mice at 4 days post-infection (d.p.i.). (C) Viral titres present in lung homogenate at 4 d.p.i. in IAV-infected mice. (D) Viral mRNA detected by qPCR from murine lung homogenate at 4 d.p.i. Dashed line indicates the detection limit. (E) Cytokines in lung homogenate at 4 d.p.i. Data are normalised to Gapdh expression and fold change was calculated using the ΔΔCt method, expressed relative to samples infected with non-asthmatic-derived virus. Statistical analysis was performed as described in Materials and methods with *p≤0.05. Data are shown as mean ± SEM.
Furthermore, upon infection of a non-asthmatic host with non-asthmatic-derived IAV, there is an observable divergence in the PB1 haplotype differences as the virus continues to adapt to its host (Figure 9A). By contrast, upon infection of a non-asthmatic host with asthmatic-derived IAV, there is a consolidation of the PB1 mutations that arose during the initial infection (Figure 9B). Together these data suggests that asthmatic hosts could facilitate the emergence of new strains that have increased pathogenicity.

Asthmatic host-derived virus shows consolidation of PB1 viral variants population across parent and daughter.
C57BL/6 mice were infected with 6 days post-infection (d.p.i.) murine lung homogenate from either non-asthmatic or asthmatic hosts. Multiple dimension scaling (MDS) analysis of the PB1 gene in (A) non-asthmatic-derived virus and (B) asthmatic-derived virus (each data point represents an individual haplotype). Data is representative of six samples per group. Parent: virus derived from either asthmatic or non-asthmatic Ager-/- host. Daughter: virus isolated from C57BL/6 mice lung following infection with host-derived virus.
Discussion
Asthma was identified as the most prevalent underlying host condition associated with hospitalisation for pH1N1 in 2009 (Nguyen-Van-Tam et al., 2010; Morris et al., 2012). However, the impact of pre-existing asthma on the severity of influenza remains unclear (O'Riordan et al., 2010; Dawood et al., 2011; Furuya et al., 2015; Samarasinghe et al., 2014; Ishikawa et al., 2012). Mouse studies modelling allergic asthma have suggested that asthma may mediate partial protection from severe influenza (Furuya et al., 2015; Samarasinghe et al., 2014; Ishikawa et al., 2012; Doorley et al., 2017). Here, using our previously established mouse model of non-allergic paucigranulocytic asthma (Arikkatt et al., 2017), we show that the impaired interferon response (characteristic of this subtype) is associated with increased influenza severity. These data suggest that further research, both clinical and pre-clinical, into the relationship between non-allergic asthma and influenza virus severity is warranted.
Recent data suggest that an impaired interferon response correlates with an increase in the emergence of influenza viral variants (Engels et al., 2017; Honce et al., 2020). Engels and colleagues demonstrated that an allogenic pregnant host environment fails to activate the innate immune response, including that of type I interferons during IAV infection (Engels et al., 2017). The failure to mount a robust immune response was associated with an increased viral load and a higher emergence of viral variants (Engels et al., 2017). More recently, Honce and colleagues demonstrated that the decreased type I interferon response within an obese micro-environment corresponded to increased viral replication and enhanced virulence compared to WT mice (Honce et al., 2020). Consistent with these findings, we observed increased viral variants in mice with the cardinal features of paucigranulocytic asthma, particularly within PB1. The absence of an early immune response in the asthmatic mice, that is, lack of type I interferon expression, may allow for the emergence of viral variants that would otherwise be suppressed in a healthy host. Further to this, non-asthmatic host-derived viruses showed host adaption and little sign of morbidity in WT mice, while asthmatic host-derived virus resulted in increased morbidity. This suggests that viral variants born from an asthmatic host could (transmission permitting) lead to the emergence of more pathogenic viral strains.
The influenza A viral RNA-dependent RNA polymerase (RdRp) forms dimers of heterotrimers consisting of PA, PB1, and PB2 (Fan et al., 2019), with areas of inter-subunit contact stabilising the trimers (Pflug et al., 2014; Sugiyama et al., 2009). Intra-host viral mutations in this complex may be deleterious, neutral, or advantageous in terms of viral pathogenesis. The β-hairpin within the PB1 thumb domain is responsible for viral replication initiation, and deletions in this region impair viral replication and growth (Oymans and te Velthuis, 2018).
PB1 S654N (a mutation identified in asthmatic clinical samples and those from asthmatic mice) and PB1 P651Q (a mutation identified in asthmatic mice) reside within this priming loop. Interestingly, the P651 residue has previously been shown to be the primary amino acid involved in coordinating terminal initiation during replication (Te Velthuis et al., 2016). Glutamine residues have been documented to interact with nucleotide bases (Luscombe et al., 2001), and given the position of PB1 P651Q, it is possible that it alters the interaction with the 3ʹ sugar-base of the vRNA or is beneficial to structural conformation during the initiation process. Our structural analysis of PB1 P651Q predicted a stabilising effect, which was consistent with an increase in polymerase activity relative to the WT. C693Y and K691R (both of which were identified in asthmatic mice) reside in the first α-helix of the C-terminus, which constitutes the major binding region with the N-terminus of PB2 and has an overall major stabilising effect on the heterotrimer (Pflug et al., 2014; Sugiyama et al., 2009; Te Velthuis and Fodor, 2016; Sugiyama et al., 2009; Poole et al., 2007). Deletion of any of the three α-helices within the C-terminus abolishes viral gene expression (Poole et al., 2007). Double and single mutations in the C-terminus of PB1 alter viral RNA yield (Sugiyama et al., 2009). Our structural analysis of PB1 C693Y predicted side chain interactions between the mutated residue and PB2 K32 which may increase the stability of the inter-subunits and, in turn, affect polymerase activity. Thus cumulatively, the data presented here suggest that the PB1 C693Y and P651Q mutations detected in asthmatic mice resulted in increased viral replication in various epithelial cell lines, and increased polymerase activity. It is well established that single nucleotide variants (SNVs) can drastically change the replication fidelity and pathogenicity of a virus (Lin et al., 2019; Pauly et al., 2017). As such, future research would benefit from more in-depth analysis of the effect of C693Y or P651Q within PB1 on the fidelity of the RdRp. Nevertheless, our results should be taken as proof of principal that an asthmatic host may be more likely to produce more virulent viral variants – a finding that holds wide consequences for the emergence of new viral strains.
It is important to recognise that there are several limitations associated with this study. Firstly, paucigranulocytic asthma remains a poorly understood and poorly characterised phenotype of asthma (Tliba and Panettieri, 2019). Accordingly, we are not aware of any other mouse model of paucigranulocytic asthma that could be used to validate these results. Similarly, there was insufficient information associated with the clinical sequences used herein to determine which asthma phenotypes each patient had. Indeed, we cannot eliminate the possibility that all patients suffered from allergic, rather than non-allergic, asthma. The clinical sequences analysed in the present study were also obtained as consensus sequences, and do not account for minority viral variant populations present in any individual. Additionally, only a single measure of genetic diversity was used (Shannon entropy) (Honce et al., 2020; Abd Raman et al., 2020) although the possibility of read depth bias was considered. Finally, it must be emphasised that these data cannot be extrapolated to other subtypes of asthma where the relationship with influenza virus, and intra-host viral variation, may be markedly different.
Nevertheless, despite these limitations these data suggest that in a mouse model of paucigranulocytic asthma there is an increased diversity of influenza viral variants some of which are associated with increased polymerase activity. Additionally, we have demonstrated that virus derived from an asthmatic host can cause more severe disease when introduced to a non-asthmatic host. Given these findings, additional studies monitoring individuals with asthma (and in particular those with non-allergic asthma) for increased influenza severity and the emergence of influenza virus variants may be of value.
Materials and methods
Viral strains and titrations
Request a detailed protocolPVM strain J3666 stocks were prepared from mouse lung homogenate as described previously (Garvey et al., 2005). A/Auckland/4/2009 (Auckland/09; H1N1) stocks were prepared in embryonated chicken eggs as previously described (Brauer and Chen, 2015). Viral titres were determined by plaque assays on Madin–Darby canine kidney (MDCK) cells, as previously described (Short et al., 2011).
Mouse strains
Request a detailed protocolRAGE-deficient mice (Ager-/-) were kindly provided by Prof. Ann-Marie Schmidt, (New York University Langone Medical Centre, NY, USA). Ager-/- and C57BL/6 mice were housed in individually ventilated cages, under alternating 12 hr light/dark periods with access to clean drinking water and food ad libitum. All animal experiments were approved by the University of Queensland Animal Ethics Committee (permit no. 071/17).
Establishing paucigranulocytic asthma
Request a detailed protocolTo model asthma development, we employed our previously described mouse model of paucigranulocytic asthma (Arikkatt et al., 2017). Specifically, 7-day old female and male Ager-/- mice were infected intranasally with 10 plaque forming units (PFU) of PVM in a volume of 10 μL. Thirty-five days after the primary infection, mice were reinfected intranasally with 100 PFU of PVM in a total volume of 50 μL. This model reflects the fact that asthma development is often a combination of a genetic deficiency and multiple respiratory tract infections (Arikkatt et al., 2017). ‘Non-asthmatic’ control mice were deficient in RAGE but were mock infected with 10 μL of PBS at 7 days of age.
In vivo IAV infection
Request a detailed protocolAt 48 days of age, all Ager-/- mice were infected with 100 PFU of Auckland/09(H1N1) intranasally in 50 μL phosphate buffered saline (PBS) or mock-infected with PBS. Six-week-old female C57BL/6 mice were infected intranasally with 2000 PFU of either asthmatic host-derived IAV or non-asthmatic host-derived IAV. All infections were performed intranasally under isoflurane-induced anaesthesia using Stinger Research Aesthetic Gas Machine (Darvall, AZ, USA). Mice were monitored daily for weight loss and clinical signs of disease. At the relevant time points, mice were euthanised with intraperitoneal injection of a lethal dose of phentobarbitol (270 mg/kg) (Virbac, NSW, Australia).
Blood oxygen saturation
Request a detailed protocolBlood oxygen saturation of C57BL/6 mice was measured at 4 d.p.i. using a collar sensor and Mouseox Plus pulse oximeter (Starr, PA, USA).
Histology and immunohistochemistry
Request a detailed protocolThe left lung lobe was fixed in 10% neutral buffered formalin, routine processed and embedded in paraffin, sectioned at 5 μm, and subsequently used for histology and immunohistochemistry. To identify mucus producing cells, sections were stained with Periodic Acid-Schiff (PAS) and then counterstained with Harris’ hematoxylin and mounted in DEPEX. Mucus production was scored based on the percentage of mucus secreting airway epithelial cells (AECs) relative to the total number of AECs. A minimum of five airways per sample were quantified. For immunohistochemistry, sections were pre-treated with 10% normal goat serum for 30 min. Sections were then probed with anti-α-SM actin (Sigma-Aldrich, MO, USA) overnight at 4°C. Following incubation with secondary antibody, immunoreactivity was developed with Fast Red (Sigma-Aldrich) and counterstained with Mayer’s haematoxylin. Airway smooth muscle (ASM) was quantified essentially as described previously (Loh et al., 2020). ASM was calculated as area per µm of basement membrane of small airways (defined as circumference <800 µm). Lungs were assessed for vascular changes, bronchitis, interstitial inflammation, alveolar inflammation, pneumocyte hypertrophy, and pleuritis (represented by a total score) by a veterinary pathologist who was blind to the study design. All stained slides were scanned using a digital slide scanner (Scanscope XT, Aperio Technologies, CA, USA) at 200× magnification.
Cytokine levels
Request a detailed protocolCytokine concentrations were assessed using LegendPlex mouse anti-viral panel (Biolegend, CA, USA) according to the manufacturer’s instructions.
Quantification of influenza virus titres
Request a detailed protocolMouse lung tissues were collected in Dulbecco’s minimum essential medium (DMEM; Gibco, NY, USA) and homogenised using a Qiagen Tissuelyser II (Qiagen, Hilden, Germany). The homogenate was centrifuged, and the supernatant was collected and stored at −80°C. Cell culture supernatant was collected at indicated time points and stored at −80°C. Viral titres were measured by plaque assay on MDCK cells, as described previously (Short et al., 2011).
RNA extraction, cDNA synthesis, and polymerase chain reaction (PCR)
Request a detailed protocolViral and host RNA was extracted from mouse lung samples, cDNA was synthesised, and real-time PCR was performed as described previously (Hulme et al., 2020). Forward and reverse primers used in the present study are supplied in Supplementary file 1. Gene expression was normalised relative to glyceraldehyde 3-phosphate dehydrogenase (Gapdh) expression and fold change was calculated using the ΔΔCt method (Schmittgen and Livak, 2008). Viral copy number was determined using IAV strain A/Puerto Rico/8/1934 H1N1 (PR8) virus matrix (M) gene cloned into pHW2000 plasmid and viral copy number was determined as described previously (Short et al., 2013).
Viral RNA for sequencing was extracted from mouse lung, and in vitro culture supernatant samples using the High Pure Viral RNA Kit (Roche, Basel, Switzerland) according to the manufacturer’s protocol and eluted in UltraPure DNase/RNase-Free Distilled Water (ThermoFisher Scientific, MA, USA). Viral cDNA was synthesised using the Transcriptor First Strand cDNA Synthesis Kit (Roche) essentially as described previously (Mancera Gracia et al., 2017). Briefly, two reactions were performed per sample, using either 2.5 µM CommonUni12A primer (GCCGGAGCTCTGCAGATATCAGCAAAAGCAGG) or 2.5 µM CommonUni12G primer (GCCGGAGCTCTGCAGATATCAGCGAAAGCAGG). Both reactions were incubated for 10 min at 65°C and cooled immediately on ice. RNA was reverse transcribed using Mastercycler Pro 6325 (Eppendorf, Hamburg, Germany). cDNA was pooled at a 1:1 ratio and amplified using the Phusion High-Fidelity DNA Polymerase kit (New England Biolabs, MA, USA) with 10 µM CommonUni13 primer (GCCGGAGCTCTGCAGATATCAGTAGAAACAAGG). cDNA was amplified using Mastercycler Pro 6325 (Eppendorf). MinElute PCR Purification Kit (Qiagen) was used to purify the PCR products. To elute, unheated UltraPure DNase/RNase-Free distilled water was used instead of elution buffer. The presence of the correct size PCR products was confirmed using a 1% agarose gel containing SYBR Safe DNA Gel Stain (ThermoFisher Scientific) at 100 Volt and visualised by the Amersham Imager 600 (GE Healthcare, IL, USA).
Next generation sequencing (NGS) for viral genes
Request a detailed protocolLibraries were prepared from viral cDNA amplicons using Nextera XT DNA library prep kit (Illumina, cat no. FC-131–1024) and were deep sequenced (0.5 Gb) on an Illumina Nextseq 2000 Platform at the Australian Centre of Ecogenomics (the University of Queensland, QLD, Australia).
NGS analysis and Shannon entropy
Request a detailed protocolThe haplotypes for each sample were reconstructed for each gene segment using a previously published pipeline (Cacciabue et al., 2020). In brief, FastQC (Andrews, 2010) was used for quality assurance of the NGS paired-end raw reads followed by BBtools (Bushnell, 2014), for removing and filtering adapters and low-quality reads. Bowtie2 (Langmead and Salzberg, 2012), an aligner tool to align the trimmed reads to the selected reference of the influenza strain (i.e. the inoculum), was then used. Samtools suite (Li et al., 2009) was used to sort, index, and generate depth and coverage statistics for read alignment files. Next, CliqueSNV (Knyazev, 2020) was used to infer the haplotypes and frequencies for all eight gene segments for each sample.
Shannon entropy (abundance-based diversity) was calculated using QSutils, an R package (Guerrero-Murillo and Font, 2020). During analysis the following assumptions were made: each paired-end read present in an alignment comes from a true viral haplotype from the original population, the occurrence of variants in a given gene segment is independent of the rest of the genes, and each haplotype in the actual population has an equal chance to be sampled.
Multiple dimension scaling (MDS) (LMDS V1.0 – R package) was utilised to visualise distance matrices to understand the dynamics of haplotype distance in a given quasi species population across donor and acceptor hosts (i.e. C57BL/6 mice infected with either 6 d.p.i. murine lung homogenate from asthmatic or non-asthmatic host) (Jombart, 2016). We used IRMA (iterative refinement meta-genome assembler) (Shepard et al., 2016) for variant calling and each of the variant in the coding sequence of PB1 gene was translated to the corresponding amino acid residue compared to the reference.
All scripts generated for this study can be found at https://github.com/akaraw/Hulme_et_al (Karawita, 2021) and the workflow for haplotype reconstruction can be found in Supplementary file 3.
Within-host alpha diversity was measured using Shannon’s entropy (H):
Identification of asthma associated mutations of IAV in clinical samples
Request a detailed protocolH1N1 PB1sequences were obtained from the NIAID IRD through the website at http://www.fludb.org. The following settings were used for selecting human surveillance data following the 2009 pandemic: date range ≥2009, subtype H1N1, all countries, all sexes, and age ranges; exclude laboratory strains and duplicate sequences. Only records containing clinical metadata were returned. Using these criteria, a total of 51 human clinical samples with asthma as a medical condition were observed. Fifty-one age and sex matched human clinical samples with no medical conditions were taken for comparison. Samples were aligned using MEGA-X (version 10.1) (Pennsylvania State University, PA, USA) and the presence or absence of the mutations identified in vivo was determined using amino acid position.
IAV mutagenesis
Request a detailed protocolThe PB1-C693Y and PB1-P651Q mutations were independently introduced into the PB1 gene segment of influenza A/Auckland/2009 (H1N1) using the Q5 Site-Directed Mutagenesis Kit (New England BioLabs) as per the manufacturer’s instructions. Mutations were confirmed using the PureYield Plasmid Miniprep System Kit (Promega, WI, USA) as per the manufacturer’s instructions and sequenced at the AEGRC Genotyping and Sequencing Facility at the University of Queensland. DNA was isolated using an EndoFree Plasmid Maxi Kit (Qiagen) as per the manufacturer’s instructions.
Cell culture
Request a detailed protocolHuman embryonic kidney 293 with large T antigen (HEK 293T) cells, MDCK cells, adenocarcinomic human alveolar basal epithelial cells (A549 cells), and Immortalized Mouse Mammary Epithelial Cells (iMMEC) were cultured in DMEM (Gibco) containing 10% foetal bovine serum (FBS; Sigma-Aldrich) and 1% penicillin streptomycin (Lonza, Basel, Switzerland) in a humidified atmosphere of 5% CO2 at 37°C. All cell lines were obtained from American Type Culture Collection (ATCC; Virginia, USA) and have no mycoplasma contamination to report.
In vitro viral replication kinetics
Request a detailed protocolMulticycle viral replication kinetics were determined by inoculating a confluent monolayer of cells with the lung homogenate of infected mice at a multiplicity of infection (MOI) of 0.01 for 1 hr. Cells were then washed with PBS three times and fresh medium (DMEM; Gibco) was added containing TPCK-trypsin (Thermofisher Scientific). At 0, 6, 24, 48, and 72 hr post inoculation, cell culture supernatant was collected and used to determine virus titres, as described above.
Influenza virus polymerase assay
Request a detailed protocolReporter plasmids, Fluc (0.25 µg) and pRL (0.01 µg) were co-transfected with expression plasmids encoding the AUCK/09 PA, PB1, PB2, PA, and NP (0.25 µg per plasmid) into HEK 293 T cells seeded on a 24-well plate using the calcium phosphate mediated method (Kwon and Firestein, 2013). Luciferase assays were performed using a dual-specific luciferase assay kit (Promega). A Renilla luciferase reporter vector was used as the internal control. Firefly luciferase activity was normalised to Renilla luciferase activity for each well. Luminescence was measured on a CLARIOstar Plus microplate reader (BMG Labtech, Ortenberg, Germany) in flat 96-white well plates (Corning, Amsterdam, the Netherlands).
Modelling of the mutations
Request a detailed protocolFor K691R and C693Y the crystal structure of the PB1–PB2 subunits from IAV (A/Puerto Rico/8/1934(H1N1)) was used as a model (PDB: 3A1G) (Sugiyama et al., 2009). For P651Q and S654N the structure of the influenza A bat polymerase complex (PDB: 4WSB) was used as a model (Pflug et al., 2014). The effects of conformational stability were estimated using FoldX (version 5) (Schymkowitz et al., 2005) and each mutation was engineered separately. For the FoldX analysis the subroutines RepairPDB, BuildModel, and PrintNetworks were used. Images were created with PyMOL (The PyMOL Molecular Graphics System, Version 2.4.0 Schrödinger, LLC).
Statistical analysis
Request a detailed protocolStatistical analyses were performed using Graph Pad Prism software (version 9.0.1) for Windows (GraphPad Software, CA, USA). Data were tested for normality using the Shapiro–Wilk test. Where data was normally distributed, data was analysed using a two-way ANOVA (with Tukey’s multiple comparison test), a one-way ANOVA (with Tukey’s multiple comparison test), an unpaired Student’s t-test, or a one sample t test as appropriate. Where data was not normally distributed, data was analysed using a Kruskal–Wallis test (with Dunn’s multiple comparisons test), a Mann–Whitney U-test or a Wilcoxon signed-rank test as appropriate. Clinical samples were analysed using a two-sided fisher's exact test.
Data availability
All methods and scripts used for sequencing analysis are outlined in the manuscript and supporting files. The raw sequencing data generated in this study is available on the European Nucleotide Archive with accession number PRJEB40027.
-
European Nucleotide ArchiveID PRJEB40027. A paucigranulocytic asthma host environment promotes the emergence of virulent influenza viral variants.
References
-
SoftwareFastQC: A Quality Control Tool for High Throughput Sequence Data, version V3FastQC: A Quality Control Tool for High Throughput Sequence Data.
-
Influenza virus propagation in embryonated chicken eggsJournal of Visualized Experiments 1:52421.https://doi.org/10.3791/52421
-
Influenza vaccination in children with asthma: randomized double-blind placebo-controlled trialAmerican Journal of Respiratory and Critical Care Medicine 169:488–493.https://doi.org/10.1164/rccm.200309-1251OC
-
ReportBBMap: A Fast, Accurate, Splice-Aware AlignerLawrence Berkeley National Laboratory.
-
Role of deficient type III interferon-lambda production in asthma exacerbationsNature Medicine 12:1023–1026.https://doi.org/10.1038/nm1462
-
Inflammatory responses to Pneumovirus infection in IFN-alpha beta R gene-deleted miceThe Journal of Immunology 175:4735–4744.https://doi.org/10.4049/jimmunol.175.7.4735
-
BookGenetic Data Analysis Using R: Introduction to PhylogeneticsTutorial using the R software.
-
DNA transfection: calcium phosphate methodMethods in Molecular Biology 1018:107–110.https://doi.org/10.1007/978-1-62703-444-9_10
-
Fast gapped-read alignment with bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
The sequence alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
Pediatric hospitalizations associated with 2009 pandemic influenza A (H1N1) in ArgentinaNew England Journal of Medicine 362:45–55.https://doi.org/10.1056/NEJMoa0907673
-
The burden of influenza illness in children with asthma and other chronic medical conditionsThe Journal of Pediatrics 137:856–864.https://doi.org/10.1067/mpd.2000.110445
-
Risk factors and outcomes among children admitted to hospital with pandemic H1N1 influenzaCanadian Medical Association Journal 182:39–44.https://doi.org/10.1503/cmaj.091724
-
Influenza A(H1N1)pdm09 virus and asthmaFrontiers in Microbiology 4:307.https://doi.org/10.3389/fmicb.2013.00307
-
A mechanism for priming and realignment during influenza A virus replicationJournal of Virology 92:e01773–17.https://doi.org/10.1128/JVI.01773-17
-
The immune profile associated with acute allergic asthma accelerates clearance of influenza virusImmunology & Cell Biology 92:449–459.https://doi.org/10.1038/icb.2013.113
-
Analyzing real-time PCR data by the comparative C(T) methodNature Protocols 3:1101–1108.https://doi.org/10.1038/nprot.2008.73
-
The FoldX web server: an online force fieldNucleic Acids Research 33:W382–W388.https://doi.org/10.1093/nar/gki387
-
Using bioluminescent imaging to investigate synergism between Streptococcus pneumoniae and influenza A virus in infant miceJournal of Visualized Experiments 1:2357.https://doi.org/10.3791/2357
-
Influenza-induced inflammation drives pneumococcal otitis mediaInfection and Immunity 81:645–652.https://doi.org/10.1128/IAI.01278-12
-
Asthma endotypes and an overview of targeted therapy for asthmaFrontiers in Medicine 4:158.https://doi.org/10.3389/fmed.2017.00158
-
The role of the priming loop in influenza A virus RNA synthesisNature Microbiology 1:16029.https://doi.org/10.1038/nmicrobiol.2016.29
-
Influenza virus RNA polymerase: insights into the mechanisms of viral RNA synthesisNature Reviews Microbiology 14:479–493.https://doi.org/10.1038/nrmicro.2016.87
-
Paucigranulocytic asthma: Uncoupling of airway obstruction from inflammationJournal of Allergy and Clinical Immunology 143:1287–1294.https://doi.org/10.1016/j.jaci.2018.06.008
-
Asthmatic bronchial epithelial cells have a deficient innate immune response to infection with rhinovirusJournal of Experimental Medicine 201:937–947.https://doi.org/10.1084/jem.20041901
Decision letter
-
Miles P DavenportSenior and Reviewing Editor; University of New South Wales, Australia
-
Cody AllisonReviewer; WEHI, Melbourne, Australia
-
Vijaykrishna DhanasekaranReviewer; Monash University, Australia
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
This manuscript investigates the potential for evolution of viruses in susceptible hosts. Using a mouse model of asthma, the authors show that asthmatic animals suffer from more severe infection, but also that the virus appears to show greater evidence of evolution in asthmatic animals (based on variation in the viral sequence). These changes in the virus support more vigorous viral replication in vitro, suggesting an increased virulence. Crucially, infection of healthy control mice with virus derived from the asthma-prone animals also demonstrated that viral strains emerging from asthmatic animals showed increased virulence when passed back to healthy animals.
This study shows that viral evolution in compromised hosts can lead to an increase in virulence that persists when transmitted back to healthy animals. The implications of this for viral spread and evolution at a population level are important, although clearly further work is required to extend this beyond the current animal model.
Decision letter after peer review:
Thank you for submitting your article "A paucigranulocytic asthma host environment promotes the emergence of virulent influenza viral variants" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by Miles Davenport as the Senior and Reviewing Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Cody Allison (Reviewer #1); Vijaykrishna Dhanasekaran (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
As the editors have judged that your manuscript is of interest, but as described below that additional experiments are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at eLife". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)
Summary:
This manuscript addresses an interesting angle of the concept that varying conditions of immunosuppression can facilitate the emergence of virulent mutant strains of Influenza. Data presented partially support this hypothesis, however there are concerns regarding the virulence of strains, deep sequence analysis, the downstream analysis of mutations, and the interpretation on influenza evolution. Therefore additional experiments are necessary to substantiate the statements in the title, Abstract and Discussion.
Essential revisions:
1) Authors state that more virulent strains emerge. While it is clear that mutational penetrance is higher in asthmatic mice, more evidence is required to establish it is more virulent.
a) weight loss and viral load are not enough to indicate virulence. This could simply be a result of the asthmatic phenotype, being unable to control virus. The increased replication in vitro is insufficient to ascribe to a mutant with higher replicative capacity.
b) what are the predominant mutations growing out of the cell lines following incubation with lung homogenates? Does it vary between the mice? If they're replicating faster they should be the predominant clones. Were the 3 asthmatic mice used for this experiment selected based upon the frequency of mutations identified in them?
c) should have histology and lung cytokines from asthmatic and non-asthmatic mice following challenge to corroborate increased virulence.
d) what happens to non-asthmatic mice when you challenge them with virus isolated from homogenised lung from the asthmatic group? Non-asthmatic animals should develop symptoms (weight loss, cytokines, oedema, infiltrates) if indeed you are generating strains of greater virulence.
2) Shannon's entropy scores in Figure 3 are non-parametric, with many asthmatic animals giving a score of zero. For PB1 (or any of the segments), can entropy scores be correlated with individual mice across the ~2-log spread of viraemia in this group (in Figure 2B) (i.e. does increased viral diversity correlate with high or low viraemia in those animals)?
3) There is concern with the use of a single measure of genetic diversity (Shannon entropy), which has been shown repeatedly to produce biased estimates (Please see. Zhao and Illingworth, Virus Evolution 5:vey041 https://doi.org/10.1093/ve/vey041 ). It is also not clear what the effect of the initial template load on measures of genetic diversity at each site/gene. Did the authors check if there was sufficient template diversity as a low diversity can bias all estimates? Additionally, it would be important to show any variation between samples / days as supporting information.
4) As before, since the rest of the paper relies heavily on single mutations (e.g. PB1 693), it would be important to use robust methods, or multiple methods to check whether they are consistent.
5) A major goal of this study was to address the importance of co-morbidities is on influenza evolution, but this hypothesis is not directly tested in this paper. Using C693Y as example, it is not clear what their prevalence is among naturally circulation H1N1pdm09 viruses. By filtering publicly available sequences, the authors have identified sequences from 51 asthma cases and compared these to 51 cases matching age and sex with no co-morbities (randomly I assume), but I feel these number are quite low to derive meaningful comparisons. To make inferences about the effect on influenza evolution, I suggest to map important mutations onto gene phylogenies to identify the dominance of mutations in naturally circulating populations. e.g. has Y693 or other PB1 mutations fixed in any of the circulation sub-lineages of H1N1.
It is also interesting this mutation is predicted to increase stability, and that they potentially show increased polymerase activity. If correct, shouldn't this mutation be dominant among naturally circulating strains?
6) As a part of the RdRp, mutations within PB1 gene by itself can increase mutational frequency, but this is not explored/discussed in light of results presented.
https://doi.org/10.7554/eLife.61803.sa1Author response
Essential revisions:
1) Authors state that more virulent strains emerge. While it is clear that mutational penetrance is higher in asthmatic mice, more evidence is required to establish it is more virulent.
a) weight loss and viral load are not enough to indicate virulence. This could simply be a result of the asthmatic phenotype, being unable to control virus. The increased replication in vitro is insufficient to ascribe to a mutant with higher replicative capacity.
We acknowledge that weight loss and viral load alone are not enough to indicate virulence. We have now updated the figure legend of Figure 2 to reflect a move away from virulence and instead focus on disease severity “RAGE-deficient asthmatic mice experience more severe influenza than non-asthmatic mice”. Furthermore, we have included day 4 viral titre data, histology scores of all time points, and cytokine and inflammatory marker data in the updated Figure 2. We have subsequently updated the following sections of the manuscript:
Materials and methods:
“Viral and host RNA was extracted from mouse lung samples, cDNA was synthesised, and real time PCR was performed as described previously (Hulme et al., 2020). […] Viral copy number was determined using IAV strain A/Puerto Rico/8/1934 H1N1 (PR8) virus matrix (M) gene cloned into pHW2000 plasmid and viral copy number was determined as described previously (Mancera Gracia et al., 2017).”
Materials and methods:
“Lungs were assessed for vascular changes, bronchitis, interstitial inflammation, alveolar inflammation, pneumocyte hypertrophy, and pleuritis (represented by a total score) by a veterinary pathologist who was blind to the study design.”
Results:
“While there was an increase in the histopathology score of mock infected asthmatic mice compared to non-asthmatic, there was no observable difference following IAV infection (Figure 2C). Furthermore, asthmatic mice showed a significantly impaired type I interferon response with a down regulation of both IFN-α and interferon inducible MX1 following IAV infection (Figure 2D).”
b) what are the predominant mutations growing out of the cell lines following incubation with lung homogenates? Does it vary between the mice? If they're replicating faster they should be the predominant clones. Were the 3 asthmatic mice used for this experiment selected based upon the frequency of mutations identified in them?
The question of mutations arising in vitro during the replication kinetics experiment is an interesting one. It is important to note that the asthmatic and non-asthmatic mice used for the replication kinetics were chosen at random to avoid any bias. In response to the reviewer’s suggestion, we have since deep sequenced the PB1 genome segment produced from the in vitro MDCK experiments and compared them back to their inoculum parent sequence to determine the predominant mutations above a 5% threshold.
One of the parent isolates was missing sufficient reads in regions of interest and was excluded from the analysis. Of the two remaining isolates, there was no change in the consensus sequence of the daughter population after in vitro replication compared the consensus sequence of the parent. In terms of minority viral populations, in one daughter population no viral variants (with non-synonymous mutations) were detected that constituted between 5-49% of the viral population. In the second daughter population, one non-synonymous mutation (PB1 I205M) was detected that constituted 28% of the viral population. However, this mutation was present in consensus sequence of the parent isolate, albeit at a higher percentage (95%). Together, these data suggest that the increased replication of viral isolates from asthmatic mice was not associated with a large number of additional in vitro adaptations. We have subsequently updated the manuscript to include these results the following section:
Results:
“Interestingly, comparing two of the viral samples pre and post in vitro replication in MDCKs we did not detect a change in the consensus sequence of the daughter population after in vitro replication compared the consensus sequence of the parent. In terms of minority viral populations, in one daughter population no viral variants (with non-synonymous mutations) were detected that constituted between 5-49% of the viral population. […] Together, these data suggest that the increased replication of viral isolates from asthmatic mice was not associated with a large number of additional in vitro adaptations”.
c) should have histology and lung cytokines from asthmatic and non-asthmatic mice following challenge to corroborate increased virulence.
Please see our response to reviewer query 1a above, we have now included histology and lung cytokine data.
d) what happens to non-asthmatic mice when you challenge them with virus isolated from homogenised lung from the asthmatic group? Non-asthmatic animals should develop symptoms (weight loss, cytokines, oedema, infiltrates) if indeed you are generating strains of greater virulence.
In response to the reviewer’s query, we have now performed experiments whereby we inoculated non-asthmatic C57BL/6 mice with lung homogenate from either asthmatic or non-asthmatic mice at a PFU of 2000 PFU per mouse (Figure 8). Interestingly but unsurprisingly, we saw that asthmatic host-derived IAV was associated with significantly more weight loss (Figure 8A), lower blood oxygen saturation (Figure 8B), and higher viral load (Figure 8C) than non-asthmatic host-derived IAV. Furthermore, mice infected with asthmatic host derived IAV had a significant upregulation in various anti-viral and pro-inflammatory cytokines (Figure 8E). Furthermore, upon infection of a non-asthmatic host with non-asthmatic derived IAV, there is an observable divergence in the PB1 haplotype differences as the virus continues to adapt to its host (Figure 9A). By contrast, upon infection of a non-asthmatic host with asthmatic derived IAV, there is a consolidation of the PB1 mutations that arose during the initial infection (Figure 9B). Together this data suggests that asthmatic hosts could facilitate the emergence of new strains that have increased pathogenicity.
We have subsequently updated the Materials and methods, Results and Discussion:
Materials and methods:
“Mouse strains
RAGE-deficient mice (Ager-/-) were kindly provided by Prof. Ann-Marie Schmidt, (New York University Langone Medical Centre, USA). […] Blood oxygen saturation of C57BL/6 was measured at 4 days post-infection using a collar sensor and Mouseox Plus pulse oximeter (Starr, Oakmont, PA, USA)”.
Materials and methods:
“Multiple dimension scaling (MDS) (LMDS V1.0 – R package) was utilised to visualise distance matrices to understand the dynamics of haplotype distance in a given quasi species population across donor and acceptor hosts (i.e. C57BL/6 mice infected with either 6 d.p.i murine lung homogenate from asthmatic or non-asthmatic host) [Shepard et al., 2016]. […] All scripts used for this study can be found at https://github.com/akaraw/Hulme_et_al and the workflow for haplotype reconstruction can be found in Supplementary file 3.”
Results:
“Asthmatic host-derived IAV exhibits greater disease severity and viral diversity upon infection of a wild-type host
Vulnerable patient groups may be a source for the emergence of new viral strains [Honce et al., 2020]. […] Together this data suggests that asthmatic hosts could facilitate the emergence of new strains that have increased pathogenicity”.
Discussion:
“Further to this, non-asthmatic host-derived viruses showed host adaption and little sign of morbidity in wild-type mice, while asthmatic-host derived virus resulted in increased morbidity. This suggests that viral variants born from an asthmatic host could (transmission permitting) lead to the emergence of more pathogenic viral strains”.
2) Shannon's entropy scores in Figure 3 are non-parametric, with many asthmatic animals giving a score of zero. For PB1 (or any of the segments), can entropy scores be correlated with individual mice across the ~2-log spread of viraemia in this group (in Figure 2B) (i.e. does increased viral diversity correlate with high or low viraemia in those animals)?
This is an interesting query and we thank the reviewer for suggesting it. To answer this, we assume the reviewer is referring to the viral load in the lung (PFU) when they are speaking of “viraemia” and we will answer accordingly. For the correlation analysis, we focused on the PB1 genome segment as this is the only gene that showed a difference between asthmatic and non-asthmatic in our results (Figure 3). On an individual animal level, we can see a positive relationship within the D6 asthmatic group, between the weight-loss (as a measure of disease severity) and both the Shannon Entropy (Author response image 1A), and the total number of mutations in the PB1 genome segment (Author response image 1B). However, no correlation between the viral load and Shannon Entropy (Author response image 1C), or the total number of mutations in the PB1 genome segment (Author response image 1D) was detected. Together this suggests that at day 6, asthmatic mice that experienced more severe influenza (as determined by weight-loss) produced more viral variants, but that this did not correlate with an increased viral load.
We are happy to include this data as supplementary material if the reviewers believe it adds value to the paper.

Asthmatic day 6 samples show a positive relationship between weight loss and measures of viral diversity in the PB1.
Influenza viral RNA isolated from the lung tissue of asthmatic (n=27) and non-asthmatic (n=25) mice was analysed for viral variants. Weight-loss as a function of (A) Shannon entropy of PB1 genome segment and of (B) SNVs detected in the PB genome segment. Viral titre (PFU/mL) as a function of (C) Shannon entropy of PB1 genome segment and of (D) SNVs detected in the PB genome segment. Samples with a Shannon Entropy value of 0 were excluded. A simple linear regression was fit, and statistical significance was determined using Prism. * p ≤ 0.05.
3) There is concern with the use of a single measure of genetic diversity (Shannon entropy), which has been shown repeatedly to produce biased estimates (Please see. Zhao and Illingworth, Virus Evolution 5:vey041 https://doi.org/10.1093/ve/vey041 ). It is also not clear what the effect of the initial template load on measures of genetic diversity at each site/gene. Did the authors check if there was sufficient template diversity as a low diversity can bias all estimates? Additional it would be important to show any variation between samples / days as supporting information.
We agree with the reviewers that one measure of genetic diversity is limiting. However, we selected this approach to be consistent with recently published works in this field, which have relied heavily on Shannon entropy without any other measures of diversity (please see Honce et al., 2020). In light of the reviewer’s concerns, we have now compared the read depth between samples and can conclude that no bias based on difference in read depth are present. We have now included the read depth results as Figure 3—figure supplement 1, and have included the use of a single measure of genetic diversity as a limitations of the study. These additions correspond to changes in the following sections:
Results: “To account for any bias from difference in sequencing coverage, the read count for each segment was analysed and there were no significant differences between groups (Figure 3—figure supplement 1)”.
Discussion:
“Additionally, only a single measure of genetic diversity was used (Shannon entropy), however we took into account any read depth bias, and within the field of influenza diversity this method alone has been routinely used [Honce et al., 2020; Abd Raman et al., 2020].”
We would also like to remind the reviewers that all the viral sequencing data generated in this study is available on the European Nucleotide Archive with accession number PRJEB40027 for anyone who wishes to look at our data further.
4) As before, since the rest of the paper relies heavily on single mutations (e.g. PB1 693), it would be important to use robust methods, or multiple methods to check whether they are consistent.
We would like to stress to the reviewers that our paper does not rely on any one mutation for our conclusion. Rather, we are suggesting that there is an observable increase in diversity of virus isolated from an asthmatic host, compared to a non-asthmatic host. Additionally, we have now demonstrated that virus derived from an asthmatic host can cause more severe disease when introduced to a non-asthmatic host (see response to reviewer query 1d above). Given these findings, additional studies monitoring individuals with asthma (and in particular those with non-allergic asthma) for increased influenza severity and the emergence of influenza virus variants may be of value.
The purpose of creating the C693Y mutation was to demonstrate that mutations found in key areas of the PB1 protein structure, exclusive to the asthmatic group, could affect (specifically, increase) the replication capability of the virus. However, to alleviate concerns that the paper relies heavily on a single mutation (C693Y), we have now performed mutagenesis and introduced the P651Q mutation into the PB1 plasmid and observed polymerase activity. These results confirm the crystal structure findings and indicate that the P651Q mutation increases polymerase activity compared to the wild type.
We have included the data in Figure 7 and a discussion of these points in the following sections:
Abstract:
“Specifically, C693Y and P651Q mutations in the PB1 of influenza virus was exclusively detected in mice with the cardinal features of paucigranulocytic asthma and these were associated with increased viral replication and polymerase activity”.
Results:
“PB1 C693Y and P651Q increases viral polymerase activity
We next sought to investigate the individual effect of the PB1 C693Y and P651Q mutations on viral polymerase activity using a minigenome assay, given its their predicted increase in stability. […] Together these data suggest that the C693Y and P651Q mutations identified exclusively in asthmatic hosts, is altering the interaction between PB1 and PB2 in a stabilising manner, resulting in increased polymerase activity”.
Discussion:
“Our structural analysis of PB1 P651Q predicted a stabilising effect, which was consistent with an increase in polymerase activity relative to the wild type”.
Materials and methods:
“The PB1-C693Y and PB1-P651Q mutations were independently as introduced into the PB1 gene segments of influenza A/Auckland/2009 (H1N1) using the Q5 Site-Directed Mutagenesis Kit from New England BioLabs as per the manufacturer’s instructions”.
5) A major goal of this study was to address the importance of co-morbidities is on influenza evolution, but this hypothesis is not directly tested in this paper. Using C693Y as example, it is not clear what their prevalence is among naturally circulation H1N1pdm09 viruses. By filtering publicly available sequences, the authors have identified sequences from 51 asthma cases and compared these to 51 cases matching age and sex with no co-morbities (randomly I assume), but I feel these number are quite low to derive meaningful comparisons. To make inferences about the effect on influenza evolution, I suggest to map important mutations onto gene phylogenies to identify the dominance of mutations in naturally circulating populations. e.g. has Y693 or other PB1 mutations fixed in any of the circulation sub-lineages of H1N1.
It is also interesting this mutation is predicted to increase stability, and that they potentially show increased polymerase activity. If correct, shouldn't this mutation be dominant among naturally circulating strains?
The reviewer is correct in the assumption that the 51 cases were randomly selected. We also agree with the reviewers that the sample size of 51 clinical cases is quite small. Unfortunately, at the time of manuscript preparation, there were only 51 clinical samples that listed asthma as a co-morbidity, thus limiting our analysis.
We agree with the reviewers that in order to make inferences about the effect on influenza evolution we would need to map the mutations onto phylogenies. However, as we used a mouse model, we do not expect to see these exact asthmatic mouse-specific mutations in the human population. Rather than suggesting that the mutations found within this study could arise in naturally circulating population, we are suggesting that asthmatic hosts give rise to more pathogenic mutations and using the clinical samples to highlight a similar trend.
6) As a part of the RdRp, mutations within PB1 gene by itself can increase mutational frequency, but this is not explored/discussed in light of results presented.
We thank the reviewer for this comment. Due to biosafety limitations, we cannot introduce the mutations that were identified into an infectious virus as this would be classified by our institution as gain-of-function work. As a result, we cannot experimentally investigate if the C693Y or P651Q mutations lead to increased mutation rates. However, previous studies have indeed shown that specific mutations in the PB1 gene segment can contribute to altering the fidelity of the RdRp (Lin, et al., 2019) and as such the emergence of mutations.
We have now included a discussion of this possibility in the Discussion:
“It is well established that SNVs can drastically change the replication fidelity and pathogenicity of a virus [Lin et al., 2019; Pauly et al., 2017]. […] Nevertheless, our results should be taken as proof of principal that an asthmatic host is more likely to produce more virulent viral variants – a finding that holds wide consequences for the emergence of new viral strains”.
https://doi.org/10.7554/eLife.61803.sa2Article and author information
Author details
Funding
Australian Research Council (DE180100512)
- Kirsty Renfree Short
Fonds Wetenschappelijk Onderzoek (12U7118N)
- Bram Vrancken
Fonds Wetenschappelijk Onderzoek
- Silvie Van den Hoecke
National Health and Medical Research Council (APP1130815)
- Kyle R Upton
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
SVDH was a Ph.D. fellow supported by FWO-Vlaanderen. B.V was supported by postdoctoral grant (12U7118N) of the FWO (Fonds Wetenschappelijk Onderzoek – Vlaanderen). KRU was supported by a NHMRC Career Development Fellowship (APP1130815). KRS was supported by an Australian Research Council DECRA (DE180100512). We would like to thank Professor Phillip Hugenholtz for his support in this project.
Ethics
Animal experimentation: All animal experiments were approved by the University of Queensland Animal Ethics Committee (permit no. 071/17).
Senior and Reviewing Editor
- Miles P Davenport, University of New South Wales, Australia
Reviewers
- Cody Allison, WEHI, Melbourne, Australia
- Vijaykrishna Dhanasekaran, Monash University, Australia
Version history
- Received: August 5, 2020
- Accepted: January 23, 2021
- Version of Record published: February 16, 2021 (version 1)
Copyright
© 2021, Hulme 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,071
- Page views
-
- 97
- Downloads
-
- 4
- Citations
Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.
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
-
- Genetics and Genomics
- Immunology and Inflammation
Age-associated DNA methylation in blood cells convey information on health status. However, the mechanisms that drive these changes in circulating cells and their relationships to gene regulation are unknown. We identified age-associated DNA methylation sites in six purified blood-borne immune cell types (naive B, naive CD4+ and CD8+ T cells, granulocytes, monocytes, and NK cells) collected from healthy individuals interspersed over a wide age range. Of the thousands of age-associated sites, only 350 sites were differentially methylated in the same direction in all cell types and validated in an independent longitudinal cohort. Genes close to age-associated hypomethylated sites were enriched for collagen biosynthesis and complement cascade pathways, while genes close to hypermethylated sites mapped to neuronal pathways. In silico analyses showed that in most cell types, the age-associated hypo- and hypermethylated sites were enriched for ARNT (HIF1β) and REST transcription factor (TF) motifs, respectively, which are both master regulators of hypoxia response. To conclude, despite spatial heterogeneity, there is a commonality in the putative regulatory role with respect to TF motifs and histone modifications at and around these sites. These features suggest that DNA methylation changes in healthy aging may be adaptive responses to fluctuations of oxygen availability.
-
- Immunology and Inflammation
- Neuroscience
Infection with Influenza A virus (IAV) causes the well-known symptoms of the flu, including fever, loss of appetite, and excessive sleepiness. These responses, mediated by the brain, will normally disappear once the virus is cleared from the system, but a severe respiratory virus infection may cause long-lasting neurological disturbances. These include encephalitis lethargica and narcolepsy. The mechanisms behind such long lasting changes are unknown. The hypothalamus is a central regulator of the homeostatic response during a viral challenge. To gain insight into the neuronal and non-neuronal molecular changes during an IAV infection, we intranasally infected mice with an H1N1 virus and extracted the brain at different time points. Using single-nucleus RNA sequencing (snRNA-seq) of the hypothalamus, we identify transcriptional effects in all identified cell populations. The snRNA-seq data showed the most pronounced transcriptional response at 3 days past infection, with a strong downregulation of genes across all cell types. General immune processes were mainly impacted in microglia, the brain resident immune cells, where we found increased numbers of cells expressing pro-inflammatory gene networks. In addition, we found that most neuronal cell populations downregulated genes contributing to the energy homeostasis in mitochondria and protein translation in the cytosol, indicating potential reduced cellular and neuronal activity. This might be a preventive mechanism in neuronal cells to avoid intracellular viral replication and attack by phagocytosing cells. The change of microglia gene activity suggest that this is complemented by a shift in microglia activity to provide increased surveillance of their surroundings.