We characterized influenza A(H1N1)pdm09 isolates from large-scale outbreaks that occurred in Nepal and India in early 2015. Although no specific viral features, which may have caused the outbreaks, were identified, an S84N substitution in hemagglutinin was frequently observed. Chronological phylogenetic analysis revealed that these Nepalese and Indian viruses possessing the S84N substitution constitute potential ancestors of the novel genetic subclade 6B.1 virus that spread globally in the following (2015/16) influenza season. Thus, active surveillance of circulating influenza viruses in the Southern Asia region, including Nepal and India, would be beneficial for detecting novel variant viruses prior to their worldwide spread.This work was supported, in part, by a Grant-in-Aid for Emerging and Re-emerging Infectious Diseases from the Ministry of Health, Labor and Welfare, Japan, as well as by JSPS KAKENHI grant JP26460816. DFB and DJS wish to acknowledge the support received from the Center for Research on Influenza Pathogenesis (CRIP) funded by the NIAID and BARDA contracts HHSN272201400008C and HHSO100201500033C, respectively

Burke, DF

Influenza and Other Respiratory Viruses

We characterized influenza A(H1N1)pdm09 isolates from large-scale outbreaks that occurred in Nepal and India in early 2015. Although no specific viral features, which may have caused the outbreaks, were identified, an S84N substitution in hemagglutinin was frequently observed. Chronological phylogenetic analysis revealed that these Nepalese and Indian viruses possessing the S84N substitution constitute potential ancestors of the novel genetic subclade 6B.1 virus that spread globally in the following (2015/16) influenza season. Thus, active surveillance of circulating influenza viruses in the Southern Asia region, including Nepal and India, would be beneficial for detecting novel variant viruses prior to their worldwide spread

Apollo (Cambridge)

Influenza Other Respi Viruses. 2017;11:399–403.	 	 	 | 	399wileyonlinelibrary.com/journal/irv Accepted: 11 July 2017DOI: 10.1111/irv.12469S H O R T  A R T I C L ECharacterization of influenza A(H1N1)pdm09 viruses isolated from Nepalese and Indian outbreak patients in early 2015Kazuya Nakamura1  | Masayuki Shirakura1 | Seiichiro Fujisaki1 | Noriko Kishida1 |  David F. Burke2 | Derek J. Smith2 | Tomoko Kuwahara1 | Emi Takashita1 |  Ikuyo Takayama1 | Mina Nakauchi1 | Mandeep Chadha3 | Varsha Potdar3 |  Arvind Bhushan3 | Bishnu Prasad Upadhyay4 | Geeta Shakya4 | Takato Odagiri1 |  Tsutomu Kageyama1 | Shinji Watanabe1This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.© 2017 The Authors. Influenza and Other Respiratory Viruses Published by John Wiley & Sons Ltd.1Influenza Virus Research Center, National Institute of Infectious Diseases, Tokyo, Japan2Center for Pathogen Evolution, University of Cambridge, Cambridge, UK3Influenza group, National Institute of Virology, Indian Council of Medical Research, Pune, India4Department of Health Services, National Public Health Laboratory, Kathmandu, NepalCorrespondenceShinji Watanabe, Influenza Virus Research Center, National Institute of Infectious Diseases, Tokyo, Japan.Email: sw@nih.go.jpFunding informationMinistry of Health, Labor and Welfare, Japan; JSPS KAKENHI, Grant/Award Number: JP26460816; NIAID, Grant/Award Number: HHSN272201400008C; BARDA, Grant/Award Number: HHSO100201500033CWe characterized influenza A(H1N1)pdm09 isolates from large- scale outbreaks that occurred in Nepal and India in early 2015. Although no specific viral features, which may have caused the outbreaks, were identified, an S84N substitution in hemaggluti-nin was frequently observed. Chronological phylogenetic analysis revealed that these Nepalese and Indian viruses possessing the S84N substitution constitute potential an-cestors of the novel genetic subclade 6B.1 virus that spread globally in the following (2015/16) influenza season. Thus, active surveillance of circulating influenza viruses in the Southern Asia region, including Nepal and India, would be beneficial for detecting novel variant viruses prior to their worldwide spread.K E Y W O R D SIndia, influenza A(H1N1)pdm09 virus, Nepal, outbreakDespite the worldwide predominance of the A(H3N2) virus in the 2014/15 season, India reported a huge outbreak of influenza A(H1N1)pdm09 virus (H1N1pdm) infection with severe clinical outcomes at the beginning of 2015; over 30 000 cases were reported with over 2000 deaths, compared with 5044, 5250 and <1000 reported cases in 2012, 2103, and 2014, respectively.1-3 Subsequently, Nepal also expe-rienced a regional outbreak with an increased number of patients in-fected with H1N1pdm (personal communication).4 Although a number of Indian case reports have proposed several possible causes for the outbreak and characterized viruses isolated from the patients,2,5,6 the viral properties of isolates from outbreak patients remain largely un-known. Thus, in this study, the antigenicity, genetic properties, and an-tiviral susceptibility of 43 Nepalese isolates from 50 clinical specimens of Nepalese outbreak patients (patient information is summarized in Table S1) and 10 Indian isolates were characterized to identify particu-lar features related to outbreaks with severe and fatal cases.The antigenicity of these Nepalese and Indian isolates was charac-terized using the hemagglutination inhibition (HI) test. The hemagglu-tination activity of all isolates tested was inhibited by ferret antisera raised against several reference viruses including the 2014/15 sea-son vaccine strain, egg- propagated A/California/07/2009, the high- growth reassortant strain, A/California/07/2009 (X- 179A), rep-resentative field isolates, cell- propagated A/Narita/1/2009, and A/Wakayama/153/2013 or A/Yokohama/50/2015. The HI titers of these antisera against the Nepalese and Indian isolates were equiv-alent to or within a fourfold reduction compared to the homologous 400  |     NAKAMURA et Al.titer of each reference virus (data not shown), indicating that all Nepalese and Indian isolates tested were antigenically similar to the contemporary circulating H1N1pdm.The susceptibility of the isolates to four neuraminidase (NA) inhibi-tors, oseltamivir, peramivir, zanamivir, and laninamivir, was determined using a fluorescent NA inhibition assay with the NA- Fluor influenza neuraminidase assay kit (Applied Biosystems, Foster City, CA, USA), as previously described.7 The results showed that the Nepalese and Indian isolates were susceptible to all four inhibitors (data not shown), suggesting that no antiviral- resistant viruses had spread in these outbreaks.Taken together, our viral antigenicity and antiviral susceptibility analyses did not identify any unusual features in the Nepalese and Indian outbreak isolates.Genetic sequencing and phylogenetic analysis were conducted for all Nepalese and Indian isolates. As shown in Figure 1, the isolates contain a hemagglutinin (HA) gene belonging to genetic clade 6B represented by globally circulating H1N1pdm.8 Interestingly, 70% of the Nepalese viruses and all Indian viruses, except for one (A/India/P1510025/2015), formed a subclade branch with a high bootstrap value (92/100) due to an S84N amino acid substitution in their HA (designated as subclade 1 in Figure 1). Of the remaining viruses an-alyzed (15 Nepalese viruses and one Indian virus), nine Nepalese vi-ruses and one Indian virus formed a small cluster (94/100 bootstrap value). Within this small cluster, five Nepalese viruses formed a sub-cluster with an E506D amino acid substitution (88/100 bootstrap value; designated as subclade 2 in Figure 1). The other six Nepalese viruses belonged to a different cluster possessing an E491G substi-tution; within this cluster, five viruses formed a subcluster exhibit-ing an N129D substitution (87/100 of bootstrap value; designated as subclade 3 in Figure 1). Analysis of the NA gene showed that the Nepalese and Indian viruses possessing an S84N substitution in HA formed a subclade exhibiting V13I and I314M substitutions in NA (designated as subclade 1 in Fig. S1). The remaining viruses formed NA subclades related to the HA subclades (designated as subclade 2 and subclade 3, respectively, in Fig. S1). Although the biological sig-nificance of these amino acid substitutions in HA and NA has yet to be precisely elucidated, they seem to be correlated with a certain ex-tent in terms of phylogenetic tree (sub)clade structure, suggesting the coevolution of HA and NA with viral fitness.9 The Indian isolate, A/India/P152122/2015, possessed a D222G/N amino acid substitution in its HA. This polymorphism in HA has been occasionally associated with viral pathogenesis.10 However, because this is the only isolate analyzed in this study that had this substitution, it is unlikely that this substitution was solely responsible for the outbreaks. The nucleotide sequences of the internal genes (PA, PB1, PB2, NP, M, and NS) of all Indian and seven randomly selected Nepalese isolates were deter-mined using next- generation sequencing with the Illumina MiSeq sys-tem (Illumina K. K., Tokyo, Japan). Genetic and phylogenetic analysis of these genes revealed that most of the Nepalese and Indian isolates were clustered into a genetic subclade separated from other subclades formed by globally circulating H1N1pdm (Fig. S2). The amino acid sub-stitutions identified in the internal genes analyzed are summarized in Table S2. Although several amino acid substitutions accumulated in these isolates compared to the globally circulating H1N1pdm, no known amino acid substitutions related to viral pathogenicity were identified. In addition, there were no characteristic differences in the genetic features of the internal genes between the patients with severe and mild outcomes (Table S2). Sequence data determined in this study were deposited with the GISAID database under accession numbers EPI610408 to EPI610535, EPI630406 to EPI630425, and EPI635481 to EPI635540.In the following influenza season (2015/16), H1N1pdm infections were more prevalent than A/H3N2 and type B virus infections in the Northern Hemisphere temperate regions, as reported by GISRS.11 Most of these H1N1pdm isolates reacted well with the ferret antisera described above (data not shown), indicating that no apparent anti-genic change had occurred in the circulating H1N1pdm strains. To de-termine the genetic correlation between these predominate 2015/16 season H1N1pdm strains and the Nepalese and Indian isolates, further phylogenic analysis with a large number of H1N1pdm strains was con-ducted (Figure 2). The vast majority of the 2015/16 season H1N1pdm strains (designated by a blue vertical line in Figure 2) formed part of a large cluster possessing an S84N substitution in HA (isolates col-ored in light blue in Figure 2). These isolates were descended from the Nepalese and Indian isolates described above (designated by a green vertical line in Figure 2). Most of these 2015/16 season isolates con-tained S162N (on the Sa antigenic site of HA resulting in a potential glycosylation site) and I216T amino acid substitutions in addition to the S84N substitution and formed a subclade subsequently desig-nated as 6B.1.12 Intriguingly, this set of S84N variants was preceded by the emergence of isolates possessing an S84N substitution in Nepal in March and April 2014 (designated by a red rectangle in Figure 2 and enlarged in Fig. S3). Although isolates with an S84N substitution were occasionally observed earlier than 2014 in other parts of the world (colored in light blue from 2009 to 2013 in Figure 2), they disap-peared shortly thereafter for unknown reasons. The 2014/15 season Nepalese and Indian isolates possessing an S84N substitution origi-nated from isolates with S84N that emerged in Nepal in early 2014 and contributed to further worldwide spread (2015/16 season), sug-gesting the importance of these areas as a source for the evolution and global spread of novel variant viruses.F IGURE  1 Phylogenetic analysis of the hemagglutinin gene of H1N1pdm strains. Nucleotide sequence data of the former representative H1N1pdm strains and recent 2014/15 season H1N1pdm isolates from around the world were randomly selected from the GISAID database and used as the references. The phylogenetic trees were constructed using the MEGA6 software with the neighbor- joining method. Nepalese and Indian strains analyzed in this study are colored in blue and red, respectively. Open squares indicate that only genetic sequence data were available in this study because of virus isolation failure. Previously reported8 strain A/Fukuoka/SDC1/2015 is underlined. Amino acid substitutions relative to the A/California/07/2009 virus and high bootstrap values (>80) are shown to the left of the nodes. The scale bar indicates the nucleotide substitutions per site. The GISAID EpiFlu database accession numbers for the sequence data of each reference strain are shown     |  401NAKAMURA et Al.402  |     NAKAMURA et Al.     |  403NAKAMURA et Al.Although no viral features related to expanded infection with severe and fatal outcomes, including antigenicity, genetic properties, or suscepti-bility to antivirals, were identified in this study, it is important to note that many Nepalese and Indian H1N1pdm isolates possessing a novel genetic feature (S84N substitution in HA) preceded strains with additional sub-stitutions observed in the subsequent annual influenza season (2015/16 season) in the Northern Hemisphere temperate regions. These novel H1N1pdm variants have also caused widespread outbreaks with severe or fatal outcomes.11 The phylogenetic analyses in our current study sug-gest that South Asian countries, such as Nepal and India, might serve as an important source of variant viruses, and improve our understanding of the evolutionary dynamics of circulating influenza viruses. Thus, the future establishment of systematic measures for monitoring circulating influenza viruses in these countries would be of great importance.ACKNOWLEDGEMENTSWe wish to thank H. Sugawara, A. Sato, R. Ogawa, H. Miura, and M. Akimoto for technical assistance with this study; D. T. Mourya (National Institute of Virology, Pune) for providing support through-out the experiments; and P. Koul (Sher- i- Kashmir Institute of Medical Sciences, Sri Nagar), L. Dar (All India Institute of Medical Sciences, New Delhi), and M. Chawla Sarkar (National Institute of Cholera and Enteric Diseases, Kolkata) for submitting influenza isolates/clinical samples to the National Institute of Virology, NIC India. In addition, we wish to thank the various laboratories that compiled and submitted (GISAID EpiFlu database) the sequences used in this study and the Influenza Virus Surveillance Group of Japan for providing Japanese isolates. Moreover, we wish to thank Editage (www.editage.jp) for English language editing. This work was supported, in part, by a Grant- in- Aid for Emerging and Re- emerging Infectious Diseases from the Ministry of Health, Labor and Welfare, Japan, as well as by JSPS KAKENHI grant JP26460816. DFB and DJS wish to acknowledge the support re-ceived from the Center for Research on Influenza Pathogenesis (CRIP) funded by the NIAID and BARDA contracts HHSN272201400008C and HHSO100201500033C, respectively.CONFLICT OF INTERESTThe authors declare that no conflict of interest exists.REFERENCES 1. Bagchi S. India tackles H1N1 influenza outbreak. Lancet. 2015;385:e21. 2. Mishra B. 2015 resurgence of influenza A (H1N1) 09: smoldering pan-demic in India? J Glob Infect Dis. 2015;7:56-59. 3. D’Silva J. Swine flu: how well did India respond? BMJ. 2015;350:h2286. 4. World Health Organization. Influenza surveillance report. https://extranet.who.int/sree/Reports?op=vs&path=/WHO_HQ_Reports/G5/PROD/EXT/Influenza%20Surveillance+Report+by+Country. Accessed December 10, 2015. 5. Tharakaraman K, Sasisekharan R. Influenza surveillance: 2014- 2015 H1N1 “swine”- derived influenza viruses from India. Cell Host Microbe. 2015;17:279-282. 6. Takashita E, Fujisaki S, Shirakura M, et al. Characterization of an A(H1N1)pdm09 virus imported from India in March 2015. Jpn J Infect Dis. 2016;69:83-86. 7. Takashita E, Kiso M, Fujisaki S, et al. Characterization of a large clus-ter of influenza A(H1N1)pdm09 viruses cross- resistant to oseltami-vir and peramivir during the 2013–2014 influenza season in Japan. Antimicrob Agents Chemother. 2015;59:2607-2617. 8. World Health Organization. Recommended composition of influenza virus vaccines for use in the 2015–2016 northern hemisphere influ-enza season. Wkly Epidemiol Rec. 2015;90:97-108. 9. Wagner R, Matrosovich M, Klenk H-D. Functional balance between haemagglutinin and neuraminidase in influenza virus infections. Rev Med Virol. 2002;12:159-166. 10. Resende PC, Motta FC, Oliveira M, et al. Polymorphisms at residue 222 of the hemagglutinin of pandemic influenza A(H1N1)pdm09: as-sociation of quasi- species to morbidity and mortality in different risk categories. PLoS ONE 2014;9:e92789. 11. World Health Organization. Risk assessment - seasonal influenza A(H1N1)pdm09. http://www.who.int/influenza/publications/riskas-sessment_AH1N1pdm09_201602/en/. Accessed May 1, 2016. 12. Broberg E, Melidou A, Prosenc K, Bragstad K, Hungnes O, WHO European Region and the European Influenza Surveillance Network members of the reporting countries. Predominance of influenza A(H1N1)pdm09 virus genetic subclade 6B.1 and influenza B/Victoria lineage viruses at the start of the 2015/16 influenza season in Europe. Euro Surveill. 2016;21:pii=30184. 13. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post- analysis of large phylogenies. Bioinformatics. 2014;30: 1312-1313.SUPPORTING INFORMATIONAdditional Supporting Information may be found online in the supporting information tab for this article.     How to cite this article: Nakamura K., Shirakura M., Fujisaki S., et al. (2017), Characterization of influenza A(H1N1)pdm09 viruses isolated from Nepalese and Indian outbreak patients in early 2015. Influenza Other Respi Viruses 11, 399-403.  https://doi.org/10.1111/irv.12469F IGURE  2 Phylogenetic relationships of the hemagglutinin (HA) genes of Nepalese and Indian isolates and H1N1pdm isolates since 2009. Nucleotide sequence data of H1N1pdm isolates from around the world were obtained from the GISAID database and aligned using MAFFT (available from http://www.ebi.ac.uk/mafft/). MODELTEST was run on the aligned sequences, and GTR+I+gamma4 was determined as the most appropriate evolutionary model for phylogenetic tree construction. The tree was constructed using RAxML v813 according to the GTR+I+gamma4 model. Global optimization of branch topology was performed on the tree with the best likelihood score from RAxML using Garli v0.951 with a run of 150 000 generations. Isolates possessing an S84N substitution in HA are colored in light blue. Isolates possessing S84, S84G, and a mixed population of amino acids are colored in orange, gray, and black, respectively

English

Characterization of influenza A(H1N1)pdm09 viruses isolated from Nepalese and Indian outbreak patients in early 2015

https://www.repository.cam.ac.uk/bitstream/1810/267403/3/Nakamura_et_al-2017-Influenza_and_Other_Respiratory_Viruses.pdf

Characterization of influenza A(H1N1)pdm09 viruses isolated from Nepalese and Indian outbreak patients in early 2015

Abstract

Similar works

Full text

Available Versions

Apollo (Cambridge)