|Home | About | Journals | Submit | Contact Us | Français|
The genotype of the host is one of several factors involved in the pathogenesis of an infectious disease and may be a key parameter in the epidemiology of highly pathogenic H5N1 influenza virus infection in humans. Gene polymorphisms may affect the viral replication rate or alter the host’s immune response to the virus. In humans, it is unclear which aspect dictates the severity of H5N1 virus disease. To identify the mechanism underlying differential responses to H5N1 virus infection in a genetically diverse population, we assessed the host responses and lung viral loads in 21 inbred mouse strains upon intranasal inoculation with A/Hong Kong/213/03 (H5N1). Resistant mouse strains survived large inocula while susceptible strains succumbed to infection with 1,000- to 10,000-fold-lower doses. Quantitative analysis of the viral load after inoculation with an intermediate dose found significant associations with lethality as early as 2 days postinoculation, earlier than any other disease indicator. The increased viral titers in the highly susceptible strains mediated a hyperinflamed environment, indicated by the distinct expression profiles and increased production of inflammatory mediators on day 3. Supporting the hypothesis that viral load rather than an inappropriate response to the virus was the key severity-determining factor, we performed quantitative real-time PCR measuring the cytokine/viral RNA ratio. No significant differences between susceptible and resistant mouse strains were detected, confirming that it is the host genetic component controlling viral load, and therefore replication dynamics, that is primarily responsible for a host’s susceptibility to a given H5N1 virus.
Highly pathogenic H5N1 influenza virus has circulated in Southeast Asia since 2003 but has been confirmed in relatively few individuals. It has been postulated that host genetic polymorphisms increase the susceptibility to infection and severe disease. The mechanisms and host proteins affected during severe disease are unknown. Inbred mouse strains vary considerably in their ability to resist H5N1 virus and were used to identify the primary mechanism determining disease severity. After inoculation with H5N1, resistant mouse strains had reduced amounts of virus in their lungs, which subsequently resulted in lower production of proinflammatory mediators and less pathology. We therefore conclude that the host genetic component controlling disease severity is primarily influencing viral replication. This is an important concept, as it emphasizes the need to limit virus replication through antiviral therapies and it shows that the hyperinflammatory environment is simply a reflection of more viral genetic material inducing a response.
Genetic polymorphisms within the genome of the infected host play an important role in the response to and outcome of microbial infections. In humans, several polymorphisms that affect HIV, hepatitis C virus (HCV), and herpes simplex virus (HSV) pathogenesis (1–6) have been described. However, very few data are available on how host genetic differences influence the course of influenza A virus infection, and what is available is at best inconclusive. Studies on highly pathogenic H5N1 influenza virus in Asia support the notion of host genetic polymorphisms affecting susceptibility to H5N1 virus (7, 8), with a high proportion of infected contacts being genetically related to the index case within family clusters. Also, during the 2009 H1N1 virus pandemic, individuals of indigenous descent in the Americas, Australia, and the Pacific had higher hospitalization rates than did the overall population (9, 10). Although the increased incidence of risk factors such as diabetes and asthma may have contributed to this difference, gene polymorphisms likely played an important role in susceptibility to influenza virus.
Despite the evidence favoring a role for host genetic polymorphisms in susceptibility to influenza virus infection, no human polymorphisms have been discovered. More importantly, we have yet to understand the mechanistic basis for the difference in pathogenesis: does the genetic change in susceptible hosts allow for increased replication dynamics of the virus or an altered immune environment, or does it completely change the hosts’ response to the pathogen independently of viral burden? Existing literature reporting large differences in the host responses after infection with low- or highly pathogenic viruses (11–17) would suggest that individuals with severe disease have an aberrant immune response to the virus. Alternatively, the genetic differences between the hosts allow for an increased virus replication rate, resulting in the induction of a more vigorous immune response. In humans, H5N1 virus pathogenesis is characterized by high viral load and increased production of proinflammatory cytokines culminating in severe disease and often death (i.e., 60% mortality rate in humans) (12, 15). Due to the correlative nature of the 2 parameters, it is unclear how they are linked and how each affects survival after infection. The importance of replication dynamics was recently highlighted in a paper by Hatta et al. demonstrating that the more virulent H5N1 virus also had the highest replication rate of the viruses studied (18).
To increase our understanding of H5N1 virus pathogenesis and to study the effects of host genetic variation, we initiated a large systematic effort to define the underlying mechanism involved in susceptibility or resistance to influenza virus infection. This knowledge is crucial for the identification of genes or gene networks associated with differences in pathogenesis. In summary, H5N1 virus pathogenesis in genetically diverse hosts is largely determined by the host factors controlling viral load in the lungs. Increased viral loads initiate the production of excessive amounts of proinflammatory cytokines, causing tissue damage and ultimately mortality of the susceptible host.
We experimentally inoculated 21 mouse strains with the highly pathogenic H5N1 influenza A virus A/Hong Kong/213/03 (HK213) and monitored the animals for 30 days thereafter for signs of morbidity and mortality. The 50% mouse lethal dose (MLD50) values varied from 40 50% egg infective doses (EID50) for the influenza virus-susceptible strain DBA/2S (susceptibility indicated by “S”) to more than 106 EID50 for the influenza virus-resistant strains BALB/cR and BALB/cByR (resistance indicated by “R”) (Fig. 1).
Several strains were resistant to the high-dose challenge and had MLD50 values between 105 and 106 EID50. Others were considered either susceptible (MLD50 < 104 EID50) or highly susceptible (MLD50 < 103 EID50). The diversity in pathogenic potential of a single highly pathogenic H5N1 virus prompted us to define a dose of virus that was lethal to the susceptible mouse strains but induced only mild morbidity without mortality in the resistant strains. Thus, we set the dose of 104 EID50 as the cutoff for H5N1 virus resistance/susceptibility and subsequently used it in all host response experiments.
An appropriate immune response mediated by hematopoietic cells is essential for clearance of influenza viruses and subsequent survival of the infected host. To test whether susceptible mouse strains have an inferior immune response due to one or more genetic polymorphisms, we generated chimeric DBA/2S mice whose bone marrow (BM) was replaced with that of C57BL/6R (DBA/2C57), BALB/cR (DBA/2BALB), or DBA/2S (DBA/2DBA/2) as a control. Twelve weeks post-BM transfer, successful engraftment was determined using strain-specific cell surface markers. The chimeric mice were subsequently infected with 1 of 2 doses of HK213 virus. At a dose of 104 EID50, the DBA/2C57 and DBA/2BALB mice succumbed to infection approximately 7 to 8 days postinoculation (dpi). This finding was similar to that seen in age-matched DBA/2 control mice and DBA/2DBA/2 mice (Fig. 2A). Intranasal inoculation with a 10-fold-lower dose (103 EID50) produced similar results: DBA/2BALB mice died around the same time (postinoculation day 8) as did the DBA/2DBA/2 control mice. DBA/2C57 mice survived a few days longer; however, all of the mice had died by day 12 postinoculation (Fig. 2B). Thus, the immune response of hematopoietic origin from resistant mouse strains was unable to rescue the phenotype of a susceptible strain, consistent with the hypothesis that the genetic polymorphism(s) affecting H5N1 virus pathogenesis in susceptible mouse strains promotes viral replication in the epithelial cells of the respiratory tract.
To identify differences in the molecular responses to HK213 virus infection in the susceptible and resistant mouse strains, we performed large-scale expression analysis of lung tissue before and at various time points after inoculation in 3 susceptible strains (DBA/2S, 129/SvImS, and A/JS) and 3 resistant strains (SMR, C57BL/6R, and BALB/cR). Probes (26,000) that were detected in one or more strains at any time point after infection were included in subsequent analyses. Prior to detailed gene and pathway expression analysis, we performed an unsupervised principal component analysis (PCA) to identify similarities or differences in overall gene expression between strains of mice before inoculation and 1, 3 and 7 days postinoculation (dpi). If all biological replicates of a particular condition (mouse strain) form a cluster and are located elsewhere on the plot compared to the rest of the data, that particular mouse strain is considered to have a distinct expression profile. Our goal was to identify similarities and differences among the strains and correlate those with our disease phenotype. Prior to infection or shortly thereafter (day 1 postinfection), no clustering was observed, and individual data were scattered throughout the plot (Fig. 3). At day 3 postinfection, the expression profiles of the mice clustered into 3 groups. Group 1 (red circle) contained only highly susceptible strains (DBA/2S and 129/SvImS). Group 2 (green circle) was composed of resistant strains (C57BL/6R and BALB/cR) only. Group 3 (blue circle) was composed of a susceptible strain and a resistant strain (A/JS and SMR, respectively) and was situated between groups 1 and 2 on the second component. By day 7, the clustering of susceptible (group 1) and resistant (group 2) strains remained intact. More importantly, the susceptible strain in group 3, A/JS, had joined group 1, and the resistant strain (SMR) was found in close proximity to the other resistant strains of group 2. These plots clearly indicate that a gene expression signature is associated with death or survival in these mice but also that these signatures are not necessarily fixed until later in the infectious process.
To identify qualitative and quantitative differences in the host responses associated with survival, we probed the transcriptomes of susceptible and resistant mouse strains before and after inoculation with H5N1 virus. Because our phenotype occurred after H5N1 virus infection, we focused on genes whose expression increased or decreased ≥2-fold compared to that in uninfected controls in at least two of the three susceptible or resistant mouse strains; this approach yielded a comprehensive set of 2,038 genes.
At day 1 postinoculation, 76 genes were differentially expressed, of which the majority (64) were shared by the susceptible and resistant strains (see Fig. S1A in the supplemental material). By day 3, the number of differentially expressed genes increased dramatically to 829, and 569 (68.6%) of those were shared by the two groups. Surprisingly, only 28 (3.4%) genes were uniquely up- or downregulated in the resistant strains and 232 (70.5%) were found only in the susceptible strains. These included several genes that encode cytokines associated with increased inflammation, e.g., Il28b, Ifna4, Ifna5, and Tnfrsf1b. Similar results were found at day 7, with 582 of 1,084 (53.7%) genes differentially expressed in the susceptible strains only. Overall, these data suggest that the resistant mouse strains do not express a unique set of genes or a pathway(s) controlling replication or disease. Rather, the difference in expression profiles identified by PCA (Fig. 3) was caused by altered or increased gene expression in the susceptible mouse strains.
Next, we performed pathway analysis of the up- or downregulated genes within each mouse strain. At day 1 postinoculation, gene sets in 5 of 6 strains were enriched for the glutathione pathway (SMR was the exception). By day 3, the identified pathways were similar among all 6 strains (see Fig. S1B in the supplemental material), with significant P values for all the well-known pathways, including the cytokine-cytokine receptor interaction pathway, Toll-like receptor pathway, NOD-like receptor pathway, and chemokine signaling pathway. The increased expression of many of these pathway-associated genes in all strains suggested that the difference in resistance to H5N1 virus pathology is not due to a defect in innate or antiviral immune signaling in the susceptible strains. On the contrary, our analysis of the total number of genes identified within a particular pathway (Fig. S1C) showed many more upregulated genes in the DBA/2S and 129/SvImS strains than in other strains. At day 7 postinfection, we found an enriched gene expression profile in the resistant mouse strains that is associated with systemic lupus erythematosus, cell adhesion molecules, the antigen processing and presentation pathway, and the graft-versus-host disease pathway. This profile is indicative of adequate virus-specific T- and B-cell responses in these mice. In contrast, the 3 susceptible strains continued to have many cytokine genes upregulated, as evidenced by the low P value for the cytokine-cytokine receptor interaction pathway (Fig. S1B and S1C; Table S1).
Finally, we performed statistical analysis on the array data to identify genes whose expression values correlated with the PCA clustering observed on days 3 and 7 postinoculation (see Table S2 in the supplemental material). The majority of the 85 identified genes belonged to proinflammatory pathways, including cytokine-cytokine receptor pathways, the Jak-Stat pathway, and innate sensing pathways like the Toll-like receptor pathway. To summarize, the difference between resistant and susceptible strains is predominantly quantitative, with the excessive production of proinflammatory cytokines in the susceptible strains most likely affecting downstream induction of robust adaptive immune responses. The overrepresentation of proinflammatory genes in the susceptible strains is indicative of ongoing viral replication and suggests higher replication rates and higher viral loads in these strains.
Many of the genes associated with severe disease that were identified by expression analysis are considered proinflammatory. Although inflammation is important for the induction of an innate and adaptive immune response, too much inflammation can be pathological and exacerbate disease. To validate the expression data and demonstrate increased production of proinflammatory mediators in the susceptible strains, we measured the production of CCL2, alpha interferon (IFN-α), IFN-β, tumor necrosis factor alpha (TNF-α), CXCL2, and CSF3 after inoculation with 104 EID50 in 3 susceptible mouse strains (DBA/2S, 129/SvImS, and A/JS) and 3 resistant strains (SMR, C57BL/6R, and BALB/cR). At day 3 postinfection, lung homogenates of susceptible strains contained significantly higher concentrations of proinflammatory cytokines than did those of most of the resistant strains (Fig. 4); however, intriguing patterns emerged. The concentration of CSF3 was significantly higher in all susceptible strains and in the SMR strain, which, based on expression analysis, clustered with the susceptible A/JS strain. A similar pattern was found for IFN-α and IFN-β. In contrast, the production levels of CCL2 and CXCL2 were significantly higher in the DBA/2S and 129/SvImS strains than in the other 4 strains. Finally, we found significantly higher levels of TNF-α in the susceptible strains than in the resistant strains. The production of certain cytokines, i.e., CCL2, TNF-α, and CXCL2, was significantly correlated with the MLD50 (P < 0.01; R2 = 0.77, 0.96, and 0.72, respectively). These results were similar to those identified earlier by array analysis.
We hypothesized that increased pathogenicity coupled with enhanced inflammation in the susceptible strains may be the result of higher viral titers. To test this hypothesis, we quantified the amount of infectious virus at days 2, 4, and 7 postinoculation with 104 EID50 of HK213 in the lungs of the 6 mouse strains used for expression and cytokine analysis. At day 2, the virus titer ranged from 104.9 50% tissue culture infective doses (TCID50)/ml in C57BL/6R and SMR strains to 106.5 TCID50/ml in the DBA/2S strain (Table 1). The virus titer differed significantly between susceptible and resistant strains (P < 0.05), and it correlated well with the MLD50 (R2 = 0.70, P < 0.05). In addition, we found significant associations between virus titer and MLD50 on day 4 (R2 = 0.82, P < 0.05) and day 7 (R2 = 0.81, P < 0.05). To estimate the viral burden over a period of 7 days, we calculated the area under the curve (AUC) for all 6 strains; again, we found a significant correlation with MLD50 (R2 = 0.87, P < 0.01).
To confirm the association between high viral load and increased susceptibility, we applied computer-assisted immunohistochemical analysis of lung tissue from 3 susceptible strains (DBA/2S, 129/SvImS, and A/JS) and 2 resistant strains (C57BL/6R and BALB/cR) to detect the frequency of influenza A virus nucleoprotein-positive (NP+) nuclei in cross sections of an entire formalin-fixed lung. Despite the small number of mice used in this experiment, it was clear that the susceptible strains had more NP+ nuclei than did the resistant strains on days 2, 4, and 7 postinfection (see Table S3 in the supplemental material).
To demonstrate the effect of increased viral load on proinflammatory cytokine production and survival, we inoculated DBA/2S and C57BL/6R mice with a 100-fold-higher dose (106 EID50). As predicted, the production of CCL2 and TNF-α 3 dpi was significantly higher (P < 0.01) (see Fig. S4 in the supplemental material) in the mice inoculated with 106 EID50 than in the mice inoculated with 104 EID50. More importantly, the amounts of CCL2 and TNF-α produced in the lethally infected C57BL/6R mice were similar to those found in DBA/2S mice infected with a much lower, but for this strain similarly lethal, dose of 104 EID50.
The newfound associations between viral loads, proinflammatory cytokine production, and outcome after infection were confirmed using quantitative real-time PCR on a select set of proinflammatory mediators in conjunction with influenza A virus matrix RNA species. Twenty-four hours postinoculation, the amount of influenza virus RNA was significantly higher (analysis of variance [ANOVA], P < 0.001) in the three susceptible strains DBA/2S, 129SvImS, and A/JS than in the two most resistant strains, C57BL/6R and BALB/cR (Fig. 5A). Not surprisingly the reduced viral load in the resistant strains minimized the production of proinflammatory cytokines on day 3 (Fig. 4) and allowed these two mouse strains to cluster at the transcriptional level (Fig. 3), displaying a “normal” inflammatory response to the invading pathogen. The third resistant strain, SMR, had a high viral RNA load at 1 dpi similar to those of the susceptible strains and therefore clustered with a susceptible (A/JS) strain at the transcriptional level and also produced more proinflammatory mediators. Three days postinoculation, the resistant C57BL/6R and BALB/cR strains continued to have smaller amounts of viral RNA (ANOVA, P < 0.0001) than did the susceptible strains. In addition, the viral load in the third resistant strain, SMR, was lower than those in all three susceptible strains, providing further support for our hypothesis that viral load dictates the pathogenic response to infection in genetically diverse mice. Interestingly, the increase in amount of viral RNA between 1 and 3 dpi is significantly lower (P < 0.01) in the resistant strains than in the susceptible strains (Fig. 5B), confirming that lower viral load and early control of viral replication are key to minimizing disease and promoting survival.
To exclude the possibility that excessive production of cytokines promotes virus replication and therefore increases the viral load, we determined the mRNA levels of several proinflammatory mediators after inoculation with H5N1 influenza A virus. On day 1 postinoculation, the increases in expression of ccl2, ifnb1, and tnfa (threshold cycle [CT] value) were similar between all 6 mouse strains (Fig. 5C; see also Fig. S3 in the supplemental material), indicating that the production of proinflammatory cytokines is not responsible for the observed difference in viral RNA at 3 dpi. At 3 dpi, the fold increase in mRNA associated with ccl2, ifnb1, and tnfa was significantly lower (P < 0.001) in the resistant strains C57BL/6R and BALB/cR than in the three susceptible strains plus SM/JR, providing further support that early viral load determines the inflammatory response and as such the outcome after infection.
To demonstrate that the genetic differences do not affect the production of proinflammatory cytokines per unit of viral RNA, we computed the ratio between cytokine RNA and viral RNA CT values (see Fig. S3 in the supplemental material). At 1 and 3 dpi, the ratios for IFN-β1 were similar across all strains tested. For the other cytokines, we observed significant differences in the expression ratios among the different mouse strains; however, limited evidence for an association of these differences with disease severity or susceptibility to H5N1 virus infection was found. Combined, these analyses suggest not that the increased susceptibility is caused by an intrinsic difference in the response to infection but rather that it is caused by a higher viral load inducing excessive production of proinflammatory cytokines and hence increased morbidity and mortality in the susceptible mouse strains.
Polymorphisms in the genome of the host play an important role in the severity of infection with a highly pathogenic H5N1 influenza A virus. The mechanism dictating the difference in H5N1 pathogenicities among genetically diverse hosts is, however, unknown. Detailed analysis of 6 H5N1 virus-infected inbred mouse strains demonstrated that the viral load is responsible for the heightened inflammation in susceptible hosts. This produces a pathogenic environment as early as 3 days postinoculation, ultimately causing the host to succumb.
Several studies have compared high- and low-pathogenic influenza A viruses in various host species and, without exception, identified a strong association between pathogenicity and excessive production of proinflammatory mediators (16, 18–25). Two possible mechanisms explain this difference: first, increased replication dynamics of the more pathogenic viruses may result in the infection of more epithelial cells and subsequent production of higher levels of cytokines. Second, altered intracellular signaling properties of viral proteins (e.g., PB1-F2, NS1, or hemagglutinin [HA]) of the pathogenic virus may trigger cells to produce more or different cytokines. Although both of these mechanisms are likely involved when comparing subtypes of influenza A viruses in a single homogenous host species, it is not well understood how a single highly pathogenic H5N1 virus produces such varied pathologies in genetically diverse hosts. Does the genetic change in susceptible hosts allow for increased replication dynamics or an altered immune environment, or does it completely change the hosts’ response to the pathogen independently of viral burden?
Based on the data presented, we conclude that the susceptible host provides a milieu that allows for increased viral replication, which subsequently induces the production of more proinflammatory cytokines and tissue damage and ultimately results in the death of the animal. Although this study is the first of its kind, detailed analysis of H5N1 virus-infected humans in Vietnam also demonstrated a significant correlation between viral load, cytokine concentration, and severity of disease (12), thereby supporting the conclusion that differences in pathogenesis within a single host species are the result of variable replication dynamics.
The difference in viral load between susceptible and resistant strains was presented very early after inoculation (24 to 48 h), suggesting that the virus replication rate in the susceptible mouse strains is higher than that in the resistant strains. At this point, it is unclear whether this increase in replication rate is due to an inefficient antiviral mechanism or to an increase in cell metabolism and/or the availability of nutrients required for replication. The observation that the viral load increased nearly 8-fold in the susceptible strains between 1 and 3 dpi compared to <4-fold in the resistant strains suggests that the increased replication rate is caused by one or more defective antiviral effector molecules induced by the production of type I and II interferons shortly after infection. Alternatively, there could be a difference in infectious dose (50% infectious dose [ID50]) between the resistant and susceptible mouse strains due to variations such as the sialic acid receptor content of the respiratory tract. Previously, we had shown a 100% infection rate in a resistant strain and a susceptible strain following inoculation with a very low dose of 100 EID50 of HK213, suggesting that the ID50s are similar between resistant and susceptible mice and that this is unlikely the reason for the differences observed.
The importance of replication kinetics or polymerase activity in pathogenicity is well established (18, 26–30). The PB2E627K and PB2D701N mutations are associated with increased polymerase activity and faster replication; thus, viruses containing these mutations are often more pathogenic in mice (18, 27, 31). Molecular characterization of a variant PR8 virus adapted to induce severe disease in Mx1-positive mice found several mutations in the polymerase complex that allowed the virus to grow faster and therefore be more pathogenic than its wild-type counterpart (32). A recent comparison of A/Hong Kong/483/97 and A/Hong Kong/486/97 H5N1 viruses also surmised that replication rate regulates disease severity and the magnitude of the virus-specific CD8+ T-cell response (18). After low-dose infection, the HK483 virus harboring a PB2 protein with a lysine residue at position 627 grew to significantly higher titers early during infection than did the variant containing a glutamic acid at this position. This difference in viral loads resulted in quantitative and qualitative differences in the cytotoxic T-cell responses that could be reversed upon administration of oseltamivir, a drug that limits viral replication. Combined, these studies suggest that the replication rate of the virus and hence viral load are responsible for the increase in inflammation and therefore pathogenesis of the virus. Reducing viral load or replication rate will most likely reduce morbidity and mortality upon infection.
As a result of the heightened viral load during the early stages of infection, the immune response in susceptible mice is excessive and can be detected within 3 days of infection. This relatively brief period of time has implications for potential intervention strategies for severe cases of influenza virus infection. It also explains the relatively short time period during which current commercially available antiviral therapy has proven effective. Similar observations have been made with monoclonal antibodies, whose treatment window is 3 to 4 days postinfection (33, 34). The opportunity to treat susceptible mice is much shorter (less than 1 day postinfection) than that for resistant strains of mice (3 days postinfection [unpublished data]). In our study, one mouse strain that displayed an intermediate pathogenic profile (SMR) recovered from the infection. This finding indicates that with the right treatment a pathogenic profile can be reversed and the animal can survive.
Clues about what biological trait is important for survival of influenza virus infection may be found in the A/JS and SMR strains. Although the strains had similar pathogenic expression profiles, the SMR mice survived the infection while the A/JS mice did not. The absence of glutathione metabolism pathway-associated gene expression in the SMR mice may provide us with a mechanism for the difference in pathogenesis. Glutathione metabolism gene transcription is initiated early after infection to reduce harmful reactive oxygen species (ROS). Mice lacking a functional NADPH oxidase (Cybb-deficient mice) can no longer produce ROS, and intranasal infection with influenza virus was shown to increase the production of proinflammatory cytokines and promote cellular infiltration into the lungs (35, 36). Despite the increase in inflammation, Cybb-deficient mice had lower viral loads and improved resolution of the infection, similar to what we observed in SMR mice. As such, further studies to determine the production of ROS and proinflammatory cytokines and to measure the induction of the glutathione metabolism response pathway in activated macrophages of SMR and C57BL/6R mice would be of interest. A lack of ROS production, minimal expression of glutathione response genes, and reduced proinflammatory cytokine production in the SM mouse macrophages will form the basis for future studies into the role of ROS in influenza virus pathogenesis. An alternative explanation is the presence of one or more genetic changes affecting viral clearance. An example of this is the hemolytic complement (Hc) gene, which affects T-cell responses and therefore viral load at later time points during infection. The SMR strain expresses a fully functional Hc gene; thus, viral load is reduced, inflammation is limited, and the animals are more likely to survive infection than are A/JS mice, which do not express a fully functional Hc gene.
Although the susceptible strains and SMR produced large amounts of proinflammatory cytokines (e.g., IFN-β and CSF3), certain cytokines were predominantly produced in the extremely susceptible strains such as DBA/2S and 129/SvImS. The excessive production of certain proinflammatory mediators such as CSF3 and type I interferons may cause a mild, controllable form of inflammation, though the addition of a second signal such as CCL2 and/or TNF-α may create an uncontrollable cascade of events resulting in tissue damage and diminished adaptive immune responses in these mice. Results from in vitro experiments suggest that IFN-β and TNF-α act synergistically to induce a response quantitatively and qualitatively different from that of either cytokine alone (37, 38). Also, TNF-α affects morbidity after H5N1 virus infection (17, 39).
While the current study focused on identifying the mechanism of severe disease, e.g., virus replication, there is obvious merit in identifying the specific genetic polymorphisms responsible for this phenotype. The MLD50 curve for highly pathogenic H5N1 virus (Fig. 1) suggests that several genetic polymorphisms are involved. A similar complexity was previously identified with the recombinant inbred BXD strain family, in which 3 genetic loci contribute to the response to influenza virus infection (40). To identify the responsible polymorphisms, genetic tools like the BXD RI strains, consomic strain sets, or the Collaborative Cross will be extremely useful (41–44). Several of the candidate genes identified in our earlier BXD study were evaluated for their association with disease severity in our larger mouse panel using the available single nucleotide polymorphism (SNP) data. Hemolytic complement, on Qivr2, did not differentiate the resistant and susceptible mouse strains. Qivr7 contains several potential antiviral genes (Trim genes); however, there are not enough public SNP data available to do the analysis. The candidate genes on Qivr11 (Grn, Ifi35, and Dhx58) all associated significantly (P < 0.01) with disease severity and are currently under investigation for their role in influenza virus pathogenesis. Finally, none of the candidate genes on Qivr17 (Xdh, Eif2ak2, Emilin2, and Nlrc4) separated with our phenotype.
In conclusion, gene polymorphisms in the genome of the infected host can have a profound impact on the course of an H5N1 influenza virus infection. Although featured as too much inflammation with excessive production of proinflammatory cytokines, the pathogenicity of influenza virus infection is rooted in increased virus titers in the lungs of the susceptible hosts.
Female mice (6 to 8 weeks old) from the following strains were purchased from Jackson Laboratories (Bar Harbor, ME) and housed in the Animal Resource Center at St. Jude Children’s Research Hospital: C57BL/6J (stock no. 664), BALB/cJ (stock no. 651), SM/J (stock no. 687), DBA/2J (stock no. 671), FVB/NJ (stock no. 1800), LP/J (stock no. 676), NZB/BlnJ (stock no. 684), C3H/HeJ (stock no. 659), BALB/cByJ (stock no. 1026), AKR/J (stock no. 648), B10.D2-H2/oSnJ (stock no. 461), A/J (stock no. 646), A/HeJ (stock no. 645), 129/SvImJ (stock no. 2448), MRL/MpJ (stock no. 486), NZW/LacJ (stock no. 1058), BTBR T+ tf/J (stock no. 2282), KK/HlJ (stock no. 2106), 129X1/SvJ (stock no. 691), DBA/1J (stock no. 670), and NOD/ShiLtJ (stock no. 1976). Strain selection was based on representation in large mouse genetic databases like Haplotype Mapping by Roche Pharmaceuticals, Ancestry Mapping by Perlegen, and full-genome sequencing by the Sanger Institute (45–47). The mice received water and food ad libitum, and all experimental procedures were approved by the Animal Care and Use Committee of St. Jude Children’s Research Hospital.
A reverse genetics variant of A/Hong Kong/213/03, a highly pathogenic H5N1 influenza A virus, was propagated in 10-day-old embryonated chicken eggs. The allantoic fluid containing the infectious virus was harvested, and the infectious virus titer was determined in eggs. The A/Hong Kong/213/03 virus (H5N1) variant contains 7 segments of A/Hong Kong/213/03 and 1 segment of A/Chicken/Hong Kong/52/03 H5N1 virus (40). This virus is referred to as HK213 throughout the paper.
Experimental inoculation of mice with HK213 virus was done intranasally in 30 µl phosphate-buffered saline (PBS) as described previously (40, 48). MLD50 experiments were repeated, and the reported results were calculated from the cumulative results.
CD4+ and CD8+ cells were depleted in vivo after intraperitoneal (i.p.) injection of DBA/2S, C57BL/6R, and BALB/cR donor mice with anti-CD4 and anti-CD8 antibodies. BM was then collected from the femurs of these mice, and approximately 6 million cells were injected intravenously into 9-Gy-irradiated DBA/2S mice. The mice were treated with Sulfatrim for 3.5 weeks to prevent any opportunistic infection.
At 12 weeks postengraftment, a small volume of peripheral blood was collected and subjected to flow cytometric analysis to test for the expression of cell surface markers of the respective donor mouse strains. In the C57BL/6R BM-recipient mice, we measured the percentage of H-2b+ (C57BL/6) and H-2k+ (DBA/2) cells within the CD3+ T-cell population. In the BALB/cR BM-recipient mice, we measured the percentage of CD22.1+ (DBA/2) and CD22.2+ (BALB/c) cells within the CD19+ T-cell population. Animals in which more than 90% of the T or B cells were donor-derived cells were considered successful BM chimeras.
Following intranasal inoculation with 104 EID50 of HK213 virus, lungs were harvested on days 2, 4, and 7 postinfection, and virus titers were determined on Madin-Darby canine kidney cells (40). At each time point, 4 to 6 animals from each mouse strain were used. The average virus titer was calculated, and the AUC was calculated between days 0 and 7 postinfection.
Immunohistochemical analysis was performed on formalin-fixed lung tissue harvested from HK213 virus-infected or uninfected DBA/2S, 129/SvImS, A/JS, C57BL/6R, and BALB/cR mice, as described previously (40). Digital images were obtained with a ScanScope (Aperio, Vista, CA), and the percentage of NP-positive nuclei was determined with ImageScope (Aperio, Vista, CA) using a nucleus-based algorithm.
Lungs from 6 strains of mice inoculated with 104 EID50 of HK213 virus were harvested on day 3 postinfection and homogenized in 1.0 ml PBS for two 30-s intervals at 30 Hz (TissueLyser II; Qiagen, Valencia, CA). After centrifugation for 30 s at 16,000 × g, the supernatant was collected, divided into aliquots, and stored at −80°C. The homogenates were subsequently thawed to quantify the amounts of CCL2, TNF-α, CSF3, and CXCL2 in the lungs of infected animals by using Quantikine enzyme-linked immunosorbent assay (ELISA) kits (R&D Systems, Minneapolis, MN). The amounts of IFN-α and -β in the lung were determined with the respective ELISA kits (PBL Laboratories, Piscataway, NJ). For each cytokine at a given time point, 4 to 6 animals were tested, and the average concentration ± standard error of the mean (SEM) is reported.
Six mouse strains, DBA/2S, 129SvImS, A/JS, SMR, C57BL/6R, and BALB/cR, were selected for an extensive RNA expression analysis before and after inoculation with 104 EID50 of HK213 virus. Lungs were isolated in 2 batches from uninfected mice and infected mice at days 1, 3, and 7 postinfection. RNA was extracted using TRIzol (Invitrogen, Carlsbad, CA), as previously reported (40) and submitted to a DNA cleanup protocol (Qiagen). Next, microarray analysis was done on RNA obtained from DBA/2S, 129SvImS, A/JS, SMR, C57BL/6R, and BALB/cR mice, using Illumina MouseWG-6 v1.1 Expression BeadChips (Illumina, San Diego, CA). PCA analysis identified a large batch effect; however, both batches included uninfected samples, which allowed us to calculate fold differences in gene expression. Also, batches were done according to the day postinfection, and all 6 strains were assessed simultaneously for each day.
Before analysis, probes without a signal (i.e., detection P value of >0.05) in any of the 104 samples were removed from the data set. Signal values were ln-start transformed, ln(signal + 20), to stabilize variance across the range of intensities and approximate normalities in preparation for parametric tests. Using Partek Genomics Suite 6.3, we visualized the structure of the data via PCA. As a result of the complexity of the disease, we compared expression patterns of genes up- or downregulated in most of the susceptible or resistant strains. For a gene to qualify for an individual strain, we required a 2-fold increase or decrease in expression in the majority of the individual samples for that strain on that day. Once a gene qualified for a particular strain, we next assessed whether that gene was shared among the majority of the resistant or susceptible (2/3 or 3/4) mouse strains. Pathway analysis was performed using the online bioinformatics tool DAVID (49).
Samples of cDNA were prepared from lung tissue RNA (200 ng) by using Superscript III (Invitrogen) and random hexamers and used for quantitative real-time PCR. The CT values for murine interleukin-6 (IL-6) (Mm00446190_m1), CCL2 (Mm00441242_m1), CCL4 (Mm00443111_m1), IFN-β1 (Mm00439552_s1), TNF-α (Mm99999068_m1), CSF3 (Mm00438334_m1), and β-actin were determined using commercially available primer-probe pairs from Applied Biosystems (Foster City, CA). To quantify the amount of virus in each RNA sample, we used the cDNA prepared with random hexamers in combination with a primer probe mix that was specific for the matrix gene (50). Baseline levels of each cytokine (ΔCT value) were determined after normalization against β-actin, and the increase in expression (ΔΔCT value) was calculated. No RNA species for IFN-β1, CSF3, or influenza virus matrix RNA were detected prior to infection. To calculate the increase, we used a CT value of 40 for the uninfected samples. RNA samples (n = 3 to 4) were tested for each strain and time point prior to or after infection with HK213 virus.
Statistical analysis of differences in lung virus titers was determined by a Student t test following ln transformation of the data. Cytokine and chemokine production and quantitative PCR results were analyzed using a Student t test or ANOVA when comparing more than two groups. For Fig. 4, statistical significance between the six mouse strains is represented by letters above each column, with different letters signifying distinct statistical groups. Linear regression analyses were done to analyze the viral load (log10 TCID50/ml and ΔΔCT value) with MLD50 (log10 EID50). P values less than 0.05 were considered significant.
Pathway analysis on differentially expressed genes in 6 mouse strains. (A) Venn diagram displaying differentially expressed genes that are shared (overlap areas) or unique to the majority (2/3) of the susceptible (red) or resistant (green) mouse strains. The genes showing at least a 2-fold up- or downregulation in expression are included. (B) Differentially expressed genes were tested for enrichment of specific pathways. The P values for particular pathways are indicated by different colors (white, P > 0.01; yellow, P < 10−2; orange, P < 10−4; red, P < 10−6; dark red, P < 10−8). (C) The total numbers of differentially expressed genes within the cytokine-cytokine receptor, Toll-like receptor, and Jak-Stat pathways are shown for each strain on days 3 and 7 postinoculation. For a complete list of genes and their fold differences in expression from those of uninfected controls, see Table S1. The color scheme is a representation for the number of differentially expressed genes within a pathway with above-average values (>average + 1× SEM) in red and low values (<average − 1× SEM) in yellow. Intermediate or average values are indicated by orange. The total numbers of genes in each individual KEGG pathway are 244, 99, and 152 for the cytokine-cytokine receptor, Toll-like receptor, and Jak-Stat pathways, respectively. Download Figure S1, PDF file, 0.1 MB.
Increased expression of proinflammatory cytokines in susceptible mouse strains. Quantitative real-time PCR was used to measure the expression of multiple cytokines and chemokines on days 1, 3, and 7 after inoculation with 104 EID50 of A/Hong Kong/213/03 virus. The bars indicate the fold increases (ΔΔCT value) over uninfected samples. Red and green bars indicate susceptible and resistant mouse strains, respectively. The mouse order from top to bottom is BALB/cR, C57BL/6R, SMR, A/JS, 129/SvImS, and DBA/2S. Download Figure S2, PDF file, 0.1 MB.
Higher viral load in susceptible strains increases the production of proinflammatory mediators. The ratio between viral and inflammatory RNAs (ΔΔCT value) is determined for 1 and 3 days postinoculation in three susceptible strains (DBA/2S, 129/SvImS, and A/JS) and three resistant mouse strains (SMR, C57BL/6R, and BALB/cR). The similarity in RNA ratios among the 6 strains suggests that the excessive production of cytokines and chemokines is merely a reflection of increased viral loads in the susceptible mouse strains. Download Figure S3, PDF file, 0.1 MB.
Infection with larger doses of H5N1 virus results in increased production of TNF-α and CCL2 in susceptible and resistant mouse strains. Groups of DBA/2S and C57BL/6R mice were inoculated with 104 (gray bars) or 106 EID50 (black bars) of A/Hong Kong/213/03 virus, and 3 dpi the production of TNF-α and CCL2 in the lungs of these animals was quantified by ELISA. Bars represent the average production levels of TNF-α or CCL2 (pg/ml) + SEMs. Download Figure S4, PDF file, 0.1 MB.
List of genes within the cytokine-cytokine receptor pathway, Jak-Stat pathway, and Toll-like receptor pathway (KEGG) and their respective fold differences (log2 values) on days 3 and 7 after inoculation with highly pathogenic H5N1 influenza virus.
List of genes whose expression values are significantly different between the three groups (DBA/2S + 129SvImS, A/JS + SMR, and C57BL/6R + BALB/cR) on day 3 and resistant and susceptible strains on day 7 postinoculation with H5N1 influenza A virus. Genes whose expression values were significantly different prior to inoculation (baseline) were excluded.
Frequency of influenza A virus nucleoprotein-positive cells in a cross section of H5N1 virus-infected lung tissue of different inbred mouse strains.
We thank Stacey Schultz-Cherry and Hui-Ling Yen for critically reviewing the manuscript and Paul Thomas for the anti-CD4 and -CD8 antibodies. We thank David Carey and Scott Krauss for their help in the ABSL3+ facility.
This project was funded, in part, by grants from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under contract no. HHSN266200700005C and by the American Lebanese Syrian Associated Charities (ALSAC).
Citation Boon ACM, et al. 2011. H5N1 influenza virus pathogenesis in genetically diverse mice is mediated at the level of viral load. mBio 2(5):e00171-11. doi:10.1128/mBio.00171-11.