Re-Invasion of H5N8 High Pathogenicity Avian Influenza Virus Clade 2.3.4.4b in Hokkaido, Japan, 2020

Abstract

Global dispersion of high pathogenicity avian influenza (HPAI), especially that caused by H5 clade 2.3.4.4, has threatened poultry industries and, potentially, human health. An HPAI virus, A/northern pintail/Hokkaido/M13/2020 (H5N8) (NP/Hok/20) belonging to clade 2.3.4.4b, was isolated from a fecal sample collected at a lake in Hokkaido, Japan where migratory birds rested, October 2020. In the phylogenetic trees of all eight gene segments, NP/Hok/20 fell into in the cluster of European isolates in 2020, but was distinct from the isolates in eastern Asia and Europe during the winter season of 2017–2018. The antigenic cartography indicates that the antigenicity of NP/Hok/20 was almost the same as that of previous isolates of H5 clade 2.3.4.4b, whereas the antigenic distances from NP/Hok/20 to the representative strains in clade 2.3.4.4e and to a strain in 2.3.4 were apparently distant. These data imply that HPAI virus clade 2.3.4.4b should have been delivered by bird migration despite the intercontinental distance, although it was not defined whether NP/Hok/20 was transported from Europe via Siberia where migratory birds nest in the summer season. Given the probability of perpetuation of transmission in the northern territory, periodic updates of intensive surveys on avian influenza at the global level are essential to prepare for future outbreaks of the HPAI virus.

1. Introduction

Avian influenza is a highly contagious disease that is caused by infection with influenza A viruses, categorized into highly or low pathogenic influenza viruses according to its pathogenicity in chickens [1]. Influenza A viruses originating from wild birds are subtyped by reactivity to antisera against two glycoproteins, hemagglutinin (HA) and neuraminidase (NA), on the virus surface, to H1-16 and N1-9, respectively [2]. High pathogenicity avian influenza (HPAI) viruses are limited to either the H5 or H7 subtypes and possess multiple basic amino acid residues at the site of HA proteolytic cleavage, leading to a systemic infection with ubiquitous protease in the avian hosts [3].

In 1997, a large-scale HPAI outbreak occurred in Hong Kong leading to 18 human infections with 3 deaths and culling of 1.2 million poultry [4]. The causal agent of the 1997 HPAI outbreak was the H5N1 HPAI virus (HPAIV), which originated from the prevalent strain in domestic geese of Guangdong Province, China, in 1996 [5]. Following continuous HPAI outbreaks in Asia, the introduction of H5N1 HPAIV by migratory birds, which are natural reservoirs of the avian influenza virus, could contribute to an interregional spread of infectious viruses [6]. Through the worldwide spread of H5 HPAIV, the H5 HA strains have evolved from an isolated strain in Guangdong, 1996, into 10 genetically distinct virus clades and multiple subclades [7], in which the antigenicity of H5 HPAIVs has varied accordingly [8,9]. In particular, the emergence of HPAIV H5 HA clade 2.3.4.4, which evolved from the prevalent strain within China (clade 2.3.4) was recognized not only in eastern Asia due to genetic reassortment of other H5 HPAIV clades and local low pathogenic avian influenza viruses (LPAIV), but also in western Asia, Europe, and North America [10,11,12,13,14,15,16,17]. As of end of 2019, nomenclature of the H5 HA was updated to divide the group a to h due to further genetic evolvement confirmed among the recent isolates [18].

Circulation of various groups of HPAIV necessitates an intensive surveillance of avian influenza virus among wild and migratory birds to provide an update on the virus circulation situation so that proper measures can be implemented to mitigate its invasion. In October 2020, we isolated an HPAIV H5N8 clade 2.3.4.4b from a fecal sample of a migratory bird that usually nests in the northern territory, collected at Lake Komuke in the northeastern part of Hokkaido, Japan. The objective of the present study was to characterize the genetic and antigenic properties of the novel isolate of A/northern pintail/Hokkaido/M13/2020 (H5N8) (NP/Hok/20). 2. Materials and Methods

A total of 100 duck fecal samples were collected at the lakeside of Lake Komuke (latitude: 44°25′47″ north; longitude: 143°50′25″ east) in the eastern part of Hokkaido, Japan on 24 October 2020 (Figure 1). Each of the fecal samples mixed with Minimum Essential Medium (Nissui, Tokyo, Japan) containing antibiotics was inoculated into 10-day-old embryonated chicken eggs for virus isolation [19]. After virus propagation in the allantoic fluid, the influenza virus subtype was determined by hemagglutination inhibition (HI) and neuraminidase inhibition (NI) tests using antisera against the referenced influenza virus strain for subtyping [20].

DNA was extracted from the fecal sample using a DNA Mini Kit (Qiagen, Hilden, Germany), according to the manual instructions. A 749 bp region near the 5’ terminus of the cytochrome c oxidase I (COI) gene of a host bird was amplified using primers described in the previous study [21]. DNA sequences of the products of polymerase chain reaction were compared with the sequence data of the COI gene obtained from the Barcode of Life Data Systems (http://www.barcodinglife.com).

Viral RNA was extracted from the allantoic fluid of embryonated chicken eggs using TRIzol LS Reagent (Life Technologies, Carlsbad, CA, USA), and then was amplified by PCR reaction using each of 8 segment gene-specific primer sets for each of the 8 segments after reverse transcription [22]. Direct sequencing of each gene was performed using BigDye Terminator version 3.1, a Cycle Sequencing Kit (Life Technologies), and a 3500 Genetic Analyzer (Life Technologies). The nucleotide sequence datasets of 8 gene segments from the representative H5N8 clade 2.3.4.4 strains, or HxN8, including outgroup viruses downloaded from the Global Initiative on Sharing All Influenza and GenBank, were prepared with adequate alignment for constructing phylogenetic trees. The maximum-likelihood method was applied to construct the phylogenetic tree of each gene segment using the best-fit general time-reversible model of the nucleotide substitution with gamma-distribution rate variation among sites (with four rate categories, Γ) in MEGA 10 software [23]. The definition of the H5 HA clade in the present study followed the update of nomenclature for phylogenetic relationships, which had been scheduled at the end of 2019 [18].

Cross-reactivity of the antisera and their corresponding antigens was assessed using the HI test (Table S1). The A/chicken/Kumamoto/1-7/2014 (H5N8) (Ck/Kum/1-7/14) was isolated from a chicken in Japan. The A/Muscovy duck/DR Congo/KAF1/2017 (H5N8) (MDk/CD/KAF1/17) was isolated from a Muscovy duck in DR Congo. The A/black swan/Akita/1/2016 (H5N6) (BS/Aki/1/16) was isolated from a black swan in Japan. The A/duck/Vietnam/HU1-1151/2016 (H5N6) (Dk/VTN/1151/16) was isolated from a Pekin duck in Vietnam. The A/peregrine falcon/Hong Kong/810/2009 (H5N1) (PF/HK/810/09) was isolated from a peregrine falcon in Hong Kong. We also included in the panel antiserum against A/mallard duck/Hokkaido/24/2009 (H5N1) (Mal/Hok/24/09), which is an unclassified LPAIV isolated from wild birds in Japan [19]. An antigenic cartograph based on the results of the HI test was estimated using web-based software [24]. The data of the results of cross-HI titers were loaded to obtain the x/y coordinates of each antiserum and antigen.

3. Results

Among the 100 fecal samples collected at Lake Komuke, two influenza A viruses were isolated, one of which was H14N2, the other was subtyped as H5N8 by the HI and NI tests and designated as A/northern pintail/Hokkaido/M13/2020 (H5N8) (NP/Hok/20). The sequence analyses revealed that the HA of NP/Hok/20 had multiple basic amino residues at the proteolytic cleavage site, REKRRKR/G, which is defined as a molecular marker of HPAIV [25]. In the COI gene database, the sequence of mitochondrial DNA extracted from the fecal sample showed 100% homology to the COI gene of Anas acuta (northern pintail). The embryonated chicken eggs and two chickens intranasally inoculated with NP/Hok/20 died within 48 to 72 h post-inoculation, respectively (data not shown), suggesting that NP/Hok/20 is an HPAIV.

The full-length nucleotide sequences for all eight gene segments of NP/Hok/20 were analyzed with those of other H5 HPAIVs and some N8 LPAIVs for a phylogenetic analysis of the NA (Figure 2 and Figure S1). All eight gene segments of NP/Hok/20 were completely sequenced and submitted to GenBank (accession numbers: MW228062-MW228069) and were confirmed to be almost identical to strains isolated in Korea, October 2020, or isolated in the southern part of Japan (Kagoshima), November 2020. The phylogenetic tree of the HA segment indicated that the HAs of NP/Hok/20 and the Korean strain were classified in clade 2.3.4.4b and further classified into the same cluster of European isolates in the winter season of 2019−2020, including A/chicken/Slovakia-Pah 14/2020 (H5N8), with which NP/Hok/20 showed the third highest homology (Figure 2A). Within 2.3.4.4b, on the other hand, the cluster containing NP/Hok/20 was apparently distinct both from the H5N6 isolates of eastern Asia and Europe in the winter season of 2017–2018, as well as H5Nx Europe strains isolated in October and November 2020. The phylogenetic tree of other gene segments except N8 NA indicated similar results as that of H5 HA; NP/Hok/20 was classified into the same cluster of European isolates in 2020 (Figure S1). In the phylogenetic tree of the NA gene segment, the NA of NP/Hok/20 was classified into the cluster of H5N8 clade 2.3.4.4b strains isolated in Europe in the winter season of 2019–2020, although it was distinct from the cluster of clade 2.3.4.4c (Figure 2B).

The antigenicity of NP/Hok/20 was characterized and compared with that of other groups of clade 2.3.4.4, using an HI test (Table S1). The HI titer of NP/Hok/20 against the antisera of Ck/Kum/1-7/14, MDk/CD/KAF1/17, Dk/VTN/1151/14, and PF/HK/810/09 had the same antigenicity as the HI titers of each homologous antigen and antiserum combination. On the other hand, NP/Hok/20 did not react well to the antisera of BS/Aki/1/16 and Mal/Hok/24/09. The antigenic cartograph developed using the results of the cross-HI tests of Table S1 demonstrated that the clade 2.3.4.4b, including NP/Hok/20, was antigenically close to clade 2.3.4.4c (Figure 3). The isolates of clade 2.3.4.4b would partially react to clade 2.3.4.4a, but were more than two antigenic units away, which were estimated to be the equivalent of a 4-fold HI difference from that of clade 2.3.4.4e.

4. Discussion

In October 2020, an HPAIV, A/northern pintail/Hokkaido/M13/2020 (H5N8) was isolated from a fecal sample of the migratory bird in the northern part of Japan, and a series of HPAI outbreaks in poultry farms have been confirmed in Japan since November 2020 [26]. The results of the sequence analysis indicated that H5 HA of NP/Hok/20 was classified into clade 2.3.4.4b, but into a separate cluster of isolates in eastern Asia during the winter season of 2017–2018. Moreover, the high sequence homology of NP/Hok/20 to the isolates in Europe during the winter season of 2019–2020 implies that NP/Hok/20 was delivered from Europe via bird migration, rather than as a product of genetic reassortment between the HPAI H5 clade 2.3.4.4b strain, which was potentially maintained in East Asia, and other local LPAIVs. Given that dissemination of identical viruses in two far-distinct regions was confirmed, as well as genetically distinct HPAIVs from NP/Hok/20 being prevalent in Europe at the end of 2020, several kinds of highly contagious viruses could have already been perpetuated in northern territories, like Siberia, though this assumption should be highly influenced by the isolation of a virus possessing similar genetic characteristics as NP/Hok/20 in there. It has been understood that influenza viruses perpetuate in aquatic birds including ducks in northern territories like Siberia or Alaska where they usually nest in summer [27,28]. Eastern Asia has been a hot spot for HPAIV outbreaks caused by H5 HPAIV infection for more than 20 years, since most of the new lineage strains and genetic reassortments were produced within or near this region. Japan, like other nations, has tolerated the introduction of the H5 clade 2.3.4.4 by bird migration [29,30,31]. The waterfowl migration pathways satellite-tracked between eastern Asia and eastern Russia support the assumption of HPAIV introduction from eastern Siberia in the fall−winter season and from eastern Asia in the winter−spring season [32,33]. Moreover, a satellite-tracked migration investigation could not confirm the direct migration between central Asia and Hokkaido, Japan, but could verify migration between the eastern part of Hokkaido Japan and the Kamchatka Peninsula [33]. It may indicate that contagious pathogens that were originally isolated in Europe could be transferred via the northern territory to eastern Asia, and could also be carried to the northern part of North America. These speculations would warn us of the likely perpetuation of HPAIV in the north and continuous invasion of HPAIV through bird migration to the south.

Antigenicity variance among different HA lineages is well known. The antigenicity of the African isolate in clade 2.3.4.4b was neither different from the isolate in clade 2.3.4.4a nor that in 4c, whereas a moderate difference observed via antigenic cartography was confirmed between isolates in clades 2.3.4.4b and 4e [34]. The results obtained in the present study are similar, and the antigenicity of between two clade 2.3.4.4b isolates, one from the Democratic Republic of Congo and one from Japan, would almost be the same (Figure 2). On the other hand, expansion of the clusters, including European and Japanese isolates in 2020, from the clusters of past isolates in the H5 HA phylogenetic tree reveals another concern. This finding implied that further genetic variants were selected from strains currently circulating among migratory birds. Further investigation on the H5 clade 2.3.4.4 is required to monitor the potential genetical and antigenic changes.

An HPAIV H5N8 clade 2.3.4.4b virus was isolated from the fecal sample of migratory bird in Hokkaido, Japan, in October 2020, and was almost identical to European isolates in the winter season of 2019–2020 but was distinct from the European isolates in the end of 2020. As the perpetuation of contagious viruses in northern territories is likely, intensive preparedness for their invasion should be implemented. The possibility of intercontinental HPAIV transportation could admonish us to strengthen the preparedness for its invasion. Enforcement of biosecurity measures on poultry farms and of the exclusion of wild birds from the vicinity, and the monitoring of avian influenza virus circulation in both domestic poultry and wild migratory birds, even in previously unaffected areas, are essential to minimize HPAI threats. Indeed, in Japan, between 30 October and 8 December 2020, a total of 19 poultry outbreaks and eight wild bird cases of HPAI were reported in six and three prefectures, respectively [35]. In parallel, strengthening the function of the global network for related influenza viruses is an urgent agenda so that the genetic and antigenic properties of H5N8 isolates can be used to develop a more effective control strategies against variable strains of H5N8 clade 2.3.4.4.

Figure 1. Location of a fecal sample collection. Fecal samples were collected at Lake Komuke, northeastern Hokkaido, in Japan in October 2020.

Figure 2. Phylogenetic tree analysis on the hemagglutinin (HA) and neuraminidase (NA) gene segments. Phylogenetic tree analysis was based on (A) the HA gene segment containing a total of 67 H5 clade 2.3.4.4, other clades, and outgroups, and (B) the NA gene segment containing a total of 49 N8 NA, including highly and low pathogenic avian influenza viruses. The red-colored strain indicates H5N8 high pathogenicity avian influenza viruses (HPAIVs) isolated in Japan and Korea, 2000 of which A/northern pintail/Hokkaido/M13/2020 (H5M8), an isolate in Hokkaido, 2020 is underlined. The green-, yellow-, and blue-colored strains indicate the H5N8 HPAIVs isolated in Europe in the winter of 2019-2020, H5N6 HPAIVs isolated in Japan and Korea in the winter season of 2017-2018, and H5Nx HPAIVs isolated in Europe in October and November, 2020, respectively. The numbers below or above the node indicate bootstrap values greater than 60%.

Figure 2. Phylogenetic tree analysis on the hemagglutinin (HA) and neuraminidase (NA) gene segments. Phylogenetic tree analysis was based on (A) the HA gene segment containing a total of 67 H5 clade 2.3.4.4, other clades, and outgroups, and (B) the NA gene segment containing a total of 49 N8 NA, including highly and low pathogenic avian influenza viruses. The red-colored strain indicates H5N8 high pathogenicity avian influenza viruses (HPAIVs) isolated in Japan and Korea, 2000 of which A/northern pintail/Hokkaido/M13/2020 (H5M8), an isolate in Hokkaido, 2020 is underlined. The green-, yellow-, and blue-colored strains indicate the H5N8 HPAIVs isolated in Europe in the winter of 2019-2020, H5N6 HPAIVs isolated in Japan and Korea in the winter season of 2017-2018, and H5Nx HPAIVs isolated in Europe in October and November, 2020, respectively. The numbers below or above the node indicate bootstrap values greater than 60%.

Figure 3. An antigenic map based on the cross-hemagglutination inhibition (HI) tests on H5 clade 2.3.4.4 viruses and the sera of different genetic groups. The antigenic relationship among H5 clade 2.3.4.4 viruses, including NP/Hok/20, was visualized using antigenic cartography based on the results of the cross-HI tests, as indicated in Table S1. In an antigenic map, both vertical and horizontal axes represent the antigenic distance, and the spacing between grid lines represents a distance of 1 antigenic-unit distance, which corresponds to a 2-fold dilution in the HI assay.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/12/1439/s1, Figure S1. Phylogenetic tree analysis on the internal gene segments, Table S1. Cross-reactivity of H5 Highly pathogenic avian influenza virus clade 2.3.4.4.

Author Contributions

Conceptualization, N.I., K.S. and Y.S.; sample collection, K.O., Y.W. and K.S.; methodology, N.I. and Y.S.; validation, N.I., A.T.T., E.B. and Y.S.; formal analysis, A.T.T., E.B. and Y.S. investigation, A.T.T., E.B., H.H., Z.-J.W. and D.K.; resources, Y.S.; data curation, N.I., A.T.T., E.B. and H.H.; writing-original draft preparation, N.I., Z.-J.W., and Y.S.; writing-review and editing, Y.S. and H.K.; visualization, N.I.; supervision, Y.S. and H.K.; project administration, Y.S.; funding acquisition, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was mainly performed by the Environment Research and Technology Development, grant number JPMEERF18S20103. This work was partially supported by the e-ASIA joint research program, funded by the Japan Agency for Medical Research and Development (AMED) grant number JP20jm0210054h0004. The AMED and Japan International Cooperation Agency (JICA) also partially supported this work within the framework of the Science and Technology Research Partnership for Sustainable Development (SATREPS) grant number JP20jm0110019.

Acknowledgments

We are greatly thankful to Kazue OKA, Laboratory of Microbiology, Graduate School of Veterinary Medicine, Hokkaido University, and Hibiki MORIMOTO, Hida Livestock Hygiene Service Center, Gifu Prefecture, for kind contribution to obtain the sequence data presented in the present manuscript. The present work was partially supported by Doctoral Program for World-leading Innovative & Smart Education (WISE), by the Ministry of Education, Culture, Sports, Science and Technology (MEXT). Our great appreciation is extended to Youn-Jeong Lee, Animal and Plant Quarantine Agency, Korea, for providing the beneficial information regarding the HPAIV isolated in Korea. We are greatly thankful to Makoto Ozawa, Joint Faculty of Veterinary Medicine, Kagoshima University for providing the beneficial information regarding the HPAIV isolated in Kagoshima as well as the method of mitochondrial DNA sequencing.

Conflicts of Interest

The funders did not have any influence on the analysis or interpretation of the results in this present study. Authors declare no conflict of interest.

References

1. Alexander, D.J. A review of avian influenza in different bird species. Vet. Microbiol. 2000, 74, 3-13.

2. Howley, M.; Knipe, D. Orthomyxoviruses. In Fields Virology: Emerging Viruses, 7th ed.; Whelan, S., Kluwer, W., Eds.; Wolters Kluwer Health Inc.: Philadelphia, PA, USA, 2021; Volume 1, p. 650.

3. Klenk, H.D.; Garten, W. Host cell proteases controlling virus pathogenicity. Trends Microbiol. 1994, 2, 39-43.

4. Shortridge, K.F. Poultry and the influenza H5N1 outbreak in Hong Kong, 1997: Abridged chronology and virus isolation. Vaccine 1999, 17 (Suppl. 1), S26-S29.

5. Xu, X.; Subbarao, K.; Cox, N.J.; Guo, Y. Genetic characterization of the pathogenic influenza A/Goose/Guangdong/1/96 (H5N1) virus: Similarity of its hemagglutinin gene to those of H5N1 viruses from the 1997 outbreaks in Hong Kong. Virology 1999, 261, 15-19.

6. Wang, G.; Zhan, D.; Li, L.; Lei, F.; Liu, B.; Liu, D.; Xiao, H.; Feng, Y.; Li, J.; Yang, B.; et al. H5N1 avian influenza re-emergence of Lake Qinghai: Phylogenetic and antigenic analyses of the newly isolated viruses and roles of migratory birds in virus circulation. J. Gen. Virol. 2008, 89, 697-702.

7. WHO; OIE; FAO; Evolution Working Group. Toward a unified nomenclature system for highly pathogenic avian influenza virus (H5N1). Emerg. Infect. Dis. 2008, 14, e1.

8. Ohkawara, A.; Okamatsu, M.; Ozawa, M.; Chu, D.H.; Nguyen, L.T.; Hiono, T.; Matsuno, K.; Kida, H.; Sakoda, Y. Antigenic diversity of H5 highly pathogenic avian influenza viruses of clade 2.3.4.4 isolated in Asia. Microbiol. Immunol. 2017, 61, 149-158.

9. Yang, L.; Zhu, W.; Li, X.; Bo, H.; Zhang, Y.; Zou, S.; Gao, R.; Dong, J.; Zhao, X.; Chen, W.; et al. Genesis and Dissemination of Highly Pathogenic H5N6 Avian Influenza Viruses. J. Virol. 2017, 91.

10. The Global Consortium for H5N8 and Related Influenza Viruses. Role for migratory wild birds in the global spread of avian influenza H5N8. Science 2016, 354, 213-217.

11. Gu, M.; Liu, W.; Cao, Y.; Peng, D.; Wang, X.; Wan, H.; Zhao, G.; Xu, Q.; Zhang, W.; Song, Q.; et al. Novel reassortant highly pathogenic avian influenza (H5N5) viruses in domestic ducks, China. Emerg. Infect. Dis. 2011, 17, 1060-1063.

12. Lee, D.H.; Bahl, J.; Torchetti, M.K.; Killian, M.L.; Ip, H.S.; DeLiberto, T.J.; Swayne, D.E. Highly Pathogenic Avian Influenza Viruses and Generation of Novel Reassortants, United States, 2014-2015. Emerg. Infect. Dis. 2016, 22, 1283-1285.

13. Lee, D.H.; Torchetti, M.K.; Winker, K.; Ip, H.S.; Song, C.S.; Swayne, D.E. Intercontinental Spread of Asian-Origin H5N8 to North America through Beringia by Migratory Birds. J. Virol. 2015, 89, 6521-6524.

14. Lee, Y.J.; Kang, H.M.; Lee, E.K.; Song, B.M.; Jeong, J.; Kwon, Y.K.; Kim, H.R.; Lee, K.J.; Hong, M.S.; Jang, I.; et al. Novel reassortant influenza A(H5N8) viruses, South Korea, 2014. Emerg. Infect. Dis. 2014, 20, 1087-1089.

15. Swayne, D.E.; Suarez, D.L.; Spackman, E.; Jadhao, S.; Dauphin, G.; Kim-Torchetti, M.; McGrane, J.; Weaver, J.; Daniels, P.; Wong, F.; et al. Antibody titer has positive predictive value for vaccine protection against challenge with natural antigenic-drift variants of H5N1 high-pathogenicity avian influenza viruses from Indonesia. J. Virol. 2015, 89, 3746-3762.

16. Wu, H.; Peng, X.; Xu, L.; Jin, C.; Cheng, L.; Lu, X.; Xie, T.; Yao, H.; Wu, N. Novel reassortant influenza A(H5N8) viruses in domestic ducks, eastern China. Emerg. Infect. Dis. 2014, 20, 1315-1318.

17. Zhao, G.; Gu, X.; Lu, X.; Pan, J.; Duan, Z.; Zhao, K.; Gu, M.; Liu, Q.; He, L.; Chen, J.; et al. Novel reassortant highly pathogenic H5N2 avian influenza viruses in poultry in China. PLoS ONE 2012, 7, e46183.

18. WHO. Antigenic and Genetic Characteristics of Zoonotic Influenza A Viruses and Development of Candidate Vaccine Viruses for Pandemic Preparedness. Available online: https://www.who.int/influenza/vaccines/virus/201909_zoonotic_vaccinevirusupdate.pdf?ua=1 (accessed on 14 December 2020).

19. Yamamoto, N.; Sakoda, Y.; Motoshima, M.; Yoshino, F.; Soda, K.; Okamatsu, M.; Kida, H. Characterization of a non-pathogenic H5N1 influenza virus isolated from a migratory duck flying from Siberia in Hokkaido, Japan, in October 2009. Virol. J. 2011, 8, 65.

20. Kida, H.; Yanagawa, R. Isolation and characerization of influenza a viruses from wild free-flying ducks in Hokkaido Japan. Zentralbl. Bakteriol. Orig. 1979, 1979, 135-143.

21. Hebert, P.D.; Stoeckle, M.Y.; Zemlak, T.S.; Francis, C.M. Identification of Birds through DNA Barcodes. PLoS Biol. 2004, 2, e312.

22. Hoffmann, E.; Stech, J.; Guan, Y.; Webster, R.G.; Perez, D.R. Universal primer set for the full-length amplification of all influenza A viruses. Arch. Virol. 2001, 146, 2275-2289.

23. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870-1874.

24. Smith, D.J.; Lapedes, A.S.; de Jong, J.C.; Bestebroer, T.M.; Rimmelzwaan, G.F.; Osterhaus, A.D.; Fouchier, R.A. Mapping the antigenic and genetic evolution of influenza virus. Science 2004, 305, 371-376.

25. Offlu. Influenza A Cleavage Sites. Available online: http://www.offlu.net/fileadmin/home/en/resource-centre/pdf/Influenza_A_Cleavage_Sites.pdf (accessed on 12 November 2020).

26. World Organisation for Animal Health. Weekly Disease Information in World Animal Health Information Database Interface. Available online: https://www.oie.int/wahis_2/public/wahid.php/Diseaseinformation/WI (accessed on 12 November 2020).

27. Ito, T.; Okazaki, K.; Kawaoka, Y.; Takada, A.; Webster, R.G.; Kida, H. Perpetuation of influenza A viruses in Alaskan waterfowl reservoirs. Arch. Virol. 1995, 140, 1163-1172.

28. Okazaki, K.; Takada, A.; Ito, T.; Imai, M.; Takakuwa, H.; Hatta, M.; Ozaki, H.; Tanizaki, T.; Nagano, T.; Ninomiya, A.; et al. Precursor genes of future pandemic influenza viruses are perpetuated in ducks nesting in Siberia. Arch. Virol. 2000, 145, 885-893.

29. Mine, J.; Uchida, Y.; Nakayama, M.; Tanikawa, T.; Tsunekuni, R.; Sharshov, K.; Takemae, N.; Sobolev, I.; Shestpalov, A.; Saito, T. Genetics and pathogenicity of H5N6 highly pathogenic avian influenza viruses isolated from wild birds and a chicken in Japan during winter 2017-2018. Virology 2019, 533, 1-11.

30. Okamatsu, M.; Ozawa, M.; Soda, K.; Takakuwa, H.; Haga, A.; Hiono, T.; Matsuu, A.; Uchida, Y.; Iwata, R.; Matsuno, K.; et al. Characterization of Highly Pathogenic Avian Influenza Virus A(H5N6), Japan, November 2016. Emerg. Infect. Dis. 2017, 23, 691-695.

31. Tanikawa, T.; Kanehira, K.; Tsunekuni, R.; Uchida, Y.; Takemae, N.; Saito, T. Pathogenicity of H5N8 highly pathogenic avian influenza viruses isolated from a wild bird fecal specimen and a chicken in Japan in 2014. Microbiol. Immunol. 2016, 60, 243-252.

32. Chen, W.; Doko, T.; Fujita, G.; Hijikata, N.; Tokita, K.; Uchida, K.; Konishi, K.; Hiraoka, E.; Higuchi, H. Migration of Tundra Swans (Cygnus columbianus) Wintering in Japan Using Satellite Tracking: Identification of the Eastern Palearctic Flyway. Zoolog. Sci. 2016, 33, 63-72.

33. Yamaguchi, N.; Hiraoka, E.; Fujita, M.; Hijikata, N.; Ueta, M.; Takagi, K.; Konno, S.; Okuyama, M.; Watanabe, Y.; Osa, Y.; et al. Spring migration routes of mallards (Anas platyrhynchos) that winter in Japan, determined from satellite telemetry. Zoolog. Sci. 2008, 25, 875-881.

34. Twabela, A.T.; Okamatsu, M.; Tshilenge, G.M.; Mpiana, S.; Masumu, J.; Nguyen, L.T.; Matsuno, K.; Monne, I.; Zecchin, B.; Sakoda, Y. Molecular, antigenic, and pathogenic characterization of H5N8 highly pathogenic avian influenza viruses isolated in the Democratic Republic of Congo in 2017. Arch. Virol. 2020, 165, 87-96.

35. MAFF. Report of High Pathogenicity Avian Influenza. Available online: https://www.maff.go.jp/j/syouan/douei/tori/r2_hpai_kokunai.html (accessed on 10 December 2020).