To study the virulence and pathogenesis of the 1918 influenza virus, mice were infected intranasally with a contemporary human influenza A/Texas/36/91 H1N1 virus (Tx91), Tx91 virus containing the 1918 HA and NA genes (2:6 1918), the 1918 HA, NA, M, NP and NS genes (5:3 1918), or the reconstructed 1918 (H1N1) recombinant virus (r1918). Mice infected with r1918 showed the most severe weight loss and the earliest mean time of death of all viruses tested (Supplementary Information). Moreover, titration of virus recovered from lung tissue showed that mice infected with r1918 had about tenfold higher levels of virus on days 1 and 3 after infection than those infected with 2:6 1918 or 5:3 1918. By day 5 after infection there was a slight decrease in the lung titres of r1918-infected mice, and mice infected with 2:6 1918, 5:3 1918 and r1918 showed similar lung titres of about 107
50% embryo-infective doses per ml. Our results were consistent with those in previous reports3
: most mice infected with r1918 were dead by day 5 after infection.
Previous studies showed increased lung pathology in mice infected with influenza viruses bearing genes from the 1918 virus3–5
. A more detailed histological analysis of mice infected with 5:3 1918, 2:6 1918 or r1918 viruses revealed differences in pathology during the course of infection. Whereas mice infected with Tx91 reassortant virus had minimal to mild diffuse histiocytic alveolitis without neutrophils or lesions in the bronchi or bronchioles (), mice infected with 5:3 1918, 2:6 1918 or r1918 viruses had pulmonary lesions that consisted of mild to severe necrotizing bronchitis with accompanying moderate to severe alveolitis and associated moderate to severe alveolar oedema (). However, the character of the alveolitis varied between the 1918 viruses and the sampling periods. At day 3, the r1918-infected mice had acute moderate to severe diffuse alveolitis with a prominent neutrophilic component and severe alveolar oedema (). For some lungs, alveolitis was multifocal with intense neutrophilic infiltration that distended the alveoli (). By contrast, mice infected with 5:3 1918 and 2:6 1918 viruses had moderate to severe alveolitis, but histiocytes predominated over neutrophils and the alveolar oedema was only moderately severe (). At day 5 the pulmonary response in r1918-infected mice had changed to subacute and was more severe with diffuse alveolitis retaining the severe oedema, but changing to a predominance of histiocytes with fewer neutrophils (). The bronchial and bronchiolar lesions were similar to those at day 3. However, in some r1918-infected lungs, an acute inflammatory response continued, as was evident from multiple scattered foci of intense neutrophilic alveolitis. In contrast, some mice infected with 5:3 1918 and 2:6 1918 viruses showed slightly less alveolitis, alveolar oedema, bronchitis and peribronchitis on day 5 than on day 3. These results show that the pathology induced by r1918 was distinct from that observed in mice infected with 5:3 1918 and 2:6 1918 viruses, indicating that cooperation between the individual 1918 viral proteins might be required for the induction of maximum pathology.
Photomicrographs of lung sections from mice infected with influenza A virus, stained with haematoxylin and eosin
To understand why 1918 influenza virus showed extreme virulence in mice, we used gene expression analysis to provide a global view of the host response in lungs of infected mice. Expression microarray analysis was performed by taking equal masses of total RNA isolated from lungs of mice infected with strain Tx91, 2:5 1918, 5:3 1918 or r1918 and comparing them with a common reference sample prepared from lungs of mock-infected mice. Lungs from three mice were harvested for each virus on days 1, 3 and 5 after infection (total n
= 36). Initially, our analysis was directed to genes associated with the activation of key immune cells known to be important for the response to influenza virus infection, including T helper 1 (TH
1), natural killer, macrophage and neutrophil cells (shown in ). This analysis revealed a significant increase in the expression of genes associated with these immune cell populations in lungs of mice infected with r1918. Expression patterns of cytokines, chemokines and apoptosis-related genes showed the most significant and earliest activation during r1918 infection (Supplementary Information). The increased expression of immune-cell-related genes was correlated with virulence, because strain r1918 induced the most significant and earliest expression of immune-response-related genes, followed in magnitude and occurrence by 2:6 1918, 5:3 1918 and Tx91. Because many innate response-related genes showed saturation of the fluorescence signal by expression microarray, we measured the messenger RNA and protein expression of several key cytokines, including tumour necrosis factor (TNF), by quantitative real-time polymerase chain reaction (PCR) and enzyme-linked immunosorbent assay (ELISA). This analysis demonstrated that r1918 infections consistently resulted in a higher and prolonged expression of these mRNAs and proteins (Supplementary Information). It was shown recently6
that increased production of cytokines and chemokines during infection with 2:6 1918 could be attributed to the presence of neutrophils and alveolar macrophages in the lungs of infected mice. Our studies showed that r1918 infections also resulted in increased neutrophil and macrophage recruitment, accompanied by a significant increase in the expression of cytokine and chemokine genes. Expression of genes related to T helper 2 (TH
2) cell and nitric oxide metabolism showed similar patterns in all viral infections on days 1, 3 and 5, indicating that these responses might not have been primary mediators of r1918 virulence (Supplementary Information). The cellular origins of the genes identified by microarray are unknown, but taken together these data show that r1918 produced the most severe and rapid disease in mice, which could be correlated with an increased expression of inflammatory response genes similar to that reported for highly pathogenic avian influenza viruses7,8
Increased and earlier expression of genes associated with the activation of key immune cells in mouse lung infected with r1918 influenza virus
To characterize the functional consequences of gene expression changes associated with infection with r1918, we performed pathway analysis of the gene expression data with Ingenuity Pathways Analysis (see Methods). As shown in , pathway analysis of the expression array data revealed that infection with r1918 virus resulted in the most significant activation of death receptor, interleukin-6, type I interferon and Toll-like receptor responses. These data show that the aggregate host response to r1918 virus infection results in an increase in pro-inflammatory and cell death responses. We observed an inverse correlation between influenza virulence and G2/M checkpoint control and glutathione metabolism. The increased expression of cell-cycle and glutathione metabolism genes in Tx91 infections could be important components in the successful clearance of influenza virus infection and tissue remodelling without severe immunopathology. We also suggest that activation of these responses reflects the importance of protecting cells from damage by oxygen radicals and the need for the proliferation and repopulation of damaged or killed cells in the respiratory epithelium.
Increased activation of inflammatory and death receptor gene expression by r1918 influenza virus
To explore the host responses associated with the severe immunopathology of infections with r1918 virus, we characterized the functional relationship between genes that were induced more than twofold (P
< 0.01) only in r1918 virus infections. As shown in , many of these genes that were uniquely induced more than twofold in r1918-infected mouse lung compared with controls were involved in cell death and inflammatory responses and have been reported to have direct and/or indirect interactions. Infections with r1918 resulted in a statistically significant coordinated activation of multiple cell-death-related genes, particularly those related to death receptor responses. This is significant because activation of TNF-related pathways has been shown to be associated with H5N1 influenza viral infections of macrophages9
. Moreover, r1918 infection-induced expression of mRNAs for Fas, caspase-8 and caspase-9 is significant because of the central roles of these proteins in the induction of apoptosis through the death-receptor and mitochondrial apoptosis pathways. Caspase-9 is activated by the release of cytochrome c
from permeabilized mitochondria, and it has been reported that an influenza viral protein expressed from an alternative reading frame of the PB1 gene segment (PB1-F2) is trafficked to mitochondria and can induce apoptosis10
. It has recently been reported11
that PB1-F2 associates with components of the permeability transition pore complex and sensitizes cells to apoptotic stimuli, including TNF-α
, resulting in mitochondrial permeabilization and the release of cytochrome c
. We believe our data showing increased cell death responses in r1918-infected mice may be related to the effects of the 1918 PB1-F2 protein. However, these observations do not rule out additional functions of the 1918 polymerase proteins, because these genes are also important virulence factors for highly pathogenic avian influenza viruses12–17
Functional relationships of activated of cell death responses during infection with r1918 influenza virus
Taken together, our results indicate that enhanced inflammatory and cell death responses might be contributors to severe immunopathology associated with r1918 influenza virus infections, although it is not possible to completely exclude the possibility that increased immune and cell death responses may be a consequence rather than the cause of enhanced pathology of r1918 infection. Our study also shows the synergistic effect of the 1918 gene segments and that exceptional virulence and exaggerated host global gene expression response may be dependent on the complete constellation of 1918 genes. We acknowledge that a more complete understanding of r1918 virulence will require experiments in additional animal models, including ferrets, pigs and non-human primates. Characterization of the global host response to other highly virulent viruses will be important in understanding the role of the host in influenza virus pathology, and we are currently comparing these results with findings from similar studies of animals infected with highly pathogenic avian influenza viruses.