Volume 13, Number 9—September 2007
Research
Detecting HumantoHuman Transmission of Avian Influenza A (H5N1)
On This Page
Figures
Tables
Technical Appendices
Downloads
Highlight and copy the desired format.
EID  Yang Y, Halloran M, Sugimoto JD, Longini IM. Detecting HumantoHuman Transmission of Avian Influenza A (H5N1). Emerg Infect Dis. 2007;13(9):1348. https://dx.doi.org/10.3201/eid1309.070111 

AMA  Yang Y, Halloran M, Sugimoto JD, et al. Detecting HumantoHuman Transmission of Avian Influenza A (H5N1). Emerging Infectious Diseases. 2007;13(9):1348. doi:10.3201/eid1309.070111. 
APA  Yang, Y., Halloran, M., Sugimoto, J. D., & Longini, I. M. (2007). Detecting HumantoHuman Transmission of Avian Influenza A (H5N1). Emerging Infectious Diseases, 13(9), 1348. https://dx.doi.org/10.3201/eid1309.070111. 
Abstract
Highly pathogenic avian influenza A (HPAI) subtype H5N1 has caused family case clusters, mostly in Southeast Asia, that could be due to humantohuman transmission. Should this virus, or another zoonotic influenza virus, gain the ability of sustained humantohuman transmission, an influenza pandemic could result. We used statistical methods to test whether observed clusters of HPAI (H5N1) illnesses in families in northern Sumatra, Indonesia, and eastern Turkey were due to humantohuman transmission. Given that humantohuman transmission occurs, we estimate the infection secondary attack rates (SARs) and the local basic reproductive number, R_{0}. We find statistical evidence of humantohuman transmission (p = 0.009) in Sumatra but not in Turkey (p = 0.114). For Sumatra, the estimated household SAR was 29% (95% confidence interval [CI] 15%–51%). The estimated lower limit on the local R_{0} was 1.14 (95% CI 0.61–2.14). Effective HPAI (H5N1) surveillance, containment response, and field evaluation are essential to monitor and contain potential pandemic strains.
Highly pathogenic avian influenza A (HPAI) subtype H5N1 is repeatedly crossing the species barrier to humans. Since December 2003, a total of 291 cases of HPAI (H5N1) have been reported in humans, resulting in 172 deaths (i.e., 59% casefatality ratio) in 12 countries, mostly in Southeast Asia (1). Among these cases, 31 family clusters have been documented, ranging in size from 2 to 8 family members. How many of these clusters are due to a common avian source and how many are due to humantohuman transmission are important facts to determine. Should one of these HPAI (H5N1) strains gain the capacity for sustained humantohuman transmission, the resulting outbreak, if not contained, would spread worldwide through the global transportation network more rapidly than adequate supplies of vaccine matched to the new variant could be manufactured and distributed (2,3). We analyzed data from 2 of the largest of the familial clusters to ascertain if humantohuman transmission took place, and if so, how transmissible the strain was.
May 2006 Human Avian Influenza Family Cluster, Indonesia
During late April and early May 2006, a cluster of 8 cases of HPAI (H5N1) was detected and investigated by the Indonesian public health surveillance system in northern Sumatra (4–6). All casepatients were members of the same extended family. Seven of them resided within 3 adjacent houses in the village of Kubu Sembilang. The remaining patient resided with his immediate family in the village of Kabanjahe (≈10 km away).
The index patient was a 37yearold woman, thought to have been exposed to dead poultry and chicken fecal material before onset of illness. She also reportedly maintained a market stall that sold live chickens. Although her illness was not confirmed to have been caused by the (H5N1) avian influenza virus, her death on May 5, 2006, is suspected to be the result of HPAI (H5N1) infection because of her reported symptoms, illness progression, and prior contact with diseased or dead poultry.
Appendix Figure 1
Twenty members of her extended family are suspected to have been in contact with her, many during a family gathering on April 29, 2006 (7). At that time, she was manifesting symptoms (i.e., she had a heavy cough, was severely ill, and was prostrate). That night, 9 of these members slept in the same small room as she did (indicated by a black triangle in Appendix Figure 1). Of these 9 family members, 2 of her sons (15 and 17 years of age) and her 25yearold brother, who lived in Kabanjhe, became ill in the next 3 weeks. The sons died. The brother was the only person from this family cluster to recover.
Of the remaining 11 family members, 4 became ill and died. The 29yearold sister of the index patient, who lived in an adjacent house, became ill after she provided direct personal care to her ill sister (7). The 18monthold daughter of this sister also became ill after she was in the presence of the index patient with her mother. The 10yearold nephew of the index patient, who lived in the other house adjacent to hers, became ill after he attended the family gathering and frequently visited his aunt’s house. The nephew’s father became ill after he personally cared for his son. The possibility that HPAI (H5N1) was transmitted from the nephew to his father is also supported by genetic sequencing data (4). Though symptoms did not develop in the mother of the nephew, she was directly exposed to her husband during his illness. All casepatients, except for the index patient, were confirmed as influenza (H5N1) positive by PCR. The nephew’s mother was confirmed as influenza (H5N1) negative. As an intervention, 54 surviving relatives and close contacts were identified and placed under voluntary quarantine (7). All of these persons, except for pregnant women and infants, received oseltamivir prophylactically.
December 2005 Human Avian Influenza Family Cluster, Eastern Turkey
Appendix Figure 2
From December 18, 2005, (8) to January 15, 2006 (9), a cluster of 8 confirmed (H5N1) influenza cases was detected in Dogubayazit District in eastern Turkey (Appendix Figure 2) (10–13). These casepatients were among 21 members of 3 households located within 1.5 km of each other (14). All confirmed casepatients were hospitalized after onset of symptoms (9). Four of the confirmed casepatients died; the other 4 recovered (9). Ten of the remaining 14 household residents were hospitalized with avian influenzalike symptoms but were never confirmed to be (H5N1) infected (9). All but one of the hospitalized residents were children (6–15 years of age) (9).
Before onset of symptoms, 4 children from 1 household, 3 of whom had confirmed cases (including the index patient), were reported to have had close contact with the dead bodies of sick chickens (15). The 2 confirmed casepatients in the second household reportedly slaughtered a duck together on January 1, 2006, at the beginning of a dieoff in the household’s flock (14). Two of the remaining confirmed casepatients lived in the third household and had no history of contact with sick or dying poultry. The remaining confirmed case occurred in a fourth residence located near the first household (10), but because we lacked information on the number of household members and the casepatient’s exposure history, we excluded it from these analyses. Most, if not all, of the 21 residents attended a dinner hosted by the family of the index patient on December 24, 2006, while he was symptomatic (8).
Statistical Methods
We used a previously developed statistical transmission model (16,17) to test whether humantohuman transmission occurred, and if it did, to estimate transmission parameters. In the model, persons mix with one another in households and between households. In addition, we include a common source of infection due to zoonotic exposure. Mathematical and statistical details are given in the Technical Appendix [PDF  39 KB  2 pages].
Model of Probability of Transmission
Figure
We define p_{1} as the probability that an infectious household member infects another household member in 1 day. If the distribution of the infectious period is known, we can obtain the household secondary attack rate (SAR_{1}) from p_{1}, defined as the probability that an infectious household member infects another household member over his or her infectious period. Similarly, we define the daily transmission probability (p_{2}) and the community SAR (SAR_{2}) for between household spread. Finally, we define the daily probability (b) that any person is infected from a zoonotic source. The contact structure used for parameter estimation is shown in the Figure. We assume that the distributions of the incubation and infectious periods are predetermined by the investigator.
We establish the likelihood function for each person and then for the whole population for statistical inference. The likelihood function for a person is equivalent to the probability of observing the realized data on that person throughout the outbreak. The likelihood function for a person labeled i is built with the following steps: 1) Obtain the probability that person i is infected by an infectious source labeled j on day t, given person i is not infected up to day t – 1. If source j is a person, this probability is p_{1}, for the same household, or p_{2} for exposure in the community, multiplied by the probability of person j being infectious on day t. The probability of person j being infectious on day t is derived from the symptomonset day of person j and the distribution of the infectious period. If source j is zoonotic, the infection probability is b. The probability of escaping infection is simply 1 minus the corresponding probability of infection. 2) Take the product of the probabilities obtained in step 1 over all humans and zoonotic sources j to obtain the probability of person i escaping infection by any infectious source on day t. 3) Take the product of the probabilities obtained in step 2 over all days before and including day t to obtain the probability of person i escaping infection up to day t. 4) If person i is not infected by the end of the outbreak, the likelihood function for person i is the product of the probabilities of person i escaping infection up to the last day of observation. 5) If person i is observed to have symptom onset on day
The likelihood function for the whole population is the product of all the individual likelihood functions. In the event that humantohuman transmission occurs, SAR estimates are used to estimate the local basic reproductive number (R_{0}), which is defined as the average number of secondary cases infected by a typical index casepatient in the beginning of the outbreak (Technical Appendix [PDF  39 KB  2 pages]). There is potential for sustained transmission if R_{0} is >1. If humantohuman transmission is determined to be occurring, then the above parameters are estimated from the symptom dates and contact information from the population under study. Data on exposed persons who do not become ill form an important component of the inference procedure.
Statistical Test
We set up a statistical test with the null hypothesis being that no humantohuman transmission occurs, that is, p_{1} = p_{2} = 0. The alternative hypothesis is either p_{1} or p_{2} is not equal to 0, or both are not equal to zero. The test statistic we use is proportional to the ratio of the maximum value of the likelihood function assuming the null hypothesis is true (null likelihood) and the maximum value of the likelihood function at the estimated parameter values (full likelihood).
Specifically, we define the likelihood ratio test statistic as –2 log (the null likelihood function divided by the full likelihood function). If no humantohuman transmission occurs, the 2 likelihood functions would be roughly equal, and we expect to see a likelihood ratio close to 1, and, thus, a likelihood ratio statistic close to 0. A large value of the likelihood ratio statistic is evidence of deviation from the null hypothesis. The question is how to obtain a reference set of the likelihood ratio statistic values that we would see under the null hypothesis. Given no humantohuman transmission, all the observed casepatients must have been infected by the zoonotic source. Since the exposure to the zoonotic source is assumed constant for each person on each day, the null likelihood function will not change if we reassign the infection and symptom status of the observed casepatients to a different group of people in the population. By performing such reassignment many times, we obtained a collection of datasets that were each equally likely to have been observed had there been no humantohuman transmission. The values of the likelihood ratio statistic calculated from these datasets form the null distribution for statistical testing. This method is referred to as a permutation test. The p value is given by the proportion of the reference values that are equal to or larger than the observed likelihood ratio statistic value. More technical details are given in the online appendix.
The probability of infection by the zoonotic source may not be estimable together with SAR_{1} or SAR_{2} from an observed cluster. In such a situation, a statistical test of the occurrence of humantohuman transmission is still meaningful because the likelihood ratio test statistic is still estimable from the permuted datasets.
Data Required
A list of the inputs that are required for estimation and statistical testing are listed in the Table. Three categories of input parameters are required for this estimation model: outbreakwide, individual level, and analysis parameters. The duration of the outbreak, the duration of the incubation period for the pathogen, and the minimum and maximum durations of the infectious period for the pathogen are the required outbreakwide inputs. For each person, their residential location (neighborhood and household), their demographic characteristics (sex and age), and whether they were a casepatient or not are required input parameters. Casepatients require additional input of their illnessonset dates, types of outcome, outcome dates, and whether or not they are the index patient in the outbreak. Hospitalization and treatment dates (considered prophylactic for nonpatients) are optional input parameters for each person. For each person who visits another residence during the outbreak period, his or her identifiers, the neighborhood and household visited, and the start and end dates of the visit are required inputs. Analysisrelated inputs include the last date of community exposure to potential common sources of infection, the last date of observation, and inputs for R_{0} estimation (mean number of residents per household and mean number of outofresidence contacts per person per day). An expanded version of the model will require the input of other exposure information such as from schools or hospitals.
For the outbreak in Indonesia, Figure shows that the incubation period had a probable range of 3–7 days and the infectious period, a probable range of 5–13 days. Thus, we let the incubation period have a uniform distribution of 3–7 days (mean 5 days) and the infectious period a uniform distribution of 5–13 days (mean 9 days). For the data shown in Figure, only the household SAR (SAR_{1}) can be estimated. We determined that humantohuman spread did occur by rejecting the null hypothesis of no humantohuman transmission (p = 0.009). The estimated household SAR is 0.29 (95% confidence interval [CI] 0.15–0.51). Thus, a single infected person in a household infected another household member with the probability of 0.29. The average household size for rural Indonesia is ≈5 people. Because we do not have an estimate of the community SAR, we have an estimate of the lower limit of the local R_{0},_{,} i.e., 1.14 with a 95% CI of 0.61–2.14. A sensitivity analysis on the distribution of the incubation and infectious period shows that the test and estimates for SAR_{1} and R_{0} are insensitive to uncertainty about these distributions within plausible ranges.
For the outbreak in Turkey, all the parameters are estimable, but we do not reject the null hypothesis of no humantohuman transmission (p = 0.114). Our estimate of the daily probability of infection from the common source is 0.011 (95% CI 0.005–0.025).
We have presented statistical evidence that the strain of HPAI (H5N1) that caused the family cluster of human cases in northern Sumatra was spread from human to human and that the household SAR was 29%. This household SAR is similar to statistical estimates for interpandemic influenza A in the United States (12.7%–30.6%) (18,19). The mean incubation period of this strain appears to have been ≈5 days, nearly twice as long as for past pandemic strains and current interpandemic strains of influenza. The CI for the estimated lower bound for the local R_{0} covers 1. Therefore, even though we determined that humantohuman transmission probably occurred, whether the virus was capable of sustained humantohuman transmission is not clear. This virus may have required very close human contact to be transmitted. Even with no intervention, the finding that R_{0} = 1.14 indicates that the chance that a single introduction would result in any further spread is ≈12%. In addition, the reported prophylactic use of oseltamivir may have played some role in limiting further spread. We did not find statistical evidence of humantohuman spread for the outbreak in eastern Turkey. This does not mean that no lowlevel humantohuman spread occurred in this outbreak, only that we lack statistical evidence of such spread. The power would be too low to detect such spread for an outbreak with 7 total cases and small SARs (17).
We did not consider the role of heterogeneity—such as age, sex, treatment status, or quarantine—in transmission. The parameters could be made to be functions of timedependent covariates, as we have done with similar models (16,19,20). We can easily extend the model used here for covariates; however, we must have sufficient data to support such models.
Computer simulations have shown that the targeted use of influenza antiviral agents could be effective in containing a potential pandemic strain of influenza at the source (21,22), if initiated within 3 weeks of the initial case in the community, and if the R_{0} is <1.8. This strategy, known as targeted antiviral prophylaxis, involves treating identified index patients in a mixing group and offering a single course of prophylaxis to the contacts of these index patients in predefined close contact groups, i.e., households at a minimum but also possibly neighborhood clusters, preschool groups, schools, and workplaces. In addition, the voluntary household quarantine of suspected close contacts of casepatients was recommended. Targeted antiviral prophylaxis at the household and neighborhood cluster level was carried out for the outbreak in Sumatra.
Ascertaining whether a potential pandemic strain of influenza is capable of sustained humantohuman transmission and estimating key transmission parameters are important. To estimate more than the household SAR, more detailed community data need to be collected. This would include a complete census of potentially exposed households and persons in the area where immediate transmission could occur from both potential zoonotic and human sources. Such data would enable estimation of important parameters and a more complete estimate of the R_{0} rather than just the lower limit.
We have developed a software application, TRANSTAT, for implementing these analyses. This application provides a standalone environment for the entry, storage, and analysis of data from outbreaks of acute infectious diseases. A partial list of the input information is given in the Table. The statistical methods presented here can be applied to the data along with several standard epidemiologic tools. This information system would allow for realtime analysis and evaluation of control measures for an outbreak. We would encourage outbreak investigators to use this tool, taking care to input data on the exposed nonpatients as well as casepatients. The authors will provide a link to this software upon request.
Acknowledgments
This work was supported by the National Institute of General Medical Sciences MIDAS grant U01GM070749 and National Institute of Allergy and Infectious Diseases grant R01AI32042.
Dr Yang is a staff scientist in the Biostatistics and Biomathematics Program in the Division of Public Health Sciences at the Fred Hutchinson Cancer Research Center, Seattle, Washington. His primary research interest is in the statistical and mathematical analysis of infectious disease data and intervention studies.
References
 World Health Organization. Confirmed human cases of avian influenza A (H5N1). Epidemic and pandemic alert and response. 2007 Apr 11. [cited 2007 Apr 25]. Available from http://www.who.int/csr/disease/avian_influenza/country/cases_table_2007_04_11/en/index.html
 Webby RJ, Webster RG. Are we ready for pandemic influenza? Science. 2003;302:1519–22. DOIPubMed
 Stohr K. Avian influenza and pandemics—research needs and opportunities. N Engl J Med. 2005;352:405–7. DOIPubMed
 Butler D. Family tragedy spotlights flu mutations. Nature. 2006;442:114–5. DOIPubMed
 World Health Organization. Avian influenza—situation in Indonesia—update 12. Jakarta, Indonesia. Report no. 12. 2006 May 19. [cited 2006 Sep 22]. Available from http://www.who.int/csr/don/2006_05_16b/en/index.html
 Soebandrio A. Indonesia: avian flu from Indonesia’s experience and perspectives. Pandemic Preparedness and Infection Control. 2006 July 13–14.
 World Health Organization. Avian influenza—situation in Indonesia—update 16. Jakarta, Indonesia. Report no. 16. 2006 Jun 2. [cited 2006 Sep 16]. Available from http://www.who.int/csr/don/2006_05_31/en/index.html
 Recombinomics. H5N1 bird flu pandemic phase evolution. 2006 Jun 5. [cited 2007 Mar 16]. Available from http://www.recombinomics.com/news/06050601/h5n1_phase_evolution.html
 Recombinomics. Timeline for Kocyigit Ozcan family clusters in Dogubeyazit. January 22, 2006. [cited 2006 Mar 16]. Available from http://www.recombinomics.com/news/01220601/h5n1_kocyigit_ozcan_timeline.html
 World Health Organization. Avian influenza, situation in Turkey, update 2. Ankara, Turkey. Report no. 2. Epidemic and pandemic alert and response. 2006 Jan 9. [cited 2007 Feb 27]. Available from http://www.who.int/csr/don/2006_01_09/en/index.html
 World Health Organization. Avian influenza, situation in Turkey, update 4: sequencing of human virus. Ankara, Turkey. Report no. 4. Epidemic and pandemic alert and response. 2006 Jan 12. [cited 2007 Feb 27]. Available from http://www.who.int/csr/don/2006_01_12/en/index.html
 World Health Organization. Avian influenza, situation in Turkey, update 5. Ankara, Turkey. Report no. 5. Epidemic and pandemic alert and response. 2006 Jan 16. [cited 2007 Feb 27]. Available from http://www.who.int/csr/don/2006_01_16/en/index.html
 Lauer C. Turkey’s bird flu outbreak: one year later. Turkish Daily News. 2007 Jan 12. [cited 2007 Mar 20]. Available from http://www.turkishdailynews.com.tr/article.php?enewsid=63680
 World Health Organization. Shindo N. Avian influenza outbreak response in Turkey, 2006. Epidemic and pandemic alert and response. [cited 2007 Mar 20]. Available from http://www.col.opsoms.org/servicios/influenza/reunion/docs/trad/10_turkey_summary.pdf
 World Health Organization. Avian influenza, situation in Turkey. Ankara, Turkey. Epidemic and pandemic alert and response. Report no. 1. 2006 Jan 5. [cited 2007 Feb 27]. Available from http://www.who.int/csr/don/2006_01_05/en/index.html
 Yang Y, Longini IM, Halloran ME. Design and evaluation of prophylactic interventions using infectious disease incidence data from close contact groups. Appl Stat. 2006;55:317–30. DOI
 Yang Y, Longini IM, Halloran ME. A resamplingbased test to detect persontoperson transmission of infectious disease. Annals of Applied Statistics. 2007;1:211–28. [cited 2007 Jul 19]. Available from http://projecteuclid.org/dpubs/repository/1.0/disseminate?handle=euclid.aoas/1183143736&view=body&contenttype=pdfview_1
 Longini IM, Koopman JS, Monto AS, Fox JP. Estimating household and community transmission parameters for influenza. Am J Epidemiol. 1982;115:736–51.PubMed
 Longini IM, Koopman JS, Haber MJ, Cotsonis GA. Statistical inference for infectious diseases: riskspecific household and community transmission parameters. Am J Epidemiol. 1988;128:845–59.PubMed
 Rampey AH, Longini IM, Haber MJ, Monto AS. A discretetime model for the statistical analysis of infectious disease incidence data. Biometrics. 1992;48:117–28. DOIPubMed
 Longini IM, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings DA, Containing pandemic influenza at the source. Science. 2005;309:1083–7. DOIPubMed
 Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005;437:209–14. DOIPubMed
Figures
 Figure. Schematic of estimation method. An infectious person (in red) infects a susceptible person (in green) in the same household with probability of household secondary attack rate (SAR_{1}) and infects a...
 Appendix Figure 1. Exposure and disease events for each member of the family cluster in northern Sumatra, Indonesia. Dark boxes, duration of illness; white boxes without text, recovery period; thick dark...
 Appendix Figure 2. Exposure and disease events for each member of the family cluster in Eastern Turkey. Dark boxes, period of illness; white boxes without text, recovery period; thick dark vertical...
Table
Technical Appendix
Cite This ArticleDOI: 10.3201/eid1309.070111
Table of Contents – Volume 13, Number 9—September 2007
Comments
Please use the form below to submit correspondence to the authors or contact them at the following address:
Ira M. Longini, Jr, Program in Biostatistics and Biomathematics, Fred Hutchinson Cancer Research Center, Seattle, WA 981091024, USA;
Comment submitted successfully, thank you for your feedback.
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
 Page created: July 01, 2010
 Page last updated: July 01, 2010
 Page last reviewed: July 01, 2010
 Centers for Disease Control and Prevention,

National Center for Emerging and Zoonotic Infectious Diseases (NCEZID)
Office of the Director (OD)