Next Article in Journal
Decrypting the Origin and Pathogenesis in Pregnant Ewes of a New Ovine Pestivirus Closely Related to Classical Swine Fever Virus
Next Article in Special Issue
Putative Role of Arthropod Vectors in African Swine Fever Virus Transmission in Relation to Their Bio-Ecological Properties
Previous Article in Journal
Aphid Transmission of Potyvirus: The Largest Plant-Infecting RNA Virus Genus
Previous Article in Special Issue
Genotypic Characterization of Epstein Barr Virus in Blood of Patients with Suspected Nasopharyngeal Carcinoma in Ghana
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 7
10.3390/v12070774
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Characterization of a Novel North American-Origin Avian Influenza A (H6N5) Virus Isolated from Bean Goose of South Korea in 2018

by
Ngoc Minh Nguyen
1,
Haan Woo Sung
2,
Ki-Jung Yun
3,
Hyun Park
1,* and
Seon-Ju Yeo
1,*
1
Zoonosis Research Center, Department of Infection Biology, School of Medicine, Wonkwang University, Iksan 570-749, Korea
2
College of Veterinary Medicine, Kangwon National University, Chuncheon 200-701, Korea
3
Department of Pathology, School of Medicine, Wonkwang University, Iksan 570-749, Korea
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(7), 774; https://0-doi-org.brum.beds.ac.uk/10.3390/v12070774
Submission received: 7 June 2020 / Revised: 13 July 2020 / Accepted: 15 July 2020 / Published: 17 July 2020
(This article belongs to the Special Issue Emerging Viruses 2020: Surveillance, Prevention, Evolution and Control)
Browse Figures
Versions Notes

Abstract

:
The complex overlap in waterfowl migratory pathways across the world has established numerous occurrences of genetic reassortment and intercontinental spread of avian influenza virus (AIV) over long distances, thereby calling for huge efforts and targeted surveillance for infection control. During annual surveillance in South Korea in 2018, a novel avian influenza H6N5 (K6) subtype was isolated from the fecal sample of wild bird. Genomic characterization using a phylogenetic tree indicated the K6 virus to be of North American-origin, with partial homology to an H6N5 strain, A/Aix galericulata/South Korea/K17-1638-5/2017 (K17). A monobasic residue at the HA cleavage site and absence of a notable mutation at the HA receptor-binding site suggested the isolate to be of low pathogenicity. However, molecular analysis revealed the E119V mutation in the NA gene and a human host marker mutation E382D in the polymerase acidic (PA) gene, implying their susceptibility to neuraminidase inhibitors and potential infectivity in humans, respectively. For comparison, K6 and K17 were found to be dissimilar for various mutations, such as A274T of PB2, S375N/T of PB1, or V105M of NP, each concerning the increased virulence of K6 in mammalian system. Moreover, kinetic data presented the highest viral titer of this H6N5 isolate at 106.37 log10TCID50 after 48 h of infection, thus proving efficient adaptability for replication in a mammalian system in vitro. The mouse virus challenge study showed insignificant influence on the total body weight, while viral load shedding in lungs peaked at 1.88 ± 0.21 log10 TICD50/mL, six days post infection. The intercontinental transmission of viruses from North America may continuously be present in Korea, thereby providing constant opportunities for virus reassortment with local resident AIVs; these results hint at the increased potential risk of host jumping capabilities of the new isolates. Our findings reinforce the demand for regular surveillance, not only in Korea but also along the flyways in Alaska.

1. Introduction

Based on surface glycoproteins, avian influenza virus (AIV) has been divided into subtypes including 18 HA and 11 NA [1], which are further distributed based on their geographical separation, into the Eurasian and American lineages [2]. However, despite the physical barriers, intercontinental circulation and reassortment of these ancestries have been frequently observed.
The main reason underlying viral dispersal over long distances is likely the overlap of flyways of wild migratory birds, which happen to be the major hosts for AIV. Particularly, four flyways (East Asian–Australasian, West Pacific, Pacific Americas, and Mississippi Americas) intersect at the boundary of the Eastern and Western hemispheres, near Alaska, causing an extensive diversity of Eurasian and North American bird species to be found in this area during spring and fall migrations [3]. Certain key vector bird species, such as green-winged teal (Anas crecca), mallard (Anas platyrhynchos), and northern pintails (Anas acuta), have been reported to carry viruses containing genetic divergence of both Asian and North American species, showing a high prevalence of group HA-1, 3, 4, 5, 6 and NA-1, 2, 6, 8 during the 2006–2008 surveillance in Alaska [4,5,6].
South Korea has recently recorded converse intercontinental exchanges from the North American lineage of the nonstructural (NS) gene segment in several subtypes (H1N1, H1N8, H4N6, H6N1, H6N2, and H6N8) during the surveillance from 2012 to 2017, as a consequence of wild bird (greater white-fronted geese) movement during summer and winter, thereby suggesting the introduction of a gene from American-origin AIV in Korea [7].
In a recent publication, an H6N5 strain, A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5) (designated as K17) was isolated from the feces of a mandarin duck during July 2017–March 2018. This virus carried the full gene segment from a North American ancestor [8]. However, mandarin ducks (Aix galericulata) belong to the Asian population that is primarily found in Southeast China, and, therefore, does not regularly participate in the transfer between Eurasia and North America; hence, there would likely be alternative wild birds involved in the transmission of H6N5 to the mandarin ducks in Korea [8]. Frequently, the H6 subtype has been recorded as the cause of widespread outbreaks in poultry and may also sporadically infect humans [9,10,11]. Previously, an unconventional H6N5 strain in Korea had been characterized as a novel pathogenic AIV, capable of infecting mammalian species, with no evidence of adaptation [12]; this suggested the potential risk of cross-species, although the isolate was confirmed to contain all gene segments from Eurasian ancestry.
The impact of foreign genes imported to the local population of AIV H6 in Korea has been rarely described. Therefore, we aimed to emphasize the need for continued surveillance and characterization of H6 outbreaks. In this study, a novel H6N5 strain introduced into Korea from North America was investigated for replicability, pathogenicity, and toxicity in mice. Specifically, the analysis based on satellite tracking devices attached to host species provided valuable information to understand how North American-origin H6N5 could spread in Korea.

2. Materials and Methods

2.1. Sample Collection

Fresh feces of wild birds in South Korea were collected using sterile swabs in March 2018. The collected samples were stored at 2–8 °C, before being shipped to the laboratory within 12 h, for further analysis.

2.2. Virus Isolation from the Feces

As described previously [13], the fecal samples were resuspended in phosphate-buffered saline (pH 7.4; supplemented with 100 U/μL penicillin and 100 mg/μL streptomycin), thoroughly mixed by vortexing, and centrifugated at 3000 rpm for 10 min at 4 °C. The supernatants were reserved for subsequent filtering (0.45-µm pore size; GVS Syringe, Nova-Tech, Kingwood, TX, USA), and then inoculated into the allantoic cavities of specific pathogen-free 10-day-old embryonated chicken eggs (Seng-Jin Inc., Eumsung, Korea) for 72 h at 37 °C, under humidified conditions, eventually chilling them at 4 °C overnight. The presence of influenza virus in allantoic fluids of the inoculated eggs was first verified using an HA assay [14].

2.3. RNA Extraction and Bird Identification

Viral RNA was extracted directly from confirmed-HA-positive samples using the RNeasy Mini kit (Macherey-Nagel, Düren, Germany) according to the manufacturer’s instructions. RNA was eluted in 60 μL of RNase-free water, supplied with 20 units of RNase inhibitor for storage at −80 °C. To verify that the samples contained Influenza A viruses, reverse transcription polymerase chain reaction (RT-PCR) was performed to amplify the matrix gene, following the World Health Organization guideline, using universal primers and probes [15]. Bird identification was conducted via barcoding technique, using primers for identifying the mitochondrial cytochrome oxidase (COI) gene in DNA isolated from the fecal samples, as described previously [16,17].

2.4. Next-Generation Sequencing (NGS) by Illumina HiSeq X Method

NGS was conducted by GnCBio (Daejeon, Korea), following the HiSeq X method, as reported previously [13,18]. Briefly, influenza viral RNA was detected using an RNA 6000 pico kit (Agilent Technologies, Santa Clara, CA, USA), and RNA concentration was measured with a spectrophotometer. A cDNA library of influenza viral RNA was prepared using the QIAseq FX single cell RNA library kit (Qiagen, Hilden, Germany). The library concentration was measured by LightCycler qPCR (Roche, Basel, Switzerland), and library size was checked using TapeStation HS D5000 ScreenTape (Agilent Technologies, Santa Clara, CA, USA. For cluster generation, the library was loaded into a flow cell, in which fragments were captured on a lawn of surface-bound oligos complementary to the library adapters. Each fragment was then amplified into distinct clonal clusters through bridge amplification. When cluster generation was complete, the templates were used for sequencing. Illumina SBS technology utilizes a proprietary reversible terminator-based method that detects single bases as they become incorporated into DNA template strands. Since all four reversible, terminator-bound dNTPs were present during each sequencing cycle, natural competition minimized their incorporation bias and greatly reduced raw error rates, compared to that in other technologies. Results provided highly accurate base-by-base sequencing that virtually eliminated sequence-context-specific errors, even within repetitive sequence regions and homopolymers. Sequencing data was then converted into raw data for analysis.

2.5. Phylogenetic Analysis

Reference sequences were downloaded from both the NCBI-Influenza Virus Resource (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genbank/) and Global Initiative on Sharing All Influenza Data (GISAID; http://www.gisaid.org). Nucleotide homology analysis was compared by BLAST search of the NCBI and/or GISAID databases. The phylogenetic relationships with reference sequences of eight gene segments (polymerase basic 2 (PB2): 2313 nucleotides (nt), polymerase basic 1 (PB1): 2288 nt, polymerase acidic (PA): 2199 nt, HA: 1705 nt, nucleoprotein (NP): 1532 nt, neuraminidase (NA): 1425 nt, matrix (M): 1001 nt, and NS: 864 nt) were determined by alignment using MUSCLE and edited with Molecular Evolutionary Genetics Analysis Version 6.0 (MEGA6) [19]. Rooted phylograms were prepared by the maximum likelihood method with 1000 bootstrap replicates for each gene using MEGA6 software.

2.6. Viral Replication Kinetics in MDCK Cells

To examine the viral replication capacity in mammalian cells, an isolated H6N5 (A/Bean goose/South Korea/KNU18-6/2018(H6N5)) virus, along with reference strains: human-origin H1N1 (A/California/04/2009 (H1N1)), H7N1 (Youngsoori-12-2F4(2012)), and chicken-origin H9N2 (A/chicken/Korea/KNUSWR09/2009 (H9N2)) were infected into Madin-Darby canine kidney (MDCK; obtained from American Type Culture Collection) cells. Briefly, the selected viruses were inoculated into MDCK monolayers at a multiplicity of infection (MOI) of 0.01. The inoculated cells were incubated with DMEM containing 1 µg/mL TPCK-treated trypsin (Sigma-Aldrich, St Louis, MO, USA) at 37 °C. The supernatant was collected at different time points: 12, 24, 36, 48, 60, and 72 h postinfection (hpi) [9]. Viral replication titers were determined by TCID50 assay.

2.7. TCID50—50% Tissue Culture Infectious Dose and EID50—50% Egg Infectious Dose

The TCID50 titers were determined by enzyme-linked immunosorbent assay (ELISA), as reported previously [20]. Briefly, MDCK cells were grown on flat-bottom 96-well plates and washed with phosphate-buffered saline. Thereafter, the cells were inoculated with 10-fold serial dilutions of virus samples. Inoculated cells were incubated at 37 °C and 5% CO2. After 3 days, the medium was removed from microtiter plates, and the cells were washed with PBS, followed by fixation with 80% cold acetone. ELISA was performed and TCID50 titers were calculated via the Reed and Muench method [21]. To measure the EID50 titer, 10-fold serial dilutions of the virus were prepared, and 100 μL of each dilution was inoculated into the chorio-allantoic cavities of 10-day-old specific-pathogen-free embryonated chicken eggs (Seong-Jin Inc., Eumsung, Korea). The eggs were incubated at 37 °C for 3 days. Five eggs were infected with each viral dilution. Harvested allantoic fluid was then tested for hemagglutination (HA) activity and EID50 was determined [22].

2.8. Animal Study

Six-week-old female BALB/c mice were purchased from Orient Bio (Seongnam, Gyeonggi, Korea). Each group of 16 mice was inoculated intranasally (i.n.) with 105 EID50/50 μL of K6-H6N5, H1N1, or H7N1 strains [12]. For the K6-H6N5 group, 105, 106, and 107 EID50/50 μL were inoculated. Anesthesia was conducted with isoflurane, following manufacturer’s instruction (Hana Pharm Co. Ltd., Hwasung, Korea) under the guidelines of vertebrate animal research, the University of Iowa [23]. Three mice per group were sacrificed at 3, 6, and 14 days post infection (dpi), and the virus titer in lung was determined thereafter. Ten mice per group were monitored for weight change, morbidity, and mortality, each day for 14 days post infection. This study was approved by the Animal Ethics Committee of the Wonkwang University (WKU19-64), December 19, 2019 and all methods were performed in accordance with the relevant guidelines and regulations. At 6 dpi, lungs from three mice were placed in 10% formalin/saline for studying pathology. Tissues in 10% formalin/PBS were processed and embedded in paraffin. Lung sections were cut at 4–5 μm from paraffin embedded lung tissues and mounted on glass slides. Histological characterization was analyzed using standard hematoxylin and eosin (H&E) staining and light microscopy (magnification ×100).

2.9. Statistics

The mean, SD, and Student’s t-test were conducted using GraphPad Prism version 5.0. Results were presented as the mean ± SD. A value of p < 0.05 was considered significant. Two-way or one-way analysis of variance (ANOVA) was used to analyze viral replication kinetics in MDCK cells and compare lung weight of viruses infected mouse, respectively.

3. Results

3.1. Genetic Characterization of Novel Avian Influenza A (H6N5) in 2018

In March 2018, fresh fecal samples from wild birds were collected from Gyeonggi-do (37° 45′28.35″N, 126° 40′ 8.60″E), South Korea. Full-length genomes of the isolates were investigated by NGS; consequently, the H6N5 isolate, A/Bean goose/South Korea/KNU18-6/2018(H6N5) (abbreviated as K6), was selected for genomic characterization and deposited in GenBank (Table 1). Bean goose was identified as the bird species carrying the H6N5 virus, based on the COI gene sequence, as shown in Figure S1.
All gene segments of K6, isolated from bean goose in March 2018, shared ≥ 99% nucleotide identity with those of ≥ 1 isolates recovered from wild and domestic birds sampled in North America during 2007–2018 (Table S1). According to BLAST results, five of the eight gene segments (PB2, PA, HA, NA, and NS) for the K6 isolate shared ≥ 99% nucleotide identity with those of K17 (Table 1). In contrast, three internal genes (PB1, NP, and MP) of K6 and K17 shared ≤ 99% nucleotide identity and possessed ≥ 99% nucleotide identity with other strains in North America, implying K6 and K17 to be partially identical (Table 1).
The phylogenetic tree strongly supported the structuring of major clades by continental affiliation of reference strains. Sequence information for all eight gene segments of K6 clustered within the clades composed of reference sequences for AIV originating from samples collected in North America (Figure 1 and Figure S2 and S3).
Based on the phylogenetic analysis, both K6 and K17 belonged to North American strains; specifically, the HA and NA genes may have originated from the same AIV ancestor (A/northern pintail/California/2015 (H6N5)). We, therefore, provided evidence to support the North American ancestry of K6, along with a close relationship with K17. Our results may also provide evidence for the introduction of H6N5, of North American-origin, into Korea via bean geese, based on the previous evidence for the introduction of K17 into Korea by mandarin ducks.

3.2. Molecular Characterization of H6N5 Isolate

Acquisition of multiple basic amino acids at the HA cleavage site and/or additional changes in amino acids in the HA receptor-binding site, or in other gene segments, have been involved in the evolution of an LPAI virus into HPAI [24,25]. Table 2 shows the comparison of several well-known mutations in gene segments between our isolate (K6-2018) and reference strains isolated in Korea: K17-2017 (A/Aix galericulata/South Korea/K17-1638-5/2017 (H6N5)), CN5-2009 (A/aquatic bird/Korea/CN5/2009 (H6N5)), and W69-2005 (A/aquatic bird/Korea/W69/2005(H6N5)), together with an H6 infected-human isolate: TW02/2013: A/Taiwan/2/2013 (H6N1) human infected isolate.
The monobasic residue at the HA cleavage site, with no notable mutation at the HA receptor-binding site, indicated our H6N5 isolate and reference strains to be LPAI (Table 2). Additionally, the NAs of the H6N5 isolates had no amino acid deletion in the stalk region [26]. However, the E119V mutation in NA was detected in all H6N5 isolates, suggesting their susceptibility to neuraminidase inhibitors, such as Oseltamivir, Peramivir, and Zanamivir [27]. Mutations in internal genes, also presented in Table 2, might imply an increase in viral replication efficiency, along with virulence, in mammals. Particularly, one human host marker mutation, E382D, in the PA gene was observed in two new strains imported from North America (K6 and K17), although they were not displayed in two old H6N5 isolates from Korea (CN5 and W69), which had originated from Eurasia. K6 and K17 were variable in certain molecular aspects, such as mutation at A274T of PB2, S375N/T of PB1, or V105M of NP, each contributing to the potentially increased virulence of K6 in mammals (Table 2).
Table
Table 2. Molecular characteristics of influenza A (H6N5) isolates used in this study.
Table 2. Molecular characteristics of influenza A (H6N5) isolates used in this study.
Viral ProteinAmino AcidK6-2018K17-2017CN5-2009W69-2005TW02/2013CommentsReference
HA *Cleavage sitesPQIETR↓GLFPQIETR↓GLFPQIETR↓GLFPQIETR↓GLFPQIATR↓GLFMonobasic-LPAI[24]
HA *-Receptor-binding siteA138SAAAAASwitch the preference from avian- to human-type receptors[25]
E190DEEEEVSwitch the preference from avian- to human-type receptors[25]
G225DGGGGGSwitch the preference from avian- to human-type receptors[25]
Q226LQQQQQSwitch the preference from avian- to human-type receptors[25]
G228SGGGGSSwitch the preference from avian- to human-type receptors[25]
HA*L102FLLLLLIncreased pathogenicity in mice[28]
T/E173G/D/VTTNNTIncreased virus binding to α-2,6-linked sialic acid[28]
N379TNNNNNIncreased virulence in mammals[29]
NAAmino acid deletion at stalk regionNoNoNoNoYesIncreased virulence in mice[26]
E119VVVVVIReduced inhibition by OS (Oseltamivir) and PER (Peramivir) and highly reduced inhibition by ZA (Zanamivir)[21]
PB2T63I (with PB1 M677T)IIIIIPathogenic in mice[30]
L89VVVVVVEnhanced polymerase activity, increased virulence in mice[31]
K251RRRRRKIncreased virulence in mice[32]
A274TTAAAAIncreased polymerase activity, increased virulence in mammals and birds[33]
G309DDDDDDEnhanced polymerase activity, increased virulence in mice[31]
T339KKKKKKEnhanced polymerase activity, increased virulence in mice[31]
Q368RRRRRRIncreased polymerase activity, increased virulence in mammals[34,35]
H447QQQQQQIncreased polymerase activity, increased virulence in mammals[34,35]
R477GGGGGGEnhanced polymerase activity, increased virulence in mice[31]
I495VVVVVVEnhanced polymerase activity, increased virulence in mice[31]
A676TTTTTVEnhanced polymerase active, increased virulence in mice[31]
D701NDDNDDIncreased polymerase activity, increased virulence in mammals, mammalian host marker[36,37,38]
PB1D/A3VVVVVVIncreased polymerase activity, increased virulence in mammals[34,35]
L13PPPPPPIncreased polymerase activity, increased virulence in mammals, mammalian host marker[36,37]
R207KKKKKKIncreased polymerase activity in mammalian cells[39]
K328NNNNNSIncreased polymerase activity, increased virulence in mammals[34,35]
S375N/TNSNNNIncreased polymerase activity, increased virulence in mammals, human host marker[34,35,40]
H436YYYYYYIncreased polymerase activity and virulence in mallards, ferrets, and mice[39]
A469T (with NS1 N205K; NEP T48N)TTTTTConferred in contact transmissibility in guinea pigs[41]
L473VVVVVVIncreased polymerase activity and replication efficiency[42]
V652AAAAAAIncreased virulence in mice[32]
M677T (with PB2 T63I)TTTTTPathogenic in mice[30]
PB1-F2N66SNSNNTruncatedIncreased virulence in mammals[43,44]
PAS37AAAAAASignificantly increased viral growth and polymerase activity in mammalian cells[45]
H266RRRRRRIncreased polymerase activity, increased virulence in mammals and birds[33]
F277SSSSSSAdapt to mammalian hosts[46]
C278QQQQQQAdapt to mammalian hosts[46]
E382DDDEEEHuman host marker[40,47]
S/A515TTTTTTIncreased polymerase activity, increased virulence in mammals and birds[33,39]
L653PSPPPPAdapt to mammalian hosts[46]
NPV41IIIIIIMight contribute to viral transmissibility[48]
V105MMVMIIContribute to the increased virulence[49]
F253IIIIIIResults in attenuated pathogenicity of the virus in mice[50]
I353VVVVVVIncreased virulence in mice[32]
MPN30DDDDDDIncreased virulence in mammals[51]
T215AAAAAAIncreased virulence in mammals[52]
NS1A/P42SSSSSSIncreased virulence in mammals, antagonism of IFN induction[52,53]
T/D/V/R/A127NNNNNNIncreased virulence in mammals[54]
V149AAAAAAPathogenicity in mice, antagonism of IFN induction[55]
* H3 numbering. K6-2018: A/Bean goose/South Korea/KNU18-6/2018(H6N5) K17-2017: A/Aix galericulata/South Korea/K17-1638-5/2017 (H6N5) CN5-2009: A/aquatic bird/Korea/CN5/2009 (H6N5) W69-2005: A/aquatic bird/Korea/W69/2005 (H6N5) TW02/2013: A/Taiwan/2/2013 (H6N1) infected-human isolates. The letters are amino acid abbreviation. A, Alanine; E, Glutamate; G, Glycine; Q, Glutamine; L, Leucine; T, Threonine; N, Asparagine; I, Isoleucine; V, Valine.

3.3. Viral Replication Kinetics in MDCK Cells

To address the potentially increased virulence of the novel H6N5 originating from North America, the viral replication capability of our isolate in mammalian cells was investigated relative to that of H1N1, H7N1, and H9N2 strains in vitro (Figure 2).
Results showed that the K6 virus grew efficiently and reached a peak at 48 hpi. The same kinetic curve pattern was observed in the control H7N1 strain. However, the titer of the H7N1 virus was approximately 10–30-fold higher than that of H6N5 at 36 hpi. In contrast, the H9N2 virus grew slowly and showed the lowest titer from the beginning to 12, 24, and 36 hpi. Raw data for the TCID50 assay are presented in Figure S4.

3.4. Animal Study

To determine the pathogenic potential of the new isolate in mammals, 105,106, and 107 EID50/ 50 µl of H6N5 viruses were inoculated into six-week-old female BALB/c mice. Overall, the infected mice displayed no severe clinical signs, including ruffled fur, depression, labored breathing, and/or severe weight loss during the 14 days post infection (dpi). Although increasing titers of the H6N5 virus infected to the mice from 105 to 106 to 107 EID50, it did cause insignificant changes in the total body weight during the 14 dpi period (Figure 3A).
For comparison, control H1N1 and H7N1 strains were inoculated at 105 EID50 /mouse. The body weight of the mice infected with these viruses was found to decrease regularly, with the lowest weight (86.81% and 86.08%) observed at approximately day 7–8 post challenge for H1N1 and H7N1, respectively. Meanwhile, the group infected with the K6-H6N5 isolate, retained stable bodyweight throughout the 14 days. Moreover, from day 9 the mice infected with all three strains recovered their body weight (Figure 3B). As a result, no significant differences were observed between the groups of mice in survival rate (Figure 3C). The viral load shedding in lungs after 3, 6, and 14 dpi has been represented in Figure 3D and Figure S5. K6-H6N5 replicated in the lung with the lowest titer (1.54 ± 0.08 log10 TICD50/mL) at 3 dpi, which then increased slightly to 1.88 ± 0.21 log10 TICD50/mL at 6 dpi, and was absent at 14 dpi; whereas H1N1 and H7N1 were maintained at a comparatively high titer at 3 dpi (3.96 ± 0.31; 3.79 ± 0.26 log10 TICD50/mL, respectively) and at 6 dpi (4.25 ± 0.37; 4.84 ± 0.19 log10 TICD50/mL, respectively). On day 14, there was no viral shedding in the lungs of all groups.
The mouse infected with K6 did not show any considerable change in the body weight or the survival rate; however, it exhibits the viral shedding in mice lungs of those at 3; 6 dpi, the lung histopathology of infected mice were examined. Hematoxylin and eosin (H&E) stained section of the K6, H7N1, H1N1-infected lung at 6 dpi showing an extensive penetrate of neutrophils that fill the alveolar air spaces (Figure 3E). K6-H6N5 virus-infected lungs presented cellular consolidation compared with those of the normal (uninfected) mice at 6 dpi (Figure S6) although lungs of K6-H6N5-infected mice did not show a significant increase in weight whereas those of H7N1- or H1N1-infection increased a significant increase of lung weight (Figure 3F). Taken together, it indicated that K6-H6N5 has a relatively lower adaptation to mice compared to H7N1 and H1N1.

3.5. Hypothesis for Reassortment Event of A/Bean goose/South Korea/KNU18-6/2018(H6N5)

K17 was identified during September 2017–March 2018 when a K6 strain was isolated in March 2018, and the locations of K17 (37°43′37″ N, 126°39′54″ E) (Figure S7) [8] and K6 were placed within a radius of 4.27 km (Figure 4).
During 2015–2017 in North America, the HA and NA of A/California/H6N5 may have reassorted with other internal gene segments of local AIV residents, thereby transmitting it to Alaska to recruit PB2 for K17, and gathering both PB1 and PB2 to K6 from a donor virus A/Mallard/Alaska/Ah0029066S.2.A/2016 H12N5.
The intercontinental transmission may have been facilitated by certain key vector species, such as the northern pintail and mallard, which have wide geographic distribution across Europe, Palearctic, and North America. Somewhere in Russia, near Alaska, the North American-origin H6N5 may have spread to bean geese before being introduced to Korea as K6 (Figure 5B). However, Aix galericulata as a host of K17 is an Asian population that breeds, feeds, and overwinters in lowland eastern China and southern Japan. Therefore, K17 has been suggested to be transferred from migratory birds to other native species before spreading to mandarin ducks in Korea (Figure 5A).

4. Discussion

The H5 HPAI viruses has been spreading widely across several countries, providing an increasing opportunity for the Asian-origin viruses to be reassorted in the Alaska area before subsequently being transferred to Canada and other states in the USA, rather than being created by multiple reassortments within the local resident LPAI strains. Therefore, the emergence of intercontinental distribution of AI has been tightly monitored to explain the high variations in AIV.
Alaska is a feasible location for intercontinental exchanges, since various AIV of Eurasian and North American lineages have been identified there. Several subtypes, such as H8N4, H9N2, H13, H16, H3, and N7, found in Alaska, contain mixed Eurasian-origin gene segments [56,57,58,59].
The HPAI H5N8 epidemic that occurred during January 2014 in South Korea, infecting 161 commercial poultry flocks, led to the culling of 14,000,000 birds; subsequently, the virus was detected continuously in Germany, England, and the Netherlands in early November 2014. Soon afterward, in early 2015, HPAI H5 had threatened Canada, along with various states in the USA, with more than 200 confirmed cases and 48,000,000 birds affected [60]. Certain HPAI H5 clade 2.3.4.4 viruses isolated from North America were proven to possess the HA gene originating from Asian lineages [61,62].
To date, intercontinental reassortment of several AIV genes in Asian waterfowl has been described in Japan (H4N6) [63] and Korea (H3N8 and H4N6), resulting in the eventual transmission to live animal markets [64]. Moreover, the imported single NS gene segment of AIV from North American to South Korea was recorded between 2012 and 2017 to have occurred via a greater white-fronted goose (Anser albifrons) [7], whereas the full K17-H6N5 genome was introduced in 2018 [8].
Phylogenetic tree analysis indicated eight K17 gene segments to share full sequence homology with North American AIV strains, during 2016–2017, rather than with those of Eurasian ancestry (Figure S2). However, the host for K17 was determined to be the mandarin duck [8], which shows narrow distribution in the Far East and Southeast Asia, within the East Asian–Australasian flyway, and is, therefore, not likely to move between the two continents of Eurasia and North America [65].
In contrast, the host of K6 was the bean goose (Anser fabalis), which is widely distributed across Asia, Europe, and Russia for both breeding and feeding [65]. Moreover, an overlap area has been identified in Russia, near Alaska, where the migratory birds from North America, such as northern pintail or mallard ducks, share the same distribution as bean goose. Indeed, northern pintails breed across Poland, Russia, Mongolia, Canada, Alaska, and Midwestern United States, and winter primarily at the equator, in Panama, northern sub-Saharan Africa, and tropical South Asia [65]. They have also been detected as the major hosts of several AIV subtypes (H1–8, 10–12, and N1–9), isolated over three years (2006–2008) in Alaska, with seven of their eight gene segments acquired as a result of intercontinental transfer from Asia [6]. Therefore, considering the close habitats of bean geese (host of K6) and mandarin ducks (host of K17) over the same seasonal period in Korea (Figure 4), we hypothesized that the introduction of K17 might be associated with that of K6 to Korea. Although it is possible that K6 and K17 were independently transported, due to the intersection of the K6 and K17 host near Alaska, reassortment may occur during groupings in these similar habitats as well.
This hypothesis was supported by records acquired by the National Institute of Biological Resources from two different satellite tracking devices (br1691 and br1692) monitoring the movement of bean geese during January to November 2017 in Korea. According to the data, in the summer, the bean goose migrated to Russia (near Alaska) from Korea and returned to Korea in November 2017. Furthermore, the migration pattern of the bean goose, with the br1691 device attached, corresponded well with the isolation site identified for K17 during October–November 2017 (Figure S8A), while that of the K6 virus corresponded to the migratory path of the bean goose with br1692 attached, in November 2017, i.e., four months prior to the isolation of K6 in this study (Figure S8B). This data, therefore, supported the theory for possible reassortment of North American-origin virus in Russia, near Alaska, and potential transmission of the K17 virus from a bean goose to mandarin duck during the summer months, before they are returned to Korea, or possible reassortment via cross grouping within Korea (Figure 5B).
To date, gene segments from the novel Eurasian, highly pathogenic, H5 virus in Washington, USA, in 2014, have been found to be genetically similar to that of A/bean goose/Korea/H40/2014 [62]. Therefore, bean geese warrant intensive surveillance for intercontinental spread of Asian- or North American-origin AIV through Beringia in Alaska.
Considering the lack of satellite tracking data for bean geese in Alaska, their migration from North America directly to Asia may still be possible. Bean geese could join other migratory birds and transmit AIV to Eurasia and North America via flyway overlap with wild birds (East Asia, Pacific Americas, and Mississippi Americas), thus indicating a demand for enhancing surveillance in these regions.
The H6 subtype reportedly infects various mammals, including humans. A study examining the antigenic and genetic relatedness of 14 H6 influenza viruses reported efficient replication of an H6 subtype in mice and ferrets without preadaptation [66]. The first report of human infection with a wild H6N1 was from Taiwan [67]. This virus has since been circulating throughout both North America and Eurasia [68]. However, an important G228S (H3 numbering) substitution at the receptor-binding site of HA, as well as a 14-amino acid deletion in the stalk region of NA noted in a human-infecting H6N1 strain, were not identified in the K6 or K17 isolates (Table 2). These results suggest the low pathogenicity of K6 compared to that of the H6 subtype in Taiwan.
Meanwhile, gene reassortment has permitted the H6 AIV to serve as a source of novel HPAI, facilitating the circulation of strains that had previously threatened humans. For instance, during 2014–2016 in China, the novel H5N6 HPAI was reported as lethal in chickens, and moderately pathogenic in mice, leading to several cases of infection and deaths in humans [69,70]. The whole genome of these H5N6 strains exhibited high similarity with those from human isolates and their NA genes originating from H6N6. Thus, it demands the phylogenetic tree analysis, as well as the study replication ability in mammalian systems of new isolates.
According to phylogenetic tree analyses, the previously reported K17 and our current K6 possess a highly similar genome sequence; however, the adaptation in mammals or mice was not characterized for the K17 strain. In our study, K6 showed the highest viral titer of 106.37 log10 TCID50 at 36–48 hpi, which corresponded with the previous results that had shown H6 viruses to grow efficiently and reach a peak of 105 log10 TCID50 at 36 hpi in MDCK cells [9]. Moreover, the K6 virus was found to replicate efficiently in lungs without prior adaptation; however, it was not lethal in mice, which is consistent with other publications that examined the toxicity in mice, induced by H6-AIV [71]. The mice challenged intranasally with 106 EID50 exhibited slight loss in body weight, and viral load in the lung was found to be 4.82 log10 TCID50 /100 mg. Alternatively, the A/aquatic bird/Korea/CN5/2009 (H6N5) strain causes illness and lethality in mice, with severe clinical signs, including significant weight loss, ruffled fur, lethargy, and ataxia, with 100% mortality noted within 7 days [12]. According to the molecular analysis of gene segments between K6 and A/aquatic bird/Korea/CN5/2009 (H6N5), the D701N mutation in PB2 may contribute to the increased virulence observed for H6N5 in the mammalian host (Table 2).
Although the K6 or K17 are low pathogenic H6N5 strains imported to Korea with a full genome from North America, they pose a potential risk for increasing pathology in mice and humans via participating in reassortment with local AIV residents. Particularly, the PB2 and/or PA gene of the CN5-2009 strain harboring D701N and L653P substitutions, respectively, can become reassorted with other K6 genes, resulting in potentially increased virulence in mammals. Additionally, K6 may increase the possibility of host jumping by providing the PB2 with A274T mutation or/and the PA gene with E382D alteration to CN5-2009 or other Korea strains reassortment (Table 2).
For comparison, H7N1 and H1N1 were used as positive controls of mouse-adapted strains while H9N2 was a non-mouse-adapted strain. Results showed a kinetic peak for H7N1 at 107.7 log10 TCID50 at 48 hpi, whereas another publication showed a peak at 105 log10 TCID50 [72]. Besides, mice challenged with the H7N1 strain displayed 30% weight loss at 8 dpi and only 60% survival by day 15. Meanwhile, our data showed 100% survival in mice with H7N1 at 15 dpi, with only 20% loss in body weight. This difference may be derived from genomic variation of H7 strains.
Overall, our results contributed to furthering the current understanding regarding how North American-origin H6N5 in mandarin duck was introduced into Korea, based on the genetic identity and migration pattern of bean geese during the isolation season of K6 and K17. However, a limitation of our study was that we only suggested the potential import of H6N5 by bean geese to Korean mandarin ducks, considering the database for the spatial migratory patterns of bean geese. Additional influenza surveillance in bean geese, along the flyways in Alaska, might provide useful inferences for early introduction of novel strains into Korea. Further studies should, therefore, focus on collecting direct AIV samples from bean geese with satellite tracking devices attached.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/7/774/s1. Figure S1, host information by COI gene sequence (A) and information of bird species identification (B). Figure S2, phylogenetic tree of whole genome of the K6 avian influenza virus (AIV) compared to influenza sequences available in NCBI and GISAID. Figure S3, maximum likelihood phylogenetic trees showing inferred relationships among nucleotide sequences for the complete coding regions of K6 and K17 gene segments. Figure S4, raw ELISA data to conduct TCID50 assay. Figure S5, raw ELISA data to conduct TCID50 assay to measure virus K6-H6N5, H1N1, H7N1 titer in lung at day 3 (A); day 6 (B); day 14 (C) post infected. Figure S6, lungs from (A) uninfected control; (B) K6-H6N5; (C) H7N1; (D) H7N1 -infected mouse at 6 dpi. Scale bar: 0.5cm. Figure S7, location map of mandarin duck marked with satellite transmitters in Korea in October in 2017. Figure S8, spatial distribution patterns of Bean geese during Jan 2017- Nov 2017. Table S1, information of strains with more than 99% genetic homology.

Author Contributions

Conceptualization, H.P. and S.-J.Y.; Methodology, N.M.N.; Software, N.M.N.; Validation, N.M.N.; Resources, H.W.S.; Investigation, K.-J.Y.; Writing–Original Draft Preparation, N.M.N.; Writing–Review and Editing, S.-J.Y.; Visualization, N.M.N.; Supervision, H.P. and S.-J.Y.; Project Administration, S.-J.Y.; Funding Acquisition, H.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Priority Research Centers Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, (NRF-2015R1A6A1A03032236) and the Bio and Medical Technology Development Program of the National Research Foundation (NRF) and funded by the Korean government (MSIT) (NRF-2018M3A9H4055193).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Prevention, C.D.C. Avian Influenza Type A Viruses. Available online: https://www.cdc.gov/flu/avianflu/influenza-a-virus-subtypes.htm (accessed on 19 April 2017).
  2. Olsen, B.; Munster, V.; Wallensten, A.; Waldenström, J.; Osterhaus, A.; Fouchier, R.A.M. Global patterns of influenza a virus in wild birds. Science 2006, 312, 384–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lam, T.T.-Y.; Ip, H.S.; Ghedin, E.; Wentworth, D.; Halpin, R.A.; Stockwell, T.B.; Spiro, D.J.; Dusek, R.J.; Bortner, J.B.; Hoskins, J.; et al. Migratory flyway and geographical distance are barriers to the gene flow of influenza virus among North American birds. Ecol. Lett. 2012, 15, 24–33. [Google Scholar] [CrossRef] [Green Version]
  4. Winker, K.; McCracken, K.G.; Gibson, D.D.; Pruett, C.L.; Meier, R.; Huettmann, F.; Wege, M.; Kulikova, I.V.; Zhuravlev, Y.N.; Perdue, M.L.; et al. Movements of birds and avian influenza from Asia into Alaska. Emerg. Infect. Dis. 2007, 13, 547–552. [Google Scholar] [CrossRef]
  5. Pearce, J.M.; Ramey, A.M.; Flint, P.L.; Koehler, A.; Fleskes, J.P.; Franson, J.; Hall, J.S.; Derksen, D.V.; Ip, H.S. Avian influenza at both ends of a migratory flyway: Characterizing viral genomic diversity to optimize surveillance plans for North America. Evol. Appl. 2009, 2, 457–468. [Google Scholar] [CrossRef] [PubMed]
  6. Ramey, A.M.; Pearce, J.M.; Flint, P.L.; Ip, H.S.; Derksen, D.V.; Franson, J.; Petrula, M.J.; Scotton, B.D.; Sowl, K.M.; Wege, M.L.; et al. Intercontinental reassortment and genomic variation of low pathogenic avian influenza viruses isolated from northern pintails (Anas acuta) in Alaska: Examining the evidence through space and time. Virology 2010, 401, 179–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Cheon, S.-H.; Lee, Y.-N.; Kang, S.-I.; Kye, S.-J.; Lee, E.-K.; Heo, G.-B.; Lee, M.-H.; Kim, J.W.; Lee, K.-N.; Son, H.-M.; et al. Genetic evidence for the intercontinental movement of avian influenza viruses possessing North American-origin nonstructural gene allele B into South Korea. Infect. Genet. Evol. 2018, 66, 18–25. [Google Scholar] [CrossRef]
  8. Jeong, S.; Lee, D.-H.; Kim, Y.-J.; Lee, S.-H.; Cho, A.Y.; Noh, J.-Y.; Tseren-Ochir, E.-O.; Jeong, J.-H.; Song, C.-S. Introduction of Avian Influenza A (H6N5) Virus into Asia from North America by Wild Birds. Emerg. Infect. Dis. 2019, 25, 2138–2140. [Google Scholar] [CrossRef]
  9. Zou, S.; Gao, R.; Zhang, Y.; Li, X.; Chen, W.; Bai, T.; Dong, L.; Wang, D.; Shu, Y. Molecular characterization of H6 subtype influenza viruses in southern China from 2009 to 2011. Emerg. Microbes Infect. 2016, 5, e73. [Google Scholar] [CrossRef] [Green Version]
  10. Peng, C.; Sun, H.; Li, J.; Hou, G.; Wang, S.; Liu, S.; Zhuang, Q.; Cheng, S.; Chen, J.; Jiang, W. Molecular epidemiological survey and complete genomic phylogenetic analysis of H6 subtype avian influenza viruses in poultry in China from 2011 to 2016. Infect. Genet. Evol. 2018, 65, 91–95. [Google Scholar] [CrossRef]
  11. Li, J.; Quan, C.; Xie, Y.; Ke, C.; Nie, Y.; Chen, Q.; Hu, T.; Chen, J.; Wong, G.; Wang, Q.; et al. Continued reassortment of avian H6 influenza viruses from Southern China, 2014–2016. Transbound. Emerg. Dis. 2019, 66, 592–598. [Google Scholar] [CrossRef] [Green Version]
  12. Nam, J.-H.; Kim, E.-H.; Song, D.; Choi, Y.K.; Kim, J.-K.; Poo, H. Emergence of mammalian species-infectious and -pathogenic avian influenza H6N5 virus with no evidence of adaptation. J. Virol. 2011, 85, 13271–13277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tuong, H.T.; Nguyen, N.M.; Sung, H.W.; Park, H.; Yeo, S.-J. Genetic Characterization of Avian Influenza A (H11N9) Virus Isolated from Mandarin Ducks in South Korea in 2018. Viruses 2020, 12, 203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Killian, M.L. Hemagglutination assay for influenza virus. Methods Mol. Biol. 2014, 1161, 3–9. [Google Scholar] [PubMed]
  15. World Health Organization. Geneva Switzerland. Available online: https://www.who.int/influenza/gisrs_laboratory/WHO_information_for_the_molecular_detection_of_influenza_viruses_20171023_Final (accessed on 23 November 2017).
  16. Hebert, P.D.; Stoeckle, M.Y.; Zemlak, T.S.; Francis, C.M. Identification of birds through DNA barcodes. PLoS Biol. 2004, 2, e312. [Google Scholar] [CrossRef] [Green Version]
  17. Yeo, S.-J.; Than, D.-D.; Park, H.-S.; Sung, H.W.; Park, H. Molecular Characterization of a Novel Avian Influenza A (H2N9) Strain Isolated from Wild Duck in Korea in 2018. Viruses 2019, 11, 1046. [Google Scholar] [CrossRef] [Green Version]
  18. Ambardar, S.; Gupta, R.; Trakroo, D.; Lal, R.; Vakhlu, J. High Throughput Sequencing: An. Overview of Sequencing Chemistry. Indian J. Microbiol. 2016, 56, 394–404. [Google Scholar] [CrossRef] [Green Version]
  19. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  20. WHO. Manual for the Laboratory Diagnosis and Virological Surveillance of Influenza; World Health Organization Global Influenza Surveillance Network: Geneva, Switzerland, 2011. [Google Scholar]
  21. Reed, L.; Muench, H. A simple method of estimating fifty per cent endpoints. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  22. Jimenez-Bluhm, P.; Karlsson, E.A.; Ciuoderis, K.A.; Cortez, V.; Marvin, S.A.; Hamilton-West, C.; Schultz-Cherry, S.; Osorio, J.E. Avian H11 influenza virus isolated from domestic poultry in a Colombian live animal market. Emerg. Microbes Infect. 2016, 5, e121. [Google Scholar] [CrossRef]
  23. The University of IOWA. 2018. Available online: https://animal.research.uiowa.edu/iacuc-guidelines-anesthesia (accessed on 6 September 2018).
  24. Diederich, S.; Berhane, Y.; Embury-Hyatt, C.; Hisanaga, T.; Handel, K.; Cottam-Birt, C.; Ranadheera, C.; Kobasa, D.; Pasick, J. Hemagglutinin-Neuraminidase Balance Influences the Virulence Phenotype of a Recombinant H5N3 Influenza A Virus Possessing a Polybasic HA0 Cleavage Site. J. Virol. 2015, 89, 10724–10734. [Google Scholar] [CrossRef] [Green Version]
  25. Lazniewski, M.; Dawson, W.K.; Szczepińska, T.; Plewczynski, D. The structural variability of the influenza A hemagglutinin receptor-binding site. Brief. Funct. Genom. 2018, 17, 415–427. [Google Scholar]
  26. Matsuoka, Y.; Swayne, D.E.; Thomas, C.; Rameix-Welti, M.-A.; Naffakh, N.; Warnes, C.; Altholtz, M.; Donis, R.; Subbarao, K. Neuraminidase Stalk Length and Additional Glycosylation of the Hemagglutinin Influence the Virulence of Influenza H5N1 Viruses for Mice. J. Virol. 2009, 83, 4704–4708. [Google Scholar] [CrossRef] [Green Version]
  27. Choi, W.-S.; Jeong, J.H.; Kwon, J.J.; Ahn, S.J.; Lloren, K.K.S.; Kwon, H.-I.; Chae, H.B.; Hwang, J.; Kim, M.H.; Kim, C.-J.; et al. Screening for Neuraminidase Inhibitor Resistance Markers among Avian Influenza Viruses of the N4, N5, N6, and N8 Neuraminidase Subtypes. J. Virol. 2017, 92, e01580-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tan, L.; Su, S.; Smith, D.K.; He, S.; Zheng, Y.; Shao, Z.; Ma, J.; Zhu, H.; Zhang, G. A Combination of HA and PA Mutations Enhances Virulence in a Mouse-Adapted H6N6 Influenza A Virus. J. Virol. 2014, 88, 14116–14125. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, Z.; Cheng, K.; Xin, Y.; Sun, W.; Li, X.; Huang, J.; Zhang, K.; Yang, S.; Wang, T.; Zheng, X.; et al. Multiple amino acid substitutions involved in the adaptation of H6N1 avian influenza virus in mice. Vet. Microbiol. 2014, 174, 316–321. [Google Scholar] [CrossRef] [PubMed]
  30. Li, J.; Li, Y.; Hu, Y.; Chang, G.; Sun, W.; Yang, Y.; Kang, X.; Wu, X.; Zhu, Q. PB1-mediated virulence attenuation of H5N1 influenza virus in mice is associated with PB2. J. Gen. Virol. 2011, 92 Pt 6, 1435–1444. [Google Scholar] [CrossRef]
  31. Li, J.; Ishaq, M.; Prudence, M.; Xi, X.; Hu, T.; Liu, Q.; Guo, D. Single mutation at the amino acid position 627 of PB2 that leads to increased virulence of an H5N1 avian influenza virus during adaptation in mice can be compensated by multiple mutations at other sites of PB2. Virus Res. 2009, 144, 123–129. [Google Scholar] [CrossRef] [PubMed]
  32. Prokopyeva, E.A.; Sobolev, I.A.; Prokopyev, M.V.; Shestopalov, A.M. Adaptation of influenza A (H1N1)pdm09 virus in experimental mouse models. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2016, 39, 265–271. [Google Scholar] [CrossRef]
  33. Leung, B.W.; Chen, H.; Brownlee, G.G. Correlation between polymerase activity and pathogenicity in two duck H5N1 influenza viruses suggests that the polymerase contributes to pathogenicity. Virology 2010, 401, 96–106. [Google Scholar] [CrossRef] [Green Version]
  34. Salomon, R.; Franks, J.; Govorkova, E.A.; Ilyushina, N.A.; Yen, H.L.; Hulse-Post, D.J.; Humberd, J.; Trichet, M.; Rehg, J.E.; Webby, R.J.; et al. The polymerase complex genes contribute to the high virulence of the human H5N1 influenza virus isolate A/Vietnam/1203/04. J. Exp. Med. 2006, 203, 689–697. [Google Scholar] [CrossRef] [Green Version]
  35. Govorkova, E.A.; Rehg, J.E.; Krauss, S.; Yen, H.-L.; Guan, Y.; Peiris, M.; Nguyen, T.D.; Hanh, T.H.; Puthavathana, P.; Long, H.T.; et al. Lethality to ferrets of H5N1 influenza viruses isolated from humans and poultry in 2004. J. Virol. 2005, 79, 2191–2198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gabriel, G.; Abram, M.; Keiner, B.; Wagner, R.; Klenk, H.D.; Stech, J. Differential polymerase activity in avian and mammalian cells determines host range of influenza virus. J. Virol. 2007, 81, 9601–9604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Gabriel, G.; Dauber, B.; Wolff, T.; Planz, O.; Klenk, H.-D.; Stech, J. The viral polymerase mediates adaptation of an avian influenza virus to a mammalian host. Proc. Natl. Acad. Sci. USA 2005, 102, 18590–18595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Ágnes, F.; Loucaides, E.M.; Mura, M.; Wise, H.M.; Barclay, W.S.; Digard, P. Influence of PB2 host-range determinants on the intranuclear mobility of the influenza A virus polymerase. J. Gen. Virol. 2011, 92 Pt 7, 1650–1661. [Google Scholar]
  39. Hulse-Post, D.J.; Franks, J.; Boyd, K.; Salomon, R.; Hoffmann, E.; Yen, H.-L.; Webby, R.J.; Walker, D.; Nguyen, T.D.; Webster, R.G. Molecular changes in the polymerase genes (PA and PB1) associated with high pathogenicity of H5N1 influenza virus in mallard ducks. J. Virol. 2007, 81, 8515–8524. [Google Scholar] [CrossRef] [Green Version]
  40. Taubenberger, J.K.; Reid, A.H.; Lourens, R.M.; Wang, R.; Jin, G.; Fanning, T.G. Characterization of the 1918 influenza virus polymerase genes. Nature 2005, 437, 889–893. [Google Scholar] [CrossRef]
  41. Wei, K.; Sun, H.; Sun, Z.; Sun, Y.; Kong, W.; Pu, J.; Ma, G.; Yin, Y.; Yang, H.; Guo, X.; et al. Influenza A virus acquires enhanced pathogenicity and transmissibility after serial passages in swine. J. Virol. 2014, 88, 11981–11994. [Google Scholar] [CrossRef] [Green Version]
  42. Xu, C.; Hu, W.-B.; Xu, K.; He, Y.-X.; Wang, T.-Y.; Chen, Z.; Li, T.-X.; Liu, J.-H.; Buchy, P.; Sun, B. Amino acids 473V and 598P of PB1 from an avian-origin influenza A virus contribute to polymerase activity, especially in mammalian cells. J. Gen. Virol. 2012, 93 Pt 3, 531–540. [Google Scholar] [CrossRef]
  43. Schmolke, M.; Manicassamy, B.; Pena, L.; Sutton, T.; Hai, R.; Varga, Z.T.; Hale, B.G.; Steel, J.; Perez, D.R.; García-Sastre, A. Differential contribution of PB1-F2 to the virulence of highly pathogenic H5N1 influenza A virus in mammalian and avian species. PLoS Pathog. 2011, 7, e1002186. [Google Scholar] [CrossRef] [Green Version]
  44. Ozawa, M.; Basnet, S.; Burley, L.M.; Neumann, G.; Hatta, M.; Kawaoka, Y. Impact of amino acid mutations in PB2, PB1-F2, and NS1 on the replication and pathogenicity of pandemic (H1N1) 2009 influenza viruses. J. Virol. 2011, 85, 4596–4601. [Google Scholar] [CrossRef] [Green Version]
  45. Yamayoshi, S.; Yamada, S.; Fukuyama, S.; Murakami, S.; Zhao, D.; Uraki, R.; Watanabe, T.; Tomita, Y.; Macken, C.; Neumann, G.; et al. Virulence-affecting amino acid changes in the PA protein of H7N9 influenza A viruses. J. Virol. 2014, 88, 3127–3134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Mei, K.; Liu, G.; Chen, Z.; Gao, Z.; Zhao, L.; Jin, T.; Yu, X.; Chen, Q. Deep sequencing reveals the viral adaptation process of environment-derived H10N8 in mice. Infect. Genet. Evol. 2016, 37, 8–13. [Google Scholar] [CrossRef] [PubMed]
  47. Shaw, M.; Cooper, L.; Xu, X.; Thompson, W.; Krauss, S.; Guan, Y.; Zhou, N.; Klimov, A.; Cox, N.; Webster, R.; et al. Molecular changes associated with the transmission of avian influenza a H5N1 and H9N2 viruses to humans. J. Med. Virol. 2002, 66, 107–114. [Google Scholar] [CrossRef] [PubMed]
  48. Zhu, W.; Zou, X.; Zhou, J.; Tang, J.; Shu, Y. Residues 41V and/or 210D in the NP protein enhance polymerase activities and potential replication of novel influenza (H7N9) viruses at low temperature. Virol. J. 2015, 12, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Liu, Q.; Chen, H.; Huang, J.; Chen, Y.; Gu, M.; Wang, X.; Hu, S.; Liu, X.; Liu, X. A nonpathogenic duck-origin H9N2 influenza A virus adapts to high pathogenicity in mice. Arch. Virol. 2014, 159, 2243–2252. [Google Scholar] [CrossRef] [PubMed]
  50. Li, J.; Zheng, W.; Hou, L.; Chen, C.; Fan, W.; Qu, H.; Jiang, J.; Liu, J.; Gao, G.F.; Zhou, J.; et al. Differential nucleocytoplasmic shuttling of the nucleoprotein of influenza a viruses and association with host tropism. Cell. Microbiol. 2017, 19, e12692. [Google Scholar] [CrossRef] [Green Version]
  51. Ou, X.; Chen, F.; Zhang, R.; Chen, J.; Liu, R.; Sun, B. Analysis of the full-length genome of a novel strain of the H7N9 avian influenza virus. Exp. Ther. Med. 2014, 7, 1369–1375. [Google Scholar] [CrossRef] [Green Version]
  52. Fan, S.; Deng, G.; Song, J.; Tian, G.; Suo, Y.; Jiang, Y.; Guan, Y.; Bu, Z.; Kawaoka, Y.; Chen, H. Two amino acid residues in the matrix protein M1 contribute to the virulence difference of H5N1 avian influenza viruses in mice. Virology 2009, 384, 28–32. [Google Scholar] [CrossRef] [Green Version]
  53. Jiao, P.; Tian, G.; Li, Y.; Deng, G.; Jiang, Y.; Liu, C.; Liu, W.; Bu, Z.; Kawaoka, Y.; Chen, H. A single-amino-acid substitution in the NS1 protein changes the pathogenicity of H5N1 avian influenza viruses in mice. J. Virol. 2008, 82, 1146–1154. [Google Scholar] [CrossRef] [Green Version]
  54. Min, J.-Y.; Li, S.; Sen, G.C.; Krug, R.M. A site on the influenza A virus NS1 protein mediates both inhibition of PKR activation and temporal regulation of viral RNA synthesis. Virology 2007, 363, 236–243. [Google Scholar] [CrossRef] [Green Version]
  55. Li, Z.; Jiang, Y.; Jiao, P.; Wang, A.; Zhao, F.; Tian, G.; Wang, X.; Yu, K.; Bu, Z.; Chen, H. The NS1 gene contributes to the virulence of H5N1 avian influenza viruses. J. Virol. 2006, 80, 11115–11123. [Google Scholar] [CrossRef] [Green Version]
  56. Ramey, A.M.; Reeves, A.B.; Donnelly, T.; Poulson, R.L.; Stallknecht, D.E. Introduction of Eurasian-Origin Influenza A(H8N4) Virus into North America by Migratory Birds. Emerg. Infect. Dis. 2018, 24, 1950–1953. [Google Scholar] [CrossRef] [PubMed]
  57. Ramey, A.M.; Reeves, A.B.; Sonsthagen, S.A.; TeSlaa, J.; Nashold, S.; Donnelly, T.; Casler, B.; Hall, J.S. Dispersal of H9N2 influenza A viruses between East Asia and North America by wild birds. Virology 2015, 482, 79–83. [Google Scholar] [CrossRef] [Green Version]
  58. Wille, M.; Robertson, G.J.; Whitney, H.; Bishop, M.A.; Runstadler, J.A.; Lang, A.S. Extensive geographic mosaicism in avian influenza viruses from gulls in the northern hemisphere. PLoS ONE 2011, 6, e20664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Ramey, A.M.; Pearce, J.M.; Reeves, A.B.; Poulson, R.L.; Dobson, J.; Lefferts, B.; Spragens, K.; Stallknecht, D.E. Surveillance for Eurasian-origin and intercontinental reassortant highly pathogenic influenza A viruses in Alaska, spring and summer 2015. Virol. J. 2016, 13, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Nonthabenjawan, N.; Cardona, C.J.; Amonsin, A.; Sreevatsan, S. Time-space analysis of highly pathogenic avian influenza H5N2 outbreak in the US. Virol. J. 2016, 13, 147. [Google Scholar] [CrossRef] [Green Version]
  61. Saito, T.; Tanikawa, T.; Uchida, Y.; Takemae, N.; Kanehira, K.; Tsunekuni, R. Intracontinental and intercontinental dissemination of Asian H5 highly pathogenic avian influenza virus (clade 2.3.4.4) in the winter of 2014-2015. Rev. Med. Virol. 2015, 25, 388–405. [Google Scholar] [CrossRef]
  62. Ip, H.S.; Kim-Torchetti, M.; Crespo, R.; Kohrs, P.; Debruyn, P.; Mansfield, K.G.; Baszler, T.; Badcoe, L.; Bodenstein, B.; Shearn-Bochsler, V.; et al. Novel Eurasian highly pathogenic avian influenza A H5 viruses in wild birds, Washington, USA, 2014. Emerg. Infect. Dis. 2015, 21, 886–890. [Google Scholar] [CrossRef]
  63. Okuya, K.; Kanazawa, N.; Kanda, T.; Kuwahara, M.; Matsuu, A.; Horie, M.; Toda, S.; Ozawa, M.; Masatani, T. Genetic characterization of an avian H4N6 influenza virus isolated from the Izumi plain, Japan. Microbiol. Immunol. 2017, 61, 513–518. [Google Scholar] [CrossRef] [Green Version]
  64. Lee, D.-H.; Lee, H.-J.; Lee, Y.-N.; Park, J.-K.; Lim, T.-H.; Kim, M.-S.; Youn, H.-N.; Lee, J.-B.; Park, S.-Y.; Choi, I.-S.; et al. Evidence of intercontinental transfer of North American lineage avian influenza virus into Korea. Infect. Genet. Evol. 2011, 11, 232–236. [Google Scholar] [CrossRef] [PubMed]
  65. The Global Register of Migratory Species. Available online: http://www.groms.de/groms/JPGs/Small_Ones/index.html (accessed on 30 April 2019).
  66. Gillim-Ross, L.; Santos, C.; Chen, Z.; Aspelund, A.; Yang, C.-F.; Ye, D.; Jin, H.; Kemble, G.; Subbarao, K. Avian influenza h6 viruses productively infect and cause illness in mice and ferrets. J. Virol. 2008, 82, 10854–10863. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Wei, S.-H.; Yang, J.-R.; Wu, H.-S.; Chang, M.-C.; Lin, J.-S.; Lin, C.-Y.; Liu, Y.-L.; Lo, Y.-C.; Yang, C.-H.; Chuang, J.-H.; et al. Human infection with avian influenza A H6N1 virus: An epidemiological analysis. Lancet Respir. Med. 2013, 1, 771–778. [Google Scholar] [CrossRef]
  68. Shi, W.; Shi, Y.; Wu, Y.; Liu, D.; Gao, G.F. Origin and molecular characterization of the human-infecting H6N1 influenza virus in Taiwan. Protein Cell 2013, 4, 846–853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Sun, W.; Li, J.; Hu, J.; Jiang, D.; Xing, C.; Zhan, T.; Liu, X. Genetic analysis and biological characteristics of different internal gene origin H5N6 reassortment avian influenza virus in China in 2016. Vet. Microbiol. 2018, 219, 200–211. [Google Scholar] [CrossRef]
  70. Jiao, P.; Cui, J.; Song, Y.; Song, H.; Zhao, Z.; Wu, S.; Qu, N.; Wang, N.; Ouyang, G.; Jiao, P. New Reassortant H5N6 Highly Pathogenic Avian Influenza Viruses in Southern China, 2014. Front. Microbiol. 2016, 7, 754. [Google Scholar] [CrossRef]
  71. Hu, C.; Li, X.; Zhu, C.; Zhou, F.; Tang, W.; Wu, D.; Li, Z.; Zhou, L.; Liu, J.; Wei, X.; et al. Co-circulation of multiple reassortant H6 subtype avian influenza viruses in wild birds in eastern China, 2016–2017. Virol. J. 2020, 17, 62. [Google Scholar] [CrossRef]
  72. Yu, Z.; Sun, W.; Li, X.; Chen, Q.; Chai, H.; Gao, X.; Guo, J.; Zhang, K.; Wang, T.; Feng, N.; et al. Adaptive amino acid substitutions enhance the virulence of a reassortant H7N1 avian influenza virus isolated from wild waterfowl in mice. Virology 2015, 476, 233–239. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00774 g001 550
Figure 1. Maximum likelihood phylogenetic trees showing inferred relationships among nucleotide sequences for the complete coding regions of K6 and K17 gene segments. Red open circle indicates K6-H6N5 and blue open circle indicates K17. These are highlighted with red and blue arrows used for K6 and K17, respectively. Closed and open circles indicate lineages of North America and Eurasia, respectively.
Figure 1. Maximum likelihood phylogenetic trees showing inferred relationships among nucleotide sequences for the complete coding regions of K6 and K17 gene segments. Red open circle indicates K6-H6N5 and blue open circle indicates K17. These are highlighted with red and blue arrows used for K6 and K17, respectively. Closed and open circles indicate lineages of North America and Eurasia, respectively.
Viruses 12 00774 g001
Viruses 12 00774 g002 550
Figure 2. In vitro growth properties of H6N5 in MDCK cells. Viruses were inoculated into MDCK monolayers at a multiplicity of infection (MOI) of 0.01. The supernatants were collected at different time points: 12, 24, 36, 48, 60, and 72 hours post infection (hpi). Viral replication titers were determined by TCID50 assay. *, p < 0.05; **, p < 0.01; ***, p < 0.001.
Figure 2. In vitro growth properties of H6N5 in MDCK cells. Viruses were inoculated into MDCK monolayers at a multiplicity of infection (MOI) of 0.01. The supernatants were collected at different time points: 12, 24, 36, 48, 60, and 72 hours post infection (hpi). Viral replication titers were determined by TCID50 assay. *, p < 0.05; **, p < 0.01; ***, p < 0.001.
Viruses 12 00774 g002
Viruses 12 00774 g003 550
Figure 3. Pathogenicity of K6-H6N5 virus in vivo. (A) Balb/c mice were intranasally challenged with different titers of the K6-H6N5 isolate. (B–F) H1N1 and H7N1 strains were used as controls, administered intranasally with 105 EID50. Mouse weight (B) and survival rate (C) were observed for 14 dpi. Body weight is presented as % of those of original mice (n = 10). (D) Mean viral titers in the lungs of mice (n = 3) were measured at 3, 6, and 14 dpi. (E,F) Histology of lung inflammation following virus infection was observed. (E) H&E staining of uninfected control (normal); H6N5; H1N1; H7N1-infected mouse lungs at 6 dpi (scale bar, 100 μm; original magnification, 100 ×). (F) Mean of lungs weight from an uninfected control (normal) and H1N1-, H7N1-, and H6N5-infected mice at 6 dpi (n = 3). *, p < 0.05; ***, p < 0.001
Figure 3. Pathogenicity of K6-H6N5 virus in vivo. (A) Balb/c mice were intranasally challenged with different titers of the K6-H6N5 isolate. (B–F) H1N1 and H7N1 strains were used as controls, administered intranasally with 105 EID50. Mouse weight (B) and survival rate (C) were observed for 14 dpi. Body weight is presented as % of those of original mice (n = 10). (D) Mean viral titers in the lungs of mice (n = 3) were measured at 3, 6, and 14 dpi. (E,F) Histology of lung inflammation following virus infection was observed. (E) H&E staining of uninfected control (normal); H6N5; H1N1; H7N1-infected mouse lungs at 6 dpi (scale bar, 100 μm; original magnification, 100 ×). (F) Mean of lungs weight from an uninfected control (normal) and H1N1-, H7N1-, and H6N5-infected mice at 6 dpi (n = 3). *, p < 0.05; ***, p < 0.001
Viruses 12 00774 g003
Viruses 12 00774 g004 550
Figure 4. Location of H6N5 isolates in Korea: A/Bean goose/South Korea/KNU18-6/2018(H6N5) and A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5). The map was retrieved from Google Earth.
Figure 4. Location of H6N5 isolates in Korea: A/Bean goose/South Korea/KNU18-6/2018(H6N5) and A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5). The map was retrieved from Google Earth.
Viruses 12 00774 g004
Viruses 12 00774 g005 550
Figure 5. Evolutionary hypothesis for the ancestor of each gene segment of A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5) (A) and A/Bean goose/South Korea/KNU18-6/2018(H6N5) (B).
Figure 5. Evolutionary hypothesis for the ancestor of each gene segment of A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5) (A) and A/Bean goose/South Korea/KNU18-6/2018(H6N5) (B).
Viruses 12 00774 g005
Table
Table 1. Sequence homology of the whole A/Bean goose/South Korea/KNU18-6/2018(H6N5) genome.
Table 1. Sequence homology of the whole A/Bean goose/South Korea/KNU18-6/2018(H6N5) genome.
GeneGenBank Accession # of K6Strain Having Highest Nucleotide Identity with K6Genetic Identity between
K6 and K17
StrainGenBank Accession #
(% Genetic Identity)
PB2MN749559A/Mallard/Alaska/AH0029066S.1.A/2016(H12N5)MN254630.1
(99.34%)		99.21%
PB1MN749560A/Mallard /Alaska/AH0029066S.2.A /2016(H12N5)MN254483.1
(99.56%)		94.64%
PAMN749561A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5)MK830102.1
(99.49%)		99.49%
HAMN749562A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5)EPI1526581
99.30%		99.30%
NPMN749563A/lesser scaup/Wisconsin/17OS5811/2017(H6N1)EPI1537365
99.67%		96.32%
NAMN749564A/Aix galericulata/South Korea/K17-1638-5/2017(H6N5)EPI1526586
99.16%		99.16%
MPMN749565A/mallard/Kentucky/AH00729344/2016(H7N2)EPI952771
99.59%		97.66%
NSMN749566A/blue-wingedteal/Missouri/17OS3211/2017(A/H3N1)EPI1537413
99.52%		99.40%

Share and Cite

MDPI and ACS Style

Nguyen, N.M.; Sung, H.W.; Yun, K.-J.; Park, H.; 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. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070774

AMA Style

Nguyen NM, Sung HW, Yun K-J, Park H, 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(7):774. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070774

Chicago/Turabian Style

Nguyen, Ngoc Minh, Haan Woo Sung, Ki-Jung Yun, Hyun Park, and Seon-Ju Yeo. 2020. "Genetic Characterization of a Novel North American-Origin Avian Influenza A (H6N5) Virus Isolated from Bean Goose of South Korea in 2018" Viruses 12, no. 7: 774. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070774

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop