- 1College of Veterinary Medicine, Kyungpook National University, Daegu, Republic of Korea
- 2Avian Disease Laboratory, College of Veterinary Medicine, Konkuk University, Seoul, Republic of Korea
- 3Wildlife Health Laboratory, College of Veterinary Medicine, Konkuk University, Seoul, Republic of Korea
Avian influenza viruses (AIVs) in wild birds are phylogenetically separated in Eurasian and North American lineages due to the separated distribution and migration of wild birds. However, AIVs are occasionally dispersed between two continents by migratory wild birds flying across the Bering Strait. In this study, we isolated three AIVs from wild bird feces collected in South Korea that contain gene segments derived from American lineage AIVs, including an H6N2 isolated in 2015 and two H6N1 in 2017. Phylogenetic analysis suggests that the H6N2 virus had American lineage matrix gene and the H6N1 viruses had American lineage nucleoprotein and non-structural genes. These results highlight that novel AIVs have continuously emerged by reassortment between viruses from the two continents. Therefore, continuous monitoring for the emergence and intercontinental spread of novel reassortant AIV is required to prepare for a possible future outbreak.
1. Introduction
Wild waterfowl are the natural hosts of avian influenza viruses (AIVs) (1) Because of the geographical barrier, AIVs are separated into two phylogenetic lineages, the Eurasian and American lineages (2). However, some migratory waterfowl (e.g., the Northern Pintail duck; Anas acuta and Greater White-fronted goose; Anser albifrons) move intercontinentally, causing genetic mixing between the two lineages of AIVs (3–8).
Previous surveillance studies on AIVs in wild birds provide evidence for the intercontinental exchange of AIVs. For example, the Eurasian H6 subtype AIV has been reported in North America since the 1990s and it replaced the prevailing North American H6 AIV (8). AIVs subtype H9N2, which contains six genes originating from North America, were simultaneously isolated in South Korea, China, and Alaska (5–7). In 2019, subtype H6N5 AIV carrying all eight gene segments from North American ancestors was detected in Mandarin duck in South Korea (4). In addition, the intercontinental spread of the highly pathogenic avian influenza (HPAI) virus from Eurasia to North America was detected in 2014 and 2021 (9, 10). Western Alaska, where migratory flyways of waterfowl overlap, is the location to encounter AIVs from Eurasian and North American lineages (3).
In the current study, we report three reassortant H6 viruses containing gene segments originating from North America. Complete genome sequencing and phylogenetic analysis were used to find the origin of each gene segment.
2. Materials and methods
We collected fresh fecal samples from wild bird habitats in South Korea for routine AIV surveillance during the fall migration and wintering seasons of wild waterfowl. In the 2015–2016 winter, we collected total of 1,896 fecal samples and in the 2016–2017 winter, we collected total of 8,096 fecal samples. Following as the previous study, fecal samples were examined for influenza A virus by real-time reverse transcription polymerase chain reaction (rRT-PCR) targeting the matrix (M) gene (11). Virus isolation was done using embryonated SPF chicken eggs. We found three reassortant H6 viruses containing gene segments originating from the American continent: A/Mandarin duck/Korea/K15-68/2015 (H6N2), A/Greater White-fronted goose/Korea/K16-727-5/2017 (H6N1), and A/Greater White-fronted goose/Korea/K16-738/2017 (H6N1) (designated as K15-68, K16-727-5, and K16-738, respectively). The date and location information of the isolated viruses is shown in Table 1. We sequenced full-length genomes of the isolates using the Illumina MiSeq system. We deposited the nucleotide sequences of each virus into the Global Initiative for Sharing All Influenza Data (GISAID) database (accession nos. EPI_ISL_11110143, EPI_ISL_11112483, and EPI_ISL_11112543, respectively). The host of each fecal sample positive for AIV was identified using a DNA barcoding technique, as previously described (12).
Comparative phylogenetic analysis of each gene was conducted to trace their origin. For phylogenetic analysis, all available sequences collected between 2010 and 2021 were downloaded from GISAID. To prevent the omission of intercontinental spread viruses during the subsampling process, we classified all downloaded sequences into two groups before subsampling: one containing the viruses isolated from Asia, Africa, and Europe and the other containing the viruses isolated from South and North America. We reduced the number of sequences in each gene segment of each group based on sequence identities of 97 ~ 99% using the program CD-HIT (13). Maximum-likelihood (ML) phylogenetic trees were constructed using a general time reversible (GTR) substitution model with 500 bootstrap replications in RAxML version 8.2.11 (14). We used BLAST1 to search for sequences with the highest identity to each virus for each gene segment.
To verify the result of the phylogenetic tree, we constructed time-scaled phylogenetic trees using BEAST version 1.10.4 (15). The GTR nucleotide substitution model and uncorrelated lognormal relaxed molecular clock model were used for constructing a time-scaled maximum clade credibility (MCC) tree. MCC trees were visualized using FigTree 1.4.2.2
To evaluate the pathogenicity of the viruses in chickens, a total of 15 three-week-old SPF chickens (Namdeok SPF, Korea) were divided into 3 groups (5 chickens/group). The chickens were inoculated with H6 LPAI viruses at 106.0 EID50 in a volume of 100ul by intranasal route. At 2-, 3-, 5-, and 7-days post-infection (dpi), oropharyngeal (OP) and cloacal (CL) swabs were collected from all chickens and examined for virus shedding using a quantitative real-time reverse transcription polymerase chain reaction (rRT-PCR) targeting the matrix gene, as described previously (11). For each virus, the standard curve was used to convert the Ct values into equivalent EID50 titer. All chickens were monitored daily for clinical signs and mortality. Serum samples were collected for serological investigations including anti-NP ELISA (Bionote, Inc., Korea) and Hemagglutination inhibition (HI) testing for homologous HA-specific antibodies.
3. Results
We successfully obtained complete genome sequences of the K15-68, K16-727-5, and K16-738 viruses. Host species of K15-68 was identified as Mandarin duck, and K16-727-5 and K16-738 viruses’ host were identified as Greater White-fronted goose using DNA barcoding. The genome sequences of K16-727-5 and K16-738 which were isolated on the same sample collection day were almost identical (NP: 99.936%, NS: 99.944%, and other 6 genes: 100%).
BLAST research indicated that the matrix (M) gene of K15-68 strain shared >99% nucleotide identity with the Guatemalan origin H14N3 subtype AIV. The nucleoprotein (NP) gene of K16-727-5 and K16-738 strains sharing 99.94% nucleotide identity with H9N6 subtype AIV isolated from Missouri, United States, and the nonstructural (NS) gene shared 98.31% nucleotide identity with Ohio, USA isolate (Table 2). All other segments showed >98% identity with the low-pathogenicity AIVs (LPAI) identified in South Korea, Japan, China, and Mongolia. Corresponding to the BLAST results, the ML and Bayesian phylogenetic analysis showed that the M gene of K15-68 and the NP and NS gene of K16-727-5 and K16-738 clustered with the American lineage wild bird AIVs and the rest of the gene segments were clustered with the Eurasian lineage wild bird AIVs (Figure 1; Supplementary Figures S1, S2). These results demonstrate that the three viruses were generated by reassortment between Eurasian and American AIVs.
Figure 1. Maximum-likelihood phylogenetic trees of each gene segment of avian influenza virus including Eurasian lineage and American lineage genome sequences. Red and blue circles indicate lineages of Eurasia and America, respectively. Green circles indicate genome sequences of novel viruses which isolated in this study; K15-68, K16-727-5, K16-738. Maximum-likelihood phylogenetic trees for each gene segment are shown in Supplementary Figure S1.
The HA genes from H6N1 and H6N2 AIVs have some phylogenetic distance. Eurasian lineage H6 genes are divided into seven groups; ST339-like, W312-like, ST2853-like, HN573-like, South African-like, Taiwan-like, and Europe (Figure 1) (16). HN573-like is divided into three subgroups, subgroup 1 (mixed group), 2 (Eurasian group), and 3 (American group). All three isolates cluster with HN573-like, but with different subgroups (Supplementary Figure S1D). H6N2 belongs to subgroup 1 and H6N1 belongs to subgroup 2.
The amino acid changes related to the mammalian adaptation, including Q591K, E627K, and D701N mutations in PB2 (17), N137/ E190V/ G228S triad, and Q226L mutation in HA (H3 numbering) (18, 19), were not detected in the viruses isolated in this study. The NA stalk deletion, which has been associated with adaptation to gallinaceous hosts (20), was also not detected.
Three-week-old SPF chickens were inoculated with the H6 viruses to study their infectivity and virulence. During the 14-day experiment, no obvious clinical signs and no mortality were observed. Virus shedding through the OP and CL routes was detected during 3–7 dpi, but the detected amount of viruses was very low (<103.5 EID50) (Table 2). For the chickens inoculated with K15-68 and K16-727-5 virus, mean shedding titers for both routes did not exceed 101.0EID50 throughout the whole period of the experiment. Chickens inoculated with K16-738 showed a peak of OP shedding at 3 dpi and a peak of CL shedding at 7 dpi, with mean shedding titers of 101.85 EID50 and 101.84 EID50, respectively. Both HI and NP-specific antibodies were detected in two of five, one of five, and one of five chickens inoculated with K15-68, K16-727-5, and K16-768, respectively (Table 3).
Table 3. Viral shedding and antibody responses of chickens inoculated with H6 avian influenza viruses.
4. Discussion
In the current study, we report American-Eurasian reassortants H6N1 and H6N2 viruses isolated from two wild bird species, the Greater White-fronted goose and Mandarin duck. Since the first report of North American-Eurasian reassortant AIVs detected from wild birds in South Korea in 2010 (21), North American lineage AIV gene segments have been detected continuously in South Korea (4, 21–23). Previous studies also detected multiple Eurasian lineage AIV genes in Alaska, indicating the bi-directional flow of AIV exchange between East Asia and North America (24).
Recently, intercontinental reassortants AIVs have been frequently detected (4, 6, 9, 21–23, 25, 26). The reason for the frequent discovery of intercontinental reassortants is not fully determined yet. We assume that climate change can be a factor in contributing to frequent exchange of AIVs between two continents. A previous study indicated that an abnormal climate in Africa might have contributed to the trans-continental introduction of HPAI H5Nx in Africa by migratory birds (27). The temperature of Alaska has increased since 2014 and it might affect the migration pattern of migratory birds and the ecology of AIVs in wild birds (28). On the other hand, Next-Generation Sequencing techniques have been widely used for AIV sequencing (23), and this high-throughput sequencing system would also contribute to the frequent detections of novel reassortants.
The reassortment of AIVs causes rapid changes in the biological characteristics of viruses and it could be a considerable threat to the poultry industry and public health (29). In 2014, in China, for example, human infection of AIV subtype H7N9 occurred by reassortment of AIVs in poultry and wild birds (30). In 2013, a case of human infection with subtype H6N1, AIV generated by reassortment carrying N137/E190V/G228S triad amino acid changes in HA, was reported in Taiwan (31). Although no mammalian adaptation mutation was detected in this study, the emergence of novel AIVs reassortants and the detection of mutation related to mammalian adaptation highlights the importance of wild bird surveillance in terms of the ‘One Health’ concept.
The viruses isolated in this study showed limited infectivity in chickens similar to other wild bird-origin H6 viruses tested in previous studies (32–35). These H6 viruses may be poorly adapted to chickens, but we assume that they possibly replicate in domestic waterfowl without prior adaptation which could lead to the spread and maintenance of H6 viruses in land-based poultry. Bahl et al. showed that the introduction and establishment of Eurasian H6 viruses in North America dramatically changed the evolutionary dynamics of the influenza virus in wild birds (36). Thus, genomic surveillance and in vivo pathobiology studies should be continued to monitor the evolution of reassortant H6 viruses and their host ranges.
Constant monitoring for AIVs in wild birds is essential to detect the introduction of new viruses and trace the dispersion path. Particularly, as a major wintering site for various wild birds (37), South Korea would be an important location for monitoring the intercontinental exchange of AIVs. Combined with recently developed sequencing technology, continuing monitoring of the intercontinental dispersion of AIVs would expand our knowledge of AIV ecology and epidemiology.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Konkuk University.
Author contributions
C-SS, D-HL, and J-HK contributed to conception and design of the study. S-HL collected and analyzed samples and conducted an animal experiment. J-YK, D-WK, and D-WL organized data and performed genetic analysis. J-YK visualized data and wrote the first draft of the manuscript. S-HL wrote animal experiment sections. D-HL and J-HK revised final manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by Kyungpook National University Research Fund, 2020.
Acknowledgments
The authors gratefully acknowledge the authors and submitting laboratories of the sequences from GISAID’s EpiFlu (https://www.gisaid.org/) database. The GISAID acknowledgement tables for laboratory contributions are shown in Supplementary Table S1.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1157984/full#supplementary-material
Footnotes
References
2. Olsen, B, Munster, VJ, Wallensten, A, Waldenström, J, Osterhaus, AD, and Fouchier, RA. Global patterns of influenza a virus in wild birds. Science. (2006) 312:384–8. doi: 10.1126/science.1122438
3. Pearce, JM, Reeves, AB, Ramey, AM, Hupp, JW, Ip, HS, Bertram, M, et al. Interspecific exchange of avian influenza virus genes in Alaska: the influence of trans-hemispheric migratory tendency and breeding ground sympatry. Mol Ecol. (2011) 20:1015–25. doi: 10.1111/j.1365-294X.2010.04908.x
4. Jeong, S, Lee, D-H, Kim, Y-J, Lee, S-H, Cho, AY, Noh, J-Y, et al. Introduction of avian influenza a (H6N5) virus into Asia from North America by wild birds. Emerg Infect Dis. (2019) 25:2138–40. doi: 10.3201/eid2511.190604
5. Lee, D-H, Park, J-K, Yuk, S-S, Erdene-Ochir, T-O, Kwon, J-H, Lee, J-B, et al. Complete genome sequence of a natural reassortant H9N2 avian influenza virus found in bean goose (Anser fabalis): direct evidence for virus exchange between Korea and China via wild birds. Infect Genet Evol. (2014) 26:250–4. doi: 10.1016/j.meegid.2014.06.007
6. Ramey, AM, Reeves, AB, Sonsthagen, SA, TeSlaa, JL, Nashold, S, Donnelly, T, et al. Dispersal of H9N2 influenza a viruses between East Asia and North America by wild birds. Virology. (2015) 482:79–83. doi: 10.1016/j.virol.2015.03.028
7. Wang, B, Chen, Q, and Chen, Z. Complete genome sequence of an H9N2 avian influenza virus isolated from egret in Lake Dongting wetland. J Virol. (2012) 86:11939. doi: 10.1128/JVI.02042-12
8. Zu Dohna, H, Li, J, Cardona, CJ, Miller, J, and Carpenter, TE. Invasions by Eurasian avian influenza virus H6 genes and replacement of its north American clade. Emerg Infect Dis. (2009) 15:1040–5. doi: 10.3201/eid1507.090245
9. Ip, HS, Torchetti, MK, Crespo, R, Kohrs, P, DeBruyn, P, Mansfield, KG, et al. Novel Eurasian highly pathogenic avian influenza a H5 viruses in wild birds, Washington, USA, 2014. Emerg Infect Dis. (2015) 21:886–90. doi: 10.3201/eid2105.142020
10. Bevins, SN, Shriner, SA, Cumbee, JC Jr, Dilione, KE, Douglass, KE, Ellis, JW, et al. Intercontinental movement of highly pathogenic avian influenza a (H5N1) clade 2.3. 4.4 virus to the United States, 2021. Emerg Infect Dis. (2022) 28:1006–11. doi: 10.3201/eid2805.220318
11. Spackman, E. Avian influenza virus detection and quantitation by real-time RT-PCR. Methods Mol Biol. (2020) 2123:137–48. doi: 10.1007/978-1-0716-0346-8_11
12. Hebert, PDN, Stoeckle, MY, Zemlak, TS, Francis, CM, and Godfray, C. Identification of birds through DNA barcodes. PLoS Biol. (2004) 2:e312. doi: 10.1371/journal.pbio.0020312
13. Li, W, and Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. (2006) 22:1658–9. doi: 10.1093/bioinformatics/btl158
14. Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. (2014) 30:1312–3. doi: 10.1093/bioinformatics/btu033
15. Suchard, MA, Lemey, P, Baele, G, Ayres, DL, Drummond, AJ, and Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. (2018) 4:vey016. doi: 10.1093/ve/vey016
16. Everest, H, Hill, SC, Daines, R, Sealy, JE, James, J, Hansen, R, et al. The evolution, spread and global threat of H6Nx avian influenza viruses. Viruses. (2020) 12:673. doi: 10.3390/v12060673
17. Yamada, S, Hatta, M, Staker, BL, Watanabe, S, Imai, M, Shinya, K, et al. Biological and structural characterization of a host-adapting amino acid in influenza virus. PLoS Pathog. (2010) 6:e1001034. doi: 10.1371/journal.ppat.1001034
18. Maines, TR, Chen, L-M, Van Hoeven, N, Tumpey, TM, Blixt, O, Belser, JA, et al. Effect of receptor binding domain mutations on receptor binding and transmissibility of avian influenza H5N1 viruses. Virology. (2011) 413:139–47. doi: 10.1016/j.virol.2011.02.015
19. Liu, WJ, Wu, Y, Bi, Y, Shi, W, Wang, D, Shi, Y, et al. Emerging HxNy influenza a viruses. Cold Spring Harb Perspect Med. (2022) 12:a038406. doi: 10.1101/cshperspect.a038406
20. Li, J, Dohna, HZ, Cardona, CJ, Miller, J, and Carpenter, TE. Emergence and genetic variation of neuraminidase stalk deletions in avian influenza viruses. PLoS One. (2011) 6:e14722. doi: 10.1371/journal.pone.0014722
21. Lee, D-H, Lee, H-J, Lee, Y-N, Park, J-K, Lim, T-H, Kim, M-S, et al. Evidence of intercontinental transfer of north American lineage avian influenza virus into Korea. Infect Genet Evol. (2011) 11:232–6. doi: 10.1016/j.meegid.2010.09.012
22. Nguyen, NM, Sung, HW, Yun, K-J, Park, H, and Yeo, S-J. Genetic characterization of a novel north American-origin avian influenza a (H6N5) virus isolated from bean goose of South Korea in 2018. Viruses. (2020) 12:774. doi: 10.3390/v12070774
23. Lee, D-H, Park, J-K, Yuk, S-S, Erdene-Ochir, T-O, Kwon, J-H, Lee, J-B, et al. Complete genome sequence of a natural recombinant H9N2 influenza virus isolated from a white-fronted goose (Anser albifrons) in South Korea. Genome Announc. (2013) 1:e00149–13. doi: 10.1128/genomeA.00149-13
24. Reeves, AB, Pearce, JM, Ramey, AM, Ely, CR, Schmutz, JA, Flint, PL, et al. Genomic analysis of avian influenza viruses from waterfowl in western Alaska, USA. J Wildl Dis. (2013) 49:600–10. doi: 10.7589/2012-04-108
25. Lee, D-H, Torchetti, MK, Winker, K, Ip, HS, Song, C-S, and Swayne, DE. Intercontinental spread of Asian-origin H5N8 to North America through Beringia by migratory birds. J Virol. (2015) 89:6521–4. doi: 10.1128/JVI.00728-15
26. González-Reiche, AS, Morales-Betoulle, ME, Alvarez, D, Betoulle, J-L, Müller, ML, Sosa, SM, et al. Influenza a viruses from wild birds in Guatemala belong to the north American lineage. PLoS One. (2012) 7:e32873. doi: 10.1371/journal.pone.0032873
27. Fusaro, A, Zecchin, B, Vrancken, B, Abolnik, C, Ademun, R, Alassane, A, et al. Disentangling the role of Africa in the global spread of H5 highly pathogenic avian influenza. Nat Commun. (2019) 10:1–13. doi: 10.1038/s41467-019-13287-y
28. Coggin, JDP. New Report Highlights Alaska's Last Five Years of Dramatic Climate Change 2019. Available at: https://www.climate.gov/news-features/understanding-climate/new-report-highlights-alaska%E2%80%99s-last-five-years-dramatic-climate
29. Su, S, Bi, Y, Wong, G, Gray, GC, Gao, GF, and Li, S. Epidemiology, evolution, and recent outbreaks of avian influenza virus in China. J Virol. (2015) 89:8671–6. doi: 10.1128/JVI.01034-15
30. Cui, L, Liu, D, Shi, W, Pan, J, Qi, X, Li, X, et al. Dynamic reassortments and genetic heterogeneity of the human-infecting influenza a (H7N9) virus. Nat Commun. (2014) 5:1–9. doi: 10.1038/ncomms4142
31. Wei, S-H, Yang, J-R, Wu, H-S, Chang, M-C, Lin, J-S, Lin, C-Y, et al. Human infection with avian influenza a H6N1 virus: an epidemiological analysis. Lancet Respir Med. (2013) 1:771–8. doi: 10.1016/S2213-2600(13)70221-2
32. Jackwood, MW, Suarez, DL, Hilt, D, Pantin-Jackwood, MJ, Spackman, E, Woolcock, P, et al. Biologic characterization of chicken-derived H6N2 low pathogenic avian influenza viruses in chickens and ducks. Avian Dis. (2010) 54:120–5. doi: 10.1637/8987-070909-ResNote.1
33. Lin, W, Cui, H, Teng, Q, Li, L, Shi, Y, Li, X, et al. Evolution and pathogenicity of H6 avian influenza viruses isolated from southern China during 2011 to 2017 in mice and chickens. Sci Rep. (2020) 10:1–12. doi: 10.1038/s41598-020-76541-0
34. Nam, J-H, Kim, E-H, Song, D, Choi, YK, Kim, J-K, and Poo, H. Emergence of mammalian species-infectious and-pathogenic avian influenza H6N5 virus with no evidence of adaptation. J Virol. (2011) 85:13271–7. doi: 10.1128/JVI.05038-11
35. Yao, Y, Wang, H, Chen, Q, Zhang, H, Zhang, T, Chen, J, et al. Characterization of low-pathogenic H6N6 avian influenza viruses in Central China. Arch Virol. (2013) 158:367–77. doi: 10.1007/s00705-012-1496-3
36. Bahl, J, Vijaykrishna, D, Holmes, EC, Smith, GJ, and Guan, Y. Gene flow and competitive exclusion of avian influenza a virus in natural reservoir hosts. Virology. (2009) 390:289–97. doi: 10.1016/j.virol.2009.05.002
Keywords: influenza virus, reassortment, phylogenetic analysis, wild bird, South Korea
Citation: Kim J-Y, Lee S-H, Kim D-W, Lee D-W, Song C-S, Lee D-H and Kwon J-H (2023) Detection of intercontinental reassortant H6 avian influenza viruses from wild birds in South Korea, 2015 and 2017. Front. Vet. Sci. 10:1157984. doi: 10.3389/fvets.2023.1157984
Edited by:
Sarah J. Edwards, Commonwealth Scientific and Industrial Research Organisation (CSIRO), AustraliaCopyright © 2023 Kim, Lee, Kim, Lee, Song, Lee and Kwon. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dong-Hun Lee, ZG9uZ2h1bmxlZUBrb25rdWsuYWMua3I=; Jung-Hoon Kwon, SnVuZ2hvb24ua3dvbkBrbnUuYWMua3I=