Fatal outcome of human influenza A (H5N1) is associated with high viral load and hypercytokinemia.

Avian influenza A (H5N1) viruses cause severe disease in humans, but the basis for their virulence remains unclear. In vitro and animal studies indicate that high and disseminated viral replication is important for disease pathogenesis. Laboratory experiments suggest that virus-induced cytokine dysregulation may contribute to disease severity. To assess the relevance of these findings for human disease, we performed virological and immunological studies in 18 individuals with H5N1 and 8 individuals infected with human influenza virus subtypes. Influenza H5N1 infection in humans is characterized by high pharyngeal virus loads and frequent detection of viral RNA in rectum and blood. Viral RNA in blood was present only in fatal H5N1 cases and was associated with higher pharyngeal viral loads. We observed low peripheral blood T-lymphocyte counts and high chemokine and cytokine levels in H5N1-infected individuals, particularly in those who died, and these correlated with pharyngeal viral loads. Genetic characterization of H5N1 viruses revealed mutations in the viral polymerase complex associated with mammalian adaptation and virulence. Our observations indicate that high viral load, and the resulting intense inflammatory responses, are central to influenza H5N1 pathogenesis. The focus of clinical management should be on preventing this intense cytokine response, by early diagnosis and effective antiviral treatment.

Main

Influenza H5N1 viruses cause severe and often fatal disease in humans that is characterized by fulminant pneumonia and multi-organ failure[1,2]. High replication efficiency, broad tissue tropism and systemic replication seem to determine the pathogenicity of H5N1 viruses in animals[3,4,5]. To examine the relevance of these viral properties in the context of human disease, we carried out virological analyses in respiratory and non-respiratory specimens of 18 previously healthy individuals with influenza H5N1 who were admitted to referral hospitals in Ho Chi Minh City during the years 2004 and 2005, of whom 13 died. (Table 1). For comparison, we studied eight patients who were hospitalized during the same period with human influenza H3N2 or H1N1. These patients presented earlier in the course of illness (Table 1), which may be explained by their origin from Ho Chi Minh City or neighboring provinces, in contrast with H5N1 patients who were mostly from more distant provinces.

Table 1

Patient characteristics

				H5N1
	H5N1	H3/H1	P	Fatal	Not fatal	P
Number of patients	18	8a		13	5
Age (years)	18 (6–35)	14(1–49)	0.71	16 (6–35)	24 (8–26)	0.33
Illness day at presentation	6 (4–8)	2 (2–3)	<0.001	6 (4–7)	6 (4–8)	0.14
Patients treated (number):
Antibioticsb	18	8		13	5
Oseltamivirc	17	8		12	5
Time to death (days since admission)	4 (1–16)	NA		4 (1–16)	–
Time to discharge (days since admission)	27 (10–38)	9 (5–15)	0.037	–	27 (10–38)

Data are given as medians (range).

a Five patients were diagnosed with influenza H3N2 and three patients with influenza H1N1.

bAll patients received third-generation cephalosporins in combination with either aminoglycosides, vancomycin or azithromycin.

c Oseltamivir was given at an oral dose of 75 mg b.i.d. in adults with weight-based adjustments for children as recommended by the manufacturer.

Despite their presentation late in the course of illness, we were able to isolate virus from pharyngeal specimens of 12 of 16 H5N1-infected individuals (Table 2). Genetic characterization and phylogenetic analysis revealed that all viral strains were of the genotype Z, H5N1 sublineage of viruses prevalent in Vietnam, Cambodia and Thailand, as previously reported[10]. Pairwise comparison of all gene segments of viruses isolated from eight fatal and four surviving cases did not reveal unique amino acid changes in either group. No viruses contained Glu92 in the NS1 protein, which is associated with increased virulence of H5N1 viruses[6], but all contained the recently reported PDZ-domain ligand ESEV[11]. An E627K substitution in the viral polymerase basic protein 2 (PB2), which is associated with adaptation and virulence of H5N1 viruses in mammals[12], was present in five of eight isolates from fatal cases and in three of four isolates from patients who survived. There was no association between the presence of Lys627 and clinical outcome (data not shown). Notably, three of four viruses without this substitution, but none of the viruses containing Lys627, contained an D701N substitution in PB2 that has been associated with adaptation of H7N7 viruses to mammalian cells[13]. This suggests that the D701N substitution may compensate for the absence of Lys627 in conferring enhanced viral polymerase activity and virulence in the mammalian host. Other amino-acid residues in the viral polymerase complex implicated in mammalian adaptation[13], including Pro13 in polymerase basic protein 1 (PB1) and Arg615 in polymerase acidic protein (PA), were present in all viruses. These observations are in agreement with recent studies demonstrating the importance of the viral polymerase complex for H5N1 virulence in mammals[14].

Table 2

Detection of influenza virus in respiratory and non-respiratory sites

					H5N1
		H5N1	H3/H1	P	Fatal	Not fatal	P
Nasopharynx
Virus isolation rate (positive/tested; %)		12/16a (75)	NA		8/12 (67)	4/4 (100)
Detectable RNA (positive/tested; %)	Nose	13/17 (76)	6/8 (75)		10/12 (83)	3/5 (60)
	Throat	18/18 (100)	8/8 (100)		13/13 (100)	5/5 (100)
Viral load (median; range)	Nose	5.5 (und.–8.1)	4.5 (und.–7.7)	0.59	5.8 (und.–8.1)	4.5 (und.–6.4)	0.20
	Throat	7.0 (4.3–8.2)	4.8 (4.2–5.8)	0.003	7.5 (4.7–8.2)	5.9 (4.3–7.0)	0.058
	P	0.001	0.87
Rectum b
Virus isolation rate (positive/tested; %)		1/7 (14)	NA		1/7 (14)	NA
Detectable RNA (positive/tested; %)		5/7 (71)	NA		5/7 (71)	NA
Viral load (median; range)		4.8 (3.6–5.8)	NA		4.8 (3.6–5.8)	NA
Blood c
Virus isolation rate (positive/tested; %)		1/6 (17)	NA		1/6 (17)	NA
Detectable RNA (positive/tested; %)		9/16 (56)	0/6 (0)	0.046	9/11 (82)	0/5 (0)	0.005
Viral load (median; range)		4.5 (3.2–-5.7)	Und.		4.5 (3.2–5.7)	Und.

Rates of virus isolation, detection of viral RNA and viral loads in nasopharyngeal, rectal and blood specimens of patients with influenza H5N1 and H3N2 or H1N1. Viral loads are given as log10 cDNA copies per ml of viral transport medium. Und., below detection limit; NA, not assessed.

aVirus isolation was not performed in two cases due to insufficient specimens.

bRectal swabs were obtained after 6–11 d (median 8) of illness.

cBlood specimens were obtained after 5–9 d (median 7) of illness.

Viral RNA could be detected at high levels in nasopharyngeal specimens of H5N1-infected individuals, without a clear decline, when viral loads were plotted against time after onset of illness (Table 2 and Supplementary Fig. 1 online). Prolonged viral shedding is not unique to H5N1 influenza but is also observed in young children with human influenza, probably reflecting the absence of pre-existing immunity.

In H5N1-infected individuals, viral RNA levels were higher in pharyngeal specimens than in nasal specimens (Table 2). Tracheal aspirates obtained from two patients during the course of illness showed more prolonged detection of viral RNA and at substantially higher levels as compared with corresponding nasopharyngeal swabs (Supplementary Fig. 2 online). Particularly in the pharynx, viral RNA levels were higher in H5N1-infected individuals than in individuals infected with human influenza H3N2 or H1N1, despite presentation later in the course of the illness (Table 2 and Supplementary Fig. 1). The highest viral RNA levels occurred in H5N1-infected individuals who died (Table 2), suggesting that the level of viral replication influences outcome. Viral RNA could also be detected in blood specimens from 9 of 16 H5N1-infected individuals and in the rectal swabs from 5 of 7 of these individuals (Table 2). It remains unclear whether the detection of virus in rectal specimens reflects genuine gastrointestinal infection. However, symptoms of diarrhea frequently occur during influenza H5N1 (refs. 1,2,15) and were present in three of the five individuals with detectable viral RNA in the rectum. No viral RNA was present in the blood of surviving H5N1-infected individuals or in those with human influenza. Individuals with detectable H5N1 RNA in blood also had higher pharyngeal viral loads (median 7.4 log10 copies per ml; range 5.1–7.4) than those without evidence of H5N1 RNA (6.0 log10 copies per ml; range 4.3–7.0; P = 0.021), suggesting that the presence of viral RNA in blood reflects an overall high viral burden.

Virus was also isolated from one rectal and one serum specimen from a patient who presented with diarrhea and coma as the only initial symptoms[16]. The failure to isolate virus from remaining non-respiratory specimens (Table 2) could suggest that the viral RNA does not necessarily reflect replicating virus, but may also reflect relatively low viral loads, loss of replication-competent virus during storage and limitations in the sensitivity of viral culture. The recent isolation of H5N1 virus from the plasma of a patient from Thailand[17] confirmed that viremia can occur. Although viremia indicates the potential of the H5N1 virus to disseminate systemically, it remains unclear whether virus replication commonly occurs in non-respiratory organs. Postmortem examinations have not shown evidence of extrapulmonary infection, but revealed reactive hemophagocytosis (which is believed to be a cytokine-driven condition) as the most prominent feature[18,19]. In vitro and animal experiments also implicate cytokine dysregulation in H5N1 pathogenesis, as characterized by evasion of the antiviral effects of interferons and tumor necrosis factor (TNF)-α6, in addition to induction of proinflammatory cytokines and chemokines by H5N1 viruses[7,8,9].

Observations in our subjects support a role for the inflammatory response in the pathogenesis of human H5N1 disease. Peripheral lymphopenia often accompanies H5N1 infection[1,2]. Measurement of peripheral blood lymphocyte counts in a subset of our subjects showed low total and CD3+ lymphocyte counts and inverted ratios of CD4+ to CD8+ cells in H5N1-infected individuals, particularly in those who died (Supplementary Table 1 online). Total and CD3+ lymphocyte numbers in H5N1-infected individuals correlated inversely with pharyngeal viral RNA load (Fig. 1a,b), suggesting an association between lymphopenia and the level of viral replication. In addition to decreased production of lymphocytes (as suggested by bone marrow hypoplasia reported in one individual[18]), or increased destruction through apoptosis (as suggested by in vitro and mouse experiments with H5N1 viruses[20,21]), peripheral blood lymphopenia may result from trafficking of lymphocytes to infected tissue. However, postmortem studies in H5N1-infected individuals have not shown a predominance of lymphocytes, but instead of macrophages, in pulmonary infiltrates[19]. The latter observation is in agreement with the chemokine profile that was measured in peripheral blood of the individuals in this study (Table 3): levels of IP-10, MIG and MCP-1 (chemoattractants of monocytes and macrophages that are produced in bronchial epithelial cells and alveolar macrophages[9,22,23]) were elevated in patients with avian and human subtypes of influenza, but were higher in H5N1-infected individuals and particularly high in those who died (Table 3). Levels of the neutrophil chemoattractant interleukin (IL)-8 were also elevated in H5N1-infected individuals, particularly in those who died (Table 3). This chemokine is produced by bronchial epithelial cells and may function in the pathogenesis of acute respiratory distress syndrome (ARDS)[24], which may be particularly relevant to H5N1 influenza, as progression to respiratory failure is associated with development of ARDS[15].

Figure 1

Correlations between virus load and immunological parameters in influenza H5N1 infection.

(a–g) Correlations between pharyngeal viral RNA load and total peripheral blood lymphocyte numbers (a) and CD3-positive lymphocyte counts (b), and plasma levels of IP-10 (c), MCP-1 (d), IL-8 (e), IL-6 (f) and IL-10 (g). Throat and blood specimens were obtained in parallel on the same day. Spearman's correlation coefficients (ρ) and P values are given in each graph. No significant correlations were observed for MIG (ρ = 0.41, P = 0.10) and IFN-γ (ρ = 0.48, P = 0.062).

Table 3

Levels of chemokines and cytokines in the peripheral blood

						H5N1
	H5N1a (n = 16)	P b	H3/H1 (n = 6)	P c	Controls (n = 15)	Fatal (n = 11)	Not fatal (n = 5)	P
IP-10	5.1 (3.5–6.3)	0.005	3.8 (3.4–4.6)	0.001	2.7 (2.4–3.8)	5.4 (3.5–6.3)	4.2 (4.0–5.0)	0.031
MCP-1	2.4 (1.5–4.0)	0.083	1.9 (und.–2.4)	0.045	1.4 (und.–2.0)	2.8 (2.0–4.0)	1.8 (1.5–2.3)	0.015
MIG	4.3 (3.1–5.2)	0.013	3.2 (2.9–3.9)	0.002	2.6 (2.2–3.3)	4.6 (3.3–5.2)	3.3 (3.1–4.2)	0.011
IL-8	2.0 (0.7–3.2)	0.001	0.8 (0.4–1.5)	0.34	0.7 (und.–1.0)	2.4 (1.1–3.2)	1.7 (0.7–1.9)	0.020
IL-10	1.5 (und.–2.8)	0.002	–1.0 (und.–0.4)	0.85	–1.0 (und.–1.0)	1.6 (und.–2.8)	0.8 (und.–2.2)	0.6
IL-6	2.1 (und.–3.7)	0.001	–0.2 (und.–0.7)	0.30	–1.0 (und.–1.0)	2.2 (1.5–3.7)	1.0 (und.–2.4)	0.054
IFN-γ	2.0 (und.–4.2)	0.029	0.1 (und.–2.4)	0.42	–1.0 (und.–1.4)	2.3 (1.0–4.2)	2.0 (und.–2.6)	0.2

Levels of chemokines and cytokines in the peripheral blood of patients with influenza H5N1 and H3N2 or H1N1. Levels are given as median log10 pg per ml (range).

a Plasma levels of chemokines and cytokines in H5N1 patients were all higher than in healthy controls at <0.001 significance levels.

bComparison between H5N1 and H3/H1 patients.

cComparison between H3/H1 patients and healthy controls. und., undetectable.

Measurement of cytokines in plasma revealed elevated levels of IL-10, IL-6 and IFN-γ in H5N1-infected individuals (Table 3); the latter is a strong inducer of IP-10 and MIG expression by bronchial epithelial cells[22]. Plasma levels of other cytokines were not higher than healthy controls, although low levels of TNF-α, IL-1b and IL-12 were observed in four H5N1-infected individuals (data not shown). Although levels of cytokines in lung tissue may be more relevant to pathogenesis than levels in blood[25], it is likely that blood cytokine levels reflect the events in the highly vascular lung. In vitro infection of bronchial epithelial cells and macrophages with H5N1 viruses is associated with hyperinduction of chemokines and cytokines, suggesting that intrinsic viral properties may contribute to the increased levels and to pathogenesis. However, plasma levels of IP10, MCP-1, IL-8, IL-6 and IL-10 in our subjects correlated with pharyngeal H5N1 load (Fig. 1c–g), which indicates that the observed hypercytokinemia and hyperchemokinemia reflect, at least partly, increased viral replication.

Our observations point to a central role for high viral burden in the pathogenesis of human H5N1 disease and suggest that timely suppression of viral replication should remain the mainstay for treatment of influenza H5N1. This is supported by our previous data demonstrating that successful control of viral replication by antivirals was associated with a good clinical outcome[26]. The seemingly limited clinical efficacy of antiviral treatment of H5N1 influenza may be due to the inability of antivirals to interfere with the sequence of events leading to the observed hypercytokinemia, and associated complications, when treatment is started late in the course of illness. Although immunomodulatory treatment has potential benefits at this stage, the focus of clinical management should be on preventing the intense cytokine response by early diagnosis and effective antiviral treatment.

Methods

Clinical specimens.

At admission, we collected pharyngeal and nasal swabs in viral transport medium (vtm: minimum essential medium Eagle with Hanks' salts, supplemented with 0.5% gelatin and antibiotics (Sigma-Aldrich)). We obtained repeat throat swabs during admission from 8 H5N1-infected individuals[26], rectal swabs from 7 such individuals, and serum or plasma from 16 H5N1-infected individuals and 6 individuals with influenza H3N2 or H1N1. We used plasma specimens from 15 healthy donors for comparison of cytokine and chemokine levels. All specimens were stored at −80 °C before analysis.

Experimental ethics policy.

This study was approved by the institutional review board of the Hospital for Tropical Diseases, Ho Chi Minh City, and the Oxford Tropical Research Ethical Committee. Informed consent was obtained from all participating patients or their parents or legal guardians.

RT–PCR.

We purified nucleic acids from 100 μl of pharyngeal, nasal and rectal swab specimens in vtm, or 100 μl of serum or plasma, using a previously decscribed method[27]. After reverse transcription, influenza A viral RNA was detected and quantified by real-time PCR targeted at a highly conserved region of the matrix gene, as previously described[16]. Quantitative analyses of influenza A viral load in clinical specimens were performed batch-wise. Influenza A virus was subtyped using H5-, H3- and H1-specific primers as previously described[1]. The preparation of reagents, the extraction of nucleic acids, and nucleic acid amplification and analysis were performed in physically separated laboratories.

Viral culture.

We isolated influenza A H5N1 viruses by cell culture in Madin-Darby canine kidney cells in biosafety level III laboratory facilities and identified influenza viruses by serotype-specific RT-PCR and hemagglutination inhibition assays, as previously described[16].

Sequence analysis.

We carried out sequence analysis of 12 H5N1 virus isolates as previously described[28,29], using BigDye Terminator v3.1 Cycle Sequencing Kit on an ABI PRISM 3700 DNA Analyzer (Applied Biosystems) following the manufacturers' instructions. We sequenced all gene segments from all viruses completely, with the exception of PA of Vietnam/CL2/04 and PB2 of Vietnam/CL17/04 (owing to technical difficulties). From the latter we obtained a partial sequence encompassing residue 627. Sequence fragments were assembled and edited with Lasergene (version 6.0; DNASTAR). Alignments and residue analyses were performed in BioEdit version 7. MrModeltest 2.2 was used to determine the appropriate DNA substitution model for constructing neighbor-joining trees in PAUP* 4.0. The sequences of all eight gene segments were compared pairwise between viruses isolated from patients who died and those who survived.

Flow cytometry and cytokine measurements.

We carried out flow cytometric analyses of whole-blood stained with fluorochrome-conjugated monoclonal antibodies on a FACScalibur flow cytometer (Becton Dickinson). Plasma concentrations of the cytokines IFN-γ, IL-2, IL-5, IL-4, IL-6, IL-10, IL-1β, TNF-α and IL-12p70 and of the chemokines MCP-1, IP-10, MIG, RANTES and IL-8 were measured using cytometric bead-array assays (Becton Dickinson) according to the manufacturers' instructions, with the exception that all samples were fixed in 2% paraformaldehyde before analysis. IFN-α plasma concentrations were measured using a commercial capture ELISA assay (Biosource International) according to the manufacturers' instructions.

Statistical analysis.

Virus load and levels of cytokines and chemokines were analyzed after log transformation. For statistical purposes, the lower detection limit of the assay (1,000 cDNA copies per ml) was used in cases of negative test results of virus load measurements, and a value of 0.1 pg ml−1 was used in cases of undetectable concentrations of cytokines or chemokines in specimens. The Mann-Whitney U or Fisher exact tests were used for group comparisons of numerical and categorical data respectively, and the Wilcoxon ranked-sum test was used for comparison of viral load in paired throat and nasal swabs. Correlations between viral load and immunological parameters in specimens obtained on the same day were calculated using Spearman's rank correlation test. For all analyses, a P value of less than 0.05 derived from a two-tailed test was considered significant. All statistical analyses were performed with SPSS 14.0 for Windows software (SPSS Inc).

Accession codes.

GenBank: DQ497719-497729, DQ535724 (HA); DQ493068-493071, DQ493073-493078, DQ250160, DQ535726 (NA); DQ492980-492990, DQ535725 (M); DQ493156-493166, DQ535727 (NP); DQ493244-493254, DQ535728 (NS); DQ493332-493341, DQ535729 (PA); DQ493418-493428, DQ535730 (PB1); DQ492895-492902, DQ535731 (PB2). Virus names and the clinical outcome of patients from whom the viruses were isolated are shown in Supplementary Table 2 online.

Note: is available on the Nature Medicine website.

Supplementary Fig. 1

Pharyngeal virus load in relation to onset of illness during H5N1 influenza and human influenza. (PDF 425 kb)

Supplementary Fig. 2

Nasal, pharyngeal and tracheal virus load during the course of illness in two H5N1 patients (PDF 650 kb)

Supplementary Table 1

Peripheral blood lymphocyte counts and CD4+/CD8+ ratios. (PDF 59 kb)

Supplementary Table 2

GenBank accession numbers of H5N1 viral sequences and clinical outcome of patients. (PDF 64 kb)