Volume 22, Number 12—December 2016
Highly Pathogenic Reassortant Avian Influenza A(H5N1) Virus Clade 188.8.131.52a in Poultry, Bhutan
In India and Bangladesh, highly pathogenic avian influenza (HPAI) A(H5N1) viruses of the 184.108.40.206a genetic lineage have circulated in poultry since 2011 (1–3). Subtype H5N1 endemicity is complicated by co-circulating subtype H9N2, G1_Mideast lineage (4,5), which derives 5 internal genes from HPAI A(H7N3) virus from Pakistan (4). A reassortant H5N1 220.127.116.11a virus, rH5N1, with an H9N2-like polymerase basic protein 1 (PB1) gene (H7N3 origin), was reported in Bangladesh (2,5,6), India, and Nepal (7). However, its virulence and transmissibility are undetermined.
In Bhutan, the poultry sector consists of free-range backyard chickens, a rising number of commercial chicken farms, and domestic waterfowl in the south (8,9). Live-bird markets do not exist, but live birds are imported from India (8,9).
Bhutan’s poultry sector was severely affected by outbreaks of HPAI A(H5N1) clade 18.104.22.168 virus infection during 2012–2013 (10). Veterinary authorities enforced strict control measures, including depopulation of poultry in affected regions and burning of related housing and equipment (11). Illegal movement of poultry was the major source of outbreaks (11). Although the introduction of HPAI A(H5N1) from neighboring H5N1-endemic countries is a constant threat, the subtype is not yet entrenched in poultry in Bhutan.
We isolated HPAI A(H5N1) viruses from samples from 36 chickens and 9 wild birds in Bhutan, all from affected backyard farms adjacent to the highway connecting India with the capital, Thimphu (Figure 1; Technical Appendix 1[PDF - 744 KB - 9 pages] Table 1). Phylogenetic analysis (Technical Appendix 1[PDF - 744 KB - 9 pages]) suggested that the 2012–2013 outbreaks in Bhutan were caused by the rH5N1 genotype (22.214.171.124a/H9-like PB1 [H7N3 origin]), reported in Bangladesh and India at that time (Technical Appendix 1[PDF - 744 KB - 9 pages]Figures 1, 2; other data not shown). PB1 phylogeny suggested that this genotype underwent 4 independent reassortment events on the Indian subcontinent (Technical Appendix 1[PDF - 744 KB - 9 pages] Figure 2).
Antigenic analysis of selected H5N1 isolates from Bhutan (Technical Appendix 1[PDF - 744 KB - 9 pages]) showed homogeneity and a reactivity pattern similar to that of H5N1 reference viruses from Bangladesh (Table). Amino acid differences were observed between strains A/chicken/Bhutan/346/2012 (Ck/Bh/346) (rH5N1) and A/chicken/Bangladesh/22478/2014 (Ck/BD/22478), representing the parental H5N1 clade 126.96.36.199a genotype (pH5N1) (Technical Appendix 1[PDF - 744 KB - 9 pages] Table 2).
To assess whether the rH5N1-PB1 gene conferred a fitness advantage over the pH5N1 genotype, we examined replication kinetics in vitro (Technical Appendix 1[PDF - 744 KB - 9 pages]). The replication patterns of rH5N1 and pH5N1 were similar in Madin-Darby canine kidney (mammalian) cells (Figure 2, panel A). However, in chicken embryo fibroblasts (CEFs), Ck/Bh/346 (rH5N1) titers were significantly higher than those of Dk/BD/21326 (rH5N1) (p< 0.05) and Ck/BD/22478 (pH5N1) (p< 0.01) at 12 hours postinoculation (hpi) and those of Ck/BD/22478 (pH5N1) (p< 0.001), and Dk/BD/19097 (pH5N1) (p< 0.01) at 24 hpi. Dk/BD/21326 (rH5N1) had significantly higher titers than did Ck/BD/22478 (pH5N1) (p< 0.01) at 24 hpi (Figure 2, panel B). These results suggest rH5N1 viruses have a selective growth advantage in avian cells at early time points.
Next, we examined whether this growth advantage reflected higher pathogenicity or transmissibility for Ck/Bh/346 (rH5N1) in chickens than did Ck/BD/22478 (pH5N1) (online Technical Appendix 1). The 50% lethal dose (LD50) for chicken was 16 EID50 (50% egg infective dose) for Ck/Bh/346 (rH5N1) and 50 EID50 for Ck/BD/22478 (pH5N1). After inoculation with 30 LD50 and cohousing with naive contacts, all donors shed virus oropharyngeally and cloacally (Figure 2, panels C, D). All Ck/Bh/346 (rH5N1) donors died within 48 hpi, whereas only 50% of chickens inoculated with Ck/BD/22478 (pH5N1) died. Naive chickens in contact with donors inoculated with Ck/Bh/346 (rH5N1) or Ck/BD/22478 (pH5N1) became infected by day 2 after contact (Figure 2, panel C), started shedding cloacally on day 3 (Figure 2, panel D), and died by day 4. On day 3 after contact, Ck/Bh/346 (rH5N1) contacts had oropharyngeal and cloacal titers > 1 log10 EID50/mL higher than those of Ck/BD/22478 (pH5N1) contacts (Figure 2, panels C, D), but the difference was not significant.
We placed Ck/Bh/346 (rH5N1) and Ck/BD/22478 (pH5N1) in direct competition by co-housing chickens inoculated with each virus with naive contacts (Technical Appendix 1[PDF - 744 KB - 9 pages]). All donors shed virus oropharyngeally and cloacally, starting at 1 day postinoculation (dpi). By day 3 after contact, real-time reverse transcription PCR to detect PB1 (Technical Appendix 1[PDF - 744 KB - 9 pages]) revealed that 7 of 8 naive contacts simultaneously exposed to both viruses were infected with Ck/Bh/346 (rH5N1) alone, none was infected with Ck/BD/22478 (pH5N1) alone, and 1 was co-infected with both viruses. Thus, despite the lower infectious dose used for 30 LD50, Ck/Bh/346 (rH5N1) killed inoculated chickens faster than did Ck/BD/22478 (pH5N1) and was transmitted faster and more efficiently to naive contacts.
We assessed the risk for human infection with rH5N1 by investigating its pathogenicity and transmissibility in ferrets (Technical Appendix 1[PDF - 744 KB - 9 pages]). Donors shed 4.5 log10 EID50/mL and 3.4 log10 EID50/mL in nasal wash samples at 2 dpi and 4 dpi, respectively, but cleared the virus by 6 dpi. No direct or aerosol contacts shed virus, suggesting that Ck/Bh/346 (rH5N1) was not transmitted (data not shown). No Ck/Bh/346 (rH5N1)–inoculated ferrets lost >5% of their body weight or showed elevated body temperature (data not shown). Histopathologic analysis showed that 1 donor, who was lethargic at 3-10 dpi, had mild meningoencephalitis at 14 dpi (Technical Appendix 2[PDF - 744 KB - 9 pages]). A nucleocapsid protein–positive cell was detected in the brain, suggesting that Ck/Bh/346 (rH5N1) is neurotropic. The other ferrets showed no clinical signs of disease. Virus replication was detected in the lung at 4 dpi (log10 2.75 EID50/g) (Technical Appendix 2[PDF - 744 KB - 9 pages]).
Our study revealed that the viruses that caused the 2012 outbreaks in Bhutan belonged to the rH5N1 genotype (188.8.131.52a/H9-like PB1 [7:1]), whereas neither H9N2 nor the pH5N1 genotype have been detected there. rH5N1 has been isolated sporadically at live-bird markets and from chickens on farms where outbreaks occurred in Bangladesh (5,6), India (12), and Nepal (7) in 2011–2013. The lack of data on the effect of the H9-like PB1 gene on the virulence of rH5N1 makes determining its pathogenicity and transmissibility a critical public-health goal for Bhutan and neighboring countries.
Ck/Bh/346 (rH5N1) killed inoculated chickens faster than did Ck/BD/22478 (pH5N1), despite the lower infectious dose used for Ck/Bh/346. In CEFs, Ck/Bh/346replicated with greater efficiency during the first 36 hpi than did Ck/BD/22478, which possibly explains why rH5N1 transmits more efficiently to naive chickens when directly competing with pH5N1. How faster replication contributes to the increased mortality rate of naive chickens might be crucial to eradicating the disease in Bhutan. In a mountainous region with widely separated villages, small-scale poultry farming, and no live-bird markets, the severity and rapid onset of the infection could lead to host-resource exhaustion and self-limitation.
To determine whether the reassortant PB1 gene accounts for the observed phenotypic properties of rH5N1, reverse genetics experiments are required. Despite its enhanced transmissibility, rH5N1 did not supplant pH5N1 in India or Bangladesh after undergoing multiple reassortment events. Possible reasons for this include the involvement of other influenza subtypes, which would complicate the competition/transmission model, especially at live-bird markets, as well as the large duck population, which is prone to inapparent HPAI infection (indicating possible underreporting).
Our ferret model results suggest that avian-to-human transmission of rH5N1 is possible. However, further adaptation to the host is necessary for rH5N1 to become transmissible among mammals. Similar results have been reported for H5N1 clade 184.108.40.206 (13), H5N1 clade 2.3.4 (14), and H5Nx clade 220.127.116.11 (15). rH5N1 is potentially neurotropic, manifesting clinically as mild meningoencephalitis with no obvious respiratory involvement. This finding has implications on early diagnosis and the use of antiviral drugs during the first 48 h of clinical onset for optimal therapeutic effect.
pH5N1 and H9N2 virus strains will likely continue to co-circulate on the Indian subcontinent, enabling further reassortment events. Our results highlight the need for active surveillance and full-genome sequencing of all H5N1 virus isolates.
Dr. Marinova-Petkova was a postdoctoral research associate at St. Jude Children’s Research Hospital, Memphis, Tennessee, USA, while this research was conducted. She is now affiliated with the Centers for Disease Control and Prevention, Atlanta, Georgia, USA, where her research interests include emerging influenza viruses at the animal–human interface, the evolution of influenza A viruses, and animal models for studying influenza pathogenesis and transmission.
We thank M. Kamrul Hasan and Sharmin Akhtar; Jasmine Turner for technical support; Lisa Kercher, Beth Little, and David Carey for help in the enhanced animal Biosafety Level 3 laboratory during animal experiments; James Knowles for administrative support; Elizabeth Stevens for generating Figure 1; Jeremy Jones for critically reviewing the manuscript; and Vani Shanker and Keith A. Laycock for scientific editing.
This work was funded by contracts HHSN266200700005C and HHSN272201400006C from the National Institute of Allergy and Infectious Diseases, the National Institutes of Health, and the Department of Health and Human Services, and by the American Lebanese Syrian Associated Charities.
- Nagarajan S, Tosh C, Smith DK, Peiris JSM, Murugkar HV, Sridevi R, et al. Avian influenza (H5N1) virus of clade 2.3.2 in domestic poultry in India.PLoS One. 2012;7:e31844.
- Gerloff NA, Khan SU, Balish A, Shanta IS, Simpson N, Berman L, et al. Multiple reassortment events among highly pathogenic avian influenza A(H5N1) viruses detected in Bangladesh. Virology. 2014;450-451:297–307.
- Marinova-Petkova A, Feeroz MM, Rabiul Alam SM, Kamrul Hasan M, Akhtar S, Jones-Engel L, et al. Multiple introductions of highly pathogenic avian influenza H5N1 into Bangladesh. Emerg Microbes Infect. 2014;3:e11.
- Shanmuganatham K, Feeroz MM, Jones-Engel L, Smith GJD, Fourment M, Walker D, et al. Antigenic and molecular characterization of avian influenza A(H9N2) viruses, Bangladesh. [Epub 2013 Sep.]. Emerg Infect Dis. 2013;19:•••.
- Marinova-Petkova A, Shanmuganatham K, Feeroz MM, Jones-Engel L, Hasan MK, Akhtar S, et al. The continuing evolution of H5N1 and H9N2 influenza viruses in Bangladesh between 2013 and 2014. Avian Dis. 2016;60(Suppl):108–17.
- Monne I, Yamage M, Dauphin G, Claes F, Ahmed G, Giasuddin M, et al. Reassortant avian influenza A(H5N1) viruses with H9N2-PB1 gene in poultry, Bangladesh. Emerg Infect Dis. 2013;19:1630–4.
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL. GenBank. Nucleic Acids Res. 2005;33:D34–8.
- Nidup K, Tshering P. Status of the family poultry production and HPAI in Bhutan. In: Proceedings of the 8th Asian Pacific Poultry Conference of the World’s Poultry Science Association; Bangkok, Thailand; March 2007. p. 78–83 [cited 2016 Mar 8]. http://cms.cnr.edu.bt/cms/files/docs/File/Karma%20Nidup/Poultry/FamilyPoultry&HPAIinBhutan.pdf
- Ministry of Agriculture and Forests of the Royal Government of Bhutan. Bhutan RNR statistics, 2015 [cited 2016 Mar 8]. http://www.moaf.gov.bt/bhutan-rnr-statistics-2015-online/
- World Organisation for Animal Health. Update on highly pathogenic avian influenza in animals (type H5 and H7). Immediate notification report, Bhutan. 2012 [cited 2016 Feb 8]. https://web.oie.int/wahis/reports/en_imm_0000011465_20120110_120756.pdf
- Ministry of Agriculture and Forests of Bhutan. RNR newsletter. 2013;I(3) [cited 2016 Feb 8]. http://www.moaf.gov.bt/rnr-newsletter-september-2013/
- Bhat S, Bhatia S, Pillai AS, Sood R, Singh VK, Shrivas OP, et al. Genetic and antigenic characterization of H5N1 viruses of clade 18.104.22.168 isolated in India. Microb Pathog. 2015;88:87–93.
- Xu L, Bao L, Yuan J, Li F, Lv Q, Deng W, et al. Antigenicity and transmissibility of a novel clade 22.214.171.124 avian influenza H5N1 virus. J Gen Virol. 2013;94:2616–26.
- Sun H, Pu J, Wei YD, Sun Y, Hu J, Liu L, et al. Highly pathogenic avian influenza H5N6 viruses exhibit enhanced affinity for human type sialic acid receptor and in-contact transmission in model ferrets. J Virol. 2016. Epub 2016 Apr 27.
- Pulit-Penaloza JA, Sun X, Creager HM, Zeng H, Belser JA, Maines TR, et al. Pathogenesis and transmission of novel highly pathogenic avian influenza H5N2 and H5N8 viruses in ferrets and mice. J Virol. 2015;89:10286–93.
Suggested citation for this article: Marinova-Petkova A, Franks J, Tenzin S, Dahal N, Dukpa K, Dorjee J, et al. Highly pathogenic reassortant avian influenza A(H5N1) virus clade 126.96.36.199a in poultry, Bhutan. Emerg Infect Dis. 2016 Dec [date cited]. http://dx.doi.org/10.3201/eid2212.160611
1These authors contributed equally to this article.
2Current affiliation: Centers for Disease Control and Prevention, Atlanta, Georgia, USA.